Direct numerical simulation of turbulence modulation by particles in isotropic turbulence

1998 ◽  
Vol 375 ◽  
pp. 235-263 ◽  
Author(s):  
MARC BOIVIN ◽  
OLIVIER SIMONIN ◽  
KYLE D. SQUIRES
Keyword(s):  
Numerical Simulation ◽  
Kinetic Energy ◽  
Direct Numerical Simulation ◽  
Isotropic Turbulence ◽  
Stokes Equations ◽  
Relaxation Times ◽  
Relative Increase ◽  
Fluid Particle ◽  
Small Scale ◽  
Small Particles

The modulation of isotropic turbulence by particles has been investigated using direct numerical simulation (DNS). The particular focus of the present work is on the class of dilute flows in which particle volume fractions and inter-particle collisions are negligible. Gravitational settling is also neglected and particle motion is assumed to be governed by drag with particle relaxation times ranging from the Kolmogorov scale to the Eulerian time scale of the turbulence and particle mass loadings up to 1. The velocity field was made statistically stationary by forcing the low wavenumbers of the flow. The calculations were performed using 963 collocation points and the Taylor-scale Reynolds number for the stationary flow was 62. The effect of particles on the turbulence was included in the Navier–Stokes equations using the point-force approximation in which 963 particles were used in the calculations. DNS results show that particles increasingly dissipate fluid kinetic energy with increased loading, with the reduction in kinetic energy being relatively independent of the particle relaxation time. Viscous dissipation in the fluid decreases with increased loading and is larger for particles with smaller relaxation times. Fluid energy spectra show that there is a non-uniform distortion of the turbulence with a relative increase in small-scale energy. The non-uniform distortion significantly affects the transport of the dissipation rate, with the production and destruction of dissipation exhibiting completely different behaviours. The spectrum of the fluid–particle energy exchange rate shows that the fluid drags particles at low wavenumbers while the converse is true at high wavenumbers for small particles. A spectral analysis shows that the increase of the high-wavenumber portion of the fluid energy spectrum can be attributed to transfer of the fluid–particle covariance by the fluid turbulence. This in turn explains the relative increase of small-scale energy caused by small particles observed in the present simulations as well as those of Squires & Eaton (1990) and Elghobashi & Truesdell (1993).

scholarly journals Anisotropic Characteristics of Turbulence Dissipation in Swirling Flow: A Direct Numerical Simulation Study

2015 ◽  
Vol 2015 ◽  
pp. 1-9 ◽  
Author(s):  
Xingtuan Yang ◽  
Nan Gui ◽  
Gongnan Xie ◽  
Jie Yan ◽  
Jiyuan Tu ◽  
...  
Keyword(s):  
Numerical Simulation ◽  
Energy Dissipation ◽  
Direct Numerical Simulation ◽  
Large Scale ◽  
Isotropic Turbulence ◽  
Vortex Breakdown ◽  
Swirling Flow ◽  
Turbulent Energy ◽  
Swirling Flows ◽  
Small Scale

This study investigates the anisotropic characteristics of turbulent energy dissipation rate in a rotating jet flow via direct numerical simulation. The turbulent energy dissipation tensor, including its eigenvalues in the swirling flows with different rotating velocities, is analyzed to investigate the anisotropic characteristics of turbulence and dissipation. In addition, the probability density function of the eigenvalues of turbulence dissipation tensor is presented. The isotropic subrange of PDF always exists in swirling flows relevant to small-scale vortex structure. Thus, with remarkable large-scale vortex breakdown, the isotropic subrange of PDF is reduced in strongly swirling flows, and anisotropic energy dissipation is proven to exist in the core region of the vortex breakdown. More specifically, strong anisotropic turbulence dissipation occurs concentratively in the vortex breakdown region, whereas nearly isotropic turbulence dissipation occurs dispersively in the peripheral region of the strong swirling flows.

scholarly journals Direct Numerical Simulation and Visualization of Biswirling Jets

2014 ◽  
Vol 6 ◽  
pp. 193731 ◽  
Author(s):  
Jie Yan ◽  
Nan Gui ◽  
Gongnan Xie ◽  
Jinsen Gao
Keyword(s):  
Numerical Simulation ◽  
Kinetic Energy ◽  
Direct Numerical Simulation ◽  
Turbulent Kinetic Energy ◽  
Vortex Structures ◽  
Small Scale ◽  
Convergence Region ◽  
Swirling Jets ◽  
Vortex Interactions ◽  
Close Distance

Two parallel swirling/rotating jets with a distance between them are termed biswirling jets here, which have important and complicated vortex structures different from the single swirling jet due to the negligible vortex-vortex interactions. The visualization of vortex-vortex interaction between the biswirling jets is accomplished by using direct numerical simulation. The evolution of vortex structures of the biswirling jets is found rather complicated. The turbulent kinetic energy and turbulence dissipation in the central convergence region are augmented locally and rather strongly. The modulation of turbulent kinetic energy by jet-jet interaction upon different scales of vortices is dominated by the swirling levels and the distance between the jets. The turbulent kinetic energy upon intermediate and small scale vortices in bijets with not very high swirling level and at a very close distance is smaller than that in single swirling jets, whereas the opposite is true under a far distance, and so forth.

Direct numerical simulation of inertio-elastic turbulent Taylor–Couette flow

2021 ◽  
Vol 926 ◽  
Author(s):  
Jiaxing Song ◽  
Fenghui Lin ◽  
Nansheng Liu ◽  
Xi-Yun Lu ◽  
Bamin Khomami
Keyword(s):  
Numerical Simulation ◽  
Kinetic Energy ◽  
Direct Numerical Simulation ◽  
Couette Flow ◽  
Flow Structures ◽  
Small Scale ◽  
Wall Region ◽  
Flow Physics ◽  
Taylor Couette ◽  
Near Wall

The flow physics of inertio-elastic turbulent Taylor–Couette flow for a radius ratio of $0.5$ in the Reynolds number ( $Re$ ) range of $500$ to $8000$ is investigated via direct numerical simulation. It is shown that as $Re$ is increased the turbulence dynamics can be subdivided into two distinct regimes: (i) a low $Re \leqslant 1000$ regime where the flow physics is essentially dominated by nonlinear elastic forces and the main contribution to transport and mixing of momentum, stress and energy comes from large-scale flow structures in the bulk region and (ii) a high $Re \geqslant 5000$ regime where inertial forces govern the flow physics and the flow dynamics is mainly governed by small-scale flow structures in the near-wall region. Flow–microstructure coupling analysis reveals that the elastic Görtler instability in the near-wall region is triggered via significant polymer extension and commensurately high hoop stresses. This instability gives rise to small-scale elastic vortical structures identified as elastic Görtler vortices which are present at all $Re$ considered. In fact, these vortices develop herringbone streaks near the inner wall that have a longer average life span than their Newtonian counterparts due to their elastic origin. Examination of the budgets of mean streamwise enstrophy, mean kinetic energy, turbulent kinetic energy and Reynolds shear stress demonstrates that increasing fluid inertia hinders the generation of elastic stresses, leading to a monotonic reduction of the elastic-related effects on the flow physics.

Direct numerical simulation of a turbulent jet impinging on a heated wall

2015 ◽  
Vol 764 ◽  
pp. 362-394 ◽  
Author(s):  
T. Dairay ◽  
V. Fortuné ◽  
E. Lamballais ◽  
L.-E. Brizzi
Keyword(s):  
Heat Transfer ◽  
Numerical Simulation ◽  
Nusselt Number ◽  
Direct Numerical Simulation ◽  
Temperature Fields ◽  
Heated Wall ◽  
Small Scale ◽  
Turbulent Regime ◽  
High Heat Transfer ◽  
Radial Evolution

AbstractDirect numerical simulation (DNS) of an impinging jet flow with a nozzle-to-plate distance of two jet diameters and a Reynolds number of 10 000 is carried out at high spatial resolution using high-order numerical methods. The flow configuration is designed to enable the development of a fully turbulent regime with the appearance of a well-marked secondary maximum in the radial distribution of the mean heat transfer. The velocity and temperature statistics are validated with documented experiments. The DNS database is then analysed focusing on the role of unsteady processes to explain the spatial distribution of the heat transfer coefficient at the wall. A phenomenological scenario is proposed on the basis of instantaneous flow visualisations in order to explain the non-monotonic radial evolution of the Nusselt number in the stagnation region. This scenario is then assessed by analysing the wall temperature and the wall shear stress distributions and also through the use of conditional averaging of velocity and temperature fields. On one hand, the heat transfer is primarily driven by the large-scale toroidal primary and secondary vortices emitted periodically. On the other hand, these vortices are subjected to azimuthal distortions associated with the production of radially elongated structures at small scale. These distortions are responsible for the appearance of very high heat transfer zones organised as cold fluid spots on the heated wall. These cold spots are shaped by the radial structures through a filament propagation of the heat transfer. The analysis of probability density functions shows that these strong events are highly intermittent in time and space while contributing essentially to the secondary peak observed in the radial evolution of the Nusselt number.

On Direct Numerical Simulation of Turbulent Flows

2011 ◽  
Vol 64 (2) ◽  
Author(s):  
Giancarlo Alfonsi
Keyword(s):  
Numerical Simulation ◽  
Direct Numerical Simulation ◽  
High Performance Computing ◽  
Turbulent Flows ◽  
High Performance ◽  
Stokes Equations ◽  
Navier Stokes ◽  
Turbulent Shear ◽  
Navier Stokes Equations ◽  
Performance Computing

The direct numerical simulation of turbulence (DNS) has become a method of outmost importance for the investigation of turbulence physics, and its relevance is constantly growing due to the increasing popularity of high-performance-computing techniques. In the present work, the DNS approach is discussed mainly with regard to turbulent shear flows of incompressible fluids with constant properties. A body of literature is reviewed, dealing with the numerical integration of the Navier-Stokes equations, results obtained from the simulations, and appropriate use of the numerical databases for a better understanding of turbulence physics. Overall, it appears that high-performance computing is the only way to advance in turbulence research through the front of the direct numerical simulation.

scholarly journals Flow topologies in bubble-induced turbulence: a direct numerical simulation analysis

2018 ◽  
Vol 857 ◽  
pp. 270-290 ◽  
Author(s):  
Josef Hasslberger ◽  
Markus Klein ◽  
Nilanjan Chakraborty
Keyword(s):  
Numerical Simulation ◽  
Direct Numerical Simulation ◽  
Turbulent Flows ◽  
Volume Fraction ◽  
Fluid Velocity ◽  
Small Scale ◽  
Flow Topology ◽  
Two Phase ◽  
Local Flow ◽  
Interface Curvature

This paper presents a detailed investigation of flow topologies in bubble-induced two-phase turbulence. Two freely moving and deforming air bubbles that have been suspended in liquid water under counterflow conditions have been considered for this analysis. The direct numerical simulation data considered here are based on the one-fluid formulation of the two-phase flow governing equations. To study the development of coherent structures, a local flow topology analysis is performed. Using the invariants of the velocity gradient tensor, all possible small-scale flow structures can be categorized into two nodal and two focal topologies for incompressible turbulent flows. The volume fraction of focal topologies in the gaseous phase is consistently higher than in the surrounding liquid phase. This observation has been argued to be linked to a strong vorticity production at the regions of simultaneous high fluid velocity and high interface curvature. Depending on the regime (steady/laminar or unsteady/turbulent), additional effects related to the density and viscosity jump at the interface influence the behaviour. The analysis also points to a specific term of the vorticity transport equation as being responsible for the induction of vortical motion at the interface. Besides the known mechanisms, this term, related to surface tension and gradients of interface curvature, represents another potential source of turbulence production that lends itself to further investigation.

scholarly journals Numerical stability analysis of a vortex ring with swirl

2019 ◽  
Vol 878 ◽  
pp. 5-36 ◽  
Author(s):  
Yuji Hattori ◽  
Francisco J. Blanco-Rodríguez ◽  
Stéphane Le Dizès
Keyword(s):  
Numerical Simulation ◽  
Growth Rate ◽  
Direct Numerical Simulation ◽  
Vortex Ring ◽  
Stokes Equations ◽  
Base Flow ◽  
Critical Layer ◽  
Navier Stokes ◽  
Instability Mode ◽  
Navier Stokes Equations

The linear instability of a vortex ring with swirl with Gaussian distributions of azimuthal vorticity and velocity in its core is studied by direct numerical simulation. The numerical study is carried out in two steps: first, an axisymmetric simulation of the Navier–Stokes equations is performed to obtain the quasi-steady state that forms a base flow; then, the equations are linearized around this base flow and integrated for a sufficiently long time to obtain the characteristics of the most unstable mode. It is shown that the vortex rings are subjected to curvature instability as predicted analytically by Blanco-Rodríguez & Le Dizès (J. Fluid Mech., vol. 814, 2017, pp. 397–415). Both the structure and the growth rate of the unstable modes obtained numerically are in good agreement with the analytical results. However, a small overestimation (e.g. 22 % for a curvature instability mode) by the theory of the numerical growth rate is found for some instability modes. This is most likely due to evaluation of the critical layer damping which is performed for the waves on axisymmetric line vortices in the analysis. The actual position of the critical layer is affected by deformation of the core due to the curvature effect; as a result, the damping rate changes since it is sensitive to the position of the critical layer. Competition between the curvature and elliptic instabilities is also investigated. Without swirl, only the elliptic instability is observed in agreement with previous numerical and experimental results. In the presence of swirl, sharp bands of both curvature and elliptic instabilities are obtained for $\unicode[STIX]{x1D700}=a/R=0.1$, where $a$ is the vortex core radius and $R$ the ring radius, while the elliptic instability dominates for $\unicode[STIX]{x1D700}=0.18$. New types of instability mode are also obtained: a special curvature mode composed of three waves is observed and spiral modes that do not seem to be related to any wave resonance. The curvature instability is also confirmed by direct numerical simulation of the full Navier–Stokes equations. Weakly nonlinear saturation and subsequent decay of the curvature instability are also observed.

scholarly journals Flow topologies in primary atomization of liquid jets: a direct numerical simulation analysis

2018 ◽  
Vol 859 ◽  
pp. 819-838 ◽  
Author(s):  
Josef Hasslberger ◽  
Sebastian Ketterl ◽  
Markus Klein ◽  
Nilanjan Chakraborty
Keyword(s):  
Numerical Simulation ◽  
Reynolds Number ◽  
Direct Numerical Simulation ◽  
Incompressible Flows ◽  
Flow Structures ◽  
Liquid Jets ◽  
Small Scale ◽  
Flow Topology ◽  
Local Flow ◽  
Turbulent Length Scale

The local flow topology analysis of the primary atomization of liquid jets has been conducted using the invariants of the velocity-gradient tensor. All possible small-scale flow structures are categorized into two focal and two nodal topologies for incompressible flows in both liquid and gaseous phases. The underlying direct numerical simulation database was generated by the one-fluid formulation of the two-phase flow governing equations including a high-fidelity volume-of-fluid method for accurate interface propagation. The ratio of liquid-to-gas fluid properties corresponds to a diesel jet exhausting into air. Variation of the inflow-based Reynolds number as well as Weber number showed that both these non-dimensional numbers play a pivotal role in determining the nature of the jet break-up, but the flow topology behaviour appears to be dominated by the Reynolds number. Furthermore, the flow dynamics in the gaseous phase is generally less homogeneous than in the liquid phase because some flow regions resemble a laminar-to-turbulent transition state rather than fully developed turbulence. Two theoretical models are proposed to estimate the topology volume fractions and to describe the size distribution of the flow structures, respectively. In the latter case, a simple power law seems to be a reasonable approximation of the measured topology spectrum. According to that observation, only the integral turbulent length scale would be required as an input for the a priori prediction of the topology size spectrum.

Sign in / Sign up

Export Citation Format

Share Document