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

看了《力学学报》第26卷第2期（1994年3月）发表的胡光初[1]等同志的论文“M阶迎风TVD数值格式在非平衡高超声速钝体绕流中的应用”一文以后，关于计算精度和结果可靠性问题。我们提出以下不同的意见．在文问中多处讲到关于格式的精度问题，并一再强调以前这项研究工作别人用的     

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

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

对流项占优问题的MLPG/SUPG方法数值模拟

在计算对流项占优问题时易产生假扩散,本文把流线型迎风格式应用于MLPG方法中可以减少对流项的影响,通过两个典型例子（旋转流场问题和Brezzi问题）验证该格式的精度与有效性,并与文献中的迎风格式的计算结果进行比较,计算结果表明,该方法能有效地克服假扩散现象,有较好的稳定性和较高的计算精度。     

变量的两种独立性和哈密顿原理的无穷多样性

变量的两种独立性和哈密顿原理的无穷多样性朱如曾，向程（中国科学院力学研究所，北京１０００８０）最近《力学与实践》发表了“对广义动量的探讨”一文￣［１２］，该文批评了９本力学著作￣［１－２，４－１０］的有关论述￣１）．如果文［１２］的批评是正确的，那么...     

变量的两种独立性和哈密顿原理的无穷多样性

变量的两种独立性和哈密顿原理的无穷多样性朱如曾，向程（中国科学院力学研究所，北京１０００８０）最近《力学与实践》发表了“对广义动量的探讨”一文￣［１２］，该文批评了９本力学著作￣［１－２，４－１０］的有关论述￣１）．如果文［１２］的批评是正确的，那么...     

屈曲杆最大挠度近似公式的再改进

屈曲杆最大挠度近似公式的再改进朱华满（大连铁道学院大连１１６０２２）两端铰支细长压杆最大挠度近似公式已有多种形式￣［１－３］,精度较高、适用范围较大的当数文献［３］推荐的公式。本文导出的近似公式在精度和适用范围方面又有新的改进。欧拉从曲率精确式出发,...     

也论弹性理论中半无限平面问题的一些解答

也论弹性理论中半无限平面问题的一些解答薛福林（哈尔滨工业大学，哈尔滨１５０００１）文［１］认为一些弹性半无限平面问题的解是不正确的，理由是它们与用叠加法得到的解不同事实上，文［１］认为不正确的那两个解都是正确的．文献［２］１１２页的一段话具有指导意义...     

也论弹性理论中半无限平面问题的一些解答

也论弹性理论中半无限平面问题的一些解答薛福林（哈尔滨工业大学，哈尔滨１５０００１）文［１］认为一些弹性半无限平面问题的解是不正确的，理由是它们与用叠加法得到的解不同事实上，文［１］认为不正确的那两个解都是正确的．文献［２］１１２页的一段话具有指导意义...     

无限元与相同子结构串的一致性及其算法

讨论了无限元方法￣［１］与相同子结构串分析￣［２］的一致性，综合文［１］和［２］得到一种通用算法，依据文［３］尚可应用于ＬＱ控制问题，且提出一种直接求出转移矩阵的算法，具有简便和较高精度的优点．     

流动显示图像的处理和分析

流动显示图像的处理和分析李欣周建和周兴华（天津大学力学系，天津３０００７２）文［１］和［２］曾介绍过天津大学流体力学实验室所开设的两个流体力学新教学实验，文中提到用摄录像设备与计算机连接，将图像信息存入计算机内进行计算和处理，快捷方便．本文介绍天津大...     

TVD三阶和四阶格式求解渗流方程新方法

在使用显式格式计算含水率曲线过程中,如果空间步长太大,则会造成差分格式的不稳定.为了减弱空间步长带来的影响,建立了TVD三阶和四阶格式来计算两相流问题的含水率,与一点上游加权格式和两点上游加权格式进行比较,可以看出,用TVD格式计算的结果与实际结果更接近;该格式在稳定性方面也比上游加权格式稳定,Courant数的大小对格式的稳定性起决定作用.     

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

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

对流扩散方程的绝对稳定高阶中心差分格式

将作者提出的数值摄动算法改进为区分离散单元内上游和下游并分别对通量进行高精度重构的双重数值摄动算法,与原(单重)摄动算法相比,双重摄动算法既提高了格式精度又明显扩大了格式的稳定域范围.利用双重摄动算法,即分别利用上游和下游基点变量的摄动重构将高阶流体力学关系及迎风机制耦合进二阶中心格式之中,由此构建了对流扩散方程的对网格Reynolds数的任意值均稳定(绝对稳定)高精度(四阶和八阶精度)三基点中心TVD差分格式,通过解析分析以及3个算例计算证实了构建格式的优良性能;3个算例包括一维线性、非线性(Burgers方程)和二维变系数对流扩散方程.数值计算表明:构建的格式在粗网格下不振荡,构建格式在粗网格时的最大误差L_∞和均方误差L_2与二阶中心格式在细网格时的相应误差一致,对线性方程,构建格式在细网格下可达到L_2精度阶.     

Sine-Gordon方程的截断系统的同宿轨道

研究Ｓｉｎｅ Ｇｏｒｄｏｎ方程的广义渐近惯性流形上的常微分方程组，证实了在一定参数条件下存在Ｗｉｇｇｉｎｓ［１］意义下的同宿轨道．计算表明，与Ｂｉｓｈｏｐ［２］用数值计算得到的Ｓｉｎｅ Ｇｏｒｄｏｎ方程产生混沌的参数值尚有差别，考虑到同宿出现参数值往往低于混沌出现参数值，故结果在定性上正确，而且改进了文［１］中的结果．     

求解温盐双扩散系统的一种高精度方法

引入边界拟合坐标系来研究温盐双扩散系统。为了提高求解的精确性 ,对流项采用四阶精度的迎风格式 ,扩散项和涡量方程的浮力项则采用四阶精度的中心差 ,因此本文的方法是高精度的方法。首先针对温度占优( Rρ=0 .32 )和盐度占优 ( Rρ=1 .6 8)的情形进行了验证性计算 ,得到了与前人一致的结果。进一步 ,本文系统给出了不同的盐度通量强度下的流动形态 ,包括对称结构 ,不对称结构 ,反转结构等 ,结果与前人的吻合。     

四种陶瓷材料与SUS304不锈钢的高温摩擦学特性研究

为了探索陶瓷与金属组合作为高温润滑材料的可能性，利用端面摩擦磨损试验机测定了４种陶瓷ＳｉＣ、Ｓｉ３Ｎ４、Ａｌ２Ｏ３和ＺｒＯ２与ＳＵＳ３０４不锈钢在室温至５００℃下的摩擦学性能．摩擦试验结果表明，ＳｉＣ、Ｓｉ３Ｎ４和Ａｌ２Ｏ３在低于２００℃时的摩擦系数稳定且都低于０．２，但在２００℃以上时的摩擦性能却都不稳定，摩擦系数在０．２－０．４之间；ＺｒＯ２在２００℃以下时的摩擦性能不稳定，而在２００℃以上时的摩擦系数低于６．２且较稳定。磨损试验结果表明，在４种陶瓷中ＺｒＯ２的磨损率最低［－２．６０×１０－９ｍｍ３／（Ｎ·ｍ）］，ＳｉＣ和Ｓｉ３Ｎ４的磨损率居中［分别为１．８０×１０－６ｍｍ３／（Ｎ·ｍ）和４．４０×１０－６ｍｍ３／（Ｎ·ｍ）］，Ａｌ２Ｏ３的磨损率最高［３．６４×１０－５ｍｍ３／（Ｎ·ｍ）］；分别与这４种陶瓷对磨的不锈钢的磨损率都高［１．４０×１０－５－４．５２×１０－５ｍｍ３／（Ｎ·ｍ）］．     

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

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

基于比拟理论的翼型扰流声场数值模拟

从二维模型方程的全离散形式出发,重点分析了差分格式的色散特性和各向异性效应,证实迎风紧致格式比对称格式有更好的色散和各向同性特性,故有利于声场的数值模拟,并采用三阶迎风紧致格式（ＵＣＤ３）和四阶对称紧致格式（ＳＣＤ４）计算了绕ＮＡＣＡ００１２翼型的可压缩非定常流场,并将此流场作为近场声源,运用声学比拟理论对气动声进行模拟。     

求解双曲型守恒律的半离散三阶中心迎风格式

给出了求解一维双曲型守恒律的一种半离散三阶中心迎风格式,并利用逐维进行计算的方法将格式推广到二维守恒律。构造格式时利用了波传播的单侧局部速度,三阶重构方法的引入保证了格式的精度。时间方向的离散采用三阶TVD Runge—Kutta方法。本文格式保持了中心差分格式简单的优点,即不需用Riemann解算器,避免了进行特征分解过程。用该格式对一维和二维守恒律进行了大量的数值试验,结果表明本文格式是高精度、高分辨率的。     

层状半空间中单桩对任意入射地震波的响应

本文以任意入射地震波作为初始运动输入，采用层状半空间中显式表达形式的动力Ｇｒｅｅｎ函数形成阶－土相互作用分析所需的土介质柔度矩阵，用梁单元模拟桩的运动，在文献［５］的基础上，建立了层状半空间中单桩对任意入射地震波响应的实用计算模型，研究了层状半空间中单桩对不同角度入射ＳＨ波、ＳＶ波和Ｐ波的动力响应，得到了一些新的结果。