基于非结构四边形网格的WENO格式

基于非结构四边形网格发展求解双曲守恒律的三阶加权基本无振荡（WENO）格式.针对任意非结构四边形网格选取重构模板,并给出基于线性多项式的三阶线性重构.但对于一般的非结构四边形网格,会出现非常大的线性权和负权,使得非线性重构的WENO格式对光滑问题也不稳定.本文给出一个处理非常大的线性权的优化重构方法,对优化后得到的负线性权采用分裂方法进行处理.对于非线性权,提出一种考虑局部网格和物理量间断的新光滑度量因子.采用优化重构方法和新的非线性权,当前的三阶WENO格式在质量很差的网格上也具有很好的稳定性.理论的三阶精度在数值精度测试算例中得到验证,同时一范数和无穷范数的误差绝对值不依赖于网格质量;具有强间断的数值结果证明了当前格式的有效性.     

求解辐射传递的非结构混合有限体积／有限元法

本文给了一种适用于任意非结构网格的有限体积/有限元法的混合算法用于求解多维半透明吸收、发射、散射性灰矩形介质内的辐射传递.该方法使用有限元法进行角度离散,有限体积法进行空间离散.与基于辐射传递离散坐标方程的方法不同的是,该方法在迭代求解的过程中,针对每一个空间体元,所有角度方向的辐射强度同时耦合求出.通过两个算例验证了该解法的正确性.     

求解多维欧拉方程的二阶非结构网格混合旋转Riemann求解器

将基于旋转近似Riemann求解器的二阶精度迎风型有限体积方法推广到非结构网格,采用基于网格中心的有限体积法,梯度的计算采用基于节点的方法引入更多的控制体模板,限制器的构造采用与非结构化网格相适应的形式.在求解Riemann问题时,沿具有一定物理意义的两个迎风方向,即控制体界面两侧速度差矢量方向及与之正交的方向.能够完全消除基于Riemann求解器的通量差分裂格式存在的激波不稳定或红斑现象.为减小计算量,采用HLL和Roe FDS混合旋转格式.     

非结构混合网格上的NS方程求解方法

提出了一套较为通用的,完全自动化的非结构混合网格生成方法.在物面粘性作用区,采用一种改进的推进层方法生成三棱柱形和金字塔形网格;在其他流动区域采用阵面推进方法生成四面体网格.采用一种改进精度的格心有限体积法对三维NS方程进行了求解,在加速收敛措施方面,提出了一种新的当地时间步长取定方法来减小质量较差的网格单元对流场计算稳定性和收敛速度的不利影响.以M6机翼和DLR/F4翼身组合体外形的粘性流场作为数值算例,验证了上述网格生成和流场求解方法的正确性和实用性.     

基于非结构网格的不可压流动耦合求解器

采用非结构化网格有限容积法求解了不可压N-S方程组,对流项采用GAMMA格式,扩散项采用二阶中心差分格式建立离散方程,用SOAR算法处理压力与速度的耦合关系,得到了一种求解不可压N-S方程的非结构网格耦合求解器。通过方腔顶盖驱动流、后台阶绕流以及方腔自然对流等几个典型的算例,考察了求解器的计算精度及收敛特性,并与SIMPLE算法进行了比较,结果表明该求解器是有效可行的。     

多重网格法在非结构网格中的应用

将多重网格法运用于非结构网格.网格是通过聚合法得到的,网格之间是相互关联的.方程的求解采用Jamson的有限体积法.给出了二维、三维情况的数值算例.     

二维非结构网格生成及Euler方程计算的方法研究

用离翼型表面最小作为阵面推进法中的参数选择依据。生成二维问题的非结构网格。这种方法戏了传统的背景网格观念,直接提供网格生成过程中所需的背景信息。在求解Ｅｌｕｅｒ方程时,用格心格式的有限体积法作空间离散,用四步Ｒｕｎｇｅ－Ｋｕｔｔａ作时间推进,采用不同的加速收敛措施获得定常流动。提出了两这界条件的构造办法,并 不同边界条件对结果的影响。     

二维非结构自适应多重网格的Euler方程解

将原来在结构网格中提出的AUSM⁺(Advection Upstream Splitting Method)格式,推广应用于二维非结构网格中。利用格林-高斯公式对控制体内的变量进行线性重构,获得空间二阶精度,并采用Barth型限制器以抑制数值解的振荡。为提高计算效率,采用了自适应多重网格法。在生成非结构网格时,采用了一种新颖的数据结构——数组链接表。最后给出了NACA0012翼型和带10%鼓包的管道流动的几个Euler方程的算例,并进行了分析比较。     

自然对流换热的高精度非结构网格有限体积法

用非结构网格有限体积法求解自然对流换热时,传统的对流项离散格式难以兼顾数值精度与计算效率,我们发展了一种耦合高精度格式的延迟修正方法,用于对流项的离散.高Re数下方腔驱动流数值计算验证了该方法具有较高的计算精度和较好的稳定性.Boussinesq流体的自然对流换热数值模拟,表明该方法能有效克服高Ra数时数值计算发散,可准确捕捉自然对流换热问题中不同偏心率下的等温线和流线分布特征.     

无结构网格有限体积方法及在热对流中的应用研究

首次将无结构三角网格的有限体积方法和压强连接半隐式算法相结合,用于求解非平行壁管道中的热对流问题。并由此分析了化学汽相淀积薄膜生长的均匀性问题。计算结果对于分析一类管道中热和动量输运现象均有普遍指导意义。     

基于WENO重构的高阶移动网格动理学格式

提出一种基于WENO重构的高阶(至少三阶)移动网格动理学格式.利用流体力学方程的积分形式得到移动网格上离散格式,再利用自适应移动网格方法移动网格,进而得到网格速度,利用WENO重构得到高阶插值多项式,最后使用时间方向上精确的动理学数值方法构造数值通量,得到移动网格单元上新的物理量.数值实验表明这种格式同时具有高精度、高分辨率的特点.     

Van der Waals状态方程可压缩多流体流动的PPM方法

基于流体体积分数的混合型多流体数值模型,将Piecewise Parabolic Method（PPM）方法应用于可压缩多流体流动的数值模拟,采用双波近似求解多流体van der Waals状态方程的Riemann问题．模拟高密度比且含有激波的可压缩多流体流动,典型的纯界面平移问题模拟结果表明,在接触间断的界面附近,压力和速度没有任何的振荡且界面数值耗散都被控制在2—3个网格之内;一维和二维算例表明,该数值方法可以有效地处理接触间断、激波和多维滑移线等物理问题,并能够比其它多流体数值方法更精细地模拟多流体交界面．     

A Parameter-Free Generalized Moment Limiter for High-Order Methods on Unstructured Grids

A parameter-free limiting technique is developed for high-order unstruc- tured-grid methods to capture discontinuities when solving hyperbolic conservation laws. The technique is based on a "troubled-cell" approach, in which cells requiring limiting are first marked, and then a limiter is applied to these marked cells. A parameter-free accuracy-preserving TVD marker based on the cell-averaged solutions and solution derivatives in a local stencil is compared to several other markers in the literature in identifying "troubled cells". This marker is shown to be reliable and efficient to consistently mark the discontinuities. Then a compact high-order hierarchical moment limiter is developed for arbitrary unstructured grids. The limiter preserves a degree $p$ polynomial on an arbitrary mesh. As a result, the solution accuracy near smooth local extrema is preserved. Numerical results for the high-order spectral difference methods are provided to illustrate the accuracy, effectiveness, and robustness of the present limiting technique.     

一维多介质可压缩流动数值方法

应用高精度界面追踪方法计算一般状态方程的多介质可压缩流动问题;应用LevelSet技术捕捉界面位置,在界面附近采用守恒数值离散,用双波近似求解一般状态方程Riemann问题,并采用统一高阶PPM格式进行内点和交界面点的计算.一维算例表明,该方法对于光滑区域以及多介质交界面具有二阶精度,能准确地模拟交界面的位置,交界面计算无数值振荡和数值耗散,并能处理一般状态方程的多介质可压缩流动问题.     

三维非结构网格Euler方程的LU-SGS算法及其改进

将LU-SGS隐式时间推进格式运用到非结构网格Euler方程的求解中,并对传统LU-SGS格式进行改进,结合网格重排序,发展了一套效率更高的三维Euler方程求解器.以M6机翼及超临界LANN机翼的跨音速无粘流场为算例,将改进LU-SGS格式与四步龙格-库塔显式格式及传统LU-SGS格式进行了比较.计算结果表明:所有格式的计算结果与实验结果都符合很好;传统LU-SGS格式计算效率为显式格式的3倍多,而改进LU-SGS格式计算效率为显式格式的7倍多.     

基于Riemann解的二维流体力学Lagrange有限点无网格方法

在高维流体力学计算中,对于多介质大变形等一类问题,采用有网格方法常遇到较大的困难.针对二维问题,研究了一种无网格方法——Lagrange有限点方法:在求解区域上设置适当的离散点集,视其中每一点为流体力学Lagrange点;对于点集的任一点,确定邻点集合,并基于该点同邻点集合的联系,应用Godunov方法将流体力学Lagrange方程进行离散;考虑到算法的稳健性,方法中可设置较多邻点并采用最小二乘法.将该方法应用于典型的数值算例,取得了良好效果.     

可压缩非守恒型两相流模型的GKS方法

提出一种求解BN非守恒型两相流模型的方法.从微观角度出发,构造与BN模型相匹配的Gas-kineticscheme(GKS),通过对粒子速度分布函数积分求矩得到通量,把非守恒项直接包含到数值通量的演化和构造中,较好地处理了方程中包含的非守恒项.数值试验说明本方法的有效性.     

双曲型守恒律的一种五阶半离散中心迎风格式

给出一种求解双曲型守恒律的五阶半离散中心迎风格式.对一维问题,该格式以五阶中心WENO重构为基础;对二维问题,用逐维计算的方法将五阶中心WENO重构进行推广.时间方向的离散采用Runge-Kutta方法.格式保持了中心差分格式简单的优点,即不用求解Riemann问题,避免进行特征分解.用该格式对一维和二维Euler方程进行数值试验,结果表明该格式是高精度、高分辨率的.