求解浅水波方程的熵相容格式

提出了一种求解浅水波方程组的熵相容格式．在熵稳定通量中添加特征速度差分绝对值的项来抵消解在跨过激波时所产生的熵增,从而实现熵相容．新的数值差分格式具有形式简单、计算效率高、无需添加任何的人工数值粘性的特点．数值算例充分说明了其显著的优点．利用新格式成功地模拟了不同类型溃坝问题的激波、稀疏波传播及溃坝两侧旋涡的形成,是求解浅水波方程组较为理想的方法.     

求解二维Euler方程的旋转通量混合格式

为提高求解二维Euler方程数值结果的分辨率,提出了一种旋转通量混合格式.该算法采用旋转通量法的类一维处理思想,通量函数选用满足热力学第二定律的熵稳定数值通量和具有良好鲁棒性的HLL数值通量耦合的混合格式,时间方向采用三阶强稳定Runge-Kutta方法进行推进.该旋转通量混合格式具有结构简单、分辨率高的优点,数值结果表明了该算法的良好特性.     

数值求解Euler方程的UCGVC差分格式的注记

１．引言 近年来高精度差分格式的研究引起国内外的普遍重视,目的是更准确地模拟复杂流场的流动．众所周知,传统的二阶ＴＶＤ类格式虽然能较好地捕捉激波,但却存在局部极值点降阶的问题,而且由于一些格式的数值粘性过大,当用该格式计算粘性流特别是高雷诺数问题时,格式本身的数值粘性可能掩盖了流场的物理粘性,从而降低了格式对边界层的分辨率,因而无法正确计算热流值。文献［３］指出,采用高精度格式可适当放松对网格雷诺数的要求,因此发展三阶或三阶以上的格式是需要的。近年来,人们已经发展了一些无伪振荡的高阶格式,如ＥＮ…     

粘性机制下一类熵相容格式的对比研究

研究了一种人工和物理耗散机制下的离散熵相容格式,探讨数值粘性和物理粘性的大小以及它们所起的作用.所得结论是:在激波捕捉的过程中,粘性系数越大,则无需加入人工粘性项;粘性系数较小时,除了物理粘性项,还需要加入人工粘性项来得到熵相容格式.首先研究了一维粘性Burgers方程离散熵相容格式,再将其推广至Navier-Stokes方程.数值算例采用空间半离散格式,并结合显式三步三阶Runge-Kutta(RK3)方法进行时间推进.这两类方程的数值结果表明,最终选取的熵相容格式能够准确地捕捉到激波.     

求解理想磁流体方程的四阶WENO型熵稳定格式

构造了一种用于求解理想磁流体方程的四阶熵稳定半离散有限体积格式.该格式空间方向上将高阶熵守恒通量与采用WENO重构的耗散项结合,得到高阶熵稳定通量.通过在耗散项中添加开关函数,使得数值通量具有更低的耗散并且高阶WENO重构满足符号性质.对用来控制磁场散度的源项采用中心格式离散,最终得到与熵守恒通量一致的高阶精度.几个一维、二维算例表明该格式无振荡,鲁棒性强,可以精确捕捉间断.     

求解多维Euler方程的二阶旋转混合型格式

提出了一个基于旋转Riemann求解器的二阶精度的Euler（欧拉）通量函数.不同于“网格相关”的有限体积方法或者维数分裂的有限差分方法,本格式是基于旋转Riemann求解器将HLLC格式与HLL格式进行特定结合而得到的一类混合型数值格式.在激波法向采用HLL格式从而抑制红斑现象,在激波方向采用HLLC格式从而避免产生过多的耗散.新的旋转混合型格式具有结构简单、无红斑、高分辨率等优点.数值算例充分说明了新格式消除Euler方程激波不稳定现象的有效性和鲁棒性.     

一种求解Burgers方程的高精度紧致差分格式

针对Burgers方程,采用余项修正法和欧拉公式,推导了一种新的四层高精度紧致差分隐格式,其截断误差为O(τ~2+τh~2+h~4),即当τ=O(h~2)时,格式空间具有四阶精度;然后通过数值实验验证了格式的精确性和可靠性.     

求解一维对流方程的高精度紧致差分格式

本文提出数值求解一维对流方程的一种两层隐式紧致差分格式,采用泰勒级数展开法以及对截断误差余项中的三阶导数进行修正的方法对时间和空间导数进行离散.格式的截断误差为O(τ~4+τ~2h~2+h~4),即该格式在时间和空间上均可达到四阶精度.利用von Neumann方法分析得到该格式是无条件稳定的.通过数值实验验证了本文格式的精确性和稳定性.     

一个求解Euler方程的特殊矩阵分裂格式

§1.引言 自[1]提出矢通量分裂格式以来,在求解气动方程方面得到广泛应用。矢通量分裂格式是一种求解守恒型双曲方程组的方法,它将方程中代表质量、动量和能量的矢通量按照矢通量Jacobian矩阵正负特征值分裂为两个亚矢通量项,目的在于改进显式格式和隐式格式的计算效率和提高求解时的稳定性。在求解方法上,对于二维问题,需要求解以4×4块矩阵为矩阵元的上三角矩阵和下三角矩阵,比中心差分格式需要求解两个块三     

抛物型方程的一个高精度隐格式

利用待定参数法,对一维抛物型方程构造出了一个截断误差为Ｏ（△ｘ＾４＋△ｘ＾４）的隐式差分格式,格式的稳定性条件为ｒ＝ａ△ｔ／△ｘ＾２≤１／√２,可用追赶法求解。     

抛物型方程的一个新的高精度恒稳定的隐式差分格式

本文用待定参数法对一维抛物型方程构造出一个截断误差为 0 (△ t3+△ x6)的隐式差分格式 ,格式绝对稳定且可用追赶法求解 .     

High Accuracy Difference Schemes for the Cylindrical Heat Conduction Equation

High accuracy implicit difference methods are derived for thecylindrical heat conduction equation. Some unconditionally stableimplicit formulas are derived. The utility of the new schemesare shown by testing the schemes on two examples.     

High Accuracy Spectral Method for the Space-Fractional Diffusion Equations

In this paper, a high order accurate spectral method is presented for the space-fractional diffusion equations. Based on Fourier spectral method in space and Chebyshev collocation method in time, three high order accuracy schemes are proposed. The main advantages of this method are that it yields a fully diagonal representation of the fractional operator, with increased accuracy and efficiency compared with low-order counterparts, and a completely straightforward extension to high spatial dimensions. Some numerical examples, including Allen-Cahn equation, are conducted to verify the effectiveness of this method.     

高精度特征解及其外推求解位势方程

根据位势理论,基本边界特征值问题可转化为具有对数奇性的边界积分方程.利用机械求积方法求解特征值和特征向量,以及利用这些特征解求解Laplace方程.特征解和Laplace方程的解具有高精度和低的计算复杂度.利用Anselone聚紧和渐近紧理论,证明了方法的收敛性和稳定性.此外,还给出了误差的奇数阶渐近展开.利用h3-Richardson外推,不仅误差近似的精度阶大为提高,而且,得到的后验误差估计可以构造自适应算法.具体的数值例子说明了算法的有效性.     

Generalized Euler Equations

Zusammenfassung Es sei vorausgesetzt, dass in einem optimalen Steuerprozess die Anzahl der Steuervariabeln s–1 ist, wos die Anzahl der Zustandsvariabeln bedeutet. Es wird gezeigt, dass dann die adjungierten Variabeln eliminiert werden können, so dass ein System von Differentialgleichungen entsteht, das als eine Verallgemeinerung der Eulerschen Gleichungen der Variationsrechnung angesehen werden kann. Dieses Resultat gilt sowohl für diskrete als auch für kontinuierliche Steuersysteme, falls die Steuervariabeln keinen einschränkenden Nebenbedingungen unterworfen werden.     

Block Difference Schemes of High Order for Stiff Linear Differential-Algebraic Equations

Computational Mathematics and Mathematical Physics - The initial value problem for stiff linear differential-algebraic equations is considered. A block variant of multistep difference...     

A Survey of High Order Schemes for the Shallow Water Equations

In this paper, we survey our recent work on designing high order positivity-preserving well-balanced finite difference and finite volume WENO (weighted essentially non-oscillatory) schemes, and discontinuous Galerkin finite element schemes for solving the shallow water equations with a non-flat bottom topography. These schemes are genuinely high order accurate in smooth regions for general solutions, are essentially non-oscillatory for general solutions with discontinuities, and at the same time they preserve exactly the water at rest or the more general moving water steady state solutions. A simple positivity-preserving limiter, valid under suitable CFL condition, has been introduced in one dimension and reformulated to two dimensions with triangular meshes, and we prove that the resulting schemes guarantee the positivity of the water depth.     

High Accuracy Solution of Three-Dimensional Biharmonic Equations

In this paper, we consider several finite-difference approximations for the three-dimensional biharmonic equation. A symbolic algebra package is utilized to derive a family of finite-difference approximations for the biharmonic equation on a 27 point compact stencil. The unknown solution and its first derivatives are carried as unknowns at selected grid points. This formulation allows us to incorporate the Dirichlet boundary conditions automatically and there is no need to define special formulas near the boundaries, as is the case with the standard discretizations of biharmonic equations. We exhibit the standard second-order, finite-difference approximation that requires 25 grid points. We also exhibit two compact formulations of the 3D biharmonic equations; these compact formulas are defined on a 27 point cubic grid. The fourth-order approximations are used to solve a set of test problems and produce high accuracy numerical solutions. The system of linear equations is solved using a variety of iterative methods. We employ multigrid and preconditioned Krylov iterative methods to solve the system of equations. Test results from two test problems are reported. In these experiments, the multigrid method gives excellent results. The multigrid preconditioning also gives good results using Krylov methods.     

Multiparameter Schemes for Evolutionary Equations

Multiparameter extensions (MP) of (linear and nonlinear) descent methods have been proposed for the solution of finite dimensional time independent problems; these new methods are based on a different treatment of several blocks of components of the solution, basically via the substitution of a scalar relaxation by a (suitable) matricial relaxation. Similarly, the Nonlinear Galerkin Method (NLG), that stems from the dynamical system theory, propose to apply distinct temporal integration schemes to different sets of data scales when solving dissipative PDEs. In this paper, the algebraic similarity of Richardson iteration and Forward-Euler time integration is extended to new grounds through the expansion of the realm of MP methods to the field of the numerical integration of dissipative PDEs. The separation of the structures is realized by the utilization of hierarchical preconditioners in finite differences, which are conjugated to a MP temporal integration steeming from NLG theory. Numerical examples of fluid dynamics problems show the improved temporal stability of these new methods as compared to the classical ones.     

On the Construction of Combined Finite-Difference Schemes of High Accuracy

Doklady Mathematics - A method is proposed for constructing combined shock-capturing finite-difference schemes that localize shock fronts with high accuracy and preserve the high order of...