A PRESSURE-BASED EULER SCHEME FOR TRANSONIC INTERNAL AND EXTERNAL FLOW SIMULATION

SUMMARY

The development of a pressure-based Euler scheme, based on a collocated grid arrangement and solving for Cartesian mass flux components is described. An implicit numerical dissipation model, which includes second and fourth difference terms expressed in pressure, is demonstrated. A smoothing function, which includes a two-level filter, is used to adjust the second order-dissipation. The successful calculations for the inviscid transonic aerodynamic flow around airfoil NACA0012 are presented. The transonic quasi-one-dimensional flow calculation in a converging-diverging nozzle is chosen as an example of an internal flow simulation.     

A Hybrid Domain Decomposition and Multigrid Method for the Acceleration of Compressible Viscous Flow Calculations on Unstructured Triangular Meshes

This paper is concerned with the formulation and the evaluation of a hybrid solution method that makes use of domain decomposition and multigrid principles for the calculation of two-dimensional compressible viscous flows on unstructured triangular meshes. More precisely, a non-overlapping additive domain decomposition method is used to coordinate concurrent subdomain solutions with a multigrid method. This hybrid method is developed in the context of a flow solver for the Navier-Stokes equations which is based on a combined finite element/finite volume formulation on unstructured triangular meshes. Time integration of the resulting semi-discrete equations is performed using a linearized backward Euler implicit scheme. As a result, each pseudo time step requires the solution of a sparse linear system. In this study, a non-overlapping domain decomposition algorithm is used for advancing the solution at each implicit time step. Algebraically, the Schwarz algorithm is equivalent to a Jacobi iteration on a linear system whose matrix has a block structure. A substructuring technique can be applied to this matrix in order to obtain a fully implicit scheme in terms of interface unknowns. In the present approach, the interface unknowns are numerical fluxes. The interface system is solved by means of a full GMRES method. Here, the local system solves that are induced by matrix-vector products with the interface operator, are performed using a multigrid by volume agglomeration method. The resulting hybrid domain decomposition and multigrid solver is applied to the computation of several steady flows around a geometry of NACA0012 airfoil.     

Chimera technique for transporting disturbances

The Chimera technique for moving grids is used to take into account nonhomogeneous unsteady inflow conditions in the simulation of aerodynamic flows. The method is applied to simulate the transport of a large‐scale vortex by a mean velocity field over a large distance, where it finally interacts with an airfoil. The Chimera approach allows one to resolve the vortex on a fine grid, whereas the unstructured background grid covering most of the computational domain can be much coarser. This method shows the same low numerical dissipation as a simulation on a globally fine grid. Several precursor tests are performed with a finite modified analytical Lamb–Oseen type vortex to study the influence of spatial and temporal resolution and the employed numerical scheme. Then, the interaction of an analytical vortex with a NACA0012 airfoil and with an ONERA‐A airfoil near stall is studied. Finally, a realistic vortex is generated by a ramping airfoil and is transported on a moving Chimera block and then interacts with a two‐element airfoil, which allows one to simulate a typical setup for a gust generator in aerodynamic facilities. Copyright © 2012 John Wiley & Sons, Ltd.     

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

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

翼型、机翼非定常运动模式下动态数值模拟

基于RANS方程,通过刚性动网格技术实现对翼型和机翼典型运动模式的描述,采用双时间推进方法和Roe空间离散格式对流场求解,构建了一个非定常气动计算平台;以NACA0012翼型为算倒进行了动态数值模拟可信度验证。数值模拟结果与试验数据吻合较好,升力和俯仰力矩的最大计算误差分别为3％和10％,表明了该平台的可靠性。另外,还数值模拟了M6机翼的动态非定常流场,并分析了两种湍流模型对非定常流场激波的捕捉能力。结果表明非定常流动中S-A湍流模型对激波的捕捉较B-L模型更敏感。文中开发的非定常计算平台对进一步解决三维复杂流场的流动问题有很高的工程应用价值。     

THE APPLICATION OF A FINITE DIFFERENCE/FINITE VOLUME ALGORITHM TO SOLVE MAXWELL'S EQUATIONS IN THE TIME DOMAIN

Abstract

A finite volume/finite difference method based on Ni's multigrid formulation is introduced for the solution of Maxwell's equations. The scheme is presented for the cases of transverse magnetic scattering from two-dimensional circular and square cylinders, as well as from NACA 0012 airfoil. The codes are validated against the traditional Method of Moments, which is analogous to a panel method in CFD. The circular cylinder scattering is compared to the analytical series solution for better understanding how the roles of numerical dispersion and dissipation errors affect the solution. The reflecting boundary conditions are modeled by the idea of inducing fields inside the conductor and a method of modeling the singularities that arise at a sharp corner is presented. Absorbing boundary conditions are modeled by integrating along the characteristic compatibility equations in the direction of the outgoing wave.     

基于Newton/Gauss-Seidel迭代的DGM隐式方法

在Newton迭代方法的基础上,对高阶精度间断Galerkin有限元方法(DGM)的时间隐式格式进行了研究. Newton迭代 法的优势在于收敛效率高效,并且定常和非定常问题能够统一处理,对于非定常问题无需引入双时间步策略. 为了避免大型矩阵的求逆,采用一步Gauss-Seidel迭代和Matrix-free技术消去残值Jacobi矩阵的上、下三角矩阵,从而只需计算和存储对角(块)矩阵. 对角(块)矩阵采用数值方法计算. 空间离散采用Taylor基,其优势在于对于任意形状的网格,基函数的形式是一致的,有利于在混合网格上推广. 利用该方法,数值模拟了Bump绕流和NACA0012翼型绕流. 计算结果表明,与显式的Runge-Kutta时间格式相比,隐式格式所需的迭代步数和CPU时间均在很大程度上得到减少,计算效率能够提高1～ 2个量级.     

Application of hybrid Cartesian grid and gridless approach to moving boundary flow problems

The paper presents a hybrid Cartesian grid and gridless approach to solve unsteady moving boundary flow problems. Unlike the Chimera clouds of points approach, the hybrid approach uses a Cartesian grid to cover most of the computational domain and a gridless method to calculate a relatively small region adjacent to the body surface, making use of the flexibility of the gridless method in handling surface grid with complicated geometry and the computational efficiency of the Cartesian grid. Four cases were conducted to examine the applicability, accuracy and robustness of the hybrid approach. Steady flows over a single NACA0012 airfoil and dual NACA0012 airfoils at different Mach numbers and angles of attack were simulated. Moreover, by implementing a dynamic hole cutting, node identification and information communication between the Cartesian grid and the gridless regions, unsteady flows over a pitching NACA0012 airfoil (small displacement) and two‐dimensional airfoil/store separation (large displacement) were performed. The computational results were found to agree well with earlier experimental data as well as computational results. Shock waves were accurately captured. The computational results show that the hybrid approach is of potential to solve the moving boundary flow problems. Copyright © 2013 John Wiley & Sons, Ltd.     

一种改进混合迎风格式在翼型流场激波捕捉中的应用

波阻是飞行器超音速飞行的关键设计因素,精确捕捉激波在流场中的位置,是数值模拟含激波流场和精确计算波阻的一个重要研究内容.本文基于网格节点有限体积空间离散方法,采用AUSM格式与FVS格式的混合格式(MAUSM方法)计算对流通量,从而抑制在数值模拟流场出现的激波处振荡和过冲现象,确保AUSM准确捕获接触间断的特性和FVS格式捕捉激波的能力.本文使用MAUSM方法分别计算了在跨声速和超声速条件下的NACA0012翼型流场,并与中心差分格式的计算结果进行比较.结果表明,对于存在激波的翼型流场,MAUSM方法是有效的.     

Numerical solution of the unsteady Euler equations for airfoils using approximate boundary conditions

This paper presents an efficient numerical method for solving the unsteady Euler equations on stationary rectilinear grids. Boundary conditions on the surface of an airfoil are implemented by using their first-order expansions on the mean chord line. The method is not restricted to flows with small disturbances since there are no restrictions on the mean angle of attack of the airfoil. The mathematical formulation and the numerical implementation of the wall boundary conditions in a fully implicit time-accurate finite-volume Euler scheme are described. Unsteady transonic flows about an oscillating NACA 0012 airfoil are calculated. Computational results compare well with Euler solutions by the full boundary conditions on a body-fitted curvilinear grid and published experimental data. This study establishes the feasibility for computing unsteady fluid-structure interaction problems, where the use of a stationary rectilinear grid offers substantial advantages in saving computer time and program design since it does not require the generation and implementation of time-dependent body-fitted grids.     

A sensitivity equation method for fast evaluation of nearby flows and uncertainty analysis for shape parameters

This paper presents the application of the continuous sensitivity equation method (CSEM) to fast evaluation of nearby flows and to uncertainty analysis for shape parameters. The flow and sensitivity fields are solved using an adaptive finite-element method. A new approach is presented to extract accurate flow derivatives at the boundary, which are needed in the shape sensitivity boundary conditions. Boundary derivatives are evaluated via high order Taylor series expansions used in a constrained least-squares procedure. The proposed method is first applied to fast evaluation of nearby flows: the baseline flow and sensitivity fields around a NACA 0012 airfoil are used to predict the flow around airfoils with nearby shapes obtained by modifications of the thickness (NACA 0015), the angle of attack and the camber (NACA 4512). The method is then applied to evaluate the influence of geometrical uncertainties on the flow around a NACA 0012 airfoil.     

动态混合网格生成及隐式非定常计算方法

建立了一种基于动态混合网格的非定常数值计算方法. 混合网格由贴体的四边形网格、外场 的多层次矩形网格和中间的三角形网格构成. 当物体运动时，贴体四边形网格随物体运动而 运动，而外场的矩形网格保持静止，中间的三角形网格随之变形；当物体运动位移较大，导 致三角形网格的质量降低，甚至导致网格相交时，在局部重新生成网格. 新网格上的物理量 由旧网格上的物理量插值而得. 为了提高计算效率，采用了双时间步和子迭代相结合的隐式 有限体积格式计算非定常Navier-Stokes方程. 子迭代采用高效的块LU-SGS方法. 利用该 方法数值模拟了NACA0012振荡翼型的无黏和黏性绕流，得到了与实验和他人计算相当一致 的结果.     

Experimental and numerical investigation on sound generation from airfoil-flow interaction

Aerodynamic noise due to interaction between incoming turbulence and rotating blades is an important component in the wind turbine noise. The rod-airfoil configuration is used to investigate the interactive phenomenon experimentally and numerically. Distribution of unsteady pressure on the airfoil surface is measured for different rod positions and airfoil attack angles. Two National Advisory Committee for Aeronatics (NACA) airfoils, NACA0012 and NACA0018, and two wind turbine airfoils, S809 and S825 are investigated. In addition, for low angles of attack, the flow field around the airfoil's leading edge is investigated with the particle image velocimetry (PIV). The experimental results indicate that unsteady pressure disturbances on the airfoil surface are related to the rod vortex shedding. Meanwhile, the interaction flow field of the rod and NACA0012 airfoil is simulated with the unsteady Reynolds averaged Navier-Stokes method (URANS), and the obtained pressure spectra are compared with the experimental results.     

基于自适应直角网格的二维全速势方程有限体积解法A finite volume method for 2D Full-Potential equation on adaptive Cartesian grids

为满足亚声速和跨声速飞机概念设计中快速气动计算的需求,研究和发展一种基于自适应直角网格的非线性全速势方程有限体积解法。要点如下。（1）在几何自适应直角网格的基础上,使用结合单元融合的网格切割算法处理物面边界,提出一种修正非贴体切割网格的方法。（2）采用隐式格式结合GM RES算法求解该非线性位流方程,针对流场的自适应来捕捉激波。（3）采用镜像法处理物面边界处的无穿透条件,并提出解析的方法来修正镜像单元的值。（4）针对直角网格的特点,提出在库塔线上插入库塔单元的方法施加库塔条件。NACA0012翼型绕流的算例结果表明,该方法用于亚声速和跨声速气动计算能得到令人满意的结果,且自动化程度高、收敛速度快。     

Numerical simulation of compressible turbulent flow via improved gas‐kinetic BGK scheme

In order to solve compressible turbulent flow problems, this study focuses on incorporating the Spalart–Allmaras turbulence model into gas‐kinetic BGK (Bhatnagar–Gross–Krook) scheme. The Spalart–Allmaras turbulence model is solved using finite difference discretization. The variables on the cell interface are interpolated via the van Leer limiter in the reconstruction stage. Simulation of subsonic and transonic flow over a NACA0012 airfoil has been implemented using two‐dimensional body‐fitted grids. The numerical results obtained appear in good agreement with the AGARD results, demonstrating the effectiveness and usefulness of the strategy of coupling the Spalart–Allmaras turbulence model with the BGK scheme for compressible turbulent flow simulation. Copyright © 2010 John Wiley & Sons, Ltd.     

EXPERIMENTAL INVESTIGATION ON LIFT INCREMENT OF A PLASMA CIRCULATION CONTROL AIRFOIL

利用等离子体激励器发展了新型的环量增升技术,并对二维NACA0012翼型绕流实施控制。由于NACA0012翼型为尖后缘构型,环量增升装置由2个非对称型介质阻挡放电等离子体激励器构成。一个等离子体激励器贴附于翼型吸力面靠近后缘处,其诱导的壁面射流沿来流方向指向下游;另一个等离子体激励器贴附于翼型压力面靠近后缘处,其诱导的壁面射流与来流方向相反指向上游。在风洞中通过时间解析二维PIV系统对翼型绕流流场进行了测量,基于翼型弦长的雷诺数Re=20 000。结果表明在等离子体激励器的控制下,翼型压力面靠近后缘处可以形成一个定常回流区,从而起到虚拟气动外形的作用,因此翼型吸力面的流场得到加速,压力面的流场得到减速,使得翼型压力面的吸力以及压力面的压力都得到增加,进而增加了翼型的环量。风洞天平测力实验进一步验证了该环量增升技术的有效性。在整个攻角范围内,施加控制的翼型的升力系数相比没有控制的工况有明显的提高。     

非结构网格上p型多重网格法流场数值模拟

在二维、三维非结构网榕上,针对间断Galerkin方法计算量大、收敛慢的缺点将p型多重网格方法应用于该方法求解跨音速Euler方程,提高计算效率。p型多重网格方法是通过对不同阶次多项式近似解进行递归迭代求解,来达到加速收敛。文中对高阶近似(p＞0)使用显式格式,最低阶近似(p=0)采用隐式格式。NACA0012翼型和O...     

Low-diffusion preconditioning scheme with linear cut-off parameter

A low-diffusion preconditioning Roe scheme with an adjustable parameter to control the numerical dissipation is proposed. This scheme reflects the real physical dissipation in the extremely low-speed region. The preconditioning parameter in scheme is improved by linear cut-off and correction factor. The numerical results of low-Mach-number/low-Reynolds-number steady solutions of viscous flows past a circular cylinder and past a NACA0012 airfoil show the efficiency of the new scheme.     

Design optimization of unsteady airfoils with continuous adjoint method

In this paper, a new unsteady aerodynamic design method is presented based on the Navier-Stokes equations and a continuous adjoint approach. A basic framework of time-accurate unsteady airfoil optimization which adopts time-averaged aerodynamic coefficients as objective functions is presented. The time-accurate continuous adjoint equation and its boundary conditions are derived. The flow field and the adjoint equation are simulated numerically by the finite volume method (FVM). Feasibility and accuracy of the approach are perfectly validated by the design optimization results of the plunging NACA0012 airfoil.     

Large-eddy simulation of wing tip vortex in the near field

Large-eddy simulation with filtered-structure-function subgrid model and implicit large-eddy simulation (ILES without explicit subgrid model) using high-order accuracy and high resolution compact scheme have been performed on the tip vortex shedding from a rectangular half-wing with a NACA 0012 airfoil section and a rounded wing tip. The formation of the tip vortex and its initial development in the boundary layer and the near field wake are investigated and analysed in detail. The physics, why the tip vortex, which is originally turbulent in the boundary layer, is re-laminarised and becomes stable and laminar rapidly after shedding in the near field, is revealed by this simulation. The computation also shows the widely used second-order subgrid model is not consistent to six-order compact scheme and would degenerate the six-order LES results to second-order. Therefore, high-order schemes, grid refinement and six-order subgrid models are critical to LES approaches.