On a fast direct elliptic solver by a modified Fourier method

We describe high order numerical algorithms for the solution of second order elliptic equations in rectangular domains. These algorithms are based on the Fourier method in combination with a subtraction procedure. The singularities at the corner points, arising due to non-smoothness of the boundaries, are treated explicitly using properly constructed singular corner functions. The present algorithm is a generalization of the Fast Poisson Solver developed in our previous paper. This revised version was published online in August 2006 with corrections to the Cover Date.     

A discontinuous Galerkin method/HLLC solver for the Euler equations

This paper proposes a fully three‐dimensional non‐linear Euler methodology for solving aerodynamic and acoustic problems in the presence of strong shocks and rarefactions. It uses a discontinuous Galerkin method (DGM) within the element, and a Riemann solver (HLLC) at the boundaries to propagate rarefactions while preserving the entropy condition and capturing shocks with no spurious oscillations. This approach is thought to marry the best aspects of finite element and finite volume methods, achieving conservation while not requiring the solution of a large matrix. Examples in which shock and rarefaction waves are well captured are presented and the propagation of acoustic pulses is well demonstrated. Copyright © 2003 John Wiley & Sons, Ltd.     

Numerical simulation on an interaction of a vortex street with an elliptical leading edge using an unstructured grid

An algorithm for a time accurate incompressible Navier–Stokes solver on an unstructured grid is presented. The algorithm uses a second order, three‐point, backward difference formula for the physical time marching. For each time step, a divergence free flow field is obtained based on an artificial compressibility method. An implicit method with a local time step is used to accelerate the convergence for the pseudotime iteration. To validate the code, an unsteady laminar flow over a circular cylinder at a Reynolds number of 200 is calculated. The results are compared with available experimental and numerical data and good agreements are achieved. Using the developed unsteady code, an interaction of a Karman vortex street with an elliptical leading edge is simulated. The incident Karman vortex street is generated by a circular cylinder located upstream. A clustering to the path of the vortices is achieved easily due to flexibility of an unstructured grid. Details of the interaction mechanism are analysed by investigating evolutions of vortices. Characteristics of the interactions are compared for large‐ and small‐scale vortex streets. Different patterns of the interaction are observed for those two vortex streets and the observation is in agreement with experiment. Copyright © 2004 John Wiley & Sons, Ltd.     

Development of a 3D Eulerian/Lagrangian Aircraft Icing Simulation Solver Based on OpenFOAM

A 3D icing simulation code is developed in the open-source CFD toolbox OpenFOAM. A hybrid Cartesian/body-fitted meshing method is used to generate high-quality meshes around complex ice shapes. Steady-state 3D Reynolds-averaged Navier-Stokes (RANS) equations are solved to provide the ensemble-averaged flow around the airfoil. Considering the multi-scale nature of droplet size distribution, and more importantly, to represent the less uniform nature of the Super-cooled Large Droplets (SLD), two droplet tracking methods are realized: the Eulerian method is used to track the small-size droplets (below 50 μm) for the sake of efficiency; the Lagrangian method with random sampling is used to track the large droplets (above 50 μm); the heat transfer of the surface overflow is solved on a virtual surface mesh; the ice accumulation is estimated via the Myers model; finally, the final ice shape is predicted by time marching. Limited by the availability of experimental data, validations are performed on 3D simulations of 2D geometries using the Eulerian and Lagrangian methods, respectively. The code proves to be feasible and accurate enough in predicting ice shapes. Finally, an icing simulation result of the M6 wing is presented to illustrate the full 3D capability.     

一种适用于磁流体力学切向间断的HLLC黎曼算子

磁场的存在使得磁流体力学特征波模不同于流体力学, 因此直接由流体力学HLLC黎曼算子导出的HLLC双中间态在交界的间断面会出现不守恒的问题。通常降级采用HLL单磁场中间态代替HLLC双磁场中间态以实现守恒和计算稳定, 代价是切向间断的模拟精度不足。本文对此进行改进, 在模拟切向间断时仍然保留原有的HLLC双磁场中间态, 同时各守恒量仍然能够满足Toro相容条件; 改进型HLLC算子在间断两侧的磁场分量存在差异, 因此能够更精确还原切向间断面。基于数值测试, 包括一维激波管和切向间断的时变模拟, 以及地球磁层三维数值模拟, 将模拟结果进行对比, 结果表明: 相比于已发展的HLLC算子, 改进型HLLC算子对切向间断具有更好的捕捉精度, 能够达到或接近耗时更多的HLLD算子的模拟精度。     

Analysis of a meshless solver for high Reynolds number flow

In this paper, the extension of an upwind least‐square based meshless solver to high Reynolds number flow is explored, and the properties of the meshless solver are analyzed both theoretically and numerically. Existing works have verified the meshless solver mostly with inviscid flows and low Reynolds number flows, and in this work, we are interested in the behavior of the meshless solver for high Reynolds number flow, especially in the near‐wall region. With both theoretical and numerical analysis, the effects of two parameters on the meshless solver are identified. The first one is the misalignment effect caused by the significantly skewed supporting points, and it is found that the meshless solver still yields accurate prediction. It is a very interesting property and is opposite to the median‐dual control volume based vertex‐centered finite volume method, which is known to give degraded result with stretched triangular/tetrahedral cells in the near‐wall region. The second parameter is the curvature, and according to theoretical analysis, it is found in the region with both large aspect ratio and curvature, and the streamwise residual is less affected; however, the wall‐normal counterpart suffers from accuracy degradation. In this paper, an improved method that uses a meshless solver for the streamwise residual and finite difference for wall‐normal residual is developed. This method is proved to be less sensitive to the curvature and provides improved accuracy. This work presents an understanding of the meshless solver for high Reynolds number flow computation, and the analysis in this paper is verified with a series of numerical experiments. Copyright © 2012 John Wiley & Sons, Ltd.     

A simplified adaptive Cartesian grid system for solving the 2D shallow water equations

This paper presents a new simplified grid system that provides local refinement and dynamic adaptation for solving the 2D shallow water equations (SWEs). Local refinement is realized by simply specifying different subdivision levels to the cells on a background uniform coarse grid that covers the computational domain. On such a non‐uniform grid, the structured property of a regular Cartesian mesh is maintained and neighbor information is determined by simple algebraic relationships, i.e. data structure becomes unnecessary. Dynamic grid adaptation is achieved by changing the subdivision level of a background cell. Therefore, grid generation and adaptation is greatly simplified and straightforward to implement. The new adaptive grid‐based SWE solver is tested by applying it to simulate three idealized test cases and promising results are obtained. The new grid system offers a simplified alternative to the existing approaches for providing adaptive mesh refinement in computational fluid dynamics. Copyright © 2011 John Wiley & Sons, Ltd.     

Roe‐type Riemann solvers for general hyperbolic systems

We present a Roe‐type weak formulation Riemann solver where the average coefficient matrix is computed numerically. The novelty of this approach is that it is general enough that can be applied to any hyperbolic system while retaining the accuracy of the original Roe solver. We show applications to the compressible Euler equations with general equation of state. An alternative version of the method uses directly the eigenvectors in the averaging process, simplifying the algorithm. These new solvers are applied in conservative and path‐conservative schemes with high‐order accuracy and on unstructured meshes. Copyright © 2014 John Wiley & Sons, Ltd.     

MFCAV近似Riemann解法器在相容拉氏方法中的熵条件分析

 对基于MFCAV(Multi Fluid Channel on Averaged Volume)近似Riemann解法器的相容拉氏方法的熵条件进行了分析. 结果表明与满足声学形式Riemann解法器的熵不同, 前者只能在每个网格边界左、右两侧网格的熵随时间变化的和保证大于零, 即能保证整体熵增, 但不保证传统意义上的在每个网格中的熵增;而后者不仅保证整体熵增, 而且还满足传统意义上的熵增. 因此MFCAV的熵增相对声学形式解法器而言要弱一些, 由此表明其熵增可能要小些, 使得格式的耗散可能要小些.数值算例也验证了分析的正确性.     

适用于等熵流动的交错拉氏Godunov方法

为消除传统单元中心型Godunov方法在求解稀疏波问题时的非物理过热现象,发展一种适用于等熵流动的交错拉氏Godunov方法.主要的特征是采用速度与热力学变量交错分布的形式,避免在单元内进行速度平均,从而消除由于动量平均过程导致的动能耗散.与传统的von Neumann型交错网格方法相比,网格的边界通量由节点处的多维黎曼求解器提供,克服了多维人工粘性选取带来的困难.为减少多维黎曼求解器在求解稀疏波问题时的非物理熵增,给出稀疏波出现的合理判据,从而保证了热力学关系式的满足.数值实验表明：该方法能很好地消除稀疏波的过热现象,同时在求解激波问题时又能保持与传统单元中心型拉氏方法相同的激波捕捉能力.