Simulation of flow through porous anode in MFC at higher power density

Microbial fuel cell (MFC) is a novel environmental friendly energy device which has received great attention due to its technology for producing electricity directly from organic or inorganic matter by using bacteria as catalyst. To date, many experiments have been carried out to achieve the maximum power output with advective flow through porous anode to the cathode in the MFC. However, the precise mechanical mechanism of flow through anode and the quantified relationship between electrode spacing and MFC performance are not yet clearly understood. It has been found experimentally that the power output can be increased apparently at certain electrode spacing configuration. Based on these available experimental data, this paper investigates the effect of spacing between electrodes, the Darcy number of porous anode and the Reynolds number on the power production performance of MFC by using lattice Boltzmann method. The numerical simulation results present that the distance between electrodes significantly influences the flow velocity and residence time of the organic matter attached to the anode in the MFC. Moreover, it is found that the Darcy number of porous anode and the Reynolds number can regulate the output efficiency of MFC. These results perform better understanding of the complex phenomena of MFC and will be helpful to optimize MFC design.     

A unified mesh refinement method with applications to porous media flow

A unified algorithm is presented for the refinement of finite element meshes consisting of tensor product Lagrange elements in any number of space dimensions. The method leads to repeatedly refined n-irregular grids with associated constraint equations. Through an object-oriented implementation existing solvers can be extended to handle mesh refinements without modifying the implementation of the finite element equations. Various versions of the refinement procedure are investigated in a porous media flow problem involving singularities around wells. A domain decomposition-type finite element method is also proposed based on the refinement technique. This method is applied to flow in heterogeneous porous media. © 1998 John Wiley & Sons, Ltd.     

A hybrid FVM–LBM method for single and multi‐fluid compressible flow problems

The lattice Boltzmann method (LBM) has established itself as an alternative approach to solve the fluid flow equations. In this work we combine LBM with the conventional finite volume method (FVM), and propose a non‐iterative hybrid method for the simulation of compressible flows. LBM is used to calculate the inter‐cell face fluxes and FVM is used to calculate the node parameters. The hybrid method is benchmarked for several one‐dimensional and two‐dimensional test cases. The results obtained by the hybrid method show a steeper and more accurate shock profile as compared with the results obtained by the widely used Godunov scheme or by a representative flux vector splitting scheme. Additional features of the proposed scheme are that it can be implemented on a non‐uniform grid, study of multi‐fluid problems is possible, and it is easily extendable to multi‐dimensions. These features have been demonstrated in this work. The proposed method is therefore robust and can possibly be applied to a variety of compressible flow situations. Copyright © 2009 John Wiley & Sons, Ltd.     

缝洞型多孔介质中多相流的有限体积法数值模拟

本文研究的碳酸盐岩油藏储集体属于缝洞型多孔介质.这类缝洞型多孔介质由裂缝、溶蚀孔洞和低孔隙度低渗透率的基岩组成.裂缝是空隙流体流动的主要通道;溶蚀孔洞大小从几厘米到数米不等,渗透率和孔隙度都很高,是流体主要的储集空间.由于缝洞型多孔介质空隙空间的复杂性和强非均质性,数值计算中基本控制方程的空间离散应采用非结构化网格的计算模型.本文采用有限体积法模拟缝洞型多孔介质中多相流体的流动,并给出了相应的单元中心格式有限体积法的计算公式.裂缝介质和溶洞介质中单元间多相流体的流动考虑为高速非达西流,其质量通量采用Forchheimer定律计算.非线性方程的离散选取全隐式格式,并采用Newton-Raphson迭代进行求解.通过两个二维模型注水驱油的数值模拟,验证了本文方法的有效性.     

The approximate solution of the self-similar problem for radial flow of non-newtonian fluids through porous media

In this paper, we study the approximate solution of the self-simikar problem for radial flow of non-Newtonian fluids through porous media. Assuming that the fluids obey the exponential function law, we obtain an exact solution for the exponent n=0 and compare it with the approximate solution in ref. [1]. For n>1 and n<1, we obtain respectively approximate solutions. Some exampls are presented.     

Determination of permeability in fibrous porous media using the lattice Boltzmann method with application to PEM fuel cells

The lattice Boltzmann method (LBM) is used to simulate the flow through an idealized proton exchange membrane fuel cell (PEMFC) porous transport layer (PTL) geometry generated using a Monte Carlo method. Using the calculated flow field, Darcy's law is applied and the permeability is calculated. This process is applied in both through‐ and in‐plane directions of the paper as both of these permeability values are important in computational fluid dynamics models of PEMFCs. It is shown that the LBM can be used to determine permeability in a random porous media by solving the flow in the microstructure of the material. The permeability in the through‐ and in‐plane directions is shown to be different and the anisotropic nature of the geometry creates anisotropic permeability. It is also found that fiber arrangement plays a large role in the permeability of the PTL. New correlations are presented for in‐ and though‐plane permeabilities of fibrous porous media with (0.6<ε<0.8). Copyright © 2008 John Wiley & Sons, Ltd.     

Application of the lattice‐Boltzmann method to study flow and dispersion in channels with and without expansion and contraction geometry

The lattice‐Boltzmann (LB) method, derived from lattice gas automata, is a relatively new technique for studying transport problems. The LB method is investigated for its accuracy to study fluid dynamics and dispersion problems. Two problems of relevance to flow and dispersion in porous media are addressed: (i) Poiseuille flow between parallel plates (which is analogous to flow in pore throats in two‐dimensional porous networks), and (ii) flow through an expansion–contraction geometry (which is analogous to flow in pore bodies in two‐dimensional porous networks). The results obtained from the LB simulations are compared with analytical solutions when available, and with solutions obtained from a finite element code (FIDAP) when analytical results are not available. Excellent agreement is found between the LB results and the analytical/FIDAP solutions in most cases, indicating the utility of the lattice‐Boltzmann method for solving fluid dynamics and dispersion problems. Copyright © 1999 John Wiley & Sons, Ltd.     

LBM simulation of fluid flow in fractured porous media with permeable matrix

To analyze and depict complicated fluid behaviors in fractured porous media with variably permeable matrix, an integrated discrete computational algorithm is proposed based on lattice Boltzmann method (LBM). This paper combines with the external force model and statistical material physics to effectively describe the feature changes while the fluid passes through the fractures within the permeable matrix. As an application example, a two dimensional rock sample is reconstructed using the digital image and characterized with different feature values at each LBM grid to distinguish pores, impermeable and permeable matrix by stating its local physical property. Compared with the conventional LBM, the results demonstrate the advantages of proposed algorithm in modeling fluid flow phenomenon in fractured porous media with variably permeable matrix.     

Numerical simulation of the flow around a porous covering square cylinder in a channel via lattice Boltzmann method

In this paper, a detailed investigation on the flow past a porous covering cylinder is presented through the lattice Boltzmann method. The Brinkman‐Forchheimer‐extended Darcy model is adopted for the entire flow field with the solid, fluid, and porous medium. The effects of several parameters, such as porous layer thickness, Darcy number, porosity, and Reynolds number on flow field are discussed. Compared with the case of a solid cylinder, the present work shows that the porous layer may play an important role on the flow, the lift and drag force exerted on the cylinder. The numerical results indicate that the maximal drag coefficient C_d and maximal amplitude of lift coefficient C_l exist at certain Darcy number which is in the range of 10^?6–10^?2. Copyright © 2010 John Wiley & Sons, Ltd.     

New studying of Lattice Boltzmann method for two-phase driven in porous media

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｔｗｏ＿ｐｈａｓｅｄｒｉｖｅｎｉｎｐｏｒｏｕｓｍｅｄｉａｉｓａｎｉｍｐｏｒｔａｎｔｓｕｂｊｅｃｔｏｆｗｉｄｅｉｎｔｅｒｄｉｓｃｉｐｌｉｎａｒｙｃｏｎｃｅｒｎ ,ｓｕｃｈａｓｐｅｔｒｏｌｅｕｍｉｎｄｕｓｔｒｙ ,ｓｕｂｔｅｒｒａｎｅａｎ ,ｃｈｅｍｉｃａｌｐｒｏｊｅｃｔｅｔｃ..Ｉｔ’ｓｖｅｒｙｄｉｆｆｉｃｕｌｔｉｆｗｅａｄｏｐｔｎｕｍｅｒｉｃａｌｉｎｔｅｇｒａｌｍｅｔｈｏｄｆｏｒｔｈｅｃａｓｅｏｆｍｕｌｔｉｐｈａｓｅａｎｄｍｕｌｔｉｃｏｍｐｏｓｉｔｉｏｎ ,ｔｈｅｃｒｕｘｏｆ…     

A lattice Boltzmann model for simulation of compressible flows

This paper presents a new model of lattice Boltzmann method for full compressible flows. On the basis of multi‐speed model, an extra potential energy distribution function is introduced to recover the full compressible Navier–Stokes equations with a flexible specific‐heat ratio and Prandtl number. The Chapman–Enskog expansion of the kinetic equations is performed, and the two‐dimension‐seventeen‐velocity density equilibrium distribution functions are obtained. The governing equations are discretized using the third order monotone upwind scheme for scalar conservation laws finite volume scheme. The van Albada limiter is used to avoid spurious oscillations. In order to verify the accuracy of this double‐distribution‐function model, the Riemann problems, Couette flows, and flows around a NACA0012 airfoil are simulated. It is found that the proposed lattice Boltzmann model is suitable for compressible flows, even for strong shock wave problem, which has an extremely large pressure ratio, 100,000. Copyright © 2014 John Wiley & Sons, Ltd.     

Numerical study of Prandtl effect on MHD flow at a lid‐driven porous cavity

In this paper, the lattice Boltzmann method is used to study the Prandtl number effect on flow structure and heat transfer rates in a magnetohydrodynamic flow mixed convection in a lid‐driven cavity filled with a porous medium. The right and left walls are at constant but different temperatures (θ_h and θ_c), while the other walls are adiabatic. Gallium and salt water (0.02 < Pr < 13.4) are used as samples of the electroconducting fluids in the cavity. Typical sets of streamlines and isotherms are presented to analyze the flow patterns set up by the competition among the forced flow created by the lid‐driven wall, the buoyancy force of the fluid and the magnetic force of the applied magnetic field. Mathematical formulations in the porous media were constructed based on the Brinkman–Forchheimer model, while the multidistribution‐function model was used for the magnetic field effect. Numerical results were obtained and the effects of the Prandtl number and the other effective parameters such as Richardson, Hartman, and Darcy numbers were investigated. It was found that the fluid fluctuations within the cavity were reduced by increasing the Hartman number. A similar pattern was observed for the Darcy number reduction. Heat transfer was essentially dominated by the conduction for the low Prandtl number and forced convection dominated as the Prandtl number increased. Also, the average Nusselt number was raised by increasing the Prandtl number. It was discovered that a remarkable heat transfer enhancement of up to 28% could be reached by increasing the Prandtl number (from 0.02 to 13.4) at constant Richardson and Darcy numbers. Copyright © 2011 John Wiley & Sons, Ltd.     

Lattice Boltzmann simulations for flow and heat/mass transfer problems in a three‐dimensional porous structure

The lattice Boltzmann method (LBM) for a binary miscible fluid mixture is applied to problems of transport phenomena in a three‐dimensional porous structure. Boundary conditions for the particle distribution function of a diffusing component are described in detail. Flow characteristics and concentration profiles of diffusing species at a pore scale in the structure are obtained at various Reynolds numbers. At high Reynolds numbers, the concentration profiles are highly affected by the flow convection and become completely different from those at low Reynolds numbers. The Sherwood numbers are calculated and compared in good agreement with available experimental data. The results indicate that the present method is useful for the investigation of transport phenomena in porous structures. Copyright © 2003 John Wiley & Sons, Ltd.     

Lattice Boltzmann simulations of turbulent shear flow between parallel porous walls

The effects of two parallel porous walls are investigated, consisting of the Darcy number and the porosity of a porous medium, on the behavior of turbulent shear flows as well as skin-friction drag. The turbulent channel flow with a porous surface is directly simulated by the lattice Boltzmann method （LBM）. The Darcy-Brinkman- Forcheimer （DBF） acting force term is added in the lattice Boltzmann equation to simu- late the turbulent flow bounded by porous walls. It is found that there are two opposite trends （enhancement or reduction） for the porous medium to modify the intensities of the velocity fluctuations and the Reynolds stresses in the near wall region. The parametric study shows that flow modification depends on the Darcy number and the porosity of the porous medium. The results show that, with respect to the conventional impermeable wall, the degree of turbulence modification does not depend on any simple set of param- eters obviously. Moreover, the drag in porous wall-bounded turbulent flow decreases if the Darcy number is smaller than the order of O（10-4） and the porosity of porous walls is up to 0.4.     

Numerical simulation of laminar jet-forced flow using lattice Boltzmann method

In the paper, a numerical study on symmetrical and asymmetrical laminar jet-forced flows is carried out by using a lattice Boltzmann method （LBM） with a special boundary treatment. The simulation results are in very good agreement with the available numerical prediction. It is shown that the LBM is a competitive method for the laminar jet-forced flow in terms of computational efficiency and stability.     

A new stochastic method of reconstructing porous media

We present a new stochastic method of reconstructing porous medium from limited morphological information obtained from two-dimensional micro- images of real porous medium. The method is similar to simulated annealing method in the capability of reconstructing both isotropic and anisotropic structures of multi-phase but differs from the latter in that voxels for exchange are not selected completely randomly as their neighborhood will also be checked and this new method is much simpler to implement and program. We applied it to reconstruct real sandstone utilizing morphological information contained in porosity, two-point probability function and linear-path function. Good agreement of those references verifies our developed method’s powerful capability. The existing isolated regions of both pore phase and matrix phase do quite minor harm to their good connectivity. The lattice Boltzmann method (LBM) is used to compute the permeability of the reconstructed system and the results show its good isotropy and conductivity. However, due to the disadvantage of this method that the connectivity of the reconstructed system’s pore space will decrease when porosity becomes small, we suggest the porosity of the system to be reconstructed be no less than 0.2 to ensure its connectivity and conductivity.     

Local discontinuous galerkin method for radial porous flow with dispersion and adsorption

Based on the local discontinuous Galerkin methods for time-dependent convection-diffusion systems newly developed by Corkburn and Shu, according to the form of the generalized convection-diffusion equations which model the radial porous flow with dispersion and adsorption, a local discontinuous Galerkin method for radial porous flow with dispersion and adsorption was developed, a high order accurary new scheme for radial porous flow is obtained. The presented method was applied to the numerical tests of two cases of radial porous, i. e. , the convection-dispersion flow and the convection-dispersion-adsorption flow, the corresponding parts of the numerical results are in good agreement with the published solutions, so the presented method is reliable. Reckoning of the computational cost also shows that the method is practicable.     

Intercomparison of boundary schemes in Lattice Boltzmann method for flow simulation in porous media

The Lattice Boltzmann method has been widely adopted to simulate flow in porous media. The choice of appropriate boundary schemes is essential to achieve simulation accuracy; however, the criteria for the most suitable boundary treatment in the simulation of flow in porous media flow remain unresolved. Here, three types of the most commonly used boundary conditions are tested: interpolation bounce back (IBB), partial saturated method (PSM), and immersed boundary method (IBM). The dimensionless drag of face-centered cubic (FCC) sphere array and the dimensionless permeability of a random closely packed (RCP) sphere array are calculated and compared at different viscosities and resolutions. In the FCC sphere array case where spheres are not contacted, the IBB and PSM exhibit the same accuracy and both are of the second-order convergence rate. The IBM is less accurate and is of the first-order convergence rate. In the RCP sphere array case where the spheres are contacted, the IBB shows finer results and a second-order convergence rate. PSM underestimates the dimensionless permeability and increases resolution only slightly improved the results of PSM. The IBM overestimates the dimensionless permeability. These results indicate that among the three methods, the IBB is the most accurate. The PSM has the same accuracy as the IBB when sediments are not contacted; however, it loses its accuracy in the simulation of flow in closely packed porous media. This work could serve as a benchmark for further research in choosing the most appropriate method in the simulation of flow in porous media.     

A domain decomposition method for modelling Stokes flow in porous materials

An algorithm is presented for solving the Stokes equation in large disordered two‐dimensional porous domains. In this work, it is applied to random packings of discs, but the geometry can be essentially arbitrary. The approach includes the subdivision of the domain and a subsequent application of boundary integral equations to the subdomains. This gives a block diagonal matrix with sparse off‐block components that arise from shared variables on internal subdomain boundaries. The global problem is solved using a biconjugate gradient routine with preconditioning. Results show that the effectiveness of the preconditioner is strongly affected by the subdomain structure, from which a methodology is proposed for the domain decomposition step. A minimum is observed in the solution time versus subdomain size, which is governed by the time required for preconditioning, the time for vector multiplications in the biconjugate gradient routine, the iterative convergence rate and issues related to memory allocation. The method is demonstrated on various domains including a random 1000‐particle domain. The solution can be used for efficient recovery of point velocities, which is discussed in the context of stochastic modelling of solute transport. Copyright © 2002 John Wiley & Sons, Ltd.     

Theoretical and numerical study of Couette-Taylor flow with an axial flow using lattice Boltzmann method

In this study, the numerical models for swirling flows developed by Li et al and Zhou for lattice Boltzmann method (LBM) are chosen. These models were firstly validated using the Couette-Taylor flow between two concentric cylinders simulations. Numerical results showed the efficiency of the Zhou's model. Numerical simulation results using LBM are in good agreement for the steady and unsteady regimes compared to the literature review. In a second step, the Zhou model was then adopted to our study to determine the Couette-Taylor instabilities with an axial flow. Two protocols are tested. The first one (direct protocol) starts with an azimuthal flow without any axial flow (Re = 0). Once the regime is established, an axial flow is then superposed to the Couette-Taylor flow (with a sudden or a progressive manner). The second protocol (inverse protocol) starts with an axial flow at a given Reynolds number (Poiseuille flow). Once the regime is established, an azimuthal flow is the executed (with a sudden or a progressive manner). The effect of various parameters controlling the physical situation is also discussed. The increase of the azimuthal velocity mainly led to the emergence and development of Taylor vortices. Its influence decreases when the axial Reynolds number increases. The relevant result for this study is the change of the critical axial Reynolds number Re_c (total disappearance of instabilities) with both protocols and both manners.