结合前沿推进的Delaunay三角化网格生成及应用

采用一种新的混合网格生成方法,生成复杂区域的非结构化网格.结合前沿推进法和Delaunay三角化两种非结构网格生成方法的特点,在边界处采用前沿推进法进行三角形初始网格的生成,在边界区域内部采用Delaunay三角化方法自动生成内部节点.分析表明,该算法简化网格生成过程,能够快速有效地生成非结构化网格.在计算时间以及网格的均匀性方面与其他方法相比具有一定的优势.最后,用混合网格生成方法生成方柱绕流的计算域网格,并运用基于特征线方程的分离算法进行流场计算.     

基于自适应非结构化网格的VOF方法

本文提出了一种基于自适应非结构化网格的VOF算法,根据相界面网格的相函数值对相界面网格进行自适应细化与合并,通过基于非结构化网格的界面构造方法构造相界面,在适量增加网格单元数量的情况下提高了计算的精度。该方法随着时间及相界面的变化无需重新整体生成网格,算法效率较高。经典算例的验证结果表明,本文自适应网格方法计算得到的结...     

基于结构化动网格的涡旋压缩机内部流场模拟

为了增加涡旋压缩机啮合间隙处的动网格数量。本文根据涡旋压缩机的流体区域变化规律,提出一种结构化动网格的生成方法,实现了啮合间隙处的网格加密。进而得到一种基于结构化动网格的涡旋压缩机非定常流动的数值模拟方法;得到了任意时刻流体区域内的流场分布。与现有的非结构化动网格进行比较,结构化动网格解决了非结构化动网格在啮合间隙处易产生负体积的问题,能够在啮合间隙处生成足够的网格数量,提高了计算精度;且具有网格质量高、网格数量和疏密易控制等优点,更适合于涡旋压缩机工作过程的数值模拟计算。     

压气机变形叶片的网格生成及区域分解方法

针对压气机优化设计和多学科分析对网格自动可靠生成方法的需求,建立了一套集叶片几何参数化建模、多块结构化网格生成、区域自动分解、叶片自由变形于一体的多功能网格生成系统。基于椭圆型微分方程法和O4H型拓扑结构,建立了压气机单叶道流域的三维分块结构化网格生成方法,改进了适合于变形叶片及网格的自由变形方法,自主开发了相应的计算程序,并用NASA Stage35对方法及程序进行了验证。结果表明网格生成质量较好,特别是近壁网格正交性和光滑性俱佳;改进的自由变形方法能够保持网格周期性和轮毂轮盖几何形状;流场计算结果与实验数据吻合良好。研究工作对发展多级高速轴流压气机优化设计和多学科分析系统具有重要参考价值。     

一种非结构化同位网格算法

本文介绍了任意网格系统中的SIMPLE算法的推导过程，并将此算法应用于二维非结构化网格（三角形网格）之中，网格由DELAUNAY三角化方法生成．对顶盖驱动空腔流的计算表明该算法是令人满意的．     

短突扩压器流场计算

本文在三维贴体坐标系下对短环突扩扩压器进行流声计算．利用偏微分方程法和曲面论生成三维贴体网格，在按非交错网格系统求解各守恒方程时，采用抑制措施，避免压力振荡产生．考虑到扩压器是多通道区域，本文采用区域法进行数值分析，计算结果表明本计算方法合理可行．     

结构网格与非结构网格联合技术及在内爆问题中的应用

对纯拉氏方法，格式的健壮性、精度和效率强烈地依赖于初始网格的生成。结构网格与非结构网格的联合技术，不仅可适应于复杂区域的网格生成，而且可根据实际问题的特殊性，在初始自适应的生成网格。本文在研究结构网格和非结构网格生成的基础上，研究了结构网格与非结构网格联合使用的自适应网格生成技术。首先是根据实际问题的特征，把整个计算区域划分成若干子区域，其次是把每个子区域分成若干子块，最后根据每个子块的特征生成网格，再把子块拼接起来。     

离心叶轮及无叶扩器内湍流的非结构化网格数值解法

利用非结构化网格上的有限体积法对三维湍流进行了数值研究,提出了不用求解单元顶点处变量值的二阶混合差分格式,既避免了复杂耗时的单元顶点处变量值的计算,又避免了普通差分格式在非结构化网格上易于引起网格界面方向相关性问题。最后利用非结构化网格计算程序,数值求解了一台离心风机内叶轮及无叶扩压器通道内三维相对定常湍流,计算结果与实验值的比较表明速度的计算值与实验值吻合较好。     

基于非结构化网格的多相流场求解与表面张力优化

本文推导了基于非结构化网格的投影法,将其与基于非结构化网格的PLIC方法进行耦合,实现了多相流动相界面的精确构造。从CSF模型出发,对相界面表面张力在非结构化网格上进行了优化,实现了表面张力的精确计算。对静态液滴算例进行了模拟,分别验证了方法的静态和动态准确性。最后模拟和分析了方腔顶盖反向驱动剪切流场中的液滴变形和破裂过程,发现两相流体的黏度比与液滴破裂的临界毛细数呈分段线性关系;液滴半径与临界毛细数之间呈线性关系。     

一种非结构/结构多层混合网格方法及其应用

对于粘性绕流的数值模拟,在自适应直角网格基础上,结合三角形非结构网格和结构化网格,利用其各自的优势和特点,提出一种生成混合杂交网格的思路和方法.在物面附近生成适合粘性流计算的大长宽比结构化网格,在远场分布自适应直角网格,快速离散计算空间.对于复杂的多体问题,采用三角形网格来连接各体网格,并运用网格合并的方法,保证各网格之间的光滑过渡与连接,提高网格质量.针对一些二维、三维外形的绕流问题,在上述网格基础上,采用B-L代数湍流模型和中心有限体积法,完成Navier-Stokes和Euler方程数值模拟的对比计算,结果表明网格生成和流场计算是正确的.     

应用笛卡尔非结构切割网格进行外挂物投放的数值模拟

描述了一种新的网格生成技术,即笛卡尔非结构切割网格技术,采用叉树数据结构,完成了几种单段和多段翼型以及三维机翼的网格生成.应用中心有限体积法,对其绕流问题进行Euler方程数值模拟,并将计算结果与实验数据进行对比.在机翼绕流数值模拟的基础上,求解出机翼带外挂物的分离投放的流场计算问题.     

非结构化网格的扩散通量计算方法评价

本文分析对比了5种用于非结构化网格的扩散项的离散格式的精度和运算量。结果表明,除简单投影法外,其它方法的精度相差不大,但计算量差别较大。在此基础上,改进了Kobayashi的法向梯度法,克服了它对网格质量依赖较强和易于发散的缺点。总的结果表明在非结构化网格上提高扩散通量的计算精度仍然是一个值得继续研究的问题。     

含运动物体流场计算的自适应非结构二维网格生成方法

采用自动插点的Delaunay方法和局部网格重新生成方法,对含运动边界自适应非结构二维网格生成方法进行了系统研究,实现了带有非定常运动物体流场的数值计算,并进行了数值实验.     

用墨浸法自动生成三维无结构网格

提出了为用Delaunay法自动生成无结构网格而发展的"墨浸式"加点法。采用"墨浸式"加入边界点的方法,使得Delaunay方法的固有缺陷,即可能出现边界穿透的问题,得到简化。"墨浸式"加点法成功地在一些复杂边界域内生成了可用于流场计算的无结构网格。     

子网格在DSMC中的作用及三角形子网格

本文采用直接模拟蒙特卡罗方法对二维微通道内Poiseuille流动进行了数值模拟,通过对不同网格密度下计算结果与"收敛"结果的比较,得到了计算网格特别是子网格系统对计算结果的影响;在此基础上,发展了基于Delaunay三角形非结构网格的子网格系统,讨论了采用该子网格的非结构化网格对模拟结果的影响,结果表明:该方法能有效提高模拟的准确性,而常规无子网格的非结构网格将"虚假"地增加气体流动的稀薄性。     

多点择优推进阵面法生成曲面三角形网格

在现有曲面非结构网格生成法的基础上,提出了一种新的曲面网格生成法——多点择优推进阵面法.它可在曲面上直接进行三角形网格划分,克服了映射法的网格变形问题,并且可以在网格生成结束后,对曲面网格直接进行Laplace格点松弛光顺.该方法使用简单,不受曲面块类型的限制,且网格质量高,可以为三维非结构网格生成提供高质量的初始阵面,并给出了若干个算例.     

八叉树网格到非结构混合网格的转换

在数值模拟中, 非结构网格的优势是可以采用相同的数值格式统一处理任意复杂的计算区域, 但在网格生成过程中难度大, 并且不容易控制网格质量。树结构网格可以认为是介于结构网格和非结构网格之间的一种网格, 目前已经有相对成熟的方法快速在复杂区域内生成二维四叉树网格和三维八叉树网格。在实际应用中, 数值方法往往需要在连接协调的非结构网格上做离散, 树结构网格中不同尺寸的网格之间连接不是协调的, 在应用上会受到很多限制。文章实现了树结构网格到非结构混合网格的转换, 这种转换在二维情况下就是将四叉树网格转换为非结构三角形和四边形的混合网格, 三维情况下则将八叉树网格转换为非结构混合网格。这一转换过程的难点在于需要考虑数千种不同的八叉树单元, 并给出能实现连接协调的非结构混合网格划分。可以出现的网格单元包括六面体、三棱柱、金字塔和四面体这4种不同情况。通过特别的分类, 实现了程序的自动生成, 这种程序自动生成技术一方面可以避免人工编写大量程序时的失误, 另一方面也使得对数以千计的不同情况的处理成为可能。通过对几个简单网格的测试, 对网格数据转换方法做了初步的验证。      

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

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

Solution of N-S equations based on the quadtree cut cell method

With the characteristic of the quadtree data structure, a new mesh generation method, which adopts square meshes to decompose a background domain and a cut cell approach to express arbitrary boundaries, is proposed to keep the grids generated with a good orthogonality easily. The solution of N-S equations via finite volume method for this kind of unstructured meshes is derived. The mesh generator and N-S solver are implemented to study two benchmark cases, i.e. a lid driven flow within an inclined square and a natural convection heat transfer flow in a square duct with an inner hot circular face. The simulation results are in agreement with the benchmark values, verifying that the present methodology is valid and will be a strong tool for two-dimensional flow and heat transfer simulations, especially in the case of complex boundaries. Supported by Japan Society for the Promotion of Science (Grant No. C20560175) and the National Natural Science Foundation of China (Grant Nos. 10872159 and 40675011)