二维矢量边界推进非结构网格生成方法

提出二维矢量边界推进生成非结构三角形网格方法并证明其可行性.根据流场边界尺度布置边界节点并运用矢量边界推进方法生成背景网格,运用符号面积函数和概率筛选方法布置初始点阵,提出Spring-Laplace方法优化节点位置,同时利用边交换技术优化网格结构.该方法可包含任意点源、线源和内嵌边界,可自由进行局部自适应加密或稀疏,实现任意平面域内与尺度要求一致的高效光滑三角网格剖分.     

Anti-diffusion method for interface steepening in two-phase incompressible flow

In this paper, we present a method for obtaining sharp interfaces in two-phase incompressible flows by an anti-diffusion correction, that is applicable in a straight-forward fashion for the improvement of two-phase flow solution schemes typically employed in practical applications. The underlying discretization is based on the volume-of-fluid (VOF) interface-capturing method on unstructured meshes. The key idea is to steepen the interface, independently of the underlying volume-fraction transport equation, by solving a diffusion equation with reverse time, i.e. an anti-diffusion equation, after each advection time step of the volume fraction. As the solution of the anti-diffusion equation requires regularization, a limiter based on the directional derivative is developed for calculating the gradient of the volume fraction. This limiter ensures the boundedness of the volume fraction. In order to control the amount of anti-diffusion introduced by the correction algorithm we propose a suitable stopping criterion for interface steepening. The formulation of the limiter and the algorithm for solving the anti-diffusion equation are applicable to 3-dimensional unstructured meshes. Validation computations are performed for passive advection of an interface, for 2-dimensional and 3-dimensional rising-bubbles, and for a rising drop in a periodically constricted channel. The results demonstrate that sharp interfaces can be recovered reliably. They show that the accuracy is similar to or even better than that of level-set methods using comparable discretizations for the flow and the level-set evolution. Also, we observe a good agreement with experimental results for the rising drop where proper interface evolution requires accurate mass conservation.     

Hybrid mesh generation with embedded surfaces using a multiple marching direction approach

This paper proposes a method for the creation of hybrid meshes with embedded surfaces for viscous flow simulations as an extension of the multiple marching direction approach (AIAA J. 2007; 45 (1):162–167). The multiple marching direction approach enables to place semi‐structured elements around singular points, where valid semi‐structured elements cannot be placed using conventional hybrid mesh generation methods. This feature is discussed first with a couple of examples. Elements sometimes need to be clustered inside a computational domain to obtain more accurate results. For example, solution features, such as shocks, vortex cores and wake regions, can be extracted during the process of adaptive mesh generation. These features can be represented as surface meshes embedded in a computational domain. Semi‐structured elements can be placed around the embedded surface meshes using the multiple marching direction approach with a pretreatment method. Tetrahedral elements can be placed easily instead. A couple of results are presented to demonstrate the capability of the mesh generation method. Copyright © 2008 John Wiley & Sons, Ltd.     

Surface Extension in Layered Structures with the Use of Metal Meshes for Heat-Transfer Enhancement

The properties of multilayer metal-fiber coverings superimposed on plane base layers, which consist of various combinations of a copper mesh and different-length copper wires, are examined. Coverings from discrete copper microwires, two-layer coverings including a mesh and microwires, and three-layer coverings consisting of a mesh and two layers of microwire are investigated. The contacts between the fibers, the copper mesh, and the base layer are analyzed by using the metallography technique. Quantitative estimates for these contacts are presented. The side effects connected with the hydrogen disease of copper are discussed. The thermal properties of the porous structures are analyzed with the use of a termographic chamber. It is shown that the porous metal-fiber structures have good mechanical and heat-transfer properties.     

钝体后湍流预混燃烧的PDF模拟

本文采用PDF方法对矩形燃烧室内钝体后的湍流预混火焰进行了数值模拟。脉动速度-频率-标量联合的PDF输运方程用Monte Carlo方法求解,质量、动量和能量的平均值由基于无结构网格的有限体积法求解,压力通过状态方程获得。PDF方程中所需的平均密度、平均速度和压力由有限体积法提供,并将用Monte Carlo方法求出的雷诺应力、化学反应源项和比热比传递给有限体积法。本文对丙烷和空气燃烧的不同简化化学反应机理进行了研究,并与实验结果进行比较,获得满意的结果。     

A conservative unstructured scheme for rapidly varied flows

This article introduces a new semi‐implicit, staggered finite volume scheme on unstructured meshes for the modelling of rapidly varied shallow water flows. Rapidly varied flows occur in the inundation of dry land during flooding situations. They typically involve bores and hydraulic jumps after obstacles such as road banks. Near such sudden flow transitions, the grid resolution is often low compared with the gradients of the bathymetry. Locally the hydrostatic pressure assumption may become invalid. In these situations, it is crucial to apply the correct conservation properties to obtain accurate results. An important feature of this scheme is therefore its ability to conserve momentum locally or, by choice, preserve constant energy head along a streamline. This is achieved using a special interpolation method and control volumes for momentum. The efficiency of inundation calculations with locally very high velocities, and in the case of unstructured meshes locally very small grid distances, is severely hampered by the Courant condition. This article provides a solution in the form of a locally implicit time integration for the advective terms that allows for an explicit calculation in most of the domain, while maintaining unconditional stability by implicit calculations only where necessary. The complex geometry of flooded urban areas asks for the flexibility of unstructured meshes. The efficient calculation of the pressure gradient in this, and other semi‐implicit staggered schemes, requires, however, an orthogonality condition to be put on the grid. In this article a simple method is introduced to generate unstructured hybrid meshes that fulfil this requirement. Copyright © 2008 John Wiley & Sons, Ltd.     

Two splitting schemes for nonstationary convection-diffusion problems on tetrahedral meshes

Two splitting schemes are proposed for the numerical solution of three-dimensional nonstationary convection-diffusion problems on unstructured meshes in the case of a full diffusion tensor. An advantage of the first scheme is that splitting is generated by the properties of the approximation spaces and does not reduce the order of accuracy. An advantage of the second scheme is that the resulting numerical solutions are nonnegative. A numerical study is conducted to compare the splitting schemes with classical methods, such as finite elements and mixed finite elements. The numerical results show that the splitting schemes are characterized by low dissipation, high-order accuracy, and versatility.     

Evaluation of adaptive mesh refinement and coarsening for the computation of compressible flows on unstructured meshes

The compressible gas flows of interest to aerospace applications often involve situations where shock and expansion waves are present. Decreasing the characteristic dimension of the computational cells in the vicinity of shock waves improves the quality of the computed flows. This reduction in size may be accomplished by the use of mesh adaption procedures. In this paper an analysis is presented of an adaptive mesh scheme developed for an unstructured mesh finite volume upwind computer code. This scheme is tailored to refine or coarsen the computational mesh where gradients of the flow properties are respectively high or low. The refinement and coarsening procedures are applied to the classical gas dynamic problems of the stabilization of shock waves by solid bodies. In particular, situations where oblique shock waves interact with an expansion fan and where bow shocks arise around solid bodies are considered. The effectiveness of the scheme in reducing the computational time, while increasing the solution accuracy, is assessed. It is shown that the refinement procedure alone leads to a number of computational cells which is 20% larger than when alternate passes of refinement and coarsening are used. Accordingly, a reduction of computational time of the same order of magnitude is obtained. Copyright © 2005 John Wiley & Sons, Ltd.     

美式期权定价问题的变网格差分方法

本文提出一种求解美式期权定价自由边值问题的变网格差分方法.通过建立一个自由边界所满足的方程,利用变网格技术可同时求出期权的差分解和最佳执行边界.本文分别讨论了显式和隐式变网格差分格式,并给出了差分解的收敛性和稳定性分析.数值实验表明本文算法是一个非常有效的期权定价算法.