| 查看: 6646 | 回复: 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加剧腐蚀的预期相反。 我想请问，对于这种分子动力学模拟装置而言，这个结果在物理上是否合理？或者是否存在其他问题？  |
回复此楼» 猜你喜欢
应该是93bebmhtak前后十一个字符比较关键
已经有10人回复
-
这样的filecode谁见过
已经有21人回复
-
filecode与中标关系的预测
已经有5人回复
-
关于代码变化问题，想知道的进来
已经有17人回复
-
fileCode有新解读？
已经有18人回复
-
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
-
2026国自然放榜时间
已经有3人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有7人回复
» 本主题相关商家推荐: (我也要在这里推广)
HER12025
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 122.3
- 散金: 35
- 帖子: 193
- 在线: 102.6小时
- 虫号: 36722173
- 注册: 2025-12-26
- 专业: 高分子材料结构与性能
2楼2026-05-13 18:24:08
100