非结构动网格生成的弹簧-插值联合方法

 提出了基于弹簧法的两点改进,以解决边界发生大变形时的非结构网格变形问题。为了使边界运动引起的网格变形能更好地由边界传递到内部网格中,提出了一种新的基于Delaunay网格插值的弹簧倔强系数逐层改进方法。为了提高弹簧法的计算效率,引入背景网格和直接插值方法,提出了弹簧-插值法。弹簧-插值法首先生成计算域的背景网格(粗网格),然后由弹簧法求解边界运动引起的背景网格变形,最后利用变形后的背景网格直接插值得到计算网格的变形。算例结果表明：改进后的方法一方面有效地提高了动网格的变形能力和变形后的网格质量;另一方面通过降低弹簧法的求解规模,显著地提高了动网格的变形效率。     

基于非结构背景网格的动网格方法研究

借鉴弹簧动网格方法和Delaunay背景图方法,提出一种基于非结构背景网格的动网格生成方法。先用少量的物面控制点生成相对较粗的非结构背景网格,再用弹簧法将动边界的变化量作用到该背景网格,然后利用计算网格与背景网格的映射关系插值生成对应的动网格,最后对物面网格点及其附近受影响的密网格点进行局部网格修正,获得最终计算所需的动网格。所提方法与直接弹簧法比较,由于弹簧法作用的背景网格远粗于实际的计算网格,提高了计算效率,此外,由于网格的运动受控于受弹簧法作用的背景网格,边界变形或运动能较为均匀地传递到内部网格上,使网格单元过渡光滑,提高了动网格的质量。给出若干算例,展示了方法的效率和可行性。     

非结构网格生成技术方法研究

利用点源或线源所影响的势函数分布规律,通过建立势函数与网格背景信息之间的关系,提供非结构网格生成过程中所需的网格背景参数;在物体边界节点已知的情况下,适当布置点源和线源个数,用最小二乘法确定各点源和线源的强度;通过对部分点源或线源的影响势函数施加方向控制开关,使所生成的网格在期望的区域得到加密;算例表明本文方法是一种有效、快速的处理方法,与传统用背景网格提供背景信息的方法相比,在计算量、程序量、光滑性和灵活性方面都能表现出良好的性质。     

两种典型非结构网格变形方法特性对比研究

非定常流动仿真中常常遇到运动边界的情形,基于物理模型的线性弹簧法和基于数学方法的径向基函数(RBFs)插值法是实现计算网格随运动边界变化的两种主要的非结构网格变形方法。以NACA0012翼型的俯仰运动为例,对上述两种方法的变形特性进行较详细地定量分析研究,比较线性弹簧法和RBFs方法的CPU计算时间和最大变形能力,分析两种方法的网格质量随俯仰角度和运动步数的关系,并给出典型状态下网格单元质量云图和流场计算结果。结果表明:基于RBFs插值的网格变形方法计算效率较高,变形能力较强且能保持较高的网格质量,是一种高效的网格变形方法。     

非结构化网格生成中背景网格对源的响应

背景网格对源的响应方法对非结构网格的生成具有重要意义,如果背景网格参数能恰当地反映出源的变化,则所生成的非结构网格就能满足网格设计的要求。将背景网格的三角形单元分为7个子域,根据源在单元中的位置对网格进行局部加密或变形,使每一个源点都与一个背景网格结点相对应,再通过求解热传导方程,使背景网格参数分布光滑,基本反映网格设计的要求,从而提高了非结构化网格的生成质量。     

结冰计算中翼面网格重构研究

对结冰数值计算中的网格生成进行了深入研究,针对结冰后几何边界的复杂性及边界依赖于时间步长变化对网格生成的影响展开了分析.网格生成是结冰计算的基础,文中基于求解椭圆形偏微分方程(PDEs)生成了未结冰时翼面周围的背景网格,采用Hilgenstock源项构造方法,合理控制了结冰边界周围的网格间距和正交性.模型表面结冰后,为...     

三维非结构自适应多重网格技术

发展了一种基于几何外形的非结构自适应多重网格技术。根据流场参数的变化梯度确定加密边,由加密准则对粗网格进行自适应剖分得到分布合理的较密网格。通过预先生成的初始极密表面网格来将边界网格的加密点投影到边界上,使得边界保持初始外形。为了提高流场求解速度,将自适应加密后的一系列网格作为多重网格,采用了多重网格算法进行计算。M6机翼的流场计算表明该方法可以准确捕捉激波位置,并且表面的压力系数分布与实验值相吻合。     

基于径向基函数与混合背景网格的动态网格变形方法

发展了一种基于径向基函数(RBF)插值和混合背景网格映射的动态网格变形方法。混合背景网格由运动物面附近的各向异性网格和偏离运动物面的各向同性网格组成。物面变形或运动后,包含在各向异性背景网格中的所有物面网格点的位移由指定的物面运动规律得到,而各向同性背景网格点的位移通过径向基函数插值获得。然后变形后的计算网格利用计算网格与混合背景网格之间固定的映射关系代数插值得到。各向同性网格距离物面一定的距离,消除了物面保形对径向基函数插值的要求,从而可以减少径向基函数基点的数量,以提高插值的效率和可靠性。最后,通过二维和三维算例对该方法的变形网格质量、参数影响、网格拓扑依赖性、复杂外形变形能力及网格变形效率进行了测试。测试结果表明,该方法能够高效生成大变形下高质量的变形网格。     

复杂外形的动态混合网格生成方法

本文将层推进方法、四分树方法和阵面推进法相结合,建立了二维复杂外形的动态混合网格生成方法.首先用层推进方法在运动物体附近生成贴体的四边形网格,然后用四分树方法在外场生成自适应的矩形网格,而中间由三角形网格过渡.三角形网格由阵面推进法生成,并由节点松弛和"对角线交换"方法进行优化.当物体运动时,贴体的四边形网格随物体运动而运动,而外场的矩形网格保持静止.中间过渡层的三角形网格随物体运动而变形.当物体运动位移较大导致出现质量较差的单元甚至导致网格相交时,则在局部重新生成网格.利用上述动态混合网格生成方法,生成了多种含相对运动的复杂外形的动态混合网格,网格质量令人满意.     

任意平面域的三角形网格和混合网格生成

本文以阵面推进法为基础,给出了一种可对任意平面域进行三角形网格和混合网格自动生成的方法。该方法包含了易于生成与调控的矩形背景网格技术,简便的外形输入和初始阵面划分技术,加速查寻的数据结构,以及节点松弛和基于Ｄｅｌａｕｎａｙ准则的对角线交换技术。网格生成的实例表明,方法能快速可靠地生成任意平面域的高质量三角形网格和适于粘性计算的三角形／四边形混合网格,具有较好的通用性与实用性。     

空间曲面的非结构网格生成

介绍一种空间曲面的 Delaunay三角形网格生成技术。采用参数表示法将空间曲面变换为平面参数域,引入相关的度量矩阵后将平面参数域进行三角剖分,最后得到空间曲面的三角剖分。空间曲面的边界点作为已知点被给定,内点采用一类自动生成的方法计算,例如按边插点法,Voronoi边插点法等。该方法可以很好地控制曲面网格的质量,使用方便,可以直接用于物理问题数值模拟中的曲面网格生成,也可以用于三维非结构网格生成中提供边界面三角剖分。文中给出的算例说明了本方法的应用情况。     

三维Euler方程的自适应非结构网络计算

介绍一种三维自适应非结构网格的Euler方程求解技术。采用有限体积中心差分的格点格式对Euler方程进行空间离散,多步龙格-库塔时间推进,并采用当地时间步长、焓阻尼修正等加速收敛。采用Delaunay三角剖分技术生成三维流场的四面体网格。网格自适应技术中的网格局部加密也采用Delaunay三角剖分技术。最后用ONERAM6机翼的跨音速绕流计算说明本文方法的有效性。     

非结构多重网格流场数值模拟研究

本文在非结构网格上应用多重网格技术来加速Euler方程的收敛过程。此方法在一系列粗网格和细网格之间进行值的传递，并且粗网格和细网格之间结构上无关。在每一重网格上用有限体积法进行计算。采用本文的方法数值求解了NACA0012翼型的绕流。     

三维复杂域的Delaunay三角剖分

介绍一类自动生成内点的Delaunay三角剖分算法 ,其中包括重心插点法 ,外心插点法和Voronoi边插点法。该类方法适用于二维和三维复杂计算域的非结构网格生成 ,所需人工干预少 ,使用方便 ,计算效率高 ,生成的网格质量也很好。文中并给出几个算例说明本方法的应用。     

三维非结构网格的生成及优化

叙述了一种新的生成三维非结构网格的Delaunay方法,该方法能够高效地生成优质四面体网格。由于采用了局部剖分的方法来保证边界的完整性,因而不需要添加新的节点到原始表面三角形中。该方法在四面体外接球的球心处插入新的节点,在插入全部内点后,采用删除边、边/面交换、网格光顺等方法对生成的网格进行了优化,基本能删除质量差的四面体。     

非结构网格多区域动态加密及分布式并行计算

针对分布式并行系统,提出一种二维非结构网格多区域动态并行策略,可在计算过程中对各计算区域的网格进行局部自适应网格加密,可在不过多消耗计算机资源的情况下,快速,有效地实现多台同构计算机上的计算网格的动态负载平衡.并针对二维Euler方程,利用已有的高精度MmB有限体积格式,应用提出的并行计算策略,进行了多个数值试验,得到令人满意的效果.     

非结构网格上欧拉方程的自适应逆风隐式有限元解法

本文发展了一种解二维欧拉方程的隐式逆风有限元格式。空间近似应用的是二阶精度逆风方法；时间近似上根据数值通量函数的线性化处理导出了有限元的隐格式。为减少数值振荡，引入了基于特征变量的限制因子。在经自适应处理的非结构网络上，本文对跨、超、高超声速绕流进行了数值模拟。计算实践表明发展的隐式有限元格式与显格式相比具有较高的计算效率和较强的适应性。     

各向异性网格自适应计算在超燃模拟中的应用

采用各向异性网格自适应求解技术,将其应用于DLR超燃冲压发动机燃烧室中的超燃模拟.开展了3个算例,包括采用滑移壁面条件的燃烧室冷流场模拟、采用无滑移壁面条件的冷流场模拟及采用无滑移壁面条件的反应流模拟.模拟中,各向异性网格自适应计算捕捉到了如激波、射流、边界层、火焰面等具有各向异性特征的大梯度区域,并利用各向异性网格进行了很好的加密.对比利用各向同性网格的初始流场计算,各向异性网格自适应计算使基于滑移条件、无滑移条件的冷流计算及反应流计算的网格单元数量分别下降了36.2%,36.4%和36.8%,有效降低了计算规模,而且流场大梯度区域的计算结果更准确,辨析度更好.结果表明:对于像超燃这类具有各向异性特征的问题,各向异性网格系统比各向同性网格系统有更好的计算效率及准确性,同时也表明基于Mach数场构造的各向异性网格系统可以有效应用于超燃计算.      

Delaunay graph-based moving mesh method with damping functions

The fluid–structure interaction and aerodynamic shape optimization usually involve the moving or deforming boundaries, thus the dynamic mesh techniques are the key techniques to cope with such deformation. A novel dynamic mesh method was developed based on the Delaunay graph in this paper. According to the Delaunay graph, the mesh points were divided into groups. In each group, a factor ranging from 0 to 1 was calculated based on the area/volume ratio. By introducing a proper function for this factor, this method can control the mesh quality with high efficiency. Several test cases were compared with other dynamic mesh methods regarding mesh quality and CPU time, such as radial basis function method and Delaunay graph mapping method.     

叶片参数化有限元建模与罩量优化设计

基于网格变形技术提出了一种叶片参数化有限元建模方法,该方法具有快速实现叶片参数化有限元建模的优点,同时能够在保证叶片有限元网格拓扑形式不变的前提下,有效解决叶片优化设计过程中的网格自动更新问题.在此基础上,结合叶片罩量优化设计的灵敏度分析技术,建立了叶片罩量自动优化一般设计流程,开展了叶片实例优化设计应用研究.优化结果表明,相对于初始设计,罩量优化设计可以有效地改善应力分布,降低叶片应力水平.