复杂无粘流场数值模拟的矩形/三角形混合网格技术

建立了一套模拟复杂无粘流场的矩形／三角形混合网格技术，其中三角形仅限于物面附近，发挥非结构网格的几何灵活性，以少量的网格模拟复杂外型；同时在以外的区域采用矩形结构网格，发挥矩形网格计算简单快速的优势，有效地克服全非结构网格计算方法需要较大内存量和较长ＣＰＵ时间的不足．混合网格系统由修正的四分树方法生成．将ＮＮＤ有限差分与ＮＮＤ有限体积格式有机地融合于混合网格计算，消除了全矩形网格模拟曲边界的台阶效应，同时保证了网格间的通量守恒．数值实验表明本方法在模拟复杂无粘流场方面的灵活性和高效性．     

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

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

一个高分辨率的矢通量分裂—TVD杂交新格式

本文提出了一个新的杂交格式,它将Steger-Warming的矢通量分裂与Harten的TVD格式紧密结合在一起,构造了一个高分辨率的新格式,用于计算跨声速流场和捕获激波。典型的定常跨声速叶栅流算例表明:当Courant数取6～100时,一般在90步内残差的二范数下降三个数量级,这样的收敛率要比Beam-Warming格式快得多;观察残差的收敛历史发现:收敛曲线并无大的波动;分析激波附近的数值结果,没出现“低亏,过跳”、伪振荡现象,叶盆和叶背面上激波前或激波后的参数无波动;在60×15网格下捕获的激波过渡区不超过2个网格,表明了该格式具有较高的分辨率,能在不人为附加耗散项的条件下给出高质量无数值波动的激波流场解。     

HLLE++格式在高马赫数空腔流动模拟中的应用

空腔流动存在剪切层运动、涡脱落与破裂,以及激波与激波、激波与剪切层、激波与膨胀波和激波/涡/剪切层相互干扰等现象,流动非常复杂,特别是高马赫数(M>2)时,剪切层和激波更强,激波与激波干扰更严重,对数值格式的要求更高,既需要格式耗散小,对分离涡等有很高的模拟精度,又需要格式在激波附近具有较大的耗散,可以很好地捕捉激波,防止非物理解的出现。Roe和HLLC等近似Riemann解格式在高马赫数强激波处可能会出现红玉现象,而HLLE++格式大大改善了这种缺陷,在捕捉高超声速激波时避免了红玉现象的发生,同时还保持在光滑区域的低数值耗散特性。本文在结构网格下HLLE++格式的基础上,通过改进激波探测的求解,建立了基于非结构混合网格的HLLE++计算方法,通过无粘斜坡算例,验证了HLLE++格式模拟高马赫数流动的能力,并应用于高马赫数空腔流动的数值模拟,开展了网格和湍流模型影响研究,验证了方法模拟高马赫数空腔流动的可靠性和有效性。     

关于无振荡、无自由参数有限元格式的研究

利用双曲守恒律方程的Ｔａｙｌｏｒ弱解表达式，建立了有限元法修正方程，选择合适的展开式系数能得到一系列数值格式．通过稳定性分析研究了格式的稳定性、色散误差与有限元修正方程导数项系数之间的关系，该关系与差分法的ＮＮＤ格式一致．在选定格式下，通过ＣＦＬ数可控制有限元离散解的振荡而使格式不含自由参数．最后，用数值算例验证了这一关系，并在二、三维欧拉方程作了推广应用．     

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

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

高超声速层流尾迹的数值模拟

本文利用无波动、无自由参数、耗散的差分格式(NND格式),通过求解NS方程,数值模拟了高超声速层流尾迹的流动,清晰地给出了主激波、拐角膨胀波、迹激波及自由剪切层,所得流场物理量的分布与实验结果甚为一致。计算发现了底部迴流区由起始向定常的发展中,在瞬时流线图上经历了极限环形成、胀大、缩小、再胀大最后消失的演变过程。     

网格与高精度差分计算问题

研究NS方程差分求解时来流雷诺数、计算格式精度和计算网格之间的关系．给出了判定空间三个方向上的粘性贡献在给定雷诺数、格式精度和网格下是否能够正确计入的估计方法．指出在NS方程的二阶差分方法的数值模拟中,由于物面法向采用了压缩网格技术,物面附近的网格间距很小,该方向上的粘性贡献可被计入．但是如果流向和周向的网格较粗,相应的差分方程中的粘性贡献可能落入截断误差相同的量级,因此在精度上等于仍是求解略去流向和周向粘性项的薄层近似方程．指出,高阶精度的差分计算格式,可以避免对网格要求苛刻的困难．并进一步讨论了建立高阶精度格式的问题,提出了建立高阶精度格式应该满足的原则：耗散控制原则以及色散控制原则．为了避免激波附近可能出现的微小非物理振荡,建议发展混合高阶精度格式,即在激波区,采用网格自适应的NND格式,在激波以外的区域,采用按上述原则发展的高阶格式．     

跨音速压气机级的三维周期性非定常流动计算

对跨音速压气机级动静叶排相干形成的三维非定常流场进行了数值研究，利用时间推进ＬＵ－ＳＧＳ稳式迭代法求解三维非定常欧拉方程，对流项采用高分辨率ＮＮＤ格式离散。对某压气机第一级动静叶排相干非定常流场的计算结果表明，本文方法不仅在收敛速度上明显地优于一般显式方法，而且保持了流场中激波的高分辨率，适于推广到计算量巨大的的多级轴流压气机三维非定常流场的数值分析问题     

近水面水下爆炸的数值研究

近水面水下爆炸后形成的冲击波很快与空气一水的自由界面相互作用，产生折射激波和反射稀疏波，这些波在形成初期是规则的，但很快转变为不规则波。伴随冲击波与自由界面的相互作用，水面要上升直至被炸开。这些复杂的物理现象对数值计算提出了很高的要求。本文采用近期发展起来的位标函数方法以及高精度的ＮＮＤ格式来数值模拟冲击波、自由界面的运动及其相互作用，获得了成功，近水面水下爆炸时的一系列复杂的物理现象均得到了合理的再现。与现有的跟踪自由界面运动的方法相比较，本文方法更简洁，更实用。     

Unsteady oblique shock waves

A spatially and temporally local analysis is provided for unsteady, oblique shock waves, in which the flow is assumed to be two-dimensional or axisymmetric. Three unsteady parameters, in a laboratory frame, are viewed as the known independent variables. These are the upstream Mach number, the shock Mach number, and the angle of the shock relative to the instantaneous upstream velocity. Other steady and unsteady parameters, such as the velocity turn angles and downstream Mach numbers, are evaluated in closed form, in terms of these three quantities. Trends are assessed, and a sensitivity analysis is provided. It is suggested that the theory may find application in converting a shock capturing algorithm, at an early time during the computational process, into a shock fitting algorithm. Received 30 April 1999 / Accepted 29 November 1999     

Numerical solution of the Navier-Stokes equations using a finite element method

A finite element algorithm for solving the Navier-Stokes equations is presented for the analysis of high-speed viscous flows. The algorithm uses triangular elements. The unsteady equations are integrated to steady state with a Runge-Kutta time-marching scheme. A postprocessing artificial dissipation term is introduced to stabilize the computations and to dampen dissipation errors. Numerical results are compared with the calculation of uniform flow on a rectangular region which encounters an embedded oblique shock. A shock/turbulent boundary layer problem is also solved and results are compared with experimental data. It is shown that the postprocessing smoothing term and boundary conditions similar to the finite difference method work well in the present numerical studies.     

从流场干涉图定量计算冲击波传播速度和压力

研究了含冲击波流场干涉图的图像处理方法。首先利用快速条纹细化算法细化干涉条纹,并利用条纹追踪算法提取冲击波波阵面,然后根据配准算法把序列冲击波波阵面配准在同一幅画面上。根据上述算法,在PC-VISION100图像处理机上开发了应用软件。实验结果表明,该方法不但可用于从流场干涉图中提取冲击波波阵面,更重要的是可用于从含冲击波流场干涉图中定量计算冲击波的传播速度和压力分布。     

Modeling wall effects in a micro-scale shock tube using hybrid MD–DSMC algorithm

Wall effects in a micro-scale shock tube are investigated using the Direct Simulation Monte Carlo method as well as a hybrid Molecular Dynamics–Direct Simulation Monte Carlo algorithm. In the Direct Simulation Monte Carlo simulations, the Cercignani–Lampis–Lord model of gas–surface interactions is employed to incorporate the wall effects, and it is shown that the shock attenuation is significantly affected by the choice of the values of tangential momentum accommodation coefficient. A loosely coupled Molecular Dynamics–Direct Simulation Monte Carlo approach is then employed to demonstrate incomplete accommodation in micro-scale shock tube flows. This approach uses fixed values of the accommodation coefficients in the gas–surface interaction model, with their values determined from a separate dynamically similar Molecular Dynamics simulation. Finally, a completely coupled Molecular Dynamics–Direct Simulation Monte Carlo algorithm is used, wherein the bulk of the flow is modeled using Direct Simulation Monte Carlo, while the interaction of gas molecules with the shock tube walls is modeled using Molecular Dynamics. The two regions are separate and coupled both ways using buffer zones and a bootstrap coupling algorithm that accounts for the mismatch of the number of molecules in both regions. It is shown that the hybrid method captures the effect of local properties that cannot be captured using a single value of accommodation coefficient for the entire domain.     

Focusing and Scattering of Plane Shock Waves at an Interface between Anisotropic Elastic Media

The paper is focused on the problem of constructing evolving fronts of quasilongitudinal and quasitransverse shock waves formed by incidence of an initial plane shock wave on a curvilinear interface between elastic transverse isotropic media with different physical properties. The parameter continuation method and the Newton algorithm are used to solve nonlinear Snell's equations. A method for calculating discontinuities of field functions is proposed. Shockwave scattering and focusing as a particular case of bifurcation of shock fronts and formation of caustics are considered. A numerical example is given.     

Unsteady shock‐fitting for unstructured grids

An unstructured, shock‐fitting algorithm, originally developed to simulate steady flows, has being further developed to make it capable of dealing with unsteady flows. The present paper discusses and analyses the additional features required to extend to unsteady flows, the steady algorithm. The properties of the unsteady version of this novel, unstructured shock‐fitting technique, are tested by reference to the inviscid interaction between a vortex and a planar shock: a comparative assessment of shock‐capturing and shock‐fitting is made for the same test problem. Copyright © 2015 John Wiley & Sons, Ltd.     

Performance Improvement of a Supercritical Airfoil by a Multi-point Optimized Shock Control Channel

A shock control channel (SCC) is a flow control method introduced here to control the shock wave/boundarylayer interaction (SWBLI) in order to reduce the resulting wave drag in transonic flows. An SCC transfers an appropriate amount of mass and momentum from downstream of the shock wave location to its upstream to decrease the pressure gradient across the shock wave and as a result the shock-wave strength is reduced. Here, a multi-point optimization method under a constant-lift-coefficient constraint is used to find the optimum design of the SCC. This flow control method is implemented on a RAE-2822 supercritical airfoil for a wide range of off-design transonic Mach numbers. The RANS flow equations are solved using Roe’s averages scheme and a gradient-based adjoint algorithm is used to find the optimum location and shape of the SCC. The solver is validated against experimental works studying effect of cavities in transonic airfoils. It is shown that the application of an SCC improves the average aerodynamic efficiency in a range of off-design conditions by 13.2% in comparison with the original airfoil. The SCC is shown to be an effective tool also for higher angle of attack at transonic flows. We have also studied the SWBLI and how the optimization algorithm makes the flow wave structure and interactions of the shock wave with the boundary layer favorable.     

Simulation of supersonic and hypersonic flows

A simple convection algorithm for simulation of time-dependent supersonic and hypersonic flows of a perfect but viscous gas is described. The algorithm is based on conservation and convection of mass, momentum and energy in a grid of rectangular cells. Examples are given for starting flow in a shock tube and oblique shocks generated by a wedge at Mach numbers up to 30·4. Good comparisons are achieved with well-known perfect gas flows.     

圆柱形汇聚激波诱导 Richtmyer-Meshkov不稳定的 SPH 模拟

汇聚激波诱导不同物质界面的Richtmyer-Meshkov(RM)不稳定现象在惯性约束核聚变领域有重要的学术意义和工程背景.基于网格离散的宏观流体力学方法由于数值扩散问题往往需要高阶精度算法才能准确追踪界面演化,且对大变形和破碎合并等复杂界面追踪也极为困难.光滑粒子流体动力学(smoothed particlehydrodynamics,SPH)方法采用纯拉格朗日算法,可以有效克服上述难点.但经典SPH算法需采用人工黏性处理强间断,在激波间断处往往会出现严重的非物理振荡,对于涉及强冲击不稳定性问题,很难达到理想的模拟效果.本文采用基于HLL黎曼求解器的SPH算法,实现了对强激波和大密度比物质界面的有效分辨和追踪.一维数值校核证明了代码的可靠性、健壮性,并进一步模拟了二维圆柱形汇聚冲击波冲击四边形轻/重气界面诱导的RM不稳定性问题,与已有实验结果进行了对比,发现模拟结果与实验结果吻合.通过分析界面演化过程中的密度及压力变化,发现本文所采用的方法可准确地追踪激波与界面作用的复杂界面和波系演化规律.研究结果为进一步理解和解释汇聚冲击条件下的RM不稳定性机理奠定了基础.      

On the performance of adaptive mesh refinement computation

A finite difference scheme for the unsteady Euler equations using an adaptive mesh refinement (AMR) algorithm was applied to the time-dependent flowfield of shock diffraction problems. The effectiveness of the AMR algorithm was evaluated against a uniform mesh algorithm. Computational results showed that to obtain solutions with equivalent resolution, the AMR algorithm requires much less processing time, when compared with a uniform mesh algorithm.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1.0 and the AMS fonts, developed by the American Mathematical Society.