Generalized Wall Functions for High-speed Separated Flows Over Adiabatic and Isothermal Walls

We present the extension of our wall-laws developed for low-speed flows to super-and hyper-sonic configurations. In particular, we are interested in flows over isothermal walls and in the modeling of heat transfer. We recall the main steps of the development:

?Obtaining generalized wall functions for low-speed fluids, valid for all y ⁺.

?Taking into account transversal effects in wall-laws.

?Accounting for the compressible feature of the flow on adiabatic walls without using information on the local boundary layer structure, for instance its thickness, but only using information available at the fictitious wall.

?Taking into account thermal effects on isothermal walls. In particular, the heat flux at the real wall is obtained by an a posteriori evaluation using information at the fictitious one.

?Only using information available on unstructured meshes and avoiding the information coming from a Cartesian hypothesis for the mesh in near-wall regions.

These ingredients are validated on hyper-sonic configurations on adiabatic and isothermal walls for expansion and compression ramps.     

自然对流换热的高精度非结构网格有限体积法

用非结构网格有限体积法求解自然对流换热时,传统的对流项离散格式难以兼顾数值精度与计算效率,我们发展了一种耦合高精度格式的延迟修正方法,用于对流项的离散.高Re数下方腔驱动流数值计算验证了该方法具有较高的计算精度和较好的稳定性.Boussinesq流体的自然对流换热数值模拟,表明该方法能有效克服高Ra数时数值计算发散,可准确捕捉自然对流换热问题中不同偏心率下的等温线和流线分布特征.     

基于混合网格的低声爆反设计方法研究

低声爆设计方法已成为新一代军民用超声速飞机研制过程中必须解决的关键难题之一。针对传统SGD低声爆外形反设计方法无法对声爆近场非线性效应进行描述和分析的缺点,提出了利用CFD方法求解得到的声爆近场压力分布代替F函数进行低声爆反设计的方法。声爆近场预测采用点-点对接的结构/非结构混合网格,充分利用非结构网格对复杂外形适应性强和结构化网格计算效率高的优点。结果分析表明,基于改进后的低声爆反设计方法得到的方案在声爆超压以及感觉噪声级等方面都比基于原始SGD方法得到的方案有较大改善。     

三角形网格上多介质流动问题界面处理方法

多介质流动问题的求解一般是在结构网格上实现,而三角形网格对于复杂计算区域具有更好的适应性,本文结合rGFM方法,给出三角形网格上多介质流动问题界面处理方法.利用level-set方法跟踪界面,在界面处构造Riemann问题,得到界面处流体准确的流动状态.通过定义界面边界条件,将多介质流动问题转化为单介质流动问题,利用高精度RKDG方法求解.采用多个算例验证该方法的稳健性和有效性,结果表明该方法能准确捕捉界面和激波的位置,保持界面清晰.     

基于背景网格变形的动态网格移动方法

基于Delaunay背景网格插值技术的动态网格生成方法无需迭代计算,效率较高。但对复杂构形大幅运动的动边界问题,尤其当边界大幅转动时,背景网格极易交叉重叠。重新生成背景网格和重新定位网格节点信息不仅费时而且会导致网格质量的严重下降。本文提出改进的基于背景网格的动态网格变形方法,通过在初始Delaunay背景网格中添加辅助点,生成一层新的背景网格和新的映射关系;采用ball-vertex弹簧法驱动新背景网格的变形,进而牵动目标网格的变形。算例表明,本文提出的动态网格变形方法对所关心区域的网格具有良好保形性,边界可转动更大角度而不会出现网格交叉重叠问题,总体上提高了动态网格更新的效率和质量。     

二阶椭圆问题的各向异性非协调Crouzeix-Raviart型有限元方法

对满足最大角条件和坐标系条件的二维区域中的各向异性一般三角形网格,研究了二阶椭圆问题的非协调Crouzeix-Raviart型线性三角形有限元逼近,得到了最优的能量模和L2-模误差估计结果．     

Comparisons of some difference forms for compressible flow in cylindrical geometry on arbitrary Lagrangian and Eulerian framework

The study of cylindrically symmetric compressible fluid is interesting from both theoretical and numerical points of view. In this paper, the typical spherical symmetry properties of the numerical schemes are discussed, and an area weighted scheme is extended from a Lagrangian method to an arbitrary Lagrangian and Eulerian (ALE) method. Numerical results are presented to compare three discrete configurations, i.e., the control volume scheme, the area weighted scheme, and the plane scheme with the addition of a geometrical source. The fact that the singularity arises from the geometrical source term in the plane scheme is illustrated. A suggestion for choosing the discrete formulation is given when the strong shock wave problems are simulated.     

A posteriori error estimates for optimal control problems constrained by convection-diffusion equations

We propose a characteristic finite element discretization of evolutionary type convection-diffusion optimal control problems. Nondivergence-free velocity fields and bilateral inequality control constraints are handled. Then some residual type a posteriori error estimates are analyzed for the approximations of the control, the state, and the adjoint state. Based on the derived error estimators, we use them as error indicators in developing efficient multi-set adaptive meshes characteristic finite element algorithm for such optimal control problems. Finally, one numerical example is given to check the feasibility and validity of multi-set adaptive meshes refinements.     

Assessment of a discontinuous Galerkin method for the simulation of vortical flows at high Reynolds number

This paper focuses on the assessment of a discontinuous Galerkin method for the simulation of vortical flows at high Reynolds number. The Taylor–Green vortex at Re = 1600 is considered. The results are compared with those obtained using a pseudo‐spectral solver, converged on a 512³ grid and taken as the reference. The temporal evolution of the dissipation rate, visualisations of the vortical structures and the kinetic energy spectrum at the instant of maximal dissipation are compared to assess the results. At an effective resolution of 288³, the fourth‐order accurate discontinuous Galerkin method (DGM) solution (p = 3) is already very close to the pseudo‐spectral reference; the error on the dissipation rate is then essentially less than a percent, and the vorticity contours at times around the dissipation peak overlap everywhere. At a resolution of 384³, the solutions are indistinguishable. Then, an order convergence study is performed on the slightly under‐resolved grid (resolution of 192³). From the fourth order, the decrease of the error is no longer significant when going to a higher order. The fourth‐order DGM is also compared with an energy conserving fourth‐order finite difference method (FD4). The results show that, for the same number of DOF and the same order of accuracy, the errors of the DGM computation are significantly smaller. In particular, it takes 768³ DOF to converge the FD4 solution. Finally, the method is also successfully applied on unstructured high quality meshes. It is found that the dissipation rate captured is not significantly impacted by the element type. However, the element type impacts the energy spectrum in the large wavenumber range and thus the small vortical structures. In particular, at the same resolution, the results obtained using a tetrahedral mesh are much noisier than those obtained using a hexahedral mesh. Those obtained using a prismatic mesh are already much better, yet still slightly noisier. Copyright © 2013 John Wiley & Sons, Ltd.