TY - GEN A1 - Vallem, Rishindra A1 - Klein, Marten A1 - Schmidt, Heiko T1 - Capabilities and limitations of smoothed particle hydrodynamics for the simulation of two‐phase flow instabilities T2 - Proceedings in Applied Mathematics and Mechanics N2 - Smoothed particle hydrodynamics (SPH) is a mesh‐free, Lagrangian particle‐based method that is able to simulate multiphase flows in an economical manner. However, its ability to capture the flow regimes and regime transitions in two phase (liquid‐gas) internal flows, such as pipe or channel flows is not yet generally established. To address this lack in understanding, we first examine a laminar rising bubble case in order to evaluate the fluid‐fluid interface representation and transient interface evolution by the solver. With a focus towards the transition mechanism from a stratified flow regime to a slug flow regime, we investigate the Kelvin–Helmholtz instability (KHI) both qualitatively and quantitatively, initially focusing on a low density ratio () and then extending it to a high density ratio (). For the low density ratio, we conduct an analysis of the temporal evolution and demonstrate that the SPH solver captures the initial exponential growth in qualitative agreement with inviscid linear stability theory (LST) and reference numerical data for shear‐dominated flow with Richardson number . By conducting eight additional simulations for various for the high density ratio, we demonstrate that the numerically obtained parameter value for instability is around , which is in reasonable agreement with the theoretically expected value of . Based on the SPH results obtained for the range , we suggest a simple parameterization of the reduction of the effective growth rate proportional to . KW - Smoothed Particle Hydrodynamics KW - two-phase flow KW - Kelvin-Helmholtz instability KW - rising bubble KW - growth rate Y1 - 2024 U6 - https://doi.org/10.1002/pamm.202400206 SN - 1617-7061 VL - 24/2024 PB - Wiley CY - Weinheim ER - TY - GEN A1 - Klein, Marten A1 - Kessler, Robert William A1 - Schmidt, Heiko T1 - Numerical investigation of drag reduction effects on a track bicycle fork using wings with a wavy leading edge T2 - Proceedings in Applied Mathematics and Mechanics N2 - Reynolds‐averaged Navier–Stokes (RANS) and large‐eddy simulations (LES) of the flow around wings with a wavy leading edge (WLE) are conducted in order to assess the capabilities of a passive flow control strategy for drag reduction. The intended application is indoor track cycling with controlled flow conditions. A section of a single fork rod is investigated in order to make the numerical simulations feasible. The present study reveals that net drag reduction is possible by a nonsinusoidal modification of the leading edge of the wing. However, the drag reduction effect remains limited to a few percent. While RANS and LES yield the same drag coefficient for a reference case, RANS underestimates the drag reduction effect for a longer wing and the WLE cases, but exhibits otherwise a qualitatively similar trend as the LES. With the aid of RANS, an optimal geometry is obtained defined by the wavelength‐to‐chord length ratio of and the amplitude‐to‐chord length ratio of . Corresponding LES results give an indication of the origin of drag reduction by a hampered vortex shedding. The generation of smaller and more streamwise oriented vortical flow structures at the trailing edge and behind the WLE wing is correlated with significantly reduced lift fluctuations and drag reduction. KW - sports aerodynamics KW - passive flow control KW - aerodynamic drag KW - drag reduction KW - computational fluid dynamics KW - turbulence modeling Y1 - 2024 U6 - https://doi.org/10.1002/pamm.202400178 SN - 1617-7061 VL - 24 SP - 1 EP - 14 PB - Wiley ER - TY - GEN A1 - Klein, Marten A1 - Schmidt, Heiko T1 - Capturing features of transient boundary layers with a map-based stochastic modeling approach N2 - Atmospheric boundary layers (ABLs) exhibit transient processes on various time and length scales, with a scale separation between the large-scale forcing and the small-scale response. Some crucial but standing challenges in modeling and simulation of ABL flows lie in the detailed representation of boundary layer turbulence (e.g. [1]). This includes intermittent and transient processes and the resulting turbulent and laminar response mechanisms. State-of-the-art subgrid-scale models utilize statistical closures for an averaged resolved flow state on the basis of the Monin-Obhukov similarity theory (MOST) to represent scalar fluxes and momentum fluxes (e.g. [2]). Fluctuations are not resolved in MOST. Instead, their ensemble effect is parameterized by the resolved large scales, neglecting backscatter from the unresolved small scales. Data-driven stochastic approaches aim to incorporate fluctuations and the spontaneous occurrence of instabilities, but at the expense of ad hoc forcings (e.g. [3]). The mentioned limitations can be removed by a physically compatible representation of turbulent fluctuations. This is addressed here by utilization of a map-based stochastic approach that is based on the one-dimensional turbulence (ODT) model [4]. ODT autonomously evolves vertical flow profiles for prescribed initial and boundary conditions, and physical forcings. The model captures turbulent cascade phenomenology and aims to resolve all relevant turbulent scales along a physical coordinate. Turbulent advection is modeled by a stochastically sampled sequence of spatial mapping events that punctuate the deterministic advancement due to viscous and Coriolis forces. The offered dynamical complexity removes the need for artificial forcings. In the contribution, key results from recent and ongoing studies related to the reduced-order modeling of ABL flows will be presented. First, surface scalar and momentum fluxes in turbulent channels are discussed emphasizing the correctly predicted inapplicability of the Reynolds analogy [5]. Second, the influence of system rotation and stratification is discussed for low-order velocity statistics and the participating turbulent scales [6,7]. Third, results for nonequilibrium conditions are presented for a transient ABL that exhibits turbulent bursts in response to an oscillatory geostrophic forcing [8]. Last, some preliminary results on the stochastic deconvolution of averaged data [9] will be presented focusing on the additional physical insight that is offered by the model. References [1] L. Mahrt. Annu. Rev. Fluid Mech. 46:23–45, 2014. [2] I. Stiperski, and M. Calaf. Phys. Rev. Lett. 130:124001, 2023. [3] V. Boyko, and N. Vercauteren. Q. J. R. Meteorol. Soc. 149(755):2125–2145, 2023. [4] A. R. Kerstein, and S. Wunsch. Bound.-Lay. Meteorol. 118:325–356, 2006. [5] M. Klein, H. Schmidt, and D. Lignell. Int. J. Heat Fluid Flow 93:108889, 2022. [6] M. Klein, and H. Schmidt. Adv. Sci. Res. 19:117–136, 2022. [7] L. S. Freire. Bound.-Lay. Meteorol. 184:25–43, 2022. [8] M. Klein, and H. Schmidt. Adv. Sci. Res. 20:55–64, 2023. [9] C. Glawe, M. Klein, and H. Schmidt. Proc. Appl. Math. Mech. 23:e202300055, 202 KW - atmospheric boundary layer KW - stochastic modeling KW - one-dimensional turbulence KW - transient effects Y1 - 2024 U6 - https://doi.org/10.5194/egusphere-egu24-15560 PB - Copernicus GmbH ER - TY - GEN A1 - Klein, Marten A1 - Medina Méndez, Juan Alí A1 - Schmidt, Heiko T1 - Resolving the electrostatic boundary layer in a turbulent electrohydrodynamic flow with a map-based stochastic modeling approach T2 - 1st European Fluid Dynamics Conference - Daily Scientific Program N2 - Electrohydrodynamically (EHD) enhanced wall-bounded turbulent flows are encountered in various technical applications ranging from air-cleaning devices (like precipitators) to electrolyte flows (like redox flow batteries). The multi-physical processes governing the flow properties, however, are not yet very well understood. This is due the inaccessibility to and limitations of measurement equipment and numerical resolution requirements imposed by the electrostatic charge boundary layer that interacts on different time scales with the turbulent boundary layer. Recent advances in measurement techniques allow to resolve the exponential electrostatic charge boundary layer in a charged particle-ladden, weakly turbulent gas flow [1], which is qualitatively compatible with direct numerical simulation (DNS) results for a weakly turbulent flow of an electrolyte [2] at low Reynolds (Re) number. The challenge that remains is the extrapolation to highly turbulent flow conditions. Charged particles are heavy compared with fluid molecules and can be treated as a high Schmidt (Sc) number scalar, which is a burden for DNS. In this contribution, it is demonstrated that this burden can be overcome for the transient boundary layer evolution by utilizing a radically reduced, map-based stochastic one-dimensional turbulence (ODT) model. The model predicts a significant enhancement of the skin friction drag due to turbulence-induced screening layer depletion [3], as summarized in Fig. 1. In the talk, features of the instantaneous and mean velocity and electrostatic boundary layer will be presented. The plausibility of the model prediction is evaluated on a physical basis, encompassing details of the model formulation and the emerging hydrodynamic and electrokinetic properties of the boundary layer. References [1] W. Xu, S. Jantaˇc, T. Matsuyama, and H. Grosshans. arXiv:2306.06970, 2023. (Accepted for publication by Exp. Fluids.) [2] R. Ostilla-M´onico, and A. A. Lee. Faraday Discuss., 199:159–173, 2017. [3] M. Klein, J. A. Medina M´endez, and H. Schmidt. Tech. Mech., 43:111–127, 2023. KW - EHD turbulence KW - electric double layer KW - boundary layer KW - Schmidt number effects KW - one-dimensional turbulence KW - stochastic modeling KW - skin friction drag Y1 - 2024 UR - https://www.conftool.org/efdc1/index.php?page=browseSessions&presentations=show&search=822+klein UR - https://www.aia.rwth-aachen.de/fileadmin/user_upload/Daily_Scientific_Program.pdf CY - Aachen ER - TY - GEN A1 - Parekh, Parshva Atulbhai A1 - Gschwander, Stefan A1 - Klein, Marten A1 - Gamisch, Sebastian T1 - CFD-based analysis and minimization of mixing during the charging phase of a thermal energy storage tank T2 - 1. Jahresfachtagung des Energie-Innovationszentrums (EIZ) Cottbus N2 - Transient numerical simulations are performed for a cuboidal storage tank in order to resolve the transient features of the charging phase in a feasible manner. A diffuser is utilized for the inflow of water in order avoid large-scale overturning fluid motions in order to establish a thermal stratification. It is demonstrated that the thermal stratification can be furrther enhanced by introducing a layer of a porous medium at the top of the storage tank based on computational fluid dynamics (CFD) simulations using COMSOL Multiphysics®. Initially, a storage tank configuration without a porous medium is simulated in order to establish a baseline understanding. Subsequent simulations are performed systematically varying various parameters of the porous medium, such as porosity, location, and inclination. The inclined placement of the porous sheet yields a reduction of the thermocline thickness by approximately 38% compared to the reference case, thereby significantly enhancing the thermal stratification. As next step, the results obtained will be verified in an experimental apparatus at Fraunhofer ISE. In the talk, the set-up of the numerical model, including the treatment of the porous sheet, the thermocline evolution together with the governing fluid flow, and the effect of an additionally installed porous sheet will be discussed KW - thermal energy storage tank KW - mixing KW - porous media KW - computational fluid dynamics Y1 - 2024 UR - https://www-docs.b-tu.de/fg-stroemungsmodellierung/public/Klein_2024_EIZ1_abstract.pdf CY - Cottbus ER - TY - GEN A1 - Polasanapalli, Sai Ravi Gupta A1 - Klein, Marten A1 - Schmidt, Heiko T1 - Investigation of the impact of transient pressure gradients on turbulent channel flow dynamics T2 - Proceedings in Applied Mathematics and Mechanics N2 - The effect of transient pressure gradients, or transient pumping in turbulent channel flow configuration, is investigated. Employing a cost-efficient reduced-order stochastic method known as one-dimensional turbulence (ODT) modeling, simulations explore different signal shapes for the modulation of the prescribed pressure gradient forcing, including sinusoidal modulation, step-like modulation, and piecewise sinusoidal beating. Various cycle periods and active pumping times are investigated. The simulations are conducted at a frictional Reynolds number of $Re_\tau = 395$ and a molecular Prandtl number of $Pr = 1$. The study adopts a passive scalar formulation to investigate heat transfer properties. The study quantifies the effects of transient pressure gradients on heat transfer rate, drag, and pumping power. Preliminary ODT predictions suggest that all transient cases exhibit lower heat transfer rates and a higher pumping power requirement than the constant pressure gradient case, with the step-like modulation yields superior skin-friction drag and heat transfer rate reductions relative to other signals. KW - turbulent heat transfer KW - turbulent drag KW - pumping power KW - one-dimensional turbulence KW - pulsating flow Y1 - 2024 U6 - https://doi.org/10.1002/pamm.202400183 SN - 1617-7061 VL - 24 SP - 1 EP - 12 PB - Wiley-VCH GmbH CY - Weinheim ER - TY - GEN A1 - Polasanapalli, Sai Ravi Gupta A1 - Klein, Marten A1 - Schmidt, Heiko T1 - Investigating modifications of the heat transfer by velocity boundary conditions in turbulent thermal convection using an off-lattice Boltzmann method T2 - STAB Jahresbericht 2024 N2 - In our contribution to the STAB workshop, we will present the effect of different surfaces on flow and heat transfer characteristics, taking into account the no-slip, free-slip, and finite-slip lengths of the walls. Following that, the impact of different walls, such as horizontal or side walls, will be discussed separately. Finally, the influence of finite slip length on turbulent characteristics will also be addressed. The work was supported by the North- German Supercomputing Alliance (HLRN) and numerical simulations were carried out on HLRN high-performance computing facilities at Berlin and Göttingen (project ID: bbi00022). KW - lattice Boltzmann method KW - thermal convection KW - hydrophilic and hydrophobic surfaces KW - Navier slip boundary condition Y1 - 2024 UR - https://bfm.rcbe.de/files/2024/11/2024_STAB-Jahresbericht.pdf SP - 40 EP - 41 PB - Deutsche Strömungsmechanische Arbeitsgemeinschaft, STAB CY - Regensburg [et al.] ER - TY - CHAP A1 - Klein, Marten A1 - Tsai, Pei-Yun A1 - Schmidt, Heiko ED - Dillmann, Andreas ED - Heller, Gerd ED - Krämer, Ewald ED - Wagner, Claus ED - Weiss, Julien T1 - Stochastic Modeling and Large-Eddy Simulation of Heated Concentric Coaxial Pipes T2 - New Results in Numerical and Experimental Fluid Mechanics XIV, STAB/DGLR Symposium 2022 N2 - 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. KW - Heat and mass transfer KW - Stochastic turbulence modeling KW - Spanwise curvature effects KW - Pipe flow Y1 - 2024 SN - 978-3-031-40482-5 SN - 978-3-031-40481-8 U6 - https://doi.org/10.1007/978-3-031-40482-5_41 SN - 1612-2909 SN - 1860-0824 N1 - This research is supported by the German Federal Government, the Federal Ministry of Education and Research and the State of Brandenburg within the framework of the joint project EIZ: Energy Innovation Center (project numbers 85056897 and 03SF0693A) with funds from the Structural Development Act (Strukturstärkungsgesetz) for coal-mining regions. M.K. acknowledges support by the BTU Graduate Research School (Conference Travel Grant). SP - 435 EP - 444 PB - Springer CY - Cham ER - TY - GEN A1 - Klein, Marten T1 - Map-based stochastic simulation data of a transient Ekman boundary layer N2 - A journal paper in Advances in Science and Research details the numerical modeling approach used to create the data. Here, the model input files, the raw data, processed data, and plot scripts are provided that support the research. The code used to generate the data is an extended version of the one-dimensional turbulence (ODT) model. The current model implementation utilizes an adaptive grid that further increases numerical efficiency. A reduced version of the adaptive ODT code used in this work is available free of charge at: https://github.com/BYUignite/ODT The bash script makePlot.sh is the top-level driver and contains all additional information about the cases. Some other Details are provided by low-level README files. Python-3.8 is required to run the scripts. KW - atmospheric turbulence KW - boundary layer theory KW - multiscale flow KW - reduced-order stochastic modeling Y1 - 2023 ER - TY - GEN A1 - Polasanapalli, Sai Ravi Gupta A1 - Klein, Marten A1 - Schmidt, Heiko T1 - Towards stochastic subgrid-scale modeling of turbulent thermal convection in an under-resolved off-lattice Boltzmann method T2 - Proceedings in applied mathematics and mechanics : PAMM N2 - A characteristic-based Off-Lattice Boltzmann Method (OLBM) and a stochastic One-Dimensional Turbulence (ODT) model is utilized for numerical simulation of turbulent thermal convection. Standalone ODT results for low-order statistics are compared with those from various eddy-viscosity-based subgrid-scale models utilized in Large-Eddy Simulations (LES) with OLBM. The predictive capabilities of both approaches are discussed by comparison with available reference Direct Numerical Simulation (DNS) results. All turbulence models are able to predicted the mean temperature, but fail to fully capture fluctuations. While the OLBM aims to represent large-scale structures, it misses some constitutional small-scale fluctuations. By contrast, the reduced-order ODT model captures small-scale fluctuations in the vicinity of the wall, but cannot resolve the organized bulk flow. Here, the modeling capabilities of both OLBM and ODT as standalone tools are discussed. On this basis, a strategy for the incorporation of ODT as wall model in OLBM is suggested. KW - thermal convection KW - heat transfer KW - lattice Boltzmann method KW - turbulence modeling KW - one-dimensional turbulence Y1 - 2023 UR - https://onlinelibrary.wiley.com/doi/10.1002/pamm.202300223 U6 - https://doi.org/10.1002/pamm.202300223 SN - 1617-7061 SP - 1 EP - 9 ER - TY - GEN A1 - Klein, Marten A1 - Zenker, Christian A1 - Starick, Tommy A1 - Schmidt, Heiko T1 - Stochastic modeling of multiple scalar mixing in a three-stream concentric coaxial jet based on one-dimensional turbulence T2 - International Journal of Heat and Fluid Flow N2 - 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. KW - map-based stochastic advection modeling KW - multiple passive scalars KW - one-dimensional turbulence KW - turbulent jet KW - turbulent mixing Y1 - 2023 U6 - https://doi.org/10.1016/j.ijheatfluidflow.2023.109235 SN - 0142-727X N1 - This article is part of the "TSFP12 Special Issue". VL - 104 SP - 1 EP - 17 ER - TY - GEN A1 - Tsai, Pei-Yun A1 - Schmidt, Heiko A1 - Klein, Marten T1 - Investigating Reynolds number effects in turbulent concentric coaxial pipe flow using stochastic one-dimensional turbulence modeling T2 - Proceedings in Applied Mathematics and Mechanics N2 - The present study numerically investigates turbulent momentum transfer in concentric coaxial (annular) pipe flow with small radius ratios (𝜂 = 𝑅𝑖∕𝑅𝑜 = 0.1,0.04,0.02). To model the flow, a stochastic one-dimensional turbulence (ODT) model formulated for cylindrical geometry is used that provides full-scale resolution along a representative radial coordinate. The present investigation extends the model validation by Tsai et al. (PAMM,22:e202200272, 2023), to radius ratios smaller than 0.1 and addresses boundary layers with strong span-wise curvature effects. The focus is on the assessment and analysis of statistical flow features in the vicinity of the inner cylinder wall, particularly in cases with small radius ratios. Following Boersma & Breugem (Flow Turbul. Com-bust.,86:113–127, 2011), classical boundary-layer and mixing-length theory is utilized to analyze the model predictions. The results demonstrate that the ODT model captures leading-order curvature and mixing-length effects by its physics-compatible construction. Utilizing the model for extrapolation to high Reynolds numbers inaccessible to conventional high-fidelity numerical approaches shows that curvature effects persist and nonlocally affect the entire boundary layer. The model results provide support for a spanwise-curvature-modified wall function. KW - turbulent coaxial pipe flow KW - stochastic modeling KW - one-dimensional turbulence KW - spanwise curvature KW - turbulent boundary layer KW - high Reynolds number Y1 - 2023 U6 - https://doi.org/10.1002/pamm.202300167 SN - 1617-7061 N1 - Proceedings of the GAMM Annual Meeting 2023, Dresden, Germany VL - 23 IS - 4 SP - 1 EP - 8 ER - TY - GEN A1 - Klein, Marten A1 - Starick, Tommy A1 - Zenker, Christian A1 - Medina Méndez, Juan Alí A1 - Schmidt, Heiko T1 - Reduced order stochastic modeling of turbulent mixing based on conservative baker’s maps T2 - Proceedings of the 14th International ERCOFTAC Symposium on Engineering Turbulence Modelling and Measurements (ETMM-14) N2 - The detailed numerical representation of turbulent mixing processes is a standing challenge for non-premixed chemically reacting flows. The full range of relevant flow scales needs to be captured and it is also necessary to distinguish turbulent advective from molecular diffusive processes in order to represent Reynolds and Schmidt number effects. These requirements are addressed here by utilizing two different map-based stochastic turbulence modeling strategies. The one-dimensional turbulence (ODT) model utilizes event-based turbulence modeling, whereas the hierarchical parcel-swapping (HiPS) model is a fully event-based mixing model. ODT provides full-scale resolution at affordable costs by dimensional model reduction based on the boundary-layer approximation to shear flow. HiPS is far less costly than ODT but currently limited to locally homogeneous isotropic turbulence. The physics-compatible modeling capabilities with respect to phase-space representation of turbulent mixing are demonstrated for two canonical cases using standalone model formulations. KW - turbulent mixing KW - one-dimensional turbulence (ODT) KW - hierarchical parcel swapping (HiPS) KW - stochastic turbulence modeling KW - round jet KW - passive scalar Y1 - 2023 UR - https://etmm.ercoftac.org/etmm/program/conference-program/ UR - https://drive.google.com/file/d/1q2BDOO5bXfqq0Y4z4HCGndiFI033bPyg/view?usp=drive_link SP - 613 EP - 618 PB - ERCOFTAC CY - Barcelona, Spain ER - TY - GEN A1 - Vallem, Rishindra A1 - Klein, Marten A1 - Schmidt, Heiko T1 - Numerical modeling and simulation of two-phase internal flow instabilities using Smoothed Particle Hydrodynamics (SPH) T2 - STAB Jahresbericht 2023 KW - Smoothed Particle Hydrodynamics (SPH) KW - Kelvin-Helmholtz instability KW - linear stability analysis KW - two-phase flow Y1 - 2023 UR - https://www.dlr.de/as/Portaldata/5/Resources/dokumente/veranstaltungen/stab_workshop/Jahresbericht2023.pdf UR - https://www.dlr.de/as/desktopdefault.aspx/tabid-128/268_read-1678/ VL - 21/2023 SP - 158 EP - 159 PB - Deutsche Strömungsmechanische Arbeitsgemeinschaft (STAB) CY - Göttingen, Germany ER - TY - GEN A1 - Tsai, Pei-Yun A1 - Schmidt, Heiko A1 - Klein, Marten T1 - Stochastic modeling of asymmetric turbulent boundary layers in annular pipe flow Y1 - 2023 UR - https://www-docs.b-tu.de/fg-stroemungsmodellierung/public/Tsai_2023_poster_PragueSummerSchoolX.pdf U6 - https://doi.org/10.13140/RG.2.2.18789.78567 N1 - Stochastics in Fluids Summer School and Workshop 2023, Czech Academy of Sciences Institute of Mathematics, Prague, Czech Republic ER -