LAMMPS30Sep2026：SPH密度核、Tait压力与人工黏性

解压进入lammps-sph，在分配好的CPU作业内bash run.sh；也可以sbatch job.slurm。cpu-8173m按56物理核整节点分配和计费，本例仅使用1线程，CPU oneAPI2024.2构建，不声称使用GPU加速。分析绘图sisl0.16.4。

四个实际run0算例为density、pressure、approach、recede。4个SPH粗粒子每个0.01kg，坐标为米，核半径0.03m。density实际调用sph/rhosum，包含自项与邻居项；另外三组固定rho1100kg/m3，用sph/taitwater，rho0=1000kg/m3、c0=10m/s、alpha=1。approach为v=-0.2*(r-0.05) m/s，recede相反；只有相向运动的原子对触发本人工黏性项。

这是静态密度和受力教学，不是4个水分子，也没有进行MD积分、流动、温度平衡或黏度测量。run0只建立邻居并计算当前受力；尽管输入列出fix sph和时间步，并没有推进物理时间。density的848.68–914.83kg/m3来自四点稀疏核求和，不能据此推断实际水密度偏低。压力三组的1100为人为指定并明确分开。

本版本sph/rhosum使用三维quadric核3465/(512pi*h^3)*(1-r^2/h^2)^4；sph/taitwater压力力使用Lucy核导数，两者不能混成同一个势能梯度。Tait p=rho0*c0^2/7*((rho/rho0)^7-1)。独立程序逐点核对4个密度、12个原子受力分量、总力零和相向/背向黏性对照。没有把sph/idealgas、SPH全部风格、DPD或peridynamics混称为已完成。

raw含真实LAMMPS输入和输出，results/particle-fields.csv为全粒子数据，results/full-fields.npz保留完整数组，plots为本次新图。reference-plots仅为原示例参考图。输入可读，但本包数值检查限定这4个固定教学条件；改结构、质量、单位或核半径时另建目录并重新设定物理检查。
