一种简单快速的多块网格生成方法

利用有限元中的形函数进行网格插值映射分块生成网格,构造了一种块网格拓扑结构,使得块与相邻块之间自动组合,并利用该方法对二维、三维复杂区域进行分块网格生成,说明该方法生成的网格质量较好,适合于复杂区域分块求解。     

一种部分连续复合网格生成方法

本文描述了一种构造部分连续复合网络的偏微分方程的数值生成方法，该方法以泊松方程为变换控制方程，构造特定的边界条件，对复杂的蜗尾船型的横剖面进行了变换，变换出来的复合网格与单一区域网格相比，生成的网格质量有明显提高，并计算了另外一个例子，证明了该方法的通用性。     

一种部分连续复合网格生成方法

本文描述了一种构造部分连续复合网格的偏微分方程的数值生成方法．该方法以泊松方程为变换控制方程，构造特定的边界条件，对复杂的蜗尾船型的横剖面进行了变换．变换出来的复合网格与单一区域网格相比，生成的网格质量有明显提高．并计算了另外一个例子．证明了该方法的通用性。     

CFD多块网格生成新进展

网格生成是计算流体力学的重要组成部分,多块网格在CFD实践中获得了广泛的应用．结合对网格生成技术规范和网格生成系统的讨论,综述了多块网格近年来的新进展,重点评述了网格拓扑和网格拼接技术（包括所谓的连续拼接、非结构拼接和广义拼接）,完整飞机外形多块网格生成策略,自动分块技术以及相应的块合井技术,CAD和CFD之间的数据交换技术和基于NURBS的曲面网格生成技术,网格质量分析和控制技术,若干网格生成新方法,以及多块网格在航空气动力数值模拟中的应用．      

两级分块结构化网格在SK型静态混合器中的应用

给出一种适用于SK型静态混合器流场数值模拟的两级分块结构化网格划分方案和分块圆弧区点的极坐标方程。与常用的非结构四面体网格划分对比,两种方法划分网格数量相近时,两级分块结构化网格达到的迭代精度高于非结构化网格一个等级,而非结构化网格收敛所用时间约为结构化网格收敛所用时间的1.25倍;运用控制体积容积变化值对两种方法划分网格进行质量评估,结构化网格划分方法接近99%网格质量指标在1.0~2.0范围内,而非结构化网格仅接近60%网格质量指标在1.0~2.0范围内。通过实验测量对两种方法划分网格数值模拟结果进行对比验证,结构化网格计算结果与实验值更为接近,表明其可行性和准确性均优于非结构四面体网格划分。     

无网格算法在多段翼型流动计算中的应用

研究了一种求解欧拉方程的无网格算法，发展出了一套布点及点云自动生成的方法；在点云离散的基础上，采用最小二乘法求解矛盾方程的方法来求取空间导数，进而获得数值通量；采用四步龙格-库塔方法进行时间推进，并引入当地时间步长和残值光顺等加速收敛措施。通过对NA-CA0012翼型的跨音速流动和多段翼型复杂绕流的数值模拟，验证了上述无网格算法的正确性和实用性。     

基于聚合多重网格方法的跨音速颤振的数值模拟

采用流固耦合方法对跨音速颤振进行了数值模拟。流体方面在非结构网格上用有限体积方法求解了Euler方程;结构方面则求解了后掠机翼典型剖面的结构模态方程。时间推进采用双时间步长：对每一真实时间步,都通过基于聚合多重网格方法的伪时间步推进,对流体和结构方程交替迭代．得到一个稳态的流固耦合的解。文章最后给出了NACA64A010翼型剖面的跨音速颤振边界．与相关文献的计算结果符合良好。     

一种快速稳健的并行多块结构动网格方法

为解决传统网格处理方法不能满足复杂外形在大设计空间内进行优化时对网格质量的要求的问题,提出了一种并行多块动网格方法,该方法基于初始外形的多块结构网格,根据优化过程中个体外形与初始外形拓扑结构相近的特性,利用体样条插值方法来拟合多块结构网格各块顶点的位移,得到几何外形变化后的拓扑结构,再利用无限插值方法并行地移动初始外形多块结构网格的边、面和块内的网格点,进行光顺处理后得到变形后几何外形的空间网格;该方法在保证网格质量的同时,可以极大地提高网格生成效率,本文以某翼身组合体为例结果表明,该方法在大设计空间的复杂外形设计问题中具有很强的实用性。     

一种有效的广义特征值分析方法

提出了一种适合于自适应有限元分析中求解广义特征值问题的多重网格方法．这种方法充分利用了初始网格下的结果，通过插值或最小二乘拟合技术来得到网格变化后的新的近似特征向量，然后由多重网格迭代过程实现对结构广义特征值问题的求解．在多重网格迭代的光滑步中，选择了收敛梯度法以提高其收敛率；在粗网格校正步中，则导出了一种近似求解特征向量误差的方程．这种方法将网格离散过程和数值求解过程很好地相结合，建立了一个网格细分后广义特征值问题的快速重分析方法，与传统有限元方法相比较，具有计算简便、计算量少等特点，可以作为结构动力问题自适应有限元分析的一种十分有效的工具．     

静气动弹性计算方法研究

对基于结构网格的Euler方程及N-S方程求解器和基于非结构网格的Euler方程求解器,采用结构模态分析方法和柔度矩阵方法,对无人机大展弦比机翼在Ma=0.6,α=2?, 飞行高度20km的巡航状态下的静气动弹性特性进行了数值模拟. 验证了两种求解器对静气动弹性模拟的准确性. 同时,对模态分析方法和柔度矩阵方法进行了对比研究,发现柔度矩阵方法更适用于静气动弹性数值模拟. 另外,对应用物面法向偏转方法替代网格变形技术模拟静气动弹性进行了研究,计算表明物面法向偏转方法可以大大提高静气动弹性计算效率和克服机翼结构变形过大时动网格技术无法处理的不足.      

翼吊式双发民机机体/动力装置一体化数值分析

介绍了多块网格技术与流场分区求解方法在翼吊式双发民机机体／动力装置一体化研究中的应用。数值求解Ｅｕｌｅｒ方程模拟复杂组合体绕流。采用边界层方程／Ｅｕｌｅｒ方程耦合迭代技术进行翼面粘性修正。为保持Ｅｕｌｅｒ注解中计算网格固定,粘流／无粘流耦合迭采用表面源模型。该方法对某民用飞机模型跨音速绕流流场进行了数值模拟,机翼表面计算压力分布与实验吻合良好。     

An algorithm of finite element grid generation for three-dimensional complex connected domain

An automated quasi three-dimensional finite element grid generation method is presented for particular three-dimensional complex connected domain, across which some are simply connected two-dimensional.regions and some are multiply connected two-dimensional regions. A subdivision algorithm based on the variational principle has been developed to ascertain the smoothness and orthogonality of the generated grid in any cross sections. Smooth transition between the simply and multiply connected regions is maintained. For illustration, the method is applied to generate a finite element three-dimensional grid for human above knee stump.     

An adaptive Cartesian grid generation method for ‘Dirty’ geometry

Traditional structured and unstructured grid generation methods need a ‘water‐tight’ boundary surface grid to start. Therefore, these methods are named boundary to interior (B2I) approaches. Although these methods have achieved great success in fluid flow simulations, the grid generation process can still be very time consuming if ‘non‐water‐tight’ geometries are given. Significant user time can be taken to repair or clean a ‘dirty’ geometry with cracks, overlaps or invalid manifolds before grid generation can take place. In this paper, we advocate a different approach in grid generation, namely the interior to boundary (I2B) approach. With an I2B approach, the computational grid is first generated inside the computational domain. Then this grid is intelligently ‘connected’ to the boundary, and the boundary grid is a result of this ‘connection’. A significant advantage of the I2B approach is that ‘dirty’ geometries can be handled without cleaning or repairing, dramatically reducing grid generation time. An I2B adaptive Cartesian grid generation method is developed in this paper to handle ‘dirty’ geometries without geometry repair. Comparing with a B2I approach, the grid generation time with the I2B approach for a complex automotive engine can be reduced by three orders of magnitude. Copyright © 2002 John Wiley & Sons, Ltd.     

Generation of linear and nonlinear waves in numerical wave tank using clustering technique-volume of fluid method

A two-dimensional(2D) numerical model is developed for the wave simulation and propagation in a wave flume.The fluid flow is assumed to be viscous and incompressible,and the Navier-Stokes and continuity equations are used as the governing equations.The standard k-ε model is used to model the turbulent flow.The NavierStokes equations are discretized using the staggered grid finite difference method and solved by the simplified marker and cell(SMAC) method.Waves are generated and propagated using a piston type wave maker.An open boundary condition is used at the end of the numerical flume.Some standard tests,such as the lid-driven cavity,the constant unidirectional velocity field,the shearing flow,and the dam-break on the dry bed,are performed to valid the model.To demonstrate the capability and accuracy of the present method,the results of generated waves are compared with available wave theories.Finally,the clustering technique(CT) is used for the mesh generation,and the best condition is suggested.     

直接控制网格节点分布的曲线网格生成技术研究

指出了Thompson与Thomas曲线网格生成方法中控制网格分布的调节函数的问题所在,克服了Thomas曲线网格生成法中边界处局部线性化近似假定的缺陷,经过严格推导得出一组新的调节函数P、Q的表达式,并给出了曲线网格生成实例.实例检验表明,该调节函数能够对复杂边界的单连通域或多连通域生成理想的曲线网格,即边界处网格正交,内部网格分布能够适应物理量场变化的情形.在实际水利工程流场数值模拟中,该方法能够准确地使用边界条件,提高求解的精确度.     

A large-scale parallel hybrid grid generation technique for realistic complex geometry

High-Performance Computing (HPC) systems and Computational Fluid Dynamics (CFD) have made significant progress in recent years; however, as the basis of the large-scale parallel computing, the massive grid generation of billions of cells has become a bottleneck problem. In this study, a parallel grid generation technique is proposed to generate large-scale mixed grids with arbitrary cell types and scales. The basic idea of our method is analogous to the global mesh refinement technique. An initial coarse grid with arbitrary cell types is regarded as a background mesh which is partitioned into subzones, and subzones are assigned onto different CPU cores. After the cells and faces in each subzone are split, the inserted new points of the solid wall are projected onto the original CAD entities to preserve the geometry accurately. Finally, the tangled cells caused by the projection in the boundary layer are untangled by a local Radial Basis Function mesh deformation technique. Furthermore, a parallel partition approach and an efficient wall distance computing technique for massive grids are developed also to shorten the preprocessing time. The tests show that the preprocessing efficiency has been increased by two or three orders compared with traditional methods. Billions of grids are generated for the AIAA JSM high-lift model and the Chinese CHN-T1 transport model to test the ability of the parallel grid generation technique. The maximum scale up to 19 billion mixed elements is generated using 16 384 CPU cores in parallel, and the mesh quality is acceptable for CFD simulations.     

A NON–LINEAR ADAPTED TRI–TREE MULTIGRID SOLVER FOR FINITE ELEMENT FORMULATION OF THE NAVIER–STOKES EQUATIONS

An iterative adaptive equation multigrid solver for solving the implicit Navier–Stokes equations simultaneously with tri-tree grid generation is developed. The tri-tree grid generator builds a hierarchical grid structur e which is mapped to a finite element grid at each hierarchical level. For each hierarchical finite element multigrid the Navier–Stokes equations are solved approximately. The solution at each level is projected onto the next finer grid and used as a start vector for the iterative equation solver at the finer level. When the finest grid is reached, the equation solver is iterated until a tolerated solution is reached. The iterative multigrid equation solver is preconditioned by incomplete LU factorization with coupled node fill-in. The non-linear Navier–Stokes equations are linearized by both the Newton method and grid adaption. The efficiency and behaviour of the present adaptive method are compared with those of the previously developed iterative equation solver which is preconditioned by incomplete LU factorization with coupled node fill-in.     

HEAT FLOW CONTROL BY TRANSFORMATION METHOD WITH GRID GENERATION METHOD

In the context of the transformation method, we propose a general approach to construct numerically the mapping generated by imposing specific boundary conditions with a targeted function, and the necessary material and heat source spatial distributions are then derived with the help of transformation method. The construction of mapping by grid generation method through solving partial differential equations circumvents the limitation of device geometry, which paves the way for designing more complex heat flow control devices. Two numerical examples are also given to show how to design material properties and heat source in order to control temperature patterns.     

AN ADAPTIVE UNSTRUCTURED TRI-TREE ITERATIVE SOLVER FOR MIXED FINITE ELEMENT FORMULATION OF THE STOKES EQUATIONS

An iterative adaptive equation solver for solving the implicit Stokes equations simultaneously with tri-tree grid generation is developed. The tri-tree grid generator builds a hierarchical grid structure which is mapped to a finite element grid at each hierarchical level. For each hierarchical finite element grid the Stokes equations are solved. The approximate solution at each level is projected onto the next finer grid and used as a start vector for the iterative equation solver at the finer level. When the finest grid is reached, the equation solver is iterated until a tolerated solution is reached. In order to reduce the overall work, the element matrices are integrated analytically beforehand. The efficiency and behaviour of the present adaptive method are compared with those of the previously developed iterative equation solver which is preconditioned by incomplete LU factorization with coupled node fill-in. The efficiency of the incomplete coupled node fill-in preconditioner is shown to be largely dependent on the global node numbering. The preconditioner is therefore tested for the natural node ordering of the tri-tree grid generator and for different ways of sorting the nodes.     

采用自适应直角网格计算三维增升装置绕流

针对三维增升装置绕流，对存在剪刀叉的不连续外形，基于自适应直角网格，提出并介绍了分区和面搭接技术，采用变长宽比网格，进行了直角网格生成和流场Euler方程数值计算. 根据几何外形的特点，在直角网格生成过程中，以外形不连续面作为分区边界，对初始``根'网格实施分区处理，降低了整个网格的生成难度. 通过基于外形的自适应网格加密，详细描述了剪刀叉外形和缝道，提高了网格质量. 在分区边界面上，基于面搭接技术，构造重叠面积切割算法，实现边界两侧网格间的流场信息传递，保证流场计算中的通量守恒. 采用中心有限体积方法，结合双时间推进算法，完成了两段机翼、带增升襟翼翼身组合体绕流流场的Euler方程数值模拟，对计算结果与实验数据进行了对比，验证了所提方法、算法的合理性和实用性.