一维可压缩Euler方程组的两个模型

作者考察了一维可压缩Euler方程组的两个模型.利用特征分解和Gronwall不等式,首先得到具有几何结构且绝热指数γ=3的一维可压缩Euler方程组L~∞模的一致有界性.进一步,考虑当绝热指数γ=-1时,一维非等熵可压缩Euler方程组的Cauchy问题.在适当的假设下,得到该系统的整体经典解.     

基于不同数值流通量求解可压缩Euler方程组的Lax-Wendroff控制体积方法（英文）

基于Lax-Wendroff时间离散的控制体积间断Petrov-Galerkin方法是求解双曲守恒律的一种高精度和高分辨率数值方法.本文通过几个数值算例对8种数值流通量的数值表现作了详尽的比较,内容涉及耗时、精度、分辨率以及模拟复杂波形相互作用的能力.     

一维双曲守恒方程组的Taylor-Galerkin有限元方法

１．引言 在文章［８］中,利用双曲守恒律的Ｈａｍｉｌｔｏｎ－Ｊａｃｏｂｉ方程形式,应用 Ｇａｌｅｒｋｉｎ有限元给出了求解一维双曲守恒律的计算方法．不同于间断有限元方法［２］、［３］和 Ｔａｙｌｏｒ－Ｇａｌｅｒｋｉｎ有限元方法［１］求解双曲守恒律,文章［８］采用连续 Ｇａｌｅｒｋｉｎ有限元求解双曲守恒律． 在文章［８］中,对差分方法和有限元方法求解双曲守恒律作了较为详细的讨论．同时在文章［８］中,采用积分变换,将双曲守恒律方程变成 Ｈａｍｉｌｔｏｎ－Ｊａｃｏｂｉ方程形式．由于 Ｈａｍｉｌｔｏｎ－Ｊａｃｏ…     

一种守恒型间断跟踪法用于处理一维Euler方程组中的若干问题

1引言间断跟踪法(front-tracking)是数值求解双曲型守恒型方程(组)的一种重要的数值方法,其主要特点是把间断作为移动的内边界来处理,光滑区域中的数值解用计算光滑解有效的数值方法来求解,而间断的移动和间断两侧的数值解的修正要满足Rankine-Hugoniot条件．我们称这样的跟踪法为传统的间断跟踪法(见[3],[14])．本文的第二作者多年来研究设计了一种基于解的守恒性质的间断跟踪法(见[11],[12]),其最主要的特点是以解的守恒性作为跟踪的机制,而不是象传统的间断跟踪法那样利用     

一维多介质可压缩流动的守恒型界面追踪方法

本文针对一维问题的ProntTracking方法,提出了一种较易实现的守恒型界面追踪方法.利用双波近似求解Riemann问题来确定界面处的数值通量,在固定的网格上采用统一的有限体积格式进行内点和交界面点的计算,通过守恒插值以及守恒量的重新分配,保证数值解在全场实现一致守恒,将该方法应用于一维多介质可压缩流动的模拟,给出了满意的数值模拟结果.     

不可压磁流体力学方程组的高精度数值解法

为数值求解低雷诺数下不可压流体在电磁场作用下的流动,提出一种四阶紧致差分方法.由二维原始变量的MHD方程组出发,推导出具有较少未知量的电流密度-涡量-流函数形式MHD方程组.建立了求解二维非定常不可压MHD方程组的电流密度-涡量-流函数形式的四阶精度紧致差分格式.为验证本文提出的高精度紧致差分方法的精确性和可靠性,对有...     

二维可压缩Euler方程组C^1解整体存在的必要条件

对于具有某类初值条件的二维可压缩流体Euler方程组，给出了其C1解整体存在的必要条件，从而对[1]、[2]中的“未解决问题”提供了有意义的说明．     

二维可压缩Euler方程组轴对称解的生命区间

本文应用非线性几何光学的方法，对二维可压缩Ｅｕｌｅｒ方程组给出了轴对称光滑解的生命区间的精细估计，此外，还指出了光滑解爆破的原因     

一维椭圆方程边值问题的高精度差分方法

本文研究一维椭圆方程边值问题的差分方法,利用Lagrange插值理论与积分因子技巧,发展了一套有效的高精度算法,对非等距节点和等距节点,其精度分别可达O(h~4)和O(h~5).数值结果显示了该方法的优越性     

求解Euler方程的高精度熵相容格式

熵相容格式相比于一般的熵稳定格式进一步控制了激波处的熵增量,一维情况下能有效消除膨胀激波及间断处的伪振荡等现象.对于Euler方程,可以通过对特征变量进行WENO重构以获得高阶熵相容格式的数值粘性项,然后与高阶熵守恒格式结合得到高精度熵相容格式,在WENO重构过程中的权重关于特征变量是非线性的,这导致了大量的向量内积运算.通过用压强和熵代替特征变量来计算权重,可以显著减少重构的计算量,并且数值算例表明这种权重的计算方式能很好地保持数值格式的高阶精度和基本无振荡的效果.     

HIGH ORDER NUMERICAL METHODS TO TWO DIMENSIONAL HEAVISIDE FUNCTION INTEGRALS

In this paper we design and analyze a class of high order numerical methods to two dimensional Heaviside function integrals. Inspired by our high order numerical methods to two dimensional delta function integrals [19], the methods comprise approximating the mesh cell restrictions of the Heaviside function integral. In each mesh cell the two dimensional Heaviside function integral can be rewritten as a one dimensional ordinary integral with the integrand being a one dimensional Heaviside function integral which is smooth on several subsets of the integral interval. Thus the two dimensional Heaviside function integral is approximated by applying standard one dimensional high order numerical quadratures and high order numerical methods to one dimensional Heaviside function integrals. We establish error estimates for the method which show that the method can achieve any desired accuracy by assigning the corresponding accuracy to the sub-algorithms. Numerical examples are presented showing that the second- to fourth-order methods implemented in this paper achieve or exceed the expected accuracy.     

可压缩流体力学高精度拉格朗日格式及其保正性质

本文对可压缩流体力学高精度拉格朗日格式及其保正性质近年来的发展给出回顾与综述.文中分别介绍了一维、二维可压缩流体力学方程中心型拉格朗日格式的设计步骤,回顾了高精度拉格朗日格式以及高精度保正拉格朗日格式的研究进展.     

THE OPTIMAL CONVERGENCE ORDER OF THE DISCONTINUOUS FINITE ELEMENT METHODS FOR FIRST ORDER HYPERBOLIC SYSTEMS

In this paper, a discontinuous finite element method for the positive and symmetric, first-order hyperbolic systems （steady and nonsteady state） is constructed and analyzed by using linear triangle elements, and the O（h^2）-order optimal error estimates are derived under the assumption of strongly regular triangulation and the Ha-regularity for the exact solutions. The convergence analysis is based on some superclose estimates of the interpolation approximation. Finally, we discuss the Maxwell equations in a two-dimensional domain, and numerical experiments are given to validate the theoretical results.     

WELL-CONDITIONED FRAMES FOR HIGH ORDER FINITE ELEMENT METHODS

The purpose of this paper is to discuss representations of high order C⁰finite element spaces on simplicial meshes in any dimension.When computing with high order piecewise polynomials the conditioning of the basis is likely to be important.The main result of this paper is a construction of representations by frames such that the associated L²condition number is bounded independently of the polynomial degree.To our knowledge,such a representation has not been presented earlier.The main tools we will use for the construction is the bubble transform,introduced previously in[1],and properties of Jacobi polynomials on simplexes in higher dimensions.We also include a brief discussion of preconditioned iterative methods for the finite element systems in the setting of representations by frames.     

结合气体分子动力学方法的RKDG有限元格式求解一维Euler可压方程组

In this paper, two numerical methods are developed for solving one-dimensionl compressible ELder equations by the RKDG finite element method.The schemesare obtained based on an important relation between the Boltzmann equation andthe ELder equations.The schemes have the TVD-like property under the uniform meshes.Several numerical results also present the performance of the schemes.     

高阶Duffing方程的周期解

在本文中,我们将二阶Duffing方程周期解存在的结果推广到高阶Duffing方程 x~((2n))+g(x)=e(t), (n≥1).     

高分辨KFVS有限体积方法及其CFD应用

1．引言文中研究三维Euler方程组的数值求解·儿1）中p,（。。,。。,。z）,p和E分别表示流体密度,流体速度矢量,压力和总能．方程组（1．1）是不封闭的,除非增加一个额外的方程一状态方程p一pk句,e表示单位质量内能．本文仅限于理想气体,此时状态方程为p一（、－1加e．队2）近H十年来,涌现了许多求解方程组（1．1）的无振荡、高分辨格式,例如TVD格式问,＊＊O格式问等,它们在一定程度上促进了航空航天和造船事业的发展．其中有一类根据双曲方程组（1．l）特征值的符号建立的迎风格式尤为突出,与中心格式相比,迎风格式的耗…     

一类二阶非线性微分方程的求积问题

讨论了二阶非线性微分方程ｙ″＋ｐ（ｘ）ｙ′＋Ｗ＇（ｙ）／Ｗ（ｙ）ｙ＇＾２＝Ｑ（ｘ）１／Ｗ（ｙ）Ｆ［ｘ，（ｙ＇Ｗ（ｙ）＾ａ］的可积性问题，提供了可积的一些充分条件，在Ｆ为某些特殊类型时，给出能解公式。     

解高阶抛物型方程的高精度显式差分格式

1 引言 1960年,Saul’ev在文中讨论了如下的高阶(2m阶)抛物型方程 μ/t=(-1)~(m 1)~2mμ/x~(2m) (1)(其中m为正整数),提出了一类含极因子α的两层差分格式。当α=0时为显式格式,其稳定性条件为,r=△t/(△x)~(2m)<1/2~(2m-1),△t,△x分别为时间及空间步长。随后,文[2],[3]利用     

三步五阶迭代方法解非线性方程组

本文根据求积公式, 给出了三种求解非线性方程组的迭代方法, 并证明了所提出的三步迭代方法具有五阶收敛性. 最后给出了四个数值实例, 将本文的实验结果与现有的几种迭代方法的实验结果作了比较分析, 表明本文所提出的方法具有明显的优越性.