弱可压缩流体边界处理算法

针对流体与固体边界的交互模拟问题,提出一种基于弱可压缩光滑粒子流体动力学(SPH)的边界处理算法.首先,引入一种新的体积权重函数,解决固体边界非均匀采样区域流体密度的计算误差问题;然后,提出一种新的边界力计算模型,避免校正流体粒子位置信息,保证固体边界不可穿透;最后,提出一种改进的流体压力计算模型,保证流体的弱可压缩性.实验结果表明,所提算法可以有效地解决基于位置校正的边界处理方法在模拟弱可压缩流体与非均匀采样固体边界交互时存在的稳定性问题,且仅需边界粒子的位置信息,在节约内存的同时避免了位置校正所带来的额外计算开销.     

Lagrangian particle model for multiphase flows

A Lagrangian particle model for multiphase multicomponent fluid flow, based on smoothed particle hydrodynamics (SPH), was developed and used to simulate the flow of an emulsion consisting of bubbles of a non-wetting liquid surrounded by a wetting liquid. In SPH simulations, fluids are represented by sets of particles that are used as discretization points to solve the Navier-Stokes fluid dynamics equations. In the multiphase multicomponent SPH model, a modified van der Waals equation of state is used to close the system of flow equations. The combination of the momentum conservation equation with the van der Waals equation of state results in a particle equation of motion in which the total force acting on each particle consists of many-body repulsive and viscous forces, two-body (particle-particle) attractive forces, and body forces such as gravitational forces. Similar to molecular dynamics, for a given fluid component the combination of repulsive and attractive forces causes phase separation. The surface tension at liquid-liquid interfaces is imposed through component dependent attractive forces. The wetting behavior of the fluids is controlled by phase dependent attractive interactions between the fluid particles and stationary particles that represent the solid phase. The dynamics of fluids away from the interface is governed by purely hydrodynamic forces. Comparison with analytical solutions for static conditions and relatively simple flows demonstrates the accuracy of the SPH model.     

An Efficient Hybrid Incompressible SPH Solver with Interface Handling for Boundary Conditions

We propose a hybrid smoothed particle hydrodynamics solver for efficientlysimulating incompressible fluids using an interface handling method for boundary conditions in the pressure Poisson equation. We blend particle density computed with one smooth and one spiky kernel to improve the robustness against both fluid–fluid and fluid–solid collisions. To further improve the robustness and efficiency, we present a new interface handling method consisting of two components: free surface handling for Dirichlet boundary conditions and solid boundary handling for Neumann boundary conditions. Our free surface handling appropriately determines particles for Dirichlet boundary conditions using Jacobi‐based pressure prediction while our solid boundary handling introduces a new term to ensure the solvability of the linear system. We demonstrate that our method outperforms the state‐of‐the‐art particle‐based fluid solvers.     

Effects of surface roughness and interface wettability on nanoscale flow in a nanochannel

Non-equilibrium molecular dynamics simulations have been carried out to investigate the effect of surface roughness and interface wettability on the nanorheology and slip boundary condition of simple fluids in a nanochannel of several atomic diameters width. The solid surfaces decorated with periodic nanostrips are considered as the rough surface in this study. The simulation results showed that the interface wettability and the surface roughness are important in determining the nanorheology of the nanochannel and fluid slip at solid–fluid interface. It is observed that the presence of surface roughness always suppresses the fluid slip for hydrophilic and hydrophobic surface nanochannels. For fluids over smooth and hydrophobic surfaces, the snapshots of fluid molecules show that an air gap or nanobubble exists at the fluid–solid interface, resulting in the apparent slip velocity. For a given surface with fixed interface wettability, the fluid velocities increase by increasing the driving force, while the driving force has no significant influence on the density structure of fluid molecules. The fluid slip and the flow rate are measured for hydrophilic and hydrophobic nanochannels. The flow rates in rough surface nanochannels are smaller than those of smooth surface walls due to the increase of drag resistance at the solid–fluid interface. The dependence between fluid slip and flow rate showed that the slip length increases approximately linearly with the flow rate for both the hydrophobic and hydrophilic surface nanochannels.     

固体入水的边界压力处理及表面粒子位置的修正

基于传统光滑粒子流体动力学(SPH)方法的边界力法、虚粒子法或耦合力法处理固体入水时,流体与固体交互界面的粒子密度不连续、压力不稳定、固体边界处会发生部分流体粒子穿透或分离等现象,而流体表面因为受到力的作用,表面破碎后,液面较粗糙.针对上述问题,结合边界力和虚粒子的优点,对耦合力法进行改进,处理运动固体边界,阻止流体粒...     

对称区域边界处理方法及基于表面粒子提取的表面张力计算

流固边界处理一直是流体模拟的研究重点,边界力法和虚粒子法是研究流固边界 的常用方法。边界力法通过对铺设在边界上的粒子施加排斥力防止粒子穿透,但边界力的计算 限制了模拟速度。虚粒子法在边界处生成虚粒子,随着粒子数的增加所需的虚粒子数也随之增 加,导致计算速度下降,且会出现流体与边界分离的现象。为此,提出一种对称区域边界处理 方法,在保证逼真度的前提下满足实时性要求,随着粒子数的增加,其耗时增长也明显比其他 传统方法慢,更适合对复杂场景的模拟,同时避免了边界处流体与边界分离的现象。CSF 方法 是处理表面张力常用的方法,可将表面张力看作体积力进行计算,大大减弱了表面形状对曲率 计算的影响,而事实上曲率的计算只与表面的形状有关。为此,对CSF 方法进行了改进,提出 了一种基于表面粒子提取的表面张力计算方法,减小了传统CSF 方法计算曲率的误差,提高了 计算速度。模拟仿真的效果验证了该方法的有效性。     

Decay rates for the viscous incompressible MHD equations with and without surface tension

In this paper, we consider a layer of a viscous incompressible electrically conducting fluid interacting with the magnetic field in a horizontally periodic setting. The upper boundary is bounded by a free boundary and below bounded by a flat rigid interface. We prove the global well-posedness of the problem for both the case with and without surface tension. Moreover, we show that the global solution decays to the equilibrium exponentially in the case with surface tension, however the global solution decays to the equilibrium at an almost exponential rate in the case without surface tension.     

Molecular dynamics-based prediction of boundary slip of fluids in nanochannels

Computational modeling and simulation can provide an effective predictive capability for flow properties of the confined fluids in micro/nanoscales. In this paper, considering the boundary slip at the fluid–solid interface, the motion property of fluids confined in parallel-plate nanochannels are investigated to couple the atomistic regime to continuum. The corrected second-order slip boundary condition is used to solve the Navier–Stokes equations for confined fluids. Molecular dynamics simulations for Poiseuille flows are performed to study the influences of the strength of the solid–fluid coupling, the fluid temperature, and the density of the solid wall on the velocity slip at the fluid boundary. For weak solid–fluid coupling strength, high temperature of the confined fluid and high density of the solid wall, the large velocity slip at the fluid boundary can be obviously observed. The effectiveness of the corrected second-order slip boundary condition is demonstrated by comparing the velocity profiles of Poiseuille flows from MD simulations with that from continuum.     

A generalized continuous surface tension force formulation for phase-field models for multi-component immiscible fluid flows

We present a new phase-field method for modeling surface tension effects on multi-component immiscible fluid flows. Interfaces between fluids having different properties are represented as transition regions of finite thickness across which the phase-field varies continuously. At each point in the transition region, we define a force density which is proportional to the curvature of the interface times a smoothed Dirac delta function. We consider a vector valued phase-field, the velocity, and pressure fields which are governed by multi-component advective Cahn–Hilliard and modified Navier–Stokes equations. The new formulation makes it possible to model any combination of interfaces without any additional decision criteria. It is general, therefore it can be applied to any number of fluid components. We give computational results for the four component fluid flows to illustrate the properties of the method. The capabilities of the method are computationally demonstrated with phase separations via a spinodal decomposition in a four-component mixture, pressure field distribution for three stationary drops, and the dynamics of two droplets inside another drop embedded in the ambient liquid.     

A new algorithm for simulating flows of conducting fluids in the presence of electric fields

We propose an algorithm based on dissipative particle dynamics (DPD) for simulations of conducting fluids in the presence of an electric field. In this model, the electrostatic equations are solved in each DPD time step to determine the charge density at the fluid surfaces. These surface charges are distributed on a thin layer of fluid particles near the interface, and the corresponding interfacial electric forces are added to other DPD forces. The algorithm is applied to the electrospinning process at the Taylor cone formation stage. It is shown that, when the applied voltage is sufficiently high, the algorithm captures the formation of a Taylor cone with analytical apex angle 98.6°. Our results demonstrate the potential of the presented DPD algorithm for simulating two-phase problems in the presence of an electric field with non-periodic boundary conditions.     

A Boundary Condition Capturing Method for Multiphase Incompressible Flow

In [6], the Ghost Fluid Method (GFM) was developed to capture the boundary conditions at a contact discontinuity in the inviscid compressible Euler equations. In [11], related techniques were used to develop a boundary condition capturing approach for the variable coefficient Poisson equation on domains with an embedded interface. In this paper, these new numerical techniques are extended to treat multiphase incompressible flow including the effects of viscosity, surface tension and gravity. While the most notable finite difference techniques for multiphase incompressible flow involve numerical smearing of the equations near the interface, see, e.g., [19, 17, 1], this new approach treats the interface in a sharp fashion.     

Lattice Boltzmann simulation of viscous fingering phenomenon of immiscible fluids displacement in a channel

In this paper, the viscous fingering phenomenon of two immiscible fluids in a channel is studied by applying the lattice Boltzmann method (LBM). The fundamental physical mechanisms of a finger formation or the interface evolution between immiscible fluids are described in terms of the relative importance of viscous forces, surface tension, and gravity, which are quantifiable via the dimensionless quantities, namely, capillary number, Bond number and viscosity ratio between displaced fluid and displacing fluid. In addition, the effect of wettability on flow behaviour of fluids is investigated for the cases with and without consideration of gravity, respectively. The numerical results provide a good understanding of the mechanisms of viscous fingering phenomenon from a mesoscopic point of view and confirm that the LBM can be viewed as a promising tool for investigating fluid behaviour and other immiscible displacement problems.     

The flow induced in a three layer stratified fluid by a submerged sink or source with stagnation points on the free surfaces

A boundary integral technique is developed to study the free surface flow of a steady, two-dimensional, incompressible, irrotational and inviscid fluid which is induced in both two and three layer stratified fluids in the presence of gravity by a submerged sink or source with stagnation points on the free surfaces. A special form of the Riemann–Hilbert problem, namely the Dirichlet boundary problem, is applied in the derivation of the governing non-linear boundary integral–differential equations which have been solved for the fluid velocity on the free surfaces and this involves the use of an interpolative technique and an iterative process. Results have been obtained for the free surface flow for various Froude numbers and sink heights in both two and three layer fluids. Further, we have also studied the critical Froude numbers for which no convergent solutions are possible for any larger values of the Froude number. We have found that the free surfaces are dependent on two parameters, namely the Froude number and the ratio of sink height to the thickness of either the middle layer in a three layer system and the bottom layer in a two layer system.     

A simple stabilized finite element method for solving two phase compressible–incompressible interface flows

When computing interface flows between compressible (gas) and incompressible (liquid) fluids, one faces at least to the following difficulties: (1) transition from a gas density linked to the local temperature and pressure by an equation of state to a liquid density mainly constant in space, (2) proper approximation of the divergence constraint in incompressible regions and (3) wave transmission at the interface. The aim of the present paper is to design a global (i.e. the same for each phase) numerical method to address easily this coupling. To this end, the same set of primitive unknowns and equations is used everywhere in the flow, but with a dynamic parameterization that changes from compressible to incompressible regions. On one hand, the compressible Navier–Stokes equations are considered under weakly compressibility assumption so that a non-conservative formulation can be used. On the other hand, the incompressible non-isothermal model is retained. In addition, the level set transport equation is used to capture the interface position needed to identify the local characteristics of the fluid and to recover the adequate local modelling. For space approximation of Navier–Stokes equations, a Galerkin least-squares finite element method is used. Two essential elements for defining this numerical scheme are the stabilization and the computation of element integral of the approximated weak form. Since very different concerns motivate the need for stabilization in compressible and incompressible flows, the first difficulty is to design a stabilization operator suitable for both types of flows especially in mixed elements. In addition, some integral of discontinuous functions must be correctly computed to ensure interfacial wave transmission. To overcome these two difficulties, specific averages are computed especially near the interface. Finally, the level set transport equation is computed by a quadrature free Discontinuous Galerkin method. Numerical strategies are performed and validated for 1D and 2D applications.     

A numerical method for two phase flows with insoluble surfactants

In many practical multiphase flow problems, i.e. treatment of gas emboli and various microfluidic applications, the effect of interfacial surfactants, or surface reacting agents, on the surface tension between the fluids is important. The surfactant concentration on an interface separating the fluids can be modeled with a time dependent differential equation defined on the moving and deforming interface. The equations for the location of the interface and the surfactant concentration on the interface are coupled with the Navier–Stokes equations. These equations include the singular surface tension forces from the interface on the fluid, which depend on the interfacial surfactant concentration.A new accurate and inexpensive numerical method for simulating the evolution of insoluble surfactants is presented in this paper. It is based on an explicit yet Eulerian discretization of the interface, which for two dimensional flows allows for the use of uniform one dimensional grids to discretize the equation for the interfacial surfactant concentration. A finite difference method is used to solve the Navier–Stokes equations on a regular grid with the forces from the interface spread to this grid using a regularized delta function. The timestepping is based on a Strang splitting approach.Drop deformation in shear flows in two dimensions is considered. Specifically, the effect of surfactant concentration on the deformation of the drops is studied for different sets of flow parameters.     

A coupling algorithm for partitioned solvers applied to bubble and droplet dynamics

Bubbles and droplets both consist of a liquid in contact with a gas. In this paper, we consider the interface between the incompressible liquid and the gas as a zero thickness structure. The position of the interface is determined by the equilibrium between surface tension effects and the fluid pressure difference across the interface. So, the structure interacts with the fluids on either side. The behaviour of a limited number of bubbles and droplets can therefore be simulated as a Fluid-Structure Interaction (FSI) problem.Most existing techniques frequently used for studying bubble and droplet dynamics, such as Level Set or Volume Of Fluid, use monolithic schemes. The flow on both sides of the interface and the position of the interface are calculated in a single code. In this contribution, a partitioned approach is presented. The position of the interface is calculated with a structural solver. Given a displacement of the interface, a separate flow solver calculates the flow on the liquid side of the interface with the Arbitrary Lagrangian-Eulerian (ALE) technique. The structural solver uses a reduced order model of the flow solver to obtain implicit coupling between both solvers. This reduced order model is built up during the coupling iterations of a time step. Grid and time converged solutions of two axisymmetric problems are calculated: an oscillating water droplet in air and the growth and detachment of an air bubble from the outlet of a vertical needle, submerged in quiescent water.     

Numerical simulation of a compressible turbulent boundary layer over a conductive wall with line heat source

A conjugated problem of supersonic turbulent flow over a conductive solid wall with an embedded line heat source has been investigated as a model of a separation detector and skin friction gage. The 2-D Navier-Stokes equations for compressible fluid, including a two layer eddy viscosity model, are solved simultaneously with the heat transfer equation for the solid, written in general coordinates. The effect of the interface boundary condition on the stability of the implicit scheme of the flow field has been checked. A careful investigation of the effect of heat source strength, solid and fluid conductivity and Mach and Reynolds numbers on flow and temperature fields has been performed. The results of this investigation may be used to design an optimal gage with a minimum influence on the flow field.     

Numerical solution of a phase field model for incompressible two-phase flows based on artificial compressibility

We present an artificial compressibility based numerical method for a phase field model for simulating two-phase incompressible viscous flows. The phase model was proposed by Liu and Shen [Physica D. 179 (2003) 211–228], in which the interface between two fluids is represented by a thin transition region of fluid mixture that stores certain amount of mixing energy. The model consists of the Navier–Stokes equations coupled with the Allen–Cahn equation (phase field equation) through an extra stress term and a transport term. The extra stress in the momentum equations represents the phase-induced capillary effect for the mixture due to the surface tension. The coupled equations are cast into a conservative form suitable for implementation with the artificial compressibility method. The resulting hyperbolic system of equations are then discretized with weighted essentially non-oscillatory (WENO) finite difference scheme. The dual-time stepping technique is applied for obtaining time accuracy at each physical time step, and the approximate factorization algorithm is used to solve the discretized equations. The effectiveness of the numerical method is demonstrated in several two-phase flow problems with topological changes. Numerical results show the present method can be used to simulate incompressible two-phase flows with small interfacial width parameters and topological changes.     

基于CUDA的弱可压SPH流体建模与仿真

为了实现小尺度范围流体场景的实时、真实感模拟,采用弱可压SPH方法对水体进行建模,提出了流体计算的CPU GPU混合架构计算方法。针对邻域粒子查找算法影响流体计算效率的问题,采用三维空间网格对整个模拟区域进行均匀网格划分,利用并行前缀求和和并行计数排序实现邻域粒子的查找。最后,采用基于CUDA并行加速的Marching Cubes算法实现流体表面提取,利用环境贴图表现流体的反射和折射效果,实现流体表面着色。实验结果表明,所提出的流体建模和模拟算法能实现小尺度范围流体的实时计算和渲染,绘制出水的波动、翻卷和木块在水中晃动的动态效果,当粒子数达到1 048 576个时,GPU并行计算方法相较CPU方法的加速比为60.7。