二维双曲型守恒律的一种五阶松弛格式

松弛格式是Jin和Xin提出的无振荡有限差分方法,其主要思想是将守恒律转化为松弛方程组进行求解.本文用逐维五阶WENO重构和显隐式Runge-Kutta方法对松弛方程组的空间和时间进行离散,得到了一种求解二维双曲型守恒律五阶松弛格式.所得格式保持了松弛格式简单的优点,不用求解Riemann问题和计算通量函数的雅可比矩阵.通过二维Burgers方程和二维浅水方程的数值算例验证了格式的有效性.     

双同守恒律方程的加权本质无振荡格式新进展

近几年,在计算流体力学中,高精度、高分辨率的加权本质无振荡(weighted essentially non-oscillatory , WENO)格式得到很大的发展.WENO格式的主要思想是通过低阶的数值流通量的凸组合重构得到高阶的逼近,并且在间断附近具有本质无振荡的性质.本文综合介绍了双曲守恒律方程的有限差分和有限体积迎风型WENO,中心WENO,紧致中心WENO以及优化的WENO格式等,讨论了负权的处理和多维问题的解决方法.最后,通过一些算例证明WENO格式的高精度,本质无振荡的性质.图6参40      

求解双曲型守恒律的半离散中心差分格式

给出了一种求解双曲型守恒律的三阶半离散中心差分格式。该格式以一种推广的三阶重构为基础,同时考虑了波传播的局部速度。格式的构造方法是利用重构,先计算非一致交错网格上的均值,再将该网格均值投影回原来的非交错网格,得到新的全离散中心差分格式,该格式有半离散形式。本文半离散格式保持了中心差分格式简单的优点,即不需用R iemann解算器,避免了进行特征解耦。它具有守恒形式,数值通量满足相容性条件。数值试验结果表明该格式是高精度、高分辨率的。     

一种改进的中心-WENO混合格式

基于中心差分与WENO格式混合可以改善WENO格式耗散特性的思想,在理论推导的基础上,给出了一种用于激波捕捉计算的守恒型中心-WENO混合格式,该混合格式可视为三阶WENO格式和二阶中心差分格式的加权平均。在数值研究现有加权函数的基础上,给出了适用于该混合格式的加权函数,使其能够自适应地调整数值耗散以捕捉激波间断。数值结果表明:与三阶WENO格式相比,混合格式HY3_4能够降低数值耗散,更陡峭地捕捉间断,对复杂流场结构具有较高的分辨率;混合格式HY3_5对于包含高压比激波间断流场结构,能给出无振荡、低耗散的结果。     

基于非结构任意多边形网格的滑移面计算技术

基于任意多边形网格管理体系,针对流体多介质问题的数值模拟,发展了拉氏方法滑移面计算技术.文章给出了滑移线设置的数据结构,滑移线上主从点速度与位置的计算格式,及节点滑移后引起界面上点、相关网格邻域关系变化的算法.该滑移计算技术避免了传统算法中由于以模拟法(重叠或分离网格)代替直接法(拼接网格)而造成几何守恒律被破坏的缺陷.数值例子验证了该算法的可行性,体现了算法无缝连接的特点.     

适用于高速层流翼型的计算网格研究

高速层流翼型计算的敏感性在于转捩位置的判定与后缘压缩流的捕捉,其气动力系数和流场计算精度受计算网格影响很大。本文针对高速层流翼型流场计算的特点,考虑计算网格的关键影响,提出了影响近壁网格分布和远场网格分布的关键参数,并通过优化研究找出了能显著提高计算精度的网格分布要求。计算结果表明,对基于SST k-ω的γ-Re_(θt)转捩模型RANS数值求解方式,优化网格的网格量(9.6万)在相对于标准网格的网格量(8万)仅增加20%的情况下,显著提高了高速层流翼型LRN1015等的气动特性数值模拟精度,有利于层流翼型的设计与分析。     

一种改进混合迎风格式在翼型流场激波捕捉中的应用

波阻是飞行器超音速飞行的关键设计因素,精确捕捉激波在流场中的位置,是数值模拟含激波流场和精确计算波阻的一个重要研究内容.本文基于网格节点有限体积空间离散方法,采用AUSM格式与FVS格式的混合格式(MAUSM方法)计算对流通量,从而抑制在数值模拟流场出现的激波处振荡和过冲现象,确保AUSM准确捕获接触间断的特性和FVS格式捕捉激波的能力.本文使用MAUSM方法分别计算了在跨声速和超声速条件下的NACA0012翼型流场,并与中心差分格式的计算结果进行比较.结果表明,对于存在激波的翼型流场,MAUSM方法是有效的.     

求解双曲守恒律方程的WENO型熵相容格式

通过在单元交界面处进行高阶WENO重构,得到了一种求解双曲型守恒律方程的WENO型熵相容格式。用该格式对一维Burgers方程和Euler方程进行数值模拟,结果表明,该格式具有高精度、基本无振荡性等特点。     

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

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

一种简单的精确捕捉接触间断的黎曼求解器

基于Godunov型数值格式的有限体积法是求解双曲型守恒律系统的主流方法,其中用来计算界面数值通量的黎曼求解器在很大程度上决定了数值格式在计算中的表现。单波的Rusanov求解器和双波的HLL求解器具有简单、高效和鲁棒性好等优点,但是在捕捉接触间断时耗散太大。全波的HLLC格式能够精确捕捉接触间断,但是在计算中出现的激波不稳定现象限制了其在高马赫数流动问题中的应用。本文利用双曲正切函数和五阶WENO格式来重构界面两侧的密度值,并且结合边界变差下降算法来减小Rusanov格式耗散项中的密度差,从而提高格式对于接触间断的分辨率。研究表明,相比于全波的HLLC求解器,本文构造的黎曼求解器不仅具有更高的接触分辨率,而且还具有更好的激波稳定性。     

Numerical simulation of free surface flows with the level set method using an extremely high-order accuracy WENO advection scheme

Three-dimensional dynamic gas–liquid flow simulations that accurately track the phase interface are numerically challenging. This article presents a numerical study of the performance of the level-set phase–interface tracking method when combined with extremely high order (7th to 11th) weighted essentially non-oscillatory (WENO) advection schemes for gas–liquid free surface flows. Comparisons between simulation results and prior benchmark results suggest that such a combination of methods can be satisfactorily applied to the level-set and Navier-Stokes equations for free surface flow simulations when volume conservation is enforced at every time step, and minor numerical oscillations are suppressed through use of an artificial viscosity term. In particular, simulations of solid body rotation, the unsteady flow following an ideal dam break, tank sloshing, and the rise of a single bubble all agree with analytical or experimental results to within ± 3.12% when the level-set method is combined with an 11th order WENO scheme. Furthermore, use of an 11th order WENO advection scheme actually has a computational cost advantage because, for the same accuracy, it can be used on a coarser grid when compared with a more-common second-order advection scheme; computational savings of up to 87% are possible.     

A high‐resolution scheme for compressible multicomponent flows with shock waves

A simple methodology for a high‐resolution scheme to be applied to compressible multicomponent flows with shock waves is investigated. The method is intended for use with direct numerical simulation or large eddy simulation of compressible multicomponent flows. The method dynamically adds non‐linear artificial diffusivity locally in space to capture different types of discontinuities such as a shock wave, contact surface or material interface while a high‐order compact differencing scheme resolves a broad range of scales in flows. The method is successfully applied to several one‐dimensional and two‐dimensional compressible multicomponent flow problems with shock waves. The results are in good agreement with experiments and earlier computations qualitatively and quantitatively. The method captures unsteady shock and material discontinuities without significant spurious oscillations if initial start‐up errors are properly avoided. Comparisons between the present numerical scheme and high‐order weighted essentially non‐oscillatory (WENO) schemes illustrate the advantage of the present method for resolving a broad range of scales of turbulence while capturing shock waves and material interfaces. Also the present method is expected to require less computational cost than popular high‐order upwind‐biased schemes such as WENO schemes. The mass conservation for each species is satisfied due to the strong conservation form of governing equations employed in the method. Copyright © 2010 John Wiley & Sons, Ltd.     

加权型紧致格式与加权本质无波动格式的比较

线性紧致格式和加权本质无波动格式是两种典型的高阶精度数值格式,它们各有优缺点.线性紧致格式在具有高阶精度的同时,格式的分辨率也比较高,耗散低,是计算多尺度流场结构的较好格式,但是不能计算具有强激波的流场.加权本质无波动格式是一种高阶精度捕捉激波格式,鲁棒性好,但耗散比较高,分辨率也不理想.近年来,在莱勒的线性紧致格式基础上,采用加权本质无波动格式捕捉激波思想,发展了一系列加权型紧致格式.本文较全面地比较了加权型紧致格式和加权本质无波动格式,包括构造方法、鲁棒性、分辨率、耗散特性、收敛特性以及并行计算效率.结果表明,现有的加权型紧致格式基本保持了加权本质无波动格式的性质,对于气动力等宏观量的计算,比加权本质无波动格式没有明显的优势.      

An adaptive multilevel wavelet framework for scale‐selective WENO reconstruction schemes

We put forth a dynamic computing framework for scale‐selective adaptation of weighted essential nonoscillatory (WENO) schemes for the simulation of hyperbolic conservation laws exhibiting strong discontinuities. A multilevel wavelet‐based multiresolution procedure, embedded in a conservative finite volume formulation, is used for a twofold purpose. (i) a dynamic grid adaptation of the solution field for redistributing grid points optimally (in some sense) according to the underlying flow structures, and (ii) a dynamic minimization of the in built artificial dissipation of WENO schemes. Taking advantage of the structure detection properties of this multiresolution algorithm, the nonlinear weights of the conventional WENO implementation are selectively modified to ensure lower dissipation in smoother areas. This modification is implemented through a linear transition from the fifth‐order upwind stencil at the coarsest regions of the adaptive grid to a fully nonlinear fifth‐order WENO scheme at areas of high irregularity. Therefore, our computing algorithm consists of a dynamic grid adaptation strategy, a scale‐selective state reconstruction, a conservative flux calculation, and a total variation diminishing Runge‐Kutta scheme for time advancement. Results are presented for canonical examples drawn from the inviscid Burgers, shallow water, Euler, and magnetohydrodynamic equations. Our findings represent a novel direction for providing a scale‐selective dissipation process without a compromise on shock capturing behavior for conservation laws, which would be a strong contender for dynamic implicit large eddy simulation approaches.     

多流体界面不稳定性的守恒和非守恒高精度数值模拟方法

对多流体界面问题守恒和非守恒格式(M)WENO重构方法进行探讨,采用虚拟流动方法并用Level-Set函数捕捉界面的运动变化。数值模拟结果表明本文的数值方法具有较高的分辨率,并能有效地抑制界面附近的非物理振荡。     

A novel weighting switch function for uniformly high‐order hybrid shock‐capturing schemes

Hybrid schemes are very efficient for complex compressible flow simulation. However, for most existing hybrid schemes in literature, empirical problem‐dependent parameters are always needed to detect shock waves and hence greatly decrease the robustness and accuracy of the hybrid scheme. In this paper, based on the nonlinear weights of the weighted essentially non‐oscillatory (WENO) scheme, a novel weighting switch function is proposed. This function approaches 1 with high‐order accuracy in smooth regions and 0 near discontinuities. Then, with the new weighting switch function, a seventh‐order hybrid compact‐reconstruction WENO scheme (HCCS) is developed. The new hybrid scheme uses the same stencil as the fifth‐order WENO scheme, and it has seventh‐order accuracy in smooth regions even at critical points. Numerical tests are presented to demonstrate the accuracy and robustness of both the switch function and HCCS. Comparisons also reveal that HCCS has lower dissipation and less computational cost than the seventh‐order WENO scheme. Copyright © 2016 John Wiley & Sons, Ltd.     

An improved third-order finite difference weighted essentially nonoscillatory scheme for hyperbolic conservation laws

In this article, we present an improved third-order finite difference weighted essentially nonoscillatory (WENO) scheme to promote the order of convergence at critical points for the hyperbolic conservation laws. The improved WENO scheme is an extension of WENO-ZQ scheme. However, the global smoothness indicator has a little different from WENO-ZQ scheme. In this follow-up article, a convex combination of a second-degree polynomial with two linear polynomials in a traditional WENO fashion is used to compute the numerical flux at cell boundary. Although the same three-point information is adopted by the improved third-order WENO scheme, the truncation errors are smaller than some other third-order WENO schemes in L_∞ and L₂ norms. Especially, the convergence order is not declined at critical points, where the first and second derivatives vanish but not the third derivative. At last, the behavior of improved scheme is proved on a variety of one- and two-dimensional standard numerical examples. Numerical results demonstrate that the proposed scheme gives better performance in comparison with other third-order WENO schemes.