参考《不重新编译也能做LAMMPS二次开发吗——写插件就可以了》。 将这个pair_style命名为 sph/rhosum/norm。然后在脚本中按如下方式调用: pair_style sph/rhosum/norm Nevery pair_coeff I J h Nevery 是每多少步做一次密度重整。 改进效果 算例为圆管中的溃坝。跟昨天《SPH也要做类似“网格无关性“的验...
LAMMPS代码 首先定义维度、单位和原子类型。这是个二维算例。单位采用SI单位制。SPH中,atom_style 是sph,不是以前的 meso 了。 dimension 2 units si atom_style sph 定义区域。首先定义了一个模拟盒子。创建三类原子,分别是流体(1号)、壁面(2号)和顶盖(3号)。 注意这边先创建了晶格,然后再创建原子。如此,...
LAMMPS软件的SPH模拟教程更新滞后问题引起关注。LAMMPS-SPH中某些关键字已更名，如"meso"改为了"sph"，"meso_e"改为了"meso/e"，因此在使用教程时必须注意最新文档与示例代码。以下为顶盖驱动方腔流模拟的具体步骤。定义模拟参数，包括维度、单位和原子类型，这里使用二维模拟，单位采用国际单位制（SI）。
git-svn-id: svn://svn.icms.temple.edu/lammps-ro/trunk@6493 f3b2605a-c… Jul 1, 2011 src fix temperature Feb 14, 2012 tools add cv in meso Nov 25, 2011 LICENSE git-svn-id: svn://svn.icms.temple.edu/lammps-ro/trunk@20 f3b2605a-c51… ...
U. (2011) The implementation of Smooth Particle Hydrodynamics in LAMMPS, Retrieved from lammps.sandia.gov/doc/PDF/SPH_LAMMPS_userguide.pdf (last accessed 03/12/16).Georg, C., "The Implementation of Smooth Particle Hydrodynamics in LAMMPS", Ph.D. thesis, University of Leuven, (2011)...
FixSPH::init() 很简单,初始化了两种 \delta t。dtv 用于根据速度更新位置,因为位置是从 t 直接更新到 t+\delta t,所以 dtv 是一个完整的时间步。dtf 用于根据力更新其他变量,因为基本上都是用 \delta t/2更新,所以乘了0.5。有的用的是 \delta t/(2m),这个是 dtfm,后面会创建。 force->ftm2v 用于...
∇iW(rij,h)是SPH核函数的对颗粒i的梯度。LAMMPS在植入这个方程时,核函数采用的是Lucy's Quartic Kernel: W(q,h)=αD(1+3q)(1−q)3q={‖r‖/hr<h0r≥h, 其中αD在二维问题中取5/(πh2),在三维问题中取105/(16πh3)。 现在来求核函数的梯度。