| 查看: 735 | 回复: 5 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】请教差分得到的代数方程组如何C语言实现
|
|||
|
rt，有没有这样的模板？ 我的问题具体是这样的，如下图一个电势和电流密度的差分后的形式（由电势泊松方程及电流连续性方程推导得到），怎么样写成C描述呢？ 说明：计算区域的网格划分方式，ax、ay分别为x轴和y轴方向上的网格步长，i、j为节点坐标序号。这个网格步长在C程序中如何表示？i和j又如何处理呢？ 谢谢！ [ Last edited by 夏之篮 on 2009-12-29 at 16:56 ] |
回复此楼» 猜你喜欢
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有4人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有7人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有4人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有8人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有5人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有8人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有8人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有5人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有3人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有4人回复
hakekill
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1709.3
- 帖子: 131
- 在线: 19小时
- 虫号: 246017
- 注册: 2006-04-29
- 性别: GG
- 专业: 水环境与生态水利
2楼2009-12-29 16:08:16
3楼2009-12-29 19:48:18
hakekill
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1709.3
- 帖子: 131
- 在线: 19小时
- 虫号: 246017
- 注册: 2006-04-29
- 性别: GG
- 专业: 水环境与生态水利
★ ★ ★ ★ ★ ★ ★ ★
nono2009(金币+2,VIP+0):谢谢应助！ 12-29 20:36
夏之篮(金币+6,VIP+0):谢谢，基本完成我的提问，不好意思，留俩回答我接下来的问题吧，不排除增加悬赏哈~ 12-29 21:24
nono2009(金币+2,VIP+0):谢谢应助！ 12-29 20:36
夏之篮(金币+6,VIP+0):谢谢，基本完成我的提问，不好意思，留俩回答我接下来的问题吧，不排除增加悬赏哈~ 12-29 21:24
|
这个都是通过定义变量来完成的，如网格步长都取0.5 ax=0.5; ay=0.5; i，j分别是你的网格节点数，假如网格节点为4x4的 rows = 4; cols = 4; float u[3,3], rho[3, 3]; max_iter = 10; for(k=1, k for(j = 0, j < cols, j++){ ax2 = ax * ax; ay2 = ay * ay; epsilon = ??; u[i,j]=1/2*(ax2+ay2)*(ax2*(u[i+1,j]+u[i-1,j])+ay2*(u[i,j+1]+u[i,j-1])+ax2*ay2*rho[i,j]/epsilon); .. } } } 大致情况是这样的 [ Last edited by hakekill on 2009-12-29 at 20:07 ] |
4楼2009-12-29 20:05:31
5楼2009-12-29 21:26:24
hakekill
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1709.3
- 帖子: 131
- 在线: 19小时
- 虫号: 246017
- 注册: 2006-04-29
- 性别: GG
- 专业: 水环境与生态水利
★ ★ ★ ★
夏之篮(金币+2,VIP+0):感谢！我自定义两个变量，我就是对如何用宏实现搞不定，找不到你的联系方式，这是我的QQ(357287329)，方便的话，迫切需要你进一步的指导，谢谢！ 12-29 21:54
jjdg(金币+2,VIP+0):辛苦了！ 12-29 22:05
夏之篮(金币+2,VIP+0):感谢！我自定义两个变量，我就是对如何用宏实现搞不定，找不到你的联系方式，这是我的QQ(357287329)，方便的话，迫切需要你进一步的指导，谢谢！ 12-29 21:54
jjdg(金币+2,VIP+0):辛苦了！ 12-29 22:05
|
你如果是用udf的话，就不是我写的程序了，在fluent里面循环和我写的不一样的。 你这个要用fluent自己的宏来搞，比较复杂一点，而且你这个感觉是自己定义了uds，不知道我的理解对不对，不过这个问题已经远超过的能力了 [ Last edited by hakekill on 2009-12-29 at 21:33 ] |
6楼2009-12-29 21:31:42
