1
|
Soomro M, Ayala LF. Unrestricted component count in multiphase lattice Boltzmann: A fugacity-based approach. Phys Rev E 2023; 108:035304. [PMID: 37849190 DOI: 10.1103/physreve.108.035304] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2023] [Accepted: 08/07/2023] [Indexed: 10/19/2023]
Abstract
Studies of multiphase fluids utilizing the lattice Boltzmann method (LBM) are typically severely restricted by the number of components or chemical species being modeled. This restriction is particularly pronounced for multiphase systems exhibiting partial miscibility and significant interfacial mass exchange, which is a common occurrence in realistic multiphase systems. Modeling such systems becomes increasingly complex as the number of chemical species increases due to the increased role of molecular interactions and the types of thermodynamic behavior that become possible. The recently introduced fugacity-based LBM [Soomro et al., Phys. Rev. E 107, 015304 (2023)2470-004510.1103/PhysRevE.107.015304] has provided a thermodynamically consistent modeling platform for multicomponent, partially miscible LBM simulations. However, until now, this fugacity-based LB model had lacked a comprehensive demonstration of its ability to accurately reproduce thermodynamic behavior beyond binary mixtures and to remove any restrictions in a number of components for multiphase LBM. In this paper we closely explore these fugacity-based LBM capabilities by showcasing comprehensive, thermodynamically consistent simulations of multiphase mixtures of up to ten chemical components. The paper begins by validating the model against the Young-Laplace equation for a droplet composed of three components. The model is then applied to study mixtures with a range of component numbers from one to six, showing agreement with rigorous thermodynamic predictions and demonstrating linear scaling of computational time with the number of components. We further investigate ternary systems in detail by exploring a wide range of temperature, pressure, and overall composition conditions to produce various characteristic ternary diagrams. In addition, the model is shown to be unrestricted in the number of phases as demonstrated through simulations of a three-component three-phase equilibrium case. The paper concludes by demonstrating simulations of a ten-component, realistic hydrocarbon mixture, achieving excellent agreement with thermodynamics for both flat interface vapor-liquid equilibrium and curved interface spinodal decomposition cases. This study represents a significant expansion of the scope and capabilities of multiphase LBM simulations that encompass multiphase systems of keen interest in engineering.
Collapse
Affiliation(s)
- Muzammil Soomro
- Department of Energy and Mineral Engineering, The Pennsylvania State University, University Park, Pennsylvania 16802, USA
| | - Luis F Ayala
- Department of Energy and Mineral Engineering, The Pennsylvania State University, University Park, Pennsylvania 16802, USA
| |
Collapse
|
2
|
Soomro M, Ayala LF, Peng C, Ayala OM. Fugacity-based lattice Boltzmann method for multicomponent multiphase systems. Phys Rev E 2023; 107:015304. [PMID: 36797960 DOI: 10.1103/physreve.107.015304] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2022] [Accepted: 12/13/2022] [Indexed: 06/18/2023]
Abstract
The free-energy model can extend the lattice Boltzmann method to multiphase systems. However, there is a lack of models capable of simulating multicomponent multiphase fluids with partial miscibility. In addition, existing models cannot be generalized to honor thermodynamic information provided by any multicomponent equation of state of choice. In this paper, we introduce a free-energy lattice Boltzmann model where the forcing term is determined by the fugacity of the species, the thermodynamic property that connects species partial pressure to chemical potential calculations. By doing so, we are able to carry out multicomponent multiphase simulations of partially miscible fluids and generalize the methodology for use with any multicomponent equation of state of interest. We test this fugacity-based lattice Boltzmann method for the cases of vapor-liquid equilibrium for two- and three-component mixtures in various temperature and pressure conditions. We demonstrate that the model is able to reliably reproduce phase densities and compositions as predicted by multicomponent thermodynamics and can reproduce different characteristic pressure-composition and temperature-composition envelopes with a high degree of accuracy. We also demonstrate that the model can offer accurate predictions under dynamic conditions.
Collapse
Affiliation(s)
- Muzammil Soomro
- Department of Energy and Mineral Engineering, The Pennsylvania State University, University Park, Pennsylvania 16802, USA
| | - Luis F Ayala
- Department of Energy and Mineral Engineering, The Pennsylvania State University, University Park, Pennsylvania 16802, USA
| | - Cheng Peng
- Key Laboratory of High Efficiency and Clean Mechanical Manufacture, Ministry of Education, School of Mechanical Engineering, Shandong University, Jinan 250061, China
| | - Orlando M Ayala
- Department of Engineering Technology, Old Dominion University, Norfolk, Virginia 23529, USA
| |
Collapse
|
3
|
Lourenço RGC, Constantino PH, Tavares FW. A Unified Interaction Model for Multiphase Flows with the Lattice Boltzmann Method. CAN J CHEM ENG 2022. [DOI: 10.1002/cjce.24604] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Affiliation(s)
- Ramon G. C. Lourenço
- Program of Chemical Engineering /COPPE ‐ Federal University of Rio de Janeiro, CEP: 21949‐972 Rio de Janeiro RJ Brazil
| | - Pedro H. Constantino
- Program of Chemical Engineering /COPPE ‐ Federal University of Rio de Janeiro, CEP: 21949‐972 Rio de Janeiro RJ Brazil
| | - Frederico W. Tavares
- Program of Chemical Engineering /COPPE ‐ Federal University of Rio de Janeiro, CEP: 21949‐972 Rio de Janeiro RJ Brazil
- Program in Engineering of Chemical and Biochemical Processes ‐ Chemical School Federal University of Rio de Janeiro, CEP Rio de Janeiro RJ Brazil
| |
Collapse
|
4
|
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
|
5
|
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
|
6
|
Chiappini D, Sbragaglia M, Xue X, Falcucci G. Hydrodynamic behavior of the pseudopotential lattice Boltzmann method for interfacial flows. Phys Rev E 2019; 99:053305. [PMID: 31212544 DOI: 10.1103/physreve.99.053305] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2018] [Indexed: 06/09/2023]
Abstract
The lattice Boltzmann method (LBM) is routinely employed in the simulation of complex multiphase flows comprising bulk phases separated by nonideal interfaces. The LBM is intrinsically mesoscale with a hydrodynamic equivalence popularly set by the Chapman-Enskog analysis, requiring that fields slowly vary in space and time. The latter assumptions become questionable close to interfaces where the method is also known to be affected by spurious nonhydrodynamical contributions. This calls for quantitative hydrodynamical checks. In this paper, we analyze the hydrodynamic behavior of the LBM pseudopotential models for the problem of the breakup of a liquid ligament triggered by the Plateau-Rayleigh instability. Simulations are performed at fixed interface thickness, while increasing the ligament radius, i.e., in the "sharp interface" limit. The influence of different LBM collision operators is also assessed. We find that different distributions of spurious currents along the interface may change the outcome of the pseudopotential model simulations quite sensibly, which suggests that a proper fine-tuning of pseudopotential models in time-dependent problems is needed before the utilization in concrete applications. Taken all together, we argue that the results of the proposed paper provide a valuable insight for engineering pseudopotential model applications involving the hydrodynamics of liquid jets.
Collapse
Affiliation(s)
- Daniele Chiappini
- Department of Industrial Engineering, University of Rome "Niccolò Cusano," Via don Carlo Gnocchi 3, 00166 Rome, Italy
| | - Mauro Sbragaglia
- Department of Physics, INFN, University of Rome "Tor Vergata," Via della Ricerca Scientifica 1, 00133 Rome, Italy
| | - Xiao Xue
- Department of Physics, INFN, University of Rome "Tor Vergata," Via della Ricerca Scientifica 1, 00133 Rome, Italy and Department of Physics, Eindhoven University of Technology, 5600 MB Eindhoven, The Netherlands
| | - Giacomo Falcucci
- Department of Enterprise Engineering "Mario Lucertini," University of Rome "Tor Vergata," Via del Politecnico 1, 00133 Rome, Italy and John A. Paulson School of Engineering and Applied Physics, Harvard University, 33 Oxford Street, 02138 Cambridge, Massachusetts, USA
| |
Collapse
|
7
|
Sofonea V, Biciuşcă T, Busuioc S, Ambruş VE, Gonnella G, Lamura A. Corner-transport-upwind lattice Boltzmann model for bubble cavitation. Phys Rev E 2018; 97:023309. [PMID: 29548242 DOI: 10.1103/physreve.97.023309] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/24/2015] [Indexed: 11/07/2022]
Abstract
Aiming to study the bubble cavitation problem in quiescent and sheared liquids, a third-order isothermal lattice Boltzmann model that describes a two-dimensional (2D) fluid obeying the van der Waals equation of state, is introduced. The evolution equations for the distribution functions in this off-lattice model with 16 velocities are solved using the corner-transport-upwind (CTU) numerical scheme on large square lattices (up to 6144×6144 nodes). The numerical viscosity and the regularization of the model are discussed for first- and second-order CTU schemes finding that the latter choice allows to obtain a very accurate phase diagram of a nonideal fluid. In a quiescent liquid, the present model allows us to recover the solution of the 2D Rayleigh-Plesset equation for a growing vapor bubble. In a sheared liquid, we investigated the evolution of the total bubble area, the bubble deformation, and the bubble tilt angle, for various values of the shear rate. A linear relation between the dimensionless deformation coefficient D and the capillary number Ca is found at small Ca but with a different factor than in equilibrium liquids. A nonlinear regime is observed for Ca≳0.2.
Collapse
Affiliation(s)
- V Sofonea
- Center for Fundamental and Advanced Technical Research, Romanian Academy, Bd. Mihai Viteazul 24, 300223 Timişoara, Romania
| | - T Biciuşcă
- Center for Fundamental and Advanced Technical Research, Romanian Academy, Bd. Mihai Viteazul 24, 300223 Timişoara, Romania.,Department of Physics, West University of Timişoara, Bd. Vasile Pârvan 4, 300223 Timişoara, Romania
| | - S Busuioc
- Center for Fundamental and Advanced Technical Research, Romanian Academy, Bd. Mihai Viteazul 24, 300223 Timişoara, Romania.,Department of Physics, West University of Timişoara, Bd. Vasile Pârvan 4, 300223 Timişoara, Romania
| | - Victor E Ambruş
- Center for Fundamental and Advanced Technical Research, Romanian Academy, Bd. Mihai Viteazul 24, 300223 Timişoara, Romania.,Department of Physics, West University of Timişoara, Bd. Vasile Pârvan 4, 300223 Timişoara, Romania
| | - G Gonnella
- Dipartimento di Fisica, Università di Bari, and INFN, Sezione di Bari, Via Amendola 173, 70126 Bari, Italy
| | - A Lamura
- Istituto Applicazioni Calcolo, CNR, Via Amendola 122/D, 70126 Bari, Italy
| |
Collapse
|