An adaptive cut‐cell method for environmental fluid mechanics

In this work we present a numerical method for solving the incompressible Navier–Stokes equations in an environmental fluid mechanics context. The method is designed for the study of environmental flows that are multiscale, incompressible, variable‐density, and within arbitrarily complex and possibly anisotropic domains. The method is new because in this context we couple the embedded‐boundary (or cut‐cell) method for complex geometry with block‐structured adaptive mesh refinement (AMR) while maintaining conservation and second‐order accuracy. The accurate simulation of variable‐density fluids necessitates special care in formulating projection methods. This variable‐density formulation is well known for incompressible flows in unit‐aspect ratio domains, without AMR, and without complex geometry, but here we carefully present a new method that addresses the intersection of these issues. The methodology is based on a second‐order‐accurate projection method with high‐order‐accurate Godunov finite‐differencing, including slope limiting and a stable differencing of the nonlinear convection terms. The finite‐volume AMR discretizations are based on two‐way flux matching at refinement boundaries to obtain a conservative method that is second‐order accurate in solution error. The control volumes are formed by the intersection of the irregular embedded boundary with Cartesian grid cells. Unlike typical discretization methods, these control volumes naturally fit within parallelizable, disjoint‐block data structures, and permit dynamic AMR coarsening and refinement as the simulation progresses. We present two‐ and three‐dimensional numerical examples to illustrate the accuracy of the method. Copyright © 2008 John Wiley & Sons, Ltd.     

A moving finite element method for the population balance equation

A moving finite element algorithm has been compared against the upwind-differencing and Smolarkiewicz methods for the population balance equation of multicomponent particle growth processes. Analytical solutions and an error function have been used to test the numerical methods. The moving finite elements technique is much more accurate than other methods for a wide range of parameters. Since this method uses moving grids, it is able to model very narrow particle size distributions. It is also shown that the method can be extended to solve condensational growth problems which include particle curvature and non-continuum mass transfer effects.     

On a graded mesh method for a class of weakly singular Volterra integral equations

In this paper a class of weakly singular Volterra integral equations with an infinite set of solutions is investigated. Among the set of solutions only one particular solution is smooth and all others are singular at the origin. The numerical solution of this class of equations has been a difficult topic to analyze and has received much previous investigation. The aim of this paper is to improve the convergence rates by a graded mesh method. The convergence rates are proved and a variety of numerical examples are provided to support the theoretical findings.     

不可压缩流动的移动粒子半隐式法数值模拟

介绍了一种新的移动粒子半隐式方法。该方法采用Lagrange描述,用一系列离散粒子代替流体,通过半隐式求解和时间推进法来进行计算,弥补了传统网格法对复杂形状、大变形、高速撞击等情况下网格划分和重构过于繁琐的缺点,在模拟不可压缩流动问题中得到了较好的效果。     

FV/MC混合算法求解轴对称钝体后湍流流场

介绍一种有限容积／Monte Carlo结合求解湍流流场的相容的混合算法．有限容积法求解Reynolds平均的动量方程和能量方程，Monte Carlo方法求解模化的脉动速度—频率—标量联合的PDF方程．将该算法发展到无结构网格，探讨了在无结构网格中实现两种方法的耦合，包括颗粒定位，颗粒场和平均场之间数据交换等问题．并以二维轴对称钝体后湍流流场作为算例，比较了计算结果与实验结果．     

轴对称高速碰撞问题的拉格朗日无结构网格有限体积法的并行计算

给出计算轴对称高速碰撞问题的拉格朗日无结构三角形网格有限体积法的并行格式，并给出以小巨型机ＡｌｉａｎｔＦＸ／４０为目标计算机的算例数值模拟结果和效率分析     

用有限元自适网格控制结构应力分析的精度

本文涉及有限元结构应务分析精度控制的三个基本方面：１）自适应的有限元网格自动生成技术；２）有限元应力分析的误差估计：３）用单元尺寸场将误差分布用于形成新网格。     

Investigating droplet separation efficiency in wire-mesh mist eliminators in bubble column

Effects of design parameters on performance of wire-mesh mist eliminators were experimentally investigated in 15 cm bubble column. The demisters performances were evaluated by droplet collection efficiency as a function of wide ranges of operating and design parameters. These parameters include: droplet size exiting the demister (250–380 μm), specific surface area (236–868 m²/m³), void fraction (97–98.3%), wire diameter (0.14–0.28 mm), packing density (130–240 kg/m³), and superficial gas velocity (0.109–0.118 m/s. All demisters were 15 cm in diameter with 10 cm pad thickness, made from 316L stainless steel layered type demister pad wires. Experiments were carried out using air–water system at ambient temperature and atmospheric pressure. The experimental data on the droplet removal efficiency were obtained using Malvern Laser Droplet Sizer. The removal efficiency was found to increase with the increasing the demister specific surface area, packing density, and superficial gas velocity. In contrast, the removal efficiency was found to increase with decreasing the demister void fraction and wire diameter. The separation efficiency is correlated empirically as a function of the design parameters. A good agreement was obtained between the measured values and the correlation predictions with ±5% accuracy.     

带插值的自适应网格重构算法求解带移动热源的反应扩散方程的理论分析

本文研究一个带插值的网格重构算法求解一类带移动热源的反应扩散方程. 算法包括两步: 第一步是用旧时间网层上的计算解计算新时间层上的空间网格; 第二步是使用有限差分方法在新时间层 空间网格上离散方程, 并且将旧时间层上计算解的插值作为初始值. 对于时间, 我们获得了一阶收敛结果. 对于空间, 我们证明了使用线性插值算法的一阶收敛性和使用二次插值算法的二阶收敛性. 数值例子肯定了本文的理论结果.     

A nonlinear PDE model for reconstructing a regular surface from sampled data using a level set formulation on triangular meshes

In this paper, we propose a nonlinear PDE model for reconstructing a regular surface from sampled data. At first, we show the existence and the uniqueness of a viscosity solution to this problem. Then we propose a numerical scheme for solving the nonlinear level set equation on unstructured triangulations adapted to the data sample. We show the consistency of this scheme. In addition, we show how to compute nodewise first and second order derivatives. Some application examples of curve or surface construction are provided to illustrate the potential and to demonstrate the accuracy of this method.