Source author record

Philipp Schlatter

Philipp Schlatter 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

20works
7topics
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

20 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.

preprint2016arXiv

Aspect ratio effect on particle transport in turbulent duct flows

The dynamics of dilute micron-sized spherical inertial particles in turbulent duct flows is studied by means of direct numerical simulations of the carrier phase turbulence with one-way coupled Lagrangian particles. The geometries are a square and a rectangular duct with width-to-height aspect ratio $AR$ of 3 operating at $Re_{τ,c}=360$ . The present study is designed to determine the effect of turbulence-driven secondary motion on the particle dynamics. Our results show that a weak cross-flow secondary motion significantly changes the cross-sectional map of the particle concentration, mean velocity and fluctuations. As the geometry of the duct is widened from $AR=1$ to 3, the secondary vortex on the horizontal wall significantly expands in the spanwise direction, and although the kinetic energy of the secondary flow increases close to the corner, it decays towards the duct centreplane in the $AR=3$ case so as the turbulent carrier phase approaches the behavior in spanwise-periodic channel flows, a fact that significantly affects the particle statistics. In the square duct the particle concentration in the viscous sublayer is maximum at the duct centreplane, whereas the maximum is found closer to the corner, at a distance of $|z/h|\approx1.25$ from the centreplane, in the $AR=3$ case. Interestingly the centreplane concentration in the rectangular duct is around 3 times lower than that in the square duct. Moreover, a second peak in the accumulation distribution is found right at the corners for both ducts. At this location the concentration increases with particle inertia. The secondary motion changes also the cross-stream map of the particle velocities significantly in comparison to the fluid flow statistics. These directly affect the particle velocity fluctuations such that multiple peaks appear near the duct walls for the particle streamwise and wall-normal velocity fluctuations.

preprint2016arXiv

Bypass transition and spot nucleation in boundary layers

The spatio-temporal aspects of the transition to turbulence are considered in the case of a boundary layer flow developing above a flat plate exposed to free-stream turbulence. Combining results on the receptivity to free-stream turbulence with the nonlinear concept of a transition threshold, a physically motivated model suggests a spatial distribution of spot nucleation events. To describe the evolution of turbulent spots a probabilistic cellular automaton is introduced, with all parameters directly fitted from numerical simulations of the boundary layer. The nucleation rates are then combined with the cellular automaton model, yielding excellent quantitative agreement with the statistical characteristics for different free-stream turbulence levels. We thus show how the recent theoretical progress on transitional wall-bounded flows can be extended to the much wider class of spatially developing boundary-layer flows.

preprint2016arXiv

Edge states as mediators of bypass transition in boundary-layer flows

The concept of edge state is investigated in the asymptotic suction boundary layer in relation with the receptivity process to noisy perturbations and the nucleation of turbulent spots. Edge tracking is first performed numerically, without imposing any discrete symmetry, in a large computational domain allowing for full spatial localisation of the perturbation velocity. The edge state is a three-dimensional localised structure recurrently characterised by a single low-speed streak that experiences erratic bursts and planar shifts. This recurrent streaky structure is then compared with predecessors of individual spot nucleation events, triggered by non-localised initial noise. The present results suggest a nonlinear picture, rooted in dynamical systems theory, of the nucleation process of turbulent spots in boundary-layer flows, in which the localised edge state play the role of state-space mediator.

preprint2016arXiv

Particle transport in turbulent curved pipe flow

Direct numerical simulations (DNS) of particle-laden turbulent flow in straight, mildly curved and strongly bent pipes are performed in which the solid phase is modelled as small heavy spherical particles. A total of seven populations of dilute particles with different Stokes numbers, one-way coupled with their carrier phase, are simulated. The objective is to examine the effect of the curvature on micro-particle transport and accumulation. It is shown that even a slight non-zero curvature in the flow configuration strongly impact the particle concentration map such that the concentration of inertial particles with bulk Stokes number 0.45 (based on bulk velocity and pipe radius) at the inner-bend wall of mildly curved pipe becomes 12.8 times larger than that in the viscous sublayer of the straight pipe. Near-wall helicoidal particle streaks are observed in the curved configurations with their inclination varying with the strength of the secondary motion of the carrier phase. A reflection layer, as previously observed in particle laden turbulent S-shaped channels, is also apparent in the strongly curved pipe with heavy particles. In addition, depending on the curvature, the central regions of the mean Dean vortices appear to be completely depleted of particles, as observed also in the partially re-laminarised region at the inner bend. The turbophoretic drift of the particles is shown to be affected by weak and strong secondary motions of the carrier phase and geometry-induced centrifugal forces. The first and second-order moments of the velocity and acceleration of the particulate phase in the same configurations are addressed in a companion paper by the same authors. The current data-set will be useful for modelling particles advected in wall-bounded turbulent flows where the effects of the curvature are not negligible.

preprint2015arXiv

Turbulence transition in the asymptotic suction boundary layer

We study the transition to turbulence in the asymptotic suction boundary layer (ASBL) by direct numerical simulation. Tracking the motion of trajectories intermediate between laminar and turbulent states we can identify the invariant object inside the laminar-turbulent boundary, the edge state. In small domains, the flow behaves like a travelling wave over short time intervals. On longer times one notes that the energy shows strong bursts at regular time intervals. During the bursts the streak structure is lost, but it reforms, translated in the spanwise direction by half the domain size. Varying the suction velocity allows to embed the flow into a family of flows that interpolate between plane Couette flow and the ASBL. Near the plane Couette limit, the edge state is a travelling wave. Increasing the suction, the travelling wave and a symmetry-related copy of it undergo a saddle-node infinite-period (SNIPER) bifurcation that leads to bursting and discrete-symmetry shifts. In wider domains, the structures localize in the spanwise direction, and the flow in the active region is similar to the one in small domains. There are still periodic bursts at which the flow structures are shifted, but the shift-distance is no longer connected to a discrete symmetry of the flow geometry. Two different states are found by edge tracking techniques, one where structures are shifted to the same side at every burst and one where they are alternatingly shifted to the left and to the right.

preprint2014arXiv

Recurrent bursts via linear processes in turbulent environments

Large-scale instabilities occurring in the presence of small-scale turbulent fluctuations are frequently observed in geophysical or astrophysical contexts but are difficult to reproduce in the laboratory. Using extensive numerical simulations, we report here on intense recurrent bursts of turbulence in plane Poiseuille flow rotating about a spanwise axis. A simple model based on the linear instability of the mean flow can predict the structure and time scale of the nearly-periodic and self-sustained burst cycles. Rotating Poiseuille flow is suggested as a prototype for future studies of low-dimensional dynamics embedded in strongly turbulent environments.

preprint2013arXiv

Complexity of localised coherent structures in a boundary-layer flow

We study numerically transitional coherent structures in a boundary-layer flow with homogeneous suction at the wall (the so-called asymptotic suction boundary layer ASBL). The dynamics restricted to the laminar-turbulent separatrix is investigated in a spanwise-extended domain that allows for robust localisation of all edge states. We work at fixed Reynolds number and study the edge states as a function of the streamwise period. We demonstrate the complex spatio-temporal dynamics of these localised states, which exhibits multistability and undergoes complex bifurcations leading from periodic to chaotic regimes. It is argued that in all regimes the dynamics restricted to the edge is essentially low-dimensional and non-extensive.

preprint2011arXiv

Hairpins et al. in Turbulent Boundary Layers

A new set of three-dimensional visualisations of a large-scale direct numerical simulations (DNS) of a turbulent boundary layer is presented. The Reynolds number ranges from $Re_θ=180$ to 4300, based on the momentum-loss thickness $θ$ and the free-stream velocity $U_\infty$. The focus of the present fluid dynamics video is on analysing the coherent vortical structures in the boundary layer: It is clearly shown that the initial phases are dominated by coherent so-called hairpin vortices which are characteristic remainders of the laminar-turbulent transition at lower Reynolds numbers. At higher $Re$ (say $Re_θ>2000$), these structures are no longer seen as being dominant; the coherence is clearly lost, both in the near-wall region as well as in the outer layer of the boundary layer. Note, however, that large-scale streaks in the streamwise velocity, which have their peak energy at about half the boundary-layer thickness, are unambiguously observed. In addition to visualisation with classical three-dimensional isosurfaces, the video is also rendered using stereoscopic views using red-cyan anaglyphs.

preprint2010arXiv

Stability of a jet in crossflow

We have produced a fluid dynamics video with data from Direct Numerical Simulation (DNS) of a jet in crossflow at several low values of the velocity inflow ratio R. We show that, as the velocity ratio R increases, the flow evolves from simple periodic vortex shedding (a limit cycle) to more complicated quasi-periodic behavior, before finally exhibiting asymmetric chaotic motion. We also perform a stability analysis just above the first bifurcation, where R is the bifurcation parameter. Using the overlap of the direct and the adjoint eigenmodes, we confirm that the first instability arises in the shear layer downstream of the jet orifice on the boundary of the backflow region just behind the jet.

preprint2010arXiv

The structure of a turbulent boundary layer studied by numerical simulation

Time-dependent visualisations of large-scale direct and large-eddy simulations (DNS and LES) of a turbulent boundary layer reaching up to $Re_θ=4300$ are presented. The focus of the present fluid dynamics video is on analysing the coherent vortical structures in the boundary layer: It is clearly shown that hairpin vortices are indeed present in the simulation data, however they are characteristic remainders of the laminar-turbulent transition at lower Reynolds numbers. At higher $Re$ (say $Re_θ>2000$), these structures are no longer seen as being dominant; the coherence is clearly lost, both in the near-wall region as well as in the outer layer of the boundary layer. Note, however, that large-scale streaks in the streamwise velocity, which have their peak energy at about half the boundary-layer thickness, are unambiguously observed.