非结构动网格在可动边界问题中的应用研究

本文对用于非结构动网格生成的弹簧近似方法进行了研究。通过采用顶点弹簧方法,分析研究了弹簧倔强系数的取值,同时通过引入挤压倔强系数和边界修正,对标准弹簧近似方法进行了改进。改进后的方法可以大大提高网格变形能力和网格质量。应用本文发展的非结构动网格生成方法并通过耦合求解基于(ArbitraryLagrangian-Eulerian ALE)描述的Euler方程,模拟了谐和振动NACA0012翼型及M6机翼的跨音速绕流,计算结果与参考文献提供的结果及实验结果吻合良好。     

非结构动网格在多介质流体数值模拟中的应用

采用非结构动网格方法对含多介质的流场进行数值模拟.采用改进的弹簧方法来处理由于边界运动而产生的网格变形.采用基于格心的有限体积方法求解守恒型的ALE(Arbitrary Lagrangiall-Eulerian)方程,控制面通量的计算采用HLLC(Hartem,Lax,van Leer,Contact)方法,采用几何构造的方法使空间达到二阶精度,时间离散采用四阶Runge-Kutta方法.物质界面的处理采用虚拟流体方法.本文对含动边界的激波管、水下爆炸等流场进行数值模拟,取得较好的结果,不同时刻界面的位置和整个扩张过程被准确模拟.     

当地DFD方法在扑翼流动数值模拟中的应用

应用当地DFD（Domain-Free Discretization）方法对包含复杂运动边界的扑翼流场进行了数值模拟。该方法通过壁面法线方向的外插确定外部相关点上的流动变量值,同时强加相应的边界条件。应用这种方法,动边界流动的模拟可以在固定网格上实现,无需为了跟随物体的运动而在每个时间步上对网格进行实时更新。对三种拍动模式的双翅流场进行了数值模拟,升、阻力系数时间历程的计算结果与参考文献的实验数据吻合很好,验证了当地DFD方法处理复杂动边界问题的可靠性。最后,数值模拟了完整昆虫模型的扑翼流场,并详细分析了涡系结构和飞行机理。     

NND格式在非结构网格中的推广

在张涵信提出的无波动、无自由参数的差分格式（ＮＮＤ格式）的基础上，构造了适用于非结构网格的二阶精度ＮＮＤ有限体积格式，解决了现有非结构网格方法中为抑制激波附近的波动而必须引入含自由参数的人工粘性项的困难，并采用网格自适应技术以提高效率．通过对二维平板激波反射和前台阶在管道内的流动问题的计算，表明本方法可有效地用于Ｅｕｌｅｒ方程的求解．     

非结构网格热化学非平衡辐射流场数值计算

针对三维热化学非平衡辐射流场设计了基于非结构网格的数值计算方法. 根据原子分子光谱理论逐条计算了100$\sim$1\,500\,nm间N, O, N$^{ + }$, O$^{ + }$的谱线以及N计算流体力学; 辐射; 热化学非平衡; 非结构网格; 有限体积法针对三维热化学非平衡辐射流场设计了基于非结构网格的数值计算方法.根据原子分子光谱理论逐条计算了100～1 500nm间N,O,N+,O+的谱线以及N2,O2,NO,N+2等分子的10个谱带,特别分析了NO的β'带,γ'带,δ带和ε带的辐射特性.采用耦合辐射的双温模型计算热化学非平衡流场,辐射源项通过直接求解辐射输运方程RTE(radiative transport equation)获得.在空间和方向上分别离散后,利用有限体积法求解不同方向上的辐射输运方程.计算得出了再入飞行器前驻点的辐射强度分布.采用该数值方法计算了MUSES-C模型在速度为11.6km/s时的绕流流场及前驻点处的辐射热流密度.并通过对比分析了热辐射对流场的影响.     

非结构重叠网格方法研究及应用

航空、航天和兵器技术等领域的研究中存在大量包含运动边界的流场。非结构重叠网格方法是一种高效的处理动边界问题的新方法。围绕相对运动的每个物体单独生成非结构网格,在网格重叠区域通过搜索和插值完成网格系之间的信息传递,提出了动态八叉树搜索算法,发展了绝对坐标系和相对坐标系相结合的流场求解方式,采用二阶精度Van Leer/Hanel格式和四阶Runge-Kutta法分别进行空间和时间离散,形成了一种新的非结构重叠网格算法。对三维Riemann问题的求解结果与精确解能很好吻合,证明了本文的重叠网格算法具有较好的时空离散精度和插值精度。对7.62mm步枪射击过程进行了数值模拟,描述了弹丸离开膛口后膛口流场的发展过程,与实验结果体现的发展过程较为吻合,验证了本文提出的非结构网格算法体系具有较好的计算性能,是研究含动边界复杂流场的一种有效手段。     

N-S方程在非结构网格下的求解

在Ｒｏｅ的矢通量差分分裂的基础上，吸收了ＮＮＤ格式的优点，提出了一种非结构网格下求解Ｅｕｌｅｒ方程和Ｎ－Ｓ方程的高分辨率高精度迎风格式．这种格式具有捕捉强激波和滑移线的良好性能．在时间方向上采用了显式和隐式两种解法．文中还给出了自适应技术．最后，成功地完成了ＧＡＭＭ超音速前台阶绕流、二维平板无粘激波反射、三维Ｈｏｂｓｏｎ叶栅流动、ＶＫＩ叶栅流动、Ｃ３Ｘ叶栅流动的数值模拟，得到了满意的结果     

基于非结构网格的TTGC有限元格式的实现及在超声速流动中的应用

基于非结构网格系统,实现了时空三阶精度的TTGC有限元格式,并在三阶TTGC格式上发展了基于人工粘性的激波捕捉技术。在非结构网格下,采用这种方法对若干典型的超声速流动问题(SOD激波管、马赫数为3的前台阶流动以及马赫数为8的高超声速圆柱流动)进行了验证计算。结果表明,TTGC格式分辨率高,在粗糙网格下能够准确的模拟超声速流场中的激波、接触间断等复杂流动现象,并且能有效的控制间断附近的数值色散现象。与传统的有限体积方法相比,本文实现的TTGC有限元格式在模拟超声速流动问题方面具有格式精度高、数值耗散小等优点。     

大展弦比无人机的静气弹问题计算及分析

普遍采用大展弦比机翼的无人机,其气动弹件的问题显得尤为突出.本文主要研究基于非结构弹性网格体系的欧拉方程CFD计算及其在大展弦比无人机静气动弹性问题中的应用.针对三维非结构运动网格技术进行了研究和开发,在此基础上,利用计算流体力学的方法,发展了一套具有一定通用性的、适用于非结构网格的Euler方程流场求解器,并综合上面的技术,进一步通过耦合结构力学方程,对大展弦比无人机的静气动弹性问题进行了,计算和分析.     

GMRES方法在含动边界流场中的应用

以基于格心的有限体积法为基础，空间二阶精度，采用4阶Runge—Kutta，GMRES隐式方法求解基于ALE形式的Euler方程，网格单元边界处守恒量通量的计算采用了Hanel方法，对NACA0012翼型绕流及运动圆球绕流等问题进行数值模拟，取得了较好的结果．GMRES方法克服了以往隐式方法大量耗费内存的弱点，达到了计算耗时短和占用内存少的统一．     

ADAPTIVE SOLUTIONS FOR UNSTEADY LAMINAR FLOWS ON UNSTRUCTURED GRIDS

An adaptive finite volume method for the simulation of time-dependent, viscous flow is presented. The Navier–Stokes equations are discretized by central schemes on unstructured grids and solved by an explicit Runge–Kutta method. The essential topics of the present study are a new concept for a local Runge–Kutta time-stepping scheme, called multisequence Runge–Kutta, which reduces the severe stability restriction in unsteady problems, a common grid generation and adaptation procedure and the application of dynamic grids for capturing moving flow structures. Results are presented for laminar, separated flow around an aerofoil with a flap.     

A PRESSURE-CORRECTION METHOD FOR THE SOLUTION OF INCOMPRESSIBLE VISCOUS FLOWS ON UNSTRUCTURED GRIDS

An unstructured grid, finite volume method is presented for the solution of two-dimensional viscous, incompressible flow. The method is based on the pressure-correction concept implemented on a semi-staggered grid. The computational procedure can handle cells of arbitrary shape, although solutions presented herein have been obtained only with meshes of triangular and quadrilateral cells. The discretization of the momentum equations is effected on dual cells surrounding the vertices of primary cells, while the pressure-correction equation applies to the primary-cell centroids and represents the conservation of mass across the primary cells. A special interpolation scheme s used to suppress pressure and velocity oscillations in cases where the semi-staggered arrangement does not ensure a sufficiently strong coupling between pressure and velocity to avoid such oscillations. Computational results presented for several viscous flows are shown to be in good agreement with analytical and experimental data reported in the open literature.     

计算含动边界非定常流动的无网格算法

在无网格算法中考虑了含动边界的流动问题,研究了可以计算处理包含一定位移及扭转动边界非定常流动的算法. 创建了无网格算法的动点法则,并引入抗扭方法对弹簧方法进行改进来处理离散点运动,提高了方法的可用度及精度. 发展了求解基于无网格的ALE方程组的算法,在点云离散的基础上采用曲面逼近计算空间导数及HLLC格式计算数值通量,运用四步龙格-库塔法进行时间推进. 在跨、超音速条件下,计算模拟了典型翼型简谐振动流场,计算结果与实验结果及文献对比吻合,验证了该算法的正确性.     

Analytic study of 2D and 3D grid motion using modified Laplacian

The modified Laplacian has been used to move unstructured grids in response to changes in the surface grid for a variety of grid movement applications including store separation, aero‐elastic wing deformation and free surface flow simulations. However, the use of the modified Laplacian can result in elements with negative areas/volumes, because it has no inherent mechanism to prevent inversion of elements. In this paper, the use of a modified Laplacian is studied analytically for a two‐dimensional problem of deforming the inner circle of two concentric circles and for a three‐dimensional problem of deforming the inner sphere of two concentric spheres. By analysing the exact solution for this problem, the amount of translation and deformation of the inner circle that maintains a valid mesh is determined. A general grid movement theorem is presented which determines analytically the maximum allowable deformation before an invalid mesh results. Under certain circumstances, the inner circle and sphere can be expanded until it reaches the outer circle or sphere, while remaining a valid grid, and the inner circle and sphere can be rotated by an extreme amount before failure of the mesh occurs. By choosing the exponent to the modified Laplacian appropriately, extreme deformations for single frequency deformations is possible, although for practical applications where the grid movement has multiple frequencies, choosing the optimal exponent for the modified Laplacian may not be practical or provide much improvement. For grid movement simulations involving rigid body translation and rotation or uniform expansion, the modified Laplacian can yield excellent results, and the optimum choice of the modified Laplacian can be analytically determined for these types of motions, but when there are multiple frequencies in the deformation, the modified Laplacian does not allow much deformation before an invalid grid results. Copyright © 2006 John Wiley & Sons, Ltd.     

DYNAMIC RESPONSE OF A DOUBLE-LAYERED SUBGRADE WITH ROCK SUBSTRATUM TO A MOVING POINT LOAD

In this paper, an analytical solution for the dynamic response of a double-layered subgrade with rock substratum to a moving point load is derived. The subgrade profile is divided into two layers. The upper layer is modeled by an elastic medium and the lower layer by a fully saturated poroelastic medium governed by Biot’s theory. In the meanwhile, the subgrade is resting on the rock substratum. The analytical solutions for stress, displacement and pore pressure are derived by using the Fourier transform. Numerical results obtained by using the inverse fast Fourier transform (IFFT) are used to analyze the influence of the moving load velocity, the thickness of an elastic medium layer and a fully saturated poroelastic medium layer on the dynamic response.     

移动简谐力作用下三维多孔饱和半空间的动力问题

研究了移动荷载作用下多孔饱和地基的动力问题.应用Fourier变换求解该问题的控制偏微分方程,考虑了荷载的移动速度及振动频率对多孔饱和地基动力响应的影响,重点研究了移动速度达地基表面波速时多孔饱和半空间的振动问题(马赫效应),并与相应的弹性介质的解答进行了比较.结果显示当移动速度与多孔饱和半空间的表面波速相近时,地基会产生很大的振动;当移动速度大于表面波速时,多孔饱和半空间的动力响应与弹性半空间的动力响应有较大的差别.     

A parallel unstructured dynamic mesh adaptation algorithm for 3‐D unsteady flows

An unstructured dynamic mesh adaptation and load balancing algorithm has been developed for the efficient simulation of three‐dimensional unsteady inviscid flows on parallel machines. The numerical scheme was based on a cell‐centred finite‐volume method and the Roe's flux‐difference splitting. Second‐order accuracy was achieved in time by using an implicit Jacobi/Gauss–Seidel iteration. The resolution of time‐dependent solutions was enhanced by adopting an h‐refinement/coarsening algorithm. Parallelization and load balancing were concurrently achieved on the adaptive dynamic meshes for computational speed‐up and efficient memory redistribution. A new tree data structure for boundary faces was developed for the continuous transfer of the communication data across the parallel subdomain boundary. The parallel efficiency was validated by applying the present method to an unsteady shock‐tube problem. The flows around oscillating NACA0012 wing and F‐5 wing were also calculated for the numerical verification of the present dynamic mesh adaptation and load balancing algorithm. Copyright © 2005 John Wiley & Sons, Ltd.