分区对接网格在跨流域气体运动论统一算法的应用研究北大核心CSCD

为了模拟跨流域复杂物形绕流或多体流场,需要在位置空间建立能描述整个流场特性的网格系统,相应的数值计算方法也需要有所改进。本文采用分区对接网格处理复杂物形绕流或多体干扰流场,得到能准确描述流场特征的贴体网格系统。在此基础上数值求解考虑转动非平衡影响的Boltzmann-Rykov模型方程,在速度空间应用离散速度坐标法、位置空间采用数值差分NND格式。网格对接面上考虑分布函数的传递,包含网格点信息和网格点上速度的信息。通过计算得到复杂物形的流动特性。相关算例表明网格分区之间的流场参数光滑过渡。与文献及DSMC结果的对比证实本文基于分区对接网格的气体运动论统一算法是可靠的。     

分区对接网格在跨流域气体运动论统一算法的应用研究

为了模拟跨流域复杂物形绕流或多体流场,需要在位置空间建立能描述整个流场特性的网格系统,相应的数值计算方法也需要有所改进。本文采用分区对接网格处理复杂物形绕流或多体干扰流场,得到能准确描述流场特征的贴体网格系统。在此基础上数值求解考虑转动非平衡影响的Boltzmann-Rykov模型方程,在速度空间应用离散速度坐标法、位置空间采用数值差分NND格式。网格对接面上考虑分布函数的传递,包含网格点信息和网格点上速度的信息。通过计算得到复杂物形的流动特性。相关算例表明网格分区之间的流场参数光滑过渡。与文献及DSMC结果的对比证实本文基于分区对接网格的气体运动论统一算法是可靠的。     

桥梁断面气动导纳的数值识别方法研究北大核心CSCD

针对一种气动导纳的数值识别方法进行研究。基于二维不可压缩URANS方法,选用SSTk-ω湍流模型,通过在来流中给定单一频率的竖向谐波速度分量,计算相应的桥梁断面气动力荷载时程,识别桥梁断面的气动导纳。首先考查来流脉动特性在计算域内的自保持能力,随后再对平板和桥梁断面的气动导纳进行识别,所得结果与理论解和试验值相比较,并讨论流场初始化条件的影响。结果表明:足够小的网格尺寸和时间步长是来流脉动不发生明显衰减的必要条件;平板的气动导纳识别结果与Sears函数高度吻合;数值识别的桥梁断面升力气动导纳在低频段与Sears函数一致,在高频段略低,但与试验值较接近;力矩气动导纳与Sears函数有较大差异,但与试验值基本吻合;流场初始化条件对计算效率有影响。     

桥梁断面气动导纳的数值识别方法研究

针对一种气动导纳的数值识别方法进行研究。基于二维不可压缩URANS方法,选用SSTk-ω湍流模型,通过在来流中给定单一频率的竖向谐波速度分量,计算相应的桥梁断面气动力荷载时程,识别桥梁断面的气动导纳。首先考查来流脉动特性在计算域内的自保持能力,随后再对平板和桥梁断面的气动导纳进行识别,所得结果与理论解和试验值相比较,并讨论流场初始化条件的影响。结果表明:足够小的网格尺寸和时间步长是来流脉动不发生明显衰减的必要条件;平板的气动导纳识别结果与Sears函数高度吻合;数值识别的桥梁断面升力气动导纳在低频段与Sears函数一致,在高频段略低,但与试验值较接近;力矩气动导纳与Sears函数有较大差异,但与试验值基本吻合;流场初始化条件对计算效率有影响。     

关于电磁场的自然单元法

有限元法和无网格法在电磁场数值计算中已经得到了广泛的应用,然而有限元法存在前处理网格剖分问题,无网格法存在计算时间长、边界条件和不连续面难处理等问题。针对以上问题,本文提出了利用自然单元法求解电磁场的方法,该方法根据场域中的离散点的信息构造自然邻点插值函数,最后算出场值。通过算例证明了该方法在电磁场计算中的可行性。     

基于人工黏性的二维激波小波多尺度数值计算

基于自适应小波配点法和人工黏性，构造了二维激波问题的小波数值计算格式。利用最细尺度的小波系数构造两组激波定位函数，分别用以控制两个正交方向人工黏性的大小和分布。对三类二维激波问题进行计算和验证，结果表明：自适应小波配点法利用小波阈值滤波删除变化平缓区域大量网格点，而保留变化急剧区域的网格点，比传统方法计算效率高，且分辨率越高，计算效率更高；利用最细层小波系数构建的幂函数形式的激波定位函数能准确判断激波位置和控制X、Y方向人工黏性的大小和分布，从而捕捉不同方向和不同强度的激波。     

一种基于三阶WENO重构的无网格算法研究

提出了基于三阶WENO重构的无网格算法,以期替代传统的无网格线性逼近重构,提高无网格算法的求解精度。基于无网格点云卫星点分布特征,通过沿点云中心点与卫星点连线方向引入局部一维坐标系,并结合虚拟点的设置,构成了在点云中实施三阶WENO重构所要求的模板。算法涉及的卫星点连线中点处的左右状态值,则利用模板上各点的流场值,采用WENO重构方法计算给出。对于模板中设置的虚拟点上的流场值,本文则利用已有的无网格点云信息,基于最近节点流场值插值得到。文中给出了采用上述WENO重构方法,结合Roe的近似Riemann求解器确定通量,并采用三阶TVD Runge-Kutta方法进行时间推进求解Euler方程的具体实施过程。对模型问题的典型流动进行了计算分析,通过比较数值解和精确解,验证了算法获得的数值解能逼近三阶精度。给出的激波管流动和超声速半圆柱绕流算例展示出所发展的算法捕捉激波等间断问题具有更高分辨率,能体现出WENO重构"基本无振荡"的特性。文中最后给出了非定常激波过弯道绕双圆柱流场的数值模拟算例,展示了所发展的算法处理含非定常激波演化的复杂流动问题的效果。     

壁面函数在激波诱导分离流动中的应用

以数值模拟激波-附面层干扰引起的流动分离问题为研究背景,发展了基于有限体积方法的雷诺平均Navier—Stokes（RANS）方程的流场数值模拟方法。利用壁面函数模型得到壁面剪切应力,通过修正壁面粘性通量,构造了一种新的湍流边界处理方法,并将其耦合到RANS方程和SSTk-ω湍流模型的数值求解中;同时,针对激波诱导引起的附面层流动分离问题,提出一种附面层网格加密技术,能够自适应加密分离区内附面层网格,使得在流动分离区域也能够使用壁面函数模型。数值算例表明,壁面函数模型能够降低数值模拟结果对网格的依赖性;同时也验证了壁面函数耦合附面层网格自适应方法,在处理激波诱导引起的附面层流动分离问题时的有效性和准确性。     

桥梁断面气动导纳风场依赖特性的数值研究

基于二维不可压非定常RANS模型,分别计算了平板断面、箱梁断面、双边肋断面在简谐脉动来流和宽频湍流中的非定常气动力,并识别出各自的气动导纳函数。在此基础上,通过考察断面在三种不同脉动幅值来流作用下的气动导纳,研究了不同脉动幅值对各类断面气动导纳的影响。结果表明,在竖向简谐脉动来流下,平板断面的气动导纳与Sears函数吻合良好,证明数值识别气动导纳的可行性;湍流作用下,平板断面识别的导纳函数有着实质性的不同,表现出很强的随机跳跃性。与平板断面类似,箱梁断面的气动导纳接近Sears函数。在不同幅值的简谐脉动来流下,平板断面和流线型箱梁断面识别的气动导纳基本一致,表现出了流线型断面对来流特性的不敏感性。计算结果表明,具有显著分离流的双边肋断面,其气动导纳具有明显的来流风场依赖性,表明钝体断面的气动导纳并不能定性为断面的函数,而是由断面与来流风场特性共同确定。     

三维分块非结构化网格上的可压缩流计算方法

提出了分块非结构化网格生成方法,通过子块划分、子块网格划分、整体网格拓扑连接生成可局部加密的分块非结构化网格,既具备了非结构化网格对复杂区域的处理能力,又避免了分块结构化网格块间信息传递的复杂计算过程.通过将密度修正值引入压力修正方程,发展了非结构化网格上可压缩流动的SIMPLE算法.对孤立翼型及透平叶栅的流动进行了数值计算,压力系数、速度、折转角的数值结果与已有的实验值吻合较好,验证了数值方法的精度.      

压气机转子叶片的气动弹性数值模拟

采用了交替迭代算法,对某压气机转子叶片进行了气动弹性数值模拟.自行开发了基于有限元的结构求解器用于结构动力学求解,引用他人开发的非定常流体求解器用于气动力的求解.结构求解器提供叶片的表面位移给流体求解器以改变流场,流体求解器提供气动载荷给结构求解器来计算叶片的变形,界面处理系统在这两个求解器之间进行信息传递.算例表明,这种交替迭代算法在气动弹性数值模拟中是可行的,可以得到叶片的瞬态响应,从而判断叶片是否发生颤振.      

任意运动变形边界流场测量PIV算法的研究：薄膜涡激振动

提出了一种精确测量任意运动变形边界流场的粒子图像速度场PIV算法,能智能识别强流固耦合问题流场测量中任意运动变形的边界信息,然后生成贴体自适应的图像互相关计算窗口,计算获取运动变形边界附近的速度场数据。这一算法可大大提高任意运动边界附近流场的测量精度,为流固耦合问题的理论分析和数值计算验证提供可靠的高质量实验数据。针对数字合成给定流场的粒子图像序列,采用所发展的PIV算法对粒子图像区域中的运动变形边界进行了精确识别,高质量地复现了原始流场信息。最后,对低速闭式循环水洞中的钝体尾流柔性薄膜涡激振动现象进行了PIV实验测量,获得了柔性薄膜大变形运动状态下的瞬态流场特征。     

Computational modeling for multiphase flows with spacecraft application

Many engineering applications involve interactions between solid, gas and liquid phases under normal or micro-gravity conditions. Numerical simulations of such fluid flows need to track the location and the shape of the fluid interface as part of the solution. The merits and basic characteristics of various approaches for numerical computations of interfacial fluid dynamics are reviewed. The computational challenges include: (i) the algorithmic complexity for handling irregularly shaped moving boundaries that can experience merger and break-up; (ii) resolution refinement techniques to maintain desirable resolution of length scales, in accordance with the evolving fluid dynamics; (iii) data structure needed to support identification of the interface and satisfaction of the physical laws in the bulk fluids as well as around the phase boundaries; and (iv) efficient parallel processing techniques required for practical engineering analysis. The present review focuses on these issues related to the Lagrangian–Eulerian approach, utilizing the immersed boundary method with marker-based tracking, as the main framework for interfacial flow computations on Cartesian grids. Specifically, we offer in-depth discussion of the organization and layout of the mesh systems for both fluid and interface representations, local adaptive refinement on two-dimensional/three-dimensional (2D/3D) Cartesian grids, and multi-level domain decomposition method that utilizes Hilbert space filling curves for parallel processing strategy. The effectiveness of individual components and overall algorithm are presented using various tests such as, binary drop-collision computations to highlight grid adaptation and interface tracking algorithms to handle complex interface behavior, and bubble/droplet placed in a vortex field with various density/viscosity ratios across interfaces to address load balancing and scalability aspects of parallel computing. A time-dependent draining flow problem motivated by spacecraft applications is presented to aid spacecraft design efforts.     

热传导方程有限差分区域分裂：显—隐算法的注记

把Rodrigue和Wolitzer（1995）提出的一种新型显格式应用于热传导方程有限差分区域分裂算法的内边界点,并得到这类显—隐算法的收效性结果。新算法的逼近阶较好,而且较纯新型显式方法可以节省一倍左右的计算量。     

多块多重网格法及其跨声速转子内流并行模拟

为提高压气机内部复杂流动计算的收敛速度,发展了一种多块多重网格上的可压缩流动并行计算流体动力学(CFD)方法,数值模拟了跨声速转子NASA Rotor 35内部流动及气动性能,重点分析了对接分区方式和内界面通信模式对并行性能的影响。结果表明:对接分区方法和多块多重网格内界面数据处理方法能够保证串、并行CFD方法结果基本一致,且与实验数据吻合良好;并行精度基本不受分区数目和通信模式的影响,并行效率随着相对通信负载的增大而降低。所发展的并行计算方法对跨声速转子气动分析和设计具有一定的可靠性,为开发大型流体机械复杂流动的基础并行算法及软件提供了参考价值。     

改进虚拟流体方法在气-水界面捕捉中的应用

提出了一种改进的Ghost方法,通过求解双介质的Riemann问题,将全部解参量外推得到虚拟流体状态参量,使单一的气-水间断问题转化为两个单介质的Riemann问题。流场控制方程求解采用二阶NND格式,LevelSet方程用五阶WENO格式离散,实现了对气-水界面稳定而准确的捕捉。激波与气-水界面作用的一维和二维算例验证了程序的有效性。对气液两相均运用可压模型,计算了单个水滴被马赫数为2的空气激波扫过后的二维变形过程。     

计算地球流体力学若干新进展

本文简要介绍了有关计算地球流体力学问题研究的若干新进展,其中包括显式平方守恒格式,经济算法,辛算子法,泥沙冲积和三角洲发育的数值模拟等,从中也可以看到计算地球流体力学的发展是具有重要理论意义和实用价值的。     

Fluid structure interaction simulation of supersonic parachute inflation by an interface tracking method

An Arbitrary Lagrangian–Eulerian (ALE) approach with interface tracking is developed in this paper to simulate the supersonic parachute inflation. A two-way interaction between a nonlinear finite element method and a finite volume method is accomplished. In order to apply this interface tracking method to problems with instantaneous large deformation and self-contact, a new virtual structure contact method is proposed to leave room for the body-fitted mesh between the contact structural surfaces. In addition, the breakpoint due to the fluid mesh with negative volume is losslessly restarted by the conservative interpolation method. Based on this method, fluid and structural dynamic behaviors of a highly folded disk-gap-band parachute are obtained. Numerical results such as maximum Root Mean Square (RMS) drag, general canopy shape and the smallest canopy projected areas in the terminal descent state are in accordance with the wind tunnel test results. This analysis reveals the inflation law of the disk-gap-band parachute and provides a new numerical method for supersonic parachute design.     

亚跨超CFD软件平台研制

亚跨ＣＦＤ软件平台研制的宗旨是综合国内外计算流体动力学的最新研究成果，用计算机软件技术的新发展，建立一个面向对象的计算流体动力学的研究开发环境，加快计算流体动力学最新成果向工程应用的转换速度。本文介绍初级版本的研究进展和研究成果，主要匝数值网格生成技术，流场解算器，面向对象的界面设计等三部分组成。     

可压缩流动SIMPLE算法压力修正方程边界条件处理探析

可压缩粘性流体Nivier-Stokes方程SIMPLE算法求解中压力修正方程起着关键的作用.提出从连续性方程本质出发确定压力修正方程边界条件的方法.通过对双圆弧(凸包)通道内流场的计算验证了方法的正确性.