Lagrange柱坐标高阶中心型守恒格式

提出Lagrange柱坐标高阶中心型守恒格式.基于用对守恒律的单调迎风算法(MUSCL)构造的高阶子网格压力,引入了柱坐标高阶体权子网格力和柱坐标高阶面权子网格力,构造了柱坐标高阶体权中心型守恒格式和柱坐标高阶面权中心型格式.柱坐标高阶体权中心型守恒格式满足动量守恒、能量守恒,但不能确定保持一维球对称性.柱坐标高阶面权中心型格式满足能量守恒,保持一维球对称性.两种格式里,格点速度以与网格面的数值通量相容的方式计算.对Saltzman活塞问题等进行了数值模拟,数值结果显示Lagrange柱坐标高阶中心型守恒格式的有效性和精确性.     

薛定谔方程在非一致网格下数值模拟的界面条件

本文研究了一维线性薛定谔方程在非一致网格下数值模拟的问题.在数值模拟中,非一致网格在界面处会产生虚假反射,利用局部时间步长和界面条件的方法,成功的减小了虚假反射.改进和提高了薛定谔方程数值模拟的效率和精度.     

薛定谔方程在非一致网格下数值模拟的界面条件（英文）

本文研究了一维线性薛定谔方程在非一致网格下数值模拟的问题.在数值模拟中,非一致网格在界面处会产生虚假反射,利用局部时间步长和界面条件的方法,成功的减小了虚假反射.改进和提高了薛定谔方程数值模拟的效率和精度.     

不可压饱和多孔弹性杆动力响应的多辛方法

研究了不可压饱和多孔弹性杆的一维动力响应问题.基于多孔介质理论,在流相和固相微观不可压、固相骨架小变形的假定下,建立了不可压流体饱和多孔弹性杆一维轴向动力响应的数学模型.利用Hamilton空间体系的多辛理论,构造了不可压饱和多孔弹性杆轴向振动方程的多辛形式及其多种局部守恒律.采用中点Box离散方法得到轴向振动方程的多辛离散格式和局部能量守恒律以及局部动量守恒律的离散格式;数值模拟了不可压饱和多孔弹性杆的轴向振动过程,记录了每一时间步的局部能量数值误差和局部动量数值误差.结果表明,已构造的多辛离散格式具有很高的精确性和较长时间的数值稳定性,这为解决饱和多孔介质的动力响应问题提供了新的途径.     

基于Rosenbrock型指数积分的一维间断Galerkin有限元方法

提出基于Rosenbrock型指数积分的一维间断Galerkin有限元方法.该方法在空间上使用间断有限元方法离散,在时间上采用Rosenbrock型指数积分方法.这样不仅可以保持空间离散上的高精度,而且继承了指数时间积分方法具有显式大步长时间推进的优点.数值试验的结果表明,对于一维双曲守恒律问题,这种方法是一种有效的数值算法.     

土壤水流与溶质耦合运移问题的混合元-迎风广义差分法研究

本文研究了一维非饱和土壤水流与溶质耦合运移问题的数学模型, 建立了求其数值解的守恒混合元-迎风广义差分格式. 对非线性土壤水分入渗方程, 采用守恒混合元法进行离散模拟, 同时得到了土壤含水量和水分通量; 而对对流-扩散形式的溶质运移方程, 利用迎风的广义差分法离散求解. 且分析了解的存在唯一性, 并讨论了误差估计. 最后给出数值算例, 模拟结果表明利用本文格式来求解非饱和土壤水流与溶质耦合运移问题是可靠的, 且该格式具有稳定性和可实用性.     

不可压缩黏性流体的二维Navier-Stokes方程的间断有限元模拟

由于不可压缩Navier-Stokes方程由守恒律、扩散及约束发展方程混合构成,为测试数值方法,该文基于非结构网格,对该方程建立了DG(discontinuous Galerkin)格式,讨论了不同黏性系数ν在方腔涡流问题的数值结果,验证了该方法的有效性且不依赖于问题的维数.圆柱绕流问题的模拟结果进一步表明此方法精度高、可有效求解具有运动界面的不可压缩黏性流体问题,使得模拟边界层、剪切层及复杂涡流解十分有效,并可以成功地推广到解决复杂现象数值模拟中的激波结构.     

二维双曲守恒律的大时间步长Godunov方法(英文)

考虑了关于二维守恒律的大时间步长Godunov方法.该方法是关于一维问题的自然推广.证明了文中给出的数值流函数下,该方法是守恒的.进一步还给出了近似Riemann解应满足的条件,并且证明了利用满足这些条件的近似Riemann解的大时间步长Godunov方法守恒.最后,补充证明了满足这些条件的近似Riemann解是满足熵条件的.     

逆时热传导问题的伴随同化研究

利用伴随同化方法对一维逆时热传导问题进行了数值研究,这里在目标函数中引入了稳定项.数值模拟结果表明,利用该方法反演一维逆时热问题的初值是可行的.     

守恒混合元法在模拟污染物运移过程中的应用

研究污染物在土壤中运移的时空规律,为土壤环境质量评价及污染预测和防治提供科学的根据与途径,具有重要的理论和实际意义.通过建立土壤中污染物运移问题的全离散守恒混合元格式,讨论了守恒混合元解的存在唯一性,并给出了误差估计.最后给出了数值算例,数值模拟结果表明,用该方法模拟污染物运移问题是合理有效的.     

High-order linear compact conservative method for the nonlinear Schrödinger equation coupled with the nonlinear Klein–Gordon equation

In this paper, we design a linear-compact conservative numerical scheme which preserves the original conservative properties to solve the Klein–Gordon–Schrödinger equation. The proposed scheme is based on using the finite difference method. The scheme is three-level and linear-implicit. Priori estimate and the convergence of the finite difference approximate solutions are discussed by the discrete energy method. Numerical results demonstrate that the present scheme is conservative, efficient and of high accuracy.     

Numerical conservative solutions of the Hunter–Saxton equation

BIT Numerical Mathematics - In the article a convergent numerical method for conservative solutions of the Hunter–Saxton equation is derived. The method is based on piecewise linear...     

Splitting over physical processes as applied to the construction of a numerical method for solving the system of kinetic equations governing a hall thruster rarefied plasma jet

A numerical method is constructed for solving the system of kinetic equations describing the behavior of a rarefied plasma jet exhausted from a Hall thruster. A similar problem was previously considered in the steady case. Now the same problem is solved in the unsteady formulation. The numerical method is based on a generalization of the splitting method with respect to physical processes, which is widely used in rarefied gas dynamics. The basic difficulty faced in the natural generalization of this method as applied to thruster jets is that a boundary condition has to be taken into account when the ion distribution function is determined. This difficulty is overcome by analytically selecting the corresponding term at the stage of free-molecular motion. Another difficulty that can be coped with by the splitting method is that the ion and neutral velocity spaces have widely different scales. Techniques for making the method conservative are described, and situations in which this is necessary are discussed. Qualitative characteristics of some numerical computations are compared with experimental data. The comparison shows that the numerical method constructed adequately reproduces the behavior of the modeled object.     

A conservative parallel difference method for 2-dimension diffusion equation

In this paper, a conservative parallel difference scheme, which is based on domain decomposition method, for 2-dimension diffusion equation is proposed. In the construction of this scheme, we use the numerical solution on the previous time step to give a weighted approximation of the numerical flux. Then the sub-problems with Neumann boundary are computed by fully implicit scheme. What is more, only local message communication is needed in the program. We use the method of discrete functional analysis to give the proof of the unconditional stability and second-order convergence accuracy. Some numerical tests are given to verify the theory results.     

A conservative local discontinuous Galerkin method for the solution of nonlinear Schrdinger equation in two dimensions

In this study, we present a conservative local discontinuous Galerkin(LDG) method for numerically solving the two-dimensional nonlinear Schrdinger(NLS) equation. The NLS equation is rewritten as a firstorder system and then we construct the LDG formulation with appropriate numerical flux. The mass and energy conserving laws for the semi-discrete formulation can be proved based on different choices of numerical fluxes such as the central, alternative and upwind-based flux. We will propose two kinds of time discretization methods for the semi-discrete formulation. One is based on Crank-Nicolson method and can be proved to preserve the discrete mass and energy conservation. The other one is Krylov implicit integration factor(IIF) method which demands much less computational effort. Various numerical experiments are presented to demonstrate the conservation law of mass and energy, the optimal rates of convergence, and the blow-up phenomenon.     

A new TVD flux-limiter method for solving nonlinear hyperbolic equations

We develop a new flux-limiter method based on the Richtmyer two-step Lax-Wendroff (R2LW) method coupled with a conservative upwind method and a nonconventional flux-limiter function. The proposed method is TVD stable and preserves the linear stability condition of the R2LW method. The numerical results show that this method improves on the approximations obtained by standard first-order and second-order TVD schemes. The new method provides accurate results for nonlinear hyperbolic equations with discontinuous solutions.     

On the convergence of a conservative numerical scheme for the usual Rosenau-RLW equation

In this paper, we study the initial-boundary value problem of the usual Rosenau-RLW equation by finite difference method. We design a conservative numerical scheme which preserves the original conservative properties for the equation. The scheme is three-level and linear-implicit. The unique solvability of numerical solutions has been shown. Priori estimate and second order convergence of the finite difference approximate solutions are discussed by discrete energy method. Numerical results demonstrate that the scheme is efficient and accurate.     

Conservative schemes of matter transport in a system of vessels closed by the heart

We suggest a conservative model and the corresponding conservative numerical method for analyzing the motion of matter in a branched vessel system closed by the heart. The efficiency of the algorithm is justified by test computations.     

Error Estimate of a New Conservative Finite Difference Scheme for the Klein-Gordon-Dirac System

In this paper, we derive and analyze a conservative Crank-Nicolson-type finite difference scheme for the Klein-Gordon-Dirac (KGD) system. Differing from the derivation of the existing numerical methods given in literature where the numerical schemes are proposed by directly discretizing the KGD system, we translate the KGD equations into an equivalent system by introducing an auxiliary function, then derive a nonlinear Crank-Nicolson-type finite difference scheme for solving the equivalent system. The scheme perfectly inherits the mass and energy conservative properties possessed by the KGD, while the energy preserved by the existing conservative numerical schemes expressed by two-level's solution at each time step. By using energy method together with the 'cut-off' function technique, we establish the optimal error estimate of the numerical solution, and the convergence rate is $\mathcal{O}(τ^2 + h^2)$ in $l^∞$-norm with time step $τ$ and mesh size $h.$ Numerical experiments are carried out to support our theoretical conclusions.     

Conservative numerical method for solving the averaged Boltzmann equation

A method is proposed for averaging the Boltzmann kinetic equation with respect to transverse velocities. A system of two integro-differential equations for two desired functions depending only on the longitudinal velocity is derived. The gas particles are assumed to interact as absolutely hard spheres. The integrals in the equations are double. The reduction in the number of variables in the desired functions and the low multiplicity of the integrals ensure the computational efficiency of the averaged equations. A numerical method of discrete ordinates is developed that effectively solves the gas relaxation problem based on the averaged equations. The method is conservative, and the number of particles, momentum, and energy are conserved automatically at every time step.