Computation of Flow Through a Three-Dimensional Periodic Array of Porous Structures by a Parallel Immersed-Boundary Method

2014 ◽  
Vol 136 (4) ◽  
Author(s):  
Zhenglun Alan Wei ◽  
Zhongquan Charlie Zheng ◽  
Xiaofan Yang
Keyword(s):  
Porous Media ◽  
Parallel Implementation ◽  
Flow Behavior ◽  
Three Dimensional ◽  
Flow Simulation ◽  
Periodic Array ◽  
Immersed Boundary ◽  
Navier Stokes Equation ◽  
Flow Through ◽  
Ib Method

A parallel implementation of an immersed-boundary (IB) method is presented for low Reynolds number flow simulations in a representative elementary volume (REV) of porous media that are composed of a periodic array of regularly arranged structures. The material of the structure in the REV can be solid (impermeable) or microporous (permeable). Flows both outside and inside the microporous media are computed simultaneously by using an IB method to solve a combination of the Navier–Stokes equation (outside the microporous medium) and the Zwikker–Kosten equation (inside the microporous medium). The numerical simulation is firstly validated using flow through the REVs of impermeable structures, including square rods, circular rods, cubes, and spheres. The resultant pressure gradient over the REVs is compared with analytical solutions of the Ergun equation or Darcy–Forchheimer law. The good agreements demonstrate the validity of the numerical method to simulate the macroscopic flow behavior in porous media. In addition, with the assistance of a scientific parallel computational library, PETSc, good parallel performances are achieved. Finally, the IB method is extended to simulate species transport by coupling with the REV flow simulation. The species sorption behaviors in an REV with impermeable/solid and permeable/microporous materials are then studied.

scholarly journals Direct Numerical Simulation of Gas–Particle Flows with Particle–Wall Collisions Using the Immersed Boundary Method

2018 ◽  
Vol 8 (12) ◽  
pp. 2387 ◽  
Author(s):  
Yusuke Mizuno ◽  
Shun Takahashi ◽  
Kota Fukuda ◽  
Shigeru Obayashi
Keyword(s):  
Immersed Boundary Method ◽  
Three Dimensional ◽  
Flow Simulation ◽  
Navier Stokes ◽  
Immersed Boundary ◽  
Navier Stokes Equation ◽  
Level Set Function ◽  
Boundary Method ◽  
The One ◽  
Energy And Momentum

We investigated particulate flows by coupling simulations of the three-dimensional incompressible Navier–Stokes equation with the immersed boundary method (IBM). The results obtained from the two-way coupled simulation were compared with those of the one-way simulation, which is generally applied for clarifying the particle kinematics in industry. In the present flow simulation, the IBM was solved using a ghost–cell approach and the particles and walls were defined by a level set function. Using proposed algorithms, particle–particle and particle–wall collisions were implemented simply; the subsequent coupling simulations were conducted stably. Additionally, the wake structures of the moving, colliding and rebounding particles were comprehensively compared with previous numerical and experimental results. In simulations of 50, 100, 200 and 500 particles, particle–wall collisions were more frequent in the one–way scheme than in the two-way scheme. This difference was linked to differences in losses in energy and momentum.

High-Order Immersed-Boundary Simulation and Error Analysis for Flow Around a Porous Structure

2017 ◽  
Author(s):  
Meihua Zhang ◽  
Z. Charlie Zheng
Keyword(s):  
Porous Structure ◽  
Three Dimensional ◽  
Time Discretization ◽  
Navier Stokes ◽  
Immersed Boundary ◽  
Dimensional Sphere ◽  
Navier Stokes Equation ◽  
Flow Resistivity ◽  
Fifth Order ◽  
Ib Method

In this study, a two-dimensional cylinder or a three-dimensional sphere is used as an example for a porous structure in the flow field. Immersed-boundary (IB) methods have been used to simulate flow around a cylinder/sphere, which is a typical problem to verify the effectiveness and accuracy of the methods. The flow was previously simulated by modeling the solid obstacles as a porous medium with flow resistivity. The current study makes an extension of the previous IB method. The fifth order WENO-Z scheme and third order Runge-Kutta method are used for space and time discretization, respectively. The Navier–Stokes equation is used for the fluids region and a ZK type of source term is applied in place of the forcing term for the porous domain. These two equations are solved simultaneously. Thus, there is no need to specify the interface conditions directly. In addition, the errors are estimated in terms of the flow resistivity that is zero for fluid, finite for porous, and infinite for non-porous solid. The aim of this paper is to establish estimates of the error induced by such an IB method. Numerical tests are performed to confirm the improvement of this method.

Fabrication of a Microchannel Device With a Three-Dimensional Pore Network Using a Sacrificial Sugar Template

2020 ◽  
Author(s):  
Haipeng Zhang ◽  
Tomer Palmon ◽  
Seunghee Kim ◽  
Sangjin Ryu
Keyword(s):  
Porous Media ◽  
Fluid Flow ◽  
Pore Structure ◽  
Three Dimensional ◽  
Pore Network ◽  
Compressed Air ◽  
Microfluidic Channels ◽  
Two Phase ◽  
Flow Through ◽  
Phase Fluid

Abstract Porous media compressed air energy storage (PM-CAES) is an emerging technology that stores compressed air in an underground aquifer during the off-peak periods, to mitigate the mismatch between energy supplies and demands. Thus, PM-CAES involves repeated two-phase fluid flow in porous media, and ensuring the success of PM-CAES requires a better understanding of repetitive two-phase fluid flow through porous media. For this purpose, we previously developed microfluidic channels that retain a two-dimensional (2D) pore network. Because it was found that the geometry of the pore structure significantly affects the patterns and occupational efficiencies of a non-wetting fluid during the drainage-imbibition cycles, a more realistic microfluidic model is needed to reflect the three-dimensional (3D) nature of pore structures in the underground geologic formation. In this study, we developed an easy-to-adopt method to fabricate a microfluidic device with a 3D random pore network using a sacrificial sugar template. Instead of using a master mold made in photolithography, a sacrificial mold was made using sugar grains so that the mold could be washed away after PDMS curing. First, we made sugar templates with different levels of compaction load, and found that the thickness of the templates decreased as the compaction load increased, which suggests more packing of sugar grains and thus lower porosity in the template. Second, we fabricated PDMS porous media using the sugar template as a mold, and imaged their pore structure using micro computed tomography (micro-CT). Pores within PDSM samples appeared more tightly packed as the compacting force increased. Last, we fabricated a prototype PDMS channel device with a 3D pore network using a sugar template, and visualized flow through the pore network using colored water. The flow visualization result shows that the water was guided by the random pores and that the resultant flow pattern was three dimensional.

Application of an Immersed Boundary Method to Flow in Cerebral Aneurysms and Porous Media

2010 ◽  
Author(s):  
Julia Mikhal ◽  
David J. Lopez Penha ◽  
Steffen Stolz ◽  
Bernard J. Geurts
Keyword(s):  
Poiseuille Flow ◽  
Stokes Equations ◽  
Spatial Arrangement ◽  
Computational Method ◽  
Implicit Time ◽  
Immersed Boundary ◽  
Time Stepping ◽  
Flow Through ◽  
Masking Function ◽  
Ib Method

We present the development and application of an immersed boundary (IB) method for the simulation of incompressible flow inside and around complex geometrical shapes and cavities. The IB method is based on a volume-penalization method that is applied throughout the domain, rendering the velocity in stationary solid parts negligibly small, while the flow in the open parts of the domain is governed by the Navier-Stokes equations. The flow solver is based on a skew-symmetric finite-volume discretization in combination with explicit time-stepping for the convective and viscous fluxes, and implicit time-stepping for the IB forcing term. The complex domain is characterized in terms of a so-called ‘masking function’ which equals unity in the solid parts and zero in the open parts of the domain. The focus is on the accuracy with which gradients of the solution close to solid walls can be approximated using the IB methodology. We investigate this for flow through a model of an aneurysm as may develop in the circle of Willis in a human brain, and to flow in a structured porous medium composed of a regular spatial arrangement of square rods. The shear stress acting on the vessel walls in case of flow through an aneurysm, and the permeability of the porous material, are analyzed. The computational method converges as a first order method for Poiseuille flow, with a considerable influence derived from the precise definition of the masking function near solid-fluid interfaces. We identify the best masking function strategy and show that for plane Poiseuille flow even second order convergence may be obtained. Qualitatively reliable results are obtained already at modest resolutions of 8–16 grid cells across a characteristic opening in the flow domain, e.g., the vessel diameter or the size of the gap between individual square rods.

Numerical Analysis of Flow Behavior Inside the Blast Furnace

2009 ◽  
Author(s):  
Dong Fu ◽  
Fengguo Tian ◽  
Guoheng Chen ◽  
D. Frank Huang ◽  
Chenn Q. Zhou
Keyword(s):  
Blast Furnace ◽  
Gas Flow ◽  
Cfd Simulation ◽  
Flow Behavior ◽  
Three Dimensional ◽  
Furnace Productivity ◽  
Gas Distribution ◽  
Furnace Shaft ◽  
Parametric Studies ◽  
Flow Through

Gas and burden distributions inside a blast furnace play an important role in optimizing gas utilization versus the furnace productivity and minimizing the CO2 emission in steel industries. In this paper, a mathematical model is presented to describe the burden descent in the blast furnace shaft and gas distribution, with the alternative structure of coke and ore layers being considered. Multi-dimensional Ergun’s equation is solved with considering the turbulent compressible gas flow through the burden column. The porosity of each material will be treated as a function of three dimensional functions which will be determined by the kinetics sub-models accordingly. A detailed investigation of gas flow through the blast furnace will be conducted with the given initial burden profiles along with the effects of redistribution during burden descending. Also, parametric studies will be carried out to analyze the gas distribution cross the blast furnace under different cohesive zone (CZ) shapes, charging rate, and furnace top pressure. A good agreement was obtained between the CFD simulation and published experimental data. Based on the results, the inverse V shape is proved to be the most desirable CZ profile.

Sign in / Sign up

Export Citation Format

Share Document