-
A preconditioning for the spectral solution of incompressible variable-density flows
Authors:
L. Reynier,
Bastien Di Pierro,
Frédéric Alizard,
Anne Cadiou,
Lionel Le Penven,
Marc Buffat
Abstract:
In the present study, the efficiency of preconditioners for solving linear systems associated with the discretized variable-density incompressible Navier-Stokes equations with semiimplicit second-order accuracy in time and spectral accuracy in space is investigated. The method, in which the inverse operator for the constant-density flow system acts as preconditioner, is implemented for three itera…
▽ More
In the present study, the efficiency of preconditioners for solving linear systems associated with the discretized variable-density incompressible Navier-Stokes equations with semiimplicit second-order accuracy in time and spectral accuracy in space is investigated. The method, in which the inverse operator for the constant-density flow system acts as preconditioner, is implemented for three iterative solvers: the General Minimal Residual, the Conjugate Gradient and the Richardson Minimal Residual. We discuss the method, first, in the context of the one-dimensional flow case where a top-hat like profile for the density is used. Numerical evidence shows that the convergence is significantly improved due to the notable decrease in the condition number of the operators. Most importantly, we then validate the robustness and convergence properties of the method on two more realistic problems: the two-dimensional Rayleigh-Taylor instability problem and the three-dimensional variable-density swirling jet.
△ Less
Submitted 11 March, 2024;
originally announced March 2024.
-
Stability of variable density rotating flows: Inviscid case and viscous effects in the limit of large Reynolds numbers
Authors:
C. Jacques,
B. Di Pierro,
F. Alizard,
M. Buffat,
A. Cadiou,
L. Le Penven
Abstract:
Linear stability of solid body rotating flows with axisymmetric density variations is addressed analytically. Considering inviscid disturbances, a non trivial dispersion relation is obtained and it is shown that the instability is of Rayleigh-Taylor type in cylindrical frame. The viscous correction is derived, in the limit of large Reynolds numbers and large azimuthal wave numbers, allowing the de…
▽ More
Linear stability of solid body rotating flows with axisymmetric density variations is addressed analytically. Considering inviscid disturbances, a non trivial dispersion relation is obtained and it is shown that the instability is of Rayleigh-Taylor type in cylindrical frame. The viscous correction is derived, in the limit of large Reynolds numbers and large azimuthal wave numbers, allowing the determination of the most unstable mode. Theoretical predictions are checked by comparing them to (spectrally accurate) computed eigenvalues and direct numerical simulations.
△ Less
Submitted 23 August, 2023;
originally announced August 2023.
-
Analysis of transient growth using an orthogonal decomposition of the velocity field in the Orr-Sommerfeld Squire equations
Authors:
Marc Buffat,
Lionel Le Penven
Abstract:
Despite remarkable accomplishment, the classical hydrodynamic stability theory fails to predict transition in wall-bounded shear ow. The shortcoming of this modal approach was found 20 years ago and is linked to the non-orthogonality of the eigenmodes of the linearised problem, de noted by the Orr Sommerfeld and Squire equations. The associated eigenmodes of this linearised problem are the normal…
▽ More
Despite remarkable accomplishment, the classical hydrodynamic stability theory fails to predict transition in wall-bounded shear ow. The shortcoming of this modal approach was found 20 years ago and is linked to the non-orthogonality of the eigenmodes of the linearised problem, de noted by the Orr Sommerfeld and Squire equations. The associated eigenmodes of this linearised problem are the normal velocity and the normal vorticity eigenmodes, which are not dimensionally homogeneous quantities. Thus non-orthogonality condition between these two families of eigenmodes have not been clearly demonstrated yet. Using an orthogonal decomposition of solenoidal velocity fields, a velocity perturbation is expressed as an L2 orthogonal sum of an OrrSommerfeld velocity field (function of the perturbation normal velocity) and a Squirevelocity field (function of the perturbation normal vorticity). Using this decomposition,a variational formulation of the linearised problem is written, that is equivalent to the Orr Sommerfeld and Squire equations, but whose eigenmodes consist of two families of velocity eigenmodes (thus dimensionally homogeneous). We demonstrate that these two sets are non-orthogonal and linear combination between them can produce large transient growth. Using this new formulation, the link between optimal mode and continuous mode transition will also be clari fied, as the role of direct resonance. Numerical solutions are presented to illustrate the analysis in the case of thin boundary layers developing between two parallel walls at large Reynolds number. Characterisations of the destabilizing perturbations will be given in that case.
△ Less
Submitted 21 May, 2013;
originally announced May 2013.