1
|
Dou S, Hao L, Liu H. A mesoscopic model for simulating the physisorption process in nanopores. Chem Eng Sci 2022. [DOI: 10.1016/j.ces.2022.117988] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
2
|
Li Q, Yu Y, Huang RZ. Achieving thermodynamic consistency in a class of free-energy multiphase lattice Boltzmann models. Phys Rev E 2021; 103:013304. [PMID: 33601620 DOI: 10.1103/physreve.103.013304] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2020] [Accepted: 12/18/2020] [Indexed: 06/12/2023]
Abstract
The free-energy lattice Boltzmann (LB) model is one of the major multiphase models in the LB community. The present study is focused on a class of free-energy LB models in which the divergence of thermodynamic pressure tensor or its equivalent form expressed by the chemical potential is incorporated into the LB equation via a forcing term. Although this class of free-energy LB models may be thermodynamically consistent at the continuum level, it suffers from thermodynamic inconsistency at the discrete lattice level owing to numerical errors [Guo et al., Phys. Rev. E 83, 036707 (2010)10.1103/PhysRevE.83.036707]. The numerical error term mainly includes two parts: one comes from the discrete gradient operator and the other can be identified in a high-order Chapman-Enskog analysis. In this paper, we propose an improved scheme to eliminate the thermodynamic inconsistency of the aforementioned class of free-energy LB models. The improved scheme is constructed by modifying the equation of state of the standard LB equation, through which the discretization of ∇(ρc_{s}^{2}) is no longer involved in the force calculation and then the numerical errors can be significantly reduced. Numerical simulations are subsequently performed to validate the proposed scheme. The numerical results show that the improved scheme is capable of eliminating the thermodynamic inconsistency and can significantly reduce the spurious currents in comparison with the standard forcing-based free-energy LB model.
Collapse
Affiliation(s)
- Q Li
- School of Energy Science and Engineering, Central South University, Changsha 410083, China
| | - Y Yu
- School of Energy Science and Engineering, Central South University, Changsha 410083, China
| | - R Z Huang
- School of Energy Science and Engineering, Central South University, Changsha 410083, China
| |
Collapse
|
3
|
Huang J, Yin X, Killough J. Thermodynamic consistency of a pseudopotential lattice Boltzmann fluid with interface curvature. Phys Rev E 2019; 100:053304. [PMID: 31869878 DOI: 10.1103/physreve.100.053304] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2019] [Indexed: 11/07/2022]
Abstract
Thermodynamic consistency of pseudopotential lattice Boltzmann models is a major topic that needs comprehensive evaluations. When interface is flat, pseudopotential models can give density-pressure isotherms in excellent agreement with those from equation of state. When interface is curved, thermodynamic equilibriums are affected by interface curvature, and consistency of pseudopotential models has not been systematically evaluated. In this study, we show that the effect of Laplace pressure on phase equilibrium is quantitatively consistent with Kelvin equation at high reduced temperatures (≥0.7). At low temperatures, inconsistency that can be attributed to the effect of orientation of the interface was noted, and it can be improved by tuning of the pseudopotential. By relating interfacial tension of a simulated fluid to that of a real fluid, the lattice spacing of pseudopotential model is found to be on the order of several molecular diameters, the typical range of intermolecular interactions. Interfacial thickness at different temperatures in pseudopotential model compared well with experiments and molecular dynamics simulations, which confirms that the calculated length scale is reasonable. Evaluation of a free energy lattice Boltzmann model indicate that it is consistent with Kelvin equation at high temperatures. The free energy model, however, is not as accurate as the tested pseudopotential model, and discrepancies may come from the relative inaccuracies in the predictions of vapor densities and the thinner interfaces.
Collapse
Affiliation(s)
- Jingwei Huang
- Department of Petroleum Engineering, Texas A&M University, College Station, Texas 77840, USA
| | - Xiaolong Yin
- Department of Petroleum Engineering, Colorado School of Mines, Golden, Colorado 80401, USA
| | - John Killough
- Department of Petroleum Engineering, Texas A&M University, College Station, Texas 77840, USA
| |
Collapse
|
4
|
Zarghami A, Van den Akker HEA. Thermohydrodynamics of an evaporating droplet studied using a multiphase lattice Boltzmann method. Phys Rev E 2017; 95:043310. [PMID: 28505732 DOI: 10.1103/physreve.95.043310] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2016] [Indexed: 11/07/2022]
Abstract
In this paper, the thermohydrodynamics of an evaporating droplet is investigated by using a single-component pseudopotential lattice Boltzmann model. The phase change is applied to the model by adding source terms to the thermal lattice Boltzmann equation in such a way that the macroscopic energy equation of multiphase flows is recovered. In order to gain an exhaustive understanding of the complex hydrodynamics during evaporation, a single droplet is selected as a case study. At first, some tests for a stationary (non-)evaporating droplet are carried out to validate the method. Then the model is used to study the thermohydrodynamics of a falling evaporating droplet. The results show that the model is capable of reproducing the flow dynamics and transport phenomena of a stationary evaporating droplet quite well. Of course, a moving droplet evaporates faster than a stationary one due to the convective transport. Our study shows that our single-component model for simulating a moving evaporating droplet is limited to low Reynolds numbers.
Collapse
Affiliation(s)
- Ahad Zarghami
- Department of Process and Energy, TU Delft, Netherlands
| | - Harry E A Van den Akker
- Transport Phenomena Group, Department of Chemical Engineering, TU Delft, Netherlands.,Bernal Institute, School of Engineering, University of Limerick, Ireland
| |
Collapse
|
5
|
Li Q, Zhou P, Yan HJ. Revised Chapman-Enskog analysis for a class of forcing schemes in the lattice Boltzmann method. Phys Rev E 2016; 94:043313. [PMID: 27841508 DOI: 10.1103/physreve.94.043313] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2016] [Indexed: 11/07/2022]
Abstract
In the lattice Boltzmann (LB) method, the forcing scheme, which is used to incorporate an external or internal force into the LB equation, plays an important role. It determines whether the force of the system is correctly implemented in an LB model and affects the numerical accuracy. In this paper we aim to clarify a critical issue about the Chapman-Enskog analysis for a class of forcing schemes in the LB method in which the velocity in the equilibrium density distribution function is given by u=∑_{α}e_{α}f_{α}/ρ, while the actual fluid velocity is defined as u[over ̂]=u+δ_{t}F/(2ρ). It is shown that the usual Chapman-Enskog analysis for this class of forcing schemes should be revised so as to derive the actual macroscopic equations recovered from these forcing schemes. Three forcing schemes belonging to the above class are analyzed, among which Wagner's forcing scheme [A. J. Wagner, Phys. Rev. E 74, 056703 (2006)10.1103/PhysRevE.74.056703] is shown to be capable of reproducing the correct macroscopic equations. The theoretical analyses are examined and demonstrated with two numerical tests, including the simulation of Womersley flow and the modeling of flat and circular interfaces by the pseudopotential multiphase LB model.
Collapse
Affiliation(s)
- Q Li
- School of Energy Science and Engineering, Central South University, Changsha 410083, China
| | - P Zhou
- School of Energy Science and Engineering, Central South University, Changsha 410083, China
| | - H J Yan
- School of Energy Science and Engineering, Central South University, Changsha 410083, China
| |
Collapse
|
6
|
Zarghami A, Looije N, Van den Akker H. Assessment of interaction potential in simulating nonisothermal multiphase systems by means of lattice Boltzmann modeling. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 92:023307. [PMID: 26382546 DOI: 10.1103/physreve.92.023307] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/18/2015] [Indexed: 06/05/2023]
Abstract
The pseudopotential lattice Boltzmann model (PP-LBM) is a very popular model for simulating multiphase systems. In this model, phase separation occurs via a short-range attraction between different phases when the interaction potential term is properly chosen. Therefore, the potential term is expected to play a significant role in the model and to affect the accuracy and the stability of the computations. The original PP-LBM suffers from some drawbacks such as being capable of dealing with low density ratios only, thermodynamic inconsistency, and spurious velocities. In this paper, we aim to analyze the PP-LBM with the view to simulate single-component (non-)isothermal multiphase systems at large density ratios and in spite of the presence of spurious velocities. For this purpose, the performance of two popular potential terms and of various implementation schemes for these potential terms is examined. Furthermore, the effects of different parameters (i.e., equation of state, viscosity, etc.) on the simulations are evaluated, and, finally, recommendations for a proper simulation of (non-)isothermal multiphase systems are presented.
Collapse
Affiliation(s)
- Ahad Zarghami
- Transport Phenomena Group, Department of Chemical Engineering, TU Delft, The Netherlands
| | - Niels Looije
- Transport Phenomena Group, Department of Chemical Engineering, TU Delft, The Netherlands
| | - Harry Van den Akker
- Transport Phenomena Group, Department of Chemical Engineering, TU Delft, The Netherlands
| |
Collapse
|
7
|
Kikkinides ES, Monson PA. Dynamic density functional theory with hydrodynamic interactions: Theoretical development and application in the study of phase separation in gas-liquid systems. J Chem Phys 2015; 142:094706. [DOI: 10.1063/1.4913636] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- E. S. Kikkinides
- Department of Mechanical Engineering, University of Western Macedonia, 50100 Kozani, Greece and Chemical Process and Energy Resources Institute (CPERI), Centre for Research and Technology Hellas (CERTH), 57001 Thermi-Thessaloniki, Greece
| | - P. A. Monson
- Department of Chemical Engineering, University of Massachusetts, 159 Goessmann Laboratory, 686 North Pleasant Street, Amherst, Massachusetts 01003-9303, USA
| |
Collapse
|
8
|
Lou Q, Guo Z. Interface-capturing lattice Boltzmann equation model for two-phase flows. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 91:013302. [PMID: 25679734 DOI: 10.1103/physreve.91.013302] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/05/2013] [Indexed: 06/04/2023]
Abstract
In this work, an interface-capturing lattice Boltzmann equation (LBE) model is proposed for two-phase flows. In the model, a Lax-Wendroff propagation scheme and a properly chosen equilibrium distribution function are employed. The Lax-Wendroff scheme is used to provide an adjustable Courant-Friedrichs-Lewy (CFL) number, and the equilibrium distribution is presented to remove the dependence of the relaxation time on the CFL number. As a result, the interface can be captured accurately by decreasing the CFL number. A theoretical expression is derived for the chemical potential gradient by solving the LBE directly for a two-phase system with a flat interface. The result shows that the gradient of the chemical potential is proportional to the square of the CFL number, which explains why the proposed model is able to capture the interface naturally with a small CFL number, and why large interface error exists in the standard LBE model. Numerical tests, including a one-dimensional flat interface problem, a two-dimensional circular droplet problem, and a three-dimensional spherical droplet problem, demonstrate that the proposed LBE model performs well and can capture a sharp interface with a suitable CFL number.
Collapse
Affiliation(s)
- Qin Lou
- State Key Laboratory of Coal Combustion, Huazhong University of Science and Technology, Wuhan 430074, China
| | - Zhaoli Guo
- State Key Laboratory of Coal Combustion, Huazhong University of Science and Technology, Wuhan 430074, China
| |
Collapse
|
9
|
Siebert DN, Philippi PC, Mattila KK. Consistent lattice Boltzmann equations for phase transitions. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2014; 90:053310. [PMID: 25493907 DOI: 10.1103/physreve.90.053310] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/31/2013] [Indexed: 06/04/2023]
Abstract
Unlike conventional computational fluid dynamics methods, the lattice Boltzmann method (LBM) describes the dynamic behavior of fluids in a mesoscopic scale based on discrete forms of kinetic equations. In this scale, complex macroscopic phenomena like the formation and collapse of interfaces can be naturally described as related to source terms incorporated into the kinetic equations. In this context, a novel athermal lattice Boltzmann scheme for the simulation of phase transition is proposed. The continuous kinetic model obtained from the Liouville equation using the mean-field interaction force approach is shown to be consistent with diffuse interface model using the Helmholtz free energy. Density profiles, interface thickness, and surface tension are analytically derived for a plane liquid-vapor interface. A discrete form of the kinetic equation is then obtained by applying the quadrature method based on prescribed abscissas together with a third-order scheme for the discretization of the streaming or advection term in the Boltzmann equation. Spatial derivatives in the source terms are approximated with high-order schemes. The numerical validation of the method is performed by measuring the speed of sound as well as by retrieving the coexistence curve and the interface density profiles. The appearance of spurious currents near the interface is investigated. The simulations are performed with the equations of state of Van der Waals, Redlich-Kwong, Redlich-Kwong-Soave, Peng-Robinson, and Carnahan-Starling.
Collapse
Affiliation(s)
- D N Siebert
- Federal University of Santa Catarina, 89218-000 Joinville, SC, Brazil
| | - P C Philippi
- Mechanical Engineering Department, Federal University of Santa Catarina, 88040-900 Florianópolis, SC, Brazil
| | - K K Mattila
- Mechanical Engineering Department, Federal University of Santa Catarina, 88040-900 Florianópolis, SC, Brazil
| |
Collapse
|
10
|
Li Q, Luo KH, Li XJ. Lattice Boltzmann modeling of multiphase flows at large density ratio with an improved pseudopotential model. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2013; 87:053301. [PMID: 23767651 DOI: 10.1103/physreve.87.053301] [Citation(s) in RCA: 90] [Impact Index Per Article: 8.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/08/2012] [Revised: 01/07/2013] [Indexed: 06/02/2023]
Abstract
Owing to its conceptual simplicity and computational efficiency, the pseudopotential multiphase lattice Boltzmann (LB) model has attracted significant attention since its emergence. In this work, we aim to extend the pseudopotential LB model to simulate multiphase flows at large density ratio and relatively high Reynolds number. First, based on our recent work [Q. Li, K. H. Luo, and X. J. Li, Phys. Rev. E 86, 016709 (2012)], an improved forcing scheme is proposed for the multiple-relaxation-time pseudopotential LB model in order to achieve thermodynamic consistency and large density ratio in the model. Next, through investigating the effects of the parameter a in the Carnahan-Starling equation of state, we find that the interface thickness is approximately proportional to 1/√a. Using a smaller a will lead to a wider interface thickness, which can reduce the spurious currents and enhance the numerical stability of the pseudopotential model at large density ratio. Furthermore, it is found that a lower liquid viscosity can be gained in the pseudopotential model by increasing the kinematic viscosity ratio between the vapor and liquid phases. The improved pseudopotential LB model is numerically validated via the simulations of stationary droplet and droplet oscillation. Using the improved model as well as the above treatments, numerical simulations of droplet splashing on a thin liquid film are conducted at a density ratio in excess of 500 with Reynolds numbers ranging from 40 to 1000. The dynamics of droplet splashing is correctly reproduced and the predicted spread radius is found to obey the power law reported in the literature.
Collapse
Affiliation(s)
- Q Li
- Energy Technology Research Group, Faculty of Engineering and the Environment, University of Southampton, Southampton SO17 1BJ, United Kingdom
| | | | | |
Collapse
|
11
|
Gross M, Varnik F. Simulation of static critical phenomena in nonideal fluids with the lattice Boltzmann method. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2012; 85:056707. [PMID: 23004903 DOI: 10.1103/physreve.85.056707] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/19/2011] [Indexed: 06/01/2023]
Abstract
A fluctuating nonideal fluid at its critical point is simulated with the lattice Boltzmann method. It is demonstrated that the method, employing a Ginzburg-Landau free energy functional, correctly reproduces the static critical behavior associated with the Ising universality class. A finite-size scaling analysis is applied to determine the critical exponents related to the order parameter, compressibility and specific heat. A particular focus is put on finite-size effects and issues related to the global conservation of the order parameter.
Collapse
Affiliation(s)
- M Gross
- Interdisciplinary Centre for Advanced Materials Simulation, Ruhr-Universität Bochum, Germany.
| | | |
Collapse
|
12
|
Melchionna S, Marini Bettolo Marconi U. Stabilized lattice Boltzmann-Enskog method for compressible flows and its application to one- and two-component fluids in nanochannels. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2012; 85:036707. [PMID: 22587209 DOI: 10.1103/physreve.85.036707] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/12/2011] [Indexed: 05/31/2023]
Abstract
A numerically stable method to solve the discretized Boltzmann-Enskog equation describing the behavior of nonideal fluids under inhomogeneous conditions is presented. The algorithm employed uses a Lagrangian finite-difference scheme for the treatment of the convective term and a forcing term to account for the molecular repulsion together with a Bhatnagar-Gross-Krook relaxation term. In order to eliminate the spurious currents induced by the numerical discretization procedure, we use a trapezoidal rule for the time integration together with a version of the two-distribution method of He et al. [J. Comput. Phys. 152, 642 (1999)]. Numerical tests show that, in the case of a one-component fluid in the presence of a spherical potential well, the proposed method reduces the numerical error by several orders of magnitude. We conduct another test by considering the flow of a two-component fluid in a channel with a bottleneck and provide information about the density and velocity field in this structured geometry.
Collapse
Affiliation(s)
- Simone Melchionna
- CNR-IPCF, Consiglio Nazionale delle Ricerche, Università di Roma La Sapienza, P.le A. Moro 2, 00185 Roma, Italy
| | | |
Collapse
|
13
|
Huang H, Krafczyk M, Lu X. Forcing term in single-phase and Shan-Chen-type multiphase lattice Boltzmann models. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2011; 84:046710. [PMID: 22181310 DOI: 10.1103/physreve.84.046710] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/29/2011] [Revised: 08/05/2011] [Indexed: 05/31/2023]
Abstract
Numerous schemes have been proposed to incorporate a bulk forcing term into the lattice Boltzmann equation. In this paper we present a simple and straightforward comparative analysis of five popular schemes [Shan and Chen, Phys. Rev. E 47, 1815 (1993); Phys Rev Lett. 81, 1618 (1998); He et al., Phys. Rev. E 57, R13 (1998); Guo et al., Phys. Rev. E 65, 046308 (2002); Kupershtokh et al., Comput. Math. Appl. 58, 965 (2009)] in which their differences and similarities are identified. From the analysis we classify the schemes into two groups; the behaviors of the schemes in each group are proven to be identical up to second order. Numerical test simulating the two-dimensional unsteady Taylor-Green vortex flow problem demonstrate that all five schemes are of comparable accuracy for single-phase flow. However, for two-phase flow the situation is different, which is demonstrated by incorporating these schemes into different Shan-Chen-type multiphase models. The forcing scheme in the original Shan-Chen (SC) multiphase model turns out to be inaccurate in terms of the resulting surface tension for different density ratios and relaxation times. In the numerical tests, a typical equation of state and interparticle interactions including next-nearest neighbors were incorporated into the SC model. Our results confirm that the surface-tension values obtained from the original SC lattice Boltzmann method (LBM) simulation depend on the value of the relaxation time τ. For τ<0.7Δt, the surface tension agree well with the analytical solutions. However, when τ>0.7Δt, the surface tension turns out to be systematically larger than the analytical one, exceeding it by more than a factor of 2 for τ=2Δt. In contrast, with the application of the scheme proposed by He et al., the SC LBM produces very accurate surface tensions independent of the value of τ. We also found that the densities of the coexisting liquid and gas can be adjusted to match those at thermodynamic equilibrium if the particle interaction term includes next-nearest-neighbor contributions. The obtained results will be useful for further studies of two-phase flow with high density ratios using the SC LBM approach.
Collapse
Affiliation(s)
- Haibo Huang
- Institute for Computational Modeling in Civil Engineering, Technische Universität, D-38106 Braunschweig, Germany
| | | | | |
Collapse
|
14
|
Gan Y, Xu A, Zhang G, Li Y, Li H. Phase separation in thermal systems: a lattice Boltzmann study and morphological characterization. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2011; 84:046715. [PMID: 22181315 DOI: 10.1103/physreve.84.046715] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/14/2010] [Revised: 09/11/2011] [Indexed: 05/31/2023]
Abstract
We investigate thermal and isothermal symmetric liquid-vapor separations via a fast Fourier transform thermal lattice Boltzmann (FFT-TLB) model. Structure factor, domain size, and Minkowski functionals are employed to characterize the density and velocity fields, as well as to understand the configurations and the kinetic processes. Compared with the isothermal phase separation, the freedom in temperature prolongs the spinodal decomposition (SD) stage and induces different rheological and morphological behaviors in the thermal system. After the transient procedure, both the thermal and isothermal separations show power-law scalings in domain growth, while the exponent for thermal system is lower than that for isothermal system. With respect to the density field, the isothermal system presents more likely bicontinuous configurations with narrower interfaces, while the thermal system presents more likely configurations with scattered bubbles. Heat creation, conduction, and lower interfacial stresses are the main reasons for the differences in thermal system. Different from the isothermal case, the release of latent heat causes the changing of local temperature, which results in new local mechanical balance. When the Prandtl number becomes smaller, the system approaches thermodynamical equilibrium much more quickly. The increasing of mean temperature makes the interfacial stress lower in the following way: σ=σ(0)[(T(c)-T)/(T(c)-T(0))](3/2), where T(c) is the critical temperature and σ(0) is the interfacial stress at a reference temperature T(0), which is the main reason for the prolonged SD stage and the lower growth exponent in the thermal case. Besides thermodynamics, we probe how the local viscosities influence the morphology of the phase separating system. We find that, for both the isothermal and thermal cases, the growth exponents and local flow velocities are inversely proportional to the corresponding viscosities. Compared with the isothermal case, the local flow velocity depends not only on viscosity but also on temperature.
Collapse
Affiliation(s)
- Yanbiao Gan
- State Key Laboratory for GeoMechanics and Deep Underground Engineering, SMCE, China University of Mining and Technology (Beijing), Beijing 100083, PR China
| | | | | | | | | |
Collapse
|
15
|
|
16
|
Guo Z, Shi B, Zheng C. Chequerboard effects on spurious currents in the lattice Boltzmann equation for two-phase flows. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2011; 369:2283-2291. [PMID: 21536575 DOI: 10.1098/rsta.2011.0023] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/30/2023]
Abstract
Spurious currents near an interface between different phases are a common undesirable feature of the lattice Boltzmann equation (LBE) method for two-phase systems. In this paper, we show that the spurious currents of a kinetic theory-based LBE have a significant dependence on the parity of the grid number of the underlying lattice, which can be attributed to the chequerboard effect. A technique that uses a Lax-Wendroff streaming is proposed to overcome this anomaly, and its performance is verified numerically.
Collapse
Affiliation(s)
- Zhaoli Guo
- State Key Laboratory of Enhanced Oil Recovery, Research Institute of Petroleum Exploration and Development, Beijing 100083, People's Republic of China.
| | | | | |
Collapse
|
17
|
Gross M, Moradi N, Zikos G, Varnik F. Shear stress in nonideal fluid lattice Boltzmann simulations. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2011; 83:017701. [PMID: 21405797 DOI: 10.1103/physreve.83.017701] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/20/2010] [Indexed: 05/30/2023]
Abstract
The viscous stress in nonideal fluid lattice Boltzmann methods is investigated theoretically and by simulations. Three representative liquid-gas models are compared in a steady-state situation, where an analytical expression for the viscous stress is available. It is shown that, in the presence of nonideal fluid interactions or strong body forces, the accuracy of the computed viscous stress depends strongly on the underlying model.
Collapse
Affiliation(s)
- Markus Gross
- Interdisciplinary Centre for Advanced Materials Simulation, Ruhr-Universität Bochum, Stiepeler Strasse 129, D-44801 Bochum, Germany.
| | | | | | | |
Collapse
|
18
|
Gray WG, Miller CT. Thermodynamically Constrained Averaging Theory Approach for Modeling Flow and Transport Phenomena in Porous Medium Systems: 8. Interface and Common Curve Dynamics. ADVANCES IN WATER RESOURCES 2010; 33:1427-1443. [PMID: 21197134 PMCID: PMC3010759 DOI: 10.1016/j.advwatres.2010.07.002] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/24/2023]
Abstract
This work is the eighth in a series that develops the fundamental aspects of the thermodynamically constrained averaging theory (TCAT) that allows for a systematic increase in the scale at which multiphase transport phenomena is modeled in porous medium systems. In these systems, the explicit locations of interfaces between phases and common curves, where three or more interfaces meet, are not considered at scales above the microscale. Rather, the densities of these quantities arise as areas per volume or length per volume. Modeling of the dynamics of these measures is an important challenge for robust models of flow and transport phenomena in porous medium systems, as the extent of these regions can have important implications for mass, momentum, and energy transport between and among phases, and formulation of a capillary pressure relation with minimal hysteresis. These densities do not exist at the microscale, where the interfaces and common curves correspond to particular locations. Therefore, it is necessary for a well-developed macroscale theory to provide evolution equations that describe the dynamics of interface and common curve densities. Here we point out the challenges and pitfalls in producing such evolution equations, develop a set of such equations based on averaging theorems, and identify the terms that require particular attention in experimental and computational efforts to parameterize the equations. We use the evolution equations developed to specify a closed two-fluid-phase flow model.
Collapse
|
19
|
Kikkinides ES, Kainourgiakis ME, Yiotis AG, Stubos AK. Lattice Boltzmann method for Lennard-Jones fluids based on the gradient theory of interfaces. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2010; 82:056705. [PMID: 21230617 DOI: 10.1103/physreve.82.056705] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/12/2010] [Revised: 09/16/2010] [Indexed: 05/30/2023]
Abstract
In the present study we propose a lattice Boltzmann equation (LBE) model derived from density gradient expansions of the discrete BBGKY evolution equations. The model is based on the mechanical approach of the gradient theory of interfaces. The basic input is the radial distribution function, which is related exclusively to the molecular interaction potential, rather than semiempirical equations of state used in previous LBE models. This function can be provided from independent molecular simulations or from approximate theories. Evidently the accuracy of the interaction potential, and thus the radial distribution function, reflects on the accuracy of the thermodynamic properties and consistency of the derived LBE model. We have applied the proposed model to obtain equilibrium bulk and interfacial properties of a Lennard-Jones fluid at different temperatures, T, close to critical, T(c). The results demonstrate that the LBE model is in excellent agreement with gradient theory as well as with independent literature results based on different molecular simulation approaches. Hence the proposed LBE model can recover accurately bulk and interfacial thermodynamics for a Lennard Jones fluid at T/T(c)>0.9.
Collapse
Affiliation(s)
- E S Kikkinides
- Department of Mechanical Engineering, University of Western Macedonia, 50100 Kozani, Greece
| | | | | | | |
Collapse
|