FG Numerische Strömungs- und Gasdynamik
Refine
Year of publication
Document Type
- Scientific journal article peer-reviewed (70)
- Conference publication not peer-reviewed (67)
- Conference Proceeding (53)
- Conference publication peer-reviewed (45)
- Image (poster) (23)
- Part of a book (chapter) (7)
- Scientific journal article not peer-reviewed (4)
- Doctoral thesis (2)
- Report (1)
- Research data (1)
Way of publication
- Open Access (60)
Keywords
- one-dimensional turbulence (55)
- stochastic modeling (31)
- stochastic turbulence modeling (12)
- turbulent mixing (10)
- turbulence modeling (9)
- atmospheric boundary layer (7)
- heat transfer (7)
- passive scalar (7)
- turbulent boundary layer (7)
- high Rayleigh number (5)
Institute
- FG Numerische Strömungs- und Gasdynamik (273)
- Energie-Innovationszentrum / Scientific Computing Lab (40)
- FG Mechanische Verfahrenstechnik (12)
- FG Prozess- und Anlagentechnik (6)
- FG Technische Akustik (5)
- FG Numerische Mathematik und Wissenschaftliches Rechnen (3)
- FG Flug-Triebwerksdesign (2)
- FG Aerodynamik und Strömungslehre (1)
- FG Atmosphärische Prozesse (1)
- FG Hydrologie (1)
This work is focused on modeling the effects of homogeneous roughness on low-order velocity statistics in turbulent channel flows. Hydrodynamic effects due to the roughness are characterized on the basis of volume-averaging theory (VAT) and a discrete roughness element method. This theory exploits the homogeneous character of the roughness in order to reduce the complexity of the flow to its one-dimensional statistics. The formulated VAT-based roughness forcing is best suited for drag dominated surfaces. Turbulence modeling closure is achieved with a map-based turbulence model, the One-Dimensional Turbulence (ODT) model. This avoids the prescription of laws of the wall or other ad-hoc scalings, unlike in more traditional filter-based turbulence models. The modeling framework is applied on selected Reynolds number flows for likewise selected roughness topologies. Results are compared to direct numerical simulation (DNS) data available from the literature. Among others, model results are compared with those of a previously formulated parametric forcing approach (PFA) for roughness drag which involved a costly coefficient calibration linked to the roughness topology model. In ODT, the only calibration process required is the same one involved for the turbulence model parameters, i.e., similar to the ODT model application for smooth-wall flows. Despite all of the inherently implied shortcomings of a 1-D model, some appealing properties of ODT are discussed. Notably, the model is able to predict the roughness function, as well as the wall-normal profile of the Reynolds shear stress across the entire boundary layer thickness.
Turbulent concentric coaxial (annular) pipe flow is numerically investigated using a stochastic one-dimensional turbulence (ODT) model as standalone tool. The dimensionally reduced ODT domain enables fully resolved numerical simulations of the flow across the radial gap between the cylindrical inner wall and the cylindrical outer wall. The model is calibrated with available reference data at low bulk Reynolds number for a wide (radius ratio ) and a moderate () gap. Making use of the model’s predictive capabilities, radius ratio and Reynolds number effects are investigated, reaching bulk Reynolds numbers as large as . Despite the large values reached, spanwise wall-curvature effects remain sensible in the momentum boundary layer. The effects are more pronounced for larger wall curvature and to leading orders restricted to the convex cylindrical inner wall. Wall-curvature corrections to the law of the wall are obtained for both the viscous and Reynolds-stress dominated regions by fitting analytically derived expressions for the flow profile to the stochastic simulation data, demonstrating physical compatibility with Reynolds-averaged Navier–Stokes flow. Second-order and detailed fluctuation statistics demonstrate the permeating and nonlocal influence of spanwise wall curvature on the turbulent boundary layer. Surrogate model output in terms of conditional eddy event statistics reveals that the disparity between the near-inner and near-outer wall turbulence increases with Reynolds number for small radius ratios, suggesting that annular pipe flows require wall-curvature-aware wall models even at very large Reynolds numbers.
This study presents a numerical investigation of passive scalar mixing in homogeneous isotropic turbulence (HIT). Different volumetric forcing schemes have been used in the literature, but the side effects are rarely discussed, either because these are assumed irrelevant or because it is too costly to conduct such an analysis with a high‐fidelity model. In this study, we have used One‐Dimensional Turbulence (ODT) model to compare forcing schemes at low Reynolds numbers. Our analysis reveals critical flaws in the linear forcing model when applied to ODT. While both schemes exhibit spectral deviations from direct numerical simulation (DNS), the stochastic forcing scheme demonstrates superior dynamic fidelity, better capturing the turbulent energy cascade. In contrast, the linear forcing scheme suffers from a non‐physical energy deficit at large scales and is approximately 10 times more computationally expensive. These artefacts directly impact scalar mixing: The stochastic scheme produces classic, multi‐scale intermittency, whereas linear forcing generates extreme gradients confined only at the dissipative scales. These results demonstrate that the choice of forcing is a critical modelling decision in ODT, leading to fundamentally different model‐dependent artifacts in both turbulence dynamics and scalar mixing statistics, at least in low Reynolds number regimes.
Wind turbine predesign is challenged by the representation of site-specific wind conditions. A good deal of that challenge lies in the modeling of the inflow turbulence in the atmospheric boundary layer (ABL). A stochastic one-dimensional turbulence (ODT) model is applied to an idealized neutrally stratified ABL and evolves the instantaneous velocity profile with full-scale resolution. The model is able to reproduce the law of the wall consistently after an initial calibration with the surface drag law. Investigating turbulent time series of the horizontal velocity components it is demonstrated that the model generates physically justified intermittency features with increasing turbulence intensity.
The present study investigates the impact of various surface boundary conditions on turbulent Rayleigh–Bénard convection within a cubic cavity configuration. Simulations are conducted with a characteristic-based off-lattice Boltzmann method (LBM) solver for mildly turbulent flow of water using a direct numerical simulation (DNS) approach. The current study considers different boundary conditions such as no-slip, free-slip, and Navier-slip conditions on the walls with variations in slip length and wall-slip anisotropy. Results are evaluated through mean isotherms, streamlines, root-mean-square fluctuations, and Nusselt number. The results obtained demonstrate that the selection of wall-boundary conditions has a significant influence on the flow organization within the cavity and on the heat transfer across it.
We report on a systematic study for Reynolds–Averaged Navier-Stokes (RANS) modeling and simulations of turbulent annular pipe flow. Several simulations were performed using the most readily-available RANS models in the open-source library OpenFOAM. A customized 1-D RANS solver was also developed for ease of analysis. The focus of the study is on the reproduction of the mean velocity profile, its maximum, and maximum radial location, as well as modeled low-order fluctuation statistics. The flow in the annular gap is characterized by a radius ratio of 0.1, and a friction Reynolds number equal to 600 that is based on a mean friction velocity. Deviations from the mean velocity profile are observed for all RANS models investigated when compared with Direct Numerical Simulation (DNS) reference data. The representation of the near-wall outer cylinder flow is better than that of the near-wall inner cylinder flow.
Heated concentric coaxial (annular) pipe flows are numerically investigated by a stochastic one‐dimensional turbulence (ODT) model. The main objective of this work is to more accurately predict the heat transfer in tubular heat exchangers at low Prandtl numbers by extending the analysis for weak temperature fluctuations as recently introduced in Tsai et al., Proceedings in Applied Mathematics and Mechanics, 23:e202300167, 2023. The ODT model offers the required predictive capabilities at affordable cost by providing full‐scale resolution of viscous, conductive, and turbulent advective transport processes along a representative radial coordinate. The Prandtl numbers and are considered for which the radius ratio of the annular pipe and the Reynolds number are varied. Numerical results demonstrate that the geometry (radius ratio) has a significant influence on the thermal boundary layer that emerges over the inner and outer curved wall, respectively. Application of boundary layer theory and mixing length arguments yield an analytical expression that includes both Reynolds number and curvature effects. Unknown closure coefficients are estimated with ODT, providing a physically based correlation for the Nusselt number.
In this work, we discuss recent experiences related to the development and enhancement of a hybrid stochastic computational fluid dynamics (CFD) solver, the C++ version of the Implicit/Explicit (IMEX) time‐advancement algorithm used in the one‐dimensional turbulence‐based (ODT) large eddy simulation (LES) model, abbreviated as ODTLES. After being ported from Fortran 90, the current capabilities of the C++ code are restricted to the reproducibility of turbulent channel flow simulations with respect to the former Fortran code version that was able to achieve reasonable agreement with available reference direct numerical simulation (DNS) for low to moderate Reynolds number turbulent channel flows. This is far from satisfactory so that current efforts are centered on improving the solver code structure through comprehensive refactoring, robust unit testing, and strict adherence to code style guides, following the principles of Clean Code. We focus the discussion on a methodology to balance unit, regression, and integration testing, here for the LES component of the code. The objective is to frame a starting point that is relevant also for other CFD codes, irrespective of whether they utilize conventional or novel discretization or flow modeling approaches.
Large-eddy simulations (LESs) are known to significantly overestimate entrainment in cloud-topped boundary layers, negatively impacting predictions on cloud mass and cover. This overestimation stems from coarse model resolutions that lead to numerical broadening of the entrainment layer. While it has been shown in direct numerical simulations (DNSs) that down-to-centimetre-scale resolutions can mitigate this issue, such high resolutions are not viable for most applications in the atmospheric sciences. The one-dimensional turbulence model (ODT), introduced by Kerstein [1], offers a computationally efficient alternative that provides full-scale resolution along a 1-D vertical domain. Molecular diffusion is explicitly resolved, while turbulent advection is modelled through a stochastically sampled sequence of spatial mappings, known as eddy events. Physically plausible eddy events are selected based on their current kinetic and potential energy. This allows an accurate representation of local turbulence properties and their dynamical complexity by evolving instantaneous property profiles. This study applies ODT to investigate cloud-top turbulent mixing processes driven by radiative cooling in a smoke cloud, benchmarking the results against DNS. Building on the preliminary findings by Meiselbach [2], we demonstrate improvements in mean profiles and turbulent fluxes of buoyancy and smoke concentration, showing ODT's ability to reproduce salient features observed in DNSs. In addition, we explore convective boundary layer scalings at extended Reynolds and Richardson numbers beyond those accessible in DNS studies.
References:
[1] A. R. Kerstein, Journal of Fluid Mechanics 392, 277334 (1999).
[2] F. T. Meiselbach, Application of ODT to Turbulent Flow Problems, doctoral thesis, BTU Cottbus-Senftenberg (2015).
Investigating cutoff scales in turbulent Ekman flow with a map-based stochastic modeling approach
(2025)
The Ozmidov scale marks the cutoff scale above which overturning fluid motions in stably stratified shear flows are energetically prohibited. A recent study of a turbulent shear layer demonstrates that the Corrsin scale provides an intrinsic cutoff scale when stratification is absent [1]. For the neutral boundary layer, it is proposed by analogy to the mixing layer that the cutoff scale is linked to the Corrsin scale rather than the unbounded Ozmidov scale. The claim is numerically investigated for turbulent Ekman flow with the aid of Kerstein’s one-dimensional turbulence (ODT) model [2], utilizing the case setup described in [3]. ODT offers full-scale resolution along a vertical coordinate by autonomously evolving the instantaneous property profiles. Molecular diffusion is directly resolved, whereas turbulent advection is modeled by a stochastic process that is formulated with the aid of spatial mapping events, which are sampled based on the local available energy. In this model formulation, a cutoff scale is economically prescribed by limiting the sampling range of turbulent scales. Model results in terms of low-order and detailed turbulence statistics will be presented and compared to available reference data and theoretical analysis.
References
[1] F. G. Jacobitz and K. Schneider. Phys. Rev. Fluids 9:044602, 2024.
[2] A. R. Kerstein and S. Wunsch. Bound.-Lay. Meteorol. 118:325–356, 2006.
[3] M. Klein and H. Schmidt. Adv. Sci. Res. 19:117–136, 2022.