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.
Scheda prodotto non validato
Attenzione! I dati visualizzati non sono stati sottoposti a validazione da parte dell'ateneo
Titolo: | Stabilized Lattice Boltzmann-Enskog method for compressible flows and its application to one and two-component fluids in nanochannels | |
Autori: | ||
Data di pubblicazione: | 2012 | |
Rivista: | ||
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. | |
Handle: | http://hdl.handle.net/11581/244342 | |
Appare nelle tipologie: | Articolo |