Shan-Chen Model Research Articles

Dielectrophoresis effect on anti-icing characteristics of superhydrophobic surface

In present study, icing experiment using semiconductor refrigeration method and 3D LBM (Lattice Boltzmann method) simulation by employing improved Shan-Chen model coupling the electric field model are conducted to confirm dielectrophoresis effect on freezing characterisitic of a sessile droplet at superhydrophobic surface. Tests are performed with the same droplet diameter of d = 2.6 mm, cold wall temperature of TW = -5 °C and ambient condition (TA = 20 °C, RH = 48%). The sessile droplet freezing process is recorded at different contact angles and electric field intensities. The improved Shan-Chen model with large density ratio (ρl/ρv = 36.55), electric field model and new wettability model are applied to simulate droplet deformation. The present results show that the droplet icing time has a positive relationship with contact angle and electric field intensity, the droplet icing time is extended up to about 40 times by the couple effect of superhydrophobic coating and the electric field. The size of droplet-solid interface and electric-induced convection are the key factors of anti-icing performance.

Simulating Two-Phase Seepage in Undisturbed Soil Based on Lattice Boltzmann Method and X-ray Computed Tomography Images.

The two-phase seepage fluid (i.e., air and water) behaviors in undisturbed granite residual soil (U-GRS) have not been comprehensively studied due to a lack of accurate and representative models of its internal pore structure. By leveraging X-ray computed tomography (CT) along with the lattice Boltzmann method (LBM) enhanced by the Shan-Chen model, this study simulates the impact of internal pore characteristics of U-GRS on the water-gas two-phase seepage flow behaviors. Our findings reveal that the fluid demonstrates a preference for larger and straighter channels for seepage, and as seepage progresses, the volume fraction of the water/gas phases exhibits an initial increase/decrease trend, eventually stabilizing. The results show the dependence of two-phase seepage velocity on porosity, while the local seepage velocity is influenced by the distribution and complexity of the pore structure. This emphasizes the need to consider pore distribution and connectivity when studying two-phase flow in undisturbed soil. It is observed that the residual gas phase persists within the pore space, primarily localized at the pore margins and dead spaces. Furthermore, the study identifies that hydrophobic walls repel adjacent fluids, thereby accelerating fluid movement, whereas hydrophilic walls attract fluids, inducing a viscous effect that decelerates fluid flow. Consequently, the two-phase flow rate is found to increase with then-enhanced hydrophobicity. The apex of the water-phase volume fraction is observed under hydrophobic wall conditions, reaching up to 96.40%, with the residual gas-phase constituting 3.60%. The hydrophilic wall retains more residual gas-phase volume fraction than the neutral wall, followed by the hydrophobic wall. Conclusively, the investigations using X-ray CT and LBM demonstrate that the pore structure characteristics and the wettability of the pore walls significantly influence the two-phase seepage process.

Lattice Boltzmann simulation of droplets accelerate shedding on an inverted triangular surface with wettability gradient

The droplet combining and shedding on a gradient-wetting inverted triangle surface is investigated using the lattice Boltzmann method based on the Shan-Chen pseudo-potential model. The impacts of the tip angle, wetting gradient, and Bond number are investigated. The study shows that the smaller the tip angle, the better for the droplets to merge and shed. The droplet combining and shedding process can be sped up by the triangular surface’s wetting gradient. As the wetting gradient increases, the droplet combining and shedding becomes easier. Applying a wetting gradient to speed drainage is better when the original surface is more hydrophilic. The combination and shedding of droplets on the wet surface are significantly influenced by the Bond number as well.

Metastable and unstable hydrodynamics in multiphase lattice Boltzmann.

Metastability in liquids is at the foundation of complex phase transformation dynamics such as nucleation and cavitation. Intermolecular interaction details, beyond the equationof state, and thermal hydrodynamic fluctuations play a crucial role. However, most numerical approaches suffer from a slow time and space convergence, thus hindering the convergence to the hydrodynamic limit. This work shows that the Shan-Chen lattice Boltzmann model has the unique capability of simulating the hydrodynamics of the metastable state. The structure factor of density fluctuations is theoretically obtained and numerically verified to a high precision, for all simulated wave vectors, reduced temperatures, and pressures, deep into the metastable region. Such remarkable agreement between the theory and simulations leverages the exact implementation at the lattice level of the mechanical equilibrium condition. The static structure factor is found to consistently diverge as the temperature approaches the critical point or the density approaches the spinodal line at a subcritical temperature. Theoretically predicted critical exponents are observed in both cases. Finally, the phase separation in the unstable branch follows the same pattern, i.e., the generation of interfaces with different topology, as observed in molecular dynamics simulations.

Pore-Scale Modeling of Pressure-Driven Flow and Spontaneous Imbibition in Fracturing-Shut-In-Flowback Process of Tight Oil Reservoirs

Tight oil reservoirs are characterized by multiple pore spaces where nano-micropores and multiscale fractures coexist, and each type of medium varies in scale, implying a tight coupling of multiscale fractures with matrix and giving rise to extremely complicated flow patterns. To further investigate its flow mechanism, we first construct three two-dimensional (2-D) fracture-pore geometry models based on microfocus computed tomography (CT) imaging of a typical tight rock. A pore-scale modeling workflow is thereafter developed using the Shan-Chen lattice Boltzmann model (SC-LBM) to simulate the pressure-driven flow and spontaneous imbibition. The influence of fracture-pore geometry on the pore-scale fluid exchange dynamics in the fracturing-shut-in-flowback process is clearly clarified. Results show that, for the porous medium model without fracture, the fracturing fluid can displace and replace some crude oil by spontaneous imbibition while a large amount of crude oil droplets remains unexploited away from the oil/water contact line, resulting in a low oil imbibition recovery. The injected fracturing fluid migrates into the deep position along the fractures, only a few entering the matrix pore space near the fractures. Small pores are the main channel for fracturing fluid to imbibe into the matrix pores, and the replaced crude oil droplets flow into fractures through large pores as intermittent or continuous pipe flow. The complexity of fracture network typically exerts a significant impact on the fluid exchange dynamics during pressure-driven flow and spontaneous imbibition. The more complex the fracture network, the larger the volume of fracturing fluid injected, the easier for oil droplets replaced from the matrix pores, and the more difficult for the flowback of fracturing fluid in the fracture-pore geometry model. Our understanding will provide a basis for explaining the underlying mechanisms of oil replacement by pressure-driven flow and spontaneous imbibition in the fracturing-shut-in-flowback process.

Simulation of phase change during the freezing of unsaturated porous media by using a coupled lattice Boltzmann model

PurposeThe purpose of this paper is to develop a coupled lattice Boltzmann model for the simulation of the freezing process in unsaturated porous media.Design/methodology/approachIn the developed model, the porous structure with complexity and disorder was generated by using a stochastic growth method, and then the Shan-Chen multiphase model and enthalpy-based phase change model were coupled by introducing a freezing interface force to describe the variation of phase interface. The pore size of porous media in freezing process was considered as an influential factor to phase transition temperature, and the variation of the interfacial force formed with phase change on the interface was described.FindingsThe larger porosity (0.2 and 0.8) will enlarge the unfrozen area from 42 mm to 70 mm, and the rest space of porous medium was occupied by the solid particles. The larger specific surface area (0.168 and 0.315) has a more fluctuated volume fraction distribution.Originality/valueThe concept of interfacial force was first introduced in the solid–liquid phase transition to describe the freezing process of frozen soil, enabling the formulation of a distribution equation based on enthalpy to depict the changes in the water film. The increased interfacial force serves to diminish ice formation and effectively absorb air during the freezing process. A greater surface area enhances the ability to counteract liquid migration.

Lattice Boltzmann simulations of quasi-steady film and axisymmetric nucleate boiling

An axisymmetric multiple relaxation time lattice Boltzmann method utilizing the Shan-Chen pseudo-potential model is developed and combined with an axisymmetric finite difference approximation of the energy equation to form a hybrid method with a view toward studying axisymmetric nucleate boiling. The mechanism of phase change in the Shan–Chen model is investigated, and the model is validated by simulating a Stefan problem as well as simulations of quasi-steady film boiling with comparisons to established results. Axisymmetric quasi-steady nucleate boiling is then investigated including examining the effect of the wetting properties of surfaces by varying the wettability force to vary the dynamic contact angle.

Combining neutral scalar and isothermal Shan-Chen lattice Boltzmann method to simulate droplet placement on a wall in isothermal and non-isothermal states

BackgroundIn the current article, two-phase thermal fluxes are created by combining the thermal model of the neutral scalar model with the two-phase Shan-Chen model of the lattice Boltzmann method (LBM). MethodsThe different intermolecular powers for the isothermal Shan-Chen model show how a droplet would be placed on a wall. By raising the droplet intermolecular power parameter, the surface area increases and becomes wet. Next, the isothermal Shan-Chen method and the neutral scalar method are combined to investigate multi-phase thermal problems. The droplet placement on the hot wall is therefore done at relatively high Rayleigh numbers. By raising the Rayleigh number, the isothermal lines within the droplet's interior gradually become less ascending and less descending until they eventually achieve a uniform state when it is placed against a hot wall. Additionally, the channel's Rayleigh-Bénard convective heat transfer is enhanced by increasing the Rayleigh number. Significant findingsNatural convection in the enclosures can be used in solar collectors. As the Rayleigh number increases, the average Nusselt number (Nuavg) rises as would be expected. The results demonstrate that LBM is a practical method for simulating multi-phase thermal flows.

Comparative analysis and computational optimization of potential-based multiphase lattice Boltzmann models

The potential-based multiphase lattice Boltzmann models are widely used because they root in thermodynamics and evade the interface tracking or integrating. This paper investigates several potential-based models with the common equations of state (EOS) by the theoretical analyses and numerical computations of the thermodynamic consistency and spurious currents. Surprisingly, the Shan–Chen model presents a superior accuracy compared to the Zhang–Chen models, although they are mathematically equivalent. We find that the great improvement is attributed to the square root form of the pseudopotential model, which significantly lessens the error of numerical gradient calculation. Inspired by the improvement, a general formula φ′=n−1φ1−n∂x(φn) is introduced for calculating the gradient, and the coefficient n=0.1 yields better results than n=0.5, which is equal to the pseudopotential model. This scheme is further applied to optimize the evaluation of the chemical potential model. The improved chemical potential model displays lower numerical errors in the liquid–gas transition region and smaller spurious currents near the curved phase interface than the pseudopotential model. Additionally, the improved model is confirmed to meet the Young–Laplace law and Galilean invariance.

Lattice Boltzmann simulation of droplet coalescence in fuel cell flow channel with mechanical vibration

Fuel cells are the promising power source for future green vehicles and water management is an important issue for proton exchange membrane (PEM) fuel cells, which affects the lifetime and durability of fuel cells critically. In the available studies, the key parameters such as channel shape, airflow rate and material wettability have been extensively studied for water management. However, so far there is no open research on the effect of mechanical vibration on water management although most fuel cells work under vibration conditions. To bridge the knowledge gap, based on the Shan-Chen lattice Boltzmann model, this paper aims to reveal the water droplet coalescence process in an airflow channel under vibration conditions and several interesting observations are obtained for the first time. A large vibration amplitude or a small vibration frequency will not only hinder the movement of droplets, but also accelerate the coalescence of droplets. Also by comparing the simulation results it can be found that the droplet movement velocity decreases significantly after coalescence. Therefore, it is important to avoid large amplitude and small frequency for the operation of PEM fuel cells.

Pore-scale modelling of water sorption in nanopore systems of shale

Water vapor sorption in nanoporous media with complex pore structures, like shale, is not well understood. To address this, a pseudopotential lattice Boltzmann method (LBM) is developed to investigate water sorption behavior. The LBM model incorporates long-range molecular forces using a modified Shan-Chen model based on the Carnahan-Starling equation of state. The simulation results show that the presence of water films in nanopores can create a liquid pressure difference of up to 100 MPa between confined and free states. The adsorption theories based on simple pore shapes are not applicable to nanoporous systems with complex pore geometries. Additionally, when the relative humidity exceeds 1, water vapor condenses inside hydrophobic nanostructures while being attracted by neighboring liquid, leading to cluster formation. In water-wet nanoporous media, capillary condensation occurs progressively from small throats to large pores, and the sorption curves vary smoothly with pressure variation. In mixed-wet nanoporous media, the sorption curves exhibit no hysteresis in the early desorption stage, where the adsorption and desorption processes occur in the region around the hydrophobic particles and are reversible. Overall, this study provides valuable insights into water sorption behavior in nanopore systems of shale and sets the foundation for modelling liquid-vapor distribution in nanoporous media at the pore scale.

Droplet-driven particle motions in fluids studied using coupled lattice Boltzmann method

A numerical method for describing droplet-driven particle motion at low Weber numbers in fluids is proposed. A coupling strategy combining the Shan–Chen (SC) model and the smoothed profile method (SPM) is employed. The proposed scheme correctly resolves the momentum transfer between solid particles and fluid phases, while effectively controlling the wetting condition. An interaction zone where three phases coexist is introduced to solve the contradiction between the SC model and SPM (i.e. whether fluid particles are allowed at the solid nodes). In the interaction zone, partial solid particles are entrapped in the fluids owing to the coexistence of solid and fluid particles, causing solid particles to exchange momentum with fluid particles. Furthermore, the interaction forces between fluid and solid particles near the solid surface are considered at nodes near pure solid particles within this region. Based on the proposed scheme, the interactions between freely moving particles and freely moving droplets were investigated, along with the effects of wettability (Gads) and volume forces (Fp*) on particle–droplet interactions. When a liquid droplet comes into contact with a solid particle, the adhesion force between the droplet and the solid surface promotes the movement of solid particles toward the droplet. When a double emulsion comes into contact with solid particles, the adhesion force between the emulsion film and solid surface causes the double emulsion to break at the connecting point, and the volume force of the emulsion film Fd* causes the solid particles to penetrate the emulsion.

Numerical Study of Microstructures and Roughness Design Effects on Surface Hydrophilicity through the Lattice Boltzmann Method

Hydrophilicity is one of the most vital characteristics of titanium (Ti) implants. Surface structure design is a powerful and efficient strategy for improving the intrinsic hydrophilic ability of Ti implants. Existing research has focused on experimental exploration, and hence, a reliable numerical model is needed for surface structure design and corresponding hydrophilicity prediction. To address this challenge, we proposed a numerical model to analyze the droplet dynamics on Ti surfaces with specific microstructures designed through the lattice Boltzmann method (LBM). In this work, a Shan-Chen (SC) model was applied in the simulations. We simulated droplets spreading on smooth and micropillar surfaces with various wettability and provided a comprehensive discussion of the edge locations, contact line, droplet height, contact area, surface free energy, and forces to reveal more details and mechanisms. To better tune and control the surface hydrophilicity, we investigated the effects of micropillar geometric sizes (pillar width a, height h, and pitch b) on hydrophilicity via single factor analysis and the response surface method (RSM). The results show that the hydrophilicity initially increases and then decreases with an increasing a, increases with an increasing h, and decreases with an increasing d. In addition, the interaction effects of a-d and h-d are significant. The optimization validation of the RSM also demonstrates the accuracy of our lattice Boltzmann (LB) model with an error of 0.687%. Here, we defined a dimensionless parameter ξ to integrate the geometric parameters and denote the surface roughness. The hydrophilicity of Ti surfaces improves with an increasing surface roughness. In addition, the effect of the microstructure geometry shape was investigated under the same value of surface roughness. Surfaces with micropillars show the best hydrophilicity. Moreover, this study is expected to provide an accurate and reliable LB model for predicting and enhancing the intrinsic hydrophilicity of Ti surfaces via specific microstructure and roughness designs.

Pore-scale simulation of gas displacement after water flooding using three-phase lattice Boltzmann method

Water flooding is a commonly used technique to improve oil recovery, although the amount of oil left in reservoirs after the procedure is still significant. Gas displacement after water flooding is an effective way to recover residual oil, but the occurrence state and flow principles of multiphase fluid after gas injection are still ambiguous. Therefore, the gas displacement process after water flooding should be studied on the pore scale to provide a basis for formulating a reasonable gas injection program. Most of the current pore-scale studies focus on two-phase flow, while simulations that account for the influence of oil-gas miscibility and injected water are seldom reported. In this work, the multi-component multi-phase Shan-Chen lattice Boltzmann model is used to simulate the gas displacement after water flooding in a porous medium, and the effects of injected water, viscosity ratio, pore structure, and miscibility are analyzed. It is established that the injected water will cause gas flow path variations and lead to premature gas channeling. Under the impact of capillary pressure, the water retained in the porous medium during the water flooding stage further imbibes into the tiny pores during gas injection and displaces the remaining oil. When miscibility is considered, the oil-gas interface disappears, eliminating the influence of the capillary effect on the fluid flow and enabling the recovery of remaining oil at the corner. This study sheds light on the gas displacement mechanisms after water flooding from the pore-scale perspective and provides a potential avenue for improving oil recovery. Cited as: Wang, S., Chen, L., Feng, Q., Chen, L., Fang, C., Cui, R. Pore-scale simulation of gas displacement after water flooding using three-phase lattice Boltzmann method. Capillarity, 2023, 6(2): 19-30. https://doi.org/10.46690/capi.2023.02.01

PORE SCALE INVESTIGATION ON THE EFFECT OF TEMPERATURE ON THREE-PHASE REACTIVE FLOW IN POROUS MEDIA USING THE LATTICE BOLTZMANN METHOD

In this study, the process of carbonate rock acidification with three-phase immiscible fluid [water, oil, and supercritical carbon dioxide (scCO<sub>2</sub>)] reactive flow in nonisothermal condition at pore scale was simulated by the lattice Boltzmann method (LBM). The three immiscible fluids' flow, solute transport, and heat transfer were solved by the Shan-Chen multiphase multicomponent lattice Boltzmann model, mass transport lattice Boltzmann model, and multicomponent thermal lattice Boltzmann model, respectively. The solid phase was updated by the volume-of-pixel method. In addition, the effect of temperature on the multiphase multicomponent reactive flow was numerically investigated. The results show that the dissolution rate is determined by three-phase distribution. The increase of temperature can accelerate the dissolution rate and decrease the inhibition influence of the non-acid phase on dissolution. Furthermore, the influence of temperature decreases with the increase of inlet velocity and oil wettability.

Pre-wetting of sand for high speed oil-water separation

The search for a simple, affordable, and rapid separation method for oil-water mixtures is one of the most serious challenges that the world is facing today. Superhydrophobic, or superoleophilic, materials and surfaces are promising choices for these types of separations. In the present study, a high speed gravity-driven oil–water separation method has been realized by using pre-wetting of river sand, sea sand, and desert sand. The superoleophobicity properties of these pre-wetted sand materials have been found to be good, with a maximum oil contact angle of about 141.1° for desert sand and a minimum oil contact angel of 113.8° for sea sand. Also, the relationship between the porosity of the respective sand type and the water phase separation rate was, for the first time, obtained by using the LBM Shan-Chen model and micro-CT characterization. As a result of a simple particle size screening, the highest kerosene-water separation rate for sea sand was 24.3 ± 2.4 mm/s. In fact, this is the highest oil-water separation rate that has ever been achieved for a sand material. In addition to the effect of sand particle size, other factors affecting the rate of oil-water separation (such as number of separation cycles, thickness of sand layer, temperature, pH value, salt concentration, and different oil-water mixtures) have also been investigated.

Lattice Boltzmann simulation of wetting gradient accelerating droplets merging and shedding on a circumferential surface

A lattice Boltzmann method (LBM) based on Shan-Chen pseudo-potential model is used to investigate the process of droplets merging and shedding on a gradient wetting circular surface. The effects of wetting gradient, radius, and radius ratio on droplet merging and shedding were mainly explored. The results show that applying a wetting gradient on the circumferential surface can accelerate the process of droplet merging and shedding. The velocity of droplets merging and shedding increase with the increase of the wetting gradient. The more hydrophilic the original surface, the better the optimization effect of accelerating drainage is by applying a wetting gradient. The radius and radius ratio significantly affect droplets merging and shedding on the gradient wetting surface. The larger the radius and radius ratio, the shorter the droplets merging and shedding time.

A Quantitative Comparison of Physical Accuracy and Numerical Stability of Lattice Boltzmann Color Gradient and Pseudopotential Multicomponent Models for Microfluidic Applications

The performances of the Color-Gradient (CG) and of the Shan-Chen (SC) multicomponent Lattice Boltzmann models are quantitatively compared side-by-side on multiple physical flow problems where breakup, coalescence and contraction of fluid ligaments are important. The flow problems are relevant to microfluidic applications, jetting of microdroplets as seen in inkjet printing, as well as emulsion dynamics. A significantly wider range of parameters is shown to be accessible for CG in terms of density-ratio, viscosity-ratio and surface tension values. Numerical stability for a high density ratio $\mathcal{O}(1000)$ is required for simulating the drop formation process during inkjet printing which we show here to be achievable using the CG model but not using the SC model. In terms of physical accuracy, the CG model shows good agreement with analytical solutions for droplet oscillation and ligament contraction test-cases. The SC model is effective in accurately simulating ligament contraction, but less accurate in handling droplet oscillations. Rayleigh-Plateau instability simulations show that both the CG and SC models give physically realistic results when breakup occurs, although smaller satellite droplets quickly vanish in the SC model due to droplet mass evaporating into the ambient phase. Our results show that the CG model is a suitable choice for challenging simulations of droplet formation, due to a combination of both numerical stability and physical accuracy. We also present a novel approach to incorporate repulsion forces between interfaces for CG, with possible applications to the study of stabilized emulsions. Specifically, we show that the CG model can produce similar results to a known multirange potentials extension of the SC model for modelling a disjoining pressure, opening up its use for the study of dense stabilized emulsions.

Study of Fingering Dynamics of Two Immiscible Fluids in a Homogeneous Porous Medium with Considering Wettability Effects Using a Pore-Scale Multicomponent Lattice Boltzmann Model

In the present study, a pore-scale multicomponent lattice Boltzmann method (LBM) is employed for the investigation of the immiscible-phase fluid displacement in a homogeneous porous medium. The viscous fingering and the stable displacement regimes of the invading fluid in the medium are quantified which is beneficial for predicting flow patterns in pore-scale structures, where an experimental study is extremely difficult. Herein, the Shan-Chen (S-C) model is incorporated with an appropriate collision model for computing the interparticle interaction between the immiscible fluids and the interfacial dynamics. Firstly, the computational technique is validated by a comparison of the present results obtained for different benchmark flow problems with those reported in the literature. Then, the penetration of an invading fluid into the porous medium is studied at different flow conditions. The effect of the capillary number (Ca), dynamic viscosity ratio (M), and the surface wettability defined by the contact angle (θ) are investigated on the flow regimes and characteristics. The obtained results show that for M<1, the viscous fingering regime appears by driving the invading fluid through the pore structures due to the viscous force and capillary force. However, by increasing the dynamic viscosity ratio and the capillary number, the invading fluid penetrates even in smaller pores and the stable displacement regime occurs. By the increment of the capillary number, the pressure difference between the two sides of the porous medium increases, so that the pressure drop Δp along with the domain at θ=40∘ is more than that of computed for θ=80∘. The present study shows that the value of wetting fluid saturation Sw at θ=40∘ is larger than its value computed with θ=80∘ that is due to the more tendency of the hydrophilic medium to absorb the wetting fluid at θ=40∘. Also, it is found that the magnitude of Sw computed for both the contact angles is decreased by the increment of the viscosity ratio from Log(M)=−1 to 1. The present study demonstrates that the S-C LBM is an efficient and accurate computational method to quantitatively estimate the flow characteristics and interfacial dynamics through the porous medium.

A Pseudopotential Lattice Boltzmann Method for Simulation of Two-Phase Flow Transport in Porous Medium at High-Density and High–Viscosity Ratios

In this work, the capability of a multiphase lattice Boltzmann method (LBM) based on the pseudopotential Shan-Chen (S-C) model is investigated for simulation of two-phase flows through porous media at high-density and high–viscosity ratios. The accuracy and robustness of the S-C LBM are examined by the implementation of the single relaxation time (SRT) and multiple relaxation time (MRT) collision operators with integrating the forcing schemes of the shifted velocity method (SVM) and the exact difference method (EDM). Herein, two equations of state (EoS), namely, the standard Shan-Chen (SC) EoS and Carnahan-Starling (CS) EoS, are implemented to assay the performance of the numerical technique employed for simulation of two-phase flows at high-density ratios. An appropriate modification in the forcing schemes is also used to remove the thermodynamic inconsistency in the simulation of two-phase flow problems studied at low reduced temperatures. The comparative study of these improvements of the S-C LBM is performed by considering an equilibrium state of a droplet suspended in the vapor phase. The solver is validated against the analytical coexistence curve for the liquid-vapor system and the surface tension estimation from the Laplace Law. Then, according to the results obtained, a conclusion has been made to choose an efficient numerical algorithm, including an appropriate collision operator, a realistic EoS, and an accurate forcing scheme, for simulation of multiphase flow transport through a porous medium. The patterns of two-phase flow transport through the porous medium are predicted using the present numerical scheme in different flow conditions defined by the capillary number and the dynamic viscosity ratio. The results obtained for the nonwetting phase saturation, penetration structure of the invading fluid, and the displacement patterns of two-phase flow in the porous medium are comparable with those reported in the literature. The present study demonstrates that the S-C LBM with employing the MRT-EDM scheme, CS EoS, and the modified forcing scheme is efficient and accurate for estimation of the two-phase flow characteristics through the porous medium.