Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA third-order discrete unified gas kinetic scheme for continuum and rarefied flows: low-speed isothermal case
291
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
An efficient third-order discrete unified gas kinetic scheme (DUGKS) with efficiency is presented in this work for simulating continuum and rarefied flows. By employing two-stage time-stepping scheme and the high-order DUGKS flux reconstruction strategy, third-order of accuracy in both time and space can be achieved in the present method. It is also analytically proven that the second-order DUGKS is a special case of the present method. Compared with the high-order lattice Boltzmann equation {LBE} based methods, the present method is capable to deal with the rarefied flows by adopting the Newton-Cotes quadrature to approximate the integrals of moments. Instead of being constrained by the second-order (or lower-order) of accuracy in time splitting scheme as in the conventional high-order Runge-Kutta (RK) based kinetic methods, the present method solves the original BE, which overcomes the limitation in time accuracy. Typical benchmark tests are carried out for comprehensive evaluation of the present method. It is observed in the tests that the present method is advantageous over the original DUGKS in accuracy and capturing delicate flow structures. Moreover, the efficiency of the present third-order method is also shown in simulating rarefied flows.
rate research
Read More
Fluid dynamic equations are valid in their respective modeling scales. With a variation of the modeling scales, theoretically there should have a continuous spectrum of fluid dynamic equations. In order to study multiscale flow evolution efficiently, the dynamics in the computational fluid has to be changed with the scales. A direct modeling of flow physics with a changeable scale may become an appropriate approach. The unified gas-kinetic scheme (UGKS) is a direct modeling method in the mesh size scale, and its underlying flow physics depends on the resolution of the cell size relative to the particle mean free path. The cell size of UGKS is not limited by the particle mean free path. With the variation of the ratio between the numerical cell size and local particle mean free path, the UGKS recovers the flow dynamics from the particle transport and collision in the kinetic scale to the wave propagation in the hydrodynamic scale. The previous UGKS is mostly constructed from the evolution solution of kinetic model equations. This work is about the further development of the UGKS with the implementation of the full Boltzmann collision term in the region where it is needed. The central ingredient of the UGKS is the coupled treatment of particle transport and collision in the flux evaluation across a cell interface, where a continuous flow dynamics from kinetic to hydrodynamic scales is modeled. The newly developed UGKS has the asymptotic preserving (AP) property of recovering the NS solutions in the continuum flow regime, and the full Boltzmann solution in the rarefied regime. In the mostly unexplored transition regime, the UGKS itself provides a valuable tool for the flow study in this regime. The mathematical properties of the scheme, such as stability, accuracy, and the asymptotic preserving, will be analyzed in this paper as well.
The recently proposed discrete unified gas kinetic scheme (DUGKS) is a finite volume method for deterministic solution of the Boltzmann model equation with asymptotic preserving property. In DUGKS, the numerical flux of the distribution function is determined from a local numerical solution of the Boltzmann model equation using an unsplitting approach. The time step and mesh resolution are not restricted by the molecular collision time and mean free path. To demonstrate the capacity of DUGKS in practical problems, this paper extends the DUGKS to arbitrary unstructured meshes. Several tests of both internal and external flows are performed, which include the cavity flow ranging from continuum to free molecular regimes, a multiscale flow between two connected cavities with a pressure ratio of 10000, and a high speed flow over a cylinder in slip and transitional regimes. The numerical results demonstrate the effectiveness of the DUGKS in simulating multiscale flow problems.
Unified gas kinetic scheme (UGKS) is an asymptotic preserving scheme for the kinetic equations. It is superior for transition flow simulations, and has been validated in the past years. However, compared to the well known discrete ordinate method (DOM) which is a classical numerical method solving the kinetic equations, the UGKS needs more computational resources. In this study, we propose a simplification of the unified gas kinetic scheme. It allows almost identical numerical cost as the DOM, but predicts numerical results as accurate as the UGKS. Based on the observation that the equilibrium part of the UGKS fluxes can be evaluated analytically, the equilibrium part in the UGKS flux is not necessary to be discretized in velocity space. In the simplified scheme, the numerical flux for the velocity distribution function and the numerical flux for the macroscopic conservative quantities are evaluated separately. The simplification is equivalent to a flux hybridization of the gas kinetic scheme for the Navier-Stokes (NS) equations and conventional discrete ordinate method. Several simplification strategies are tested, through which we can identify the key ingredient of the Navier-Stokes asymptotic preserving property. Numerical tests show that, as long as the collision effect is built into the macroscopic numerical flux, the numerical scheme is Navier-Stokes asymptotic preserving, regardless the accuracy of the microscopic numerical flux for the velocity distribution function.
In this paper, a high-order gas-kinetic scheme in general curvilinear coordinate (HGKS-cur) is developed for the numerical simulation of compressible turbulence. Based on the coordinate transformation, the Bhatnagar-Gross-Krook (BGK) equation is transformed from physical space to computational space. To deal with the general mesh given by discretized points, the geometrical metrics need to be constructed by the dimension-by-dimension Lagrangian interpolation. The multidimensional weighted essentially non-oscillatory (WENO) reconstruction is adopted in the computational domain for spatial accuracy, where the reconstructed variables are the cell averaged Jacobian and the Jacobian-weighted conservative variables. The two-stage fourth-order method, which was developed for spatial-temporal coupled flow solvers, is used for temporal discretization. The numerical examples for inviscid and laminar flows validate the accuracy and geometrical conservation law of HGKS-cur. As a direct application, HGKS-cur is implemented for the implicit large eddy simulation (iLES) in compressible wall-bounded turbulent flows, including the compressible turbulent channel flow and compressible turbulent flow over periodic hills. The iLES results with HGKS-cur are in good agreement with the refereed spectral methods and high-order finite volume methods. The performance of HGKS-cur demonstrates its capability as a powerful tool for the numerical simulation of compressible wall-bounded turbulent flows and massively separated flows.
The discrete unified gas kinetic scheme (DUGKS) is a new finite volume (FV) scheme for continuum and rarefied flows which combines the benefits of both Lattice Boltzmann Method (LBM) and unified gas kinetic scheme (UGKS). By reconstruction of gas distribution function using particle velocity characteristic line, flux contains more detailed information of fluid flow and more concrete physical nature. In this work, a simplified DUGKS is proposed with reconstruction stage on a whole time step instead of half time step in original DUGKS. Using temporal/spatial integral Boltzmann Bhatnagar-Gross-Krook (BGK) equation, the transformed distribution function with inclusion of collision effect is constructed. The macro and mesoscopic fluxes of the cell on next time step is predicted by reconstruction of transformed distribution function at interfaces along particle velocity characteristic lines. According to the conservation law, the macroscopic variables of the cell on next time step can be updated through its macroscopic flux. Equilibrium distribution function on next time step can also be updated. Gas distribution function is updated by FV scheme through its predicted mesoscopic flux in a time step. Compared with the original DUGKS, the computational process of the proposed method is more concise because of the omission of half time step flux calculation. Numerical time step is only limited by the Courant-Friedrichs-Lewy (CFL) condition and relatively good stability has been preserved. Several test cases, including the Couette flow, lid-driven cavity flow, laminar flows over a flat plate, a circular cylinder, and an airfoil, as well as micro cavity flow cases are conducted to validate present scheme. The numerical simulation results agree well with the references results.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
