Abstract: A nonlinear model of the mathematical fluid dynamics which describes the motion of an incompressible viscous rotating fluid in a homogeneous gravitational field is considered. The model is a generalization of the known Navier-Stokes system with the addition of the Coriolis parameter and the equations for changeable density. An explicit algorithm for the solution is constructed, and the proof of the existence and uniqueness theorems for the strong solution of the nonlinear problem is given. For the linear case, the localization and the structure of the spectrum of inner waves are also investigated.
Abstract: This article presents a numerical method to find the
heat flux in an inhomogeneous inverse heat conduction problem with
linear boundary conditions and an extra specification at the terminal.
The method is based upon applying the satisfier function along with
the Ritz-Galerkin technique to reduce the approximate solution of the
inverse problem to the solution of a system of algebraic equations.
The instability of the problem is resolved by taking advantage of
the Landweber’s iterations as an admissible regularization strategy.
In computations, we find the stable and low-cost results which
demonstrate the efficiency of the technique.
Abstract: This paper develops a meshless approach, called Element Free Galerkin (EFG) method, which is based on the weak form Moving Least Squares (MLS) of the partial differential governing equations and employs the interpolation to construct the meshless shape functions. The variation weak form is used in the EFG where the trial and test functions are approximated bye the MLS approximation. Since the shape functions constructed by this discretization have the weight function property based on the randomly distributed points, the essential boundary conditions can be implemented easily. The local weak form of the partial differential governing equations is obtained by the weighted residual method within the simple local quadrature domain. The spline function with high continuity is used as the weight function. The presently developed EFG method is a truly meshless method, as it does not require the mesh, either for the construction of the shape functions, or for the integration of the local weak form. Several numerical examples of two-dimensional static structural analysis are presented to illustrate the performance of the present EFG method. They show that the EFG method is highly efficient for the implementation and highly accurate for the computation. The present method is used to analyze the static deflection of beams and plate hole
Abstract: In this paper, the influence of van der Waals, as well as electrostatic forces on the structural behavior of MEMS and NEMS actuators, has been investigated using of a Euler-Bernoulli beam continuous model. In the proposed nonlinear model, the electrostatic fringing-fields and the mid-plane stretching (geometric nonlinearity) effects have been considered. The nonlinear integro-differential equation governing the static structural behavior of the actuator has been derived. An original Galerkin-based reduced-order model has been developed to avoid problems arising from the nonlinearities in the differential equation. The obtained reduced-order model equations have been solved numerically using the Newton-Raphson method. The basic design parameters such as the pull-in parameters (voltage and deflection at pull-in), as well as the detachment length due to the van der Waals force of some investigated micro- and nano-actuators have been calculated. The obtained numerical results have been compared with some other existing methods (finite-elements method and finite-difference method) and the comparison showed good agreement among all assumed numerical techniques.
Abstract: This paper investigates the parametric stability of an
axially moving web subjected to non-uniform in-plane edge
excitations on two opposite, simply-supported edges. The web is
modeled as a viscoelastic plate whose constitutive relation obeys the
Kelvin-Voigt model, and the in-plane edge excitations are expressed
as the sum of a static tension and a periodical perturbation. Due to the
in-plane edge excitations, the moving plate may bring about
parametric instability under certain situations. First, the in-plane
stresses of the plate due to the non-uniform edge excitations are
determined by solving the in-plane forced vibration problem. Then,
the dependence on the spatial coordinates in the equation of transverse
motion is eliminated by the generalized Galerkin method, which
results in a set of discretized system equations in time. Finally, the
method of multiple scales is utilized to solve the set of system
equations analytically if the periodical perturbation of the in-plane
edge excitations is much smaller as compared with the static tension of
the plate, from which the stability boundaries of the moving plate are
obtained. Numerical results reveal that only combination resonances
of the summed-type appear under the in-plane edge excitations
considered in this work.
Abstract: New physical insights into the nonlinear Lorenz
equations related to flow resistance is discussed in this work. The
chaotic dynamics related to Lorenz equations has been studied in
many papers, which is due to the sensitivity of Lorenz equations to
initial conditions and parameter uncertainties. However, the physical
implication arising from Lorenz equations about convectional motion
attracts little attention in the relevant literature. Therefore, as a first
step to understand the related fluid mechanics of convectional motion,
this paper derives the Lorenz equations again with different forced
conditions in the model. Simulation work of the modified Lorenz
equations without the viscosity or buoyancy force is discussed. The
time-domain simulation results may imply that the states of the
Lorenz equations are related to certain flow speed and flow resistance.
The flow speed of the underlying fluid system increases as the flow
resistance reduces. This observation would be helpful to analyze the
coupling effects of different fluid parameters in a convectional model
in future work.
Abstract: Stator elements «Vane diffuser + crossover + return
channel» of stages with different specific speed were investigated by
CFD calculations. The regime parameter was introduced to present
efficiency and loss coefficient performance of all elements together.
Flow structure demonstrated advantages and disadvantages of design.
Flow separation in crossovers was eliminated by its shape
modification. Efficiency increased visibly. Calculated CFD
performances are in acceptable correlation with predicted ones by
engineering design method. The information obtained is useful for
design method better calibration.
Abstract: The 6th version of Universal modeling method for
centrifugal compressor stage calculation is described. Identification
of the new mathematical model was made. As a result of
identification the uniform set of empirical coefficients is received.
The efficiency definition error is 0,86 % at a design point. The
efficiency definition error at five flow rate points (except a point of
the maximum flow rate) is 1,22 %. Several variants of the stage with
3D impellers designed by 6th version program and quasi threedimensional
calculation programs were compared by their gas
dynamic performances CFD (NUMECA FINE TURBO).
Performance comparison demonstrated general principles of design
validity and leads to some design recommendations.
Abstract: Parameters of flow are calculated in vaneless diffusers
with relative width 0,014–0,10. Inlet angles of flow and similarity
criteria were varied. There is information on flow separation,
boundary layer development, configuration of streamlines.
Polytrophic efficiency, loss coefficient and recovery coefficient are
used to compare effectiveness of diffusers. The sample of
optimization of narrow diffuser with conical walls is presented. Three
wide diffusers with narrowing walls are compared. The work is made
in the R&D laboratory “Gas dynamics of turbo machines” of the TU
SPb.
Abstract: Universal modeling method well proven for industrial
compressors was applied for design of the high flow rate supersonic
stage. Results were checked by ANSYS CFX and NUMECA Fine
Turbo calculations. The impeller appeared to be very effective at
transonic flow velocities. Stator elements efficiency is acceptable at
design Mach numbers too. Their loss coefficient versus inlet flow
angle performances correlates well with Universal modeling
prediction. The impeller demonstrates ability of satisfactory operation
at design flow rate. Supersonic flow behavior in the impeller inducer
at the shroud blade to blade surface Φ des deserves additional study.
Abstract: There is an evident trend to elevate pressure ratio of a
single stage of a turbo compressors - axial compressors in particular.
Whilst there was an opinion recently that a pressure ratio 1,9 was a
reasonable limit, later appeared information on successful modeling
tested of stages with pressure ratio up to 2,8. The authors recon that
lack of information on high pressure stages makes actual a study of
rational choice of design parameters before high supersonic flow
problems solving. The computer program of an engineering type was
developed. Below is presented a sample of its application to study
possible parameters of the impeller of the stage with pressure ratio
3,0. Influence of two main design parameters on expected efficiency,
periphery blade speed and flow structure is demonstrated. The results
had lead to choose a variant for further analysis and improvement by
CFD methods.
Abstract: In this paper, a numerical algorithm using a coupled Galerkin-Differential Quadrature (DQ) method is proposed for the solution of dam-reservoir interaction problem. The governing differential equation of motion of the dam structure is discretized by the Galerkin method and the DQM is used to discretize the fluid domain. The resulting systems of ordinary differential equations are then solved by the Newmark time integration scheme. The mixed scheme combines the simplicity of the Galerkin method and high accuracy and efficiency of the DQ method. Its accuracy and efficiency are demonstrated by comparing the calculated results with those of the existing literature. It is shown that highly accurate results can be obtained using a small number of Galerkin terms and DQM sampling points. The technique presented in this investigation is general and can be used to solve various fluid-structure interaction problems.
Abstract: Two-dimensional Eulerian (volume-averaged) continuity and momentum equations governing multi-size slurry flow through pump casings are solved by applying a penalty finite element formulation. The computational strategy validated for multi-phase flow through rectangular channels is adapted to the present study. The flow fields of the carrier, mixture and each solids species, and the concentration field of each species are determined sequentially in an iterative manner. The eddy viscosity field computed using Spalart-Allmaras model for the pure carrier phase is modified for the presence of particles. Streamline upwind Petrov-Galerkin formulation is used for all the momentum equations for the carrier, mixture and each solids species and the concentration field for each species. After ensuring mesh-independence of solutions, results of multi-size particulate flow simulation are presented to bring out the effect of bulk flow rate, average inlet concentration, and inlet particle size distribution. Mono-size computations using (1) the concentration-weighted mean diameter of the slurry and (2) the D50 size of the slurry are also presented for comparison with multi-size results.
Abstract: A numerical approach for solving constant-coefficient differential equations whose solutions exhibit boundary layer structure is built by inserting Bernstein Partition of Unity into Galerkin variational weak form. Due to the reproduction capability of Bernstein basis, such implementation shows excellent accuracy at boundaries and is able to capture sharp gradients of the field variable by p-refinement using regular distributions of equi-spaced evaluation points. The approximation is subjected to convergence experimentation and a procedure to assemble the discrete equations without a background integration mesh is proposed.
Abstract: A theoretical investigation on the effects of both
steady-state and dynamic deformations of the foils on the dynamic
performance characteristics of a self-acting air foil journal bearing
operating under small harmonic vibrations is proposed. To take into
account the dynamic deformations of foils, the perturbation method is
used for determining the gas-film stiffness and damping coefficients
for given values of excitation frequency, compressibility number, and
compliance factor of the bump foil. The nonlinear stationary
Reynolds’ equation is solved by means of the Galerkins’ finite
element formulation while the finite differences method are used to
solve the first order complex dynamic equations resulting from the
perturbation of the nonlinear transient compressible Reynolds’
equation. The stiffness of a bump is uniformly distributed throughout
the bearing surface (generation I bearing). It was found that the
dynamic properties of the compliant finite length journal bearing are
significantly affected by the compliance of foils especially whenthe
dynamic deformation of foils is considered in addition to the static
one by applying the principle of superposition.
Abstract: The effect of internal heat generation is applied to the Rayleigh-Benard convection in a horizontal micropolar fluid layer. The bounding surfaces of the liquids are considered to be rigid-free, rigid-rigid and free-free with the combination of isothermal on the spin-vanishing boundaries. A linear stability analysis is used and the Galerkin method is employed to find the critical stability parameters numerically. It is shown that the critical Rayleigh number decreases as the value of internal heat generation increase and hence destabilize the system.
Abstract: This paper deals with a high-order accurate Runge
Kutta Discontinuous Galerkin (RKDG) method for the numerical
solution of the wave equation, which is one of the simple case of a
linear hyperbolic partial differential equation. Nodal DG method is
used for a finite element space discretization in 'x' by discontinuous
approximations. This method combines mainly two key ideas which
are based on the finite volume and finite element methods. The
physics of wave propagation being accounted for by means of
Riemann problems and accuracy is obtained by means of high-order
polynomial approximations within the elements. High order accurate
Low Storage Explicit Runge Kutta (LSERK) method is used for
temporal discretization in 't' that allows the method to be nonlinearly
stable regardless of its accuracy. The resulting RKDG
methods are stable and high-order accurate. The L1 ,L2 and L∞ error
norm analysis shows that the scheme is highly accurate and effective.
Hence, the method is well suited to achieve high order accurate
solution for the scalar wave equation and other hyperbolic equations.
Abstract: In this paper static and dynamic response of a
varactor of a micro-phase shifter to DC, step DC and AC
voltages have been studied. By presenting a mathematical
modeling Galerkin-based step by step linearization method
(SSLM) and Galerkin-based reduced order model have been
used to solve the governing static and dynamic equations,
respectively. The calculated static and dynamic pull-in
voltages have been validated by previous experimental and
theoretical results and a good agreement has been achieved.
Then the frequency response and phase diagram of the system
has been studied. It has been shown that applying the DC
voltage shifts down the phase diagram and frequency
response. Also increasing the damping ratio shifts up the
phase diagram.
Abstract: The simulation of extrusion process is studied widely
in order to both increase products and improve quality, with broad
application in wire coating. The annular tube-tooling extrusion was
set up by a model that is termed as Navier-Stokes equation in
addition to a rheological model of differential form based on singlemode
exponential Phan-Thien/Tanner constitutive equation in a twodimensional
cylindrical coordinate system for predicting the
contraction point of the polymer melt beyond the die. Numerical
solutions are sought through semi-implicit Taylor-Galerkin pressurecorrection
finite element scheme. The investigation was focused on
incompressible creeping flow with long relaxation time in terms of
Weissenberg numbers up to 200. The isothermal case was considered
with surface tension effect on free surface in extrudate flow and no
slip at die wall. The Stream Line Upwind Petrov-Galerkin has been
proposed to stabilize solution. The structure of mesh after die exit
was adjusted following prediction of both top and bottom free
surfaces so as to keep the location of contraction point around one
unit length which is close to experimental results. The simulation of
extrusion process is studied widely in order to both increase products
and improve quality, with broad application in wire coating. The
annular tube-tooling extrusion was set up by a model that is termed
as Navier-Stokes equation in addition to a rheological model of
differential form based on single-mode exponential Phan-
Thien/Tanner constitutive equation in a two-dimensional cylindrical
coordinate system for predicting the contraction point of the polymer
melt beyond the die. Numerical solutions are sought through semiimplicit
Taylor-Galerkin pressure-correction finite element scheme.
The investigation was focused on incompressible creeping flow with
long relaxation time in terms of Weissenberg numbers up to 200. The
isothermal case was considered with surface tension effect on free
surface in extrudate flow and no slip at die wall. The Stream Line
Upwind Petrov-Galerkin has been proposed to stabilize solution. The
structure of mesh after die exit was adjusted following prediction of
both top and bottom free surfaces so as to keep the location of
contraction point around one unit length which is close to
experimental results.
Abstract: Computational simulation of steam flow and heat transfer in power plant condensers on the basis of the threedimensional mathematical model for the flow through porous media is presented. In order to solve the mathematical model of steam flow and heat transfer in power plant condensers, the Streamline Upwind Petrov-Galerkin finite element method is applied. By comparison of the results of simulation with experimental results about an experimental condenser, it is confirmed that SUPG finite element method can be successfully applied for solving the three-dimensional mathematical model of steam flow and heat transfer in power plant condensers.