基于黎曼解的粒子间接触算法在SPH中的应用

 采用光滑粒子流体动力学（SPH）方法模拟大变形问题具有明显的优点,但传统的SPH方法在模拟冲击波与接触界面的作用时,往往会出现压力的反常跳动。采用黎曼解描述粒子间相互作用的接触算法对传统SPH方法进行修正,计算了激波管和飞片碰撞（包含接触界面）问题中波的传播,并将计算结果与解析解作比较。结果表明,与传统的光滑粒子法相比,该改进的光滑粒子法无需引入人工粘性项和人工热流项,程序结构简洁,且能较好地处理接触界面问题,从而能有效提高计算精度。     

基于光滑粒子动力学方法的液滴冲击固壁面问题数值模拟

采用改进的光滑粒子动力学(SPH)方法对液滴冲击固壁面问题进行了数值模拟. 为了提高传统SPH方法的计算精度和数值稳定性, 在传统的SPH方法的基础上对粒子方法中的密度和核梯度进行了修正, 采用了考虑黎曼解法的SPH流体控制方程, 构造了一种新型的粒子间相互作用力(IIF)模型来模拟表面张力的影响. 应用改进的SPH方法对液滴冲击固壁面问题进行了数值模拟. 计算结果表明：新型的IIF 模型能够较好地模拟表面张力的影响, 改进的SPH方法能够精细地描述液滴与固壁面相互作用过程中液滴的内部压力场演变和自由面形态变化, 液滴的铺展因子随初始韦伯数的增大而增大, 数值模拟结果与实验得到的结果基本一致. 关键词： 液滴 固壁面 光滑粒子动力学 表面张力     

Symplectic and multisymplectic numerical methods for Maxwell’s equations

In this paper, we compare the behaviour of one symplectic and three multisymplectic methods for Maxwell’s equations in a simple medium. This is a system of PDEs with symplectic and multisymplectic structures. We give a theoretical discussion of how some numerical methods preserve the discrete versions of the local and global conservation laws and verify this behaviour in numerical experiments. We also show that these numerical methods preserve the divergence. Furthermore, we extend the discussion on dispersion for (multi)symplectic methods applied to PDEs with one spatial dimension, to include anisotropy when applying (multi)symplectic methods to Maxwell’s equations in two spatial dimensions. Lastly, we demonstrate how varying the Courant–Friedrichs–Lewy (CFL) number can cause the (multi)symplectic methods in our comparison to behave differently, which can be explained by the study of backward error analysis for PDEs.     

三维SPH程序及其在斜高速碰撞问题的应用

叙述了平滑粒子流体动力学(简称SPH)的基本概念,重新推导了所用的基本方程。介绍了具体计算格式、算法以及策略等。对所编制的程序作了一定的检验。对SPH方法中的人为粘性和人为热流的作用作了数值考查。还提供了一个三维高速碰撞现象的算例。最后指出尚存在的若干问题。     

Robust particle methods for relativistic hydrodynamics

In this paper we describe how the equations of relativistic fluid dynamics can be solved numerically using the particle method SPH.     

液滴溅落问题的光滑粒子动力学模拟

对传统的光滑粒子动力学方法进行了改进, 改进的光滑粒子动力学方法对传统粒子方法中的核近似式和粒子近似式进行了修正, 采用Riemann 算法求解光滑粒子动力学流体控制方程, 添加了表面张力的计算程序, 考虑了表面张力对液滴溅落的影响. 应用改进的光滑粒子动力学方法对液滴静止状态下冲击液面的飞溅过程进行了数值模拟. 计算结果表明, 改进的光滑粒子动力学方法能够有效地描述液滴溅落液面的动力学特性和自由表面变化特征, 能够得到稳定精度的结果.     

Numerical simulations of instabilities in the implosion process of inertial confined fusion in 2D cylindrical coordinates

In this paper, we introduce a multi-material arbitrary Lagrangian and Eulerian method for the hydrodynamic radiative multi-group diffusion model in 2D cylindrical coordinates. The basic idea in the construction of the method is the following: In the Lagrangian step, a closure model of radiation-hydrodynamics is used to give the states of equations for materials in mixed cells. In the mesh rezoning step, we couple the rezoning principle with the Lagrangian interface tracking method and an Eulerian interface capturing scheme to compute interfaces sharply according to their deformation and to keep cells in good geometric quality. In the interface reconstruction step, a dual-material Moment-of-Fluid method is introduced to obtain the unique interface in mixed cells. In the remapping step, a conservative remapping algorithm of conserved quantities is presented. A number of numerical tests are carried out and the numerical results show that the new method can simulate instabilities in complex fluid field under large deformation,and are accurate and robust.     

Direct Numerical Simulation of the Rayleigh-Taylor Instability with the Spectral Element Method

A novel method is proposed to simulate Rayleigh-Taylor instabilities using a specially-developed unsteady three-dimensional high-order spectral element method code. The numerical model used consists of Navier-Stokes equations and a transport-diffusive equation. The code is first validated with the results of linear stability perturbation theory. Then several characteristics of the Rayleigh-Taylor instabilities are studied using this three-dimensional unsteady code, including instantaneous turbulent structures and statistical turbulent mixing heights under different initial wave numbers. These results indicate that turbulent structures of Rayleigh-Taylor instabilities are strongly dependent on the initial conditions. The results also suggest that a high-order numerical method should provide the capability of simulating small scale fluctuations of Rayleigh-Taylor instabilities of turbulent flows.     

基于不可压缩光滑粒子动力学的黏性液滴变形过程仿真

应用基于投影算法的不可压缩光滑粒子动力学(incompressible smoothed particle hydrodynamics, ISPH)法对黏性液滴变形过程进行了数值仿真. 对于张力失稳导致的粒子非物理簇集问题, 采用粒子移位技术加以解决. 为了验证本文ISPH 算法的精度和稳定性, 分别模拟了圆形黏性液滴的拉伸变形过程以及方形液滴的旋转变形过程, 得到了不同时刻液滴内部的压力变化特征, 准确地捕捉了液滴自由面演化过程, 并将数值计算结果与文献中的解析解进行了比较.分析结果表明, 基于投影算法的不可压缩光滑粒子动力学方法结合粒子移位技术, 能够有效地模拟黏性液滴变形过程, 可以得到精确和稳定的结果. 关键词： 不可压缩光滑粒子动力学 黏性液滴 自由面流动 数值仿真     

基于弱可压与不可压光滑粒子动力学方法的封闭方腔自然对流数值模拟及算法对比

为了对比研究弱可压光滑粒子动力学(WCSPH)方法和不可压光滑粒子动 力学(ISPH)方法在模拟封闭方腔自然对流问题时的特性, 采用粒子位移技术有效地解决了高瑞利数条件下, 拉格朗日型SPH方法模拟封闭方腔自然对流时流体域内的粒子聚集和空穴问题, 将拉格朗日型SPH 方法求解封闭方腔自然对流问题的最高瑞利数提高到了10⁶; 进而通过对比瑞利数分别为10⁴, 10⁵, 10⁶的条件下, 采用拉格朗日型WCSPH、 拉格朗日型ISPH、欧拉型ISPH三种SPH方法模拟得到的封闭方腔速度分布云图、 温度分布云图、壁面努赛尔特数分布曲线和平均努塞尔特数, 分析了三种SPH方法在模拟封闭方腔自然对流时的精度、稳定性和计算效率. 结果表明: 在低瑞利数条件下, 以上三种SPH方法都可以较好地模拟此问题, 在高瑞利数条件下, 欧拉型ISPH方法的模拟结果最为精确; 拉格朗日型WCSPH方法模拟所得结果比拉格朗日型ISPH方法模拟所得结果稍好些. 关键词： 光滑粒子动力学 不可压光滑粒子动力学 粒子位移技术 自然对流     

Instabilities developed in stratified flows over pronounced obstacles

In the present work we study numerical and experimentally the flow of a two-layer stratified fluid over a topographic obstacle. The problem reflects a wide number of oceanographic and meteorological situations, where the stratification plays an important role. We identify the different instabilities developed by studying the pycnocline deformation due to a pronounced obstacle. The numerical simulations were made using the model caffa3D.MB which works with a numerical model of Navier-Stokes equations with finite volume elements in curvilinear meshes. The experimental results are contrasted with numerical simulations. Linear stability analysis predictions are checked with particle image velocimetry (PIV) measurements.     

改进的物理粘性SPH方法及其在溃坝问题中的应用

在低雷诺数物理粘性SPH方法基础上引入再生核粒子法进行密度重构,既避免了用人工粘性所导致的数值耗散问题,又提高了低雷诺数物理粘性SPH方法的数值稳定性;以溃坝问题为例,对比分析低雷诺数物理粘性SPH方法和本文方法的仿真结果表明,本文方法可有效消除数值不稳定,压强和速度分布更加光滑,粒子秩序更好,可应用于雷诺数较高或粘性不可忽略的流动问题.     

微液滴振荡过程的光滑粒子动力学方法数值模拟

本文对传统的光滑粒子动力学方法进行了改进, 改进的光滑粒子动力学方法对传统粒子方法中的核梯度进行了修正, 采用了一种新型的核函数和离散格式, 应用改进的光滑粒子动力学方法对微液滴振荡过程进行了数值研究. 研究了不同纵横比和雷诺数(Re)下振荡阻尼与振荡的周期、振幅与Re数的关系. 研究表明: 对于纵横比λ≤ 4时的微液滴振荡过程, 其他参数恒定不变的前提下, Re数越大, 液滴形状变化越剧烈, 波的阻尼作用越弱, 液滴振荡周期变长; 在Re数一定的前提下, 随着液滴初始的纵横比的增大, 液滴振动的振幅增大, 液滴振荡的周期变长.     

激光辐照柱体铝靶的三维光滑粒子流体动力学仿真

在现有光滑粒子流体动力学方法基础上,提出了结合差分方法处理激光辐照边界条件的方法,并利用此方法自编程序模拟计算了功率密度为5.0109 W/cm2的强连续激光辐照铝靶时的热-力学效应,准确给出了与激光作用的边界粒子的分布和飞散情况,同时给出了不同时刻靶体的温度场,得到了温度分布规律,并对结果进行了分析。     

光滑质点流体动力学(SPH)方法(综述)

综合介绍了一种新的纯(Lagrangian方法——光滑质点流体动力学(SPH)方法,由于它计算空间导数时不需要使用网格,从而避免了高维拉氏差分网格法中的网格缠结和扭曲的麻烦,它对非对称的和内含真空区域的三维问题特别有效。重点介绍了该方法的理论基础,流体动力学方程组的推导,人为粘性和热传导,自引力、汇和磁场,光滑核的选取,以及SPH执行过程等有关问题     

二维管道充填过程的修正SPH模拟

通过施加一种密度初始化方法对传统光滑粒子动力学(SPH)方法进行修正,提出一种修正SPH方法.同时,为了提高边界上数值计算的准确性,提出一种新的固壁边界处理方法.通过修正SPH方法模拟液滴拉伸问题和溃坝问题,验证修正SPH方法的准确性和可靠性.随后,对研究很少的管道充填过程进行修正SPH模拟,并讨论Re对流场及涡的影响.数值结果表明,修正SPH方法能够准确模拟牛顿流体管道充填过程,且流动受Re的影响较大.     

SPH simulation of fuel drop impact on heated surfaces

The interaction of liquid drops and heated surfaces is of great importance in many applications. This paper describes a numerical method, based on smoothed particle hydrodynamics (SPH), for simulating n-heptane drop impact on a heated surface. The SPH method uses numerical Lagrangian particles, which obey the laws of fluid dynamics, to describe the fluid flows. By incorporating the Peng–Robinson equation of state, the present SPH method can directly simulate both the liquid and vapor phases and the phase change process between them. The numerical method was validated by two experiments on drop impact on heated surfaces at low impact velocities. The numerical method was then used to predict drop-wall interactions at various temperatures and velocities. The model was able to predict the different outcomes, such as rebound, spread, splash, breakup, and the Leidenfrost phenomenon, consistent with the physical understanding.     

Nonlinear simulations of combustion instabilities with a quasi-1D Navier-Stokes code

As lean premixed combustion systems are more susceptible to combustion instabilities than non-premixed systems, there is an increasing demand for improved numerical design tools that can predict the occurrence of combustion instabilities with high accuracy. The inherent nonlinearities in combustion instabilities can be of crucial importance, and we here propose an approach in which the one-dimensional (1D) Navier-Stokes and scalar transport equations are solved for geometries of variable cross-section. The focus is on attached flames, and for this purpose a new phenomenological model for the unsteady heat release from a flame front is introduced. In the attached flame method (AFM) the heat release occurs over the full length of the flame. The nonlinear code with the use of the AFM approach is validated against analytical results and against an experimental study of thermoacoustic instabilities in oxy-fuel flames by Ditaranto and Hals [Combustion and Flame 146 (2006) 493-512]. The numerical simulations are in accordance with the experimental measurements and the analytical results and both the frequencies and the amplitudes of the resonant acoustic pressure modes are reproduced with good accuracy.     

Simulation of the 3D viscoelastic free surface flow by a parallel corrected particle scheme

In this work,the behavior of the three-dimensional(3D) jet coiling based on the viscoelastic Oldroyd-B model is investigated by a corrected particle scheme,which is named the smoothed particle hydrodynamics with corrected symmetric kernel gradient and shifting particle technique(SPH_CS_SP) method.The accuracy and stability of SPH_CS_SP method is first tested by solving Poiseuille flow and Taylor-Green flow.Then the capacity for the SPH_CS_SP method to solve the viscoelastic fluid is verified by the polymer flow through a periodic array of cylinders.Moreover,the convergence of the SPH_CS_SP method is also investigated.Finally,the proposed method is further applied to the 3D viscoelastic jet coiling problem,and the influences of macroscopic parameters on the jet coiling are discussed.The numerical results show that the SPH_CS_SP method has higher accuracy and better stability than the traditional SPH method and other corrected SPH method,and can improve the tensile instability.     

Nonlinear hydrodynamic phenomena in Stokes flow regime

We investigate nonlinear phenomena in dispersed two-phase systems under creeping-flow conditions. We consider nonlinear evolution of a single deformed drop and collective dynamics of arrays of hydrodynamically coupled particles. To explore physical mechanisms of system instabilities, chaotic drop evolution, and structural transitions in particle arrays we use simple models, such as small-deformation equations and effective-medium theory. We find numerical and analytical solutions of the simplified governing equations. The small-deformation equations for drop dynamics are analyzed using results of dynamical systems theory. Our investigations shed new light on the dynamics of complex fluids, where the nonlinearity often stems from the evolving boundary conditions in Stokes flow.