三角B-B曲面在非结构网格生成中的应用

基于生成非结构化网格的Delaunay三角形化方法,应用三次三角Bernstein-Bézier曲面来控制网格点的分布情况。对于给定边界,利用三次非均匀B样条进行边界的拟合及重新离散。对初始化形成的Delaunay三角形,应用三次三角Bernstein-Bézier曲面来计算位于其中的点的长度标尺,通过对三角Bernstein-Bézier曲面边界点的定义来控制三角形内长度标尺的分布。对复杂通道网格剖分的实例表明,此方法可以很好地控制内部网格点的分布情况。     

用非结构网格与欧拉方程计算复杂区域的二维流动

提出用Delaunay三角化方法生成非结构网格的一种过程。所生成的网格可用于复杂多连通域内的可压流计算。采用Euler方程和格心有限体积法,研制出程序,给出了算例。     

三维非结构网格自动生成

采用各向异性平面非结构网格生成技术对曲面的参数平面进行三角剖分,从而得到曲面的非结构网格,作为三维非结构网格的边界网格.应用推进面法生成网格内点,增量法将生成的内点逐点插入现有网格进行网格细化,得到三维计算域的Delaunay非结构网格.讨论了非结构网格质量优化方法.给出几个算例说明方法的应用.     

基于Steiner点扰动和矢量边界推进技术的三维约束非结构网格生成方法

提出一种结合Steiner点扰动和矢量边界推进三角化技术的三维约束非结构四面体网格生成方法.在由物面节点生成的Delaunay网格基础上,利用Conforming方法恢复计算域形状,然后利用点扰动技术迫使所有Steiner点全部从约束边、面上转移,并利用矢量边界推进三角化方法重构约束面三角剖分,进而得到完全恢复所有约束条件的Constrained网格,并证明方法的收敛性和稳定性.     

一种新的非结构四边形网格生成方法及软件

提出一种间接生成非结构四边形网格的方法.采用Delaunay三角形网格生成的算法生成三角形单元,在三角形单元中定义生成四边形网格的前沿,根据三角形的特点(边、角大小)以及生成四边形前沿的方向,将两个三角形合并为一个四边形单元.此方法的特点是不增加节点,仅改变点与点之间的拓扑连接,使单元基本变为四边形占优的单元.然后研究四边形占优的网格转换为100%的四边形网格的转换模版.并应用于实际物理问题.数值实验显示了该方法生成非结构四边形网格的特点.     

非结构网格的生成及新型数据类型的应用

本文基于 Ｄｅｌａｕｎａｙ三角化方法及局部重组的阵面推进法，给出了任意形状平面区域的非结构化网格剖分方法。文中对加点方法及数据结构和存取模式进行了改进，同时应用链表结构、动态数组和派生类型的数据变量，对算法进行优化，节省了计算内存，提高了网格的生成效率。本文研究为特殊计算区域的三维非结构化网格的生成提供了基本方法。     

基于约束Delaunay三角化的二维非结构网格生成方法

给出基于局部重构和边交换技术的两种约束Delaunay三角剖分方法并证明其收敛性.采用边界指示法恢复流场形状;在预设尺度的指导下融合流场边界曲率、中轴线、梯度限制等信息修正流场尺度;运用Spring方法布置边界点,通过符号面积函数和概率筛选法布置计算区域节点;运用Spring-Laplace方法优化节点位置,伴同边交换和边吞噬技术优化网格结构.该方法可自由进行局部自适应加密或稀疏,并应用于映射曲面网格生成和移动网格技术.     

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

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

振动桥梁CFD数值模拟的网格运动算法

针对非流线型桥梁断面模型,采用Delaunay背景网格方法、棱边弹簧及线性弹簧网格方法,研究非结构网格与结构网格条件下,计算域网格运动的特点及算法优劣.研究表明,在正常的桥梁断面振动振幅内,三种网格运动方法均能保持网格变形时不会失效,其中Delaunay背景网格方法耗时最少,线性弹簧网格方法获得网格质量最高;在大振幅状态下,结构化网格容易失效,线性弹簧网格方法可最大限度保证网格有效.     

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

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

Spectral difference method for compressible flow on unstructured grids with mixed elements

This paper presents the development of a 2D solver for inviscid and viscous compressible flows using the spectral difference (SD) method for unstructured grids with mixed elements. A mixed quadrilateral and triangular grid is first refined using one-level h-refinement to generate a quadrilateral grid while keeping the curvature of boundary edges. The SD method designed for quadrilateral meshes can subsequently be applied for the refined unstructured grid. Results obtained with the SD method for both inviscid and viscous compressible flows compare well with analytical solutions and other published results.     

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

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

任意网格重映的样条逼近算法

在大变形流体力学问题的数值模拟中,任何方法都必须考虑网格重分或网格自适应,只要改动网格就涉及重分,或自适应后从旧的、扭曲的网格到新网格的守恒量重映,包括质量、动量和能量.在研究样条函数逼近的基础上,给出一种物理量重映的对结构网格和非结构网格均适应的算法,并给出了数值结果.     

Construction of Commutative Filters for LES on Unstructured Meshes

A method of constructing discrete filters for large eddy simulation of turbulent flows on unstructured meshes is presented. The commutation error between differentiation and filtering can be made arbitrarily small with these filters. The filtering method is applied to various test cases to demonstrate commutation. An extension to three dimensions and implementation into an unstructured solver for LES are discussed.     

Modified finite-volume method based on a cell vertex scheme for the solution of radiative transfer problems in complex 3D geometries

A modified finite-volume method based on a cell vertex scheme was applied to solve radiative transfer problems within a participating medium of complex three-dimensional shaped domain. The computational spatial domain of interest was divided into four-node tetrahedron elements with unstructured meshes while the adopted formulation was combined with a closure relation based on an exponential scheme. The studied medium was assumed to be grey, non-scattering and was bounded by black surfaces. Our results were then compared with those found in other articles on the subject. The approach shows a very good level of performance for wall heat transfer evaluation. Accurate results were obtained on coarse computational meshes and solution errors were found to decrease with grid refinement.     

二维非结构网格有限体积格式的壁热现象和人工热流方法

基于Noh人工热流方法的思想,推导出任意多边形网格人工热流的格式,将Noh的结构网格人工热流方法推广到任意多边形非结构网格的拉氏有限体积方法中,抑制或消除了壁热现象.几个数值实验,包括Noh激波反射问题和两块钢金属碰撞问题,验证了方法的可行性.     

多尺度模拟中网格守恒重映算法

针对耦合微观分子动力学(MD)和宏观有限元方法(FE)的多尺度模拟,提出一类新的基于贡献单元法的网格守恒重映算法.由于物理量是由有限元节点以及相应区域的原子信息通过积分重构得到的,对结构和非结构网格都能适用.对于未知量定义在顶点的情形,引入辅助网格.数值例子验证了算法的准确性和有效性.     

热传导方程的一类无网格方法

构造求解热传导方程的一类无网格方法,只要选择好每个节点的适当的邻点集合,便可利用节点信息顺利进行计算.作为特殊情形,也可在各种结构或非结构网格上进行计算.在矩形或均匀平行四边形网格上进行计算时具有二阶精度,当在任意的不规则四边形或三角形网格上计算时仍然是守恒的和相容的,且至少具有一阶精度.作为数值试验,将该方法用于在不规则四边形网格上及四边形与三角形混合网格上求解二维非线性抛物型方程,并在不规则四边形网格上求解二维三温辐射热传导方程,均获得了较为理想的数值结果.     

A Pseudo-Stokes Mesh Motion Algorithm

This work presents a moving mesh methodology based on the solution of a pseudo flow problem. The mesh motion is modeled as a pseudo Stokes problem solved by an explicit finite element projection method. The mesh quality requirements are satisfied by employing a null divergent velocity condition. This methodology is applied to triangular unstructured meshes and compared to well known approaches such as the ones based on diffusion and pseudo structural problems. One of the test cases is an airfoil with a fully meshed domain. A specific rotation velocity is imposed as the airfoil boundary condition. The other test is a set of two cylinders that move toward each other. A mesh quality criterion is employed to identify critically distorted elements and to evaluate the performance of each mesh motion approach. The results obtained for each test case show that the pseudo-flow methodology produces satisfactory meshes during the moving process.