MPS-FEM数值分析带自由面的流固耦合问题

随着计算科学的发展,研究人员为探索流固耦合问题的物理机理而提出了众多的数值方法。其中,耦合的移动粒子半隐式方法 MPS(Moving Particle Semi-Implicit method)和有限单元法FEM(Finite Element method)为流固耦合问题的数值仿真工作提供了新的途径。本文所有流场的数值模拟工作均采用课题组自主开发的无网格法求解器MLParticle-SJTU来完成。该求解器在原始的MPS法基础上,对核函数、压力梯度模型、压力泊松方程的求解和自由面判断方式等方面进行了改进。此外,在该求解器框架内,基于FEM法拓展了针对结构场进行求解的功能。首先,对MPS和FEM方法的理论模型及其耦合策略进行了介绍。然后,采用该自研MPS-FEM耦合求解器,数值模拟了溃坝流动对弹性结构的冲击及其相互作用的标准问题。通过将结构变形及自由面波型变化等结果与已发表结果进行对比,验证了该求解器在处理带自由面剧烈变化的粘性流体和柔性变形结构的耦合作用问题上的可行性。     

一种改进的MPS粒子作用模型

从分析角度出发,为移动粒子半隐式法(MPS)提供了一组新的粒子作用模型,该模型在粒子局部坐标系下进行泰勒展开,同步获取各算子的粒子近似式,确保了离散精度的一致性.方法对于随机分布的粒子和边界附近的粒子同样适用.文中利用三个算例对该方法进行了检验,计算结果与分析解符合较好.最后,运用该方法成功模拟了自由面流动的二维溃坝问题.     

MPS方法数值模拟楔形体入水问题

入水问题是船舶海洋工程中典型的流动问题。当船舶在恶劣海况中航行或海洋平台遭遇恶劣天气时,结构物和水体之间往往会出现剧烈的砰击作用。砰击发生时,伴随着结构物湿表面的变化、自由液面的翻卷和破碎等强非线性现象。本文采用本课题组自主开发的基于移动粒子半隐式法MPS(Moving Particle Semi-Implicit Method)的求解器MLParticle-SJTU对二维楔形体入水问题进行了数值模拟。选取斜升角为30°的楔形体模拟其入水过程,研究了不同粒子布置方式对于计算结果如垂向水动力和自由面变化的影响,并与MLM砰击模型(Modified Logvinovich Model)的结果进行了比较,吻合较好。在此基础上,进一步研究了不同斜升角对计算结果(垂向水动力和自由面变化)的影响,其中流动分离前的垂向水动力与MLM结果相近,表明了MPS方法能有效地模拟入水问题。     

具有表面张力的粘性射流的数值计算

本文将流线迭代法应用于求解带有自由表面的粘性不定边界流动问题。並利用这种方法对具有表面张力的粘性轴对称层流射流问题进行了数值计算。计算结果与实验结果基本上符合。     

具有壁面射流激励的圆管内气相流动大涡模拟

采用大涡模拟方法和单方程亚格子模式对小尺度量进行模拟。研究了不同强度壁面射流激励对圆管内气相流动的影响,模拟结果给出了射流对瞬态拟序结构发展、时平均流向速度分布的影响。随着射流强度的增加,射流入口附近流体的回流现象增强。射流强度足够大时可以减小管壁处的切应力值,同时会减小壁面附近流动速度,这种速度分布会导致气体夹带颗粒的能力下降,从而在实际两相流动中容易造成壁面附近的气粒返混现象。     

孔板空调风口送风射流的数值模拟

介绍Ｎ点风口模型用于数值模拟室内空气流动时描述孔板类送风口的入流边界条件．然后采用该风口模型对不同的孔板风口出流条件算例进行数值计算，并就轴心速度衰减、射流扩展角以及断面流速分布等射流特性与实验数据进行了对比．比较结果表明，Ｎ点风口模型用于描述数值模拟室内空气流动的孔板类风口入流边界条件，可以获得工程上足够满意的结果．     

基于格子Boltzmann方法的挤出胀大的数值模拟

将单相格子Boltzmann方法(lattice Boltzmann method, LBM)引入到粘弹流体的瞬态挤出胀大的数值模拟中,建立了基于双分布函数的自由面粘弹性流动格子Boltzmann模型．分析得到的流道中流动速度分布和构型张量结果与理论解十分吻合．对粘弹流体瞬态挤出胀大过程进行了模拟,并分析了运动粘度比和剪切速率对挤出胀大率的影响,得到的胀大率结果与理论分析和其它模拟结果基本一致．表明给出的LBM可以捕捉挤出胀大的瞬态效应．     

气流作用下同轴带电射流的不稳定性研究

通过对气体驱动同轴电流动聚焦的实验模型进行简化,开展了电场力和惯性力共同作用下同轴带电射流的不稳定性理论研究.在流动为无黏、不可压缩、无旋的假设下,建立了三层流体带电射流物理模型并得到了扰动在时间域内发展演化的解析形式色散关系,利用正则模方法求解色散方程发现了流动的不稳定模态,进而分析了主要控制参数对不稳定模态的影响.结果表明,在参考状态下轴对称模态的最不稳定增长率最大,因此轴对称扰动控制整个流场.外层气流速度越高,气体惯性力越大,射流的界面越容易失稳.内外层液-液同轴射流之间的速度差越大,射流越不稳定.表面张力对射流不稳定性起到促进作用.轴向电场对射流不稳定性具有双重影响:当加载电场强度较小时,射流不稳定性被抑制;当施加电压大于某一临界值时,轴向电场会促进射流失稳.临界电压的大小与界面上自由电荷密度和射流表面扰动发展关系密切.这些结果与已有的实验现象吻合,能够对实验的过程控制提供理论指导.     

风口区空气流动的有限元模拟与实验研究

通风和空调出风口的射流特性主要取决于出风口截面之前的一段管道内的流动状况,研究送风口的射流偏转特性有重要实用价值.应用湍流大涡模拟技术结合Taylor-Galerkin有限单元法对工程中常见工况下的风口区管道内的三维流场进行了数值模拟和实验研究,数值计算结果与实验结果符合较好.表明湍流大涡模拟方法适合于边界形状复杂,存在各向异性的大尺度涡的内流情况,能可靠地预报风口区空气流动的射流偏转特性.     

二元工质Marangoni对流的实验研究与数值模拟

用实验和数值模拟方法研究了正值表面张力温度系数的二元醇类特殊水溶液Marangoni对流流动. 首先通过实验测量确定正戊醇水溶液表面张力在特定的浓度分布和温度区间内具有明显随温度升高而增加的变化区域, 然后利用特定浓度配比的正戊醇水溶液,采用PIV方法实验观测了矩形液池中二元工质液层在水平温差驱动下的Marangoni对流, 发现了不同于常规的反向热毛细对流流动, 测量的表面速度分布与相同工况的数值模拟结果进行了比较, 发现二者变化趋势一致. 实验观测和理论结果的比较进一步验证了表面张力温度系数为正值时二元工质液层的热毛细流动输运特性.      

基于表面力模型的液滴冲击弹性基底SPH模拟

采用光滑粒子动力学SPH方法建立液滴冲击弹性基底的流固耦合数值模型,给出描述粘性流体和弹性固体运动的SPH离散方程和数值处理格式,引入人工耗散项来抑制标准SPH方法的数值震荡。为模拟液滴的表面张力效应,通过精确检测边界粒子,采用拉格朗日插值方法计算表面法向量和曲率,结合界面理论中的连续表面力CSF方法,建立了适用于自由表面液滴的表面力模型,方形液滴变形的模拟结果与拉普拉斯理论解吻合较好。随后,采用SPH流固耦合模型模拟1.0 mm直径水滴以不同速度(0.2 m/s～3.0 m/s)冲击两种薄板型基底,分析了基底弹性变形对液滴铺展、收缩以及回弹行为的影响。     

Development of a PLIC-VOF method for the dynamic simulation of entry region flow in a laminar falling film

The present work describes a numerical procedure to simulate the development of hydrodynamic entry region in a gravity-driven laminar liquid film flow over an inclined plane. It provides a better insight into the physics of developing film in entry region. A novel numerical approach is proposed which has the potential to provide solutions for the complex physics of liquid film spreading on solid walls. The method employs an incompressible flow algorithm to solve the governing equations, a PLIC-VOF method to capture the free surface evolution and a continuum surface force (CSF) model to include the effect of surface tension. To account for the moving contact line on the solid substrate, a precursor film model based wall treatment is implemented. Liquid film flow has been simulated for the Reynolds number range of 5 ≤ Re ≤ 37.5, and the predicted results are found to agree well with the available analytical and experimental data.     

Improvement of moving particle semi‐implicit method for simulation of progressive water waves

Precise simulation of the propagation of surface water waves, especially when involving breaking wave, takes a significant place in computational fluid dynamics. Because of the strong nonlinear properties, the treatment of large surface deformation of free surface flow has always been a challenging work in the development of numerical models. In this paper, the moving particle semi‐implicit (MPS) method, an entirely Lagrangian method, is modified to simulate wave motion in a 2‐D numerical wave flume preferably. In terms of consecutive pressure distribution, a new and simple free surface detection criterion is proposed to enhance the free surface recognition in the MPS method. In addition, a revised gradient model is deduced to diminish the effect of nonuniform particle distribution and then to reduce the numerical wave attenuation occurring in the original MPS model. The applicability and stability of the improved MPS method are firstly demonstrated by the calculation of hydrostatic problem. It is revealed that these modifications are effective to suppress the pressure oscillation, weaken the local particle clustering, and boost the stability of numerical algorithm. It is then applied to investigate the propagation of progressive waves on a flat bed and the wave breaking on a mild slope. Comparisons with the analytical solutions and experimental results indicate that the improved MPS model can give better results about the profiles and heights of surface waves in contrast with the previous MPS models. Copyright © 2017 John Wiley & Sons, Ltd.     

Surface tension and buoyancy-driven flow in a non-isothermal liquid bridge

The Navier–Stokes–Boussinesq equations governing the transport of momentum, mass and heat in a non-isothermal liquid bridge with a temperature-dependent surface tension are solved using a vorticity-stream-function formulation together with a non-orthogonal co-ordinate transformation. The equations are discretized using a pseudo-unsteady semi-implicit finite difference scheme and are solved by the ADI method. A Picard-type iteration is adopted which consists of inner and outer iterative processes. The outer iteration is used to update the shape of the free surface. Two schemes have been used for the outer iteration; both use the force balance normal to the free surface as the distinguished boundary condition. The first scheme involves successive approximation by the direct solution of the distinguished boundary condition. The second scheme uses the artificial force imbalance between the fluid pressure, viscous and capillary forces at the free surface which arises when the boundary condition for force balance normal to the surface is not satisfied. This artificial imbalance is then used to change the surface shape until the distinguished boundary condition is satisfied. These schemes have been used to examine a variety of model liquid bridge situations including purely thermocapillary-driven flow situations and mixed thermocapillary- and bouyancy-driven flow.     

贮腔类三维自由液面动力学问题数值研究

讨论了贮腔类三维自由液面动力学问题的数值研究,将任意的拉格朗日-欧拉运动学描述关系引入到系统的控制方程中,采用任意的拉格朗日-欧拉描述跟踪自由液面,推导了自由面上结点的法向矢量计算公式。采用Galerkin余量法推导了Navier-Stokes方程的空间离散有限元方程,采用三维自由液面上微分几何理论推导了表面张力计算公式。数值研究中考虑了接触角效应,最后进行了三维数值算例分析。     

Towards the numerical simulation of the horizontal slug front

Steps towards the numerical simulation of the flow behind the slug front in horizontal slug flow performed with a streamfunction-vorticity representation of the mean flow and an energy dissipation model for the turbulence are discussed. The flow field consists of two vortices, one saddle point and four stagnation regions. Attention is focused on the following boundary conditions: moving wall jet, moving wall, free jet velocity discontinuity and vertical liquid-gas open surface. A dissipation flux boundary condition is suggested to simulate the interaction of the turbulent eddies with the open surface. A method to assess the necessity to use a transport model equation for the dissipation rather than a geometric specification of a length is suggested. Three different ways to characterize the mixing zone length are proposed.     

An improved three-dimensional model for interface pressure calculations in free-surface flows

A three-dimensional method for the calculation of interface pressure in the computational modeling of free surfaces and interfaces is developed. The methodology is based on the calculation of the pressure force at the interfacial cell faces and is mainly designed for volume of fluid (VOF) interface capturing approach. The pressure forces at the interfacial cell faces are calculated according to the pressure imposed by each fluid on the portion of the cell face that is occupied by that fluid. Special formulations for the pressure in the interfacial cells are derived for different orientations of an interface. The present method, referred to as pressure calculation based on the interface location (PCIL), is applied to both static and dynamic cases. First, a three-dimensional motionless drop of liquid in an initially stagnant fluid with no gravity force is simulated as the static case and then two different small air bubbles in water are simulated as dynamic cases. A two-fluid, piecewise linear interface calculation VOF method is used for numerical simulation of the interfacial flow. For the static case, both the continuum surface force (CSF) and the continuum surface stress (CSS) methods are used for surface tension calculations. A wide range of Ohnesorge numbers and density and viscosity ratios of the two fluids are tested. It is shown that the presence of spurious currents (artificial velocities present in case of considerable capillary forces) is mainly due to the inaccurate calculation of pressure forces in the interfacial computational cells. The PCIL model reduces the spurious currents up to more than two orders of magnitude for the cases tested.

Also for the dynamic bubble rise case, it is shown that using the numerical solver employed here, without PCIL, the magnitude of spurious currents is so high that it is not possible to simulate this type of surface tension dominated flows, while using PCIL, we are able to simulate bubble rise and obtain results in close agreement with the experimental data.     

A numerical study of an annular liquid jet in a compressible gas medium

An annular liquid jet in a compressible gas medium has been examined using an Eulerian approach with mixed-fluid treatment. The governing equations have been solved by using highly accurate numerical methods. An adapted volume of fluid method combined with a continuum surface force model was used to capture the gas–liquid interface dynamics. The numerical simulations showed the existence of a recirculation zone adjacent to the nozzle exit and unsteady large vortical structures at downstream locations, which lead to significant velocity reversals in the flow field. It was found that the annular jet flow is highly unstable because of the existence of two adjacent shear layers in the annular configuration. The large vortical structures developed naturally in the flow field without external perturbations. Surface tension tends to promote the Kelvin–Helmholtz instability and the development of vortical structures that leads to an increased liquid dispersion. A decrease in the liquid sheet thickness resulted in a reduced liquid dispersion. It was identified that the liquid-to-gas density and viscosity ratios have opposite effects on the flow field with the reduced liquid-to-gas density ratio demoting the instability and the reduced liquid-to-gas viscosity ratio promoting the instability characteristics.     

Simulating gas–liquid multiphase flow using meshless Lagrangian method

A multiphase flow model has been established based on a moving particle semi‐implicit method. A surface tension model is introduced to the particle method to improve the numerical accuracy and stability. Several computational techniques are employed to simplify the numerical procedure and further improve the accuracy. A particle fraction multiphase flow model is developed and verified by a two‐phase Poiseuille flow. The multiphase surface tension model is discussed in detail, and an ethanol drop case is introduced to verify the surface tension model. A simple dam break is simulated to demonstrate the improvements with various modifications in particle method along with a new boundary condition. Finally, we simulate several bubble rising cases to show the capacity of this new model in simulating gas–liquid multiphase flow with large density ratio difference between phases. The comparisons among numerical results of mesh‐based model, experimental data, and the present model, indicate that the new multiphase particle method is acceptable in gas–liquid multiphase fluids simulation. Copyright © 2014 John Wiley & Sons, Ltd.