Vlasov-Poisson-Fokker-Planck 方程组 到不可压 Euler方程的收敛性

研究Vlasov-Poisson-Fokker-Planck (VPFP)方程组的拟中性极限问题. 通过使用紧致性理论和相对熵方法证明了VPFP方程组到不可压Euler方程的收敛性, 这个结果在不可压Euler方程光滑解存在的时间区间上都成立.     

二维三温热传导方程的并行自适应多重网格算法求解

3-T heat conduct equation including electron, ion and photon (radiation) temperatures can be used to approximately describe the energy broadcast across multimedia for radial flow dynamics and discover the energy swapping among photon,electron and ion. Owing to the strong nonlinear diffusion coefficients and energy swapping coefficients and strong discontinuous coefficients across media interfaces,this equation is difficult to be solved with high numerical resolution. Based on the parallel adaptive multigrid software framework UG on 2-D unstructured grid, this paper successfully solved such equation with high resolution by combining the finite volume implicit discretization scheme and parallel adaptive multigrid algorithm, and gained much significant results.     

基于非结构自适应网格的二维Euler方程数值求解方法研究

提出了一种基于非结构自适应网格的二维Euler方程的数值解法.采用有限体积法进行空间离散,通量计算采用Jamson中心格式,使得它适用于任意多边形计算单元.为了得到定常解,采用一种显式的四步Runge-Kutta迭代方法对时间进行积分.根据流场参数的变化梯度确定加密边,由加密准则进行自适应网格剖分,然后得到分布合理的加密过后的网格.求解二维Euler方程,对NACA0012翼型进行了数值模拟,通过对自适应前后的数值解的对比,说明所建立的方法是正确的.     

不可压缩Navier-Stokes方程最优动力系统建模和分析

研究了同时满足任意速度边界条件和速度不可压条件的Navier-Stokes方程最优动力系统的建模方法.通过对方柱绕流问题的最优动力系统的建模与分析,发现该最优动力系统的动力学特性为极限环.同时,该最优动力系统仅使用了三个最优基函数就很好地描述了所有主要的流场特征和该问题的动力学特性,故满足任意速度边界条件和速度不可压条件Navier-Stokes方程最优动力系统的建模方法,能够用最少的基函数最大限度地描述复杂流体问题及其动力学特性.     

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

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

二维奇异摄动对流扩散方程的自适应移动网格算法

针对二维奇异摄动对流扩散方程,在任意网格下给出了经典的迎风有限差分格式.利用二元多项式插值技术,推导出一阶最大范数的后验误差估计,并以此设计了一个自适应网格生成算法.数值实验表明本文构造的自适应移动网格算法是有效的.     

带插值的自适应网格重构算法求解带移动热源的反应扩散方程的理论分析

本文研究一个带插值的网格重构算法求解一类带移动热源的反应扩散方程. 算法包括两步: 第一步是用旧时间网层上的计算解计算新时间层上的空间网格; 第二步是使用有限差分方法在新时间层 空间网格上离散方程, 并且将旧时间层上计算解的插值作为初始值. 对于时间, 我们获得了一阶收敛结果. 对于空间, 我们证明了使用线性插值算法的一阶收敛性和使用二次插值算法的二阶收敛性. 数值例子肯定了本文的理论结果.     

奇异摄动反应扩散方程的后验误差估计及自适应算法

研究了一类奇异摄动半线性反应扩散方程的自适应网格方法.在任意非均匀网格上建立迎风有限差分离散格式,并推导出离散格式的后验误差界,然后用该误差界设计自适应网格移动算法.数值实验结果证明了所提出的自适应网格方法的有效性.     

基于二重网格的定常Navier-Stokes方程的局部和并行有限元算法

对二维定常的不可压缩的Navier-Stokes方程的局部和并行算法进行了研究．给出的算法是多重网格和区域分解相结合的算法,它是基于两个有限元空间:粗网格上的函数空间和子区域的细网格上的函数空间．局部算法是在粗网格上求一个非线性问题,然后在细网格上求一个线性问题,并舍掉内部边界附近的误差相对较大的解．最后,基于局部算法,通过有重叠的区域分解而构造了并行算法,并且做了算法的误差分析,得到了比标准有限元方法更好的误差估计,也对算法做了数值试验,数值结果通过比较验证了本算法的高效性和合理性．     

非定常不可压Navier-Stokes方程两重网格方法的稳定性和L2误差估计

讨论了二维非定常不可压Navier-Stokes方程的两重网格方法．此方法包括在粗网格上求解一个非线性问题,在细网格上求解一个Stokes问题．采用一种新的全离散(时间离散用Crank-Nicolson格式,空间离散用混合有限元方法)格式数值求解N-S方程．证明了该全离散格式的稳定性．给出了L2误差估计．对比标准有限元方法,在保持同样精度的前提下,TGM能节省大量的计算量．     

曲线网格下基于粘声分离方法的流噪声计算

在曲线网格下基于粘声分离方法对流场中的静止圆柱同声波和涡波的相互作用进行研究.首先推导了曲线坐标系下、适用于水流噪声的粘声分离方法(viscous acoustic splitting method,VASM)控制方程,并采用7点色散关系保持(dispersion-relation preserving,DRP)格式和四阶时间差分格式进行计算.然后将静止流场中圆柱壁面对声波反射的计算结果同理论值进行比较,验证了计算方法模拟水中物体对声波散射的准确性.进而模拟了旋涡行走发声的特性,并分析了流速等对声场特性的影响.     

Discontinuous Stable Elements for the Incompressible Flow

In this paper, we derive a discontinuous Galerkin finite element formulation for the Stokes equations and a group of stable elements associated with the formulation. We prove that these elements satisfy the new inf–sup condition and can be used to solve incompressible flow problems. Associated with these stable elements, optimal error estimates for the approximation of both velocity and pressure in L ² norm are obtained for the Stokes problems, as well as an optimal error estimate for the approximation of velocity in a mesh dependent norm.     

Optimal Boundary Control of Steady-State Flow of a Viscous Inhomogeneous Incompressible Fluid

We study the problem of optimal boundary control of two-dimensional steady-state flow of a viscous inhomogeneous incompressible fluid. The role of control is played by the values of the velocity on a part of the boundary of the domain considered. On the remaining part of the boundary, the vector of flow velocity and the fluid density are given. We seek the fluid density as a scalar function (determined by the initial data) of the stream function, study the solvability of the problem, and obtain necessary optimality conditions.     

Preconditioners for Incompressible Navier-Stokes Solvers

<正>In this paper we give an overview of the present state of fast solvers for the solution of the incompressible Navier-Stokes equations discretized by the finite element method and linearized by Newton or Picard's method.It is shown that block preconditioners form an excellent approach for the solution,however if the grids are not to fine preconditioning with a Saddle point ILU matrix(SILU) may be an attractive alternative. The applicability of all methods to stabilized elements is investigated.In case of the stand-alone Stokes equations special preconditioners increase the efficiency considerably.     

瞬态两相流的矢通量分裂解

本文用矢通量分裂法求解了瞬态两相流问题.首先导出了瞬态两相流控制方程的特征根,然后将矢通量分裂成二个子矢量,使每一个子矢量仅包含一种符号的特征值,再利用空间的一侧差分,得到差分方程并求解,最后将数值预测与实验结果进行了比较,其结果十分吻合.     

Adaptive Runge-Kutta Discontinuous Galerkin Methods with Modified Ghost Fluid Method for Simulating the Compressible Two-Medium Flow

In this paper, we investigate using the adaptive Runge-Kutta discontinuous Galerkin (RKDG) methods with the modified ghost fluid method (MGFM) in conjunction with the adaptive RKDG methods for solving the level set function to simulate the compressible two-medium flow in one and two dimensions. A shock detection technique (KXRCF method) is adopted as an indicator to identify the troubled cell, which serves for further numerical limiting procedure which uses a modified TVB limiter to reconstruct different degrees of freedom and an adaptive mesh refinement procedure. If the computational mesh should be refined or coarsened, and the detail of the implementation algorithm is presented on how to modulate the hanging nodes and redefine the numerical solutions of the two-medium flow and the level set function on such adaptive mesh. Extensive numerical tests are provided to illustrate the proposed adaptive methods may possess the capability of enhancing the resolutions nearby the discontinuities inside of the single medium flow region and material interfacial vicinities of the two-medium flow region.     

裂缝孔隙介质中二相驱动问题的交替方向有限元方法及理论分析

考虑裂缝孔隙介质中二相驱动问题的数值方法及理论分析。对压力方程采用混合有限元方法,对裂缝和岩块系统上的饱和度方程采用交替方向有限元方法,证明了交替方向有限元格式具有最优Ｌ２模和Ｈ１模误差估计。     

一种h型自适应有限单元

ｈ型自适应有限单元在网格局部细划时．会产生非常规节点，从而破坏了一般意义上的单元连续性假定．本文利用参照节点对非常规单元进行坐标和位移插值．为保证单元之间坐标和位移的连续性，本文提出了一组修正的形函数，常用的形函数是它的一个特例．本方法应用于有限元程序时，除形函数外无须做任何改动．算例表明水文的方法具有方法简单、精度高、自由度少、计算量小等优点．     

Boundary Asymptotic Analysis for an Incompressible Viscous Flow: Navier Wall Laws

We consider a new way of establishing Navier wall laws. Considering a bounded domain Ω of R ^N , N=2,3, surrounded by a thin layer Σ _ε , along a part Γ₂ of its boundary ∂Ω, we consider a Navier-Stokes flow in Ω∪∂Ω∪Σ _ε with Reynolds’ number of order 1/ε in Σ _ε . Using Γ-convergence arguments, we describe the asymptotic behaviour of the solution of this problem and get a general Navier law involving a matrix of Borel measures having the same support contained in the interface Γ₂. We then consider two special cases where we characterize this matrix of measures. As a further application, we consider an optimal control problem within this context.     

A Two‐Parameter Stabilized Finite Element Method for Incompressible Flows

Based on two‐grid discretizations, a two‐parameter stabilized finite element method for the steady incompressible Navier–Stokes equations at high Reynolds numbers is presented and studied. In this method, a stabilized Navier–Stokes problem is first solved on a coarse grid, and then a correction is calculated on a fine grid by solving a stabilized linear problem. The stabilization term for the nonlinear Navier–Stokes equations on the coarse grid is based on an elliptic projection, which projects higher‐order finite element interpolants of the velocity into a lower‐order finite element interpolation space. For the linear problem on the fine grid, either the same stabilization approach (with a different stabilization parameter) as that for the coarse grid problem or a completely different stabilization approach could be employed. Error bounds for the discrete solutions are estimated. Algorithmic parameter scalings of the method are also derived. The theoretical results show that, with suitable scalings of the algorithmic parameters, this method can yield an optimal convergence rate. Numerical results are provided to verify the theoretical predictions and demonstrate the effectiveness of the proposed method. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 425–444, 2017