科研科学计算高性能计算【免费下载链接】lammpsPublic development project of the LAMMPS MD software package项目地址https://gitcode.com/gh_mirrors/la/lammps点击查看免费下载导读fix viscous/sphere是 LAMMPS 中专门面向有限尺寸球体颗粒finite-size spheres的粘性阻尼力矩命令它为颗粒施加与角速度成正比的反向阻尼力矩用于可控地消耗转动动能。本文围绕该命令的语法、物理原理、scale缩放机制、单位注意事项、加速器变体与限制展开并结合仓库源码EXTRA-FIX 包与单元测试给出可复现的实战配置读完即可在自己的颗粒/粒状模拟中正确使用并调参。一、命令总览与定位fix viscous/sphere属于 EXTRA-FIX 包下文将详细说明安装要求其作用与平移版本fix viscous相对应fix viscous对原子的平动速度施加阻尼力F -γv而fix viscous/sphere对原子的角速度施加阻尼力矩T -γω只作用于组group内携带omega角速度与torque力矩属性的有限尺寸球体原子。在粒状granular模拟中该命令常用于在受控方式下排干系统的转动动能若不加额外热浴它会缓慢或快速地使系统冻结因此也可作为一种简易的能量最小化手段。其底层实现位于 src/EXTRA-FIX/fix_viscous_sphere.cpp类名为FixViscousSphere通过FixStyle(viscous/sphere, FixViscousSphere)注册见 src/EXTRA-FIX/fix_viscous_sphere.h 第 16 行。二、语法与参数详解fix ID group-ID viscous/sphere gamma keyword values ...参数说明ID、group-ID与其他 fix 命令一致参见 fix 命令文档viscous/sphere本命令的 style 名称gamma阻尼系数单位为力矩/角速度torque/angular velocity 单位keyword values可选的零个或多个关键字/值对目前仅支持scalescale 关键字的两种用法keyword scale scale values type ratio 或 v_name type 原子类型1-N ratio 阻尼系数缩放因子 v_name 引用原子风格变量atom-style variable名称按类型缩放scale type ratio可多次使用分别对不同原子类型设置不同的缩放因子。例如fix 1 damp viscous/sphere 0.5 scale 3 2.5表示将原子类型 3 的阻尼系数放大为0.5 × 2.5。按原子变量缩放scale v_name引用一个 atom-style 变量每个原子可独立得到例如基于自身半径计算的缩放因子。变量名以v_前缀传入实际变量名去掉该前缀。官方示例fix 1 flow viscous/sphere 0.1 fix 1 damp viscous/sphere 0.5 scale 3 2.5 fix 1 damp viscous/sphere 0.5 scale v_radscale参数解析的源码证据在 fix_viscous_sphere.cpp 的构造函数中强制校验atom-omega_flag第 41 行若原子不含角速度属性则直接报错Fix viscous/sphere requires atom attribute omegagamma通过utils::numeric解析第 45 行解析scale时区分两种模式参数以^v_开头则走变量模式第 53-64 行否则走按类型模式第 66-78 行。两种模式互斥混用会报错Must use only one style of scaling按类型模式下内部为每个原子类型维护一个scalegamma数组默认值为1.0第 68-71 行并对原子类型越界做检查第 74-75 行变量模式下要求所引用的变量确实是 atom-style 变量第 59-60 行并在init()阶段再次确认第 120-128 行。三、物理原理阻尼力矩如何施加阻尼力矩公式为T_i -γ · ω_i系数越大转动动能衰减越快。注意这里用户直接指定 γ 本身力矩/角速度单位而非指定黏度若需要按颗粒尺寸σ调整不同大小原子的阻尼使用scale关键字即可源码注释也明确写到direction is opposed to angular velocity vector方向与角速度向量相反见 fix_viscous_sphere.cpp。底层实现调用链FixViscousSphere的setmask()注册了POST_FORCE、POST_FORCE_RESPA与MIN_POST_FORCE三个回调fix_viscous_sphere.cpp意味着它分别在常规力计算之后、r-RESPA 分层积分器对应层级以及能量最小化迭代中改写力矩。核心的post_force()实现fix_viscous_sphere.cpp逻辑清晰取本进程管辖的局部原子数组omega、torque、mask、type若为变量缩放模式先调用input-variable-compute_atom(...)计算每个原子的缩放值存入scaleval第 165-168 行遍历nlocal个原子凡mask[i] groupbit命中的原子根据三种缩放模式计算有效阻尼系数dragNONEdrag gammaTYPEdrag gamma * scalegamma[type[i]]VARIABLEdrag gamma * scaleval[i]执行力矩更新torque[i][k] - drag * omega[i][k]k x, y, z。从源码可以推断该 fix 只修改力矩数组不直接修改角速度——真正把力矩积分成角速度变化的是下游的/sphere系列积分器如fix nve/sphere这解释了为何文档要求配合atom_style sphere与带/sphere后缀的积分器使用。四、单位注意事项易踩坑点文档中的警告非常关键γ 的单位是力矩/角速度并非质量/时间至少在 LAMMPS 的某些非自洽单位制如real、metal下两者并不相等。不同单位制下的力矩定义参见 units 命令文档单位制力矩单位ljε约化单位realkcal/molmetaleVsiNewton-meterscgsdyne-centimetersnano皮克-微米²/微秒²micro阿托克-纳米²/纳秒²因此设置 γ 时必须基于当前units下的力矩与角速度量纲来估算而不是直接照搬质量/时间单位的数值。若颗粒尺寸差异较大优先借助scale按类型或按变量对 γ 做等效折算以体现不同直径颗粒阻尼不同的物理效果。五、与 fix viscous平移阻尼的对比与联动fix viscous/sphere的姊妹命令是 fix viscous二者物理与语法结构几乎平行维度fix viscousfix viscous/sphere阻尼对象平动速度 v角速度 ω公式F -γvT -γωγ 单位力/速度力矩/角速度适用原子任意原子含点粒子仅有限尺寸球体所属包核心无包限制EXTRA-FIX 包在真实工程中两者经常同时启用。例如仓库自带示例 examples/bpm/metamaterial/in.bpm.metamaterial磁驱动杆状力学超材料模拟中就同时写了fix 3 all viscous 1e-3 fix 4 all viscous/sphere 1e-3即同一系统同时用fix viscous消耗平动动能、用fix viscous/sphere消耗转动动能使体系快速弛豫到低能构型。另外fix viscous/sphere文档列出的相关命令还包括 fix damping/cundall实现位于 src/GRANULAR/fix_damping_cundall.cpp它是基于 Cundall 非黏性阻尼公式的另一种旋转阻尼方案适合需要与颗粒接触力解耦的粒状模拟可按需对比选用。六、集成器交互r-RESPA、能量最小化与输出该 fix 与运行框架的交互规则与 fix viscous 完全一致如下重启文件不向二进制重启文件写任何信息fix_modify所有常规选项均不适用唯一例外是respa选项——它支持设置 r-RESPA 积分器run_style respa中修改力矩的层级默认在最外层最粗层级输出不存储任何全局量或每原子量compute/dump等输出命令无法读取run start/stop其参数不可与run的 start/stop 关键字搭配。源码中对 r-RESPA 的支持体现在构造函数设置respa_level_support 1fix_viscous_sphere.cppinit()中检测到respa积分风格后计算默认层级为最外层第 115-118 行post_force_respa()仅在目标层级触发实际阻尼第 187-190 行。能量最小化该 fix 施加的力矩在minimize过程中同样生效但只能配合允许非保守力的阻尼动力学类最小化器damped dynamics minimizers即 min_style quickmin 或 fire例如min_style fire minimize 0.0 1.0e-8 1000 10000这与阻尼耗散 无热浴 ≈ 冻结系统的机制一致quickmin/fire 本身就是带阻尼的非保守动力学方法允许这类非保守力矩参与迭代。普通共轭梯度cg、最速下降sd等保守方法不适合搭配。七、加速器变体viscous/sphere/kkKOKKOS该命令提供 KOKKOS 加速变体viscous/sphere/kk对应 src/KOKKOS/fix_viscous_sphere_kokkos.cpp可在 GPU/OpenMP 等后端并行执行适合大规模球体体系。其核心post_force通过Kokkos::parallel_for对局部原子做核外循环并将有效阻尼系数预先烘焙为k_eff_gamma双视图数组以减少设备端计算fix_viscous_sphere_kokkos.cpp。从源码结构看/kk变体目前存在两个明确限制均在init()中直接报错不支持 r-RESPACannot (yet) use respa with fix viscous/sphere/kk不支持变量缩放Cannot (yet) use variable style with fix viscous/sphere/kk仅支持scale type ratio或纯 γ。因此若使用 KOKKOS 加速请避免run_style respa与scale v_xxx写法改用按类型缩放。八、前置条件与限制Restrictions所属包属于 EXTRA-FIX 包编译时必须启用该包例如传统构建make yes-extra-fix或 CMake 构建-D PKG_EXTRA-FIXon详见 Build package 文档原子属性要求原子存储torque、omega角速度与半径即必须使用 atom_style sphere或atom_style hybrid sphere ...定义的原子粒子形态组内所有粒子必须是有限尺寸球体直径 0不能是点粒子直径 0上述任一条件不满足运行时即报错如缺omega属性时报Fix viscous/sphere requires atom attribute omega。九、单元测试验证仓库在 unittest/force-styles/tests/fix-timestep-viscous_sphere.yaml 中提供了该 fix 的数值回归测试前置条件为atom spherefix viscous/spherefix nve/sphere后置命令为fix move all nve/sphere fix test solute viscous/sphere 0.05 scale 1 2.0即对 32 个原子、以timestep 5e-04、epsilon 容差5e-14校验单步力矩与位置/速度/力矩输出的精度验证了scale type ratio路径的正确性。若需在本地复现可参考 unittest 目录 中 force-styles 测试的运行方式。十、最小可用配置示例结合上文给出一个可直接套用的最小 LAMMPS 输入片段球体 旋转阻尼 冻结式最小化units lj atom_style sphere region box block 0 10 0 10 0 10 create_box 1 box create_atoms 1 single 5 5 5 set group all diameter 1.0 set group all density 1.0 velocity group all set 0.0 0.0 0.0 omega 0.5 0.3 -0.2 units box pair_style none neighbor 0.3 bin # 平动 转动双阻尼消耗动能 fix 1 all nve/sphere fix 2 all viscous 0.1 fix 3 all viscous/sphere 0.1 scale 1 1.0 run 1000 # 或作为简易能量最小化改用阻尼动力学最小化器 min_style fire minimize 0.0 1.0e-8 1000 10000要点回顾atom_style sphere提供omega/torque/radius属性fix nve/sphere负责把阻尼力矩积分成角速度变化γ 需按当前units的力矩/角速度量纲选取需要按颗粒尺寸差异化阻尼时使用scale按类型或按 atom-style 变量。十一、小结核心能力fix viscous/sphere以T -γω对有限尺寸球体施加角速度阻尼是粒状模拟中耗散转动动能、冻结系统或作为简易能量最小化手段的标配工具关键参数γ 的单位是力矩/角速度scale type ratio与scale v_name提供按类型/按原子两种差异化缩放源码要点实现位于 src/EXTRA-FIX/fix_viscous_sphere.cpp通过post_force系列回调改写力矩数组KOKKOS 变体fix_viscous_sphere_kokkos.cpp不支持 respa 与变量缩放注意事项需 EXTRA-FIX 包、atom_style sphere、有限尺寸球体能量最小化仅适用于 quickmin/fire 等阻尼动力学方法。赞分享科研科学计算高性能计算【免费下载链接】lammpsPublic development project of the LAMMPS MD software package项目地址https://gitcode.com/gh_mirrors/la/lammps点击查看免费下载相关推荐LAMMPS fix npt/sphere 命令详解有限尺寸球形粒子的 Nose–Hoover NPT 积分LAMMPS fix npt/sphere 命令详解有限尺寸球形粒子的 Nose–Hoover NPT 积分 fix npt/sphere 是 LAMMPS科研科学计算高性能计算LAMMPS fix viscous 命令深度解析速度比例阻尼、Stokes 拖拽与布朗动力学的实现与实战LAMMPS fix viscous 命令深度解析速度比例阻尼、Stokes 拖拽与布朗动力学的实现与实战 LAMMPS 的 fix viscous 为组内原科研科学计算高性能计算LAMMPS 中 fix nvt/sphere 命令详解有限尺寸球形粒子的 Nose/Hoover 恒 NVT 系综控温LAMMPS 中 fix nvt/sphere 命令详解有限尺寸球形粒子的 Nose/Hoover 恒 NVT 系综控温 本文围绕 LAMMPS 的 fix科研科学计算高性能计算上一篇OpenVR-SpaceCalibrator安全与兼容性确保VR设备稳定运行的最佳实践下一篇gh_mirrors/tr/training-data-analyst项目中的药物发现与分子设计创作声明：本文部分内容由AI辅助生成（AIGC），仅供参考