A parallel cell‐based DSMC method on unstructured adaptive meshes

A parallel DSMC method based on a cell‐based data structure is developed for the efficient simulation of rarefied gas flows on PC‐clusters. Parallel computation is made by decomposing the computational domain into several subdomains. Dynamic load balancing between processors is achieved based on the number of simulation particles and the number of cells allocated in each subdomain. Adjustment of cell size is also made through mesh adaptation for the improvement of solution accuracy and the efficient usage of meshes. Applications were made for a two‐dimensional supersonic leading‐edge flow, the axi‐symmetric Rothe's nozzle, and the open hollow cylinder flare flow for validation. It was found that the present method is an efficient tool for the simulation of rarefied gas flows on PC‐based parallel machines. Copyright © 2004 John Wiley & Sons, Ltd.     

Solution of the hyperbolic mild‐slope equation using the finite volume method

A finite volume solver for the 2D depth‐integrated harmonic hyperbolic formulation of the mild‐slope equation for wave propagation is presented and discussed. The solver is implemented on unstructured triangular meshes and the solution methodology is based upon a Godunov‐type second‐order finite volume scheme, whereby the numerical fluxes are computed using Roe's flux function. The eigensystem of the mild‐slope equations is derived and used for the construction of Roe's matrix. A formulation that updates the unknown variables in time implicitly is presented, which produces a more accurate and reliable scheme than hitherto available. Boundary conditions for different types of boundaries are also derived. The agreement of the computed results with analytical results for a range of wave propagation/transformation problems is very good, and the model is found to be virtually paraxiality‐free. Copyright © 2003 John Wiley & Sons, Ltd.     

A posteriori H1 error estimation for a singularly perturbed reaction diffusion problem on anisotropic meshes

The paper deals with a singularly perturbed reaction diffusionmodel problem. The focus is on reliable a posteriori error estimatorsfor the H¹ seminorm that can be applied to anisotropic finiteelement meshes. A residual error estimator and a local problemerror estimator are proposed and rigorously analysed. They arelocally equivalent, and both bound the error reliably. Threemodifications of these estimators are introduced and discussed. Much attention is given to the performance of the error estimatorin numerical experiments. This helps to identify those estimatorsthat are suitable for practical applications.     

Hybrid mesh adaptation applied to industrial numerical combustion

This paper presents an anisotropic mesh adaptation method applied to industrial combustion problems. The method is based on a measure of the distance between two Riemannian metrics called metric non‐conformity. This measure, which can be used to build a cost function to adapt meshes comprising several types of mesh elements, provides the basis for a generic mesh adaptation approach applicable to various types of physical problems governed by partial differential equations. The approach is shown to be applicable to industrial combustion problems, through the specification of a target metric computed as the intersection of several Hessian matrices reconstructed from the main variables of the governing equations. Numerical results show that the approach is cost effective in that it can drastically improve the prediction of temperature and species distributions in the flame region of a combustor while reducing computational cost. The results can be used as a basis for pollutant prediction models. Copyright © 2011 John Wiley & Sons, Ltd.     

A moving mesh algorithm for aero‐thermo‐mechanical modelling in turbomachinery

This paper describes the development of a mesh deformation method used for aero‐thermo‐mechanical coupling of turbo‐engine components. The method is based on the nonlinear solution of an elastic medium analogy, solved using finite element discretisation and modified to let the boundary nodes be free to slide over the deflected surfaces. This sliding technique relies on a B‐spline reconstruction of the moving boundary and increases the robustness of the method in situations where the boundary deflection field presents significant gradients or large relative motion between two distinct boundaries. The performance of the method is illustrated with the application to an interstage cavity of a turbine assembly, subjected to the deformations computed by a coupled thermo‐mechanical analysis of the engine component. Copyright © 2011 John Wiley & Sons, Ltd.     

具有非光滑边界层函数的线性变系数抛物型方程初边值问题的差分格式

本文利用非均匀网格和指数型拟合差分方法给出了具有非光滑边界层函数的线性抛物型方程关于小参数ε一致收敛的差分格式.文章还给出了误差估计和数值结果.     

ROBUSTNESS OF AN UPWIND FINITE DIFFERENCE SCHEME FOR SEMILINEAR CONVECTION-DIFFUSION PROBLEMS WITH BOUNDARY TURNING POINTS

We consider a singularly perturbed semilinear convection-diffusion problem with a boundary layer of attractive turning-point type. It is shown that its solution can be decomposed into a regular solution component and a layer component. This decomposi-tion is used to analyse the convergence of an upwinded finite difference scheme on Shishkin meshes.     

光学透明频率选择表面的设计研究

以氟化镁为基底材料, 采用基底、金属网栅与频率选择表面一体化设计方法设计了一种雷达波与光学波段双带通的结构. 利用模式匹配法对设计结构的传输特性进行了仿真研究, 并将设计结果与制备样件的测试结果进行了对比分析. 结果发现: 采用一体化设计的方法设计光学透明频率选择表面, 不仅能够快速得到电场基函数而且还能够准确预估其谐振尺寸, 从而在提高计算效率的同时避免了模式互作用零点的出现. 采用一体化设计方法获得了具有稳定滤波特性的光学透明频率选择表面, 为雷达/红外双模制导头罩的电磁屏蔽技术和隐身技术提供了一种有效技术方案. 关键词： 复合制导 频率选择表面 金属网栅 一体化设计     

基于变分原理的二维热传导方程差分格式

研究二维热传导方程的差分数值模拟.用变分原理在不规则结构网格上建立热流通量形式的差分格式.将热流通量作为未知函数求泛函极值,并与温度函数联立求解.克服通常九点格式用插值方法计算网格边界上的热传导系数和网格结点上的温度所引入的误差.