An improved parallel compact scheme for domain‐decoupled simulation of turbulence

Fang, Jiaojiao; Gao, Feng; Moulinec, Charles; Emerson, David R.

doi:10.1002/fld.4731

Cited by 24 publications

(16 citation statements)

References 74 publications

(146 reference statements)

Supporting

Mentioning

Contrasting

Order By: Relevance

“…This approach has been recently applied to studies of various SWTBLI problems (see Fang et al 2014Fang et al , 2015Fang et al , 2017. The diffusion terms are solved using a sixth-order compact central scheme (Hirsh 1975;Lele 1992) with a domain decoupling scheme for parallel computation (Fang et al 2019). After all the spatial terms are solved, a three-step third-order total variation diminishing Runge-Kutta method, proposed by Gottlieb & Shu (1998), is used for the temporal integration.…”

Section: Computational Set-upmentioning

confidence: 99%

On the turbulence amplification in shock-wave/turbulent boundary layer interaction

et al. 2020

View full text Add to dashboard Cite

Section: Computational Set-upmentioning

confidence: 99%

On the turbulence amplification in shock-wave/turbulent boundary layer interaction

et al. 2020

View full text Add to dashboard Cite

“…All the spatial derivatives are approximated with a sixth-order pade-type compact central scheme [23] . The compact scheme is decoupled at the subdomain interface by using a recently proposed method [24] to handle the domain-decomposition based parallel computation. To remove the small-scale wiggles due to aliasing errors resulting from discrete evaluation of the non-linear convective terms, a tenth-order compact filter is adopted, which limits the filter only at high wavenumbers [25] .…”

Section: Methodsmentioning

confidence: 99%

Non-equilibrium Turbulent Phenomena in a Transitional Flat Plate Boundary Layer

Liu¹,

Fang²,

Lu³

et al. 2019

Proceedings of Global Power &Amp; Propulsion Society

View full text Add to dashboard Cite

Many recent laboratory experiments and numerical simulations supported a non-equilibrium dissipation scaling in decaying turbulence before it reaches an equilibrium state. However, these studies were mostly focus on canonical flows, such as homogeneous isotropic turbulence and gridgenerated turbulence. Beyond that, is it possible to observe the non-equilibrium phenomena in wall turbulence? In the present contribution, by analyzing a direct numerical simulation (DNS) database of a transitional boundary-layer flow, we show that the transition region and the nonequilibrium turbulence region, which are divided in the streamwise direction, presents different non-equilibrium scalings. Moreover, in the wall-normal direction, buffer layer, log layer and outer layer show different nonequilibrium phenomena which are very different from that in transitional channel flows. Thus, the turbulence models which consider the non-equilibrium phenomena should be modelled according to the dissipation scalings both in the streamwise direction and in the wall-normal direction. These findings are expected to shed light on the future modelling of various types of non-equilibrium turbulent flows.

show abstract

“…Finite difference compact schemes exhibit superior spectral resolution [40,42] than the explicit schemes of same order, which motivated researchers to use them in many DNS and large eddy simulations [12,20,35,[43][44][45][46][47][48]. The approximation of m th derivative of a function f at i th point using a compact scheme is represented as,…”

Section: Formulation Of Compressible Navier-stokes Equation For Rtimentioning

confidence: 99%

Ultrasound Triggering of Rayleigh-Taylor Instability: Solution of Compressible Navier-Stokes Equation by a Non-Overlapping Parallel Compact Scheme

Sundaram¹,

Sengupta²,

Sengupta³

2022

Preprint

View full text Add to dashboard Cite

Rayleigh-Taylor instability (RTI) occurs at the interface of two media when the heavier fluid is accelerated into the lighter fluid and is a prototypical hydrodynamic event present in many physical events. In high energy physics, this manifests itself across a wide range of length scales from nuclear confinement fusion at micron-scale to supernova explosion at terra scales. RTI can also be viewed as a baroclinic instability prevalent in engineering, geophysics, and astrophysics, a pedagogic description of which is given in Sengupta et al., Comput. Fluids, 225, 104995 (2021) with respect to the experimental results of Read, Physica D, 12 45-58 (1984). Here, a recently proposed non-overlapping parallel algorithm is used to solve this three-dimensional canonical problem, having the unique property of not distinguishing between sequential and parallel computing, using 4.19 billion points and a refined time step of 7.69 × 10 −8 sec. The problem achieves the required density gradient by considering two volumes of air at different temperatures (with a temperature difference of 200K) separated by a non-conducting, impermeable partition at the onset of the experiment, which is removed impulsively at t = 0. The resulting buoyancy force at the interface acting from top to bottom is the seed of the baroclinic instability. Present high precision computation enables one to capture the ensuing RTI triggered by ultrasonic waves created at the interface.

show abstract

An improved parallel compact scheme for domain‐decoupled simulation of turbulence

Cited by 24 publications

References 74 publications

On the turbulence amplification in shock-wave/turbulent boundary layer interaction

On the turbulence amplification in shock-wave/turbulent boundary layer interaction

Non-equilibrium Turbulent Phenomena in a Transitional Flat Plate Boundary Layer

Ultrasound Triggering of Rayleigh-Taylor Instability: Solution of Compressible Navier-Stokes Equation by a Non-Overlapping Parallel Compact Scheme

Contact Info

Product

Resources

About