FG Numerische Strömungs- und Gasdynamik
Refine
Year of publication
Document Type
- Scientific journal article peer-reviewed (58)
- Conference publication not peer-reviewed (55)
- Conference Proceeding (53)
- Conference publication peer-reviewed (32)
- Image (poster) (16)
- Scientific journal article not peer-reviewed (4)
- Part of a book (chapter) (4)
- Doctoral thesis (2)
- Report (1)
- Research data (1)
Way of publication
- Open Access (24)
Keywords
- one-dimensional turbulence (49)
- stochastic modeling (26)
- stochastic turbulence modeling (12)
- turbulence modeling (8)
- turbulent mixing (8)
- passive scalar (7)
- turbulent boundary layer (7)
- heat transfer (6)
- atmospheric boundary layer (5)
- high Rayleigh number (5)
Institute
- FG Numerische Strömungs- und Gasdynamik (226)
- Energie-Innovationszentrum / Scientific Computing Lab (12)
- FG Mechanische Verfahrenstechnik (12)
- FG Technische Akustik (5)
- FG Numerische Mathematik und Wissenschaftliches Rechnen (3)
- FG Flug-Triebwerksdesign (2)
- FG Prozess- und Anlagentechnik (2)
- FG Aerodynamik und Strömungslehre (1)
- FG Atmosphärische Prozesse (1)
- FG Hydrologie (1)
Turbulent concentric coaxial pipe flows are numerically investigated as canonical problem addressing spanwise curvature effects on heat and momentum transfer that are encountered in various engineering applications. It is demonstrated that the wall-adapting local eddy-viscosity (WALE) model within a large-eddy simulation (LES) framework, without model parameter recalibration, has limited predictive capabilities as signalized by poor representation of wall curvature effects and notable grid dependence. The identified lack in the modeling of radial transport processes is therefore addressed here by utilizing a stochastic one-dimensional turbulence (ODT) model. A standalone ODT formulation for cylindrical geometry is used in order to assess to which extent the predictability can be expected to improve by utilizing an advanced wall-modeling strategy.
Turbulent mixing is an omnipresent phenomenon that permanently affects our everyday life. Mixing processes also plays an important role in many industrial applications. The full resolution of all relevant flow scales often poses a major challenge to the numerical simulation and requires a modeling of the small-scale effects. In transported Probability Density Function (PDF) methods, the simplified modeling of the molecular mixing is a known weak point. At this place, the Hierarchical Parcel-Swapping (HiPS) model developed by A.R. Kerstein [J. Stat. Phys. 153, 142-161 (2013)] represents a computationally efficient and novel turbulent mixing model. HiPS simulates the effects of turbulence on time-evolving, 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 existing 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 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, a standalone HiPS model formulation for the simulation of passive scalar mixing is detailed first. The
generated scalar power spectra with forced turbulence shows the known scaling law of Kolmogorov turbulence. Furthermore, results for the PDF of the passive scalar, mean square displacement and scalar dissipation rate are shown and reveal a reasonable agreement with experimental findings. The described possibility to account for variable Schmidt number effects is an important next development step for the HiPS formulation. This enables the incorporation of differential diffusion, which represents an immense advantage compared to the established mixing models. Using a binary structure allows HiPS to satisfy a large number of criteria for a good mixing model. Considering the reduced order and associated computational efficiency, HiPS is an attractive
mixing model, which can contribute to an improved representation of the molecular mixing in transported PDF methods.
Modeling simultaneous momentum and passive scalar transfer in turbulent annular Poiseuille flow
(2023)
Simultaneous momentum and passive scalar transfer in weakly heated pressure-driven turbulent concentric annular pipe flow is numerically investigated using the cylindrical formulation of the stochastic one-dimensional turbulence (ODT) model,which is utilized here as standalone tool. In the present study, we focus on the model calibration for heated annular pipes based on recent reference direct numerical simulations (DNS) from Bagheri and Wang (Int. J. Heat Fluid Flow 86, 108725,2020; Phys. Fluids 33, 055131, 2021). It is shown that the model is able to individually capture scalar and momentum transfer, but not both equally well at the same time. We attribute this to less dissimilar scalar and momentum transport in the model at the low Reynolds number investigated. It is argued that the model prefers a fully developed turbulent state due to its construction. Nevertheless, it is demonstrated that ODT is able to reasonably capture the radial inner-outer asymmetry of the scalar and momentum boundary layers which yields better predictive capabilities than wall-function-based approaches.
Accurate but economical modeling of supersonic turbulent boundary layers is a standing challenge due to the intricate entanglement of temperature, density, and velocity fluctuations on top of the mean-field variation. Application of the van Driest transformation may describe well the mean state but cannot provide detailed flow information. This lack-in modeling coarse and fine-scale variability is addressed by the present study using a stochastic one-dimensional turbulence (ODT) model. ODT is a simulation methodology that represents the evolution of turbulent flow in a low-dimensional stochastic way. In this study, ODT is extended to fully compressible flows. An Eulerian framework and a conservative form of the governing equations serve as the basis of the compressible ODT model. Computational methods for statistical properties based on ODT realizations are also extended to compressible flows, and a comprehensive way of turbulent kinetic energy budget calculation based on compressible ODT is put forward for the first time. Two canonical direct numerical simulation cases of supersonic isothermal-wall channel flow at Mach numbers 1.5 and 3.0 with bulk Reynolds numbers 3000 and 4880, respectively, are used to validate the extended model. A rigorous numerical validation is presented, including the first-order mean statistics, the second-order root mean square statistics, and higher-order turbulent fluctuation statistics. In ODT results, both mean and root mean square profiles are accurately captured in the near-wall region. Near-wall temperature spectra reveal that temperature fluctuations are amplified at all turbulent scales as the effects of compressibility increase. This phenomenon is caused by intensified viscous heating at a higher Mach number, which is indicated by the steeper profiles of viscous turbulent kinetic energy budget terms in the very near-wall region. The low computational cost and predictive capabilities of ODT suggest that it is a promising approach for detailed modeling of highly turbulent compressible boundary layers. Furthermore, it is found that the ODT model requires a Mach-number-dependent increase in a viscous penalty parameter Z in wall-bounded turbulent flows to enable accurate capture of the buffer layer.
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.
Investigation of turbulent mixing using a stochastic hierarchical parcel swapping mixing model
(2023)
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.