三维溃坝流的光滑粒子动力学方法模拟

采用光滑粒子动力学SPH(Smoothed Particle Hydrodynamics)方法对三维溃坝流问题进行了数值模拟。为了逼真地模拟出坝内水体与壁面间相互作用而产生的水花飞溅、融合以及近壁面流动等现象,加入了混合长度形式的湍流模型。为了有效地防止粒子穿透固壁,提出了一种新型的适合三维数值模拟的固壁边界处理方法。应用SPH方法对三维溃坝流进行了数值模拟,并分别考虑了未添加障碍物和添加圆柱障碍物两种情形。计算结果表明,改进SPH方法能够精细地捕捉溃坝流在不同时刻的自由液面,并获得稳定而精确的数值结果。     

一种SPH应力修正算法及自由表面流中的应用

提出了一种SPH应力修正算法,即模型中的拉应力和压应力分别采用不同的插值核函数和状态方程来处理,改善应力稳定性问题。介绍了一种改进的Quintic核函数,用于改善模型中压应力的稳定性。通过增加钟型核函数的光滑长度,改善模型中拉应力的稳定性。采用该应力修正算法模拟了无重力条件下方形液滴的震荡变形过程,对比分析了不同算法的模拟结果。此外,为进一步验证算法的适用性,模拟了溃坝算例。研究表明,改进的Quintic型核函数明显改善了粒子聚集现象,该SPH应力修正方法可以使液滴具有更均匀的粒子分布以及更光滑的自由表面,有效改善了SPH方法中的压应力不稳定作用以及自由表面流的模拟精度。     

近自由面水下爆炸冲击载荷特性三维数值模拟

基于三维SPH 方法,对传统链表搜索算法进行了变光滑长度改进,并提出了具有较好稳定性的多 相物质交界面的处理方法,模拟了三维无限域水下爆炸问题,验证了改进的三维SPH 方法模拟水下爆炸问题 的可行性和有效性。在此基础上,建立了水下爆炸三维数值模型,模拟了近自由面水下爆炸过程,研究了冲击 波传播特征、自由面下压力场和能量场特性以及水柱的产生过程。结果表明:自由面可将冲击波压力峰值和 压力冲量最大衰减到1/3和1/7;爆深的增加会导致压力比和冲量比的等值柱面曲率变小,产生的水柱也逐 渐由破碎的喷柱向高而窄的水冢过渡。     

二维空间中基于黎曼解的SPH方法

应用传统的光滑粒子流体动力学(SPH)方法模拟激波问题能自主捕获波阵面,但是在接触间断界面处却不可避免地存在压力振荡。采用黎曼解修正粒子对之间的相互作用,并将基于黎曼解的SPH方法扩展到二维空间。通过对一维具有解析解的激波管问题和二维拟内爆问题的模拟,比较了基于黎曼解的SPH方法和传统的SPH方法。数值结果表明,相对于传统的SPH方法,基于黎曼解的SPH方法能够有效避免接触界面处的压力振荡,提高求解的精度。     

带有障碍物溃坝流动问题的有限元数值模拟研究

针对下游带有障碍物的溃坝流动问题,本文基于两相流动模型,在有限元算法框架下对其进行数值模拟研究。依据水平集(Level Set)方法追踪运动界面,并引入了一个简单的修正技术,保证较好的质量守恒性。为了精确表示运动界面,采用稳定和有效的间断有限元方法求解双曲型Level Set及其重新初始化方程。对于两相统一Navier-Stokes方程,首先利用分裂格式对其解耦,然后通过SUPG (Streamline Upwind Petrov Galerkin)方法进行数值求解。模拟研究了下游带有障碍物的牛顿流体溃坝流动问题,得到的数值结果与文献已有模拟结果及实验结果均吻合较好。此外,还考虑了幂律型非牛顿流体,并分析了不同特性非牛顿流体对于溃坝流动过程和界面形态等的影响。     

基于ALE,CEL和SPH方法的球形破片高速冲击充液结构对比研究

为研究ALE,CEL和SPH方法在高速冲击流固耦合动力学数值分析中的差异性,开展球形破片高速冲击充液结构数值模拟研究。建立经文献资料验证的ALE,CEL和SPH三种动力学模型,研究了流体压力变化、形成的空腔尺寸、破片速度衰减变化和充液结构变形等模拟精度,并分析相应的计算成本。结果表明,ALE,CEL和SPH三种方法均能有效模拟破片高速冲击充液结构的流固耦合动力学过程;ALE方法预测的空腔尺寸精度较高;CEL方法预测的流体压力、破片速度衰减和充液结构变形精度较高;SPH方法预测的空腔尺寸、破片速度衰减精度较高;当网格尺寸一致时, SPH方法计算时长约为ALE和CEL方法的两倍,但SPH方法前后处理更加简便。     

基于修正光滑粒子流体动力学算法的低能量耗散数值波浪水池开发

目前, 无网格光滑粒子流体动力学SPH粒子法在波浪与结构物相互作用研究方面得到广泛应用, 但该方法模拟波浪远距离传播时, 常常面临严重的能量耗散问题, 导致波高非物理性降低, 给大范围海域、长时间作用下的波-物耦合作用研究带来一定困难. 对此, 本文采用一种核函数修正算法, 在确保粒子间相互作用对称性的同时, 改进压力梯度离散项的计算精度, 设法解决SPH方法中能量非物理性耗散的难题. 相较于前人减缓能量非物理性衰减的方法, 本文的修正SPH算法避免了自由液面搜索等复杂处理过程, 并能保证动量守恒特性. 数值结果中, 采用振荡液滴、规则波、不规则波等算例, 验证本修正SPH算法的准确性和有效性. 结果表明, 该修正SPH算法能准确模拟振荡液滴形态变化, 且动能保持较好守恒性. 通过数值水池与物理水池两者规则波与不规则波结果的对比分析表明, 基于本文修正SPH算法建立的数值波浪水池具有较好的抗能量衰减效果, 能实现长时间、远距离波浪传播的准确模拟. 此外, 本算法能在低光滑长度系数条件下, 实现精确模拟, 将极大缩减三维SPH模拟的时间, 从而节约计算成本.      

混凝土中爆炸模拟的数值方法比较

应用AUTODYN商业软件对装药在混凝土靶介质中的爆炸过程进行数值模拟,对比了混凝土中爆炸模拟的不同方法。将各数值方法模拟得到的混凝土毁伤特征尺寸与相关实验数据进行比较,研究了各数值方法的准确性及描述混凝土毁伤特征的可靠性参数。结果发现：Euler算法计算的压力最大,SPH和Lagrange耦合算法计算的压力最小;Euler算法和SPH算法计算的压力衰减过程会发生比较显著的扰动。以不同的平均网格尺寸对数值模拟结果的影响分析了各种数值方法的稳定性。比较靶体上同一测点的压力时间历程表明网格变化对测点处压力影响不大。通过应用各种数值方法模拟再现混凝土中的爆炸过程,揭示了各种方法的优缺点,明确了现有商业软件模拟混凝土中爆炸的适用性。     

基于SPH方法的聚能射流侵彻混凝土靶板数值模拟

在完全变光滑长度SPH(smoothed particle hydrodynamics)方法的基础上,利用F.Ott等提出的修正SPH方法处理在求解多介质大密度问题时的数值不稳定性问题,运用Holmquist-Johnson-Cook本构模型处理混凝土在冲击载荷下的变形和损伤问题,对聚能装药射流侵彻混凝土靶板的过程进行了数值模拟,同时利用LS-DYNA非线性有限元程序进行对比,分析了2种方法得到的混凝土von Mises应力变化、射流头部特定节点处的速度变化及裂纹演变,验证了SPH方法的准确性。分析了另外2种不同尺寸的靶板在射流侵彻作用下的破坏形式,结果符合射流侵彻物理规律,表明该方法适合模拟聚爆炸与冲击等大变形破坏等问题。     

一种基于黎曼解处理大密度比多相流SPH的改进算法

多相流界面存在密度、黏性等物理场间断,直接采用传统光滑粒子水动力学(smoothedparticle hydrodynamics,SPH)方法进行数值模拟,界面附近的压力和速度存在震荡.一套基于黎曼解能够处理大密度比的多相流SPH计算模型被提出,该模型利用黎曼解在处理接触间断问题方面的优势,将黎曼解引入到SPH多相流计算模型中,为了能够准确求解多相流体物理黏性、减小黎曼耗散,对黎曼形式的SPH动量方程进行了改进,又将Adami固壁边界与黎曼单侧问题相结合来施加多相流SPH固壁边界,同时模型中考虑了表面张力对小尺度异相界面的影响,该模型没有添加任何人工黏性、人工耗散和非物理人工处理技术,能够反应多相流真实物理黏性和物理演变状态.采用该模型首先对三种不同粒子间距离散下方形液滴震荡问题进行了数值模拟,验证了该模型在处理异相界面的正确性和模型本身的收敛性;后又通过对Rayleigh--Taylor不稳定、单气泡上浮、双气泡上浮问题进行了模拟计算,结果与文献对比吻合度高,异相界面捕捉清晰,结果表明,本文改进的多相流SPH模型能够稳定、有效的模拟大密度比和黏性比的多相流问题.      

Modified incompressible SPH method for simulating free surface problems

An incompressible smoothed particle hydrodynamics (I-SPH) formulation is presented to simulate free surface incompressible fluid problems. The governing equations are mass and momentum conservation that are solved in a Lagrangian form using a two-step fractional method. In the first step, velocity field is computed without enforcing incompressibility. In the second step, a Poisson equation of pressure is used to satisfy incompressibility condition. The source term in the Poisson equation for the pressure is approximated, based on the SPH continuity equation, by an interpolation summation involving the relative velocities between a reference particle and its neighboring particles. A new form of source term for the Poisson equation is proposed and also a modified Poisson equation of pressure is used to satisfy incompressibility condition of free surface particles. By employing these corrections, the stability and accuracy of SPH method are improved. In order to show the ability of SPH method to simulate fluid mechanical problems, this method is used to simulate four test problems such as 2-D dam-break and wave propagation.     

A comparative study of truly incompressible and weakly compressible SPH methods for free surface incompressible flows

In this paper, the performance of the incompressible SPH (ISPH) method and an improved weakly compressible SPH (IWCSPH) method for free surface incompressible flows are compared and analyzed. In both methods, the Navier–Stokes equations are solved, and no artificial viscosity is used. The ISPH algorithm in this paper is based on the classical SPH projection method with common treatments on solid boundaries and free surfaces. The IWCSPH model includes some advanced corrective algorithms in density approximation and solid boundary treatment (SBT). In density approximation, the moving least squares (MLS) approach is applied to re‐initialize density every several steps to obtain smoother and more stable pressure fields. An improved coupled dynamic SBT algorithm is implemented to obtain stable pressure values near solid wall areas and, thus, to minimize possible numerical oscillations brought in by the solid boundaries. Three representative numerical examples, including a benchmark test for hydrostatic pressure, a dam breaking problem and a liquid sloshing problem, are comparatively analyzed with ISPH and IWCSPH. It is demonstrated that the present IWCSPH is more attractive than ISPH in modeling free surface incompressible flows as it is more accurate and more stable with comparable or even less computational efforts. Copyright © 2013 John Wiley & Sons, Ltd.     

A consistent and fast weakly compressible smoothed particle hydrodynamics with a new wall boundary condition

A modified weakly compressible smoothed particle hydrodynamics (WCSPH) is presented, which utilizes consistent discretization schemes for spatial derivatives in the flow equations. Here, each SPH particle is considered as a computational point that represents a specific part of the fluid. To overcome non‐physical oscillations that usually arise in standard WCSPH, we modified the mass conservation equation by using a numerical filter. This modification is based on the difference between two discretization schemes used for the term . Furthermore, a new implementation of wall boundary condition in SPH is introduced. This condition is imposed on the pressure of wall boundary particles to ensure that the acceleration of each boundary particle in normal direction to the wall is zero. Thus, no penetration through walls will occur. To examine the performance of the modified method, we solved a series of two‐dimensional incompressible internal flow benchmark problems. By comparing the result with analytical solutions and the results of the standard WCSPH, we show that the use of consistent schemes in conjunction with the proposed numerical filter improves both accuracy and speed of the numerical method. Copyright © 2011 John Wiley & Sons, Ltd.     

Incompressible SPH simulation of wave breaking and overtopping with turbulence modelling

In this paper a truly incompressible version of the smoothed particle hydrodynamics (SPH) method is presented to investigate the surface wave overtopping. SPH is a pure Lagrangian approach which can handle large deformations of the free surface with high accuracy. The governing equations are solved based on the SPH particle interaction models and the incompressible algorithm of pressure projection is implemented by enforcing the constant particle density. The two‐equation k–ε model is an effective way of dealing with the turbulence and vortices during wave breaking and overtopping and it is coupled with the incompressible SPH numerical scheme. The SPH model is employed to reproduce the experiment and computations of wave overtopping of a sloping sea wall. The computations are validated against the experimental and numerical data found in the literatures and good agreement is observed. Besides, the convergence behaviour of the numerical scheme and the effects of particle spacing refinement and turbulence modelling on the simulation results are also investigated in further detail. The sensitivity of the computed wave breaking and overtopping on these issues is discussed and clarified. Copyright © 2005 John Wiley & Sons, Ltd.     

Finite Difference Method for the Acoustic Radiation of an Elastic Plate Excited by a Turbulent Boundary Layer: A Spectral Domain Solution

A finite difference method is developed to study, on a two-dimensional model, the acoustic pressure radiated when a thin elastic plate, clamped at its boundaries, is excited by a turbulent boundary layer. Consider a homogeneous thin elastic plate clamped at its boundaries and extended to infinity by a plane, perfectly rigid, baffle. This plate closes a rectangular cavity. Both the cavity and the outside domain contain a perfect fluid. The fluid in the cavity is at rest. The fluid in the outside domain moves in the direction parallel to the system plate/baffle with a constant speed. A turbulent boundary layer develops at the interface baffle/plate. The wall pressure fluctuations in this boundary layer generates a vibration of the plate and an acoustic radiation in the two fluid domains. Modeling the wall pressure fluctuations spectrum in a turbulent boundary layer developed over a vibrating surface is a very complex and unresolved task. Ducan and Sirkis [1] proposed a model for the two-way interactions between a membrane and a turbulent flow of fluid. The excitation of the membrane is modeled by a potential flow randomly perturbed. This potential flow is modified by the displacement of the membrane. Howe [2] proposed a model for the turbulent wall pressure fluctuations power spectrum over an elastomeric material. The model presented in this article is based on a hypothesis of one-way interaction between the flow and the structure: the flow generates wall pressure fluctuations which are at the origin of the vibration of the plate, but the vibration of the plate does not modify the characteristics of the flow. A finite difference scheme that incorporates the vibration of the plate and the acoustic pressure inside the fluid cavity has been developed and coupled with a boundary element method that ensures the outside domain coupling. In this paper, we focus on the resolution of the coupled vibration/interior acoustic problem. We compare the results obtained with three numerical methods: (a) a finite difference representation for both the plate displacement and the acoustic pressure inside the cavity; (b) a coupled method involving a finite difference representation for the displacement of the plate and a boundary element method for the interior acoustic pressure; (c) a boundary element method for both the vibration of the plate and the interior acoustic pressure. A comparison of the numerical results obtained with two models of turbulent wall pressure fluctuations spectrums - the Corcos model [3] and the Chase model [4] - is proposed. A difference of 20 dB is found in the vibro-acoustic response of the structure. In [3], this difference is explained by calculating a wavenumber transfer function of the plate. In [6], coupled beam-cavity modes for similar geometry are calculated by the finite difference method. This revised version was published online in July 2006 with corrections to the Cover Date.     

Incompressible smoothed particle hydrodynamics simulation of multifluid flows

This paper presents an incompressible SPH (ISPH) technique to simulate multifluid flows. The SPH method is a mesh‐free particle modeling approach that can treat free surfaces and multi‐interfaces in a simple and efficient manner. The ISPH method employs an incompressible hydrodynamic formulation to solve the fluid pressure that ensures a stable pressure field. Two multifluid ISPH models are proposed following different interface treatments: the coupled ISPH model does not distinguish the different fluid phases and applies the standard ISPH technique across the interface, whereas the decoupled ISPH model first treats each fluid phase separately and then couples the different phases by applying pressure and shear stress continuities across the interface. The two proposed models were used to investigate a gravity underflow with a low density ratio in a Generalized Reservoir Hydrodynamics (GRH) flume and a horizontal lock exchange flow with a high density ratio. Comparisons with data and relevant numerical error analysis indicated that the decoupled model performed well in cases of both low and high density ratios because of the accurate treatment of interface boundaries. Copyright © 2011 John Wiley & Sons, Ltd.     

改进的SPH边界处理方法与土体大变形模拟

沙土滑坡往往会造成重大的人身财产损失,研究这类土体大变形问题对防灾工程具有指导意义。光滑粒子流体动力学SPH（Smoothed Particle Hydrodynamics）方法是一种拉格朗日型无网格粒子法,十分适用于模拟大变形问题。在SPH方法中,合适的边界处理方法一直是个难点,传统的边界虚粒子法或排斥力法较难模拟复杂边界。本文引入了一种能处理任意形状边界的方法——统一半解析壁面边界条件处理方法USAW（unified semi-analytical wall boundary conditions）,通过在控制方程中引入修正因子并保留边界面积分项来弥补边界缺失。为了更准确模拟问题域边界,提出无质量边界粒子的新概念。利用该方法成功模拟了土体滑坡算例,验证了方法的可靠性,并避免了边界零粒子层问题。通过数值模拟,分析了内摩擦角和黏聚力等土体物性参数对滑坡过程的影响。最后,应用该方法研究了滑坡冲击楔形体时的压力响应。     

Simulation of container filling process with two inlets by improved smoothed particle hydrodynamics (SPH) method

In this article, an improved smoothed particle hydrodynamics (SPH) method is proposed to simulate the filling process with two inlets. Improvements are achieved by deriving a corrected kernel gradient of SPH and a density re-initialisation. In addition, a new treatment of solid wall boundaries is presented. Thus, the improved SPH method has higher accuracy and better stability, and conserves both linear and angular momentums. The validity of the new boundary treatment is shown by simulating the spin-down problem. The bench tests are also presented to demonstrate the performance of the improved SPH method. Then the filling process with a single inlet is simulated to show the ability to capture complex-free surface of the proposed method. Finally, the filling process with two inlets is numerically investigated. The numerical results show that the filling patterns are affected significantly by Reynolds number, aspect ratio of the container and the location of the inlets.