Flow past disperse solid particles or bubbles induces fluctuations in carrier fluid velocity, which correlate with temperature fluctuations in non-isothermal flows resulting in the pseudo-turbulent heat flux (PTHF). In the Eulerian-Eulerian (EE) two-fluid (TF) model, the transport of PTHF is shown to be an important contributor to the overall energy budget, and is modeled using a pseudo-turbulent thermal diffusivity (PTTD). The PTHF and PTTD were originally quantified using particle-resolved direct numerical simulation (PR-DNS) data, and correlations were developed over a range of solid volume fraction (0.1 ≤ 𝜀𝑠 ≤ 0.5) and mean slip Reynolds number (1 ≤ 𝑅𝑒𝑚 ≤ 100) for a Prandtl number of 0.7. However, the original PTTD correlation diverges to infinity as the solid volume fraction goes to zero, which is physically unrealistic. This singular behavior is problematic for EE TF simulations at particle material fronts where solid volume fraction values can fall below the lower limit of existing data (𝜀𝑠 =0.1) to zero in the pure carrier phase. In this work, additional PR-DNS data are reported for 𝜀𝑠 < 0.1, and improved correlations are developed for the PTHF and PTTD. The new PTTD correlation is non- singular, and both the PTHF and PTTD decay exponentially to zero as the solid volume fraction approaches zero, which is physically reasonable. This improves prediction of PTHF transport in dilute flow using EE TF heat transfer simulations.
more »
« less
Eulerian-Eulerian Description of the Interaction of a Shock With Particles Through Godunov’s Scheme
This paper is on an Eulerian-Eulerian (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 sub-grid scale stresses, and particle sub-grid 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 two-dimensional 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 well-stablished Eulerian-Lagragian (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 the EL results with a very good accuracy.
more »
« less
- Award ID(s):
- 1115631
- PAR ID:
- 10090027
- Date Published:
- Journal Name:
- The ASME 2013 Fluids Engineering Division Summer Meeting
- Page Range / eLocation ID:
- V01CT25A006
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
We propose a novel Particle Flow Map (PFM) method to enable accurate long-range advection for incompressible fluid simulation. The foundation of our method is the observation that a particle trajectory generated in a forward simulation naturally embodies a perfect flow map. Centered on this concept, we have developed an Eulerian-Lagrangian framework comprising four essential components: Lagrangian particles for a natural and precise representation of bidirectional flow maps; a dual-scale map representation to accommodate the mapping of various flow quantities; a particle-to-grid interpolation scheme for accurate quantity transfer from particles to grid nodes; and a hybrid impulse-based solver to enforce incompressibility on the grid. The efficacy of PFM has been demonstrated through various simulation scenarios, highlighting the evolution of complex vortical structures and the details of turbulent flows. Notably, compared to NFM, PFM reduces computing time by up to 49 times and memory consumption by up to 41%, while enhancing vorticity preservation as evidenced in various tests like leapfrog, vortex tube, and turbulent flow.more » « less
-
Summary In this paper, a three‐dimensional numerical solver is developed for suspensions of rigid and soft particles and droplets in viscoelastic and elastoviscoplastic (EVP) fluids. The presented algorithm is designed to allow for the first time three‐dimensional simulations of inertial and turbulent EVP fluids with a large number particles and droplets. This is achieved by combining fast and highly scalable methods such as an FFT‐based pressure solver, with the evolution equation for non‐Newtonian (including EVP) stresses. In this flexible computational framework, the fluid can be modeled by either Oldroyd‐B, neo‐Hookean, FENE‐P, or Saramito EVP models, and the additional equations for the non‐Newtonian stresses are fully coupled with the flow. The rigid particles are discretized on a moving Lagrangian grid, whereas the flow equations are solved on a fixed Eulerian grid. The solid particles are represented by an immersed boundary method with a computationally efficient direct forcing method, allowing simulations of a large numbers of particles. The immersed boundary force is computed at the particle surface and then included in the momentum equations as a body force. The droplets and soft particles on the other hand are simulated in a fully Eulerian framework, the former with a level‐set method to capture the moving interface and the latter with an indicator function. The solver is first validated for various benchmark single‐phase and two‐phase EVP flow problems through comparison with data from the literature. Finally, we present new results on the dynamics of a buoyancy‐driven drop in an EVP fluid.more » « less
-
Abstract In this paper, a hybrid Lagrangian–Eulerian topology optimization (LETO) method is proposed to solve the elastic force equilibrium with the Material Point Method (MPM). LETO transfers density information from freely movable Lagrangian carrier particles to a fixed set of Eulerian quadrature points. This transfer is based on a smooth radial kernel involved in the compliance objective to avoid the artificial checkerboard pattern. The quadrature points act as MPM particles embedded in a lower‐resolution grid and enable a subcell multidensity resolution of intricate structures with a reduced computational cost. A quadrature‐level connectivity graph‐based method is adopted to avoid the artificial checkerboard issues commonly existing in multiresolution topology optimization methods. Numerical experiments are provided to demonstrate the efficacy of the proposed approach.more » « less
-
Formation, expansion, and breakage of bubbles in single bubble and freely bubbling fluidized beds were studied using an improved hybrid Lagrangian-Eulerian computational fluid dynamics (CFO) approach. Dense Discrete Phase Model (DDPM) is a novel approach to simulate industrial scale fluidized bed reactors with polydispersed particles. The model uses a hybrid Lagrangian-Eulerian approach to track the particle parcels (lumping several particles in one computational cell) in a Lagrangian framework according to Newton's laws of motion. The interactions between particles are estimated by the gradient of solids stress solved in Eulerian grid. In this work. a single bubble fluidized bed and a freely bubbling fluidized bed were simulated using DDPM coupled with kinetic theory of granular flows (KTGF). The solid stress was improved to include both tangential and normal forces compared to current hybrid methods with the consideration of only normal stress or solid pressure. The results showed that solid pressure (normal forces) as the only contributor in solid stress would lead to over prediction of bubble size and overlooking of bubble breakage in a single bubble bed. Also, the results showed the improved model bad a good prediction of bubble path in a freely bubbling bed compared to solid pressure-based model. It was shown that increasing the restitution coefficient increased the particle content of the bubbles and it lead to less breakage during the formation of the bubble. The probability of formation of bubbles was compared with experimental results and solid stress model showed less discrepancies compared to the solid pressure-based model.more » « less