This paper introduces a new weighting scheme for particle-grid transfers that generates hybrid Lagrangian/Eulerian fluid simulations with uniform particle distributions and precise volume control. At its core, our approach reformulates the construction of Power Particles [de Goes et al. 2015] by computing volume-constrained density kernels. We employ these optimized kernels as particle domains within the Generalized Interpolation Material Point method (GIMP) in order to incorporate Power Particles into the Particle-In-Cell framework, hence the name the Power Particle-In-Cell method. We address the construction of volume-constrained density kernels as a regularized optimal transportation problem and describe an iterative solver based on localized Gaussian convolutions that leads to a significant performance speedup compared to [de Goes et al. 2015]. We also present novel extensions for handling free surfaces and solid obstacles that bypass the need for cell clipping and ghost particles. We demonstrate the advantages of our transfer weights by improving hybrid schemes for fluid simulation such as the Fluid Implicit Particle (FLIP) method and the Affine Particle-In-Cell (APIC) method with volume preservation and robustness to varying particle-per-cell ratio, while retaining low numerical dissipation, conserving linear and angular momenta, and avoiding particle reseeding or post-process relaxations.
more »
« less
Dense Discrete Phase Model Coupled with Kinetic Theroy of Granular Flow to Improve Predictions of Bubbling Fluidized Bed Hydrodynamics
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
- Award ID(s):
- 1736173
- NSF-PAR ID:
- 10085263
- Date Published:
- Journal Name:
- Kona powder and particle journal
- ISSN:
- 2187-5537
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
We present the Moving Eulerian-Lagrangian Particles (MELP), a novel mesh-free method for simulating incompressible fluid on thin films and foams. Employing a bi-layer particle structure, MELP jointly simulates detailed, vigorous flow and large surface deformation at high stability and efficiency. In addition, we design multi-MELP: a mechanism that facilitates the physically-based interaction between multiple MELP systems, to simulate bubble clusters and foams with non-manifold topological evolution. We showcase the efficacy of our method with a broad range of challenging thin film phenomena, including the Rayleigh-Taylor instability across double-bubbles, foam fragmentation with rim surface tension, recovery of the Plateau borders, Newton black films, as well as cyclones on bubble clusters.more » « less
-
more » « less
Nonspherical particles are commonly found when processing biomass or municipal solid waste. In this study, cylindrical particles are used as generic nonspherical particles and are co‐fluidized with small spherical particles. X‐ray particle tracking velocimetry is used to track the three‐dimensional particle position and velocity of a single tagged cylindrical particle over a long time period in the binary fluidized bed. The effects of superficial gas velocity (
u f ), cylindrical particle mass fraction (α ), particle sphericity (Φ ), and bed material size on the cylindrical tracer particle location and velocity are investigated. Overall, the cylindrical particles are found in the near‐wall region more often than in the bed center region. Increasing the superficial gas velocityu f provide a slight improvement in the uniformity of the vertical and horizontal distributions. Increasing the cylindrical particle mass fractionα causes the bed mixing conditions to transition from complete mixing into partial mixing. © 2018 American Institute of Chemical EngineersAIChE J , 65: 520–535, 2019 -
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.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- The DOI is not currently available.
