| 查看: 8072 | 回复: 1 | ||
| 【悬赏金币】回答本帖问题，作者wzl252825将赠送您 15 个金币 | ||
[求助]
求助：模拟流动加速腐蚀，模拟结果与预期结果相反
|
||
|
求助：模拟流动加速腐蚀，模拟结果与预期结果相反 LAMMPS模拟Fe / 液态 LBE（Pb-Bi 共晶体）界面处的腐蚀/互扩散，并比较静态 LBE和流动 LBE条件下的行为。 模型设置： 下部：固态铁； 上部：液态铅铋合金（Pb/Bi）； 界面法线方向：z； 液态铅铋合金上方留有真空层 ；原子间势：EAM/alloy； 边界条件：P P P； 底部部分铁原子固定，其余原子定义为mobile模型如图： 在流动条件下，我向液态LBE施加一个x方向的体积力（平行于界面）以驱动流动。在生产阶段，我使用单个mobile nvt系综，并temp/profile用它来消除x方向的流动偏差： compute Tflow mobile temp/profile 1 0 0 z 20 region LBE_flow block INF INF INF INF 58.86 124 units box variable fx equal 5.0e-5 fix 3 LBE addforce v_fx 0.0 0.0 region LBE_flow fix 4 mobile nvt temp 973 973 0.1 fix_modify 4 temp Tflow 分析的量包括： 1、渗入Fe中的Pb/Bi原子数； 2、Pb/Bi渗入Fe的最大深度； 3、溶解到LBE中的Fe原子数。 与静态LBE情况相比，在流动条件下，我得到了与预期相反的结果： 1、渗透的Pb/Bi原子数量减少； 2、最大渗透深度变小； 3、溶解的Fe原子数量也减少。 在我目前的模拟中，流动状态并没有加速腐蚀/渗透。相反，它似乎抑制了渗透次数和渗透深度，这与通常认为流动的LBE会比静态LBE加剧腐蚀的预期相反。 我想请问，对于这种分子动力学模拟装置而言，这个结果在物理上是否合理？或者是否存在其他问题？  |
回复此楼» 猜你喜欢
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有21人回复
-
能否退出参与的面上项目解除限项
已经有17人回复
-
估计是周四
已经有4人回复
-
我面上完蛋了
已经有9人回复
-
2026国自然函评费到账
已经有18人回复
-
建议基金发布提前给出明确的时间点
已经有15人回复
-
范进中举一文的中心思想
已经有7人回复
-
朋友圈看到的
已经有9人回复
-
让我中一个面上吧！
已经有16人回复
-
放榜前的不淡定
已经有17人回复
» 本主题相关商家推荐: (我也要在这里推广)
HER12025
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 135.8
- 散金: 35
- 帖子: 206
- 在线: 112.4小时
- 虫号: 36722173
- 注册: 2025-12-26
- 专业: 高分子材料结构与性能
2楼2026-05-13 18:24:08
20