Smoothed-Particle Hydrodynamics Simulation of Ship Motion and Tank Sloshing under the Effect of Regular Waves

Predicting the response of liquefied natural gas(LNG)contained in vessels subjected to external waves is extremely important to ensure the safety of the transportation process.In this study,the coupled behavior due to ship motion and liquid tank sloshing has been simulated by the Smoothed-Particle Hydrodynamics(SPH)method.Firstly,the sloshing flow in a rectangular tank was simulated and the related loads were analyzed to verify and validate the accuracy of the present SPH solver.Then,a three-dimensional simplified LNG carrier model,including two prismatic liquid tanks and a wave tank,was introduced.Different conditions were examined corresponding to different wave lengths,wave heights,wave heading angles,and tank loading rates.Finally,the effects of liquid tank loading rate on LNG ship motions and sloshing loading were analyzed,thereby showing that the SPH method can effectively provide useful indications for the design of liquid cargo ships.

An Innovative Coupled Common-Node Discrete Element Method-Smoothed Particle Hydrodynamics Model Developed with LS-DYNA and Its Applications

In this study,a common-node DEM-SPH coupling model based on the shared node method is proposed,and a fluid–structure coupling method using the common-node discrete element method-smoothed particle hydrodynamics(DS-SPH)method is developed using LS-DYNA software.The DEM and SPH are established on the same node to create common-node DEM-SPH particles,allowing for fluid–structure interactions.Numerical simulations of various scenarios,including water entry of a rigid sphere,dam-break propagation over wet beds,impact on an ice plate floating on water and ice accumulation on offshore structures,are conducted.The interaction between DS particles and SPH fluid and the crack generation mechanism and expansion characteristics of the ice plate under the interaction of structure and fluid are also studied.The results are compared with available data to verify the proposed coupling method.Notably,the simulation results demonstrated that controlling the cutoff pressure of internal SPH particles could effectively control particle splashing during ice crushing failure.

A modified smoothed particle hydrodynamics method considering residual stress for simulating failure and its application in layered rock mass

Residual strength is an indispensable factor in evaluating rock fracture,yet the current Smoothed Particle Hydrodynamics(SPH)framework rarely considers its influence when simulating fracture.An improved cracking strategy considering residual stress in the base bond SPH method was proposed to simulate failures in layered rocks and slopes and verified by experimental results and other simulation methods(i.e.,the discrete element method).Modified Mohr–Coulomb failure criterion was applied to distinguish the mixed failure of tensile and shear.Bond fracture markψwas introduced to improve the kernel function after tensile damage,and the calculation of residual stress after the damage was derived after shear damage.Numerical simulations were carried out to evaluate its performance under different stress and scale conditions and to verify its effectiveness in realistically reproducing crack initiation and propagation and coalescence,even fracture and separation.The results indicate that the improved cracking strategy precisely captures the fracture and failure pattern in layered rocks and rock slopes.The residual stress of brittle tock is correctly captured by the improved SPH method.The improved SPH method that considers residual strength shows an approximately 13%improvement in accuracy for the safety factor of anti-dip layered slopes compared to the method that does not consider residual strength,as validated against analytical solutions.We infer that the improved SPH method is effective and shows promise for applications to continuous and discontinuous rock masses.

Application of Smoothed Particle Hydrodynamics(SPH)for the Simulation of Flow-Like Landslides on 3D Terrains

Flow-type landslide is one type of landslide that generally exhibits characteristics of high flow velocities,long jump distances,and poor predictability.Simulation of its propagation process can provide solutions for risk assessment and mitigation design.The smoothed particle hydrodynamics(SPH)method has been successfully applied to the simulation of two-dimensional(2D)and three-dimensional(3D)flow-like landslides.However,the influence of boundary resistance on the whole process of landslide failure is rarely discussed.In this study,a boundary condition considering friction is proposed and integrated into the SPH method,and its accuracy is verified.Moreover,the Navier-Stokes equation combined with the non-Newtonian fluid rheologymodel was utilized to solve the dynamic behavior of the flow-like landslide.To verify its performance,the Shuicheng landslide event,which occurred in Guizhou,China,was taken as a case study.In the 2D simulation,a sensitivity analysis was conducted,and the results showed that the shearing strength parameters have more influence on the computation accuracy than the coefficient of viscosity.Afterwards,the dynamic characteristics of the landslide,such as the velocity and the impact area,were analyzed in the 3D simulation.The simulation results are in good agreement with the field investigations.The simulation results demonstrate that the SPH method performs well in reproducing the landslide process,and facilitates the analysis of landslide characteristics as well as the affected areas,which provides a scientific basis for conducting the risk assessment and disaster mitigation design.

基于SPH方法的二维液舱晃荡特性研究

本文基于流固耦合动力学理论,建立二维液舱晃荡的光滑粒子流体动力学模型,并引入Dummy边界处理固定边界问题。首先,与实验结果对比验证模型的有效性;随后,分析不同横摇激励、多种载液率下自由液面的运动规律,研究液体抨击舱壁产生的压力特性;在此基础上,设计单一隔板和组合隔板等不同制荡结构,探究共振频率、30%载液率下,自由液面形态和壁面冲击压力的特性。结果表明:在共振频率下,横摇幅值对压力峰值的影响最大;随着载液率的增加,壁面抨击压力双峰特性逐渐消失;当隔板尺度较小时,T型隔板更有利于降低压力的整体水平;双T型隔板可有效降低舱底冲击。

Effects of hydrodynamics processes on phosphorus fluxes from sediment in large,shallow Taihu Lake

The turnover of phosphorus （P） in lake sediments, a major cause of eutrophication and subsequent deterioration of water quality, is in need of deep understanding. In this study, effects of resuspension on P release were studied in cylindrical microcosms with Yshape apparatus. The results indicated that there was a positive correlation between flux of suspended substance across sediment-water interface （Fss） and the wind speed, and an increasing Fss during each wind process followed by a steady state. The maximal Fss under fight, moderate, and strong wind conditions were 299.9±41.1,573.4±61.7, and 2093.8±215.7 g/m^2, respectively. However, flux of P across sediment-water interface （Fp） did not follow a similar pattern as Fss responding to wind intensity, which increased and reached the maximum in initial 120 rain for fight wind, then decreased gradually, with maximal flux of 9.4±1.9 mg/m^2. A rapid increase of Fp at the first 30 rain was observed under moderate wind, with maximal flux of 11.2±0.6 mg/m^2. Surprisingly, strong wind caused less Fp than under light and moderate wind conditions with maximal flux of 3.5±0.9 mg/m^2. Fss in water column declined obviously during the sedimentation process after winds, but Fp varied with wind regime. No obvious difference was detected on Fp after 8 h sedimentation process, compared with the initial value, which means little redundant P left in the water column after winds.

基于SPH露天台阶的爆破优化

针对露天深孔台阶爆破容易出现块度不均匀、大块率高等问题,运用经验公式法对台阶爆破参数进行计算并确定爆破方案,建立基于光滑粒子流体动力学(SPH)的台阶数值模型,对深孔台阶爆破孔排距共同作用下的岩体损伤情况及孔间应力大小进行研究。结果表明:孔排距会影响炮孔间的应力叠加作用,模拟过程中SPH粒子能良好地模拟出爆破过程中岩石运动状态和岩石损伤范围,当孔间距为8.5 m、排间距为5.0 m,爆破效果最佳。将模拟结果应用于现场试验,大块率较优化前降低了10.3%,块度平均合格率提高了5.2%,现场试验与模拟结果最优爆破方案效果高度一致,说明利用SPH的数值模拟方法研究台阶爆破是合理可行的。研究结果可为相关台阶爆破孔排距选择提供依据和参考。

基于SPH方法的空调进气口进水量分析及优化

由于雨刮臂摆动频次与导水槽及空调外循环进风口结构不匹配等因素,部分车型在淋雨工况下,有水滴甚至水流沿空调进风口进入空调系统,导致空调出现异味,甚至有水渗漏进入乘客舱,出现严重的感知质量问题。采用SPH方法,考虑雨刮臂运动、车身姿态及气流效应工况下,借助PreonLab软件仿真,复现某车型淋雨场景下的水管理质量问题,通过在空调外循环进风口位置增加挡板,显著减少了溅入空调箱的雨水量,并在实车淋雨试验中得到验证。方法能在车辆开发阶段发现并整改车辆在淋雨场景下的流通路径设计不合理等问题,可改变目前车辆淋雨场景水管理质量控制主要依托实车测试的现状,也可应用于车辆水管理的其他性能质量的控制,如涉水、雨污等。

基于SPH方法的船体海洋飞沫生成形态模拟

为了研究海洋飞沫的生成形态,采用光滑粒子法(SPH)的粒子捕捉技术,分析4种常见的船首特征对海洋飞沫轨迹的影响。研究发现,船首的设计在飞沫的生成和运动特性上起着决定性作用。球鼻艏和破冰艏通过镇压波浪减少飞沫的覆盖范围,而飞剪艏和前倾艏则利用其独特的流线型设计有效地对飞沫进行引流,但却导致飞沫的生成高度较高,覆盖范围较广。此外,研究还揭示了不同类型的飞沫(如崩波破碎、卷波和激散波破碎)在船首附近的工况和轨迹模式。研究结果可为船舶设计优化及防除冰技术提供重要理论支持。

基于SPH算法的热仿真系统开发

随着电子设备日益微型化和集成化,热仿真已成为其设计中的关键因素。电子封装模块的热仿真通常使用传统的有限元法FEM(finite element method),存在计算效率和精度之间的矛盾,在处理大变形问题和网格畸变方面也容易造成计算不收敛,从而导致结果错误。针对该问题,提出一种基于光滑粒子动力学SPH(smoothed particle hydrodynamics)算法的电子封装模块热仿真系统。该算法基于无网格拉格朗日数值方法,通过将热仿真对象离散为1组粒子的方式求解热传导方程,从而准确地预测电子封装模块的传热与散热,无需生成并处理大量的微小网格,不用担心网格失真等问题。SPH相对于FEM,仿真精度误差保持在1%~2%,仿真效率可提升近30倍,适合用于复杂和动态系统的模拟仿真。

基于SPH粒子法的汽车高压油箱碰撞CAE计算分析

文章基于汽车车体结构、高压燃油箱及其周围零部件组成的子模型,采用LS-DYNA显式求解器,利用SPH 粒子技术,针对整车正面碰撞、侧面碰撞以及后面碰撞过程中,高压油箱周围变形以及其内部油液的晃动对油箱壳体的影响进行有限元模拟和研究。通过碰撞CAE计算分析,得出碰撞过程中油液晃动对油箱的撞击而产生的应力分布及变化过程,分析油箱的变形、应力、应变分布等情况,可以评估汽车在发生高速碰撞时油箱损伤、破坏情况。最终通过在某整车高速追尾CAE模拟中的应用,再现了油液SPH 粒子的有效作用。

基于自适应FEM-SPH耦合算法的飞机典型部位破片冲击战伤的数值研究

针对飞机典型部位在遭到高速破片攻击后结构整体的战伤状态及破片的剩余行为开展数值模拟。应用LS-DYNA软件,结合有限单元方法(finite element method,FEM)和光滑粒子流体动力学(smoothed particle hydrodynamics,SPH)两者的优势,建立自适应的FEM-SPH耦合模拟方法,并构建2种飞机典型部位的计算模型,采用六面体网格局部细化方法实现了核心位置的精确模拟,并进行试验来验证数值模型;开展了一系列高速冲击战伤模拟,对比了不同工况下破片高速冲击结构后形成的碎片云和破口形貌,并对破片的剩余速度和质量进行分析,确定了破片在结构蒙皮上的临界跳飞角。结果表明:自适应FEM-SPH耦合算法的计算结果与试验结果吻合良好,能够对破片高速冲击战伤进行有效准确模拟;碎片云分布形状随破片速度增加变得狭长,冲击角度会改变碎片云和结构破口形状朝向;碎片云高度和扩散速度随破片速度或角度的变化趋势基本一致并都呈线性关系;破片的速度减少量不随初始速度变化,质量减少量则与冲击速度正相关,两者与冲击角度都负相关;破片临界跳飞角与冲击速度大小基本呈线性关系。

Influence of Hydrodynamic Pore Pressure Damage on the Performance of Hot-Mixed Renewable Asphalt Mixture

For evaluating the water stability of hot-mixed renewable asphalt mixture(HRM),the traditional methods are all tested under still water conditions.Except for damage in still water conditions,the hydrodynamic pore pressure generated by the tire driving on the surface water has a great impact.Thus,the RAP contents of the HRMs were designed at 0%,30%,45%and 60%with AC-25 gradation.Then,the self-designed evaluation methods of water stability and dynamic modulus were studied.Finally,the mechanism of the influence of hydrodynamic pore pressure damage on HRMs was studied.The results show that the water stability of HRM containing 30%RAP is equivalent to that of 45%RAP,and the water stability of HRM containing 60%RAP decreases significantly.The Contabro test after MIST treatment can be used as an evaluation method for hydrodynamic pore pressure damage on HRM.Low-speed,heavy-load traffic and larger RAP content have greater damage to the mixture after hydrodynamic pore pressure damage.The performance differences between the aged bitumen and pure bitumen,as well as the aged minerals and new minerals,are continuing to be enlarged in hydrodynamic pore pressure conditions,finally affecting the water stability and dynamic modulus of the HRMs.

Numerical and experimental study on the particle erosion and gas–particle hydrodynamics in an integral multi-jet swirling spout-fluidized bed

In this paper,using the computational fluid dynamics based on Euler Lagrange and the commercial software Barracuda VR,the gas-particle hydrodynamics and the erosion of particles on the inner wall and internal components of the spouted bed in the integrated multi-jet swirling spout-fluidized bed(IMSSFB)are studied.Erosion experiments have obtained the characterization of particle erosion on internal components and verified the relevant numerical models.The results show that:the particle distribution within the IMSSFB is uneven due to the cyclonic effect of the axial swirl vane(ASV),resulting in particle erosion for the ASV being concentrated on one side;when the gas reaches the top,too high an erosion gas velocity leads to gas backflow.As the filling height increases,there is a tendency for the erosion position of the particles on the ASV to expand upwards.However,the effect of increasing gas velocity on the erosion position is insignificant.

基于PCISPH的流体粒子飞溅改进方法

流体飞溅是自然界中最常见的流体现象。流体模拟是计算机图形学的一个重要研究分支。流体模拟已经广泛应用于电影、游戏和其他工业数值领域。由于飞溅场景中流体的密度和压力的速率变化非常大,因此模拟对离散化解的精度要求很高。对于流体粒子飞溅模拟数值不稳定,且缺乏真实效果,针对经典的预测校正不可压缩SPH(PCISPH),该文提出了基于粒子表面流体飞溅改进的泊松压力解的方法,通过替换压力源项提高了模拟精度。实验结果表明,该方法中流体的压力分布较之经典方法更接近现实的流体粒子飞溅模拟效果,并且改进后的方法经过剧烈运动后帧率仍符合实时性仿真效果的要求,保证了良好的压力稳定性和真实感。该文通过Anaconda3集成Python3和Taichi环境实现流体仿真实验。

基于SPH-FEM耦合方法的泥石流冲击输电塔基础的动力分析

泥石流是我国西南山区常见的地质灾害。架空输电杆塔在泥石流的冲击下往往发生基础破坏甚至会造成杆塔倒塌。首先采用光滑粒子流体动力学(smoothed particle hydrodynamics,简称SPH)方法和有限元方法(finite element method,简称FEM)相耦合的三维数值方法模拟了泥石流对杆塔基础的冲击作用;在与相关模型试验结果验证的基础上,开展了不同泥石流密度、黏度系数及初始速度条件下对输电塔基础的冲击力作用的参数分析;研究结果表明:随着泥石流初始速度的增加,冲击力峰值会随之增大;前排基础的冲击力峰值均大于后排基础;泥石流冲击过程特性受到泥石流密度和黏度系数影响。与稀性泥石流相比:黏性泥石流冲击基础后,基础下游真空区相对要小;此外,将数值模拟结果与Kwan冲击力公式及铁二院推荐的冲击压力设计公式预测值进行对比分析可以发现:Kwan冲击力公式能较好地预测出基础所受泥石流冲击力的平均趋势,最大预测误差低于30%,铁二院公式预测的稀性和黏性泥石流的冲击压力平均偏低分别约17%和28%。相关研究结果有望为泥石流频发区域输电塔基础的设计和风险评估提供一定的参考依据。

Collisionless shock acceleration of protons in a plasma slab produced in a gas jet by the collision of two laser-driven hydrodynamic shockwaves

We have recently proposed a new technique of plasma tailoring by laser-driven hydrodynamic shockwaves generated on both sides of a gas jet[Marquès et al.,Phys.Plasmas 28,023103(2021)].In a continuation of this numerical work,we study experimentally the influence of the tailoring on proton acceleration driven by a high-intensity picosecond laser in three cases:without tailoring,by tailoring only the entrance side of the picosecond laser,and by tailoring both sides of the gas jet.Without tailoring,the acceleration is transverse to the laser axis,with a low-energy exponential spectrum,produced by Coulomb explosion.When the front side of the gas jet is tailored,a forward acceleration appears,which is significantly enhanced when both the front and back sides of the plasma are tailored.This forward acceleration produces higher-energy protons,with a peaked spectrum,and is in good agreement with the mechanism of collisionless shock acceleration(CSA).The spatiotemporal evolution of the plasma profile is characterized by optical shadowgraphy of a probe beam.The refraction and absorption of this beam are simulated by post-processing 3D hydrodynamic simulations of the plasma tailoring.Comparison with the experimental results allows estimation of the thickness and near-critical density of the plasma slab produced by tailoring both sides of the gas jet.These parameters are in good agreement with those required for CSA.

基于SPH方法的沉箱开孔位置对消浪效果的影响

为了更深入地理解沉箱结构对波浪传播和反射的影响,对波浪与开孔沉箱相互作用进行数值模拟研究,为探究消浪效果最佳的沉箱上部开孔位置,通过基于SPH方法建立的数值波浪水槽,分析了3种不同开孔位置对反射系数Kr、波面高度变化、涡量场及流场的影响,总结了不同开孔位置对消浪效果的影响。研究发现当上部开孔的顶端在水下0.08 m时沉箱反射系数Kr最小,沉箱外的波面高度变化最小,开孔沉箱开孔处的涡量消耗的总量最大,流场变化最大是其他2种开孔位置的1~2倍,同时沉箱内波面最为平静,消浪效果最好。所得结论对开孔沉箱实际工程中消浪效果研究具有一定的参考意义。

Modeling Near-Field Impulsive Waves Generated by Deformable Landslide Using the HBP Model Based on the SPH Method

Landslide-generated impulsive waves(LGWs)in reservoir areas can seriously threaten waterway safety as well as hu-man life and properties around the two side slopes.The risk reduction and mitigation of such a hazard require the accurate prediction of near-field wave characteristics,such as wave amplitude and run-up.However,near-field LGW involves complicated fluid-solid interactions.Furthermore,the wave characteristics are closely related to various aspects,including the geometry and physical features of the slide,river,and body of water.However,the empirical or analytical methods used for rough estimation cannot derive accurate results,especially for deformable landslides,due to their significant geometry changes during the sliding process.In this study,the near-field waves generated by deformable landslides were simulated by smoothed particle hydrodynamics(SPH)based on multi-phase flow.The deformable landslides were generalized as a kind of viscous flow by adopting the Herschel-Bulkley-Papanastasiou(HBP)-based nonNewtonian rheology model.The HBP model is capable of producing deformable landslide dynamics even though the high-speed sliding process is involved.In this study,an idealized experiment case originating from Lituya LGW and a practical case of Gongjiafang LGW were reproduced for verification and demonstration.The simulation results of both cases show satisfactory consistency with the experiment/investigation data in terms of landslide movement and near-field impulsive wave characteristics,thus indicating the applicability and accuracy of the proposed method.Finally,the effects of the HBP model’s rheological parameters on the landslide dynamics and near-field wave characteristics are discussed,providing a parameter calibration method along with sug-gestions for further applications.

基于SPH方法的围油栏仿真研究

为解决计算流体动力学(computational fluid dynamics,CFD)方法在模拟围油栏围控溢油时流-固耦合、水-油液体两相流中的不稳定性、计算效率差、将自由表面视为滑移壁且不考虑围油栏移动等问题,基于光滑粒子流体动力学(smoothed particle hydrodynamics,SPH)方法建立适用于围油栏围控溢油的固-液耦合、液体两相流数值模型,将5种不同裙摆结构的围油栏导入两相液体相互作用的数值仿真模型中,通过SPH代码设定数值模拟相关参数,提高计算效率,得到较精确的围油栏滞油长度与溢油失效步长。结果表明:围油栏围控溢油能力与上、下裙摆高度比及裙摆角类型有关,上裙摆高度占比越大,围油栏的滞油表现越好;前折角型围油栏的滞油表现优于前转角型围油栏;滞油长度与油品性质有关,油品的密度和黏度越小,围油栏的滞油长度越大。