一个基于非均匀网格上的二阶基本无振荡差分格式

在基本无振荡格式的构造中，将通常的对流通量$f$的逼近方式推广到对通量导数的逼近，这一构造方法可以有效地应用到非均匀或非结构网格. 直接基于非均匀网格上，构造了一个二阶的基本无振荡(ENO)差分格式. 该格式具有形式简单，对网格的划分灵活，与传统格式相比不增加计算量等优点. 几个数值算例证明了格式的有效性.     

非均匀介质弹性波动方程的不规则网格有限差分方法

从弹性波动方程出发,提出了一种新的空间不规则网格有限差分方法,并用于求解非均匀各向异性介质中的弹性波正演问题。这种方法简单易行,对于复杂几何结构,例如低速层、套管井和非平面界面等,在较细的不规则网格上进行离散,计算时间和占用内存更少。与多重网格差分方法相比,该方法不需要粗、细网格之间的插值,所有网格差分计算在同一次空间迭代中完成。具有复杂几何交界面的模型计算,包括地下透镜体、套管井眼等,在确定弹性常数和密度后,用不规则网格的差分方法更易实现。该方法使用了Higdon吸收边界条件解决人工边界反射问题,引入了新的稳定性条件和网格频散条件,很好地消除了非物理散射波。理论模型的效值计算表明,该方法具有良好的稳定性和计算精度,在模拟非均匀介质弹性波传播时,比相同精度的规则网格有限差分方法计算速度更快。该方法易于推广到非结构网格和三维问题中。     

非线性双曲型守恒律的两步二阶TVD差分格式

通过Mac Cormack格式和Warming-Beam的结合,构造了一种非常简单的两步二阶TVD差分格式,该差分格式更适合于使用分量形式差分计算而无须对欧拉方程组进行特征解耦。通过对流体力学方程组的大量数值试验,并与二阶ENO格式进行了比较,充分显示了该格式高精度、高分辨并且极其简单的优良特性。     

加权基本无振荡格式研究进展

加权型基本无振荡WENO格式是近十年发展起来的一类高阶、高精度格式,它是在ENO格式的基础上采用加权思想构造的,对流场内的间断和细致结构具有较高的分辨率,适于求解包含激波、膨胀波以及接触间断等复杂结构的流场,目前已发展成为计算流体力学中一种重要的方法。本文针对加权型基本无振荡格式近年来的进展作一简要介绍。     

对“二阶迎风TVD数值格式在非平衡高超声速钝体绕流中的运用”一文的讨论

对“二阶迎风ＴＶＤ数值格式在非平衡高超声速钝体绕流中的运用”一文的讨论胡光初（中国空气动力研究与发展中心，四川绵阳６２１０００）１９９６１１１１收到．为了方便讨论，将我们于１９９４年３月在《力学学报》第２期上发表的“二阶迎风ＴＶＤ数值格式在非平衡...     

概率密度演化方程TVD格式的自适应时间步长技术及其初值条件改进

随机性普遍存在于实际工程问题中,而复杂结构的非线性随机响应分析是其中的一个难点,近年发展的概率密度演化方法为此类问题的求解提供了新的途径.由于实际问题的复杂性,概率密度演化方程通常采用数值方法求解,因此提高计算效率和求解精度对实际应用具有重要意义.本文基于变网格技术,推导了概率密度演化方程在非均匀时间步长上的总变差减小(total variation diminishing,TVD)差分格式,算例结果表明通过自适应插值可将迭代次数减少为原来的43.4%,当随机过程样本持续时间增大时均值估计的平均误差基本不变,而标准差估计的平均误差不断增大,但增大幅度不断减小;计算耗时随样本持续时间的增大也呈增大趋势,而由于使用了时间步长自适应插值算法导致有些情况下长持时样本的计算耗时反而比短持时样本的计算耗时短;在传统的脉冲函数型初值条件基础上,提出了一种高阶导数更稳定的余弦函数型初值条件形式.结果表明,脉冲函数型的初值条件是余弦函数型初值条件的一个特例,当参数取值适当时,余弦函数型初值条件的数值求解结果具有更高的精度.本文的工作进一步完善了概率密度演化方程的求解方法,为其在实际工程中的应用提供了基础.      

河道砂油藏的自适应非均匀网格粗化算法

以河道砂的观测深度为确定性数据，由贝叶斯理论通过随机楚模的方法楚立横截面为抛物线形状的河道砂油藏边界面，并将渗透率自适应网格技术应用于河道砂油藏的网格粗化算法中。在渗透率或孔隙度交化异常区域自动采用精细网格，用直接解法求解渗透率或孔隙度交化异常区域的压强分布，而在其他区域采用不均匀网格粗化方法计算，印在流体流速大的区域采用精细网格。用本文方法计算了河道砂油藏的压强分布，结果表明河道砂油藏的三维不均匀自适应网格粗化算法的解在渗透率或孔隙度异常区的压强分布规律更逼近采用精细网格的解，在其他区域压强分布规律非常逼近粗化算法的解，但计算的速度比采用精细网格提高了100多倍。     

两种高分辨率数值方法计算波前有非均匀流的变截面管内激波传播规律的比较

本文采用两种对间断解具有高分辨率的数值方法基于推广Riemann问题解的二阶Godunov型有限差分法和分裂算子的随机选取法,计算了微波在波前有非均匀定常流的一维变截面管道中的传播和波后流场特征,得到一致结果,用数值模拟方法揭示了这类运动的一些特殊规律。对比两种方法的计算过程和结果,可以看出,二阶Godunov型方法明显优于随机选取法。     

基于非结构/混合网格的高阶精度格式研究进展

尽管以二阶精度格式为基础的计算流体力学(CFD) 方法和软件已经在航空航天飞行器设计中发挥了重要的作用, 但是由于二阶精度格式的耗散和色散较大, 对于湍流、分离等多尺度流动现象的模拟, 现有成熟的CFD 软件仍难以给出满意的结果, 为此CFD 工作者发展了众多的高阶精度计算格式. 如果以适应的计算网格来分类, 一般可以分为基于结构网格的有限差分格式、基于非结构/混合网格的有限体积法和有限元方法,以及各种类型的混合方法. 由于非结构/混合网格具有良好的几何适应性, 基于非结构/混合网格的高阶精度格式近年来备受关注. 本文综述了近年来基于非结构/混合网格的高阶精度格式研究进展, 重点介绍了空间离散方法, 主要包括k-Exact 和ENO/WENO 等有限体积方法, 间断伽辽金(DG) 有限元方法, 有限谱体积(SV) 和有限谱差分(SD) 方法, 以及近来发展的各种DG/FV 混合算法和将各种方法统一在一个框架内的CPR (correctionprocedure via reconstruction) 方法等. 随后简要介绍了高阶精度格式应用于复杂外形流动数值模拟的一些需要关注的问题, 包括曲边界的处理方法、间断侦测和限制器、各种加速收敛技术等. 在综述过程中, 介绍了各种方法的优势与不足, 其间介绍了作者发展的基于"静动态混合重构" 的DG/FV 混合算法. 最后展望了基于非结构/混合网格的高阶精度格式的未来发展趋势及应用前景.     

基于非结构网格的TTGC有限元格式的实现及在超声速流动中的应用

基于非结构网格系统,实现了时空三阶精度的TTGC有限元格式,并在三阶TTGC格式上发展了基于人工粘性的激波捕捉技术。在非结构网格下,采用这种方法对若干典型的超声速流动问题(SOD激波管、马赫数为3的前台阶流动以及马赫数为8的高超声速圆柱流动)进行了验证计算。结果表明,TTGC格式分辨率高,在粗糙网格下能够准确的模拟超声速流场中的激波、接触间断等复杂流动现象,并且能有效的控制间断附近的数值色散现象。与传统的有限体积方法相比,本文实现的TTGC有限元格式在模拟超声速流动问题方面具有格式精度高、数值耗散小等优点。     

二阶迎风TVD数值格式在非平衡高超声速钝体绕流中的应用

将Ｈ．Ｃ．Ｙｅｅ所发展的一种二阶迎风ＴＶＤ格式运用于半球的高超声速化学非平衡黏性绕流的数值计算，详细介绍了流体力学方程组和化学组元守恒方程组全耦合求解的过程，以前这一研究领域中仅仅运用一阶格式及松散耦合的方法。化学松弛气体考虑了七组元（Ｏ，Ｏ＿２，ＮＯ，Ｎ，ＮＯ￣（＋），Ｎ＿２，ｅ￣（－））和主要的七反应，在驻点线上的计算结果与前人结果比较表明，本文有较大改进。     

AN ACCURATE FINITE DIFFERENCE SCHEME FOR SOLVING CONVECTION-DOMINATED DIFFUSION EQUATIONS

Approximating convection-dominated diffusion equations requires a very accurate scheme for the convection term. The most famous is the method of backward characteristics, which is very precise when a good interpolation procedure is used. However, this method is difficult to implement in 2D or 3D. The goal of this paper is to show that it is possible to construct finite difference schemes almost as accurate as the method of characteristics. Starting from a family of second- and third- order Lax–Wendroff-type schemes, a TVD and L^∞- stable scheme that is easy to implement in higher dimensions is constructed. Numerical tests are performed on various model problems whose solution is known and on classical problems. Comparisons with some other limiter schemes and the method of characteristics are discussed. © 1997 by John Wiley & Sons, Ltd.     

用物理黏性构建高阶不振荡对流扩散差分格式

利用数值摄动算法, 通过扩散格式数值摄动重构把对流扩散方程的2阶中心差分格式(2-CDS)重构为高精度高分辨率格式, 解析分析和模型方程计算证实了新格式的高精度不振荡性质. 新格式是把物理黏性使流动光滑化的扩散运动规律引入2-CDS 中的结果. 该法显然与构建高级离散格式的常见方法不同. 证实: 数值摄动重构中引入扩散运动规律的结果格式与引入对流运动规律(下游不影响上游的规律)的结果格式一致, 说明对离散方程的数值摄动运算, 在维持原格式结构形式不动的条件下, 不仅能提高格式精度和稳健性, 且可揭示对流离散运动规律与扩散离散运动规律之间的内在关联;同时证实, 文中提出和使用的上、下游分裂方法是构建高精度不振荡离散格式的一个有效方法.     

群桩作用机理的数值试验研究

 在分析桩土支撑体系及其相互作用关系的基础上，利用有限元分析软件ANSYS建立了群桩 体系计算模型. 通过对均质土体例题计算结果分析，可以看出利用该建模方法与传统解答有 很好的一致性. 根据辽宁工程技术大学实验馆场地土实测资料，利用该建模方法，求解出了承台与桩的 荷载分担比，对安全、经济地进行群桩基础设计具有重要指导作用和应用价值.     

A NEW HIGH-ORDER ACCURACY EXPLICIT DIFFERENCE SCHEME FOR SOLVING THREE-DIMENSIONAL PARABOLIC EQUATIONS

1.IntroductionInalotoffields,suchas'diffusion,seepageandheatconductionetc.,weoftenmeetproblemsofsolvingparabolicequations.Inthecaseofthree-dimension,themathematicalmodelissuchaninitialandboundaryvalueproblemasbelow:Thedifferenceschemesforsolvingtheproblemabove,whichhavehigh-orderaccuracyandcanbefiguredoutwithexplicitschemearesuchschemesasthoseinII]and[21.Theexplicitdifferenceschemeproposedinthispaperkeepsthehigh-orderaccuracy(it'hasthesameaccuracyasthosein[ifand[21),butitsexpressionismuchsimp…     

利用多小波自适应格式求解流体力学方程

高阶计算格式的高精度、高分辨率对提高复杂流场的计算水平有重要的意义,为了提高AUSMPW格式对流场计算中激波等间断的分辨率,减小数值振荡,在原有AUSMPW格式的基础之上,利用多小波对函数进行多尺度分解,并采取阈值的方法生成自适应网格,提出了一种新的基于多小波自适应算法的AUSMPW格式,理论上可以达到任意阶精度. 将所得的压强、密度与原格式、TVD格式及WENO格式的计算结果进行了比较分析. 结果表明改进后的AUSMPW格式较原格式具有更高的分辨率、更强的捕捉间断的能力及更低的数值耗散.      

有应力史差非匀质截面轴心受压杆的性能研究

对有应力史差、非匀质、非线弹性的轴心受压构件进行了理论分析和试验研究，初步揭示了此类构件的受力特点．理论分析和试验表明：此类构件截面承载力不是各部分材料承载力的简单叠加；承载力和截面破坏控制点因截面上材料的本构关系的变化和应力史差而发生改变．     

N-S方程在非结构网格下的求解

在Ｒｏｅ的矢通量差分分裂的基础上，吸收了ＮＮＤ格式的优点，提出了一种非结构网格下求解Ｅｕｌｅｒ方程和Ｎ－Ｓ方程的高分辨率高精度迎风格式．这种格式具有捕捉强激波和滑移线的良好性能．在时间方向上采用了显式和隐式两种解法．文中还给出了自适应技术．最后，成功地完成了ＧＡＭＭ超音速前台阶绕流、二维平板无粘激波反射、三维Ｈｏｂｓｏｎ叶栅流动、ＶＫＩ叶栅流动、Ｃ３Ｘ叶栅流动的数值模拟，得到了满意的结果