BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Parviz Moin

Parviz Moin contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.flu-dynphysics.comp-phphysics.ao-phComputational Engineering, Finance, and Sciencephysics.geo-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
14works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

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

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

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

Building this graph slice

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

Published work

14 published item(s)

preprint2023arXiv

Reynolds number dependence of length scales governing turbulent flow separation with application to wall-modeled large-eddy simulations

This article proposes a Reynolds number scaling of the required grid points to perform wall-modeled LES of turbulent flows encountering separation off a solid surface. Based on comparisons between the various time scales in a non-equilibrium (due to the action of an external pressure gradient) turbulent boundary layer, a simple definition of the near-wall ``under-equilibrium" and ``out-of-equilibrium" scales is put forward (where ``under-equilibrium" refers to scales governed by a quasi-balance between the viscous and the pressure gradient terms). It is shown that the former length scale varies with Reynolds number as lp Re^(-2/3). The same scaling is obtained from a simplified Green's function solution of the Poisson equation in the vicinity of the separation point. A-priori analysis demonstrates that the resolution required to reasonably predict the wall-shear stress (for example, errors lower than approximately 10-15% in the entire domain) in several nonequilibrium flows is at least O(10) lp irrespective of the Reynolds number and the Clauser parameter. Further, a series of a-posteriori validation studies are performed to determine the accuracy of this scaling including the flow over the Boeing speed bump, Song-Eaton diffuser, Notre-Dame Ramp, and the backward-facing step. The results suggest that for these flows, scaling the computational grids () such that / lp is independent of the Reynolds number results in accurate predictions of flow separation at the same ``nominal" grid resolution across different Reynolds numbers. Finally, it is suggested that in the vicinity of the separation and reattachment points, the grid-point requirements for wall-modeled large eddy simulations may scale as Re^4/3, which is more restrictive than the previously proposed flat-plate boundary layer-based estimates (Re1) of Choi and Moin (Phys. Fluids, 2012) and Yang and Griffin (Phys. Fluids, 2021).

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A kinetic energy--and entropy-preserving scheme for compressible two-phase flows

Accurate numerical modeling of compressible flows, particularly in the turbulent regime, requires a method that is non-dissipative and stable at high Reynolds ($Re$) numbers. For a compressible flow, it is known that discrete conservation of kinetic energy is not a sufficient condition for numerical stability, unlike in incompressible flows. In this study, we adopt the recently developed conservative diffuse-interface method (Jain, Mani $\&$ Moin, $\textit{J. Comput. Phys.}$, 2020) along with the five-equation model for the simulation of compressible two-phase flows. This method discretely conserves the mass of each phase, momentum, and total energy of the system. We here propose discrete consistency conditions between the numerical fluxes, such that any set of numerical fluxes that satisfy these conditions would not spuriously contribute to the kinetic energy and entropy of the system. We also present a set of numerical fluxes\textemdash which satisfies these consistency conditions\textemdash that results in an exact conservation of kinetic energy and approximate conservation of entropy in the absence of pressure work, viscosity, thermal diffusion effects, and time-discretization errors. Since the model consistently reduces to the single-phase Navier-Stokes system when the properties of the two phases are identical, the proposed consistency conditions and numerical fluxes are also applicable for a broader class of single-phase flows. To this end, we present coarse-grid numerical simulations of compressible single-phase and two-phase turbulent flows at infinite $Re$, to illustrate the stability of the proposed method in canonical test cases, such as an isotropic turbulence and Taylor-Green vortex flows. A higher-resolution simulation of a droplet-laden compressible decaying isotropic turbulence is also presented, and the effect of the presence of droplets on the flow is analyzed.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Prediction of aerothermal characteristics of a generic hypersonic inlet flow

The accurate prediction of aerothermal surface loading is of paramount importance for the design of high speed flight vehicles. In this work, we consider the numerical solution of hypersonic flow over a double-finned geometry, representative of the inlet of an air-breathing flight vehicle, characterized by three-dimensional intersecting shock-wave/turbulent boundary-layer interaction at Mach 8.3. High Reynolds numbers ($Re_L \approx 11.6 \times 10^6$ based on free-stream conditions) and the presence of cold walls ($T_w/T_o \approx 0.27$) leading to large near-wall temperature gradients necessitate the use of wall-modeled large-eddy simulation (WMLES) in order to make calculations computationally tractable. The comparison of the WMLES results with experimental measurements shows good agreement in the time-averaged surface heat flux and wall pressure distributions, and the WMLES predictions show reduced errors with respect to the experimental measurements than prior RANS calculations. The favorable comparisons are obtained using an LES wall model based on equilibrium boundary layer approximations despite the presence of numerous non-equilibrium conditions including three dimensionality, shock-boundary layer interactions, and flow separation. Lastly, it is also demonstrated that the use of semi-local eddy viscosity scaling (in lieu of the commonly used van Driest scaling) in the LES wall model is necessary to accurately predict the surface pressure loading and heat fluxes.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

General method for determining the boundary layer thickness in nonequilibrium flows

While the computation of the boundary-layer thickness is straightforward for canonical equilibrium flows, there are no established definitions for general non-equilibrium flows. In this work, a method is developed based on a local reconstruction of the "inviscid" velocity profile $U_I[y]$ resulting from the application of the Bernoulli equation in the wall-normal direction. The boundary-layer thickness $δ_{99}$ is then defined as the location where $U/U_I = 0.99$, which is consistent with its classical definition for the zero-pressure-gradient boundary layers (ZPGBLs). The proposed local-reconstruction method is parameter free and can be deployed for both internal and external flows without resorting to an iterative procedure, numerical integration, or numerical differentiation. The superior performance of the local-reconstruction method over various existing methods is demonstrated by applying the methods to laminar and turbulent boundary layers and two flows over airfoils. Numerical experiments reveal that the local-reconstruction method is more accurate and more robust than existing methods, and it is applicable for flows over a wide range of Reynolds numbers.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A conservative diffuse-interface method for compressible two-phase flows

In this article, we propose a novel conservative diffuse-interface method for the simulation of immiscible compressible two-phase flows. The proposed method discretely conserves the mass of each phase, momentum and total energy of the system. We use the baseline five-equation model and propose interface-regularization (diffusion--sharpening) terms in such a way that the resulting model maintains the conservative property of the underlying baseline model; and lets us use a central-difference scheme for the discretization of all the operators in the model, which leads to a non-dissipative implementation that is crucial for the simulation of turbulent flows and acoustics. Furthermore, the provable strengths of the proposed model are: (a) the model maintains the boundedness property of the volume fraction field, which is a physical realizability requirement for the simulation of two-phase flows, (b) the proposed model is such that the transport of volume fraction field inherently satisfies the total-variation-diminishing property without having to add any flux limiters that destroy the non-dissipative nature of the scheme, (c) the proposed interface-regularization terms in the model do not spuriously contribute to the kinetic energy of the system and therefore do not affect the non-linear stability of the numerical simulation, and (d) the model is consistent with the second law of thermodynamics. Finally, we present numerical simulations using the model and assess (a) the accuracy of evolution of the interface shape, (b) implementation of surface tension effects, (c) propagation of acoustics and their interaction with material interfaces, (d) the accuracy and robustness of the numerical scheme for simulation of complex high-Reynolds-number flows, and (e) performance and scalability of the method.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Identifying and tracking bubbles and drops in simulations: a toolbox for obtaining sizes, lineages, and breakup and coalescence statistics

Knowledge of bubble and drop size distributions in two-phase flows is important for characterizing a wide range of phenomena, including combustor ignition, sonar communication, and cloud formation. The physical mechanisms driving the background flow also drive the time evolution of these distributions. Accurate and robust identification and tracking algorithms for the dispersed phase are necessary to reliably measure this evolution and thereby quantify the underlying mechanisms in interface-resolving flow simulations. The identification of individual bubbles and drops traditionally relies on an algorithm used to identify connected regions. This traditional algorithm can be sensitive to the presence of spurious structures. A cost-effective refinement is proposed to maximize volume accuracy while minimizing the identification of spurious bubbles and drops. An accurate identification scheme is crucial for distinguishing bubble and drop pairs with large size ratios. The identified bubbles and drops need to be tracked in time to obtain breakup and coalescence statistics that characterize the evolution of the size distribution, including breakup and coalescence frequencies, and the probability distributions of parent and child bubble and drop sizes. An algorithm based on mass conservation is proposed to construct bubble and drop lineages using simulation snapshots that are not necessarily from consecutive time-steps. These lineages are then used to detect breakup and coalescence events, and obtain the desired statistics. Accurate identification of large-size-ratio bubble and drop pairs enables accurate detection of breakup and coalescence events over a large size range. Together, these algorithms enable insights into the mechanisms behind bubble and drop formation and evolution in flows of practical importance.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Shock-induced heating and transition to turbulence in a hypersonic boundary layer

The interaction between an incident shock wave and a Mach-6 undisturbed hypersonic laminar boundary layer over a cold wall is addressed using direct numerical simulations (DNS) and wall-modeled large-eddy simulations (WMLES) at different angles of incidence. At sufficiently high shock-incidence angles, the boundary layer transitions to turbulence via breakdown of near-wall streaks shortly downstream of the shock impingement, without the need of any inflow free-stream disturbances. The transition causes a localized significant increase in the Stanton number and skin-friction coefficient, with high incidence angles augmenting the peak thermomechanical loads in an approximately linear way. Statistical analyses of the boundary layer downstream of the interaction for each case are provided that quantify streamwise spatial variations of the Reynolds analogy factors and indicate a breakdown of the Morkovin's hypothesis near the wall, where velocity and temperature become correlated. A modified strong Reynolds analogy with a fixed turbulent Prandtl number is observed to perform best. Conventional transformations fail at collapsing the mean velocity profiles on the incompressible log law. The WMLES prompts transition and peak heating, delays separation, and advances reattachment, thereby shortening the separation bubble. When the shock leads to transition, WMLES provides predictions of DNS peak thermomechanical loads within $\pm 10\%$ at a computational cost lower than DNS by two orders of magnitude. Downstream of the interaction, in the turbulent boundary layer, WMLES agrees well with DNS results for the Reynolds analogy factor, the mean profiles of velocity and temperature, including the temperature peak, and the temperature/velocity correlation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The turbulent bubble break-up cascade. Part 1. Theoretical developments

Breaking waves entrain gas beneath the surface. The wave-breaking process energizes turbulent fluctuations that break bubbles in quick succession to generate a wide range of bubble sizes. Understanding this generation mechanism paves the way towards the development of predictive models for large-scale maritime and climate simulations. Garrett et al. (2000) suggested that super-Hinze-scale turbulent breakup transfers entrained gas from large to small bubble sizes in the manner of a cascade. We provide a theoretical basis for this bubble-mass cascade by appealing to how energy is transferred from large to small scales in the energy cascade central to single-phase turbulence theories. A bubble break-up cascade requires that break-up events predominantly transfer bubble mass from a certain bubble size to a slightly smaller size on average. This property is called locality. In this paper, we analytically quantify locality by extending the population balance equation in conservative form to derive the bubble-mass transfer rate from large to small sizes. Using our proposed measures of locality, we show that scalings relevant to turbulent bubbly flows, including those postulated by Garrett et al. (2000) and observed in breaking-wave experiments and simulations, are consistent with a strongly local transfer rate, where the influence of non-local contributions decays in a power-law fashion. These theoretical predictions are confirmed using numerical simulations in Part 2, revealing key physical aspects of the bubble break-up cascade phenomenology. Locality supports the universality of turbulent small-bubble break-up, which simplifies the development of subgrid-scale models to predict oceanic small-bubble statistics of practical importance.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access