Search

Turbulent mixing is an omnipresent phenomenon that constantly affects our everyday life and plays an important role in a variety of industrial applications. The simulation of turbulent mixing poses great challenges, since the full resolution of all relevant length and time scales is associated with an immense computational effort. This limitation can be overcome by only resolving the large‐scale effects and completely model the sub‐grid scales. The development of an accurate sub‐grid mixing model is therefore a key challenge to capture all interactions in the sub‐grid scales. At this place, the hierarchical parcel‐swapping (HiPS) model formulated by A.R. Kerstein [J. Stat. Phys. 153, 142–161 (2013)] represents a computationally efficient and scale‐resolving turbulent mixing model. HiPS mimics the effects of turbulence on time‐evolving, diffusive scalar fields. In HiPS, the diffusive scalar fields or a state space is interpreted as a binary tree structure, which is an alternative approach compared to the most common mixing models. Every level of the tree represents a specific length and time scale, which is based on turbulence inertial range scaling. The state variables are only located at the base of the tree and are treated as fluid parcels. The effects of turbulent advection are represented by stochastic swaps of sub‐trees at rates determined by turbulent time scales associated with the sub‐trees. The mixing only takes places between adjacent fluid parcels and at rates consistent with the prevailing diffusion time scales. In this work, the HiPS model formulation for the simulation of passive scalar mixing is detailed first. Preliminary results for the mean square displacement, passive scalar probability density function (PDF) and scalar dissipation rate are given and reveal the strengths of the HiPS model considering the reduced order and computational efficiency. These model investigations are an important step of further HiPS advancements. The integrated auxiliary binary tree structure allows HiPS to satisfy a large number of criteria for a good mixing model. From this point of view, HiPS is an attractive candidate for modeling the mixing in transported PDF methods. [ABSTRACT FROM AUTHOR]

One-dimensional turbulence (ODT) is a simulation methodology that represents the essential physics of three-dimensional turbulence through stochastic resolution of the full range of length and time scales on a one-dimensional domain. In the present study, full compressible modifications are incorporated into ODT methodology, based on an Eulerian framework and a conservative form of the governing equations. In the deterministic part of this approach, a shock capturing scheme is introduced for the first time. In the stochastic part, one-dimensional eddy events are modeled and sampled according to standard methods for compressible flow simulation. Time advancement adjustments are made to balance comparable time steps between the deterministic and stochastic parts in compressible flows. Canonical shock–turbulence interaction cases involving Richtmyer–Meshkov instability at Mach numbers 1.24, 1.5, and 1.98 are simulated to validate the extended model. The ODT results are compared with available reference data from large eddy simulations and laboratory experiments. The introduction of a shock capturing scheme significantly improves the performance of the ODT method, and the results for turbulent kinetic energy are qualitatively improved compared with those of a previous compressible Lagrangian ODT method [Jozefik et al., "Simulation of shock–turbulence interaction in non-reactive flow and in turbulent deflagration and detonation regimes using one-dimensional turbulence," Combust. Flame 164, 53 (2016)]. For the time evolution of profiles of the turbulent mixing zone width, ensemble-averaged density, and specific heat ratio, the new model also yields good to reasonable results. Furthermore, it is found that the viscous penalty parameter Z of the ODT model is insensitive to compressibility effects in turbulent flows without wall effects. A small value of Z is appropriate for turbulent flows with weak wall effects, and the parameter Z serves to suppress extremely small eddy events that would be dissipated instantly by viscosity. [ABSTRACT FROM AUTHOR]

Turbulent mixing plays an important role in a variety of applications ranging from astrophysics to combustion and even pollutant dispersion. The Direct Numerical Simulation (DNS) that resolves all scales is not feasible for most engineering applications since the flow has a wide range of length and time scales, which yields extremely high resolution requirements. Large Eddy Simulations (LES) overcome this limitation by modeling the sub‐grid scale effects. In transported Probability Density Function (PDF) methods, the key challenge is to develop an accurate mixing model. At this point, the Hierarchical Parcel‐Swapping (HiPS) model, introduced by A.R. Kerstein [J. Stat. Phys. 153, 142‐161 (2013)], is an attractive candidate. It is characterized by a computationally efficient representation of the effects of turbulence on a time‐evolving structure of diffusive scalar fields. The interpretation of the diffusive scalar fields or a state space as a binary tree structure is an alternative approach compared to the most common mixing models. The characteristic feature of HiPS is that every level of the tree corresponds to a specific length and time scale, which is based on turbulence inertial range scaling. The state variables only reside at the base of the tree and are understood as fluid parcels. The effects of turbulent advection are represented by stochastic swaps of sub‐trees at rates determined by prevailing turbulent time scales associated with the sub‐trees. The mixing of adjacent fluid parcels is done either instantaneously or at rates consistent with the corresponding diffusion time scales. In this work, HiPS is detailed for the simulation of passive scalar mixing first. Preliminary results for the scalar power spectra, mean square displacement and scalar dissipation rate are shown and reveal a reasonable agreement with experimental findings. Furthermore, the integrated binary tree structure allows to satisfy a large number of criteria for a good mixing model. Considering the reduced order and computational efficiency, HiPS is an attractive candidate for modeling the mixing in transported PDF methods. [ABSTRACT FROM AUTHOR]

[Display omitted] • Lower-order modeling of transient boundary-layer dynamics in turbulent channel flows. • All relevant time and length scales of the flow are resolved by the model. • State-space statistics of surface scalar and momentum fluxes are captured. • Surface layer obeys inner scaling independent of the scalar forcing used. • Regime-independent prediction of Sc, Pe, and Re dependence of the scalar transfer. Accurate and economical modeling of near-surface transport processes is a standing challenge for various engineering and atmospheric boundary-layer flows. In this paper, we address this challenge by utilizing an economical stochastic one-dimensional turbulence (ODT) model. ODT aims to resolve all relevant scales of a turbulent flow for a one-dimensional domain. Here ODT is applied to turbulent channel flow as stand-alone tool. The ODT domain is a wall-normal line that is aligned with the mean shear. The free model parameters are calibrated once for the turbulent velocity boundary layer at a fixed Reynolds number. After that, we use ODT to investigate the Schmidt (Sc), Reynolds (Re), and Peclet (Pe) number dependence of the scalar boundary-layer structure, turbulent fluctuations, transient surface fluxes, mixing, and transfer to a wall. We demonstrate that the model is able to resolve relevant wall-normal transport processes across the turbulent boundary layer and that it captures state-space statistics of the surface scalar-flux fluctuations. In addition, we show that the predicted mean scalar transfer, which is quantified by the Sherwood (Sh) number, self-consistently reproduces established scaling regimes and asymptotic relations with respect to Sc , Re , and Pe. For high asymptotic Sc and Re , ODT results fall between the Dittus–Boelter, Sh ∼ Re 4 / 5 Sc 2 / 5 , and Colburn, Sh ∼ Re 4 / 5 Sc 1 / 3 , scalings but they are closer to the former. For finite Sc and Re , the model prediction reproduces the relation proposed by Schwertfirm and Manhart (Int. J. Heat Fluid Flow , 28 , 1204–1214, 2007) that is based on boundary-layer theory and yields a locally steeper effective scaling than any of the established asymptotic relations. The model extrapolates the scalar transfer to small asymptotic Sc ≪ Re τ - 1 (diffusive limit) with a functional form that has not been previously described. [ABSTRACT FROM AUTHOR]

Hierarchical Parcel‐Swapping (HiPS) developed by A.R. Kerstein [J. Stat. Phys. 153, 142‐161 (2013)] is a computationally efficient and novel model for the effects of turbulence on time‐evolving, diffusive scalar fields. The characteristic feature of HiPS is the interpretation of the one‐dimensional flow domain or a state space as a binary tree structure. Every tree level corresponds to a specific length and time scale, which is based on a turbulence inertial range scaling. The state variables reside at the base of the tree and are interpreted as fluid parcels. The effects of turbulent advection are represented by stochastic swaps of sub‐trees at rates determined by turbulent time scales associated with the sub‐trees. The mixing of adjacent fluid parcels is done at rates consistent with the prevailing diffusion time scales. In this work, we investigate the influence of turbulent time scale variations on an isothermal series‐parallel reaction scheme. The production of a desired chemical species is evaluated by means of a defined selectivity and is strongly affected by the underlying mixing time scales. [ABSTRACT FROM AUTHOR]

• We introduce an extension of the spatial ODT formulation (S-ODT), in order to study a variable density internal air flow. • S-ODT results are compared to the results of a traditional temporal ODT formulation (T-ODT) and to reference DNS data. • For bulk quantities, both the T-ODT and S-ODT formulations have a similar performance in comparison to DNSs. • For the evaluation of radial distributions, the S-ODT formulation clearly outperforms the T-ODT formulation. A novel spatial formulation of the One-Dimensional Turbulence (ODT) model is applied to a vertical pipe-flow with heat transfer, analogous to the Direct Numerical Simulation (DNS) performed by Bae et al. [Phys. Fluids 18, (075102) (2006)]. The framework presented here is an extension for radially confined domains of the cylindrical ODT spatial formulation for low Mach number flows with variable density. The variable density simulations for air (Prandtl number P r = 0.71) are performed at an initial bulk Reynolds number R e b , 0 , D N S = 6000 and Grashof number G r 0 , D N S = 6.78 × 10 6. ODT results are presented for both the spatial formulation introduced in this work and the standard temporal formulation for cylindrical flows introduced by Lignell et al. [Theor. Comput. Fluid Dyn. 32, 4 (2018), pp. 495–20]. Streamwise bulk profiles and radial profiles at specific streamwise positions for the temporal and spatial formulations are in good agreement with the DNS results from Bae et al. For the present application, the spatial formulation yields physically better results in comparison to the temporal formulation. Overall, the findings in the original work of Bae et al. were corroborated with ODT. Although the framework proposed in this work is not a compressible framework and has some clear limitations regarding conservation properties, we suggest its use for future studies in the low Mach number variable density regime. [ABSTRACT FROM AUTHOR]

Modeling turbulent mixing is a standing challenge for nonpremixed chemically reacting flows. Key complications arise from the requirement to capture all relevant scales of the flow and the necessity to distinguish between turbulent advective transport and molecular diffusive transport processes. In addition, anisotropic mean shear, variable advection time scales, and the coexistence of turbulent and nonturbulent regions need to be represented. The fundamental issues at stake are addressed by investigating multi-scalar mixing in a three-stream coaxial jet with a map-based stochastic one-dimensional turbulence model. ODT provides full-scale resolution at affordable costs by a radical reduction of complexity compared to high-fidelity Navier–Stokes solvers. The approach is partly justified by an application of the boundary-layer approximation, but neglects fluctuating axial pressure gradients. It is demonstrated that low-order scalar statistics are reasonably but not fully captured. Despite this shortcoming, it is shown that the model is able to reproduce experimental state-space statistics of multi-stream multi-scalar mixing. The model therefore offers physics-compatible improvements in multi-stream mixing modeling despite some fundamental limitations that remain from unjustified assumptions. [Display omitted] • Map-based stochastic modeling of multi-stream mixing with full-scale resolution. • Comprehensive assessment of spatial sensitivities of scalar joint probabilities. • State-space statistics of scalar mixing in shear turbulence are reproduced. • Scalar statistics are obtained diagnostically without closure assumptions. • Local advection velocities and turbulent/nonturbulent regions are represented. [ABSTRACT FROM AUTHOR]