求解Oseen流的交替线松弛多重网格方法

利用Riemann解的通量差分分裂法——Godunov方法对Oseen流控制方程进行离散,得到了基于一阶上迎风格式的离散方程,并给出了使用多重网格方法求解该离散方程的V-循环算法和W-循环算法的收敛性分析.通过局部Fourier分析方法,对获得的离散方程的聚对称交替线GaussSeidel松弛的光滑性质进行了研究.结果表明:使用多重网格的两层网格及三层网格算法求解具有不同Reynolds数的Oseen流,即便是在高Reynolds数情况下,聚对称交替线Gauss-Seidel松弛具有很好的光滑性质,多重网格W-循环算法收敛性比V-循环算法好.     

带初始奇异性的多项时间分数阶扩散方程一类全离散数值方法的误差分析

本文研究了带有初始奇异性的多项时间分数阶扩散方程的一种全离散数值方法.首先,基于L1公式在渐变网格下离散多项Caputo时间分数阶导数,构造了多项时间分数阶扩散方程的时间半离散格式,证明了时间格式通过选取合适的网格参数r,时间方向的误差可以达到最优的收敛阶2-α_1,其中α_1(0 α_11)为多项时间分数阶导数阶数的最大值.然后,空间采用谱方法进行离散,得到了全离散格式,证明了全离散格式的无条件稳定性和收敛性.为了降低计算量和储存量,对多项时间分数阶扩散方程又构造了时间方向的快速算法,同时证明了该格式的收敛性.数值算例验证了算法的有效性,显示了快速算法的高效性.     

定常Navier-Stokes方程的一种高效稳定有限元方法

提出了二维定常Navier-Stokes(N-S)方程的一种两层稳定有限元方法.该方法基于局部高斯积分技术,通过不满足inf-sup条件的低次等阶有限元对N-S方程进行有限元求解.该方法在粗网格上解定常N-S方程,在细网格上只需解一个Stokes方程.误差分析和数值试验都表明:两层稳定有限元方法与直接在细网格上采用的传统有限元方法得到的解具有同阶的收敛性,但两层稳定有限元方法节省了大量的工作时间.     

带移动奇异源热方程的移动网格方法

研究了一个带若干奇异源热方程的数值求解,其源的移动由一个常微分方程描述.基于移动观察区域和区域分解思想提出了一个移动网格预估校正算法.网格方程可自然的通过并行高效求解,算法避免了跳跃信息[u]的计算而使物理方程的离散格式变得非常简单,且仍保持了空间上的二阶收敛性.数值例子验证了算法的收敛性和高效性,并模拟了非线性源函数带来的爆破现象.     

带漂移的单侧正规化回火分数阶扩散方程的Crank-Nicolson拟紧格式

基于已有的针对单侧正规化回火分数阶扩散方程的三阶拟紧算法,将该算法的思想应用于带漂移的单侧正规化回火分数阶扩散方程的数值模拟,并结合Crank-Nicolson方法导出数值格式.证明了数值格式的稳定性与收敛性,且数值格式的时间收敛阶和空间收敛阶分别是二阶和三阶.通过数值试验验证了数值格式的有效性和理论结果.     

关于热传导方程半离散差分格式的一个注记

本文研究了热传导方程初边值问题的半离散化差分格式直接解算法.分别从Dirichlet和Neumann边界条件出发,直接由空间差分格式导出与时间相关的一阶常微分方程组,随后通过正/余弦变换获得了原方程的半解析解,并给出了相关收敛性分析.并对中心差分格式和紧差分格式的精度差异,通过矩阵特征值理论给出了相关原因分析.另外,对于二维热传导方程初边值问题,应用矩阵张量积运算,该直接解算法可直接演变成二重正(余)弦变换.该方法由于不涉及时间上的离散,从而具有较好的计算效率.     

基于抛物线插值的Allen-Cahn方程时间两重网格有限元算法

本文运用抛物线插值的时间两重网格(TT-M)有限元(FE)算法求解非线性Allen-Cahn方程.首先,对非线性Allen-Cahn方程,在空间和时间上分别采用有限元方法以及二阶θ格式进行离散.其次,使用抛物线插值的时间两重网格有限元算法求解Allen-Cahn方程.同时,在粗细时间步长上,对数值解进行了稳定性分析和误差估计.最后,通过数值实验验证方法的有效性.     

求解定常不可压Stokes方程的两层罚函数方法

借助于两套有限元网格空间提出了一种求解定常不可压Stokes方程的两层罚函数方法.该方法只需要求解粗网格空间上的Stokes方程和细网格空间上的两个易于求解的罚参数方程（离散后的线性方程组具有相同的对称正定系数矩阵）.收敛性分析表明粗网格空间相对于细网格空间可以选择很小，并且罚参数的选取只与粗网格步长和问题的正则性有关.因此罚参数不必选择很小仍能够得到最优解.最后通过数值算例验证了上述理论结果，并且数值对比可知两层罚函数方法对于求解定常不可压Stokes方程具有很好的效果.     

带非线性源项的双侧空间分数阶扩散方程的隐式中点方法

本文用隐式中点方法离散一阶时间偏导数,并用拟紧差分算子逼近Riemann-Liouville空间分数阶偏导数,构造了求解带非线性源项的空间分数阶扩散方程的数值格式.给出了数值方法的稳定性和收敛性分析.数值试验表明数值方法是有效的.     

时间导数项含小参数的抛物方程的数值解法

本文讨论时间导数项含小参数的抛物方程.我们依Бахволов构造非均匀网格的差分格式,并证明了格式的一阶一致收敛性.给出了数值结果.     

A space‐time continuous Galerkin method with mesh modification for viscoelastic wave equations

In this article, we consider the space‐time continuous Galerkin (STCG) method for the viscoelastic wave equations. It allows variable temporal step‐sizes, and the changing of the spatial grids in two adjacent time levels. The existence, uniqueness, and stability of the approximate solutions are demonstrated and the error estimates with global and local spatial mesh sizes in norm are derived without any restrictive assumptions on the space‐time meshes. If the meshes in each time level satisfy some reasonable assumptions, then we can get the optimal order error estimates both in time and space. Finally, we give a numerical example on unstructured meshes to confirm the theoretical findings. © 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1183–1207, 2017     

An adaptive finite element algorithm with reliable and efficient error control for linear parabolic problems

An efficient and reliable a posteriori error estimate is derived for linear parabolic equations which does not depend on any regularity assumption on the underlying elliptic operator. An adaptive algorithm with variable time-step sizes and space meshes is proposed and studied which, at each time step, delays the mesh coarsening until the final iteration of the adaptive procedure, allowing only mesh and time-step size refinements before. It is proved that at each time step the adaptive algorithm is able to reduce the error indicators (and thus the error) below any given tolerance within a finite number of iteration steps. The key ingredient in the analysis is a new coarsening strategy. Numerical results are presented to show the competitive behavior of the proposed adaptive algorithm.

     

A Multigrid Block LU-SGS Algorithm for Euler Equations on Unstructured Grids

We propose an efficient and robust algorithm to solve the steady Euler equa- tions on unstructured grids.The new algorithm is a Newton-iteration method in which each iteration step is a linear multigrid method using block lower-upper symmetric Gauss-Seidel(LU-SGS)iteration as its smoother To regularize the Jacobian matrix of Newton-iteration,we adopted a local residual dependent regularization as the replace- ment of the standard time-stepping relaxation technique based on the local CFL number The proposed method can be extended to high order approximations and three spatial dimensions in a nature way.The solver was tested on a sequence of benchmark prob- lems on both quasi-uniform and local adaptive meshes.The numerical results illustrated the efficiency and robustness of our algorithm.     

Decay of Dissipative Equations and Negative Sobolev Spaces

We develop a general energy method for proving the optimal time decay rates of the solutions to the dissipative equations in the whole space. Our method is applied to classical examples such as the heat equation, the compressible Navier-Stokes equations and the Boltzmann equation. In particular, the optimal decay rates of the higher-order spatial derivatives of solutions are obtained. The negative Sobolev norms are shown to be preserved along time evolution and enhance the decay rates. We use a family of scaled energy estimates with minimum derivative counts and interpolations among them without linear decay analysis.     

Application of Richardson extrapolation to the numerical solution of partial differential equations

Richardson extrapolation is a methodology for improving the order of accuracy of numerical solutions that involve the use of a discretization size h. By combining the results from numerical solutions using a sequence of related discretization sizes, the leading order error terms can be methodically removed, resulting in higher order accurate results. Richardson extrapolation is commonly used within the numerical approximation of partial differential equations to improve certain predictive quantities such as the drag or lift of an airfoil, once these quantities are calculated on a sequence of meshes, but it is not widely used to determine the numerical solution of partial differential equations. Within this article, Richardson extrapolation is applied directly to the solution algorithm used within existing numerical solvers of partial differential equations to increase the order of accuracy of the numerical result without referring to the details of the methodology or its implementation within the numerical code. Only the order of accuracy of the existing solver and certain interpolations required to pass information between the mesh levels are needed to improve the order of accuracy and the overall solution accuracy. Using the proposed methodology, Richardson extrapolation is used to increase the order of accuracy of numerical solutions of the linear heat and wave equations and of the nonlinear St. Venant equations in one‐dimension. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

Moving mesh method for problems with blow-up on unbounded domains

This paper studies the numerical solutions of semilinear parabolic partial differential equations (PDEs) on unbounded spatial domains whose solutions blow up in finite time. There are two major difficulties usually in numerical solutions: the singularity of blow-up and the unboundedness. We propose local absorbing boundary conditions (LABCs) on the selected artificial boundaries by using the idea of unified approach (Brunner et al., SIAM J Sci Comput 31:4478–4496, (2010). Since the uniform fixed spatial meshes may be inefficient, we adopt moving mesh partial differential equation (MMPDE) method to adapt the spatial mesh as the singularity develops. Combining LABCs and MMPDE, we can effectively capture the qualitative behavior of the blow-up singularities in the unbounded domain. Moreover, the implementation of the combination consists of two independent parts. Numerical examples also illustrate the efficiency and the accuracy of the new method.     

Time step selection for the numerical solution of boundary value problems for parabolic equations

An algorithm is proposed for selecting a time step for the numerical solution of boundary value problems for parabolic equations. The solution is found by applying unconditionally stable implicit schemes, while the time step is selected using the solution produced by an explicit scheme. Explicit computational formulas are based on truncation error estimation at a new time level. Numerical results for a model parabolic boundary value problem are presented, which demonstrate the performance of the time step selection algorithm.     

多层快速多极子方法的快速插值

多层快速多极子方法(MLFMM)可用来加速迭代求解由Maxwell方程组或Helmholtz方程导出的积分方程,其复杂度理论上是O(NlogN),N为未知量个数.MLFMM依赖于快速计算每层的转移项,以及上聚和下推过程中的层间插值.本文引入计算类似N体问题的一维快速多极子方法(FMM1D).基于FMM1D的快速Lagr...     

Mesh-dependent stability for finite element approximations of parabolic equations with mass lumping

In this paper, based on some mesh-dependent estimates on the extreme eigenvalues of a general finite element system defined on a simplicial mesh, novel and sharp bounds on the permissible time step size are derived for the mass lumping finite element approximations of parabolic equations. The bounds are dependent not only on the mesh size but also on the mesh shape. These results provide guidance to the stability of numerical solutions of parabolic problems in relation to the unstructured geometric meshing. Numerical experiments on both uniform meshes and adaptive meshes are presented to validate the theoretical analysis.     

Mesh-dependent stability for finite element approximations of parabolic equations with mass lumping

In this paper, based on some mesh-dependent estimates on the extreme eigenvalues of a general finite element system defined on a simplicial mesh, novel and sharp bounds on the permissible time step size are derived for the mass lumping finite element approximations of parabolic equations. The bounds are dependent not only on the mesh size but also on the mesh shape. These results provide guidance to the stability of numerical solutions of parabolic problems in relation to the unstructured geometric meshing. Numerical experiments on both uniform meshes and adaptive meshes are presented to validate the theoretical analysis.