scholarly journals Numerical simulations of the Lagrangian averaged Navier–Stokes equations for homogeneous isotropic turbulence

2003 ◽  
Vol 15 (2) ◽  
pp. 524-544 ◽  
Author(s):  
Kamran Mohseni ◽  
Branko Kosović ◽  
Steve Shkoller ◽  
Jerrold E. Marsden
2019 ◽  
Vol 864 ◽  
pp. 244-272 ◽  
Author(s):  
L. Djenidi ◽  
R. A. Antonia ◽  
S. L. Tang

The problem of homogeneous isotropic turbulence (HIT) is revisited within the analytical framework of the Navier–Stokes equations, with a view to assessing rigorously the consequences of the scale invariance (an exact property of the Navier–Stokes equations) for any Reynolds number. The analytical development, which is independent of the 1941 (K41) and 1962 (K62) theories of Kolmogorov for HIT for infinitely large Reynolds number, is applied to the transport equations for the second- and third-order moments of the longitudinal velocity increment, $(\unicode[STIX]{x1D6FF}u)$. Once the normalised equations and the constraints required for complying with the scale-invariance property of the equations are presented, results derived from these equations and constraints are discussed and compared with measurements. It is found that the fluid viscosity, $\unicode[STIX]{x1D708}$, and the mean kinetic energy dissipation rate, $\overline{\unicode[STIX]{x1D716}}$ (the overbar denotes spatial and/or temporal averaging), are the only scaling parameters that make the equations scale-invariant. The analysis further leads to expressions for the distributions of the skewness and the flatness factor of $(\unicode[STIX]{x1D6FF}u)$ and shows that these distributions must exhibit plateaus (of different magnitudes) in the dissipative and inertial ranges, as the Taylor microscale Reynolds number $Re_{\unicode[STIX]{x1D706}}$ increases indefinitely. Also, the skewness and flatness factor of the longitudinal velocity derivative become constant as $Re_{\unicode[STIX]{x1D706}}$ increases; this is supported by experimental data. Further, the analysis, backed up by experimental evidence, shows that, beyond the dissipative range, the behaviour of $\overline{(\unicode[STIX]{x1D6FF}u)^{n}}$ with $n=2$, 3 and 4 cannot be represented by a power law of the form $r^{\unicode[STIX]{x1D701}_{n}}$ when the Reynolds number is finite. It is shown that only when $Re_{\unicode[STIX]{x1D706}}\rightarrow \infty$ can an $n$-thirds law (i.e. $\overline{(\unicode[STIX]{x1D6FF}u)^{n}}\sim r^{\unicode[STIX]{x1D701}_{n}}$, with $\unicode[STIX]{x1D701}_{n}=n/3$) emerge, which is consistent with the onset of a scaling range.


2013 ◽  
Vol 732 ◽  
pp. 316-331 ◽  
Author(s):  
Diego A. Donzis ◽  
John D. Gibbon ◽  
Anupam Gupta ◽  
Robert M. Kerr ◽  
Rahul Pandit ◽  
...  

AbstractThe issue of intermittency in numerical solutions of the 3D Navier–Stokes equations on a periodic box ${[0, L] }^{3} $ is addressed through four sets of numerical simulations that calculate a new set of variables defined by ${D}_{m} (t)= {({ \varpi }_{0}^{- 1} {\Omega }_{m} )}^{{\alpha }_{m} } $ for $1\leq m\leq \infty $ where ${\alpha }_{m} = 2m/ (4m- 3)$ and ${[{\Omega }_{m} (t)] }^{2m} = {L}^{- 3} \int \nolimits _{\mathscr{V}} {\vert \boldsymbol{\omega} \vert }^{2m} \hspace{0.167em} \mathrm{d} V$ with ${\varpi }_{0} = \nu {L}^{- 2} $. All four simulations unexpectedly show that the ${D}_{m} $ are ordered for $m= 1, \ldots , 9$ such that ${D}_{m+ 1} \lt {D}_{m} $. Moreover, the ${D}_{m} $ squeeze together such that ${D}_{m+ 1} / {D}_{m} \nearrow 1$ as $m$ increases. The values of ${D}_{1} $ lie far above the values of the rest of the ${D}_{m} $, giving rise to a suggestion that a depletion of nonlinearity is occurring which could be the cause of Navier–Stokes regularity. The first simulation is of very anisotropic decaying turbulence; the second and third are of decaying isotropic turbulence from random initial conditions and forced isotropic turbulence at fixed Grashof number respectively; the fourth is of very-high-Reynolds-number forced, stationary, isotropic turbulence at up to resolutions of $409{6}^{3} $.


2013 ◽  
Vol 2013 ◽  
pp. 1-4 ◽  
Author(s):  
Nicola de Divitiis

The purpose of this paper is to improve a hypothesis of the previous work of N. de Divitiis (2011) dealing with the finite-scale Lyapunov analysis of isotropic turbulence. There, the analytical expression of the structure function of the longitudinal velocity differenceΔuris derived through a statistical analysis of the Fourier transformed Navier-Stokes equations and by means of considerations regarding the scales of the velocity fluctuations, which arise from the Kolmogorov theory. Due to these latter considerations, this Lyapunov analysis seems to need some of the results of the Kolmogorov theory. This work proposes a more rigorous demonstration which leads to the same structure function, without using the Kolmogorov scale. This proof assumes that pair and triple longitudinal correlations are sufficient to determine the statistics ofΔurand adopts a reasonable canonical decomposition of the velocity difference in terms of proper stochastic variables which are adequate to describe the mechanism of kinetic energy cascade.


2021 ◽  
Vol 930 ◽  
Author(s):  
Kartik P. Iyer ◽  
Katepalli R. Sreenivasan ◽  
P.K. Yeung

Using direct numerical simulations performed on periodic cubes of various sizes, the largest being $8192^3$ , we examine the nonlinear advection term in the Navier–Stokes equations generating fully developed turbulence. We find significant dissipation even in flow regions where nonlinearity is locally absent. With increasing Reynolds number, the Navier–Stokes dynamics amplifies the nonlinearity in a global sense. This nonlinear amplification with increasing Reynolds number renders the vortex stretching mechanism more intermittent, with the global suppression of nonlinearity, reported previously, restricted to low Reynolds numbers. In regions where vortex stretching is absent, the angle and the ratio between the convective vorticity and solenoidal advection in three-dimensional isotropic turbulence are statistically similar to those in the two-dimensional case, despite the fundamental differences between them.


Author(s):  
Djordje Romanic ◽  
Horia Hangan

Analytical and semi-empirical models are inexpensive to run and can complement experimental and numerical simulations for risk analysis-related applications. Some models are developed by employing simplifying assumptions in the Navier-Stokes equations and searching for exact, but many times inviscid solutions occasionally complemented by boundary layer equations to take surface effects into account. Other use simple superposition of generic, canonical flows for which the individual solutions are known. These solutions are then ensembled together by empirical or semi-empirical fitting procedures. Few models address turbulent or fluctuating flow fields, and all models have a series of constants that are fitted against experiments or numerical simulations. This chapter presents the main models used to provide primarily mean flow solutions for tornadoes and downbursts. The models are organized based on the adopted solution techniques, with an emphasis on their assumptions and validity.


Author(s):  
Jaromi´r Hora´cˇek ◽  
Miloslav Feistauer ◽  
Petr Sva´cˇek

The contribution deals with the numerical simulation of the flutter of an airfoil with three degrees of freedom (3-DOF) for rotation around an elastic axis, oscillation in the vertical direction and rotation of a flap. The finite element (FE) solution of two-dimensional (2-D) incompressible Navier-Stokes equations is coupled with a system of nonlinear ordinary differential equations describing the airfoil vibrations with large amplitudes taking into account the nonlinear mass matrix. The time-dependent computational domain and a moving grid are treated by the Arbitrary Lagrangian-Eulerian (ALE) method and a suitable stabilization of the FE discretization is applied. The developed method was successfully tested by the classical flutter computation of the critical flutter velocity using NASTRAN program considering the linear model of vibrations and the double-lattice aerodynamic theory. The method was applied to the numerical simulations of the post flutter regime in time domain showing Limit Cycle Oscillations (LCO) due to nonlinearities of the flow model and vibrations with large amplitudes. Numerical experiments were performed for the airfoil NACA 0012 respecting the effect of the air space between the flap and the main airfoil.


Author(s):  
Patrick Queutey ◽  
Jeroen Wackers ◽  
Alban Leroyer ◽  
GanBo Deng ◽  
Emmanuel Guilmineau ◽  
...  

The paper focuses on the hydrodynamic flow around a ship with pods in waves and compares the results of an experimental campaign with numerical simulations conducted during the EU-funded STREAMLINE project. It was the first project for which the effect of waves on cavitation and ventilation was explored in both experimental and numerical ways for a ship with pods. The measurements were carried out in MARIN’s Depressurized Wave Basin (DWB) with a fully instrumented podded ship model, in sailing condition, in waves and depressurised conditions. In this way, the correct representation of cavitation and possible ventilation bubbles and vortices is ensured, resulting in a correct physical behaviour. The discretisation of the Reynolds-Averaged Navier-Stokes Equations (RANSE) is based on the unstructured finite-volume flow solver ISIS-CFD developed by ECN-CNRS. An essential feature for full RANSE simulations with this code is the use of a sliding grid technique to simulate the real propeller rotating behind a ship hull. The computational study in operational service conditions considered here has been conducted to evaluate the instantaneous flow distribution around the podded propellers and to analyse and to compare the unsteady behaviour of the forces induced by the rotating propeller in waves with the measurements from omnidirectional propeller loads as well as the blade forces and moments. The computational study has been done in model and full scale to evaluate the scale effects.


Sign in / Sign up

Export Citation Format

Share Document