No Arabic abstract
We study the acceleration of steady-state computation for microflow, which is modeled by the high-order moment models derived recently from the steady-state Boltzmann equation with BGK-type collision term. By using the lower-order model correction, a novel nonlinear multi-level moment solver is developed. Numerical examples verify that the resulting solver improves the convergence significantly thus is able to accelerate the steady-state computation greatly. The behavior of the solver is also numerically investigated. It is shown that the convergence rate increases, indicating the solver would be more efficient, as the total levels increases. Three order reduction strategies of the solver are considered. Numerical results show that the most efficient order reduction strategy would be $m_{l-1} = lceil m_{l} / 2 rceil$.
In [Z. Hu, R. Li, and Z. Qiao. Acceleration for microflow simulations of high-order moment models by using lower-order model correction. J. Comput. Phys., 327:225-244, 2016], it has been successfully demonstrated that using lower-order moment model correction is a promising idea to accelerate the steady-state computation of high-order moment models of the Boltzmann equation. To develop the existing solver, the following aspects are studied in this paper. First, the finite volume method with linear reconstruction is employed for high-resolution spatial discretization so that the degrees of freedom in spatial space could be reduced remarkably without loss of accuracy. Second, by introducing an appropriate parameter $tau$ in the correction step, it is found that the performance of the solver can be improved significantly, i.e., more levels would be involved in the solver, which further accelerates the convergence of the method. Third, Heuns method is employed as the smoother in each level to enhance the robustness of the solver. Numerical experiments in microflows are carried out to demonstrate the efficiency and to investigate the behavior of the new solver. In addition, several order reduction strategies for the choice of the order sequence of the solver are tested, and the strategy $m_{l-1} = lceil m_{l} / 2 rceil$ is found to be most efficient.
We develop proper correction formulas at the starting $k-1$ steps to restore the desired $k^{rm th}$-order convergence rate of the $k$-step BDF convolution quadrature for discretizing evolution equations involving a fractional-order derivative in time. The desired $k^{rm th}$-order convergence rate can be achieved even if the source term is not compatible with the initial data, which is allowed to be nonsmooth. We provide complete error estimates for the subdiffusion case $alphain (0,1)$, and sketch the proof for the diffusion-wave case $alphain(1,2)$. Extensive numerical examples are provided to illustrate the effectiveness of the proposed scheme.
We present an arbitrarily high-order, conditionally stable, partitioned spectral deferred correction (SDC) method for solving multiphysics problems using a sequence of pre-existing single-physics solvers. This method extends the work in [1, 2], which used implicit-explicit Runge-Kutta methods (IMEX) to build high-order, partitioned multiphysics solvers. We consider a generic multiphysics problem modeled as a system of coupled ordinary differential equations (ODEs), coupled through coupling terms that can depend on the state of each subsystem; therefore the method applies to both a semi-discretized system of partial differential equations (PDEs) or problems naturally modeled as coupled systems of ODEs. The sufficient conditions to build arbitrarily high-order partitioned SDC schemes are derived. Based on these conditions, various of partitioned SDC schemes are designed. The stability of the first-order partitioned SDC scheme is analyzed in detail on a coupled, linear model problem. We show that the scheme is conditionally stable, and under conditions on the coupling strength, the scheme can be unconditionally stable. We demonstrate the performance of the proposed partitioned solvers on several classes of multiphysics problems including a simple linear system of ODEs, advection-diffusion-reaction systems, and fluid-structure interaction problems with both incompressible and compressible flows, where we verify the design order of the SDC schemes and study various stability properties. We also directly compare the accuracy, stability, and cost of the proposed partitioned SDC solver with the partitioned IMEX method in [1, 2] on this suite of test problems. The results suggest that the high-order partitioned SDC solvers are more robust than the partitioned IMEX solvers for the numerical examples considered in this work, while the IMEX methods require fewer implicit solves.
We present a paradigm for developing arbitrarily high order, linear, unconditionally energy stable numerical algorithms for gradient flow models. We apply the energy quadratization (EQ) technique to reformulate the general gradient flow model into an equivalent gradient flow model with a quadratic free energy and a modified mobility. Given solutions up to $t_n=n Delta t$ with $Delta t$ the time step size, we linearize the EQ-reformulated gradient flow model in $(t_n, t_{n+1}]$ by extrapolation. Then we employ an algebraically stable Runge-Kutta method to discretize the linearized model in $(t_n, t_{n+1}]$. Then we use the Fourier pseudo-spectral method for the spatial discretization to match the order of accuracy in time. The resulting fully discrete scheme is linear, unconditionally energy stable, uniquely solvable, and may reach arbitrarily high order. Furthermore, we present a family of linear schemes based on prediction-correction methods to complement the new linear schemes. Some benchmark numerical examples are given to demonstrate the accuracy and efficiency of the schemes.
This paper extends the model reduction method by the operator projection to the three-dimensional special relativistic Boltzmann equation. The derivation of arbitrary order moment system is built on our careful study of infinite families of the complicate Grad type orthogonal polynomials depending on a parameter and the real spherical harmonics. We derive the recurrence relations of the polynomials, calculate their derivatives with respect to the independent variable and parameter respectively, and study their zeros. The recurrence relations and partial derivatives of the real spherical harmonics are also given. It is proved that our moment system is globally hyperbolic, and linearly stable. Moreover, the Lorentz covariance is also studied in the 1D space.