Researcher profile

Philipp Schlatter

Philipp Schlatter contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

10 published item(s)

preprint2022arXiv

An adverse-pressure-gradient turbulent boundary layer with nearly-constant $\boldsymbol{β\simeq 1.4}$ up to $\boldsymbol{Re_θ \simeq 8,700}$

In this study, a new well-resolved large-eddy-simulation (LES) of an incompressible near-equilibrium adverse-pressure-gradient (APG) turbulent boundary layer (TBL) over a flat plate is presented. In this simulation, we have established a near-equilibrium APG over a wide Reynolds-number range. In this so-called region of interest (ROI), the Clauser--Rotta pressure-gradient parameter $β$ exhibits an approximately constant value of around 1.4, and the Reynolds number based on momentum thickness reaches $Re_θ=8700$. To the authors' knowledge, this is to date the highest $Re_θ$ achieved for a near-equilibrium APG TBL under an approximately constant moderate APG. We evaluated the self-similarity of the outer region using two scalings, namely the Zagarola--Smits and an alternative one based on edge velocity and displacement thickness. Our results reveal that outer-layer similarity is achieved, and the viscous scaling collapses the near-wall region of the mean flow in agreement with classical theory. Spectral analysis reveals that the APG displaces some small-scale energy from the near-wall to the outer region, an effect observed for all the components of the Reynolds-stress tensor, which becomes more evident at higher Reynolds numbers. Generally, the effects of the APG are more noticeable at lower Reynolds numbers. For instance, the outer peak of turbulent-kinetic-energy (TKE) production is less prominent at higher $Re$. While the scale separation increases with $Re$ in zero-pressure-gradient (ZPG) TBLs, this effect becomes accentuated by the APG. Despite the reduction of the outer TKE production at higher Reynolds numbers, the mechanisms of energization of large scales are still present.

preprint2022arXiv

Large-Scale Direct Numerical Simulations of Turbulence Using GPUs and Modern Fortran

We present our approach to making direct numerical simulations of turbulence with applications in sustainable shipping. We use modern Fortran and the spectral element method to leverage and scale on supercomputers powered by the Nvidia A100 and the recent AMD Instinct MI250X GPUs, while still providing support for user software developed in Fortran. We demonstrate the efficiency of our approach by performing the world's first direct numerical simulation of the flow around a Flettner rotor at Re=30'000 and its interaction with a turbulent boundary layer. We present one of the first performance comparisons between the AMD Instinct MI250X and Nvidia A100 GPUs for scalable computational fluid dynamics. Our results show that one MI250X offers performance on par with two A100 GPUs and has a similar power efficiency.

preprint2022arXiv

New insight into the spectra of turbulent boundary layers with pressure gradients

With the availability of new high-Reynolds-number ($Re$) databases of turbulent boundary layers (TBLs) it has been possible to identify in detail certain regions of the boundary layer with more complex behavior. In this study we consider a unique database at moderately-high $Re$, with a near-constant adverse pressure gradient (APG) (Pozuelo {\it et al.}, {\it J. Fluid Mech.}, {\bf 939}, A34, 2022), and perform spectral analysis of the Reynolds stresses, focusing on the streamwise component. We assess different regions of the APG TBL, comparing this case with the zero-pressure-gradient (ZPG) TBL, and identify the relevant scaling parameters as well as the contribution of the scales of different sizes. The small scales in the near-wall region up to the near-wall spectral peak have been found to scale using viscous units. In APG TBLs, the largest scales close to the wall have a better scaling with the boundary-layer thickness ($δ_{99}$), and they are significantly affected by the APG. In the overlap and wake regions of the boundary layer, the small energetic scales exhibit a good scaling with the displacement thickness ($δ^*$) while the larger scales and the outer spectral peak are better scaled with the boundary-layer thickness. Also note that the wall-normal location of the spectral outer peak scales with the displacement thickness rather than the boundary layer thickness. The various scalings exhibited by the spectra in APG TBLs are reported here for the first time, and shed light on the complex phenomena present in these flows of great scientific and technological importance.

preprint2022arXiv

Predicting the temporal dynamics of turbulent channels through deep learning

The success of recurrent neural networks (RNNs) has been demonstrated in many applications related to turbulence, including flow control, optimization, turbulent features reproduction as well as turbulence prediction and modeling. With this study we aim to assess the capability of these networks to reproduce the temporal evolution of a minimal turbulent channel flow. We first obtain a data-driven model based on a modal decomposition in the Fourier domain (which we denote as FFT-POD) of the time series sampled from the flow. This particular case of turbulent flow allows us to accurately simulate the most relevant coherent structures close to the wall. Long-short-term-memory (LSTM) networks and a Koopman-based framework (KNF) are trained to predict the temporal dynamics of the minimal-channel-flow modes. Tests with different configurations highlight the limits of the KNF method compared to the LSTM, given the complexity of the flow under study. Long-term prediction for LSTM show excellent agreement from the statistical point of view, with errors below 2% for the best models with respect to the reference. Furthermore, the analysis of the chaotic behaviour through the use of the Lyapunov exponents and of the dynamic behaviour through Poincaré maps emphasizes the ability of the LSTM to reproduce the temporal dynamics of turbulence. Alternative reduced-order models (ROMs), based on the identification of different turbulent structures, are explored and they continue to show a good potential in predicting the temporal dynamics of the minimal channel.

preprint2021arXiv

Physics-informed neural networks for solving Reynolds-averaged Navier$\unicode{x2013}$Stokes equations

Physics-informed neural networks (PINNs) are successful machine-learning methods for the solution and identification of partial differential equations (PDEs). We employ PINNs for solving the Reynolds-averaged Navier$\unicode{x2013}$Stokes (RANS) equations for incompressible turbulent flows without any specific model or assumption for turbulence, and by taking only the data on the domain boundaries. We first show the applicability of PINNs for solving the Navier$\unicode{x2013}$Stokes equations for laminar flows by solving the Falkner$\unicode{x2013}$Skan boundary layer. We then apply PINNs for the simulation of four turbulent-flow cases, i.e., zero-pressure-gradient boundary layer, adverse-pressure-gradient boundary layer, and turbulent flows over a NACA4412 airfoil and the periodic hill. Our results show the excellent applicability of PINNs for laminar flows with strong pressure gradients, where predictions with less than 1% error can be obtained. For turbulent flows, we also obtain very good accuracy on simulation results even for the Reynolds-stress components.

preprint2020arXiv

An Uncertainty-Quantification Framework for Assessing Accuracy, Sensitivity, and Robustness in Computational Fluid Dynamics

A framework is developed based on different uncertainty quantification (UQ) techniques in order to assess validation and verification (V&V) metrics in computational physics problems, in general, and computational fluid dynamics (CFD), in particular. The metrics include accuracy, sensitivity and robustness of the simulator's outputs with respect to uncertain inputs and computational parameters. These parameters are divided into two groups: based on the variation of the first group, a computer experiment is designed, the data of which may become uncertain due to the parameters of the second group. To construct a surrogate model based on uncertain data, Gaussian process regression (GPR) with observation-dependent (heteroscedastic) noise structure is used. To estimate the propagated uncertainties in the simulator's outputs from first and also the combination of first and second groups of parameters, standard and probabilistic polynomial chaos expansions (PCE) are employed, respectively. Global sensitivity analysis based on Sobol decomposition is performed in connection with the computer experiment to rank the parameters based on their influence on the simulator's output. To illustrate its capabilities, the framework is applied to the scale-resolving simulations of turbulent channel flow using the open-source CFD solver Nek5000. Due to the high-order nature of Nek5000 a thorough assessment of the results' accuracy and reliability is crucial, as the code is aimed at high-fidelity simulations. The detailed analyses and the resulting conclusions can enhance our insight into the influence of different factors on physics simulations, in particular the simulations of wall-bounded turbulence.

preprint2020arXiv

On the use of recurrent neural networks for predictions of turbulent flows

In this paper, the prediction capabilities of recurrent neural networks are assessed in the low-order model of near-wall turbulence by Moehlis {\it et al.} (New J. Phys. {\bf 6}, 56, 2004). Our results show that it is possible to obtain excellent predictions of the turbulence statistics and the dynamic behavior of the flow with properly trained long short-term memory (LSTM) networks, leading to relative errors in the mean and the fluctuations below $1\%$. We also observe that using a loss function based only on the instantaneous predictions of the flow may not lead to the best predictions in terms of turbulence statistics, and it is necessary to define a stopping criterion based on the computed statistics. Furthermore, more sophisticated loss functions, including not only the instantaneous predictions but also the averaged behavior of the flow, may lead to much faster neural network training.

preprint2020arXiv

Optimization of Tensor-product Operations in Nekbone on GPUs

In the CFD solver Nek5000, the computation is dominated by the evaluation of small tensor operations. Nekbone is a proxy app for Nek5000 and has previously been ported to GPUs with a mixed OpenACC and CUDA approach. In this work, we continue this effort and optimize the main tensor-product operation in Nekbone further. Our optimization is done in CUDA and uses a different, 2D, thread structure to make the computations layer by layer. This enables us to use loop unrolling as well as utilize registers and shared memory efficiently. Our implementation is then compared on both the Pascal and Volta GPU architectures to previous GPU versions of Nekbone as well as a measured roofline. The results show that our implementation outperforms previous GPU Nekbone implementations by 6-10%. Compared to the measured roofline, we obtain 77 - 92% of the peak performance for both Nvidia P100 and V100 GPUs for inputs with 1024 - 4096 elements and polynomial degree 9.

preprint2020arXiv

The edge as a Lagrangian Coherent Structure in a high-dimensional state space

Dissipative dynamical systems characterised by two basins of attraction are found in many physical systems, notably in hydrodynamics where laminar and turbulent regimes can coexist. The state space of such systems is structured around a dividing manifold called the edge, which separates trajectories attracted by the laminar state from those reaching the turbulent state. We apply here concepts and tools from Lagrangian data analysis to investigate this edge manifold. This approach is carried out in the state space of automous arbitrarily high-dimensional dissipative systems, in which the edge manifold is re-interpreted as a Lagrangian Coherent Structure (LCS). Two different diagnostics, finite-time Lyapunov exponents and Lagrangian Descriptors, are used and compared with respect to their ability to identify the edge and to their scalability. Their properties are illustrated on several low-order models of subcritical transition of increasing dimension and complexity, as well on well-resolved simulations of the Navier-Stokes equations in the case of plane Couette flow. They allow for a mapping of the global structure of both the state space and the edge manifold based on quantitative information. Both diagnostics can also be used to generate efficient bisection algorithms to approach asymptotic edge states, which outperform classical edge tracking.

preprint2019arXiv

A critical point for bifurcation cascades and featureless turbulence

In this Letter we show that a bifurcation cascade and fully sustained turbulence can share the phase space of a fluid flow system, resulting in the presence of competing stable attractors. We analyse the toroidal pipe flow, which undergoes subcritical transition to turbulence at low pipe curvatures and supercritical transition at high curvatures, as was previously documented. We provide decisive evidence that the nature of the supercritical transition is of Ruelle--Takens type and that, in a narrow range of intermediate curvatures, its dynamics competes with that of sustained turbulence emerging through subcritical transition mechanisms.