A discrete forcing method dedicated to moving bodies in two‐phase flow

The numerical simulation of interaction between structures and two‐phase flows is a major concern for many industrial applications. To address this challenge, the motion of structures has to be tracked accurately. In this work, a discrete forcing method based on a porous medium approach is proposed to follow a nondeformable rigid body with an imposed velocity by using a finite‐volume Navier‐Stokes solver code dedicated to multiphase flows and based on a two‐fluid approach. To deal with the action reaction principle at the solid wall interfaces in a conservative way, a porosity is introduced allowing to locate the solid and insuring no diffusion of the fluid‐structure interface. The volumetric fraction equilibrium is adapted to this novelty. Mass and momentum balance equations are formulated on a fixed Cartesian grid. Interface tracking is addressed in detail going from the definition of the porosity to the changes in the discretization of the momentum balance equation. This so‐called time‐ and space‐dependent porosity method is then validated by using analytical and elementary test cases.     

Numerical simulation of a single bubble by compressible two‐phase fluids

The present work deals with the numerical investigation of a collapsing bubble in a liquid–gas fluid, which is modeled as a single compressible medium. The medium is characterized by the stiffened gas law using different material parameters for the two phases. For the discretization of the stiffened gas model, the approach of Saurel and Abgrall is employed where the flow equations, here the Euler equations, for the conserved quantities are approximated by a finite volume scheme, and an upwind discretization is used for the non‐conservative transport equations of the pressure law coefficients. The original first‐order discretization is extended to higher order applying second‐order ENO reconstruction to the primitive variables. The derivation of the non‐conservative upwind discretization for the phase indicator, here the gas fraction, is presented for arbitrary unstructured grids. The efficiency of the numerical scheme is significantly improved by employing local grid adaptation. For this purpose, multiscale‐based grid adaptation is used in combination with a multilevel time stepping strategy to avoid small time steps for coarse cells. The resulting numerical scheme is then applied to the numerical investigation of the 2‐D axisymmetric collapse of a gas bubble in a free flow field and near to a rigid wall. The numerical investigation predicts physical features such as bubble collapse, bubble splitting and the formation of a liquid jet that can be observed in experiments with laser‐induced cavitation bubbles. Opposite to the experiments, the computations reveal insight to the state inside the bubble clearly indicating that these features are caused by the acceleration of the gas due to shock wave focusing and reflection as well as wave interaction processes. While incompressible models have been used to provide useful predictions on the change of the bubble shape of a collapsing bubble near a solid boundary, we wish to study the effects of shock wave emissions into the ambient liquid on the bubble collapse, a phenomenon that may not be captured using an incompressible fluid model. Copyright © 2009 John Wiley & Sons, Ltd.     

Flow simulation on moving boundary‐fitted grids and application to fluid–structure interaction problems

We present a method for the parallel numerical simulation of transient three‐dimensional fluid–structure interaction problems. Here, we consider the interaction of incompressible flow in the fluid domain and linear elastic deformation in the solid domain. The coupled problem is tackled by an approach based on the classical alternating Schwarz method with non‐overlapping subdomains, the subproblems are solved alternatingly and the coupling conditions are realized via the exchange of boundary conditions. The elasticity problem is solved by a standard linear finite element method. A main issue is that the flow solver has to be able to handle time‐dependent domains. To this end, we present a technique to solve the incompressible Navier–Stokes equation in three‐dimensional domains with moving boundaries. This numerical method is a generalization of a finite volume discretization using curvilinear coordinates to time‐dependent coordinate transformations. It corresponds to a discretization of the arbitrary Lagrangian–Eulerian formulation of the Navier–Stokes equations. Here the grid velocity is treated in such a way that the so‐called Geometric Conservation Law is implicitly satisfied. Altogether, our approach results in a scheme which is an extension of the well‐known MAC‐method to a staggered mesh in moving boundary‐fitted coordinates which uses grid‐dependent velocity components as the primary variables. To validate our method, we present some numerical results which show that second‐order convergence in space is obtained on moving grids. Finally, we give the results of a fully coupled fluid–structure interaction problem. It turns out that already a simple explicit coupling with one iteration of the Schwarz method, i.e. one solution of the fluid problem and one solution of the elasticity problem per time step, yields a convergent, simple, yet efficient overall method for fluid–structure interaction problems. Copyright © 2005 John Wiley & Sons, Ltd.     

动脉分岔血管内膜增生过程的数值模拟

内膜增生从发生到阻塞血管是一个复杂的变化过程,在这个过程中,内膜的增生、血管腔体形状的改变和血流动力学之间是相互影响的。为了研究这些变化,本文提出一种单元填充方法数值模拟了三维颈动脉分岔血管在低切应力作用下血管内膜增生的过程。该方法既可以克服节点移动方法所不可避免的内膜增生的不连续性,也可以避免网格重划分的困难。结果发现,如果单纯以切应力阈值作为内膜增生的判据,低切应力的作用将无法导致血管完全阻塞,但内膜增生和血流动力学之间的相互影响是可以通过数值方法进行模拟的。在本数值模拟中,内膜增生的过程分为"增厚"(先)和"扩展"(后)两个阶段,最大狭窄率为34.4%,发生在距血管分岔5mm处动脉窦的外侧壁面。其发生位置和形状与临床观察吻合。     

A NUMERICAL METHOD FOR SIMULATING NONLINEAR FLUID-RIGID STRUCTURE INTERACTION PROBLEMS

A numerical method for simulating nonlinear fluid-rigid structure interaction problems is developed. The structure is assumed to undergo large rigid body motions and the fluid flow is governed by nonlinear, viscous or non-viscous, field equations with nonlinear boundary conditions applied to the free surface and fluid-solid interaction interfaces. An Arbitrary-Lagrangian-Eulerian (ALE) mesh system is used to construct the numerical model. A multi-block numerical scheme of study is adopted allowing for the relative motion between moving overset grids, which are independent of one another. This provides a convenient method to overcome the difficulties in matching fluid meshes with large solid motions. Nonlinear numerical equations describing nonlinear fluid-solid interaction dynamics are derived through a numerical discretization scheme of study. A coupling iteration process is used to solve these numerical equations. Numerical examples are presented to demonstrate applications of the model developed.     

New strategy of solid/fluid coupling during numerical simulation of thermo-mechanical processes

In this study, numerical methods are developed to simulate thermomechanical processes, taking into account both the fluid flows in the molten pool and the deformations of the solid parts. The methods are based on a new strategy of solid/fluid coupling. They allow to simulate the formation of the molten pool by taking into account the fluid flows through both effects of the surface tension (“curvature effect” and “Marangoni effect”) and the buoyancy. An ALE approach is used to follow the evolution of the free surface. The effects of the deformations in the base metal on the fluid flows in the molten pool (solid/fluid interaction) is ensured by imposing the velocities of the solid nodes during the thermo-fluid simulation. As an application, a thermo-fluid-mechanical simulation of laser welding is carried out. It is found that the solid/fluid interaction has a minor effect on simulation results.     

基于虚拟区域法的黏弹性流体中微生物游动的数值模拟和应用

微生物是自然生态系统的重要组成部分,掌握微生物在复杂流体中的运动特性可以为微型器件的设计制造提供理论指导.壁面效应是微生物游动研究中的重要问题之一,已有研究表明微生物在壁面附近存在复杂的行为特征.然而已有研究大多集中于微生物在牛顿流体中的游动模拟,仅有少数涉及黏弹性流体等非牛顿流体.本文采用直接力虚拟区域法与乔列斯基分解相结合的数值方法,引入Squirmer微生物游动模型,研究了微生物在黏弹性流体中的游动问题.首先给出求解黏弹性流体本构方程的数值格式;并将该方法应用于研究微生物游动中的壁面效应.研究结果表明,游动方向是影响微生物颗粒壁面效应的重要因素.流体弹性应力会对微生物产生一个反向转矩,影响微生物的游动方向,从而阻碍微生物逃离壁面.微生物颗粒在黏弹性流体中与壁面作用时间较长,几乎达到牛顿流体的两倍以上.     

Investigation of dense slurry suspensions with coaxial mixers: Influences of design variables through tomography and mathematical modelling

The coaxial mixers enhance the suspension of concentrated slurries in an agitated reactor. In this research work, the complex slurry suspension and dissemination behavior in a coaxial slurry mixing system (comprised of a close clearance anchor rotating with a low speed and an inner axial impeller rotating with a high speed) was analyzed employing ERT (electrical resistance tomography, a non-intrusive flow visualization technique), and computational fluid dynamics (CFD). The numerical models were validated by comparing the axial solid concentration profiles generated using the ERT data and the CFD simulation results. The influences of various important parameters such as the diameter of the inner axial impeller, the inner impeller type, and the inner impeller spacing on the hydrodynamic characteristics of the slurry suspensions in a coaxial mixing vessel were thoroughly analyzed. The radial and axial velocity profiles of solid particles were generated using the validated mathematical models. The assessment of energy loss due to the solid–solid collisions, the particle–fluid frictions, and the particle–vessel wall collisions was conducted. The evaluation of optimum inner impeller clearance and inner impeller diameter is essential to attain a high degree of solids suspension and dissemination in a coaxial slurry mixing system.     

Mathematical modelling of two‐phase non‐Newtonian flow in a helical pipe

Governing equations for a two‐phase 3D helical pipe flow of a non‐Newtonian fluid with large particles are derived in an orthogonal helical coordinate system. The Lagrangian approach is utilized to model solid particle trajectories. The interaction between solid particles and the fluid that carries them is accounted for by a source term in the momentum equation for the fluid. The force‐coupling method (FCM), developed by M.R. Maxey and his group, is adopted; in this method the momentum source term is no longer a Dirac delta function but is spread on a numerical mesh by using a finite‐sized envelop with a spherical Gaussian distribution. The influence of inter‐particle and particle–wall collisions is also taken into account. Copyright © 2005 John Wiley & Sons, Ltd.     

基于2D网格的轴对称浸没边界法

浸没边界法是处理颗粒两相流中运动边界问题的一种常用数值模拟方法. 当研究的物理问题的无量纲参数满足一定要求时, 该流场结构呈现轴对称状态. 为此本文提出了一种基于2D笛卡尔网格和柱坐标系的轴对称浸没边界法. 该算法采用有限体积法(FVM)对动量方程进行空间离散, 并通过阶梯状锐利界面替代真实的固体浸没边界来封闭控制方程. 为了提高计算效率, 本文采用自适应网格加密技术提高浸没边界附近网格分辨率. 由于柱坐标系的使用, 使得动量方程中的黏性项产生多余的源项, 我们对其作隐式处理. 此外, 在对小球匀速近壁运动进行直接数值模拟时, 由于球壁间隙很小, 间隙内的压力变化比较剧烈. 因此想要精确地解析流场需要很高的网格分辨率. 此时, 需要在一个时间步内多次实施投影步来保证计算的稳定性. 而在小球自由碰壁运动中, 我们通过引入一个润滑力模型使得低网格分辨率下也能模拟小球近壁处的运动. 最后通过小球和圆盘绕流、Stokes流小球近壁运动以及小球自由下落碰壁弹跳算例验证本算法对于轴对称流的静边界和动边界问题均是适用和准确的.      

基于孔隙弹性耦合的水力压裂数值模拟研究

系统地阐述基于多孔介质渗流-损伤耦合原理,进行水力压裂的FEM数值实现方法.基本架构为:(1)引入孔隙流体压力膨胀系数将孔隙流体压力与应力场进行耦合;(2)基于损伤局部化模型,提出裂缝张开度表达式;(3)提出水力压裂引起的多孔介质水-力学属性的各向异性表达式;(4)提出全流量加载的耦合分析方案.作为实例,模拟三维地层水压裂缝扩展形态,通过比较模型的数值解和经典理论解,验证该方法的正确性.     

求解三维流固耦合问题的一种全隐全耦合区域分解并行算法

三维流固耦合问题的非结构网格数值算法在很多工程领域都有重要应用,目前现有的数值方法主要基于分区算法,即流体和固体区域分别进行求解,因此存在收敛速度较慢以及附加质量导致的稳定性问题,此外,该类算法的并行可扩展性不高,在大规模应用计算方面也受到一定限制．本文针对三维非定常流固耦合问题,提出一种基于区域分解的全隐全耦合可扩展并行算法．首先基于任意拉格朗日-欧拉框架建立流固耦合控制方程,然后时间方向采用二阶向后差分隐式格式、空间方向采用非结构稳定化有限元方法进行离散．对于大规模非线性离散系统,构造一种结合非精确Newton法、Krylov子空间迭代法与区域分解Schwarz预条件子的Newton-Krylov-Schwarz (NKS)并行求解算法,实现流体、固体和动网格方程的一次性整体求解．采用弹性障碍物绕流的标准测试算例对数值方法的准确性进行了验证,数值性能测试结果显示本文构造的全隐全耦合算法具有良好的稳定性,在不同的物理参数下具有良好的鲁棒性,在“天河二号”超级计算机上,当并行规模从192增加到3072个处理器核时获得了91%的并行效率．性能测试结果表明本文构造的NKS算法有望应用于复杂...     

Effects of grid size on predictions of bed expansion in bubbling fluidized beds of Geldart B particles: A generalized rule for a grid-independent solution of TFM simulations

Numerical simulations of gas–solid fluidized beds based on the kinetic theory of granular flow exhibit a significant dependence on domain discretization. Bubble formation, bubble size and shape all vary greatly with the discretization, and the use of an inappropriate scale resolution leads to inaccurate predictions of fluidization hydrodynamics. In this study, grid-independent solutions of the two fluid model were examined by comparing the bed expansions obtained from numerical simulations with experimental results and empirical predictions, based on bubbling fluidized beds of Geldart B particles. Grid independence was achieved with a grid resolution equal to 18 times the particle diameter. The simulation results were compared with previously published data for verification purposes. The results of this work should provide a guideline for choosing the appropriate grid size and thereby minimize the time and expense associated with large simulations.     

Numerical modeling and investigation of viscoelastic fluid–structure interaction applying an implicit partitioned coupling algorithm

We investigate the interaction between a viscoelastic Oldroyd-B fluid and an elastic structure via simulations applying an implicit partitioned coupling algorithm. Simulations are done for a flow through a channel with a flexible wall and a lid-driven cavity flow with flexible bottom. In addition, we make use of a mass–spring–dashpot prototype model to study the dynamic interaction problem. Both the simulation results and the analysis of the prototype model show that there are obvious differences in the fluid–structure interaction if the fluids are viscoelastic instead of purely viscous. These differences appear in the deformation of the solid at stationary state and in the equilibrium position, amplitude, frequency as well as phase shift of the oscillation. Moreover, we investigate the influence of numerical and physical parameters on the implicit partitioned coupling algorithm for simulation of viscoelastic fluid–structure interaction problems.     

Moving immersed boundary method

A novel implicit immersed boundary method of high accuracy and efficiency is presented for the simulation of incompressible viscous flow over complex stationary or moving solid boundaries. A boundary force is often introduced in many immersed boundary methods to mimic the presence of solid boundary, such that the overall simulation can be performed on a simple Cartesian grid. The current method inherits this idea and considers the boundary force as a Lagrange multiplier to enforce the no‐slip constraint at the solid boundary, instead of applying constitutional relations for rigid bodies. Hence excessive constraint on the time step is circumvented, and the time step only depends on the discretization of fluid Navier‐Stokes equations, like the CFL condition in present work. To determine the boundary force, an additional moving force equation is derived. The dimension of this derived system is proportional to the number of Lagrangian points describing the solid boundaries, which makes the method very suitable for moving boundary problems since the time for matrix update and system solving is not significant. The force coefficient matrix is made symmetric and positive definite so that the conjugate gradient method can solve the system quickly. The proposed immersed boundary method is incorporated into the fluid solver with a second‐order accurate projection method as a plug‐in. The overall scheme is handled under an efficient fractional step framework, namely, prediction, forcing, and projection. Various simulations are performed to validate current method, and the results compare well with previous experimental and numerical studies.     

Pulsatile blood flow in a resistance vessel

The quasi-one-dimensional Newtonian fluid flow in an active vessel in which the inlet pressure oscillates about a mean value is considered. The vascular wall is assumed to display the Bayliss effect and a reaction to changes in blood flow. The first is characterized by the dependence of the contractility of the vascular smooth-muscle cells on the intravascular pressure and the second by the sensitivity of the endothelial cells to the wall shear stress. The dependence of the period-average radius and blood flow rate on the mean value and pulsation amplitude and frequency of the inlet pressure is investigated numerically. The characteristics are compared for the passive vessel, the vessel with the Bayliss effect, and the vessel with both reactions.     

改进虚拟边界算法在超声速流动问题求解中的应用

提出了一种改进的虚拟单元浸没边界法, 并与一种高阶格式的有限差分算法相结合, 运用于求解超声速复杂几何绕流问题.该算法的核心思想在于在固体边界的内部和外部分别施加满足边界关系的作用点, 使得几何边界离散更加细化, 起到了壁面附近网格局部加密的作用.采用源空间内流体点作为反距离插值算法的重构点, 有效避免了插值点数目过少而与作用点相重合情况.通过对二维激波反射现象 (马赫数为 2.81) 和三维超声速球体绕流问题 (马赫数为 1.2) 的数值模拟, 与实验结果对比表明, 本文改进算法相对一般的虚拟边界法来说能显著提高数值精度, 减小计算误差.计算结果揭示了球体绕流中剪切层、压缩波系和尾迹的相互作用导致自由剪切层失稳的机理.剪切层厚度和湍流雷诺脉动经历了线性增长、大幅度震荡和小幅度波动三个阶段, 导致剪切层表面褶皱因子变化呈指数规律增长.其湍流结构表现出明显的各向异性, 具体在流向雷诺正应力在湍流脉动中占主导地位, 激波的压缩作用对不同方向雷诺正应力的影响存在空间迟滞效应.