求解二维Euler方程的时一空守恒格式

将作者原来得出的一维时－空守恒格式推广到了二维情形,得到了二维Ｅｕｌｅｒ方程的时．空守恒格式,并用几个典型算例进行了检验计算,结果表明：得到的二维时一空守恒格式保留了一维格式所有的优点,格式简单,通用性强,对微波等间断具有很高的分辨率．     

求解二维Euler方程的时一空守恒格式

将作者原来得出的一维时－空守恒格式推广到了二维情形,得到了二维Ｅｕｌｅｒ方程的时．空守恒格式,并用几个典型算例进行了检验计算,结果表明：得到的二维时一空守恒格式保留了一维格式所有的优点,格式简单,通用性强,对微波等间断具有很高的分辨率．     

一般坐标系下的二维Euler方程时—空守恒格式

将文「１－３」中的一维时－空守恒格式推广到了二维情形,得到了一般坐标系下的二维Ｅｕｌｅｒ方程时－空守恒格式,并用几个典型算例进行了检验计算,结果表明：本文得到的二维时－空守恒格式保留了一维格式所有的优点,格式简单,通用性强,而且对激波等间断具有很高的分辨率。     

一般形式的二维时—空守恒格式

本文将经作者改进后的一维时－空守恒格式推广到了二维情形，得到了一个一般形式的二维Ｅｕｌｅｒ方程时－空守恒格式，该格式对各种不规则几何区域内的流动问题具有很强的适应性，同时它还保留了一维格式的优点。几个典型算例的计算结果表明，本文格式不仅精度高，通用性好，而且对激波等间断具有很高的分辨率。     

二维时-空守恒格式及其在激波传播计算中的应用

将改进后的一维时－空守恒格式推广到了二维情形,得到了一个新的一般形式的二维Ｅｕｌｅｒ方程时－空守恒格式,并用该格式对几个具有复杂波系的流场进行了数值模拟。结果表明,该格式保留了一维格式通用性好、结构简单的优点,其计算结果精度高,对激波等间断具有很强的分辨率。     

求解双曲守恒律方程的熵相容格式

为了更好地求解流体动力学中的双曲守恒律方程,本文提出了一种熵相容格式.通过分析单元内跨越激波时熵的产生情况,得到熵产的显式表达式.在熵守恒通量中加入耗散项与熵增项获得熵相容格式的通量,并在此基础上加入限制器构造出高分辨率熵相容格式.在一维浅水波方程与一维相对论力学方程基础上对新格式进行检验,数值模拟结果表明:这种新的格...     

动力学平衡方程的Euler中点辛差分求解格式

给出了动力学方程${\pmb M}\ddot {\pmb x} + {\pmb C}\dot {\pmb x} + {\pmb K \pmb x} = {\pmb R}$的二阶Euler中点隐式差分求解格式,分保守系统、无 阻尼受迫振动系统和阻尼系统3种情况, 讨论了算法中Jacobi矩阵${\pmb A}$的性质,譬 如${\pmb A}$是否为辛矩阵以及谱半径等. 对于无阻尼系统,证明了无论是否存在外 载荷,Jacobi 矩阵都是辛矩阵. 证明了辛矩阵的所有本征值的模为1,其谱半径永远 为1, 以及$\delta = 0.5$和$\alpha = 0.25$的Newmark算法就是Euler中点隐式差 分格式,对保守系统它们都是辛算法. 严格证 明了Euler中点辛格式是严格保持系统能量的. 通过算例详细讨论了保辛算法用于求解非保 守系统动态特性的优越性,如广义保结构特性等;分析了保辛算法的相位误差以及由其引起 的系统的附加能量特性;分析了保辛算法和$\delta \ne 0.5$的Newmark算法的精度随着激励频率与系统固有频率比的变化情况等     

求解双曲守恒律方程的WENO型熵相容格式

通过在单元交界面处进行高阶WENO重构,得到了一种求解双曲型守恒律方程的WENO型熵相容格式。用该格式对一维Burgers方程和Euler方程进行数值模拟,结果表明,该格式具有高精度、基本无振荡性等特点。

     

二维边界层方程的迭代求解

针对二维边界层方程,提出了分析积分迭代法.首先将该方法应用于Blasius方程和Falkner-Skan方程的求解,数值计算结果稳定,计算精度高;然后对外部有势流不能达到自相似要求、必须二维求解的二维层流边界层问题,在分析积分迭代法中加上计算力学的松弛迭代法,形成了一套有效的算法.数值结果表明该方法用于层流边界层的计算是很有效的.     

二维Euler方程无网格算法的精度分析

在几种典型的点云布点方案中,分析了文献[3]中的无网格算法的计算精度。分析指出影响无网格算法计算精度的主要因素是点云中点的分布是否平衡;增加点云中节点的数目或者采用高阶曲面拟合都不能提高精度;这里的理论分析为无网格算法的布点提供了重要参考。     

求解可压缩流的二维通量分裂格式

传统的一维通量分裂格式在计算界面数值通量时,只考虑网格界面法向的波系。采用传统的TV格式分别求解对流通量和压力通量。通过求解考虑了横向波系影响的角点数值通量来构造一种真正二维的TV通量分裂格式。在计算一维数值算例时,该格式与传统的TV格式具有相同的数值通量计算公式,因此其保留了传统的TV格式精确捕捉接触间断和膨胀激波的优点。在计算二维算例时,该格式比传统的TV格式具有更高的分辨率;在计算二维强激波问题时,消除了传统TV格式的非物理现象,表现出更好的鲁棒性;此外,该格式大大提高了稳定性CFL数,从而具有更高的计算效率。因此,本文方法是一种精确、高效并且具有强鲁棒性的数值方法,在可压缩流的数值模拟中具有广阔的应用前景。     

求解二维污染扩散方程的时空任意阶的三点紧致显格式

通过在泰勒级数展开中运用逐阶迭代的方法,推导出了空间二阶导数任意精度的三点紧致的表达式,并在半高散方程中通过二维扩散方程本身把时间导数转换为空间导数,从而推导出了时空任意阶的三点紧致显格式.数值实验表明,本文格式的精度很高,而且具有使用简单,易于编程的优点,对求解二维污染扩散方程具有很好的应用前景.     

二维泊松方程的两步预估校正格式

为避免用四阶紧致格式求解泊松方程所具有编程复杂和难以实现的困难,对传统的五点二阶中心差分格式进行改进;通过增加对残差的校正计算,提出了一种新型具有四阶精度的两步预估校正格式. 新格式虽需要增加一定的计算量,但它的格式精度高,边界条件处理极简单,易于编程实现. 通过数值实验,结果证明上述格式的确具有易于编程和计算精度高的优点. 预估校正格式很容易推广到其他复杂情形.     

通过Adomian分解法求解二维Helmholtz方程

提出基于Adomian分解法求解二维Helmholtz方程。通过Adomian分解法可以把Helmholtz微分方程和边界条件分别转换成递归代数公式和适用符号计算的简单代数公式。利用边界条件可以很容易得到方程的解析解表达式。Adomian分解法的主要特点在于计算简单快速,并且不需要进行线性化或离散化。最后给出数值实例以验证Adomian分解法求解二维Helmholtz方程的有效性。通过数值计算可以发现,基于Adomian分解法的计算结果非常接近精确解,并且该方法具有良好的收敛性。这表明Adomian分解法能够快速有效求解Helmholtz方程。     

通过Adomian分解法求解二维Helmholtz方程

提出基于Adomian分解法求解二维Helmholtz方程。通过Adomian分解法可以把Helmholtz微分方程和边界条件分别转换成递归代数公式和适用符号计算的简单代数公式。利用边界条件可以很容易得到方程的解析解表达式。Adomian分解法的主要特点在于计算简单快速,并且不需要进行线性化或离散化。最后给出数值实例以验证Adomian分解法求解二维Helmholtz方程的有效性。通过数值计算可以发现,基于Adomian分解法的计算结果非常接近精确解,并且该方法具有良好的收敛性。这表明Adomian分解法能够快速有效求解Helmholtz方程。     

二维双曲型守恒律的一种五阶松弛格式

松弛格式是Jin和Xin提出的无振荡有限差分方法,其主要思想是将守恒律转化为松弛方程组进行求解.本文用逐维五阶WENO重构和显隐式Runge-Kutta方法对松弛方程组的空间和时间进行离散,得到了一种求解二维双曲型守恒律五阶松弛格式.所得格式保持了松弛格式简单的优点,不用求解Riemann问题和计算通量函数的雅可比矩阵.通过二维Burgers方程和二维浅水方程的数值算例验证了格式的有效性.     

谱元方法求解二维不可压缩Navier-Stokes方程

利用有限元的思想并结合谱方法的精度提出求解偏微分方程的谱元方法,在元素内插值函数使用伪谱Chebyshev逼近,并将此方法应用于求解不可压Navier-Stokes方程,具体求解了二维方腔顶盖驱动流,与公认基准解对比获得了较好的结果。     

高阶紧致格式求解二维粘性不可压缩复杂流场

提出了一种求解二维不可压缩复杂流场的高精度算法．控制方程为原始变量、压力Ｐｏｉｓｓｏｎ方程提法．在任意曲线坐标下，采用四阶紧致格式求解Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组，时间推进采用交替方向隐式（ＡＤＩ）格式，在非交错网格上用松弛法求解压力Ｐｏｉｓｓｏｎ方程．对于复杂的流场，采用了区域分解方法，并在每一时间步对各子域实施松弛迭代使之能精确地反映非定常流场．利用该算法计算了二维受驱空腔流动，弯管流动和垂直平板的突然起动问题．计算结果与实验结果和其他研究者的计算结果相比较吻合良好．对于平板起动流动，成功地模拟了流场中旋涡的生成以及Ｋａｒｍａｎ涡街的形成     

用拟压缩性方法求解不可压Euler方程

用拟压缩性方法和Ｊａｍｅｓｏｎ的有限体积算法求解了二维和三维定常可可压Ｅｕｌｅｒ方程。分别采用显、隐式时间离散推进求解;分析了人工粘性的阶数对定常解收敛性的影响,应用该方法计算了单个翼型和翼身组合体的低速绕流,结果与实验吻合较好。     

二维对流扩散方程的欧拉—拉格朗日分裂格式

本文在[1]基础上发展了一种有效的处理大P_e(R_e)数、非定常二维对流扩散方程的欧拉-拉格朗日(E-L)分裂格式,由于方法本质上与区域形状无关,且不需再分网格,因此是一种无网格的E-L方法,特别对于定常流动,E.-L.分裂格式可以导致比一阶迎风格式更精确的单调、无振荡格式,文中对于常系数、变系数和非线性的二维非定常和定常对流扩散方程的(初)边值问题进行了数值计算,数值结果与精确解的比较表明,本方法具有很好的精度,解是单调无振荡的,比通常一阶迎风格式具有较少的数值扩散,最大计算网格P-e(R-e)数可达100—500。