强迫耗散非线性发展方程与完全平方守恒格式

从描述大气和海洋运动的强迫耗散非线性发展方程出发，对强迫耗散非线性大气和海洋方程组显式差分格式的计算稳定性进行了分析，构造了一类强迫耗散非线性发展方程的显式准安全平方守恒差分格式。理论分析和数值试验证明，这类显式准完全平方守恒差分格式是计算稳定的，值得推广应用。     

二维Euler方程组的高分辨率隐式差分格式的并行计算

本文首先对标量双曲型守恒律构造了一类新的二阶隐式ＴＶＤ差分格式，再利用Ｒｏｅ方法和隐式近似因子分解法推广格式到一维、二维守恒率方程组。为减少计算工作量，推导了二维Ｅｕｌｅｒ方程组隐格式的对角形式。最后，针对江南－Ⅱ并行机设计了适合于多处理机结构特征的并行计算程序，计算翼型的跨声速流场，结果令人满意。     

二维Euler方程组的高分辨率隐式差分格式的进行计算

本文首先对标量双曲型守恒律构造了一类新的二阶隐式ＴＶＤ差分格式，再利用Ｒｏｅ方法和隐式近似因子分解法推广格式到一维、二维守恒率方程组。为减少计算工作量，推导了二Ｅｕｌｅｒ方程组隐格式的对角形式。最后，针对江南－Ⅱ并行机设计了适合盱多处理机结构特征的并行计算程序，计算翼型的跨声速流场，结果令人满意。     

求解双曲守恒律初值问题的GAUSS格式

构造了一类计算双曲守恒律弱解的 3点二阶显式格式 ,这类格式在CFL条件数为 1的限制下为TVD格式 ,并将这类格式推广到方程组的情形 ,进行了数值试验 ,结果是满意的。     

一个基于通量分裂的高精度MmB差分格式

本文研究双曲型守恒律的高精度差分方法.在新的算法中,首先将计算区间划分为互不相交的小区间,再根据精度要求等分小区间;其次,根据流动方向进行通量分裂,重构小区间交界面上的正、负数值通量,并进行校正;然后,采用高阶Runge-Kutta TVD方法进行时间离散,构造了一维非线性双曲型守恒律方程的一个高精度、高分辨率的守恒型差分格式.推广到二维双曲型守恒律方程,证明了格式的MmB特性.进而推广到二维守恒型方程组情形.最后对二维Burgers方程及Euler方程进行了数值试验,数值结果令人满意.     

双曲守恒型方程的二阶摄动有限差分格式

对双曲守恒型方程，将其一阶迎风格式空间差商的常系数摄动展开为时间步长和空间步长的幂级数，通过确定幂级数系数而获得二阶精度的摄动有限差分(PFD)格式。进而从双曲守恒型方程的通量分裂型一阶迎风格式出发，通过类似的摄动展开方法，获得空间精度为二阶的通量分裂形式的摄动有限差分(FPFD)格式。这两类格式保留了一阶守恒迎风格式的简洁结构形式，使用三节点即可达到二阶精度，又避免了三点二阶格式的非物理数值振荡。并将这两类格式推广应用到双曲守恒型方程组，最后通过模型方程和一维激波管流动的数值算例验证了格式的高精度、高分辨率性质。     

一种改进的紧致WENO混合格式

给出了一种用于激波捕捉计算的守恒型混合紧致WENO格式.混合格式可视为五阶WENO格式和五阶守恒紧致格式的加权平均.在格式构造中,本文采用了间断分辨率更好的WENO权因子计算方法以及新的子格式权因子计算方法.对于Euler方程,仅对混合格式中WENO格式部分做特征投影处理,不仅获得了特征型WENO格式良好的间断分辨率,同时又保证了计算效率.数值实验验证了本文混合格式的高时间效率和高流场分辨率.     

多守恒性的差分格式的分析和比较

球面大气浅水波方程(在忽略摩擦和外源强迫的条件下)具有总能量、总质量、总位涡、总位涡拟能和总角动量守恒等五个重要的物理守恒性.然而,将大气浅水波方程差分离散化后,常常不能保持这些守恒性.本文对新发展的两种较好的差分格式进行比较;一种可称之为修正总能量守恒格式,它可以较好地保持上述五个守恒性中的四个;而另一种可称之为准辛格式,它可以使上述五种守恒性均近似保持.文中对这两类格式作了具体分析,并作相应的数值试验和比较,结果表明这两类格式均值得推广应用.     

求解双曲守恒律方程的自适应人工黏性熵稳定格式

采用熵守恒格式求解双曲守恒律方程不满足熵稳定条件.为了达到熵稳定,需要在熵守恒格式的基础上增加黏性机制,使总熵耗散.鉴于自适应人工黏性在不同的计算区域上其黏性大小的变化特点,在熵守恒通量的基础上添加经过修正的自适应人工黏性通量构造出一种新的结构简单的熵稳定格式.根据自适应人工黏性的自适应性,并经过简单调节黏性比例系数C之后,该格式可以在不使用限制器的情况下达到整体的高分辨率,表现出:在间断区域有足够的数值耗散保证稳定;在光滑区域,黏性很小,不会影响格式在此的高精度.最后的几个数值算例以熵稳定ERoe格式和高分辨率熵稳定EYee格式的数值结果作为参照,分析说明这种格式的特性.     

二维非结构网格上的高精度有限体积WENO格式

构造了非结构网格上二维双曲型守恒律的一类新的高精度有限体积WENO格式。其主要思想是:根据格式精度的要求,按照谱体积方法对三角形单元网格进行剖分,通过选取适当的子单元组成模板,利用WENO重构方法重构二阶和三阶多项式,利用有限体积公式和高阶Runge-Kutta TVD时间离散方法,构造了非结构网格上二维双曲型守恒律的一致二阶和三阶精度的有限体积WENO格式。然后,推广到二维Euler方程组。最后,给出几个数值算例,验证了格式的稳定性、高阶精度和高分辨捕捉激波等间断的能力。     

S型进气道的数值模拟研究

 本工作改进了Denton的有限体积的时间推进法,将这种对时间一阶精度对空间二阶精度的方法改进为对时间二阶精度对空间二阶精度的计算方法,并应用于S型进气道的内流及具有进口外流的进气道流场的数值模拟,进行了实例计算,取得了合理的计算结果。     

计算Cauchy和Mangler积分主值的新方法

对于Cauchy积分主值曾提出过一些数值积分法,主要是抽去奇异性方法和内插型或Gauss型方法。对于前者,要求积分核在积分区间上有连续的一阶导数;对于后者,实际上只能求出节点(某类函数零点)处的积分主值,且当积分核不是多项式时     

发动机短舱防冰性能计算方法研究综述

发动机短舱中可能出现的结冰现象对飞行安全有着严重的影响。总结了已有关于短舱防冰性能计算的文献资料,在明确短舱防冰性能计算的理论依据的基础上,分析了其与机翼防冰性能计算的异同及计算中可能遇到的问题及解决方法,并对国内发动机短舱防冰性能计算提出了合理化建议。     

数控加工零件热场分布的虚拟仿真

利用有限差分法对金属零件在数控机械加工过程产生热量的分布进行了数值分析,为计算金属零件热应力变形提供基础,是虚拟制造中机械加工工艺过程规划虚拟仿真的一种新方法。     

CAXA制造工程师在复杂模具数控加工中的应用

通过对可乐瓶底模具运用CAD／CAM技术进行数控加工这一实例，介绍了CAXA制造工程师软件在数控加工制造中应用的具体步骤：首先用CAXA制造工程师的CAD功能生成3D模型，再用其CAM功能自动生成数控加工程序及刀具轨迹，并进行仿真加工后再校核、修订，从而获得优化的实际生产结果；强调了CAD／CAM在现代制造技术中的地位及重要性，说明用CAD／CAM技术缩短设计和制造周期，提高零件加工质量的目的。     

固液相变问题的边界元法数值解

相变问题的数值计算方法有多种,本文用边界元法解应用工程中常见的相变导热问题。     

A Class of TVD Type Combined Numerical Scheme for MHD Equations With a Survey About Numerical Methods in Solar Wind Simulations

It has been believed that three-dimensional, numerical, magnetohydrodynamic (MHD) modelling must play a crucial role in a seamless forecasting system. This system refers to space weather originating on the sun; propagation of disturbances through the solar wind and interplanetary magnetic field (IMF), and thence, transmission into the magnetosphere, ionosphere, and thermosphere. This role comes as no surprise to numerical modelers that participate in the numerical modelling of atmospheric environments as well as the meteorological conditions at Earth. Space scientists have paid great attention to operational numerical space weather prediction models. To this purpose practical progress has been made in the past years. Here first is reviewed the progress of the numerical methods in solar wind modelling. Then, based on our discussion, a new numerical scheme of total variation diminishing (TVD) type for magnetohydrodynamic equations in spherical coordinates is proposed by taking into account convergence, stability and resolution. This new MHD model is established by solving the fluid equations of MHD system with a modified Lax-Friedrichs scheme and the magnetic induction equations with MacCormack II scheme for the purpose of developing a combined scheme of quick convergence as well as of TVD property. To verify the validation of the scheme, the propagation of one-dimensional MHD fast and slow shock problem is discussed with the numerical results conforming to the existing results obtained by the piece-wise parabolic method (PPM). Finally, some conclusions are made. This revised version was published online in August 2006 with corrections to the Cover Date.     

圆柱突然起动的快速涡解法

 文中将Chorin的涡方法、格子涡(VIC)方法和快速Fourier变换(FFT)相结合,发展了一套快速有效的求解高Re数下非定常N-S方程的算法(快速涡方法)。通过对二维圆柱突然起动问题的数值计算,得到流动分离及其演化过程(包括二次涡及二次涡对的生成);此外还得到回流区对称轴线上速度分布。这些结果与实验数据及其它数值解相比,是令人满意的。该程序对其它高Re数非定常流动问题也是适用的。     

超声速翼—身组合体Euler方程数值解

本文采用MacCormack二步显式差分格式,用空向推进法给出了超声速翼-身组合体Euler方程的数值解。应用薄翼假设生成简易的多区网格。采用Kentzer方案处理边界。计算中通过对横流速度型的修正改善了背风面压力分布特性,解决了推进中断问题,同时模拟出背风涡。     

压气机/风扇叶片自动优化设计的研究现状和关键技术

 着重介绍和分析了近年来快速发展的压气机叶片自动优化设计方法现状、应用特点以及关键技术。通过分析指出：自动优化设计方法设计效率高,适用于高气动性能指标压气机/风扇设计;为缩短优化时间和提高搜索效率,优化算法应首选并行遗传算法;今后在叶片参数化方法、加速流场计算方法以及提高遗传算法的搜索效率方面需进行进一步研究。