适于可压缩多尺度流动的紧致型激波捕捉格式

针对可压缩多尺度流动数值模拟特点,研究一种五阶高分辨率紧致型激波捕捉格式——紧致重构加权基本无振荡(CRWENO)格式。该格式利用非线性权系数将低阶紧致格式加权组合以达到高阶精度。在光滑区域蜕化成具有高分辨率的五阶线性紧致格式,在间断附近则能保持计算稳定无振荡。对CRWENO格式、目前广泛使用的加权基本无振荡(WENO)格式及两格式对应的线性格式(即五阶线性迎风格式和五阶紧致格式)进行数值性能研究,评估非线性权系数对格式耗散及频谱特性的影响。使用一维、二维、三维典型算例进行数值试验,探讨线性/非线性、紧致/非紧致格式在可压缩多尺度流动模拟中的优势和不足。结果表明,CRWENO格式在强压缩性流场模拟中能够稳定地捕捉激波,其紧致特性则改善了非线性格式普遍存在的耗散过大、分辨率较差的问题,使其能够清晰捕捉多尺度流动结构。因此,该格式在可压缩多尺度流动模拟中具有较大优势。     

Monotonicity-Preserving激波捕捉格式在湍流大尺度模拟中的评估

对高阶激波捕捉格式的性能进行了系统的测评,重点分析了Suresh和Huynh(1997)所提出的Monotonicity-Preserving格式的性能.结果表明Monotonicity-Preserving格式的性能显著优于原始WENO(Weighted Essentially Non-Oscillatory) 格式,和改进型WENO格式相当.对格式的分析进一步表明,迎风型的激波捕捉格式在湍流模拟方面的性能都不及高阶中心格式,其原因归结为激波捕捉格式所包含的线性和非线性耗散.因此,改进高阶激波捕捉格式的关键在于同时降低格式的线性耗散和非线性耗散,以提高格式对湍流脉动能量的保持和对小尺度脉动结构的捕捉能力.     

三阶中心无振荡格式性能分析

三阶中心无振荡格式Cn3(Centered nonoscillatory scheme of third order)使用对称模板构造具有三阶精度的插值公式.利用单调区域与精确区域修正原始插值参数,该格式能够获得间断附近无振荡、光滑区域高精度的计算结果.通过一维和二维典型算例,将Cn3格式与三阶和五阶WENO格式(Weighted Essentially Non-Oscillatory schemes)进行比较,重点分析了3种格式的间断分辨率、计算稳定性和数值耗散性.分析可见,Cn3格式能够精确、稳定地捕捉激波和接触间断,同时对光滑流动区域的小尺度流动结构保持较低的耗散,值得进一步研究及推广应用.     

无碰撞激波的理想MHD模拟

采用了WENO格式数值求解一维理想MHD方程组，模拟了行星际无碰撞激波，研究了垂直无碰撞激波与行星际反向磁场结构和高密度等离子体团的相互作用过程，并与粒子模拟的结果进行比对，两者的结果非常类似．模拟结果表明，对涉及无碰撞激波的大部分现象中理想MHD模拟是准确且可行的，同时相对于粒子模拟又有很好的计算效率，便于扩展至二维或三维的情形．     

基于CE/SE法的翼型流场激波捕捉计算

在二维四边形结构网格下,针对时空守恒元和解元(CE/SE, space-time Conservation Element and Solution Element)方法捕捉激波需要双倍空间网格点的问题,重新对守恒元和解元进行了定义,计算点同时包含单元中心和网格节点,并推导得到了新的CE/SE方法的计算公式.以此为基础,结合非敏感克朗数(CNIS, Courant Number Insensitive Scheme)计算格式和当地时间步长方法,对激波翼型流场进行数值模拟,并与原方法及AGARD报告进行对比.结果表明,对原CE/SE方法的改进是有效的,可明显提高激波分辨能力,并且激波前后无明显的数值震荡发生,新方法适合应用于翼型流场中捕捉激波.     

航空发动机一维数值仿真

分析了航空发动机一维气动热力数值仿真模型,此模型能够用于航空发动机整机数值仿真.进行了整机性能数值仿真,求解过程考虑了粘性效应以及气体流动损失对于整机性能的影响;采用了具有TVD性质的高阶精度Godunov格式,此格式具有良好的数值稳定性和激波捕捉能力;应用了时间推进法,交替使用显/隐格式.计算了某型带加力的混排涡扇发动机的速度特性以及节流特性并且与设计数据进行了比较分析,数值仿真结果与设计数据比较吻合.分析了典型工作点的计算过程,模型收敛速度较快.     

笔式草图的增量识别

针对笔式三维建模技术中手势交互的特点,提出了一种二维草图的增量识别方法:从原始笔划的在线识别与线元转换,到约束关系的自动捕捉与施加,到二维组合图元的识别与规整,到特征手势的识别,到最后的三维特征建模.整个识别过程依据了三维特征建模思想和基于上下文感知的智能推理技术,通过分层识别的方法从紊乱、非精确的二维草图中有效提取三维建模的特征信息,有效降低了草图识别过程的复杂度和积累误差.将特征手势遵循隐式映射机制映射到三维空间后,实时完成三维建模,实现了二维手绘草图与三维精确建模的有效集成.该方法已经在自主开发的原型系统PenSketch3D中得以验证,并取得良好效果.      

高效率的特征型紧致WENO混合格式

特征型紧致加权基本无振荡（WENO）混合格式HCW-R结合了迎风紧致格式CS5-P和WENO格式，具有十分优异的分辨率特性。但在求解多维方程组时，HCW-R格式需要求解块状三对角方程组，因而计算代价十分高昂。采用迎风紧致格式CS5-F代替CS5-P，构造了一个新的特征型紧致WENO混合格式HCW-E。由于HCW-E的特殊形式，其可沿迎风方向、由边界处向内推进求解，避免了处理三对角或块状三对角方程组，从而其计算代价与显式格式无异。虽然就分辨率而言，HCW-E稍逊于HCW-R，但前者的计算效率要显著高于后者。因此，当花费相同的计算代价，HCW-E格式可以获得更好的数值结果。一系列求解Euler方程组的数值试验验证了HCW-E的高分辨率特性和相比HCW-R更高的计算效率。HCW-E格式的效率优势在求解高维问题时更为明显。      

基于简单WENO-间断Galerkin的Euler方程自适应计算

为了得到Euler方程的高精度、高分辨率数值解,介绍了间断Galerkin方法、三角形单元上简单WENO限制器的基本原理以及基于自适应网格加密的激波捕捉方法。将简单WENO限制器-间断Galerkin方法应用到曲边四边形单元上,通过单元边界上高斯积分点的坐标来搜索相邻单元从而得到相邻单元的单元编号,实现了基于“问题单元”的局部网格加密自适应计算。对若干典型问题进行编程计算,结果表明,简单WENO限制器可以应用到曲边四边形单元上,且可适用于局部网格加密时具有“悬挂节点”的非结构网格上的激波捕捉。      

高精度多维限制器的性能分析

目前常用的限制器大都是基于一维构造,无法在多维情况下保证物理量的单调特性进而导致非物理振荡.为弥补传统方法的这一构造缺陷,多维限制器(MLP)通过多维修正使单元通量值介于周围相邻单元通量的最大值和最小值之间,在保证求解精度的情况下有效避免了多维振荡.基于一维激波管、无黏涡及激波边界层干扰等算例,对高精度MLP的特性进行了研究分析.结果显示:3阶MLP在连续和间断区域均可有效地避免多维振荡;与高阶WENO(Weighted Essentially Non-Oscillatory)方法相比,3阶MLP不仅算法简单、易于实现,还可显著提高求解的精度、保单调性及收敛性.因此可用于工程及科学研究的复杂流动,具有较好的应用前景.      

高超声速热流高精度数值模拟方法

通过高超声速下钝锥热流的计算,对影响热流计算精度的因素进行了综合研究.获得了壁面法向网格Re数对热流计算的影响规律.采用边界高阶插值,提高了热流计算精度.用Roe和AUSM(Advection Upstream Splitting Method)+格式搭配minmod和混合limiter两种限制器,探讨了空间离散方法及限制器的耗散性对热流结果的影响.应用五阶WENO(Weighted Essentially Non-Oscillatory)方法,对高阶格式计算热流的性能进行了细致研究,最终确立了网格依赖小,热流精度高的计算方法.研究认为:高阶格式不但可以放宽热流计算对网格Re数的限制,同时使得格式对本身计算方法的依赖性也变小.      

改进的二维三阶半离散中心迎风格式

对二维三阶半离散中心迎风格式中的权函数给出了简化改进. 在保持格式精度的基础上, 改进后的权函数在二维情况下具有更加简单直接的结构而且严格非负. 该改进方法得到的格式仍然具有半离散中心迎风格式的优点, 同时保持了重构函数的非振性. 时间离散采用保持强稳定性的三阶Runge-Kutta方法, 并利用四阶Lax-Wendroff (L-W) 格式计算磁流体算例中的磁场散度. 用该修正格式计算了二维磁流体数值算例, 得到高精度无振荡的结果, 验证了此方法的有效性.      

A CLASS OF TWO-STEP TVD MACCORMACK TYPE NUMERICAL SCHEME FOR MHD EQUATIONS

In this paper, a new numerical scheme of Total Variation Diminishing (TVD) Mac-Cormack type for MagnetoHydroDynamic (MHD) equations is proposed by taklng into account of the characteristics such as convergence, stability, resolution. This new scheme is established by solving the MHD equations with a TVD modified MacCormack scheme for the purpose of developing a scheme of quick convergence as well as of TVD property. To show the validation, simplicity and practicability of the scheme for modelling MHD problems, a self-similar Cauchy problem with the discontinuous initial data consisting of constant states, and the collision of two fast MHD shocks, and two-dimensional Orszag and Tang‘s MHD vortex problem are discussed with the numerical results conforming to the existing results obtained by the Roe type TVD, the high-order Godunov scheme, and Weighted Essentially Non-Oscillatory (WENO) scheme. The numerical tests show that this two-step TVD MacCormack numerical scheme for MHD system is of robust operation in the presence of very strong waves, thin shock fronts, thin contact and slip surface discontinuities.     

A 3rd Order WENO GLM-MHD Scheme for Magnetic Reconnection

A new numerical scheme of 3rd order Weighted Essentially Non-Oscillatory (WENO)type for 2.5D mixed GLM-MHD in Cartesian coordinates is proposed. The MHD equations are modified by combining the arguments as by Dellar and Dedner et al to couple the divergence constraint with the evolution equations using a Generalized Lagrange Multiplier (GLM). Moreover, the magnetohydrodynamic part of the GLM-MHD system is still in conservation form. Meanwhile, this method is very easy to add to an existing code since the underlying MHD solver does not have to be modified. To show the validation and capacity of its application to MHD problem modelling,interaction between a magnetosonic shock and a denser cloud and magnetic reconnection problems are used to verify this new MHD code. The numerical tests for 2D Orszag and Tang's MHD vortex,interaction between a magnetosonic shock and a denser cloud and magnetic reconnection problems show that the third order WENO MHD solvers are robust and yield reliable results by the new mixed GLM or the mixed EGLM correction here even if it can not be shown that how the divergence errors are transported as well as damped as done for one dimensional ideal MHD by Dedner et al.      

Application of ADER Scheme in MHD Simulation

The Arbitrary accuracy Derivatives Riemann problem method(ADER) scheme is a new high order numerical scheme based on the concept of finite volume integration,and it is very easy to be extended up to any order of space and time accuracy by using a Taylor time expansion at the cell interface position.So far the approach has been applied successfully to flow mechanics problems.Our objective here is to carry out the extension of multidimensional ADER schemes to multidimensional MHD systems of conservation laws by calculating several MHD problems in one and two dimensions: (ⅰ) Brio-Wu shock tube problem,(ⅱ) Dai-Woodward shock tube problem,(ⅲ) Orszag-Tang MHD vortex problem.The numerical results prove that the ADER scheme possesses the ability to solve MHD problem,remains high order accuracy both in space and time,keeps precise in capturing the shock.Meanwhile,the compared tests show that the ADER scheme can restrain the oscillation and obtain the high order non-oscillatory result.     

Euler方程的分裂型通量分裂双时间步隐式方法

传统隐式方法有格式复杂、计算量大等缺点,在Euler方程的差分离散过程中,利用算子分裂思想,结合通量分裂法、双时间步法等隐式离散方法,构造了一种更简单的分裂型隐式计算方法.通过对典型空气动力学问题的计算,检验了该方法的有效性和可靠性,并对其性能做了具体讨论.该方法具有稳定性好、时间步长约束小等隐式格式的普遍优点,同时具有格式简单、程序易实现等优点;避免了传统隐式方法单步推进时的方程组常规求解及矩阵求逆过程,计算量小;比LU-SGS方法收敛速度快.