Source author record

Stefan Hickel

Stefan Hickel appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

10works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

10 published item(s)

preprint2022arXiv

Permeability and turbulence over perforated plates

We perform direct numerical simulations of turbulent flow at friction Reynolds number $Re_τ\approx 500-2000$ grazing over perforates plates with moderate viscous-scaled orifice diameter $d^+\approx40$--$160$ and analyse the relation between permeability and added drag. Unlike previous studies of turbulent flows over permeable surfaces, we find that flow inside the orifices is dominated by inertial effects, and that the relevant permeability is the Forchheimer and not the Darcy one. We find evidence of a fully rough regime where the relevant length scale is the inverse of the Forchheimer coefficient, which can be regarded as the resistance experienced by the wall-normal flow. Moreover, we show that, for low porosities, the Forchheimer coefficient can be estimated with good accuracy using a simple analytical relation.

preprint2020arXiv

Adaptive Reduced-Order Modeling for Non-Linear Fluid-Structure Interaction

We present an adaptive reduced-order model for the efficient time-resolved simulation of fluid-structure interaction problems with complex and non-linear deformations. The model is based on repeated linearizations of the structural balance equations. Upon each linearization step, the number of unknowns is strongly decreased by using modal reduction, which leads to a substantial gain in computational efficiency. Through adaptive re-calibration and truncation augmentation whenever a non-dimensional deformation threshold is exceeded, we ensure that the reduced modal basis maintains arbitrary accuracy for small and large deformations. Our novel model is embedded into a partitioned, loosely coupled finite volume - finite element framework, in which the structural interface motion within the Eulerian fluid solver is accounted for by a conservative cut-element immersed-boundary method. Applications to the aeroelastic instability of a flat plate at supersonic speeds, to an elastic panel placed within a shock tube, and to the shock induced buckling of an inflated thin semi-sphere demonstrate the efficiency and accuracy of the method.

preprint2020arXiv

Customized data-driven RANS closures for bi-fidelity LES-RANS optimization

Multi-fidelity optimization methods promise a high-fidelity optimum at a cost only slightly greater than a low-fidelity optimization. This promise is seldom achieved in practice, due to the requirement that low- and high-fidelity models correlate well. In this article, we propose an efficient bi-fidelity shape optimization method for turbulent fluid-flow applications with Large-Eddy Simulation (LES) and Reynolds-averaged Navier-Stokes (RANS) as the high- and low-fidelity models within a hierarchical-Kriging surrogate modelling framework. Since the LES-RANS correlation is often poor, we use the full LES flow-field at a single point in the design space to derive a custom-tailored RANS closure model that reproduces the LES at that point. This is achieved with machine-learning techniques, specifically sparse regression to obtain high corrections of the turbulence anisotropy tensor and the production of turbulence kinetic energy as functions of the RANS mean-flow. The LES-RANS correlation is dramatically improved throughout the design-space. We demonstrate the effectiveness and efficiency of our method in a proof-of-concept shape optimization of the well-known periodic-hill case. Standard RANS models perform poorly in this case, whereas our method converges to the LES-optimum with only two LES samples.

preprint2020arXiv

Rapid multi-component phase-split calculations using volume functions and reduction methods

We present a new family of fast and robust methods for the calculation of the vapor-liquid equilibrium at isobaric-isothermal (PT-flash), isochoric-isothermal (VT-flash), isenthalpic-isobaric (HP-flash), and isoenergetic-isochoric (UV-flash) conditions. The framework is provided by formulating phase-equilibrium conditions for multi-component mixtures in an effectively reduced space based on the molar specific value of the recently introduced volume function derived from the Helmholtz free energy. The proposed algorithmic implementation can fully exploit the optimum quadratic convergence of a Newton method with the analytical Jacobian matrix. This paper provides all required exact analytic expressions for the general cubic equation of state. Computational results demonstrate the effectivity and efficiency of the new methods. Compared to conventional methods, the proposed reduced-space iteration leads to a considerable speed-up as well as to improved robustness and better convergence behavior near the spinodal and coexistence curves of multi-component mixtures, where the preconditioning by the reduction method is most effective.

preprint2019arXiv

Rossby Number Effects on Columnar Eddy Formation and the Energy Dissipation Law in Homogeneous Rotating Turbulence

Two aspects of homogeneous rotating turbulence are quantified through forced Direct Numerical Simulations in an elongated domain, which is in the direction of rotation about $340$ times larger than the typical initial eddy size. First, by following the time evolution of the integral length-scale along the axis of rotation $\ell_{\|}$, the growth rate of the columnar eddies and its dependency on the Rossby number $ Ro_{\varepsilon}$ is determined as $γ=4 \exp(-17 Ro_{\varepsilon})$, where $γ$ is the non-dimensional growth rate. Second, a scaling law for the energy dissipation rate $\varepsilon$ is sought. Comparison with current available scaling laws shows that the relation proposed by Baqui & Davidson (2015), i.e., $\varepsilon\sim u'^3/\ell_{\|}$, where $u's$ is the r.m.s. velocity, approximates well part of our data, more specifically the range $0.39\le Ro_{\varepsilon} \le 1.54$. However, relations proposed in the literature fail to model the data for the second and most interesting range, i.e., $0.06\le Ro_{\varepsilon} \le 0.31$, which is marked by the formation of columnar eddies. To find a similarity relation for the latter, we exploit the concept of a spectral transfer time introduced by Kraichnan (1965). Within this framework, the energy dissipation rate is considered to depend on both the non-linear time-scale and the relaxation time-scale. Thus, by analyzing our data, expressions for these different time-scales are obtained that results in~$\varepsilon\sim u'^4/(\ell_{\bot}^2 Ro_{\varepsilon}^{0.62} τ_{nl})$, where $\ell_{\bot}$ is the integral length-scale in the direction normal to the axis of rotation and $τ_{nl}$ is the non-liner time-scale of the initial homogeneous isotropic field.

preprint2016arXiv

Multi-component vapor-liquid equilibrium model for LES and application to ECN Spray A

We present and evaluate a detailed multi-species two-phase thermodynamic equilibrium model for large-eddy simulations (LES) of liquid-fuel injection and mixing at high pressure. The model can represent the coexistence of supercritical states and multi-component subcritical two-phase states. LES results for the transcritical Spray A of the Engine Combustion Network (ECN) are found to agree very well to available experimental data. We also address well-known numerical challenges of trans- and supercritical fluid mixing and compare a fully conservative formulation to a quasi conservative formulation of the governing equations. Our results prove physical and numerical consistency of both methods on fine grids and demonstrate the effects of energy conservation errors associated with the quasi conservative formulation on typical LES grids.

preprint2015arXiv

A cut-cell finite volume - finite element coupling approach for fluid-structure interaction in compressible flow

We present a loosely coupled approach for the solution of fluid-structure interaction problems between a compressible flow and a deformable structure. The method is based on staggered Dirichlet-Neumann partitioning. The interface motion in the Eulerian frame is accounted for by a conservative cut-cell Immersed Boundary method. The present approach enables sub-cell resolution by considering individual cut-elements within a single fluid cell, which guarantees an accurate representation of the time-varying solid interface. The cut-cell procedure inevitably leads to non-matching interfaces, demanding for a special treatment. A Mortar method is chosen in order to obtain a conservative and consistent load transfer. We validate our method by investigating two-dimensional test cases comprising a shock-loaded rigid cylinder and a deformable panel. Moreover, the aeroelastic instability of a thin plate structure is studied with a focus on the prediction of flutter onset. Finally, we propose a three-dimensional fluid-structure interaction test case of a flexible inflated thin shell interacting with a shock wave involving large and complex structural deformations.

preprint2015arXiv

A parametrized non-equilibrium wall-model for large-eddy simulations

Wall-models are essential for enabling large-eddy simulations (LESs) of realistic problems at high Reynolds numbers. The present study is focused on approaches that directly model the wall shear stress, specifically on filling the gap between models based on wall-normal ordinary differential equations (ODEs) that assume equilibrium and models based on full partial differential equations (PDEs) that do not. We develop ideas for how to incorporate non-equilibrium effects (most importantly, strong pressure-gradient effects) in the wall-model while still solving only wall-normal ODEs. We test these ideas using two reference databases: an adverse pressure-gradient turbulent boundary-layer and a shock/boundary-layer interaction problem, both of which lead to separation and re-attachment of the turbulent boundary layer.

preprint2014arXiv

Benchmarking in a rotating annulus: a comparative experimental and numerical study of baroclinic wave dynamics

The differentially heated rotating annulus is a widely studied tabletop-size laboratory model of the general mid-latitude atmospheric circulation. The two most relevant factors of cyclogenesis, namely rotation and meridional temperature gradient are quite well captured in this simple arrangement. The radial temperature difference in the cylindrical tank and its rotation rate can be set so that the isothermal surfaces in the bulk tilt, leading to the formation of baroclinic waves. The signatures of these waves at the free water surface have been analyzed via infrared thermography in a wide range of rotation rates (keeping the radial temperature difference constant) and under different initial conditions. In parallel to the laboratory experiments, five groups of the MetStröm collaboration have conducted numerical simulations in the same parameter regime using different approaches and solvers, and applying different initial conditions and perturbations. The experimentally and numerically obtained baroclinic wave patterns have been evaluated and compared in terms of their dominant wave modes, spatio-temporal variance properties and drift rates. Thus certain ``benchmarks'' have been created that can later be used as test cases for atmospheric numerical model validation.

preprint2014arXiv

DNS and LES of two-phase flows with cavitation

We report on recent progress in the physical and numerical modeling of compressible two-phase flows that involve phase transition between the liquid and gaseous state of the fluid. The high-speed dynamics of cavitation bubbles is studied in well-resolved simulations (DNS) with a sharp-interface numerical model on a micro scale. The underlying assumption of the employed evaporation/condensation model is that phase change occurs in thermal non-equilibrium and that the associated timescale is larger than that of the wave dynamics. Results for the collapse of a spherical vapor bubble close to a solid wall are discussed for three different bubble-wall configurations. The major challenge for such numerical investigations is to accurately reproduce the dynamics of the interface between liquid and vapor during the entire collapse process, including the high-speed dynamics of the late stages, where compressibility of both phases plays a decisive role. Direct interface resolving simulations are intractable for real world technical applications such as turbulent flows involving cavitation clouds with millions of vapor bubbles and a wide range of time and length scales. For this reason, we developed a coarse grained model for large-eddy simulation (LES) of turbulent two-phase flows with cavitation. In LES, vapor bubbles constitute sub-grid scales that have to be modeled accordingly. On the grid scale, we solve the compressible Navier-Stokes equations for a homogeneous mixture of liquid and vapor. For deriving an appropriate equation of state, we assume that the characteristic time scale of phase change is much smaller than the numerical time step and that phase change is isentropic and in mechanical equilibrium at the saturation pressure. This macroscopic model is applied to realistic technical problems and to more generic flows to study the mutual interaction of turbulence and cavitation.