This paper is on an EulerianEulerian (EE) approach that utilizes Godunov’s scheme to deal with a running shock that interacts with a cloud of particles. The EE approach treats both carrier phase (fluid phase) and dispersed phase (particle phase) in the Eulerian frame. In this work, the fluid equations are the Euler equations for the compressible gas while the particle equations are based on a recently developed model to solve for the number density, velocity, temperature, particle subgrid scale stresses, and particle subgrid scale heat fluxes. The carrier and dispersed phases exchange momentum and heat, which are modeled through incorporating source terms in their equations. Carrier and dispersed phase equation form a hyperbolic set of differential equations, which are numerically solved with Godunov’s scheme. The numerical solutions are obtained in this work for a twodimensional normal running shock interacting with a rectangular cloud of particles. The results generated by the EE approach were compared against the results that were generated by a wellstablished EulerianLagragian (EL) approach that treats the carrier phase in an Eulerian frame, while does the dispersed phase in a Lagrangian framework where individuals particles are traced and solved. For the considered configuration, the EE approach reproduced themore »
Sparse identification of multiphase turbulence closures for coupled fluid–particle flows
In this work, model closures of the multiphase Reynoldsaveraged Navier–Stokes (RANS) equations are developed for homogeneous, fully developed gas–particle flows. To date, the majority of RANS closures are based on extensions of singlephase turbulence models, which fail to capture complex twophase flow dynamics across dilute and dense regimes, especially when twoway coupling between the phases is important. In the present study, particles settle under gravity in an unbounded viscous fluid. At sufficient mass loadings, interphase momentum exchange between the phases results in the spontaneous generation of particle clusters that sustain velocity fluctuations in the fluid. Data generated from Eulerian–Lagrangian simulations are used in a sparse regression method for model closure that ensures form invariance. Particular attention is paid to modelling the unclosed terms unique to the multiphase RANS equations (drag production, drag exchange, pressure strain and viscous dissipation). A minimal set of tensors is presented that serve as the basis for modelling. It is found that sparse regression identifies compact, algebraic models that are accurate across flow conditions and robust to sparse training data.
 Award ID(s):
 1846054
 Publication Date:
 NSFPAR ID:
 10285288
 Journal Name:
 Journal of Fluid Mechanics
 Volume:
 914
 ISSN:
 00221120
 Sponsoring Org:
 National Science Foundation
More Like this


Mechanical properties of cellular structures, including the cell cytoskeleton, are increasingly used as biomarkers for disease diagnosis and fundamental studies in cell biology. Recent experiments suggest that the cell cytoskeleton and its permeating cytosol, can be described as a poroelastic (PE) material. Biot theory is the standard model used to describe PE materials. Yet, this theory does not account for the fluid viscous stress, which can lead to inaccurate predictions of the mechanics in the dilute filamentous network of the cytoskeleton. Here, we adopt a twophase model that extends Biot theory by including the fluid viscous stresses in the fluid's momentum equation. We use generalized linear viscoelastic (VE) constitutive equations to describe the permeating fluid and the network stresses and assume a constant friction coefficient that couples the fluid and network displacement fields. As the first step in developing a computational framework for solving the resulting equations, we derive closedform general solutions of the fluid and network displacement fields in spherical coordinates. To demonstrate the applicability of our results, we study the motion of a rigid sphere moving under a constant force inside a PE medium, composed of a linear elastic network and a Newtonian fluid. We find that themore »

This study presents the first 3D twoway coupled fluid structure interaction (FSI) simulation of a hybrid anechoic wind tunnel (HAWT) test section with modeling all important effects, such as turbulence, Kevlar wall porosity and deflection, and reveals for the first time the complete 3D flow structure associated with a lifting model placed into a HAWT. The Kevlar deflections are captured using finite element analysis (FEA) with shell elements operated under a membrane condition. Threedimensional RANS CFD simulations are used to resolve the flow field. Aerodynamic experimental results are available and are compared against the FSI results. Quantitatively, the pressure coefficients on the airfoil are in good agreement with experimental results. The lift coefficient was slightly underpredicted while the drag was overpredicted by the CFD simulations. The flow structure downstream of the airfoil showed good agreement with the experiments, particularly over the wind tunnel walls where the Kevlar windows interact with the flow field. A discrepancy between previous experimental observations and juncture flowinduced vortices at the ends of the airfoil is found to stem from the limited ability of turbulence models. The qualitative behavior of the flow, including airfoil pressures and crosssectional flow structure is well captured in the CFD. Frommore »

The drag force on a spherical intruder in dense granular shear flows is studied using discrete element method simulations. Three regimes of the intruder dynamics are observed depending on the magnitude of the drag force (or the corresponding intruder velocity) and the flow inertial number: a fluctuationdominated regime for small drag forces; a viscous regime for intermediate drag forces; and an inertial (cavity formation) regime for large drag forces. The transition from the viscous regime (linear forcevelocity relation) to the inertial regime (quadratic forcevelocity relation) depends further on the inertial number. Despite these distinct intruder dynamics, we find a quantitative similarity between the intruder drag in granular shear flows and the Stokesian drag on a sphere in a viscous fluid for intruder Reynolds numbers spanning five orders of magnitude. Beyond this firstorder description, a modified Stokes drag model is developed that accounts for the secondary dependence of the drag coefficient on the inertial number and the intruder size and density ratios. When the drag model is coupled with a segregation force model for intruders in dense granular flows, it is possible to predict the velocity of gravitydriven segregation of an intruder particle in shear flow simulations.

The relative velocities and positions of monodisperse highinertia particle pairs in isotropic turbulence are studied using direct numerical simulations (DNS), as well as Langevin simulations (LS) based on a probability density function (PDF) kinetic model for pair relative motion. In a prior study (Rani et al. , J. Fluid Mech. , vol. 756, 2014, pp. 870–902), the authors developed a stochastic theory that involved deriving closures in the limit of high Stokes number for the diffusivity tensor in the PDF equation for monodisperse particle pairs. The diffusivity contained the time integral of the Eulerian twotime correlation of fluid relative velocities seen by pairs that are nearly stationary. The twotime correlation was analytically resolved through the approximation that the temporal change in the fluid relative velocities seen by a pair occurs principally due to the advection of smaller eddies past the pair by largescale eddies. Accordingly, two diffusivity expressions were obtained based on whether the pair centre of mass remained fixed during flow time scales, or moved in response to integralscale eddies. In the current study, a quantitative analysis of the (Rani et al. 2014) stochastic theory is performed through a comparison of the pair statistics obtained using LS with those from DNS.more »