Showing 1–46 of 46 results for author: Lozano-Durán, A

Machine-learning wall-model large-eddy simulation accounting for isotropic roughness under local equilibrium

Authors: Rong Ma, Adrian Lozano-Duran

Abstract: We introduce a wall model (WM) for large-eddy simulation (LES) applicable to rough surfaces with Gaussian and non-Gaussian distributions for both transitionally and fully rough regimes. The model is applicable to arbitrary complex geometries where roughness elements are assumed to be underresolved. The wall model is implemented using a feedforward neural network, with the geometric properties of t… ▽ More We introduce a wall model (WM) for large-eddy simulation (LES) applicable to rough surfaces with Gaussian and non-Gaussian distributions for both transitionally and fully rough regimes. The model is applicable to arbitrary complex geometries where roughness elements are assumed to be underresolved. The wall model is implemented using a feedforward neural network, with the geometric properties of the roughness topology and near-wall flow quantities serving as input. The optimal set of non-dimensional input features is identified using information theory, selecting variables that maximize information about the output while minimizing redundancy among inputs. The model incorporates a confidence score based on Gaussian process modeling, enabling the detection of low model performance for unseen rough surfaces. The model is trained using a direct numerical simulation roughness database comprising approximately 200 cases. The roughness geometries for the database are selected from a large repository through active learning. This approach ensures that the rough surfaces incorporated into the database are the most informative. The model performance is evaluated both a-priori and a-posteriori in WMLES of turbulent channel flows with rough walls. Over 120 channel flow cases are considered, including untrained roughness geometries, roughness Reynolds numbers, and grid resolutions for both transitionally- and fully-rough regimes. The results show that the rough-wall model predicts the wall shear stress within 15% accuracy. The model is also assessed on a high-pressure turbine blade with two different rough surfaces. The WM predicts the skin friction and the mean velocity deficit within 10% accuracy except the region with shock waves. △ Less

Submitted 1 June, 2024; originally announced June 2024.

Decomposing causality into its synergistic, unique, and redundant components

Authors: Álvaro Martínez-Sánchez, Gonzalo Arranz, Adrián Lozano-Durán

Abstract: Causality lies at the heart of scientific inquiry, serving as the fundamental basis for understanding interactions among variables in physical systems. Despite its central role, current methods for causal inference face significant challenges due to nonlinear dependencies, stochastic interactions, self-causation, collider effects, and influences from exogenous factors, among others. While existing… ▽ More Causality lies at the heart of scientific inquiry, serving as the fundamental basis for understanding interactions among variables in physical systems. Despite its central role, current methods for causal inference face significant challenges due to nonlinear dependencies, stochastic interactions, self-causation, collider effects, and influences from exogenous factors, among others. While existing methods can effectively address some of these challenges, no single approach has successfully integrated all these aspects. Here, we address these challenges with SURD: Synergistic-Unique-Redundant Decomposition of causality. SURD quantifies causality as the increments of redundant, unique, and synergistic information gained about future events from past observations. The formulation is non-intrusive and applicable to both computational and experimental investigations, even when samples are scarce. We benchmark SURD in scenarios that pose significant challenges for causal inference and demonstrate that it offers a more reliable quantification of causality compared to previous methods. △ Less

Submitted 21 August, 2024; v1 submitted 20 May, 2024; originally announced May 2024.

Comments: arXiv admin note: text overlap with arXiv:2310.20544

Building-block flow model for computational fluids

Authors: Gonzalo Arranz, Yuenong Ling, Sam Costa, Konrad Goc, Adrian Lozano-Duran

Abstract: We introduce a closure model for wall-modeled large-eddy simulation (WMLES), referred to as the Building-block Flow Model (BFM). The foundation of the model rests on the premise that a finite collection of simple flows encapsulates the essential physics necessary to predict more complex scenarios. The BFM is implemented using artificial neural networks and introduces five advancements within the f… ▽ More We introduce a closure model for wall-modeled large-eddy simulation (WMLES), referred to as the Building-block Flow Model (BFM). The foundation of the model rests on the premise that a finite collection of simple flows encapsulates the essential physics necessary to predict more complex scenarios. The BFM is implemented using artificial neural networks and introduces five advancements within the framework of WMLES: (1) It is designed to predict multiple flow regimes (wall turbulence under zero, favorable, adverse mean-pressure-gradient, and separation); (2) It unifies the closure model at solid boundaries (i.e., the wall model) and the rest of the flow (i.e., the subgrid-scale model) into a single entity; (3) It ensures consistency with numerical schemes and gridding strategy by accounting for numerical errors; (4) It is directly applicable to arbitrary complex geometries; (5) It can be scaled up to model additional flow physics in the future if needed (e.g., shockwaves and laminar-to-turbulent transition). The BFM is utilized to predict key quantities of interest in turbulent channel and pipe flows, a Gaussian bump, a simplified aircraft, and a realistic aircraft in landing configuration. In all cases, the BFM demonstrates similar or superior capabilities in terms of accuracy and computational efficiency compared to previous state-of-the-art closure models. △ Less

Submitted 14 June, 2024; v1 submitted 13 March, 2024; originally announced March 2024.

Informative and non-informative decomposition of turbulent flow fields

Authors: Gonzalo Arranz, Adrian Lozano-Durán

Abstract: Not all the information in a turbulent field is relevant for understanding particular regions or components of the flow. Here, we present a method for decomposing a source field into its informative $\mathbfΦ_{I}(\mathbf{x},t)$ and residual $\mathbfΦ_R(\mathbf{x},t)$ components relative to another target field. The method is referred to as informative and non-informative decomposition (IND). All t… ▽ More Not all the information in a turbulent field is relevant for understanding particular regions or components of the flow. Here, we present a method for decomposing a source field into its informative $\mathbfΦ_{I}(\mathbf{x},t)$ and residual $\mathbfΦ_R(\mathbf{x},t)$ components relative to another target field. The method is referred to as informative and non-informative decomposition (IND). All the necessary information for physical understanding, reduced-order modeling, and control of the target variable is contained in $\mathbfΦ_{I}(\mathbf{x},t)$ whereas $\mathbfΦ_R(\mathbf{x},t)$ offers no substantial utility in these contexts. The decomposition is formulated as an optimization problem that seeks to maximize the time lagged mutual information of the source variable with the target variable while minimizing the mutual information of its non-informative component. The method is applied to extract the informative and non-informative components of the velocity field in a turbulent channel flow, using the wall shear stress as the target variable. We demonstrate the utility of IND in three scenarios: (i) understanding the effect of the velocity fluctuations on the wall-shear stress, (ii) predicting the wall-shear stress using velocities far from the wall, and (iii) developing control strategies for drag reduction in opposition control. △ Less

Submitted 17 February, 2024; originally announced February 2024.

Are the dynamics of wall turbulence in minimal channels and larger domain channels equivalent? A graph-theoretic approach

Authors: Ahmed Elnahhas, Emma Lenz, Parviz Moin, Adrián Lozano-Durán, H. Jane Bae

Abstract: This work proposes two algorithmic approaches to extract critical dynamical mechanisms in wall-bounded turbulence with minimum human bias. In both approaches, multiple types of coherent structures are spatiotemporally tracked, resulting in a complex multilayer network. Network motif analysis, i.e., extracting dominant non-random elemental patterns within these networks, is used to identify the mos… ▽ More This work proposes two algorithmic approaches to extract critical dynamical mechanisms in wall-bounded turbulence with minimum human bias. In both approaches, multiple types of coherent structures are spatiotemporally tracked, resulting in a complex multilayer network. Network motif analysis, i.e., extracting dominant non-random elemental patterns within these networks, is used to identify the most dominant dynamical mechanisms. Both approaches, combined with network motif analysis, are used to answer whether the main dynamical mechanisms of a minimal flow unit (MFU) and a larger unconstrained channel flow, labeled a full channel (FC), at $Re_τ\approx 180$, are equivalent. The first approach tracks traditional coherent structures defined as low- and high-speed streaks, ejections, and sweeps. It is found that the roll-streak pairing, consistent with the current understanding of self-sustaining processes, is the most significant and simplest dynamical mechanism in both flows. However, the MFU has a timescale for this mechanism that is approximately $2.83$ times slower than that of the FC. In the second approach, we use semi-Lagrangian wavepackets and define coherent structures from their energetic streak, roll, and small-scale phase space. This method also shows similar motifs for both the MFU and FC. It indicates that, on average, the most dominant phase-space motifs are similar between the two flows, with the significant events taking place approximately $2.21$ times slower in the MFU than in the FC. This value is more consistent with the implied timescale ratio of only the slow speed streaks taking part in the roll-streak pairing extracted using the first multi-type spatiotemporal approach, which is approximately $2.17$ slower in the MFU than the FC. △ Less

Submitted 15 January, 2024; originally announced January 2024.

Information-theoretic causality and applications to turbulence: energy cascade and inner/outer layer interactions

Authors: Adrián Lozano-Durán, Gonzalo Arranz, Yuenong Ling

Abstract: We introduce an information-theoretic method for quantifying causality in chaotic systems. The approach, referred to as IT-causality, quantifies causality by measuring the information gained about future events conditioned on the knowledge of past events. The causal interactions are classified into redundant, unique, and synergistic contributions depending on their nature. The formulation is non-i… ▽ More We introduce an information-theoretic method for quantifying causality in chaotic systems. The approach, referred to as IT-causality, quantifies causality by measuring the information gained about future events conditioned on the knowledge of past events. The causal interactions are classified into redundant, unique, and synergistic contributions depending on their nature. The formulation is non-intrusive, invariance under invertible transformations of the variables, and provides the missing causality due to unobserved variables. The method only requires pairs of past-future events of the quantities of interest, making it convenient for both computational simulations and experimental investigations. IT-causality is validated in four scenarios representing basic causal interactions among variables: mediator, confounder, redundant collider, and synergistic collider. The approach is leveraged to address two questions relevant to turbulence research: i) the scale locality of the energy cascade in isotropic turbulence, and ii) the interactions between inner and outer layer flow motions in wall-bounded turbulence. In the former case, we demonstrate that causality in the energy cascade flows sequentially from larger to smaller scales without requiring intermediate scales. Conversely, the flow of information from small to large scales is shown to be redundant. In the second problem, we observe a unidirectional causality flow, with causality predominantly originating from the outer layer and propagating towards the inner layer, but not vice versa. The decomposition of IT-causality into intensities also reveals that the causality is primarily associated with high-velocity streaks. △ Less

Submitted 31 October, 2023; originally announced October 2023.

Near-Field Wall-Modeled Large-Eddy Simulation of the NASA X-59 Low-Boom Flight Demonstrator

Authors: Emily Williams, Gonzalo Arranz, Adrián Lozano-Durán

Abstract: Wall-modeled large-eddy simulation (WMLES) is utilized to analyze the experimental aircraft X-59 Quiet SuperSonic Technology (QueSST) developed by Lockheed Martin at Skunk Works for NASA's Low-Boom Flight Demonstrator project. The simulations utilize the charLES solver and aim to assess the ability of WMLES to predict near-field noise levels under cruise conditions, considering various subgrid-sca… ▽ More Wall-modeled large-eddy simulation (WMLES) is utilized to analyze the experimental aircraft X-59 Quiet SuperSonic Technology (QueSST) developed by Lockheed Martin at Skunk Works for NASA's Low-Boom Flight Demonstrator project. The simulations utilize the charLES solver and aim to assess the ability of WMLES to predict near-field noise levels under cruise conditions, considering various subgrid-scale (SGS) models and grid resolutions. The results are compared with previous numerical studies based on the Reynolds-averaged Navier-Stokes (RANS) equations. Our findings demonstrate that WMLES produces near-field pressure predictions that are similar to those of RANS simulations at a comparable computational cost. Some mild discrepancies are observed between the WMLES and RANS predictions downstream the aircraft. These differences persist for finest grid refinement considered, suggesting that they might be attributed to underresolved interactions of shock waves and expansions waves at the trailing edge. △ Less

Submitted 5 July, 2023; originally announced July 2023.

The coherent structure of the energy cascade in isotropic turbulence

Authors: Danah Park, Adrian Lozano-Duran

Abstract: The energy cascade, i.e. the transfer of kinetic energy from large-scale to small-scale flow motions, has been the cornerstone of turbulence theories and models since the 1940s. However, understanding the spatial organization of the energy transfer has remained elusive. In this study, we utilize numerical data of isotropic turbulence to investigate the three-dimensional spatial structure of the en… ▽ More The energy cascade, i.e. the transfer of kinetic energy from large-scale to small-scale flow motions, has been the cornerstone of turbulence theories and models since the 1940s. However, understanding the spatial organization of the energy transfer has remained elusive. In this study, we utilize numerical data of isotropic turbulence to investigate the three-dimensional spatial structure of the energy cascade. Specifically, we focus on analyzing the averaged flow patterns associated with intense energy transfer events in the inertial range. Our findings indicate that forward energy-transfer events are predominantly confined in the high strain-rate region created between two distinct zones of elevated enstrophy. On average, these zones manifest in the form of two hairpin-like shapes with opposing orientations. The mean velocity field associated with the energy transfer exhibits a saddle point topology when observed in the frame of reference local to the event. The analysis also shows that the primary driving mechanism for the cascade involves strain-rate self-amplification, which is responsible for 85% of the energy transfer, whereas vortex stretching accounts for less than 15%. △ Less

Submitted 15 June, 2024; v1 submitted 14 June, 2023; originally announced June 2023.

Comments: 16 pages, 12 figures

Wall-modeled large-eddy simulation based on building-block flows

Authors: Yuenong Ling, Gonzalo Arranz, Emily Williams, Konrad Goc, Kevin Griffin, Adrián Lozano-Durán

Abstract: A unified subgrid-scale (SGS) and wall model for large-eddy simulation (LES) is proposed by devising the flow as a collection of building blocks that enables the prediction of the eddy viscosity. The core assumption of the model is that simple canonical flows contain the essential physics to provide accurate predictions of the SGS tensor in more complex flows. The model is constructed to predict z… ▽ More A unified subgrid-scale (SGS) and wall model for large-eddy simulation (LES) is proposed by devising the flow as a collection of building blocks that enables the prediction of the eddy viscosity. The core assumption of the model is that simple canonical flows contain the essential physics to provide accurate predictions of the SGS tensor in more complex flows. The model is constructed to predict zero-pressure-gradient wall-bounded turbulence, adverse pressure gradient effects, separation and laminar flow. The approach is implemented using a Bayesian classifier, which identifies the contribution of each building block in the flow, and a neural-network-based predictor, which estimates the eddy viscosity based on the building-block units. The training data are directly obtained from wall-modeled LES with an exact SGS/wall model for the mean quantities to guarantee consistency with the numerical discretization. The model is validated in canonical flows and the NASA High-Lift Common Research Model and shown to improve the predictions with respect to current modeling approaches. △ Less

Submitted 17 December, 2022; v1 submitted 9 December, 2022; originally announced December 2022.

Machine learning building-block-flow wall model for large-eddy simulation

Authors: Adrián Lozano-Durán, H. Jane Bae

Abstract: A wall model for large-eddy simulation (LES) is proposed by devising the flow as a combination of building blocks. The core assumption of the model is that a finite set of simple canonical flows contains the essential physics to predict the wall-shear stress in more complex scenarios. The model is constructed to predict zero/favourable/adverse mean pressure gradient wall turbulence, separation, st… ▽ More A wall model for large-eddy simulation (LES) is proposed by devising the flow as a combination of building blocks. The core assumption of the model is that a finite set of simple canonical flows contains the essential physics to predict the wall-shear stress in more complex scenarios. The model is constructed to predict zero/favourable/adverse mean pressure gradient wall turbulence, separation, statistically unsteady turbulence with mean flow three-dimensionality, and laminar flow. The approach is implemented using two types of artificial neural networks: a classifier, which identifies the contribution of each building block in the flow, and a predictor, which estimates the wall-shear stress via combination of the building-block flows. The training data are directly obtained from wall-modelled LES (WMLES) optimised to reproduce the correct mean quantities. This approach guarantees the consistency of the training data with the numerical discretisation and the gridding strategy of the flow solver. The output of the model is accompanied by a confidence score in the prediction that aids the detection of regions where the model underperforms. The model is validated in canonical flows (e.g. laminar/turbulent boundary layers, turbulent channels, turbulent Poiseuille-Couette flow, turbulent pipe) and two realistic aircraft configurations: the NASA Common Research Model High-lift and NASA Juncture Flow experiment. It is shown that the building-block-flow wall model outperforms (or matches) the predictions by an equilibrium wall model. It is also concluded that further improvements in WMLES should incorporate advances in subgrid-scale modelling to minimise error propagation to the wall model. △ Less

Submitted 30 April, 2023; v1 submitted 14 November, 2022; originally announced November 2022.

Causality analysis of large-scale structures in the flow around a wall-mounted square cylinder

Authors: Álvaro Martínez-Sánchez, Esteban López, Soledad Le Clainche, Adrián Lozano-Durán, Ankit Srivastava, Ricardo Vinuesa

Abstract: The aim of this work is to analyse the formation mechanisms of large-scale coherent structures in the flow around a wall-mounted square cylinder, due to their impact on pollutant transport within cities. To this end, we assess causal relations between the modes of a reduced-order model obtained by applying proper-orthogonal decomposition to high-fidelity-simulation data of the flow case under stud… ▽ More The aim of this work is to analyse the formation mechanisms of large-scale coherent structures in the flow around a wall-mounted square cylinder, due to their impact on pollutant transport within cities. To this end, we assess causal relations between the modes of a reduced-order model obtained by applying proper-orthogonal decomposition to high-fidelity-simulation data of the flow case under study. The causal relations are identified using conditional transfer entropy, which is an information-theoretical quantity that estimates the amount of information contained in the past of one variable about another. This allows for an understanding of the origins and evolution of different phenomena in the flow, with the aim of identifying the modes responsible for the formation of the main vortical structures. Our approach unveils that vortex-breaker modes are the most causal modes, in particular, over higher-order modes, and no significant causal relationships were found for vortex-generator modes. We validate this technique by determining the causal relations present in the nine-equation model of near-wall turbulence developed by Moehlis et al. (New J. Phys, vol. 6, 2004, p. 56), which are in good agreement with literature results for turbulent channel flows. △ Less

Submitted 30 September, 2022; originally announced September 2022.

Comments: 19 pages, 10 figures

Numerical and modeling error assessment of large-eddy simulation using direct-numerical-simulation-aided large-eddy simulation

Authors: H. Jane Bae, Adrian Lozano-Duran

Abstract: We study the numerical errors of large-eddy simulation (LES) in isotropic and wall-bounded turbulence. A direct-numerical-simulation (DNS)-aided LES formulation, where the subgrid-scale (SGS) term of the LES is computed by using filtered DNS data is introduced. We first verify that this formulation has zero error in the absence of commutation error between the filter and the differentiation operat… ▽ More We study the numerical errors of large-eddy simulation (LES) in isotropic and wall-bounded turbulence. A direct-numerical-simulation (DNS)-aided LES formulation, where the subgrid-scale (SGS) term of the LES is computed by using filtered DNS data is introduced. We first verify that this formulation has zero error in the absence of commutation error between the filter and the differentiation operator of the numerical algorithm. This method allows the evaluation of the time evolution of numerical errors for various numerical schemes at grid resolutions relevant to LES. The analysis shows that the numerical errors are of the same order of magnitude as the modeling errors and often cancel each other. This supports the idea that supervised machine learning algorithms trained on filtered DNS data might not be suitable for robust SGS model development, as this approach disregards the existence of numerical errors in the system that accumulates over time. The assessment of errors in turbulent channel flow also identifies that numerical errors close to the wall dominate, which has implications for the development of wall models. △ Less

Submitted 3 August, 2022; originally announced August 2022.

A Database for Reduced-Complexity Modeling of Fluid Flows

Authors: Aaron Towne, Scott T. M. Dawson, Guillaume A. Brès, Adrián Lozano-Durán, Theresa Saxton-Fox, Aadhy Parthasarathy, Anya R. Jones, Hulya Biler, Chi-An Yeh, Het D. Patel, Kunihiko Taira

Abstract: We present a publicly accessible database designed to aid in the conception, training, demonstration, evaluation, and comparison of reduced-complexity models for fluid mechanics. Availability of high-quality flow data is essential for all of these aspects of model development for both data-driven and physics-based methods. The database contains time-resolved data for six distinct datasets: a large… ▽ More We present a publicly accessible database designed to aid in the conception, training, demonstration, evaluation, and comparison of reduced-complexity models for fluid mechanics. Availability of high-quality flow data is essential for all of these aspects of model development for both data-driven and physics-based methods. The database contains time-resolved data for six distinct datasets: a large eddy simulation of a turbulent jet, direct numerical simulations of a zero-pressure-gradient turbulent boundary layer, particle-image-velocimetry measurements for the same boundary layer at several Reynolds numbers, direct numerical simulations of laminar stationary and pitching flat-plate airfoils, particle-image-velocimetry and force measurements of an airfoil encountering a gust, and a large eddy simulation of the separated, turbulent flow over an airfoil. These six cases span several key flow categories: laminar and turbulent, statistically stationary and transient, tonal and broadband spectral content, canonical and application-oriented, wall-bounded and free-shear flow, and simulation and experimental measurements. For each dataset, we describe the flow setup and computational/experimental methods, catalog the data available in the database, and provide examples of how these data can be used for reduced-complexity modeling. All data can be downloaded using a browser interface or Globus. Our vision is that the common testbed provided by this database will aid the fluid mechanics community in clarifying the distinct capabilities of new and existing methods. △ Less

Submitted 23 June, 2022; originally announced June 2022.

Information-theoretic formulation of dynamical systems: causality, modeling, and control

Authors: Adrián Lozano-Durán, Gonzalo Arranz

Abstract: The problems of causality, modeling, and control for chaotic, high-dimensional dynamical systems are formulated in the language of information theory. The central quantity of interest is the Shannon entropy, which measures the amount of information in the states of the system. Within this framework, causality is quantified by the information flux among the variables of interest in the dynamical sy… ▽ More The problems of causality, modeling, and control for chaotic, high-dimensional dynamical systems are formulated in the language of information theory. The central quantity of interest is the Shannon entropy, which measures the amount of information in the states of the system. Within this framework, causality is quantified by the information flux among the variables of interest in the dynamical system. Reduced-order modeling is posed as a problem related to the conservation of information in which models aim at preserving the maximum amount of relevant information from the original system. Similarly, control theory is cast in information-theoretic terms by envisioning the tandem sensor-actuator as a device reducing the unknown information of the state to be controlled. The new formulation is used to address three problems about the causality, modeling, and control of turbulence, which stands as a primary example of a chaotic, high-dimensional dynamical system. The applications include the causality of the energy transfer in the turbulent cascade, subgrid-scale modeling for large-eddy simulation, and flow control for drag reduction in wall-bounded turbulence. △ Less

Submitted 31 May, 2022; v1 submitted 17 November, 2021; originally announced November 2021.

POD-based study of turbulent plane Poiseuille flow: comparing structure and dynamics between quasi-linear simulations and DNS

Authors: Marios-Andreas Nikolaidis, Petros J. Ioannou, Brian F. Farrell, Adrián Lozano-Durán

Abstract: Turbulence in the restricted nonlinear (RNL) dynamics is analyzed and compared with DNS of Poiseuille turbulence at $R=1650$. The structures are obtained by POD analysis of the two components of the flow partition used in RNL dynamics: the streamwise-mean flow and fluctuations. POD analysis of the streamwise-mean flow indicates that the dominant POD modes, in both DNS and RNL, are roll-streaks har… ▽ More Turbulence in the restricted nonlinear (RNL) dynamics is analyzed and compared with DNS of Poiseuille turbulence at $R=1650$. The structures are obtained by POD analysis of the two components of the flow partition used in RNL dynamics: the streamwise-mean flow and fluctuations. POD analysis of the streamwise-mean flow indicates that the dominant POD modes, in both DNS and RNL, are roll-streaks harmonic in the spanwise. However, we conclude that these POD modes do not occur in isolation but rather are Fourier components of a coherent roll-streak structure. POD analysis of the fluctuations in DNS and RNL reveals similar complex structures consisting in part of oblique waves collocated with the streak. The origin of these structures is identified by their correspondence to POD modes predicted using a stochastic turbulence model (STM). These predicted POD modes are dominated by the optimally growing structures on the streak, which the STM predicts correctly to be of sinuous oblique wave structure. This close correspondence between the roll-streak structure and the associated fluctuations in DNS, RNL and the STM implies that the self-sustaining mechanism operating in DNS is essentially the same as that in RNL, which has been previously associated with optimal perturbation growth on the streak. △ Less

Submitted 29 March, 2023; v1 submitted 6 September, 2021; originally announced September 2021.

Comments: 31 pages, accepted by the Journal of Fluid Mechanics

Performance of wall-modeled LES with boundary-layer-conforming grids for external aerodynamics

Authors: Adrián Lozano-Durán, Sanjeeb T. Bose, Parviz Moin

Abstract: We investigate the error scaling and computational cost of wall-modeled large-eddy simulation (WMLES) for external aerodynamic applications. The NASA Juncture Flow is used as representative of an aircraft with trailing-edge smooth-body separation. Two gridding strategies are examined: i) constant-size grid, in which the near-wall grid size has a constant value and ii) boundary-layer-conforming gri… ▽ More We investigate the error scaling and computational cost of wall-modeled large-eddy simulation (WMLES) for external aerodynamic applications. The NASA Juncture Flow is used as representative of an aircraft with trailing-edge smooth-body separation. Two gridding strategies are examined: i) constant-size grid, in which the near-wall grid size has a constant value and ii) boundary-layer-conforming grid (BL-conforming grid), in which the grid size varies to accommodate the growth of the boundary-layer thickness. Our results are accompanied by a theoretical analysis of the cost and expected error scaling for the mean pressure coefficient ($C_p$) and mean velocity profiles. The prediction of $C_p$ is within less than $5\%$ error for all the grids studied, even when the boundary layers are marginally resolved. The high accuracy in the prediction of $C_p$ is attributed to the outer-layer nature of the mean pressure in attached flows. The errors in the predicted mean velocity profiles exhibit a large variability depending on the location considered, namely, fuselage, wing-body juncture, or separated trailing-edge. WMLES performs as expected in regions where the flow resembles a zero-pressure-gradient turbulent boundary layer such as the fuselage ($<5\%$ error). However, there is a decline in accuracy of WMLES predictions of mean velocities in the vicinity of wing-body junctions and, more acutely, in separated zones. The impact of the propagation of errors from the underresolved wing leading-edge is also investigated. It is shown that BL-conforming grids enable a higher accuracy in wing-body junctions and separated regions due to the more effective distribution of grid points, which in turn diminishes the streamwise propagation of errors. △ Less

Submitted 3 July, 2021; originally announced July 2021.

Comments: arXiv admin note: text overlap with arXiv:2101.00331

Toward a flow-structure-based wall-modeled large-eddy simulation paradigm

Authors: Ahmed Elnahhas, Adrián Lozano-Durán, Parviz Moin

Abstract: A promising and cost-effective method for numerical simulation of high Re wall-bounded flows is wall-modeled large-eddy simulation. Most wall models are formulated from the Reynolds-averaged Navier-Stokes equations (RANS). These RANS-based wall models are calibrated using mean turbulence data and make no use of the current vast knowledge on turbulent flow structure. Moreover, RANS-based models are… ▽ More A promising and cost-effective method for numerical simulation of high Re wall-bounded flows is wall-modeled large-eddy simulation. Most wall models are formulated from the Reynolds-averaged Navier-Stokes equations (RANS). These RANS-based wall models are calibrated using mean turbulence data and make no use of the current vast knowledge on turbulent flow structure. Moreover, RANS-based models are limited to predicting the mean velocity profile, and the mean wall shear stress. Using the knowledge of the wall-normal self-similarity of high-$Re_τ$ wall-bounded turbulent flows, we present a coupling between a near-wall patch of DNS resolution fixed in inner-units and an outer LES flow field. The near-wall patch captures the near-wall self-sustaining cycle as well as the lower portion of the self-similar hierarchy of eddies. Given that both the near-wall patch and the LES capture separate portions of the self-similar hierarchy of eddies an instantaneous top boundary condition for the patch is formulated. The near-wall model is capable of predicting subgrid-scale quantities such as the wall stress, velocity fluctuations, kinetic energy spectra, and flow structure across the entire near-wall layer. △ Less

Submitted 2 January, 2021; originally announced January 2021.

Journal ref: Annual Research Briefs 2020, Center for Turbulence Research, Stanford University

Performance of wall-modeled LES for external aerodynamics in the NASA Juncture Flow

Authors: Adrián Lozano-Durán, Sanjeeb T. Bose, Parviz Moin

Abstract: We investigate the performance of wall-modeled LES for external aerodynamics in the NASA Juncture Flow. We characterize the errors in the prediction of mean velocity profiles and pressure coefficient for three different locations over the aircraft: the upstream region of the fuselage, the wing-body juncture, and the wing-body juncture close to the trailing-edge. The last two locations are characte… ▽ More We investigate the performance of wall-modeled LES for external aerodynamics in the NASA Juncture Flow. We characterize the errors in the prediction of mean velocity profiles and pressure coefficient for three different locations over the aircraft: the upstream region of the fuselage, the wing-body juncture, and the wing-body juncture close to the trailing-edge. The last two locations are characterized by strong mean-flow three-dimensionality and separation. The message conveyed by our error analysis is that WMLES performs as expected in regions where the flow resembles a zero-pressure-gradient flat plate boundary layer. However, there is a clear decline of the current models in the presence of wing-body junctions and, more acutely, in separated zones. The slow convergence to the solution in these regions renders the brute-force grid-refinement approach to improve the accuracy of the solution unfeasible. The results reported above pertain to the mean velocity profile predicted using the typical grid resolution for external aerodynamics applications, i.e., 5--20 points per boundary-layer thickness. The prediction of the pressure coefficient is below 5\% error for all grid sizes considered, even when boundary layers were marginally resolved. The latter accuracy can be attributed to the outer-layer nature of the mean pressure, which becomes less sensitive to flow details within the turbulent boundary layer. Finally, we show that boundary-layer-conforming grids (i.e., grids maintaining a constant number of points per boundary-layer thickness) allow for a more efficient distribution of grid points and smaller errors. Our results suggest that novel modeling venues encompassing physical insights, together with numerical and gridding advancements, must be exploited to attain predictions within the tolerance required for Certification by Analysis. △ Less

Submitted 1 January, 2021; originally announced January 2021.

Journal ref: Annual Research Briefs, Center for Turbulence Research (Stanford University), 2020

Requirements and sensitivity analysis of RANS-free wall-modeled LES

Authors: Michael Whitmore, Adrián Lozano-Durán, Parviz Moin

Abstract: We study the sensitivity of wall model input variables to the modeling choices of the outer LES. This work is motivated by sensitivities observed in dynamic slip wall models. These dynamic wall models use variables from the near-wall LES solution as inputs to predict the wall stress without relying on a priori coefficients or equilibrium assumption. Mitigating the sensitivities in the wall model i… ▽ More We study the sensitivity of wall model input variables to the modeling choices of the outer LES. This work is motivated by sensitivities observed in dynamic slip wall models. These dynamic wall models use variables from the near-wall LES solution as inputs to predict the wall stress without relying on a priori coefficients or equilibrium assumption. Mitigating the sensitivities in the wall model inputs allows development of robust dynamic wall models. The effects of SGS model, boundary condition type, numerics, and mesh topology are assessed through a series of WMLES calculations of turbulent channels. Probability density functions are computed from planes of the WMLES solutions at a wall-normal sampling height and are used as a metric for sensitivity. Sensitivity to SGS model is alleviated when the fraction of total wall stress carried by the SGS model is held constant. Use of hexagonal close-packed grids mitigated numerical sensitivities. △ Less

Submitted 29 December, 2020; originally announced December 2020.

Comments: 12 pages, 9 figures

Journal ref: Annual Research Briefs 2020, Center for Turbulence Research, Stanford University

Self-critical machine-learning wall-modeled LES for external aerodynamics

Authors: Adrián Lozano-Durán, Hyunji Jane Bae

Abstract: The prediction of aircraft aerodynamic quantities of interest remains among the most pressing challenges for computational fluid dynamics. The aircraft aerodynamics are inherently turbulent with mean-flow three-dimensionality, often accompanied by laminar-to-turbulent transition, flow separation, secondary flow motions at corners, and shock wave formation, to name a few. However, the most widespre… ▽ More The prediction of aircraft aerodynamic quantities of interest remains among the most pressing challenges for computational fluid dynamics. The aircraft aerodynamics are inherently turbulent with mean-flow three-dimensionality, often accompanied by laminar-to-turbulent transition, flow separation, secondary flow motions at corners, and shock wave formation, to name a few. However, the most widespread wall models are built upon the assumption of statistically-in-equilibrium wall-bounded turbulence and do not faithfully account for the wide variety of flow conditions described above. This raises the question of how to devise models capable of accounting for such a vast and rich collection of flow physics in a feasible manner. In this work, we propose tackling the wall-modeling challenge by devising the flow as a collection of building blocks, whose information enables the prediction of the stress as the wall. The model relies on the assumption that simple canonical flows contain the essential flow physics to devise accurate models. Three types of building block units were used to train the model: turbulent channel flows, turbulent ducts and turbulent boundary layers with separation. This limited training set will be extended in future versions of the model. The approach is implemented using two interconnected artificial neural networks: a classifier, which identifies the contribution of each building block in the flow; and a predictor, which estimates the wall stress via non-linear combinations of building-block units. The output of the model is accompanied by the confidence in the prediction. The latter value aids the detection of areas where the model underperforms, such as flow regions that are not representative of the building blocks used to train the model. The model is validated in a unseen case representative of external aerodynamic applications: the NASA Juncture Flow Experiment. △ Less

Submitted 22 May, 2021; v1 submitted 17 December, 2020; originally announced December 2020.

Comments: Minor typos corrected

Journal ref: Annual Research Brief (Center for Turbulence Research, Stanford University) 2020

Information transfer between turbulent boundary layer and porous media

Authors: Wenkang Wang, Xu Chu, Adrián Lozano-Durán, Rainer Helmig, Bernhard Weigand

Abstract: The interaction between the flow above and below a permeable wall is a central topic in the study of porous media. While previous investigations have provided compelling evidence of the strong coupling between the two regions, few studies have quantitatively measured the directionality, i.e., causal-and-effect relations, of this interaction. To shed light on the problem, we use transfer entropy as… ▽ More The interaction between the flow above and below a permeable wall is a central topic in the study of porous media. While previous investigations have provided compelling evidence of the strong coupling between the two regions, few studies have quantitatively measured the directionality, i.e., causal-and-effect relations, of this interaction. To shed light on the problem, we use transfer entropy as a marker to evaluate the causal interaction between the free turbulent flow and the porous media using interface-resolved direct numerical simulation. Our results show that the porosity of the porous medium has a profound impact on the intensity, time scale and spatial extent of surface-subsurface interactions. For values of porosity equal to 0.5, top-down and bottom-up interactions strongly are asymmetric, the former being mostly influenced by small near-wall eddies. As the porosity increases, both top-down and bottom-up interactions are dominated by shear-flow instabilities. △ Less

Submitted 2 November, 2020; originally announced November 2020.

Cause-and-effect of linear mechanisms sustaining wall turbulence

Authors: Adrián Lozano-Durán, Navid C. Constantinou, Marios-Andreas Nikolaidis, Michael Karp

Abstract: Despite the nonlinear nature of turbulence, there is evidence that part of the energy-transfer mechanisms sustaining wall turbulence can be ascribed to linear processes. The different scenarios stem from linear stability theory and comprise exponential instabilities, neutral modes, transient growth from non-normal operators, and parametric instabilities from temporal mean-flow variations, among ot… ▽ More Despite the nonlinear nature of turbulence, there is evidence that part of the energy-transfer mechanisms sustaining wall turbulence can be ascribed to linear processes. The different scenarios stem from linear stability theory and comprise exponential instabilities, neutral modes, transient growth from non-normal operators, and parametric instabilities from temporal mean-flow variations, among others. These mechanisms, each potentially capable of leading to the observed turbulence structure, are rooted in theoretical and conceptual arguments. Whether the flow follows any or a combination of them remains elusive. Here, we evaluate the linear mechanisms responsible for the energy transfer from the streamwise-averaged mean-flow ($\bf U$) to the fluctuating velocities ($\bf u'$). We use cause-and-effect analysis based on interventions. This is achieved by direct numerical simulation of turbulent channel flows at low Reynolds number, in which the energy transfer from $\bf U$ to $\bf u'$ is constrained to preclude a targeted linear mechanism. We show that transient growth is sufficient for sustaining realistic wall turbulence. Self-sustaining turbulence persists when exponential instabilities, neutral modes, and parametric instabilities of the mean flow are suppressed. We further show that a key component of transient growth is the Orr/push-over mechanism induced by spanwise variations of the base flow. Finally, we demonstrate that an ensemble of simulations with various frozen-in-time $\bf U$ arranged so that only transient growth is active, can faithfully represent the energy transfer from $\bf U$ to $\bf u'$ as in realistic turbulence. Our approach provides direct cause-and-effect evaluation of the linear energy-injection mechanisms from $\bf U$ to $\bf u'$ in the fully nonlinear system and simplifies the conceptual model of self-sustaining wall turbulence. △ Less

Submitted 7 October, 2020; v1 submitted 9 May, 2020; originally announced May 2020.

Journal ref: J. Fluid Mech., vol. 914, A8, 2021

On the evolution of the velocity gradient tensor in transitional boundary layers

Authors: Ahmed Elnahhas, Perry L. Johnson, Adrián Lozano-Durán, Parviz Moin

Abstract: We study the transition to turbulence from the perspective of the velocity gradient tensor dynamics. Our work is motivated by the observation of nonlinear structures emerging during transition, as revealed by vortex identifiers such as the Q-criterion. To that end, we have derived transport equations based on several invariants of the velocity gradient tensor to obtain integral budgets spanning th… ▽ More We study the transition to turbulence from the perspective of the velocity gradient tensor dynamics. Our work is motivated by the observation of nonlinear structures emerging during transition, as revealed by vortex identifiers such as the Q-criterion. To that end, we have derived transport equations based on several invariants of the velocity gradient tensor to obtain integral budgets spanning the different stages of transition from laminar to turbulent flow. We have also discussed which quantity would be the most appropriate in our study, while keeping in mind the potential for future modeling applications. △ Less

Submitted 19 January, 2020; originally announced January 2020.

Journal ref: Annual Research Briefs 2019, Center for Turbulence Research, Stanford University

Wall-modeled large-eddy simulation of non-equilibrium turbulent boundary layers

Authors: Minjeong Cho, George I. Park, Adrián Lozano-Durán, Parviz Moin

Abstract: We conducted WMLES to examine the performance of a simple and widely used ODE-based equilibrium wall model in a spatially-developing 3D TBL inside a bent square duct (Schwarz and Bradshaw 1994) and 3D separated flows behind a skewed bump (Ching et al. 2018a,b; Ching and Eaton 2019). From the square duct simulation, the mean velocity profiles and crossflow angles in the outer region were predicted… ▽ More We conducted WMLES to examine the performance of a simple and widely used ODE-based equilibrium wall model in a spatially-developing 3D TBL inside a bent square duct (Schwarz and Bradshaw 1994) and 3D separated flows behind a skewed bump (Ching et al. 2018a,b; Ching and Eaton 2019). From the square duct simulation, the mean velocity profiles and crossflow angles in the outer region were predicted with high accuracy for all the considered mesh resolutions. Some disagreement was observed in the crossflow angles in the bend region where the non-equilibrium effect is most significant. Also, the simulation for the wall-mounted skewed bump showed that this simple ODE-based equilibrium wall model along with an adequate grid resolution around the 3D separation point resulted in reasonable predictions of 3D separating and reattaching flows, including mean velocity distributions, separation bubbles, and vortex structures in the bump wake. △ Less

Submitted 3 January, 2020; originally announced January 2020.

Journal ref: Annual Research Briefs 2019, Center for Turbulence Research, Stanford University

Computational model discovery with reinforcement learning

Authors: Maxime Bassenne, Adrián Lozano-Durán

Abstract: The motivation of this study is to leverage recent breakthroughs in artificial intelligence research to unlock novel solutions to important scientific problems encountered in computational science. To address the human intelligence limitations in discovering reduced-order models, we propose to supplement human thinking with artificial intelligence. Our three-pronged strategy consists of learning (… ▽ More The motivation of this study is to leverage recent breakthroughs in artificial intelligence research to unlock novel solutions to important scientific problems encountered in computational science. To address the human intelligence limitations in discovering reduced-order models, we propose to supplement human thinking with artificial intelligence. Our three-pronged strategy consists of learning (i) models expressed in analytical form, (ii) which are evaluated a posteriori, and iii) using exclusively integral quantities from the reference solution as prior knowledge. In point (i), we pursue interpretable models expressed symbolically as opposed to black-box neural networks, the latter only being used during learning to efficiently parameterize the large search space of possible models. In point (ii), learned models are dynamically evaluated a posteriori in the computational solver instead of based on a priori information from preprocessed high-fidelity data, thereby accounting for the specificity of the solver at hand such as its numerics. Finally in point (iii), the exploration of new models is solely guided by predefined integral quantities, e.g., averaged quantities of engineering interest in Reynolds-averaged or large-eddy simulations (LES). We use a coupled deep reinforcement learning framework and computational solver to concurrently achieve these objectives. The combination of reinforcement learning with objectives (i), (ii) and (iii) differentiate our work from previous modeling attempts based on machine learning. In this report, we provide a high-level description of the model discovery framework with reinforcement learning. The method is detailed for the application of discovering missing terms in differential equations. An elementary instantiation of the method is described that discovers missing terms in the Burgers' equation. △ Less

Submitted 29 December, 2019; originally announced January 2020.

Inter-scale energy transfer in turbulence from the viewpoint of subfilter scales

Authors: J. I. Cardesa, A. Lozano-Durán

Abstract: This brief is organized as follows. In Section 2 we introduce the filter-scale energy equations that lead to the subfilter-scale dissipation term P. We discuss the properties of the sources and sinks of kinetic energy, and identify our target features concerning an optimal inter-scale energy transfer (ISET) term. A similar analysis is performed in Section 3 for the kinetic energy associated with t… ▽ More This brief is organized as follows. In Section 2 we introduce the filter-scale energy equations that lead to the subfilter-scale dissipation term P. We discuss the properties of the sources and sinks of kinetic energy, and identify our target features concerning an optimal inter-scale energy transfer (ISET) term. A similar analysis is performed in Section 3 for the kinetic energy associated with the subfilter velocities, and a new ISET term T is introduced. Comparisons between T and P for different flows and for various filters are offered in Section 4 and Section 5. Finally, we present our conclusions in Section 6. △ Less

Submitted 25 December, 2019; v1 submitted 23 December, 2019; originally announced December 2019.

Comments: Submitted for publication as a 2019 Annual Research Brief

Alternative physics to understand wall turbulence: Navier-Stokes equations with modified linear dynamics

Authors: Adrán Lozano-Durán, Marios-Andreas Nikolaidis, Navid C. Constantinou, Michael Karp

Abstract: Despite the nonlinear nature of wall turbulence, there is evidence that the energy-injection mechanisms sustaining wall turbulence can be ascribed to linear processes. The different scenarios stem from linear stability theory and comprise exponential instabilities from mean-flow inflection points, transient growth from non-normal operators, and parametric instabilities from temporal mean-flow vari… ▽ More Despite the nonlinear nature of wall turbulence, there is evidence that the energy-injection mechanisms sustaining wall turbulence can be ascribed to linear processes. The different scenarios stem from linear stability theory and comprise exponential instabilities from mean-flow inflection points, transient growth from non-normal operators, and parametric instabilities from temporal mean-flow variations, among others. These mechanisms, each potentially capable of leading to the observed turbulence structure, are rooted in simplified theories and conceptual arguments. Whether the flow follows any or a combination of them remains unclear. In the present study, we devise a collection of numerical experiments in which the Navier-Stokes equations are sensibly modified to quantify the role of the different linear mechanisms. This is achieved by direct numerical simulation of turbulent channel flows with constrained energy extraction from the streamwise-averaged mean-flow. We demonstrate that (i) transient growth alone is not sufficient to sustain wall turbulence and (ii) the flow remains turbulent when the exponential instabilities are suppressed. On the other hand, we show that (iii) transient growth combined with the parametric instability of the time-varying mean-flow is able to sustain turbulence. △ Less

Submitted 13 December, 2019; originally announced December 2019.

Comments: arXiv admin note: substantial text overlap with arXiv:1909.05490

Effect of a weak current on wind-generated waves in the wrinkle regime

Authors: C. Nove-Josserand, S. Perrard, A. Lozano-Duran, M. Benzaquen, M. Rabaud, F. Moisy

Abstract: We investigate numerically the influence of a weak current on wind-generated surface deformations for wind velocity below the onset of regular waves. In that regime, the liquid surface is populated by small disorganised deformations elongated in the wind direction, referred to as wrinkles. These wrinkles are the superposition of incoherent wakes generated by the pressure fluctuations traveling in… ▽ More We investigate numerically the influence of a weak current on wind-generated surface deformations for wind velocity below the onset of regular waves. In that regime, the liquid surface is populated by small disorganised deformations elongated in the wind direction, referred to as wrinkles. These wrinkles are the superposition of incoherent wakes generated by the pressure fluctuations traveling in the turbulent boundary layer in the air. In this work, we account for the effect of a weak sheared current in the liquid, either longitudinal or transverse, by introducing a modified Doppler-shifted dispersion relation to lowest order in viscosity and current in the spectral theory previously derived by Perrard et al (2019). This theory describes the simplified one-way problem of surface deformations excited by a prescribed turbulent forcing, thereby neglecting the retroaction of waves on turbulence in the air. The forcing is taken from a set of direct numerical simulations of a turbulent channel flow. We determine the wrinkle properties (size and amplitude) as a function of the liquid viscosity and current properties (surface velocity, thickness and orientation). We find significant modifications of the wrinkle geometry by the currents: the wrinkles are tilted for a transverse current, and show finer scales for a longitudinal current. However, their characteristic size is weakly affected, and their amplitude remains independent of the current. We discuss the implications of these results on the onset of regular waves at larger wind velocity. In this work, we introduce a spectral interpolation method to evaluate the surface deformation fields, based on a refined meshing close to the dispersion relation of the waves. This method, which can be extended to any dispersive system excited by a random forcing, strongly reduces the discretization effects at a low computational cost. △ Less

Submitted 28 November, 2020; v1 submitted 15 December, 2019; originally announced December 2019.

Journal ref: Phys. Rev. Fluids 5, 124801 (2020)

Nonlinear mechanism of the self-sustaining process in the buffer and logarithmic layer of wall-bounded flows

Authors: Hyunji Jane Bae, Adrian Lozano-Duran, Beverley J. McKeon

Abstract: The nonlinear mechanism in the self-sustaining process (SSP) of wall-bounded turbulence is investigated. Resolvent analysis is used to identify the principal forcing mode which produces the maximum amplification of the velocities in numerical simulations of the minimal channel for the buffer layer and a modified logarithmic (log) layer. The wavenumbers targeted in this study are those of the funda… ▽ More The nonlinear mechanism in the self-sustaining process (SSP) of wall-bounded turbulence is investigated. Resolvent analysis is used to identify the principal forcing mode which produces the maximum amplification of the velocities in numerical simulations of the minimal channel for the buffer layer and a modified logarithmic (log) layer. The wavenumbers targeted in this study are those of the fundamental mode that is infinitely long in the streamwise direction and once periodic in the spanwise direction. The identified mode is then projected out from the nonlinear term of the Navier-Stokes equations at each time step from the simulation of the corresponding minimal channel. The results show that the removal of the principal forcing mode of the fundamental wavenumber can inhibit turbulence in both the buffer and log layer, with the effect being greater in the buffer layer. Removing other modes instead of the principal mode of the fundamental wavenumber only marginally affects the flow. Closer inspection of the dyadic interactions in the nonlinear term shows that contributions toward the principal forcing mode come from a limited set of wavenumber interactions. Using conditional averaging, the flow structures that are responsible for generating the nonlinear interaction to self-sustain turbulence are identified as spanwise rolls interacting with oblique streaks. This method, based on the equations of motion, validates the similarities in the SSP of the buffer and log layer, and characterises the underlying quadratic interactions in the SSP of the minimal channel. △ Less

Submitted 13 November, 2020; v1 submitted 16 November, 2019; originally announced November 2019.

Journal ref: Journal of Fluid Mechanics, 914, A3 (2021)

Wall turbulence without modal instability of the streaks

Authors: Adrián Lozano-Durán, Marios-Andreas Nikolaidis, Navid C. Constantinou, Michael Karp

Abstract: Despite the nonlinear nature of wall turbulence, there is evidence that the mechanism underlying the energy transfer from the mean flow to the turbulent fluctuations can be ascribed to linear processes. One of the most acclaimed linear instabilities for this energy transfer is the modal growth of perturbations with respect to the streamwise-averaged flow (or streaks). Here, we devise a numerical e… ▽ More Despite the nonlinear nature of wall turbulence, there is evidence that the mechanism underlying the energy transfer from the mean flow to the turbulent fluctuations can be ascribed to linear processes. One of the most acclaimed linear instabilities for this energy transfer is the modal growth of perturbations with respect to the streamwise-averaged flow (or streaks). Here, we devise a numerical experiment in which the Navier--Stokes equations are sensibly modified to suppress these modal instabilities. Our results demonstrate that wall turbulence is sustained with realistic mean and fluctuating velocities despite the absence of streak instabilities. △ Less

Submitted 12 September, 2019; originally announced September 2019.

Coherent structures and secondary flow in turbulent square duct

Authors: Marco Atzori, Ricardo Vinuesa, Adrián Lozano-Durán, Philipp Schlatter

Abstract: The aim of the present work is to investigate the role of coherent structures in the generation of the secondary flow in a turbulent square duct. The coherent structures are defined as connected regions of flow where the product of the instantaneous fluctuations of two velocity components is higher than a threshold based on the long-time turbulence statistics, in the spirit of the three-dimensiona… ▽ More The aim of the present work is to investigate the role of coherent structures in the generation of the secondary flow in a turbulent square duct. The coherent structures are defined as connected regions of flow where the product of the instantaneous fluctuations of two velocity components is higher than a threshold based on the long-time turbulence statistics, in the spirit of the three-dimensional quadrant analysis proposed by Lozano-Durán et al. (J. Fluid Mech., vol. 694, 2012, pp. 100-130). We consider both the direct contribution of the structures to the mean in-plane velocity components and their geometrical properties. The instantaneous phenomena taking place in the turbulent duct are compared with turbulent channel flow at Reynolds numbers of $Re_τ=180$ and $360$, based on friction velocity at the center-plane and channel half height. In the core region of the duct, the fractional contribution of intense events to the wall-normal component of the mean velocity is in very good agreement with that in the channel, despite the presence of the secondary flow in the former. Additionally, the shapes of the three-dimensional objects do not differ significantly in both flows. On the other hand, in the corner region of the duct, the proximity of the walls affects both the geometrical properties of the coherent structures and the contribution to the mean component of the vertical velocity, which is less relevant than that of the complementary portion of the flow not included in such objects. Our results show however that strong Reynolds shear-stress events, despite the differences observed between channel and duct, do not contribute directly to the secondary motion, and thus other phenomena need to be considered instead. △ Less

Submitted 3 June, 2019; originally announced June 2019.

Comments: Submitted to Journal of Fluid Mechanics on June 2nd, 2019

Non-equilibrium three-dimensional boundary layers at moderate Reynolds numbers

Authors: Adrián Lozano-Durán, Marco Giometto, George I. Park, Parviz Moin

Abstract: Non-equilibrium wall turbulence with mean-flow three-dimensionality is ubiquitous in geophysical and engineering flows. Under these conditions, turbulence may experience a counter-intuitive depletion of the turbulent stresses, which has important implications for modelling and control. Yet, current turbulence theories have been established mainly for statistically two-dimensional equilibrium flows… ▽ More Non-equilibrium wall turbulence with mean-flow three-dimensionality is ubiquitous in geophysical and engineering flows. Under these conditions, turbulence may experience a counter-intuitive depletion of the turbulent stresses, which has important implications for modelling and control. Yet, current turbulence theories have been established mainly for statistically two-dimensional equilibrium flows and are unable to predict the reduction in the Reynolds stress magnitude. In the present work, we propose a multiscale model which explains the response of non-equilibrium wall-bounded turbulence under the imposition of three-dimensional strain. The analysis is performed via direct numerical simulation of transient three-dimensional turbulent channels subjected to a sudden lateral pressure gradient at friction Reynolds numbers up to 1,000. We show that the flow regimes and scaling properties of the Reynolds stress are consistent with a model comprising momentum-carrying eddies with sizes and time scales proportional to their distance to the wall. We further demonstrate that the reduction in Reynolds stress follows a spatially and temporally self-similar evolution caused by the relative horizontal displacement between the core of the momentum-carrying eddies and the flow layer underneath. Inspection of the flow energetics reveals that this mechanism is associated with lower levels of pressure-strain correlation which ultimately inhibits the generation of Reynolds stress. Finally, we assess the ability of the state-of-the-art wall-modelled large-eddy simulation to predict non-equilibrium, three-dimensional flows. △ Less

Submitted 29 May, 2019; v1 submitted 24 April, 2019; originally announced April 2019.

Comments: typos corrected

The coherent structure of the kinetic energy transfer in shear turbulence

Authors: Siwei Dong, Yongxiang Huang, Xianxu Yuan, Adrián Lozano-Durán

Abstract: The cascade of energy in turbulent flows, i.e., the transfer of kinetic energy from large to small flow scales or vice versa (backward cascade), is the cornerstone of most theories and models of turbulence since the 1940s. Yet, understanding the spatial organisation of kinetic energy transfer remains an outstanding challenge in fluid mechanics. Here, we unveil the three-dimensional structure of th… ▽ More The cascade of energy in turbulent flows, i.e., the transfer of kinetic energy from large to small flow scales or vice versa (backward cascade), is the cornerstone of most theories and models of turbulence since the 1940s. Yet, understanding the spatial organisation of kinetic energy transfer remains an outstanding challenge in fluid mechanics. Here, we unveil the three-dimensional structure of the energy cascade across the shear-dominated scales using numerical data of homogeneous shear turbulence. We show that the characteristic flow structure associated with the energy transfer is a vortex shaped as an inverted hairpin followed by an upright hairpin. The asymmetry between the forward and backward cascade arises from the opposite flow circulation within the hairpins, which triggers reversed patterns in the flow. △ Less

Submitted 7 March, 2020; v1 submitted 4 April, 2019; originally announced April 2019.

Comments: improved wording, fixed figure labels, some rephrasing

Causality of energy-containing eddies in wall turbulence

Authors: Adrián Lozano-Durán, H. Jane Bae, Miguel P. Encinar

Abstract: Turbulent flows in the presence of walls may be apprehended as a collection of momentum- and energy-containing eddies (energy-eddies), whose sizes differ by many orders of magnitude. These eddies follow a self-sustaining cycle, i.e., existing eddies are seeds for the inception of new ones, and so forth. Understanding this process is critical for the modelling and control of geophysical and industr… ▽ More Turbulent flows in the presence of walls may be apprehended as a collection of momentum- and energy-containing eddies (energy-eddies), whose sizes differ by many orders of magnitude. These eddies follow a self-sustaining cycle, i.e., existing eddies are seeds for the inception of new ones, and so forth. Understanding this process is critical for the modelling and control of geophysical and industrial flows, in which a non-negligible fraction of the energy is dissipated by turbulence in the immediate vicinity of walls. In this study, we examine the causal interactions of energy-eddies in wall-bounded turbulence by quantifying how the knowledge of the past states of eddies reduces the uncertainty of their future states. The analysis is performed via direct numerical simulation (DNS) of turbulent channel flows in which time-resolved energy-eddies are isolated at a prescribed scale. Our approach unveils, in a simple manner, that causality of energy-eddies in the buffer and logarithmic layers is similar and independent of the eddy size. We further show an example of how novel flow control and modelling strategies can take advantage of such self-similar causality. △ Less

Submitted 25 September, 2019; v1 submitted 23 February, 2019; originally announced February 2019.

Comments: even more improved wording and typos corrected, format changed

Journal ref: Journal of Fluid Mechanics, 882, A2 (2019)

DNS-aided explicitly filtered LES

Authors: H. Jane Bae, Adrián Lozano-Durán

Abstract: The equations for LES are formally derived by low-pass filtering the NS equations with the effect of the small scales on the larger ones captured by a SGS model. However, it is known that the LES equations usually employed in practical applications are inconsistent with the filter operator when no explicit filter is used. Moreover, even for explicitly filtered LES, some inconsistencies remain in t… ▽ More The equations for LES are formally derived by low-pass filtering the NS equations with the effect of the small scales on the larger ones captured by a SGS model. However, it is known that the LES equations usually employed in practical applications are inconsistent with the filter operator when no explicit filter is used. Moreover, even for explicitly filtered LES, some inconsistencies remain in the wall-normal direction due to the constraining effect of the wall. A typically undesirable effect stemming from the inconsistency between LES and the filtered NS (fNS) equations is that DNS data can not be used to aid SGS modeling. Using the proposed form of the incompressible fNS equations for flows over flat walls in which the continuity equation, the filter operator, and the boundary condition at the wall are consistently formulated, we compute a set of plane turbulent channel flow explicitly filtered LES at different grid resolutions with the SGS term computed from a parallel DNS computation starting with an equivalent initial condition. The results show that the SGS term computed from DNS provide a truthful SGS model for short times until the numerical errors compound and the two parallel computations diverge. The effective time where the statistics remain similar is longer for lower-order statistics such as mean velocity profile, as expected. This demonstrates that with the consistent fNS equations, the abundant DNS data can be used to inform SGS model development for explicitly filtered LES. Moreover, the diverging of statitics at longer times shows that the SGS contribution is not necessarily homogeous as most SGS models assume, and a more precise correlation between the instantaneous LES velocity field and SGS model is required to capture the dynamics of turbulence. △ Less

Submitted 7 February, 2019; originally announced February 2019.

Journal ref: Center for Turbulence Research, Annual Research Briefs 2018

Wall turbulence with constrained energy extraction from mean flow

Authors: Adrián Lozano-Durán, Michael Karp, Navid. C. Constantinou

Abstract: We study the mechanism of energy injection from the mean flow to the fluctuating velocity necessary to maintain wall turbulence. This process is believed to be correctly represented by the linearized Navier--Stokes equations, and three potential linear mechanisms have been considered, namely, modal instability of the streamwise mean cross-flow $U(y,z,t)$, non-modal transient growth, and non-modal… ▽ More We study the mechanism of energy injection from the mean flow to the fluctuating velocity necessary to maintain wall turbulence. This process is believed to be correctly represented by the linearized Navier--Stokes equations, and three potential linear mechanisms have been considered, namely, modal instability of the streamwise mean cross-flow $U(y,z,t)$, non-modal transient growth, and non-modal transient growth supported by parametric instability. We have designed three numerical experiments of plane turbulent channel flow with additional forcing terms aiming to neutralize one or various linear mechanisms for energy extraction. From our preliminary experiments, only cases with mean cross-flows capable of supporting modal instabilities were found to sustain turbulence. However, the question whether such a new turbulence complies with the same physical mechanisms as those occurring in actual (unforced) turbulence remains unanswered. On the other hand, cases exclusively supported by transient growth decayed until laminarization. △ Less

Submitted 13 February, 2019; v1 submitted 5 February, 2019; originally announced February 2019.

Comments: Corrected typo in author name

Journal ref: Center for Turbulence Research Annual Research Briefs 2018

Resolvent-based estimation of space-time flow statistics

Authors: Aaron Towne, Adrian Lozano-Duran, Xiang Yang

Abstract: We develop a method to estimate space-time flow statistics from a limited set of known data. While previous work has focused on modeling spatial or temporal statistics independently, space-time statistics carry fundamental information about the physics and coherent motions of the flow and provide a starting point for low-order modeling and flow control efforts. The method is derived using a statis… ▽ More We develop a method to estimate space-time flow statistics from a limited set of known data. While previous work has focused on modeling spatial or temporal statistics independently, space-time statistics carry fundamental information about the physics and coherent motions of the flow and provide a starting point for low-order modeling and flow control efforts. The method is derived using a statistical interpretation of resolvent analysis. The central idea of our approach is to use known data to infer the statistics of the nonlinear terms that constitute a forcing on the linearized Navier-Stokes equations, which in turn imply values for the remaining unknown flow statistics through application of the resolvent operator. Rather than making an a priori rank-1 assumption, our method allows the known input data to select the most relevant portions of the resolvent operator for describing the data, making it well-suited for high-rank turbulent flows. We demonstrate the predictive capabilities of the method using two examples: the Ginzburg-Landau equation, which serves as a convenient model for a convectively unstable flow, and a turbulent channel flow at low Reynolds number. △ Less

Submitted 22 January, 2019; originally announced January 2019.

Characteristic scales of Townsend's wall attached eddies

Authors: Adrián Lozano-Durán, Hyunji Jane Bae

Abstract: Townsend (1976) proposed a structural model for the logarithmic layer (log-layer) of wall turbulence at high Reynolds numbers, where the dominant momentum-carrying motions are organised into a multi-scale population of eddies attached to the wall. In the attached eddy framework, the relevant length and velocity scales of the wall-attached eddies are the friction velocity and the distance to the wa… ▽ More Townsend (1976) proposed a structural model for the logarithmic layer (log-layer) of wall turbulence at high Reynolds numbers, where the dominant momentum-carrying motions are organised into a multi-scale population of eddies attached to the wall. In the attached eddy framework, the relevant length and velocity scales of the wall-attached eddies are the friction velocity and the distance to the wall. In the present work, we hypothesise that the momentum-carrying eddies are controlled by the mean momentum flux and mean shear with no explicit reference to the distance to the wall and propose new characteristic velocity, length, and time scales consistent with this argument. Our hypothesis is supported by direct numerical simulation (DNS) of turbulent channel flows driven by non-uniform body forces and modified mean velocity profiles, where the resulting outer-layer flow structures are substantially altered to accommodate the new mean momentum transfer. The proposed scaling is further corroborated by simulations where the no-slip wall is replaced by a Robin boundary condition for the three velocity components, allowing for substantial wall-normal transpiration at all lengths scales. We show that the outer-layer one-point statistics and spectra of this channel with transpiration agree quantitatively with those of its wall-bounded counterpart. The results reveal that the wall-parallel no-slip condition is not required to recover classic wall-bounded turbulence far from the wall and, more importantly, neither is the impermeability condition at the wall. △ Less

Submitted 5 March, 2019; v1 submitted 14 January, 2019; originally announced January 2019.

Comments: Section 3 revised, typos corrected

Journal ref: Journal of Fluid Mechanics, 868, 698-725 (2019)

Dynamic slip wall model for large-eddy simulation

Authors: H. Jane Bae, Adrián Lozano-Durán, Sanjeeb T. Bose, Parviz Moin

Abstract: Wall modelling in large-eddy simulation (LES) is necessary to overcome the prohibitive near-wall resolution requirements in high-Reynolds-number turbulent flows. Most existing wall models rely on assumptions about the state of the boundary layer and require a priori prescription of tunable coefficients. They also impose the predicted wall stress by replacing the no-slip boundary condition at the w… ▽ More Wall modelling in large-eddy simulation (LES) is necessary to overcome the prohibitive near-wall resolution requirements in high-Reynolds-number turbulent flows. Most existing wall models rely on assumptions about the state of the boundary layer and require a priori prescription of tunable coefficients. They also impose the predicted wall stress by replacing the no-slip boundary condition at the wall with a Neumann boundary condition in the wall-parallel directions while maintaining the no-transpiration condition in the wall-normal direction. In the present study, we first motivate and analyse the Robin (slip) boundary condition with transpiration (nonzero wall-normal velocity) in the context of wall-modelled LES. The effect of the slip boundary condition on the one-point statistics of the flow is investigated in LES of turbulent channel flow and flat-plate turbulent boundary layer. It is shown that the slip condition provides a framework to compensate for the deficit or excess of mean momentum at the wall. Moreover, the resulting nonzero stress at the wall alleviates the well-known problem of the wall-stress under-estimation by current subgrid-scale (SGS) models. Secondly, we discuss the requirements for the slip condition to be used in conjunction with wall models and derive the equation that connects the slip boundary condition with the stress at the wall. Finally, a dynamic procedure for the slip coefficients is formulated, providing a dynamic slip wall model free of a priori specified coefficients. The performance of the proposed dynamic wall model is tested in a series of LES of turbulent channel flow, non-equilibrium three-dimensional channel flow, and flat-plate turbulent boundary layer. The results show that the dynamic wall model is able to accurately predict one-point turbulence statistics for various flow configurations, Reynolds numbers, and grid resolutions. △ Less

Submitted 9 October, 2018; originally announced October 2018.

Journal ref: Journal of Fluid Mechanics, 859, 400-432 (2019)

Turbulent windprint on a liquid surface

Authors: Stéphane Perrard, Adrián Lozano-Durán, Marc Rabaud, Michael Benzaquen, Frédéric Moisy

Abstract: We investigate the effect of a light turbulent wind on a liquid surface, below the onset of wave generation. In that regime, the liquid surface is populated by small disorganised deformations elongated in the streamwise direction. Formally identified recently by Paquier et al. (2015), the deformations that occur below the wave onset were named wrinkles. We provide here a theoretical framework for… ▽ More We investigate the effect of a light turbulent wind on a liquid surface, below the onset of wave generation. In that regime, the liquid surface is populated by small disorganised deformations elongated in the streamwise direction. Formally identified recently by Paquier et al. (2015), the deformations that occur below the wave onset were named wrinkles. We provide here a theoretical framework for this wrinkle regime, using the viscous response of a free surface liquid submitted to arbitrary normal and tangential interfacial stresses at its upper boundary. We relate the spatio-temporal spectrum of the surface deformations to that of the applied interfacial pressure and shear stress fluctuations. For that, we evaluate the spatio-temporal statistics of the turbulent forcing using Direct Numerical Simulation of a turbulent air channel flow, assuming no coupling between the air and the liquid flows. Combining theory and numerical simulation, we thus obtain synthetic wrinkles that reproduce previous experimental investigations. We show that the wrinkles are a multi-scale superposition of random wakes generated by the turbulent fluctuations. They result mainly from the nearly isotropic pressure fluctuations generated in the boundary layer, rather than from the elongated shear stress fluctuations. The wrinkle regime described in this paper naturally arises as the viscous-saturated asymptotic of the inviscid growth theory of Phillips (1957). Experiments indicate that the onset of exponential wave growth depends on the liquid viscosity. Our theory suggests that the empirical criterion for the onset is satisfied when the wrinkle amplitude reaches a given fraction of the viscous sublayer thickness. It indicates that the turbulent fluctuations near the onset may play a role in the triggering of exponential wave growth. △ Less

Submitted 12 April, 2019; v1 submitted 19 September, 2018; originally announced September 2018.

Comments: J. Fluid Mech. 2019 (in press)

Journal ref: J. Fluid Mech. 873, 1020-1054 (2019)

Error scaling of large-eddy simulation in the outer region of wall-bounded turbulence

Authors: Adrián Lozano-Durán, Hyunji Jane Bae

Abstract: We study the error scaling properties of large-eddy simulation (LES) in the outer region of wall-bounded turbulence at moderately high Reynolds numbers. In order to avoid the additional complexity of wall-modeling, we perform LES of turbulent channel flows in which the no-slip condition at the wall is replaced by a Neumann condition supplying the exact mean wall-stress. The statistics investigated… ▽ More We study the error scaling properties of large-eddy simulation (LES) in the outer region of wall-bounded turbulence at moderately high Reynolds numbers. In order to avoid the additional complexity of wall-modeling, we perform LES of turbulent channel flows in which the no-slip condition at the wall is replaced by a Neumann condition supplying the exact mean wall-stress. The statistics investigated are the mean velocity profile, turbulence intensities, and kinetic energy spectra. The errors follow $(Δ/L)^αRe_τ^γ$, where $Δ$ is the characteristic grid resolution, $Re_τ$ is the friction Reynolds number, and $L$ is the meaningful length-scale to normalize $Δ$ in order to collapse the errors across the wall-normal distance. We show that $Δ$ can be expressed as the $L_2$-norm of the grid vector and that $L$ is well represented by the ratio of the friction velocity and mean shear. The exponent $α$ is estimated from theoretical arguments for each statistical quantity of interest and shown to roughly match the values computed by numerical simulations. For the mean profile and kinetic energy spectra, $α\approx1$, whereas the turbulence intensities converge at a slower rate $α<1$. The exponent $γ$ is approximately $0$, i.e. the LES solution is independent of the Reynolds number. The expected behavior of the turbulence intensities at high Reynolds numbers is also derived and shown to agree with the classic log-layer profiles for grid resolutions lying within the inertial range. Further examination of the LES turbulence intensities and spectra reveals that both quantities resemble their filtered counterparts from direct numerical simulation (DNS) data, but that the mechanism responsible for this similarity is related to the balance between the input power and dissipation rather than to filtering. △ Less

Submitted 4 April, 2019; v1 submitted 5 June, 2018; originally announced June 2018.

Comments: Improved clarity and typo correction

Journal ref: Journal of Computational Physics 392, 532-555 (2019)

Multiscale analysis of the topological invariants in the logarithmic region of turbulent channels at $Re_τ=932$

Authors: Adrián Lozano-Durán, Markus Holzner, Javier Jiménez

Abstract: The invariants of the velocity gradient tensor, $R$ and $Q$, and their enstrophy and strain components are studied in the logarithmic layer of an incompressible turbulent channel flow. The velocities are filtered in the three spatial directions and the results analyzed at different scales. We show that the $R$--$Q$ plane does not capture the changes undergone by the flow as the filter width increa… ▽ More The invariants of the velocity gradient tensor, $R$ and $Q$, and their enstrophy and strain components are studied in the logarithmic layer of an incompressible turbulent channel flow. The velocities are filtered in the three spatial directions and the results analyzed at different scales. We show that the $R$--$Q$ plane does not capture the changes undergone by the flow as the filter width increases, and that the enstrophy/enstrophy-production and strain/strain-production planes represent better choices. We also show that the conditional mean trajectories may differ significantly from the instantaneous behavior of the flow since they are the result of an averaging process where the mean is 3-5 times smaller than the corresponding standard deviation. The orbital periods in the $R$--$Q$ plane are shown to be independent of the intensity of the events, and of the same order of magnitude than those in the enstrophy/enstrophy-production and strain/strain-production planes. Our final goal is to test whether the dynamics of the flow are self-similar in the inertial range, and the answer turns out to be no. The mean shear is found to be responsible for the absence of self-similarity and progressively controls the dynamics of the eddies observed as the filter width increases. However, a self-similar behavior emerges when the calculations are repeated for the fluctuating velocity gradient tensor. Finally, the turbulent cascade in terms of vortex stretching is considered by computing the alignment of the vorticity at a given scale with the strain at a different one. These results generally support a non-negligible role of the phenomenological energy-cascade model formulated in terms of vortex stretching. △ Less

Submitted 4 September, 2017; originally announced September 2017.

Journal ref: J. Fluid Mech. (2016), vol. 803, pp. 356--394

An efficient algorithm to compute the genus of discrete surfaces and applications to turbulent flows

Authors: Adrián Lozano-Durán, Guillem Borrell

Abstract: A simple and efficient algorithm to numerically compute the genus of surfaces of three-dimensional objects using the Euler characteristic formula is presented. The algorithm applies to objects obtained by thresholding a scalar field in a structured-collocated grid, and does not require any triangulation of the data. This makes the algorithm fast, memory-efficient and suitable for large datasets. A… ▽ More A simple and efficient algorithm to numerically compute the genus of surfaces of three-dimensional objects using the Euler characteristic formula is presented. The algorithm applies to objects obtained by thresholding a scalar field in a structured-collocated grid, and does not require any triangulation of the data. This makes the algorithm fast, memory-efficient and suitable for large datasets. Applications to the characterization of complex surfaces in turbulent flows are presented to illustrate the method. △ Less

Submitted 4 September, 2017; originally announced September 2017.

ACM Class: G.4; I.3; J.2

Journal ref: ACM Trans. Math. Softw. 42, 4, Article 34 (June 2016)

A statistical state dynamics-based study of the structure and mechanism of large-scale motions in plane Poiseuille flow

Authors: Brian F. Farrell, Petros J. Ioannou, Javier Jiménez, Navid C. Constantinou, Adrián Lozano-Durán, Marios-Andreas Nikolaidis

Abstract: The perspective of statistical state dynamics (SSD) has recently been applied to the study of mechanisms underlying turbulence in various physical systems. An example implementation of SSD is the second order closure referred to as stochastic structural stability theory (S3T), which has provided insight into the dynamics of wall turbulence and specifically the emergence and maintenance of the roll… ▽ More The perspective of statistical state dynamics (SSD) has recently been applied to the study of mechanisms underlying turbulence in various physical systems. An example implementation of SSD is the second order closure referred to as stochastic structural stability theory (S3T), which has provided insight into the dynamics of wall turbulence and specifically the emergence and maintenance of the roll/streak structure. This closure eliminates nonlinear interactions among the perturbations has been removed, restricting nonlinearity in the dynamics to that of the mean equation and the interaction between the mean and perturbation covariance. Here, this quasi-linear restriction of the dynamics is used to study the structure and dynamics of turbulence in plane Poiseuille flow at moderately high Reynolds numbers in a closely related dynamical system, referred to as the restricted nonlinear (RNL) system. RNL simulations reveal that the essential features of wall-turbulence dynamics are retained. Remarkably, the RNL system spontaneously limits the support of its turbulence to a small set of streamwise Fourier components giving rise to a naturally minimal representation of its turbulence dynamics. Although greatly simplified, this RNL turbulence exhibits natural-looking structures and statistics. Surprisingly, even when further truncation of the perturbation support to a single streamwise component is imposed the RNL system continues to produce self-sustaining turbulent structure and dynamics. RNL turbulence at the Reynolds numbers studied is dominated by the roll/streak structure in the buffer layer and similar very-large-scale structure (VLSM) in the outer layer. Diagnostics of the structure, spectrum and energetics of RNL and DNS turbulence are used to demonstrate that the roll/streak dynamics supporting the turbulence in the buffer and logarithmic layer is essentially similar in RNL and DNS. △ Less

Submitted 14 November, 2016; v1 submitted 18 December, 2015; originally announced December 2015.

Comments: published in J. Fluid Mech., 26 pages, 15 figures. arXiv admin note: text overlap with arXiv:1401.7816

Journal ref: J. Fluid Mech., 2016, 809, pp. 290-315

Turbulence in the highly restricted dynamics of a closure at second order: comparison with DNS

Authors: Navid C. Constantinou, Adrian Lozano-Durán, Marios-Andreas Nikolaidis, Brian F. Farrell, Petros J. Ioannou, Javier Jiménez

Abstract: S3T (Stochastic Structural Stability Theory) employs a closure at second order to obtain the dynamics of the statistical mean turbulent state. When S3T is implemented as a coupled set of equations for the streamwise mean and perturbation states, nonlinearity in the dynamics is restricted to interaction between the mean and perturbations. The S3T statistical mean state dynamics can be approximately… ▽ More S3T (Stochastic Structural Stability Theory) employs a closure at second order to obtain the dynamics of the statistical mean turbulent state. When S3T is implemented as a coupled set of equations for the streamwise mean and perturbation states, nonlinearity in the dynamics is restricted to interaction between the mean and perturbations. The S3T statistical mean state dynamics can be approximately implemented by similarly restricting the dynamics used in a direct numerical simulation (DNS) of the full Navier-Stokes equations (referred to as the NS system). Although this restricted nonlinear system (referred to as the RNL system) is greatly simplified in its dynamics in comparison to the associated NS, it nevertheless self-sustains a turbulent state in wall-bounded shear flow with structures and dynamics comparable to that in observed turbulence. Moreover, RNL turbulence can be analyzed effectively using theoretical methods developed to study the closely related S3T system. In order to better understand RNL turbulence and its relation to NS turbulence, an extensive comparison is made of diagnostics of structure and dynamics in these systems. Although quantitative differences are found, the results show that turbulence in the RNL system closely parallels that in NS and suggest that the S3T/RNL system provides a promising reduced complexity model for studying turbulence in wall-bounded shear flows. △ Less

Submitted 30 January, 2014; originally announced January 2014.

Comments: 18 pages, 14 figures

Journal ref: J. Phys.: Conf. Ser., 2014, 506, 012004

The 3D structure of turbulent channels up to $\mathrm{Re}_τ=4000$

Authors: Adrián Lozano-Durán, Javier Jiménez

Abstract: Several time-resolved fluid dynamics videos of turbulent channels from $\mathrm{Re}_τ=180$ to $\mathrm{Re}_τ=4000$ are presented. The videos show the temporal evolution of sweeps (bluish) and ejections (reddish) in one half of the channel (only the bottom wall is shown). The color changes from dark for points close to the wall, to bright for those reaching the center of the channel. As the Reynold… ▽ More Several time-resolved fluid dynamics videos of turbulent channels from $\mathrm{Re}_τ=180$ to $\mathrm{Re}_τ=4000$ are presented. The videos show the temporal evolution of sweeps (bluish) and ejections (reddish) in one half of the channel (only the bottom wall is shown). The color changes from dark for points close to the wall, to bright for those reaching the center of the channel. As the Reynolds number increases the scale separation becomes more clear and the complexity of the dynamics observed rises. △ Less

Submitted 12 October, 2012; originally announced October 2012.

Comments: Videos are included