| 查看: 2539 | 回复: 5 | |||
xiaoyueriyue金虫 (正式写手)
|
[求助]
相场法步骤函数处理
|
|
我最近在做多相场法的共晶凝固的模拟，有一个小问题希望能够得到大家的解答，在多相场法的模拟中，有一个步骤函数（step function) 的定义si(x,t)=1,if phi(x,y,i)>0, otherwise si(x,t) =0;在程序编写的过程中，根据参考文献的意见（Phase-field modeling of eutectic solidification, Kim. S. G, Kim. W. T., et at. Journal of Crystal Growth， 2004)的，modifying the definition as si(x,t) = 1 if the sum of phi on the four nearest neighbor grids is between alpha and 4-alpha, 文献中alpha的参考值是在0.02，但是如果这样定义，会存在一个问题：在某一相中，如果某一个节点四个临近节点都是同样的固相，那么所得到的的si(x,t)的值成为了0，与前面一个定义是并不相符的。 这是我最近编程过程中遇到的困惑，希望能够有人帮忙解答。不胜感激。 |
回复此楼» 猜你喜欢
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有4人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有3人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有5人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有5人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有5人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有6人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有6人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有6人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有7人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 有人理解tvd算法和壁面函数法的求解步骤么
已经有4人回复
- 函数 f(x)=x*tan(x) 能够积分写出原函数吗 最好有过程
已经有7人回复
- 求助相场动力学模拟用什么软件比较好？主要用相场法
已经有5人回复
- 【讨论】预调件共轭梯度法（PCG）
已经有50人回复
xiaohu3956
铁杆木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 5568.7
- 散金: 66
- 红花: 2
- 帖子: 341
- 在线: 207.4小时
- 虫号: 688730
- 注册: 2009-01-06
- 性别: MM
- 专业: 自然语言理解与机器翻译
2楼2013-11-25 09:40:48
xiaoyueriyue
金虫 (正式写手)
- 应助: 6 (幼儿园)
- 金币: 964.4
- 散金: 1051
- 红花: 6
- 帖子: 381
- 在线: 128小时
- 虫号: 1966386
- 注册: 2012-08-31
- 性别: GG
- 专业: 微/纳机械系统
★ ★ ★
xiegangmai: 金币+3, 鼓励讨论交流 2013-11-30 21:56:49
xiegangmai: 金币+3, 鼓励讨论交流 2013-11-30 21:56:49
|
这个贴我发布之后，我与文章的作者取得了联系，他们的解释可以给大家参考： By puttimg sum_phi = phi(i-1,j)+phi(i+1,j)+phi(i,j-1)+phi(i,j+1) we then define s(i,j)=1 if sum_phi > 0 and s(i,j)=0 if sum_phi =0 It is a useful ingredient in simulation. With the strict definition of s(i,j)=0 for phi(i,j)=0 and s(i,j)=1 for phi(i,j)>0, there can not be any propagation of phase-field in simulation. |
3楼2013-11-28 13:55:04
xiaohu3956
铁杆木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 5568.7
- 散金: 66
- 红花: 2
- 帖子: 341
- 在线: 207.4小时
- 虫号: 688730
- 注册: 2009-01-06
- 性别: MM
- 专业: 自然语言理解与机器翻译
4楼2013-11-30 21:49:57
xiaoyueriyue
金虫 (正式写手)
- 应助: 6 (幼儿园)
- 金币: 964.4
- 散金: 1051
- 红花: 6
- 帖子: 381
- 在线: 128小时
- 虫号: 1966386
- 注册: 2012-08-31
- 性别: GG
- 专业: 微/纳机械系统
5楼2013-12-01 13:31:06
xiaohu3956
铁杆木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 5568.7
- 散金: 66
- 红花: 2
- 帖子: 341
- 在线: 207.4小时
- 虫号: 688730
- 注册: 2009-01-06
- 性别: MM
- 专业: 自然语言理解与机器翻译
6楼2013-12-01 22:00:20
10