Search

Pore-resolved direct numerical simulations of turbulent open channel flow are performed comparing the structure and dynamics of turbulence over impermeable rough and smooth walls to a porous sediment bed at permeability Reynolds number ($Re_K$) of 2.6, representative of aquatic beds. Four configurations are investigated; namely, (i) permeable bed with randomly packed sediment grains, (ii) an impermeable-wall with full layer of roughness elements matching the top layer of the sediment bed, (iii) an impermeable-wall with half layer of roughness elements , and (iv) a smooth wall. A double-averaging methodology is used to compute the mean velocity, Reynolds stresses, form-induced stresses, and turbulent kinetic energy budget. It is observed that the mean velocity, Reynolds stresses, and form-induced pressure-velocity correlations representing upwelling and downwelling fluxes are similar in magnitude for the permeable-bed and impermeable-full layer wall cases. However, for the impermeable-half layer case, the wall blocking effect results in higher streamwise and lower wall-normal stresses compared to the permeable-bed case. Bed roughness increases Reynolds shear stress whereas permeability has minimal influence. Bed permeability, in contrast to Reynolds stresses, significantly influences form-induced shear stress. Bed shear stress statistics show that probability of extreme events increases in permeable bed and rough wall cases as compared to the smooth wall. Findings suggest that bed permeability can have significant impact on modeling of hyporheic exchange and its effect on bed shear and pressure fluctuations are better captured by considering at least the top layer of sediment., 31 pages, 18 figures, journal paper

Accurate prediction of spray atomization process using an Euler-Lagrange (EL) approach is challenging because of high volume fraction of the liquid phase in dense regimes. This would in reality displace a remarkable portion of the gaseous phase which is commonly ignored in the standard EL approaches. In addition, deformation of droplet due to the interaction of aerodynamic force, surface tension and viscous forces is typically neglected in modeling dense sprays. In this work, to capture the volumetric displacement effects using an EL approach, the spatio-temporal changes in the volume fraction of the gaseous phase are taken into account. This leads to zero-Mach number, variable density equations that give rise to a source term in both momentum and continuity equations. It is shown that the continuity source term increases the velocity and dynamics of the carrier phase close to the nozzle. However, owing to the jet spread and dispersion of droplets, these effects decrease further downstream. In order to quantify the droplet deformation effects, different models are compared together with an experimental data. Different breakup regimes are studied in order to identify the best model for each regime. The shape deformation effect is isolated by performing a single droplet injected into the cross flow with flow conditions similar to the bag-type breakup. A significant deviation in the motion of droplet is observed compared to a case where deformation is neglected.

The accuracy of Euler-Lagrange point-particle models employed in particle-laden fluid flow simulations depends on accurate estimation of the particle force through closure models. Typical force closure models require computation of the slip velocity at the particle location, which in turn requires accurate estimation of the undisturbed fluid velocity. However, when the fluid and particle phases are two-way coupled the fluid velocity field is disturbed by the presence of the particle. A common practice is to use the disturbed velocity to compute the particle force which can result in errors as much as 100% in predicting the particle dynamics. In this work, a general velocity correction scheme is developed that facilitates accurate estimation of the undisturbed fluid velocity in particle-laden fluid flows with and without no-slip walls. The model can handle particles of different size, arbitrary interpolation functions, anisotropic grids with large aspect ratios, and wall-bounded flows. The present correction scheme is motivated by the recent work of Esmaily & Horwitz (JCP, 2018) on unbounded particle-laden flows. Modifications necessary for wall-bounded flows are developed such that the undisturbed fluid velocity at any wall distance is accurately recovered, asymptotically approaching the unbounded scheme for particles far away from walls. A detailed series of verification tests were conducted on settling velocity of a particle in parallel and perpendicular motions to a no-slip wall. A range of flow parameters and grid configurations; involving anisotropic grids with aspect ratios typically encountered in particle-laden turbulent channel flows, were considered in detail. When the wall effects are accounted for, the present correction scheme reduces the errors in predicting the near-wall particle motion by one order of magnitude smaller values compared to the unbounded correction schemes.

Accurate prediction of a dense spray using an Euler-Lagrange approach is challenging because of high volume fraction of the dispersed phase due to subgrid cluster of droplets. To accurately model dense sprays, one needs to capture this effect by taking into account the spatio-temporal changes in the volume fraction of the carrier phase due to the motion and presence of the dispersed phase. This leads to zero-Mach number, variable density equations which are commonly neglected in the standard two-way coupling spray simulations. Using pressure-based solvers, this gives rise to a source term in the pressure Poisson equation and a non-divergence free velocity field. To validate the predictive capability of such approach, an atomized non-evaporating dilute particulate round jet is first examined using Large Eddy Simulation coupled with Point-Particle approach and then higher volume loadings up to 38% are investigated with and without taking into account the volumetric displacement effects. It is shown that for volume loadings above 5%, the volumetric displacement effects enhance dynamics of the flow resulting in a higher stream-wise mean and r.m.s. velocities compared to the results of standard two-way coupling. This is more pronounced for the near field of the jet where local volume fraction of the dispersed phase is relatively high. This enhancement is conjectured to be due to the velocity divergence effect due to the modified continuity equation where spatio-temporal variations in volume fraction of the carrier phase increases velocity in the regions of high void fraction.

Modeling of a dense spray regime using an Euler-Lagrange approach is challenging because of local high volume loading. A cluster of droplets, that are assumed subgrid, can lead to locally low void fractions for the fluid phase. Under these conditions, spatio-temporal changes in the fluid volume fractions should be considered in an Euler-Lagrange, two-way coupling model. This leads to zero-Mach number, variable density governing equations. Using pressure-based solvers, this gives rise to a source term in the pressure Poisson equation and a non-divergence free velocity field. To test the validity and predictive capability of such an approach, a round jet laden with particles is investigated using Direct Numerical Simulation coupled with point-Particle based model and compared with available experimental data for a particulate turbulent round jet with $Re_j=5712$. Standard force closures including drag, lift, Magnus effect, pressure, added mass as well as viscous torque acting on each individual particle are employed in the Point-Particle based model. In addition, volume displacement effects due to the presence of solid particles or liquid droplets, which is commonly neglected in the standard two-way coupling, are taken into account in both continuity and inter-phase momentum transfer to accurately capture the underlying structure of particle-turbulence interactions. Prediction results are in well agreement with the corresponding experiment.