scholarly journals Consistent Velocity–Pressure Coupling for Second-Order L2-Penalty and Direct-Forcing Methods

Fluids ◽  
2020 ◽  
Vol 5 (2) ◽  
pp. 92
Author(s):  
Arthur Sarthou ◽  
Stéphane Vincent ◽  
Jean-Paul Caltagirone

The present work studies the interactions between fictitious-domain methods on structured grids and velocity–pressure coupling for the resolution of the Navier–Stokes equations. The pressure-correction approaches are widely used in this context but the corrector step is generally not modified consistently to take into account the fictitious domain. A consistent modification of the pressure-projection for a high-order penalty (or penalization) method close to the Ikeno–Kajishima modification for the Immersed Boundary Method is presented here. Compared to the first-order correction required for the L 2 -penalty methods, the small values of the penalty parameters do not lead to numerical instabilities in solving the Poisson equation. A comparison of the corrected rotational pressure-correction method with the augmented Lagrangian approach which does not require a correction is carried out.

2015 ◽  
Vol 31 (6) ◽  
pp. 771-782 ◽  
Author(s):  
R.-J. Wu ◽  
S.-Y. Lin

AbstractA modified direct-forcing immersed-boundary (IB) pressure correction method is developed to simulate the flows of a falling ellipse. The pressure correct method is used to solve the solutions of the two dimensional Navier-Stokes equations and a direct-forcing IB method is used to handle the interaction between the flow and falling ellipse. For a fixed aspect ratio of an ellipse, the types of the behavior of the falling ellipse can be classified as three pure motions: Steady falling, fluttering, tumbling, and two transition motions: Chaos, transition between steady falling and fluttering. Based on two dimensionless parameters, Reynolds number and the dimensionless moment of inertia, a Reynolds number-inertia moment phase diagram is established. The behaviors and characters of five falling regimes are described in detailed.


Author(s):  
Neil W. Harvey ◽  
Martin G. Rose ◽  
John Coupland ◽  
Terry Jones

A 3-D steady viscous finite volume pressure correction method for the solution of the Reynolds averaged Navier-Stokes equations has been used to calculate the heat transfer rates on the end walls of a modern High Pressure Turbine first stage stator. Surface heat transfer rates have been calculated at three conditions and compared with measurements made on a model of the vane tested in annular cascade in the Isentropic Light Piston Facility at DERA, Pyestock. The NGV Mach numbers, Reynolds numbers and geometry are fully representative of engine conditions. Design condition data has previously been presented by Harvey and Jones (1990). Off-design data is presented here for the first time. In the areas of highest heat transfer the calculated heat transfer rates are shown to be within 20% of the measured values at all three conditions. Particular emphasis is placed on the use of wall functions in the calculations with which relatively coarse grids (of around 140,000 nodes) can be used to keep computational run times sufficiently low for engine design purposes.


2018 ◽  
Vol 26 (2) ◽  
pp. 267-283
Author(s):  
M. Tezer-Sezgin ◽  
Merve Gürbüz

Abstract We consider the steady, laminar, convection ow in a long channel of 2D rectangular constricted cross-section under the inuence of an applied magnetic field. The Navier-Stokes equations including Lorentz and buoyancy forces are coupled with the temperature equation and are solved by using linear radial basis function (RBF) approximations in terms of the velocity, pressure and the temperature of the fluid. RBFs are used in the approximation of the particular solution which becomes also the approximate solution of the problem. Results are obtained for several values of Grashof number (Gr), Hartmann number (M) and the constriction ratios (CR) to see the effects on the ow and isotherms for fixed values of Reynolds number and Prandtl number. As M increases, the ow is flattened. An increase in Gr increases the magnitude of the ow in the channel. Isolines undergo an inversion at the center of the channel indicating convection dominance due to the strong buoyancy force, but this inversion is retarded with the increase in the strength of the applied magnetic field. When both Hartmann number and constriction ratio are increased, ow is divided into more loops symmetrically with respect to the axes.


Sign in / Sign up

Export Citation Format

Share Document