A Sharp Interface Direct Forcing Immersed Boundary Approach for Fully Resolved Simulations of Particulate Flows

2014 ◽  
Vol 136 (4) ◽  
Author(s):  
Jianming Yang ◽  
Frederick Stern
Keyword(s):  
Volume Fraction ◽  
Solid Volume Fraction ◽  
Density Ratio ◽  
Sharp Interface ◽  
Immersed Boundary ◽  
Particulate Flows ◽  
Order Accuracy ◽  
Wide Range ◽  
Direct Forcing ◽  
Coarse Grids

In recent years, the immersed boundary method has been well received as an effective approach for the fully resolved simulations of particulate flows. Most immersed boundary approaches for numerical studies of particulate flows in the literature were based on various discrete delta functions for information transfer between the Lagrangian elements of an immersed object and the underlying Eulerian grid. These approaches have some inherent limitations that restrict their wider applications. In this paper, a sharp interface direct forcing immersed boundary approach based on the method proposed by Yang and Stern (Yang and Stern, 2012, “A Simple and Efficient Direct Forcing Immersed Boundary Framework for Fluid-Structure Interactions,” J. Comput. Phys., 231(15), pp. 5029–5061) is given for the fully resolved simulations of particulate flows. This method uses a discrete forcing approach and maintains a sharp profile of the fluid-solid interface. It is not limited to low Reynolds number flows and the immersed boundary discretization can be arbitrary or totally eliminated for particles with analytical shapes. In addition, it is not required to calculate the solid volume fraction in low density ratio problems. A strong coupling scheme is employed for the fluid-solid interaction without including the fluid solver in the predictor-corrector iterative loop. The overall algorithm is highly efficient and very attractive for simulating particulate flows with a wide range of density ratios on relatively coarse grids. Several cases are examined and the results are compared with reference data to demonstrate the simplicity and robustness of our method in particulate flow simulations. These cases include settling and buoyant particles and the interaction of two settling particles showing the kissing-drafting-tumbling phenomenon. Systematic verification studies show that our method is of second-order accuracy on very coarse grids and approaches fourth-order accuracy on finer grids.

An Efficient Immersed Boundary Method Based on Penalized Direct Forcing for Simulating Flows Through Real Porous Media

2012 ◽  
Author(s):  
Wim-Paul Breugem ◽  
Vincent van Dijk ◽  
René Delfos
Keyword(s):  
Porous Medium ◽  
Ct Scan ◽  
Immersed Boundary Method ◽  
Volume Fraction ◽  
Solid Volume Fraction ◽  
Particle Diameter ◽  
Immersed Boundary ◽  
Average Particle ◽  
Boundary Method ◽  
Direct Forcing

A computationally efficient Immersed Boundary Method (IBM) based on penalized direct forcing was employed to determine the permeability of a real porous medium. The porous medium was composed of about 9000 glass beads with an average particle diameter of 1.93 mm and a porosity of 0.367. The forcing of the IBM depends on the local solid volume fraction within a computational grid cell. The latter could be obtained from a high-resolution X-ray Computed Tomography (CT) scan of the packing. An experimental facility was built to determine the permeability of the packing experimentally. Numerical simulations were performed for the same packing based on the data from the CT scan. For a scan resolution of 0.1 mm the numerical value for the permeability was nearly 70% larger than the experimental value. An error analysis indicated that the scan resolution of 0.1 mm was too coarse for this packing.

Comparison of Sharp Interface and Smoothed Profile Methods for Laminar Flow Analysis Over Stationary and Moving Boundaries

2014 ◽  
Author(s):  
Fazlolah Mohaghegh ◽  
John Mousel ◽  
H. S. Udaykumar
Keyword(s):  
Drag Coefficient ◽  
Immersed Boundary Method ◽  
Flow Analysis ◽  
Sharp Interface ◽  
Wake Flow ◽  
Immersed Boundary ◽  
Particulate Flows ◽  
Sharp Interface Method ◽  
Smoothed Profile Method ◽  
Coarse Grids

This study is a comparison of two techniques for simulation of particulate flows on fixed Cartesian grids: Sharp interface Method (SIM) (Udaykumar et al., 2001, 2002, 2003) and a modified version of Immersed Boundary Method (Peskin, 1977) (IBM) known as Smoothed Profile Method (SPM) (Nakayama and Yamamoto, 2005; Luo et. al, 2009). Different cases were studied includes flow over one or two moving and stationary particles. Predictions of the drag coefficient shows that SPM and SIM are very close to the experiments. SIM slightly under-predicts the value of the drag coefficient while SPM has a small over-estimation. Moreover, SPM is more accurate on coarse grids. However, with refinement of the grid SIM approaches the exact values very fast leading to better results on fine grids. Flow pattern and vortex structures of SPM and SIM are almost the same. Both methods are capable of analyzing the wake flow. Unlike SIM, SPM is able to simulate the flow when two particles are in contact. When two particles are in motion and are very close in a way that the two interfaces overlap, SPM shows a repulsion force between two spheres which reduces the accuracy in comparison with SIM. However, SPM can achieve the collision of two particles without problem.

Computational Study of Slurry Flow in Pipes to Determine Profiles Sensed in Near Infrared Slurry Measurements

2011 ◽  
Author(s):  
V. Pasangulapati ◽  
N. R. Kesana ◽  
G. Sharma ◽  
F. W. Chambers ◽  
M. E. McNally ◽  
...  
Keyword(s):  
Near Infrared ◽  
Volume Fraction ◽  
Solid Phase ◽  
Computational Study ◽  
Solid Volume Fraction ◽  
Density Ratio ◽  
Optical Measurements ◽  
Concentration Profiles ◽  
Horizontal Pipe ◽  
Volume Fractions

It is desired to perform accurate Near Infrared sensor measurements of slurries flowing in pipes leaving large batch reactors. A concern with these measurements is the degree to which the slurry sensed is representative of the material in the reactor and flowing through the pipe. Computational Fluid Dynamics (CFD) has been applied to the flow in the pipe to determine the flow fields and the concentration profiles seen by the sensors. The slurry was comprised of a xylene liquid phase and an ADP (2-amino-4, 6-dimethylpyrimidine) solid phase with a density ratio of 1.7. Computations were performed for a horizontal pipe with diameter 50.8 mm, length 2.032 m, and 1.76 m/s and 3.26 m/s mixture velocities. The corresponding pipe Reynolds numbers were 1.19E+05 and 2.21E+05. The flow through a slotted cylindrical probe inserted radially in the pipe also was considered. Spherical slurry particles with diameters from 10 μm to 1000 μm were considered with solid volume fractions of 12%, 24%, and 35%. Computations were performed with ANSYS FLUENT 12 software using the Realizable k-ε turbulence model and the enhanced wall treatment function. Comparisons of computed vertical profiles of solid volume fraction to results in the literature showed good agreement. Symmetric, nearly flat solid volume fraction profiles were observed for 38 μm particles for all three initial solid volume fractions. Asymmetric solid volume fraction profiles with greater values toward the bottom were observed for the larger particles. Changes in the profiles of turbulent kinetic energy also were observed. These changes are important for optical measurements which depend upon the mean concentration profiles as well as the turbulent motion of the slurry particles.

Flows Through Real Porous Media: X-Ray Computed Tomography, Experiments, and Numerical Simulations

2014 ◽  
Vol 136 (4) ◽  
Author(s):  
Wim-Paul Breugem ◽  
Vincent van Dijk ◽  
René Delfos
Keyword(s):  
Computed Tomography ◽  
Porous Medium ◽  
Ct Scan ◽  
Numerical Simulations ◽  
Volume Fraction ◽  
Solid Volume Fraction ◽  
Immersed Boundary ◽  
Immersed Boundary Methods ◽  
X Ray ◽  
X Ray Computed

Two different direct-forcing immersed boundary methods (IBMs) were applied for the purpose of simulating slow flow through a real porous medium: the volume penalization IBM and the stress IBM. The porous medium was a random close packing of about 9000 glass beads in a round tube. The packing geometry was determined from an X-ray computed tomography (CT) scan in terms of the distribution of the truncated solid volume fraction (either 0 or 1) on a three-dimensional Cartesian grid. The scan resolution corresponded to 19.3 grid cells over the mean bead diameter. A facility was built to experimentally determine the permeability of the packing. Numerical simulations were performed for the same packing based on the CT scan data. For both IBMs the numerically determined permeability based on the Richardson extrapolation was just 10% lower than the experimentally found value. As expected, at finite grid resolution the stress IBM appeared to be the most accurate IBM.

Immersed Boundary – Thermal Lattice Boltzmann Methods for Non-Newtonian Flows Over a Heated Cylinder: A Comparative Study

2015 ◽  
Vol 18 (2) ◽  
pp. 489-515 ◽  
Author(s):  
A. Amiri Delouei ◽  
M. Nazari ◽  
M. H. Kayhani ◽  
S. Succi
Keyword(s):  
Lattice Boltzmann ◽  
Computational Cost ◽  
Shear Thickening ◽  
The Body ◽  
Sharp Interface ◽  
Immersed Boundary ◽  
Heated Cylinder ◽  
Hydrodynamic Approach ◽  
Direct Forcing ◽  
Thermal Lattice Boltzmann

AbstractIn this study, we compare different diffuse and sharp interface schemes of direct-forcing immersed boundary — thermal lattice Boltzmann method (IB-TLBM) for non-Newtonian flow over a heated circular cylinder. Both effects of the discrete lattice and the body force on the momentum and energy equations are considered, by applying the split-forcing Lattice Boltzmann equations. A new technique based on predetermined parameters of direct forcing IB-TLBM is presented for computing the Nusselt number. The study covers both steady and unsteady regimes (20<Re<80) in the power-law index range of 0.6<n<1.4, encompassing both shear-thinning and shear-thickening non-Newtonian fluids. The numerical scheme, hydrodynamic approach and thermal parameters of different interface schemes are compared in both steady and unsteady cases. It is found that the sharp interface scheme is a suitable and possibly competitive method for thermal-IBM in terms of accuracy and computational cost.

Effective Geometric Algorithms for Immersed Boundary Method Using Signed Distance Field

2018 ◽  
Vol 141 (6) ◽  
Author(s):  
Chenguang Zhang ◽  
Chunliang Wu ◽  
Krishnaswamy Nandakumar
Keyword(s):  
Immersed Boundary Method ◽  
Volume Fraction ◽  
Solid Volume Fraction ◽  
Closed Surface ◽  
Immersed Boundary ◽  
Field Calculation ◽  
Distance Field ◽  
Signed Distance ◽  
Boundary Method ◽  
Signed Distance Field

We present three algorithms for robust and efficient geometric calculations in the context of immersed boundary method (IBM), including classification of mesh cells as inside/outside of a closed surface, projection of points onto a surface, and accurate calculation of the solid volume fraction field created by a closed surface overlapping with a background Cartesian mesh. The algorithms use the signed distance field (SDF) to represent the surface and remove the intersection tests, which are usually required by other algorithms developed before, no matter the surface is described in analytic or discrete form. The errors of the algorithms are analyzed. We also develop an approximate method on efficient SDF field calculation for complex geometries. We demonstrate how the algorithms can be implemented within the framework of IBM with a volume-average discrete-forcing scheme and applied to simulate fluid–structure interaction problems.

Stochastic Lagrangian model for hydrodynamic acceleration of inertial particles in gas–solid suspensions

2016 ◽  
Vol 788 ◽  
pp. 695-729 ◽  
Author(s):  
Sudheer Tenneti ◽  
Mohammad Mehrabadi ◽  
Shankar Subramaniam
Keyword(s):  
Particle Acceleration ◽  
Particle Velocity ◽  
Volume Fraction ◽  
Solid Volume Fraction ◽  
Density Ratio ◽  
Granular Temperature ◽  
Solid Flow ◽  
Acceleration Model ◽  
The Mean ◽  
Solid Suspensions

The acceleration of an inertial particle in a gas–solid flow arises from the particle’s interaction with the gas and from interparticle interactions such as collisions. Analytical treatments to derive a particle acceleration model are difficult outside the Stokes flow regime, but for moderate Reynolds numbers (based on the mean slip velocity between gas and particles) particle-resolved direct numerical simulation (PR-DNS) is a viable tool for model development. In this study, PR-DNS of freely-evolving gas–solid suspensions are performed using the particle-resolved uncontaminated-fluid reconcilable immersed-boundary method (PUReIBM) that has been extensively validated in previous studies. Analysis of the particle velocity variance (granular temperature) equation in statistically homogeneous gas–solid flow shows that a straightforward extension of a class of mean particle acceleration models (drag laws) to their corresponding instantaneous versions, by replacing the mean particle velocity with the instantaneous particle velocity, predicts a granular temperature that decays to zero, which is at variance with the steady particle granular temperature that is obtained from PR-DNS. Fluctuations in particle velocity and particle acceleration (and their correlation) are important because the particle acceleration–velocity covariance governs the evolution of the particle velocity variance (characterized by the particle granular temperature), which plays an important role in the prediction of the core annular structure in riser flows. The acceleration–velocity covariance arising from hydrodynamic forces can be decomposed into source and dissipation terms that appear in the granular temperature evolution equation, and these have already been quantified in the Stokes flow regime using a combination of kinetic theory closure and multipole expansion simulations. From PR-DNS data we show that the fluctuations in the particle acceleration that are aligned with fluctuations in the particle velocity give rise to a source term in the granular temperature evolution equation. This approach is used to quantify the hydrodynamic source and dissipation terms of granular temperature from PR-DNS results for freely-evolving gas–solid suspensions that are performed over a wide range of solid volume fraction ($0.1\leqslant {\it\phi}\leqslant 0.4$), Reynolds number based on the slip velocity between the solid and the fluid phase ($10\leqslant \mathit{Re}_{m}\leqslant 100$) and solid-to-fluid density ratio ($100\leqslant {\it\rho}_{p}/{\it\rho}_{f}\leqslant 2000$). The straightforward extension of drag law models does not give rise to any source in the granular temperature due to hydrodynamic effects. This motivates the development of better Lagrangian particle acceleration models that can be used in Lagrangian–Eulerian formulations of gas–solid flow. It is found that a Langevin equation for the increment in the particle velocity reproduces PR-DNS results for the stationary particle velocity autocorrelation in freely-evolving suspensions. Based on the data obtained from the simulations, the functional dependence of the Langevin model coefficients on solid volume fraction, Reynolds number and solid-to-fluid density ratio is obtained. This new Lagrangian particle acceleration model reproduces the correct steady granular temperature and can also be adapted to gas–solid flow computations using Eulerian moment equations.

Fully Resolved Simulation of Particulate Flow Using a Sharp Interface Direct Forcing Immersed Boundary Method

2013 ◽  
Author(s):  
Jianming Yang ◽  
Frederick Stern
Keyword(s):  
Immersed Boundary Method ◽  
Solid Particles ◽  
Sharp Interface ◽  
Particulate Flow ◽  
Immersed Boundary ◽  
Fluid Structure ◽  
Boundary Method ◽  
Direct Forcing ◽  
Fully Resolved Simulation ◽  
Predictor Corrector

In this paper, the sharp interface, direct forcing immersed boundary method developed by Yang and Stern (A simple and efficient direct forcing immersed boundary framework for fluid-structure interactions, J. Comput. Phys. 231 (2012) 5029–5061) is applied to the fully resolved simulation of particulate flow. This method uses a discrete forcing approach and maintains a sharp profile of the fluid/solid interface. Also, it employs a strong coupling scheme for fluid-structure interaction through a predictor-corrector algorithm. The fluid flow solver is not included in the predictor-corrector iterative loop thanks to the direct forcing idea, which makes the overall algorithm highly efficient and very attractive for the fully resolved simulation of particulate flow with numerous solid particles. Several cases including sedimenting and buoyant particles and the interaction of two sedimenting particles showing kissing, drafting, and tumbling phenomenon are examined and compared with the reference results to demonstrate the simplicity and applicability of our method in particulate flow simulations.

DNS Study of Collision and Coalescence Over a Wide Range of Volume Fraction

2010 ◽  
Author(s):  
G. Luret ◽  
T. Me´nard ◽  
J. Re´veillon ◽  
A. Berlemont ◽  
F. X. Demoulin
Keyword(s):  
Kinetic Energy ◽  
Weber Number ◽  
Volume Fraction ◽  
Density Ratio ◽  
Liquid Volume ◽  
Kinetic Energy Density ◽  
Two Phase ◽  
Liquid Volume Fraction ◽  
Gas Phases ◽  
Wide Range

Among the different processes that play a role during the atomization process, collisions are addressed in this work. Collisions can be very important in dense two-phase flows. Recently, the Eulerian Lagrangian Spray Atomization (ELSA) model has been developed. It represents the atomization by taking into account the dense zone of the spray. Thus in this context, collisions modeling are of the utmost importance. In this model results of collisions are controlled by the value of an equilibrium Weber number, We*. It is defined as the ratio between the kinetic energy to the surface energy. Such a value of We* has been studied in the past using Lagrangian collision models with various complexity. These models are based on analysis of collisions between droplets that have surface at rest. This ideal situation can be obtained only if droplet agitation created during a collision has enough time to vanish before the next collision. For a spray, this requirement is not always fulfill depending for instance on the mean liquid volume fraction. If there is not enough time, collisions will occur between agitated droplets changing the issue of the collision with respect to the ideal case. To study this effect, a DNS simulation with a stationary turbulence levels has been conducted for different liquid volume fractions in a cubic box with periodic condition in all directions. For liquid volume fraction close to zero the spray is diluted and collisions between spherical droplets can be identified. For a volume fraction close to one, collisions between bubbles are found. For a middle value of the volume fraction no discrete phase can be observed, instead a strong interaction between both liquid and gas phases is taking place. In all this case the equilibrium value of the Weber number We* can be determined. First propositions to determine We* as a function of the kinetic energy, density ratio, surface tension coefficient and the volume fraction will be proposed.

Sign in / Sign up

Export Citation Format

Share Document