24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 699  |  回复: 5
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

夏之篮

木虫 (正式写手)

[交流] 【求助】请教差分得到的代数方程组如何C语言实现

rt,有没有这样的模板?
我的问题具体是这样的,如下图一个电势和电流密度的差分后的形式(由电势泊松方程及电流连续性方程推导得到),怎么样写成C描述呢?

说明:计算区域的网格划分方式,ax、ay分别为x轴和y轴方向上的网格步长,i、j为节点坐标序号。这个网格步长在C程序中如何表示?i和j又如何处理呢?

谢谢!

[ Last edited by 夏之篮 on 2009-12-29 at 16:56 ]
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

夏之篮

木虫 (正式写手)

谢谢思路,接着回答我的问题?

引用回帖:
Originally posted by hakekill at 2009-12-29 20:05:
这个都是通过定义变量来完成的,如网格步长都取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< ...

是这样的,我这方程式编好后需要使用UDF加载到FLUENT中的,因为模型的网格划分为非结构网格,x轴和y轴方向上的网格步长都不为定值,我们是不知道的,我觉得应该采用某一个宏去循环遍历线上的所有网格,
此外,这个网格节点(你说的rows和cols)也不知道,也该通过宏读取吧,还有这个max_iter也不知道,那该如何处理?

呵呵,恰逢你也使用FLUENT就好,不知我这个想法,这个实现方式是否可行?

谢谢你的关注,并期待各位专家的高见!
5楼2009-12-29 21:26:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

hakekill

木虫 (小有名气)

★ ★ ★ ★
余泽成(金币+2,VIP+0):谢谢应助! 12-29 19:02
夏之篮(金币+2,VIP+0):谢谢,辛苦了~~ 12-29 19:48
for(k=1,k<=max_iter, k++){
u[i,j]=1/(2*ax*ax+ay*ay)*(ax*ax*(u[i+1,j]+u[i-1,j])+ay*ay*(u[i,j+1]+u[i,j-1])+ax*ax*ay*ay*rho[i,j]/epsilon);
...
}

剩下你自己参考这个写

[ Last edited by hakekill on 2009-12-29 at 16:14 ]
2楼2009-12-29 16:08:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

夏之篮

木虫 (正式写手)

上面还搞不定呢!不是我想要的啊!
我的ax、ay分别为x轴和y轴方向上的网格步长,i、j为节点坐标序号。这个网格步长在C程序中如何表示?i和j又如何处理呢?
3楼2009-12-29 19:48:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hakekill

木虫 (小有名气)

★ ★ ★ ★ ★ ★ ★ ★
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(i = 0, i < rows, i++){
       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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿211 初试270分 求调剂 +8 谷雨上岸 2026-03-23 9/450 2026-03-29 01:00 by 我是小康
[考研] 332求调剂 +4 @MZB382400 2026-03-28 4/200 2026-03-28 21:02 by 唐沐儿
[考研] 343求调剂 +5 爱羁绊 2026-03-28 5/250 2026-03-28 20:53 by 唐沐儿
[考研] 一志愿北化085600材料专硕275|有文章专利|求调剂 +7 Micky11223 2026-03-25 7/350 2026-03-28 18:34 by 无际的草原
[考研] 317求调剂 +6 十闲wx 2026-03-24 6/300 2026-03-28 13:27 by Iveryant
[考研] 求调剂 +6 芦lty 2026-03-25 7/350 2026-03-28 13:13 by 唐沐儿
[考研] 291求调剂 +15 hhhhxn.. 2026-03-23 21/1050 2026-03-28 11:26 by self2008
[考研] 277跪求调剂 +5 1915668 2026-03-27 9/450 2026-03-28 09:58 by zhshch
[考研] 07化学280分求调剂 +10 722865 2026-03-23 10/500 2026-03-27 15:51 by Plutoqq
[考研] 274求调剂 +17 顾九笙要谦虚 2026-03-24 23/1150 2026-03-27 15:16 by caszguilin
[考研] 08开头275求调剂 +4 拉谁不重要 2026-03-26 4/200 2026-03-27 14:12 by Delta2012
[考研] 一志愿华东理工大学081700,初试分数271 +6 kotoko_ik 2026-03-23 7/350 2026-03-27 12:29 by 惠州彭于晏
[考研] 316求调剂 +5 Pigcasso 2026-03-24 5/250 2026-03-27 12:10 by zhshch
[考研] 调剂求收留 +7 果然有我 2026-03-26 7/350 2026-03-27 00:26 by wxiongid
[考研] 材料调剂 5+4 想要一壶桃花水 2026-03-25 10/500 2026-03-26 19:56 by 不吃魚的貓
[考研] 334分 一志愿武理-080500 材料求调剂 +4 李李不服输 2026-03-25 4/200 2026-03-25 21:26 by 星空星月
[考研] 290分调剂求助 +3 吉祥止止陈 2026-03-25 3/150 2026-03-25 19:58 by barlinike
[考研] 0854人工智能方向招收调剂 +4 章小鱼567 2026-03-24 4/200 2026-03-25 13:29 by 2177681040
[考研] 材料/农业专业,07/08开头均可,过线就行 +3 呵唔哦豁 2026-03-23 4/200 2026-03-23 22:30 by 汪!?!
[考研] 接收2026硕士调剂(学硕+专硕) +4 allen-yin 2026-03-23 6/300 2026-03-23 15:04 by 汪!?!
信息提示
请填处理意见