Numerical Modelling of Hydrodynamic Instabilities in Supercritical Fluids

Author(s):  
Sakir Amiroudine

The case of a supercritical fluid heated from below (Rayleigh-Bénard) in a rectangular cavity is first presented. The stability of the two boundary layers (hot and cold) is analyzed by numerically solving the Navier-Stokes equations with a van der Waals gas and stability diagrams are derived. The very large compressibility and the very low heat diffusivity of near critical pure fluids induce very large density gradients which lead to a Rayleigh–Taylor-like gravitational instability of the heat diffusion layer and results in terms of growth rates and wave numbers are presented. Depending on the relative direction of the interface or the boundary layer with respect to vibration, vibrational forces can destabilize a thermal boundary layer, resulting in parametric/Rayleigh vibrational instabilities. This has recently been achieved by using a numerical model which does not require any equation of state and directly calculates properties from NIST data base, for instance.

Author(s):  
Sakir Amiroudine

The case of a supercritical fluid heated from below (Rayleigh-Bénard) in a rectangular cavity is first presented. The stability of the two boundary layers (hot and cold) is analyzed by numerically solving the Navier-Stokes equations with a van der Waals gas and stability diagrams are derived. The very large compressibility and the very low heat diffusivity of near critical pure fluids induce very large density gradients which lead to a Rayleigh–Taylor-like gravitational instability of the heat diffusion layer and results in terms of growth rates and wave numbers are presented. Depending on the relative direction of the interface or the boundary layer with respect to vibration, vibrational forces can destabilize a thermal boundary layer, resulting in parametric / Rayleigh vibrational instabilities. This has recently been achieved by using a numerical model which does not require any equation of state and directly calculates properties from NIST data base (NIST, 2000) for instance.


Analysis ◽  
2015 ◽  
Vol 35 (4) ◽  
Author(s):  
Emmanuel Grenier ◽  
Yan Guo ◽  
Toan T. Nguyen

AbstractIn this paper we show how the stability of Prandtl boundary layers is linked to the stability of shear flows in the incompressible Navier–Stokes equations. We then recall classical physical instability results, and give a short educational presentation of the construction of unstable modes for Orr–Sommerfeld equations. We end the paper with a conjecture concerning the validity of Prandtl boundary layer asymptotic expansions.


1985 ◽  
Vol 40 (8) ◽  
pp. 789-799 ◽  
Author(s):  
A. F. Borghesani

The Navier-Stokes equations for the fluid motion induced by a disk rotating inside a cylindrical cavity have been integrated for several values of the boundary layer thickness d. The equivalence of such a device to a rotating disk immersed in an infinite medium has been shown in the limit as d → 0. From that solution and taking into account edge effect corrections an equation for the viscous torque acting on the disk has been derived, which depends only on d. Moreover, these results justify the use of a rotating disk to perform accurate viscosity measurements.


1998 ◽  
Vol 371 ◽  
pp. 207-232 ◽  
Author(s):  
G. VITTORI ◽  
R. VERZICCO

Numerical simulations of Navier–Stokes equations are performed to study the flow originated by an oscillating pressure gradient close to a wall characterized by small imperfections. The scenario of transition from the laminar to the turbulent regime is investigated and the results are interpreted in the light of existing analytical theories. The ‘disturbed-laminar’ and the ‘intermittently turbulent’ regimes detected experimentally are reproduced by the present simulations. Moreover it is found that imperfections of the wall are of fundamental importance in causing the growth of two-dimensional disturbances which in turn trigger turbulence in the Stokes boundary layer. Finally, in the intermittently turbulent regime, a description is given of the temporal development of turbulence characteristics.


2013 ◽  
Vol 136 (2) ◽  
Author(s):  
Graham Ashcroft ◽  
Christian Frey ◽  
Kathrin Heitkamp ◽  
Christian Weckmüller

This is the first part of a series of two papers on unsteady computational fluid dynamics (CFD) methods for the numerical simulation of aerodynamic noise generation and propagation. In this part, the stability, accuracy, and efficiency of implicit Runge–Kutta schemes for the temporal integration of the compressible Navier–Stokes equations are investigated in the context of a CFD code for turbomachinery applications. Using two model academic problems, the properties of two explicit first stage, singly diagonally implicit Runge–Kutta (ESDIRK) schemes of second- and third-order accuracy are quantified and compared with more conventional second-order multistep methods. Finally, to assess the ESDIRK schemes in the context of an industrially relevant configuration, the schemes are applied to predict the tonal noise generation and transmission in a modern high bypass ratio fan stage and comparisons with the corresponding experimental data are provided.


1976 ◽  
Vol 78 (2) ◽  
pp. 355-383 ◽  
Author(s):  
H. Fasel

The stability of incompressible boundary-layer flows on a semi-infinite flat plate and the growth of disturbances in such flows are investigated by numerical integration of the complete Navier–;Stokes equations for laminar two-dimensional flows. Forced time-dependent disturbances are introduced into the flow field and the reaction of the flow to such disturbances is studied by directly solving the Navier–Stokes equations using a finite-difference method. An implicit finitedifference scheme was developed for the calculation of the extremely unsteady flow fields which arose from the forced time-dependent disturbances. The problem of the numerical stability of the method called for special attention in order to avoid possible distortions of the results caused by the interaction of unstable numerical oscillations with physically meaningful perturbations. A demonstration of the suitability of the numerical method for the investigation of stability and the initial growth of disturbances is presented for small periodic perturbations. For this particular case the numerical results can be compared with linear stability theory and experimental measurements. In this paper a number of numerical calculations for small periodic disturbances are discussed in detail. The results are generally in fairly close agreement with linear stability theory or experimental measurements.


2021 ◽  
Vol 11 (04) ◽  
pp. 1-11
Author(s):  
Wanwan Li

In mechanical engineering educations, simulating fluid thermodynamics is rather helpful for students to understand the fluid’s natural behaviors. However, rendering both high-quality and realtime simulations for fluid dynamics are rather challenging tasks due to their intensive computations. So, in order to speed up the simulations, we have taken advantage of GPU acceleration techniques to simulate interactive fluid thermodynamics in real-time. In this paper, we present an elegant, basic, but practical OpenGL/SL framework for fluid simulation with a heat map rendering. By solving Navier-Stokes equations coupled with the heat diffusion equation, we validate our framework through some real-case studies of the smoke-like fluid rendering such as their interactions with moving obstacles and their heat diffusion effects. As shown in Fig. 1, a group of experimental results demonstrates that our GPU-accelerated solver of Navier-Stokes equations with heat transfer could give the observers impressive real-time and realistic rendering results.


2013 ◽  
Vol 2013 ◽  
pp. 1-6
Author(s):  
Wen-Juan Wang ◽  
Yan Jia

We study the stability issue of the generalized 3D Navier-Stokes equations. It is shown that if the weak solutionuof the Navier-Stokes equations lies in the regular class∇u∈Lp(0,∞;Bq,∞0(ℝ3)),(2α/p)+(3/q)=2α,2<q<∞,0<α<1, then every weak solutionv(x,t)of the perturbed system converges asymptotically tou(x,t)asvt-utL2→0,t→∞.


Author(s):  
Kazuomi Yamamoto ◽  
Yoshimichi Tanida

A self-excited oscillation of transonic flow in a simplified cascade model was investigated experimentally, theoretically and numerically. The measurements of the shock wave and wake motions, and unsteady static pressure field predict a closed loop mechanism, in which the pressure disturbance, that is generated by the oscillation of boundary layer separation, propagates upstream in the main flow and forces the shock wave to oscillate, and then the shock oscillation disturbs the boundary layer separation again. A one-dimensional analysis confirms that the self-excited oscillation occurs in the proposed mechanism. Finally, a numerical simulation of the Navier-Stokes equations reveals the unsteady flow structure of the reversed flow region around the trailing edge, which induces the large flow separation to bring about the anti-phase oscillation.


2014 ◽  
Vol 752 ◽  
pp. 602-625 ◽  
Author(s):  
Kengo Deguchi ◽  
Philip Hall

AbstractOur concern in this paper is with high-Reynolds-number nonlinear equilibrium solutions of the Navier–Stokes equations for boundary-layer flows. Here we consider the asymptotic suction boundary layer (ASBL) which we take as a prototype parallel boundary layer. Solutions of the equations of motion are obtained using a homotopy continuation from two known types of solutions for plane Couette flow. At high Reynolds numbers, it is shown that the first type of solution takes the form of a vortex–wave interaction (VWI) state, see Hall & Smith (J. Fluid Mech., vol. 227, 1991, pp. 641–666), and is located in the main part of the boundary layer. On the other hand, here the second type is found to support an equilibrium solution of the unit-Reynolds-number Navier–Stokes equations in a layer located a distance of $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}O(\ln \mathit{Re})$ from the wall. Here $\mathit{Re}$ is the Reynolds number based on the free-stream speed and the unperturbed boundary-layer thickness. The streaky field produced by the interaction grows exponentially below the layer and takes its maximum size within the unperturbed boundary layer. The results suggest the possibility of two distinct types of streaky coherent structures existing, possibly simultaneously, in disturbed boundary layers.


Sign in / Sign up

Export Citation Format

Share Document