scholarly journals PARALLEL COMPUTATION OF ROTATING FLOWS

1999 ◽  
Vol 4 (1) ◽  
pp. 124-134
Author(s):  
L. K. Lundin ◽  
V. A. Barker ◽  
J. N. Sørensen
Keyword(s):  
Stokes Equations ◽  
Rotating Flows ◽  
Navier Stokes ◽  
Second Step ◽  
Governing Equations ◽  
Time Level ◽  
Navier Stokes Equations ◽  
Linear System Of Equations ◽  
Wide Range ◽  
New Time

This paper deals with the simulation of 3‐D rotating flows based on the velocity‐vorticity formulation of the Navier‐Stokes equations in cylindrical coordinates. The governing equations are discretized by a finite difference method. The solution is advanced to a new time level by a two‐step process. In the first step, the vorticity at the new time level is computed using the velocity at the previous time level. In the second step, the velocity at the new time level is computed using the new vorticity. We discuss here the second part which is by far the most time‐consuming. The numerical problem is that of solving a singular, large, sparse, over‐determined linear system of equations, and the iterative method CGLS is applied for this purpose. We discuss some of the mathematical and numerical aspects of this procedure and report on the performance of our software on a wide range of parallel computers.

The Effect of Electrostatic Forces on the Distribution of Drops in a Channel

2002 ◽  
Author(s):  
Arturo Ferna´ndez ◽  
Jiacai Lu ◽  
Asghar Esmaeeli ◽  
Gre´tar Tryggvason
Keyword(s):  
Electric Field ◽  
Electric Fields ◽  
Hydrodynamic Interaction ◽  
Stokes Equations ◽  
Navier Stokes ◽  
Governing Equations ◽  
Navier Stokes Equations ◽  
Dielectric Fluids ◽  
Wide Range ◽  
Leaky Dielectric Model

Direct numerical simulations are used to examine the effect of electric fields on the behavior of suspension of drops in dielectric fluids. The effect of electric field is modeled using the “leaky dielectric” model, coupled with the full Navier-Stokes equations. The governing equations are solved using a front-tracking/finite volume technique. The interaction of the drops is strongly dependant on the conductivity and the permittivity ratio, but fibration, where drops line up into long columns, takes place over a wide range of these parameters. The hydrodynamic interaction due to fluid circulation induced by the electric field has a strong influence on the drop distribution and the rate of fibration.

scholarly journals Thermochemical Nonequilibrium in Rapidly Expanding Flows of High-Temperature Air

1997 ◽  
Vol 52 (4) ◽  
pp. 358-368 ◽  
Author(s):  
Michio Nishida ◽  
Masashi Matsumotob
Keyword(s):  
High Temperature ◽  
Vibrational Energy ◽  
Stokes Equations ◽  
Computational Study ◽  
Mass Conservation ◽  
Navier Stokes ◽  
Tvd Scheme ◽  
Governing Equations ◽  
Navier Stokes Equations ◽  
Free Jet

Abstract • This paper describes a computational study of the thermal and chemical nonequilibrium occuring in a rapidly expanding flow of high-temperature air transported as a free jet from an orifice into low-density stationary air. Translational, rotational, vibrational and electron temperatures are treated separately, and in particular the vibrational temperatures are individually treated; a multi-vibrational temperature model is adopted. The governing equations are axisymmetric Navier-Stokes equations coupled with species vibrational energy, electron energy and species mass conservation equations. These equations are numerically solved, using the second order upwind TVD scheme of the Harten-Yee type. The calculations were carried out for two different orifice temperatures and also two different orifice diameters to investigate the effects of such parameters on the structure of a nonequilibrium free jet.

Numerical Investigation of Blade Flutter at or Near Stall in Axial Flow Turbomachines

2001 ◽  
Author(s):  
Wolfgang Höhn
Keyword(s):  
Viscous Flow ◽  
Stokes Equations ◽  
Separated Flow ◽  
Parallel Computer ◽  
Navier Stokes ◽  
Computational Domain ◽  
Compressor Cascade ◽  
Governing Equations ◽  
Navier Stokes Equations ◽  
Blade Flutter

During the design of the compressor and turbine stages of today’s aeroengines, aerodynamically induced vibrations become increasingly important since higher blade load and better efficiency are desired. In this paper the development of a method based on the unsteady, compressible Navier-Stokes equations in two dimensions is described in order to study the physics of flutter for unsteady viscous flow around cascaded vibrating blades at stall. The governing equations are solved by a finite difference technique in boundary fitted coordinates. The numerical scheme uses the Advection Upstream Splitting Method to discretize the convective terms and central differences discretizing the viscous terms of the fully non-linear Navier-Stokes equations on a moving H-type mesh. The unsteady governing equations are explicitly and implicitly marched in time in a time-accurate way using a four stage Runge-Kutta scheme on a parallel computer or an implicit scheme of the Beam-Warming type on a single processor. Turbulence is modelled using the Baldwin-Lomax turbulence model. The blade flutter phenomenon is simulated by imposing a harmonic motion on the blade, which consists of harmonic body translation in two directions and a rotation, allowing an interblade phase angle between neighboring blades. Non-reflecting boundary conditions are used for the unsteady analysis at inlet and outlet of the computational domain. The computations are performed on multiple blade passages in order to account for nonlinear effects. A subsonic massively stalled unsteady flow case in a compressor cascade is studied. The results, compared with experiments and the predictions of other researchers, show reasonable agreement for inviscid and viscous flow cases for the investigated flow situations with respect to the Steady and unsteady pressure distribution on the blade in separated flow areas as well as the aeroelastic damping. The results show the applicability of the scheme for stalled flow around cascaded blades. As expected the viscous and inviscid computations show different results in regions where viscous effects are important, i.e. in separated flow areas. In particular, different predictions for inviscid and viscous flow for the aerodynamic damping for the investigated flow cases are found.

On the flow near the trailing edge of a flat plate

1968 ◽  
Vol 306 (1486) ◽  
pp. 275-290 ◽  
Keyword(s):  
Boundary Layer ◽  
Flat Plate ◽  
Stokes Equations ◽  
Navier Stokes ◽  
Governing Equations ◽  
Trailing Edge ◽  
Navier Stokes Equations ◽  
Multiplicative Constant ◽  
Boundary Layer Approximation ◽  
Uniform Shear

It is shown that the boundary layer approximation to the flow of a viscous fluid past a flat plate of length l , generally valid near the plate when the Reynolds number Re is large, fails within a distance O( lRe -3/4 ) of the trailing edge. The appropriate governing equations in this neighbourhood are the full Navier- Stokes equations. On the basis of Imai (1966) these equations are linearized with respect to a uniform shear and are then completely solved by means of a Wiener-Hopf integral equation. The solution so obtained joins smoothly on to that of the boundary layer for a flat plate upstream of the trailing edge and for a wake downstream of the trailing edge. The contribution to the drag coefficient is found to be O ( Re -3/4 ) and the multiplicative constant is explicitly worked out for the linearized equations.

scholarly journals Gyrotactic swimmers in turbulence: shape effects and role of the large-scale flow

2018 ◽  
Vol 856 ◽  
Author(s):  
M. Borgnino ◽  
G. Boffetta ◽  
F. De Lillo ◽  
M. Cencini
Keyword(s):  
Large Scale ◽  
Stokes Equations ◽  
Navier Stokes ◽  
Small Scale ◽  
Navier Stokes Equations ◽  
Preferential Sampling ◽  
Wide Range ◽  
Shape Effects ◽  
Dilute Suspensions ◽  
Large Scale Flow

We study the dynamics and the statistics of dilute suspensions of gyrotactic swimmers, a model for many aquatic motile microorganisms. By means of extensive numerical simulations of the Navier–Stokes equations at different Reynolds numbers, we investigate preferential sampling and small-scale clustering as a function of the swimming (stability and speed) and shape parameters, considering in particular the limits of spherical and rod-like particles. While spherical swimmers preferentially sample local downwelling flow, for elongated swimmers we observe a transition from downwelling to upwelling regions at sufficiently high swimming speed. The spatial distribution of both spherical and elongated swimmers is found to be fractal at small scales in a wide range of swimming parameters. The direct comparison between the different shapes shows that spherical swimmers are more clusterized at small stability and speed numbers, while for large values of the parameters elongated cells concentrate more. The relevance of our results for phytoplankton swimming in the ocean is briefly discussed.

Complex dynamics in a stratified lid-driven square cavity flow

2018 ◽  
Vol 855 ◽  
pp. 43-66 ◽  
Author(s):  
Ke Wu ◽  
Bruno D. Welfert ◽  
Juan M. Lopez
Keyword(s):  
Stokes Equations ◽  
Complex Dynamics ◽  
Boussinesq Approximation ◽  
Navier Stokes ◽  
Square Cavity ◽  
Navier Stokes Equations ◽  
Stable Temperature ◽  
Wide Range ◽  
Side Walls ◽  
Steady State Solutions

The dynamic response to shear of a fluid-filled square cavity with stable temperature stratification is investigated numerically. The shear is imposed by the constant translation of the top lid, and is quantified by the associated Reynolds number. The stratification, quantified by a Richardson number, is imposed by maintaining the temperature of the top lid at a higher constant temperature than that of the bottom, and the side walls are insulating. The Navier–Stokes equations under the Boussinesq approximation are solved, using a pseudospectral approximation, over a wide range of Reynolds and Richardson numbers. Particular attention is paid to the dynamical mechanisms associated with the onset of instability of steady state solutions, and to the complex and rich dynamics occurring beyond.

scholarly journals Investigating the Effect of the Closure in Partially-Averaged Navier–Stokes Equations

2019 ◽  
Vol 141 (12) ◽  
Author(s):  
Filipe S. Pereira ◽  
Luís Eça ◽  
Guilherme Vaz
Keyword(s):  
Reynolds Number ◽  
Shear Stress ◽  
Coherent Structures ◽  
Stokes Equations ◽  
The Other ◽  
Navier Stokes ◽  
Governing Equations ◽  
Navier Stokes Equations ◽  
Auxiliary Functions ◽  
Modeling Accuracy

The importance of the turbulence closure to the modeling accuracy of the partially-averaged Navier–Stokes equations (PANS) is investigated in prediction of the flow around a circular cylinder at Reynolds number of 3900. A series of PANS calculations at various degrees of physical resolution is conducted using three Reynolds-averaged Navier–Stokes equations (RANS)-based closures: the standard, shear-stress transport (SST), and turbulent/nonturbulent (TNT) k–ω models. The latter is proposed in this work. The results illustrate the dependence of PANS on the closure. At coarse physical resolutions, a narrower range of scales is resolved so that the influence of the closure on the simulations accuracy increases significantly. Among all closures, PANS–TNT achieves the lowest comparison errors. The reduced sensitivity of this closure to freestream turbulence quantities and the absence of auxiliary functions from its governing equations are certainly contributing to this result. It is demonstrated that the use of partial turbulence quantities in such auxiliary functions calibrated for total turbulent (RANS) quantities affects their behavior. On the other hand, the successive increase of physical resolution reduces the relevance of the closure, causing the convergence of the three models toward the same solution. This outcome is achieved once the physical resolution and closure guarantee the precise replication of the spatial development of the key coherent structures of the flow.

scholarly journals Numerical Simulation of Micromixing of Particles and Fluids with Galloping Cylinder

Symmetry ◽  
2020 ◽  
Vol 12 (4) ◽  
pp. 580 ◽  
Author(s):  
Zahra Abdelmalek ◽  
Mohammad Yaghoub Abdollahzadeh Jamalabadi
Keyword(s):  
Continuity Equation ◽  
Stokes Equations ◽  
The Self ◽  
Moving Mesh ◽  
Navier Stokes ◽  
Parameter Study ◽  
Mems Devices ◽  
Governing Equations ◽  
Navier Stokes Equations ◽  
Moving Cylinder

Micromixers are significant segments inside miniaturized scale biomedical frameworks. Numerical investigation of the effects of galloping cylinder characteristics inside a microchannel Newtonian, incompressible fluid in nonstationary condition is performed. Governing equations of the system include the continuity equation, and Navier–Stokes equations are solved within a moving mesh domain. The symmetry of laminar entering the channel is broken by the self-sustained motion of the cylinder. A parameter study on the amplitude and frequency of passive moving cylinder on the mixing of tiny particles in the fluid is performed. The results show a significant increase to the index of mixing uses of the galloping body in biomedical frameworks in the course of micro-electromechanical systems (MEMS) devices.

Heat Transfer From a Circular Cylinder at Low Reynolds Numbers

1998 ◽  
Vol 120 (1) ◽  
pp. 72-75 ◽  
Author(s):  
V. N. Kurdyumov ◽  
E. Ferna´ndez
Keyword(s):  
Circular Cylinder ◽  
Energy Equation ◽  
Numerical Solutions ◽  
Stokes Equations ◽  
Navier Stokes ◽  
Order Unity ◽  
Navier Stokes Equations ◽  
Wide Range ◽  
Prandtl Numbers ◽  
High Degree

A correlation formula, Nu = W0(Re)Pr1/3 + W1(Re), that is valid in a wide range of Reynolds and Prandtl numbers has been developed based on the asymptotic expansion for Pr → ∞ for the forced heat convection from a circular cylinder. For large Prandtl numbers, the boundary layer theory for the energy equation is applied and compared with the numerical solutions of the full Navier Stokes equations for the flow field and energy equation. It is shown that the two-terms asymptotic approximation can be used to calculate the Nusselt number even for Prandtl numbers of order unity to a high degree of accuracy. The formulas for coefficients W0 and W1, are provided.

Lift, drag and added mass of a hemispherical bubble sliding and growing on a wall in a viscous linear shear flow

2008 ◽  
Vol 366 (1873) ◽  
pp. 2233-2248 ◽  
Author(s):  
Dominique Legendre ◽  
Catherine Colin ◽  
Typhaine Coquard
Keyword(s):  
Shear Flow ◽  
Stokes Equations ◽  
Three Dimensional ◽  
Added Mass ◽  
Unsteady Motion ◽  
Navier Stokes ◽  
Mass Force ◽  
Navier Stokes Equations ◽  
Wide Range ◽  
Linear Shear

The three-dimensional flow around a hemispherical bubble sliding and growing on a wall in a viscous linear shear flow is studied numerically by solving the full Navier–Stokes equations in a boundary-fitted domain. The main goal of the present study is to provide a complete description of the forces experienced by the bubble (drag, lift and added mass) over a wide range of sliding and shear Reynolds numbers (0.01≤ Re b , Re α ≤2000) and shear rate (0≤ Sr ≤5). The drag and lift forces are computed successively for the following situations: an immobile bubble in a linear shear flow; a bubble sliding on the wall in a fluid at rest; and a bubble sliding in a linear shear flow. The added-mass force is studied by considering an unsteady motion relative to the wall or a time-dependent radius.

Sign in / Sign up

Export Citation Format

Share Document