顺序输送管道水力瞬变模拟的有限体积法

针对顺序输送水力瞬变中输送介质不连续的情况,提出采用有限体积法取代传统的特征线法进行求解,以减少插值误差。首先得到水击基本方程基于压力-流速的非守恒有限体积离散格式,通过Rieman求解器得到控制体界面的Rieman解。采用MUSCL-Hancock Primitive方法进行界面数值重构与时间推进,构建时空二阶精度的TVD格式。边界的处理采用Rieman不变量与构造虚拟边界节点相结合的方法,使其与整体精度保持一致。数值计算与实验对比,证明了本文算法能有效抑制非物理振荡,具有精度高、鲁棒性好的特点。     

Numerical modeling of wave interactions with coastal structures by a constrained interpolation profile/immersed boundary method

A high‐order difference method based multiphase model is proposed to simulate nonlinear interactions between water wave and submerged coastal structures. The model is based on the Navier–Stokes equations using a constrained interpolation profile (CIP) method for the flow solver, and employs an immersed boundary method (IBM) for the treatment of wave–structure interactions. A more accurate interface capturing scheme, the volume of fluid/weighed line interface calculation (VOF/WLIC) scheme, is adopted as the interface capturing method. A series of computations are performed to verify the application of the model for simulations of fluid interaction with various structures. These problems include flow over a fixed cylinder, water entry of a circular cylinder and solitary waves passing various submerged coastal structures. Computations are compared with the available analytical, experimental and other numerical results and good agreement is obtained. The results of this study demonstrate the accuracy and applications of the proposed model to simulate the nonlinear flow phenomena and capture the complex free surface flow. Copyright © 2015 John Wiley & Sons, Ltd.     

A projection scheme for incompressible multiphase flow using adaptive Eulerian grid

This paper presents a finite element method for incompressible multiphase flows with capillary interfaces based on a (formally) second‐order projection scheme. The discretization is on a fixed Eulerian grid. The fluid phases are identified and advected using a level set function. The grid is temporarily adapted around the interfaces in order to maintain optimal interpolations accounting for the pressure jump and the discontinuity of the normal velocity derivatives. The least‐squares method for computing the curvature is used, combined with piecewise linear approximation to the interface. The time integration is based on a formally second order splitting scheme. The convection substep is integrated over an Eulerian grid using an explicit scheme. The remaining generalized Stokes problem is solved by means of a formally second order pressure‐stabilized projection scheme. The pressure boundary condition on the free interface is imposed in a strong form (pointwise) at the pressure‐computation substep. This allows capturing significant pressure jumps across the interface without creating spurious instabilities. This method is simple and efficient, as demonstrated by the numerical experiments on a wide range of free‐surface problems. Copyright © 2004 John Wiley & Sons, Ltd.     

Vortex‐in‐cell method combined with a boundary element method for incompressible viscous flow analysis

In this study, an immersed boundary vortex‐in‐cell (VIC) method for simulating the incompressible flow external to two‐dimensional and three‐dimensional bodies is presented. The vorticity transport equation, which is the governing equation of the VIC method, is represented in a Lagrangian form and solved by the vortex blob representation of the flow field. In the present scheme, the treatment of convection and diffusion is based on the classical fractional step algorithm. The rotational component of the velocity is obtained by solving Poisson's equation using an FFT method on a regular Cartesian grid, and the solenoidal component is determined from solving an integral equation using the panel method for the convection term, and the diffusion term is implemented by a particle strength exchange scheme. Both the no‐slip and no‐through flow conditions associated with the surface boundary condition are satisfied by diffusing vortex sheet and distributing singularities on the body, respectively. The present method is distinguished from other methods by the use of the panel method for the enforcement of the no‐through flow condition. The panel method completes making use of the immersed boundary nature inherent in the VIC method and can be also adopted for the calculation of the pressure field. The overall process is parallelized using message passing interface to manage the extensive computational load in the three‐dimensional flow simulations. Copyright © 2011 John Wiley & Sons, Ltd.     

Finite element analysis of flow and heat transfer with moving free surface using fixed grid system

A numerical study was performed on flow and heat transfer involving moving free surfaces that occurs in mold filling processes such as casting and injection molding. In these problems, the calculation domain changes continuously and the numerical treatment of the moving interface tends to cause artificial diffusion. Among the solution algorithms based on the Eulerian method, the volume-of-fluid (VOF) method was used because the method is simple and efficient in handling the complex flow patterns inside the cavity. To solve the transport equation of free surface without artificial smearing of the interface the baby-cell method was employed in the geometric reconstruction of the free surface. Furthermore, a predictor–corrector method was adopted in the time integration of volume-of-fluid (VOF) transport equation to increase the accuracy. The proposed scheme was verified through several benchmark problems. In order to show the capability of the proposed method, several three-dimensional mold filling processes were solved. The current algorithm was applied to the floating body problem. Three-dimensional floating body problems were tested.     

A fully conservative high‐order upwind multi‐moment method using moments in both upwind and downwind cells

We propose a fully conservative high‐order upwind multi‐moment method for the conservation equation. The proposed method is based on a third‐order polynomial interpolation function and semi‐Lagrangian formulation and is a variant of the constrained interpolation profile conservative semi‐Lagrangian scheme with third‐order polynomial function method. The third‐order interpolation function is constructed based on three constraints in the upwind cell (two boundary values and a cell average) and a constraint in the downwind cell (a cell center value). The proposed method shows fourth‐order accuracy in a benchmark problem (sine wave propagation). We also propose a less oscillatory formulation of the proposed method. The less oscillatory formulation can minimize numerical oscillations. These methods were validated through scalar transport problems, and compressible flow problems (shock tube and 2D explosion problems). Copyright © 2016 John Wiley & Sons, Ltd.     

Mathematical model and numerical method for spontaneous potential log in heterogeneous formations

This paper introduces a new spontaneous potential log model for the case in which formation resistivity is not piecewise constant. The spontaneous potential satisfies an elliptic boundary value problem with jump conditions on the interfaces. It has beer/ shown that the elliptic interface problem has a unique weak solution. Furthermore, a jump condition capturing finite difference scheme is proposed and applied to solve such elliptic problems. Numerical results show validity and effectiveness of the proposed method.     

Volume preserving immersed boundary methods for two‐phase fluid flows

In this article, we propose a simple area‐preserving correction scheme for two‐phase immiscible incompressible flows with an immersed boundary method (IBM). The IBM was originally developed to model blood flow in the heart and has been widely applied to biofluid dynamics problems with complex geometries and immersed elastic membranes. The main idea of the IBM is to use a regular Eulerian computational grid for the fluid mechanics along with a Lagrangian representation of the immersed boundary. Using the discrete Dirac delta function and the indicator function, we can include the surface tension force, variable viscosity and mass density, and gravitational force effects. The principal advantage of the IBM for two‐phase fluid flows is its inherent accuracy due in part to its ability to use a large number of interfacial marker points on the interface. However, because the interface between two fluids is moved in a discrete manner, this can result in a lack of volume conservation. The idea of an area preserving correction scheme is to correct the interface location normally to the interface so that the area remains constant. Various numerical experiments are presented to illustrate the efficiency and accuracy of the proposed conservative IBM for two‐phase fluid flows. Copyright © 2011 John Wiley & Sons, Ltd.     

A DOMAIN DECOMPOSITION ALGORITHM WITH FINITE ELEMENT-BOUNDARY ELEMENT COUPLING

A domain decomposition algorithm coupling the finite element and the boundary element was presented. It essentially involves subdivision of the analyzed domain into sub-regions being independently modeled by two methods, i.e., the finite element method (FEM) and the boundary element method (BEM). The original problem was restored with continuity and equilibrium conditions being satisfied on the interface of the two sub-regions using an iterative algorithm. To speed up the convergence rate of the iterative algorithm, a dynamically changing relaxation parameter during iteration was introduced. An advantage of the proposed algorithm is that the locations of the nodes on the interface of the two sub-domains can be inconsistent. The validity of the algorithm is demonstrated by the consistence of the results of a numerical example obtained by the proposed method and those by the FEM, the BEM and a present finite element-boundary element (FE-BE) coupling method.     

Solution of the Navier–Stokes equations in velocity–vorticity form using a Eulerian–Lagrangian boundary element method

This paper describes the Eulerian–Lagrangian boundary element model for the solution of incompressible viscous flow problems using velocity–vorticity variables. A Eulerian–Lagrangian boundary element method (ELBEM) is proposed by the combination of the Eulerian–Lagrangian method and the boundary element method (BEM). ELBEM overcomes the limitation of the traditional BEM, which is incapable of dealing with the arbitrary velocity field in advection‐dominated flow problems. The present ELBEM model involves the solution of the vorticity transport equation for vorticity whose solenoidal vorticity components are obtained iteratively by solving velocity Poisson equations involving the velocity and vorticity components. The velocity Poisson equations are solved using a boundary integral scheme and the vorticity transport equation is solved using the ELBEM. Here the results of two‐dimensional Navier–Stokes problems with low–medium Reynolds numbers in a typical cavity flow are presented and compared with a series solution and other numerical models. The ELBEM model has been found to be feasible and satisfactory. Copyright © 2000 John Wiley & Sons, Ltd.     

粘性流体中弹性板振动的有限元耦合问题

流体－结构耦合作用问题是工程中比较常见的问题，具有重要意义，由于流体计算的复杂性，迄今为止，大部分的流体－结构耦合分析都是建立在对流体充分简化的基础上，尤其是将流体视为无粘、无旋的理想流体。该文在近年来前人工作的基础上，发展了一种流体－弹性结构耦合计算模式。流体视为不可压、有粘的介质，流场没有自由表面。该文采用ＳＵ／ＰＧ方法形成流体的有限元方程，采用ＡＬＥ格式处理流体和结构之间的移动界面。采用预估     

An improved Rhie–Chow interpolation scheme for the smoothed‐interface immersed boundary method

Rhie–Chow interpolation is a commonly used method in CFD calculations on a co‐located mesh in order to suppress non‐physical pressure oscillations arising from chequerboard effects. A fully parallelized smoothed‐interface immersed boundary method on a co‐located grid is described in this paper. We discuss the necessity of modifications to the original Rhie–Chow interpolation in order to deal with a locally refined mesh. Numerical simulation with the modified scheme of Choi shows that numerical dissipation due to Rhie–Chow interpolation introduces significant errors at the immersed boundary. To address this issue, we develop an improved Rhie–Chow interpolation scheme that is shown to increase the accuracy in resolving the flow near the immersed boundary. We compare our improved scheme with the modified scheme of Choi by parallel simulations of benchmark flows: (i) flow past a stationary cylinder; (ii) flow past an oscillating cylinder; and (iii) flow past a stationary elliptical cylinder, where Reynolds numbers are tested in the range 10–200. Our improved scheme is significantly more accurate and compares favourably with a staggered grid algorithm. We also develop a scheme to compute the boundary force for the direct‐forcing immersed boundary method efficiently. Copyright © 2016 John Wiley & Sons, Ltd.     

A conservative local interface sharpening scheme for the constrained interpolation profile method

A conservative local interface sharpening scheme has been developed for the constrained interpolation profile method with the conservative semi‐Lagrangian scheme, because the conservative semi‐Lagrangian scheme does not feature a mechanism to control the interface thickness, thus causing an increase of numerical error with the advance of the time step. The proposed sharpening scheme is based on the conservative level set method proposed by Olsson and Kreiss. However, because their method can cause excessive deformation of the free‐surface in certain circumstances, we propose an improvement of the method by developing a local sharpening technique. Several advection tests are presented to assess the correctness of the advection and the improved interface sharpening scheme. This is followed by the validations of dam‐breaking flow and the rising bubble flows. The mass of the fluid is exactly conserved and the computed terminal velocity of the rising bubble agrees well with the experiments compared with other numerical methods such as the volume of fluid method (VOF), the front tracking method, and the level set method. Copyright © 2011 John Wiley & Sons, Ltd.     

模拟相变问题的精细积分边界元法

将精细积分边界元法和界面追踪法相结合求解相变问题。因为边界元法只需要将待求解空间域的边界离散,方便连续追踪移动界面位置和重构网格,所以边界元法适合应用于移动边界问题的模拟。首先,利用精细积分边界元法在固相区域和液相区域分别求解相应的瞬态热传导控制方程,从而求得温度场和边界热流密度。然后,根据固-液相变界面上的能量平衡方程,利用热流密度求得相变界面的移动速度,再采用界面追踪法预测移动相变界面的位置变化。最后,给出了几个数值算例,并通过与参考解的对比验证本文方法的准确性。     

A two‐phase flow interface capturing finite element method

A new interface capturing algorithm is proposed for the finite element simulation of two‐phase flows. It relies on the solution of an advection equation for the interface between the two phases by a streamline upwind Petrov–Galerkin (SUPG) scheme combined with an adaptive mesh refinement procedure and a filtering technique. This method is illustrated in the case of a Rayleigh–Taylor two‐phase flow problem governed by the Stokes equations. Copyright © 2006 John Wiley & Sons, Ltd.     

一种强耦合预估-校正浸入边界法

为克服传统浸入边界法的质量不守恒缺陷,提出了一种用于可压缩流固耦合问题的强耦合预估-校正浸入边界法。通过阐述一般流固耦合系统的矩阵表示,推导了流固耦合系统的强耦合Gauss-Seidel迭代格式,进一步导出预估-校正格式,提出了预估-校正浸入边界法。该方法使用无耦合边界模型对流体进行预估,将流固耦合边界视为自由面,固体原本占据的空间初始化为零质量的单元,允许流体自由穿过耦合边界。对于流体的计算,使用带有minmod限制器的二阶MUSCL有限体积格式和基于Zha-Bilgen分裂的AUSM+-up方法,配合三阶Runge-Kutta格式推进时间步。在校正步骤中,通过一组质量守恒的输运规则来实现输运过程。输运算法可概括为将边界内侧的流体进行标记,根据标记顺序以均匀方式分割和移动流体,产生一个指向边界外侧的流动,最后在边界附近施加速度校正保证无滑移条件。标记和输运算法避免了繁琐的对截断单元的几何处理,确保了算法易于实现。对于固体的计算,分别采用一阶差分格式和隐式动力学有限元格式求解刚体和线弹性体,并利用高斯积分获得固体表面的耦合力。使用预估-校正浸入边界法计算了一维问题和二维问题。在一维活塞问题中,获得了压力分布、相对质量历史和误差曲线,并与其他方法进行了对比。在二维的激波冲击平板问题中,获得了数值模拟纹影和平板结构的挠度历史,并与实验结果进行了对比。研究表明,该方法区别于传统的虚拟网格方法和截断单元方法,能够精确地维持流场的质量守恒并易于实现,且具有一阶收敛精度,能够较准确地预测激波绕射后的流场以及平板在激波作用下的挠度,为开发流固耦合算法提供了一种新的思路。     

A combined vortex and panel method for numerical simulations of viscous flows: a comparative study of a vortex particle method and a finite volume method

This paper describes and compares two vorticity‐based integral approaches for the solution of the incompressible Navier–Stokes equations. Either a Lagrangian vortex particle method or an Eulerian finite volume scheme is implemented to solve the vorticity transport equation with a vorticity boundary condition. The Biot–Savart integral is used to compute the velocity field from a vorticity distribution over a fluid domain. The vorticity boundary condition is improved by the use of an iteration scheme connected with the well‐established panel method. In the early stages of development of flows around an impulsively started circular cylinder, and past an impulsively started foil with varying angles of attack, the computational results obtained by the Lagrangian vortex method are compared with those obtained by the Eulerian finite volume method. The comparison is performed separately for the pressure fields as well. The results obtained by the two methods are in good agreement, and give a better understanding of the vorticity‐based methods. Copyright © 2005 John Wiley & Sons, Ltd.     

无单元Galerkin法和边界元耦合法

无网格方法与有限元法或边界元法耦合是无网格方法处理边界条件的方法之一,在无网格方法中研究无网格方法与有限元法或边界元法耦合的研究显得非常重要.本文在无单元Galerkin法和边界元法的基础上,基于无单元Galerkin法子域和边界元法子域的界面上位移连续和面力平衡条件,提出了一种新的无单元Galerkin法和边界元法的直接耦合方法,对弹性力学问题详细推导了在整个求解域上的耦合公式.与以往的耦合法相比,这种方法简单直观,不需要增加新的耦合区域,也不需要建立新的逼近函数来保证界面位移的连续性.算例结果表明,该方法具有较好的计算精度.     

Three‐dimensional vortex simulation of unsteady flow in hydraulic turbines

On the basis of the Helmholtz decomposition, a grid‐free numerical scheme is provided for the solution of unsteady flow in hydraulic turbines. The Lagrangian vortex method is utilized to evaluate the convection and stretch of the vorticity, and the BEM is used to solve the Neumann problem to define the potential flow. The no‐slip boundary condition is satisfied by generating vortex sticks at the solid surface. A semi‐analytical regularization technique is applied to evaluate the singular boundary surface integrals of the potential velocity and its gradients accurately. The fast multipole method was extended to evaluate the velocity and velocity gradients induced by the discretized vortex blobs in the Lagrangian vortex method. The successful simulation for the unsteady flow through a hydraulic turbine's runner has manifested the effectiveness of the proposed method. Copyright © 2011 John Wiley & Sons, Ltd.