共查询到20条相似文献,搜索用时 15 毫秒
1.
Jaw-Yen Yang Bagus Putra Muljadi Zhi-Hui Li & Han-Xin Zhang 《Communications In Computational Physics》2013,14(1):242-264
An accurate and direct algorithm for solving the semiclassical Boltzmann
equation with relaxation time approximation in phase space is presented for parallel
treatment of rarefied gas flows of particles of three statistics. The discrete ordinate
method is first applied to discretize the velocity space of the distribution function to
render a set of scalar conservation laws with source term. The high order weighted essentially non-oscillatory scheme is then implemented to capture the time evolution of
the discretized velocity distribution function in physical space and time. The method
is developed for two space dimensions and implemented on gas particles that obey
the Maxwell-Boltzmann, Bose-Einstein and Fermi-Dirac statistics. Computational examples in one- and two-dimensional initial value problems of rarefied gas flows are
presented and the results indicating good resolution of the main flow features can be
achieved. Flows of wide range of relaxation times and Knudsen numbers covering
different flow regimes are computed to validate the robustness of the method. The
recovery of quantum statistics to the classical limit is also tested for small fugacity
values. 相似文献
2.
V. A. Titarev 《Communications In Computational Physics》2012,12(1):162-192
The paper is devoted to the development of an efficient deterministic framework for modelling of three-dimensional rarefied gas flows on the basis of the numerical solution of the Boltzmann kinetic equation with the model collision integrals. The
framework consists of a high-order accurate implicit advection scheme on arbitrary
unstructured meshes, the conservative procedure for the calculation of the model collision integral and efficient implementation on parallel machines. The main application
area of the suggested methods is micro-scale flows. Performance of the proposed approach is demonstrated on a rarefied gas flow through the finite-length circular pipe.
The results show good accuracy of the proposed algorithm across all flow regimes and
its high efficiency and excellent parallel scalability for up to 512 cores. 相似文献
3.
High-Order Conservative Asymptotic-Preserving Schemes for Modeling Rarefied Gas Dynamical Flows with Boltzmann-BGK Equation
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
Manuel A. Diaz Min-Hung Chen & Jaw-Yen Yang 《Communications In Computational Physics》2015,18(4):1012-1049
High-order and conservative phase space direct solvers that preserve the
Euler asymptotic limit of the Boltzmann-BGK equation for modelling rarefied gas
flows are explored and studied. The approach is based on the conservative discrete ordinate
method for velocity space by using Gauss Hermite or Simpsons quadrature rule
and conservation of macroscopic properties are enforced on the BGK collision operator.
High-order asymptotic-preserving time integration is adopted and the spatial evolution
is performed by high-order schemes including a finite difference weighted essentially
non-oscillatory method and correction procedure via reconstruction schemes.
An artificial viscosity dissipative model is introduced into the Boltzmann-BGK equation
when the correction procedure via reconstruction scheme is used. The effects of
the discrete velocity conservative property and accuracy of high-order formulations of
kinetic schemes based on BGK model methods are provided. Extensive comparative
tests with one-dimensional and two-dimensional problems in rarefied gas flows have
been carried out to validate and illustrate the schemes presented. Potentially advantageous
schemes in terms of stable large time step allowed and higher-order of accuracy
are suggested. 相似文献
4.
With discretized particle velocity space, a multi-scale unified gas-kinetic
scheme for entire Knudsen number flows has been constructed based on the kinetic
model in one-dimensional case [J. Comput. Phys., vol. 229 (2010), pp. 7747-7764]. For
the kinetic equation, to extend a one-dimensional scheme to multidimensional flow is
not so straightforward. The major factor is that addition of one dimension in physical
space causes the distribution function to become two-dimensional, rather than axially
symmetric, in velocity space. In this paper, a unified gas-kinetic scheme based on the
Shakhov model in two-dimensional space will be presented. Instead of particle-based
modeling for the rarefied flow, such as the direct simulation Monte Carlo (DSMC)
method, the philosophical principal underlying the current study is a partial-differential-equation (PDE)-based modeling. Since the valid scale of the kinetic equation and
the scale of mesh size and time step may be significantly different, the gas evolution in a discretized space is modeled with the help of kinetic equation, instead of
directly solving the partial differential equation. Due to the use of both hydrodynamic and kinetic scales flow physics in a gas evolution model at the cell interface,
the unified scheme can basically present accurate solution in all flow regimes from
the free molecule to the Navier-Stokes solutions. In comparison with the DSMC and
Navier-Stokes flow solvers, the current method is much more efficient than DSMC in
low speed transition and continuum flow regimes, and it has better capability than
NS solver in capturing of non-equilibrium flow physics in the transition and rarefied
flow regimes. As a result, the current method can be useful in the flow simulation
where both continuum and rarefied flow physics needs to be resolved in a single computation. This paper will extensively evaluate the performance of the unified scheme
from free molecule to continuum NS solutions, and from low speed micro-flow to high
speed non-equilibrium aerodynamics. The test cases clearly demonstrate that the unified scheme is a reliable method for the rarefied flow computations, and the scheme
provides an important tool in the study of non-equilibrium flow. 相似文献
5.
S. C. Fu R. M. C. So & W. W. F. Leung 《Communications In Computational Physics》2011,9(5):1257-1283
The objective of this paper is to seek an alternative to the numerical simulation
of the Navier-Stokes equations by a method similar to solving the BGK-type
modeled lattice Boltzmann equation. The proposed method is valid for both gas and
liquid flows. A discrete flux scheme (DFS) is used to derive the governing equations
for two distribution functions; one for mass and another for thermal energy. These
equations are derived by considering an infinitesimally small control volume with a
velocity lattice representation for the distribution functions. The zero-order moment
equation of the mass distribution function is used to recover the continuity equation,
while the first-order moment equation recovers the linear momentum equation. The
recovered equations are correct to the first order of the Knudsen number (Kn); thus,
satisfying the continuum assumption. Similarly, the zero-order moment equation of
the thermal energy distribution function is used to recover the thermal energy equation.
For aerodynamic flows, it is shown that the finite difference solution of the DFS
is equivalent to solving the lattice Boltzmann equation (LBE) with a BGK-type model
and a specified equation of state. Thus formulated, the DFS can be used to simulate a
variety of aerodynamic and hydrodynamic flows. Examples of classical aeroacoustics,
compressible flow with shocks, incompressible isothermal and non-isothermal Couette
flows, stratified flow in a cavity, and double diffusive flow inside a rectangle are used
to demonstrate the validity and extent of the DFS. Very good to excellent agreement
with known analytical and/or numerical solutions is obtained; thus lending evidence
to the DFS approach as an alternative to solving the Navier-Stokes equations for fluid
flow simulations. 相似文献
6.
Jaw-Yen Yang Li-Hsin Hung & Yao-Tien Kuo 《Communications In Computational Physics》2011,10(2):405-421
Computations of microscopic circular pipe flow in a rarefied quantum gas
are presented using a semiclassical axisymmetric lattice Boltzmann method. The
method is first derived by directly projecting the Uehling-Uhlenbeck Boltzmann-BGK
equations in two-dimensional rectangular coordinates onto the tensor Hermite polynomials
using moment expansion method and then the forcing strategy of Halliday
et al. [Phys. Rev. E., 64 (2001), 011208] is adopted by adding forcing terms into the
resulting microdynamic evolution equation. The determination of the forcing terms
is dictated by yielding the emergent macroscopic equations toward a particular target
form. The correct macroscopic equations of the incompressible axisymmetric viscous
flows are recovered through the Chapman-Enskog expansion. The velocity profiles
and the mass flow rates of pipe flows with several Knudsen numbers covering different
flow regimes are presented. It is found the Knudsen minimum can be captured in
all three statistics studied. The results also indicate distinct characteristics of the effects
of quantum statistics. 相似文献
7.
A Novel Dynamic Quadrature Scheme for Solving Boltzmann Equation with Discrete Ordinate and Lattice Boltzmann Methods
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
The Boltzmann equation (BE) for gas flows is a time-dependent nonlinear
differential-integral equation in 6 dimensions. The current simplified practice is to linearize the collision integral in BE by the BGK model using Maxwellian equilibrium
distribution and to approximate the moment integrals by the discrete ordinate method
(DOM) using a finite set of velocity quadrature points. Such simplification reduces
the dimensions from 6 to 3, and leads to a set of linearized discrete BEs. The main
difficulty of the currently used (conventional) numerical procedures occurs when the
mean velocity and the variation of temperature are large that requires an extremely
large number of quadrature points. In this paper, a novel dynamic scheme that requires only a small number of quadrature points is proposed. This is achieved by
a velocity-coordinate transformation consisting of Galilean translation and thermal
normalization so that the transformed velocity space is independent of mean velocity and temperature. This enables the efficient implementation of Gaussian-Hermite
quadrature. The velocity quadrature points in the new velocity space are fixed while
the correspondent quadrature points in the physical space change from time to time
and from position to position. By this dynamic nature in the physical space, this new
quadrature scheme is termed as the dynamic quadrature scheme (DQS). The DQS was
implemented to the DOM and the lattice Boltzmann method (LBM). These new methods with DQS are therefore termed as the dynamic discrete ordinate method (DDOM)
and the dynamic lattice Boltzmann method (DLBM), respectively. The new DDOM
and DLBM have been tested and validated with several testing problems. Of the same
accuracy in numerical results, the proposed schemes are much faster than the conventional schemes. Furthermore, the new DLBM have effectively removed the incompressible and isothermal restrictions encountered by the conventional LBM. 相似文献
8.
Due to the rapid advances in micro-electro-mechanical systems (MEMS), the
study of microflows becomes increasingly important. Currently, the molecular-based
simulation techniques are the most reliable methods for rarefied flow computation,
even though these methods face statistical scattering problem in the low speed limit.
With discretized particle velocity space, a unified gas-kinetic scheme (UGKS) for entire Knudsen number flow has been constructed recently for flow computation. Contrary to the particle-based direct simulation Monte Carlo (DSMC) method, the unified
scheme is a partial differential equation-based modeling method, where the statistical
noise is totally removed. But the common point between the DSMC and UGKS is that
both methods are constructed through direct modeling in the discretized space. Due
to the multiscale modeling in the unified method, i.e., the update of both macroscopic
flow variables and microscopic gas distribution function, the conventional constraint
of time step being less than the particle collision time in many direct Boltzmann solvers
is released here. The numerical tests show that the unified scheme is more efficient
than the particle-based methods in the low speed rarefied flow computation. The main
purpose of the current study is to validate the accuracy of the unified scheme in the
capturing of non-equilibrium flow phenomena. In the continuum and free molecular
limits, the gas distribution function used in the unified scheme for the flux evaluation
at a cell interface goes to the corresponding Navier-Stokes and free molecular solutions. In the transition regime, the DSMC solution will be used for the validation of
UGKS results. This study shows that the unified scheme is indeed a reliable and accurate flow solver for low speed non-equilibrium flows. It not only recovers the DSMC
results whenever available, but also provides high resolution results in cases where
the DSMC can hardly afford the computational cost. In thermal creep flow simulation,
surprising solution, such as the gas flowing from hot to cold regions along the wallsurface, is observed for the first time by the unified scheme, which is confirmed later
through intensive DSMC computation. 相似文献
9.
Peng Wang Lianhua Zhu Zhaoli Guo & Kun Xu 《Communications In Computational Physics》2015,17(3):657-681
The lattice Boltzmann equation (LBE) methods (both LBGK and MRT) and
the discrete unified gas-kinetic scheme (DUGKS) are both derived from the Boltzmann
equation, but with different consideration in their algorithm construction. With the
same numerical discretization in the particle velocity space, the distinctive modeling
of these methods in the update of gas distribution function may introduce differences
in the computational results. In order to quantitatively evaluate the performance of
these methods in terms of accuracy, stability, and efficiency, in this paper we test LBGK,
MRT, and DUGKS in two-dimensional cavity flow and the flow over a square cylinder,
respectively. The results for both cases are validated against benchmark solutions.
The numerical comparison shows that, with sufficient mesh resolution, the LBE and
DUGKS methods yield qualitatively similar results in both test cases. With identical
mesh resolutions in both physical and particle velocity space, the LBE methods are
more efficient than the DUGKS due to the additional particle collision modeling in
DUGKS. But the DUGKS is more robust and accurate than the LBE methods in most
test conditions. Particularly, for the unsteady flow over a square cylinder at Reynolds
number 300, with the same mesh resolution it is surprisingly observed that the DUGKS
can capture the physical multi-frequency vortex shedding phenomena while the LBGK
and MRT fail to get that. Furthermore, the DUGKS is a finite volume method and
its computational efficiency can be much improved when a non-uniform mesh in the
physical space is adopted. The comparison in this paper clearly demonstrates the progressive
improvement of the lattice Boltzmann methods from LBGK, to MRT, up to the
current DUGKS, along with the inclusion of more reliable physical process in their algorithm
development. Besides presenting the Navier-Stokes solution, the DUGKS can
capture the rarefied flow phenomena as well with the increasing of Knudsen number. 相似文献
10.
G H. Tang G. X. Zhai W. Q. Tao X. J. Gu & D. R. Emerson 《Communications In Computational Physics》2013,13(5):1330-1356
Gases in microfluidic structures or devices are often in a non-equilibrium
state. The conventional thermodynamic models for fluids and heat transfer break
down and the Navier-Stokes-Fourier equations are no longer accurate or valid. In
this paper, the extended thermodynamic approach is employed to study the rarefied
gas flow in microstructures, including the heat transfer between a parallel channel and
pressure-driven Poiseuille flows through a parallel microchannel and circular microtube. The gas flow characteristics are studied and it is shown that the heat transfer in
the non-equilibrium state no longer obeys the Fourier gradient transport law. In addition, the bimodal distribution of streamwise and spanwise velocity and temperature
through a long circular microtube is captured for the first time. 相似文献
11.
This paper develops a high-order accurate gas-kinetic scheme in the framework
of the finite volume method for the one- and two-dimensional flow simulations,
which is an extension of the third-order accurate gas-kinetic scheme [Q.B. Li, K. Xu, and
S. Fu, J. Comput. Phys., 229(2010), 6715-6731] and the second-order accurate gas-kinetic
scheme [K. Xu, J. Comput. Phys., 171(2001), 289-335]. It is formed by two parts: quartic
polynomial reconstruction of the macroscopic variables and fourth-order accurate flux
evolution. The first part reconstructs a piecewise cell-center based quartic polynomial
and a cell-vertex based quartic polynomial according to the "initial" cell average approximation
of macroscopic variables to recover locally the non-equilibrium and equilibrium
single particle velocity distribution functions around the cell interface. It is in
view of the fact that all macroscopic variables become moments of a single particle velocity
distribution function in the gas-kinetic theory. The generalized moment limiter
is employed there to suppress the possible numerical oscillation. In the second part,
the macroscopic flux at the cell interface is evolved in fourth-order accuracy by means
of the simple particle transport mechanism in the microscopic level, i.e. free transport
and the Bhatnagar-Gross-Krook (BGK) collisions. In other words, the fourth-order
flux evolution is based on the solution (i.e. the particle velocity distribution function)
of the BGK model for the Boltzmann equation. Several 1D and 2D test problems are
numerically solved by using the proposed high-order accurate gas-kinetic scheme. By
comparing with the exact solutions or the numerical solutions obtained the second-order
or third-order accurate gas-kinetic scheme, the computations demonstrate that
our scheme is effective and accurate for simulating invisid and viscous fluid flows,
and the accuracy of the high-order GKS depends on the choice of the (numerical) collision
time. 相似文献
12.
13.
A Conservative Numerical Method for the Cahn–Hilliard Equation with Generalized Mobilities on Curved Surfaces in Three-Dimensional Space
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
Darae Jeong Yibao Li Chaeyoung Lee Junxiang Yang & Junseok Kim 《Communications In Computational Physics》2020,27(2):412-430
In this paper, we develop a conservative numerical method for the Cahn–
Hilliard equation with generalized mobilities on curved surfaces in three-dimensional
space. We use an unconditionally gradient stable nonlinear splitting numerical scheme
and solve the resulting system of implicit discrete equations on a discrete narrow band
domain by using a Jacobi-type iteration. For the domain boundary cells, we use the
trilinear interpolation using the closest point method. The proposing numerical algorithm is computationally efficient because we can use the standard finite difference
Laplacian scheme on three-dimensional Cartesian narrow band mesh instead of discrete Laplace–Beltrami operator on triangulated curved surfaces. In particular, we employ a mass conserving correction scheme, which enforces conservation of total mass.
We perform numerical experiments on the various curved surfaces such as sphere,
torus, bunny, cube, and cylinder to demonstrate the performance and effectiveness of
the proposed method. We also present the dynamics of the CH equation with constant
and space-dependent mobilities on the curved surfaces. 相似文献
14.
Extension of Near-Wall Domain Decomposition to Modeling Flows with Laminar-Turbulent Transition
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
M. Petrov S. Utyuzhnikov A. Chikitkin & N. Smirnova 《Communications In Computational Physics》2022,31(2):645-668
The near-wall domain decomposition method (NDD) has proved to be very
efficient for modeling near-wall fully turbulent flows. In this paper the NDD is extended to non-equilibrium regimes with laminar-turbulent transition (LTT) for the first
time. The LTT is identified with the use of the $e^N$-method which is applied to both incompressible and compressible flows. The NDD is modified to take into account LTT in
an efficient way. In addition, implementation of the intermittency expands the capabilities of NDD to model non-equilibrium turbulent flows with transition. Performance
of the modified NDD approach is demonstrated on various test problems of subsonic
and supersonic flows past a flat plate, a supersonic flow over a compression corner and
a planar shock wave impinging on a turbulent boundary layer. The results of modeling with and without decomposition are compared in terms of wall friction and show
good agreement with each other while NDD significantly reducing computational resources needed. It turns out that the NDD can reduce the computational time as much
as three times while retaining practically the same accuracy of prediction. 相似文献
15.
Simulation of Power-Law Fluid Flows in Two-Dimensional Square Cavity Using Multi-Relaxation-Time Lattice Boltzmann Method
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
Qiuxiang Li Ning Hong Baochang Shi & Zhenhua Chai 《Communications In Computational Physics》2014,15(1):265-284
In this paper, the power-law fluid flows in a two-dimensional square cavity
are investigated in detail with multi-relaxation-time lattice Boltzmann method (MRT-LBM). The influence of the Reynolds number (Re) and the power-law index (n) on the
vortex strength, vortex position and velocity distribution are extensively studied. In
our numerical simulations, Re is varied from 100 to 10000, and n is ranged from 0.25 to
1.75, covering both cases of shear-thinning and shear-thickening. Compared with the
Newtonian fluid, numerical results show that the flow structure and number of vortex
of power-law fluid are not only dependent on the Reynolds number, but also related
to power-law index. 相似文献
16.
In this paper the pressure distribution of the gaseous flow in a microchannel is studied via a lattice Boltzmann equation (LBE) method. With effective relaxation times and a generalized second order slip boundary condition, the LBE can be used to simulate rarefied gas flows from slip to transition regimes. The Knudsen minimum phenomena of mass flow rate in the pressure driven flow is also investigated. The effects of Knudsen number (rarefaction effect), pressure ratio and aspect ratio (compression effect) on the pressure distribution are analyzed. It is found the rarefaction effect tends to the curvature of the nonlinear pressure distribution, while the compression effect tends to enhance its nonlinearity. The combined effects lead to a local minimum of the pressure deviation. Furthermore, it is also found that the relationship between the pressure deviation and the aspect ratio follows a pow-law. 相似文献
17.
Are extensions to continuum formulations for solving fluid dynamic problems in the transition-to-rarefied regimes viable alternatives to particle methods? It
is well known that for increasingly rarefied flow fields, the predictions from continuum
formulation, such as the Navier-Stokes equations lose accuracy. These inaccuracies are
attributed primarily to the linear approximations of the stress and heat flux terms in the
Navier-Stokes equations. The inclusion of higher-order terms, such as Burnett or high-order moment equations, could improve the predictive capabilities of such continuum
formulations, but there has been limited success in the shock structure calculations, especially for the high Mach number case. Here, after reformulating the viscosity and heat
conduction coefficients appropriate for the rarefied flow regime, we will show that the
Navier-Stokes-type continuum formulation may still be properly used. The equations
with generalization of the dissipative coefficients based on the closed solution of the
Bhatnagar-Gross-Krook (BGK) model of the Boltzmann equation, are solved using the
gas-kinetic numerical scheme. This paper concentrates on the non-equilibrium shock
structure calculations for both monatomic and diatomic gases. The Landau-Teller-Jeans
relaxation model for the rotational energy is used to evaluate the quantitative difference
between the translational and rotational temperatures inside the shock layer. Variations
of shear stress, heat flux, temperatures, and densities in the internal structure of the
shock waves are compared with, (a) existing theoretical solutions of the Boltzmann solution, (b) existing numerical predictions of the direct simulation Monte Carlo (DSMC)
method, and (c) available experimental measurements. The present continuum formulation for calculating the shock structures for monatomic and diatomic gases in the
Mach number range of 1.2 to 12.9 is found to be satisfactory. 相似文献
18.
19.
In this paper, a gas-kinetic scheme (GKS) method coupled with a three temperature kinetic model is presented and applied in numerical study of the Edney-type IV shock/shock interaction which could cause serious problems in hypersonic
vehicles. The results showed very good agreement with the experimental data in
predicting the heat flux on the surface. It could be obviously seen that the current
method can accurately describe the position and features of supersonic jets structure
and clearly capture the thermal non-equilibrium in this case. The three temperature
kinetic model includes three different models of temperatures which are translational,
rotational and vibrational temperatures. The thermal non-equilibrium model is used
to better simulate the aerodynamic and thermodynamic phenomenon. Current results
were compared with the experimental data, computational fluid dynamics (CFD) results, and the Direct Simulation Monte Carlo (DSMC) results. Other CFD methods
include the original GKS method without considering thermal non-equilibrium, the
GKS method with a two temperature kinetic model and the Navier-Stokes equations
with a three temperature kinetic model, which is the same as the multiple temperature kinetic model in current GKS method. Comparisons were made for the surface
heat flux, pressure loads, Mach number contours and flow field values, including rotational temperature and density. By Comparing with other CFD method, the current
GKS method showed a lot of improvement in predicting the magnitude and position
of heat flux peak on the surface. This demonstrated the good potential of the current GKS method in solving thermodynamic non-equilibrium problems in hypersonic
flows. The good performance of predicting the heat flux could bring a lot of benefit for
the designing of the thermal protection system (TPS) for the hypersonic vehicles. By
comparing with the original GKS method and the two temperature kinetic model, the
three temperature kinetic model showed its importance and accuracy in this case. 相似文献
20.
Ruo Li Tiao Lu Yanli Wang & Wenqi Yao 《Communications In Computational Physics》2014,15(3):569-595
A globally hyperbolic moment system up to arbitrary order for the Wigner
equation was derived in [6]. For numerically solving the high order hyperbolic moment system therein, we in this paper develop a preliminary numerical method for
this system following the NR$xx$ method recently proposed in [8], to validate the moment system of the Wigner equation. The developedmethod can keep both mass and
momentum conserved, and the variation of the total energy under control though it
is not strictly conservative. We systematically study the numerical convergence of the
solution to the moment system both in the size of spatial mesh and in the order of
the moment expansion, and the convergence of the numerical solution of the moment
system to the numerical solution of the Wigner equation using the discrete velocity
method. The numerical results indicate that the high order moment system in [6] is
a valid model for the Wigner equation, and the proposed numerical method for the
moment system is quite promising to carry out the simulation of the Wigner equation. 相似文献