24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 696  |  回复: 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的回帖

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的回帖

夏之篮

木虫 (正式写手)

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

引用回帖:
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的回帖

hakekill

木虫 (小有名气)

★ ★ ★ ★
夏之篮(金币+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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 夏之篮 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 本科新能源科学与工程,一志愿华理能动285求调剂 +3 AZMK 2026-03-28 6/300 2026-03-28 17:29 by AZMK
[考研] 张芳铭-中国农业大学-环境工程专硕-298 +4 手机用户 2026-03-26 4/200 2026-03-28 07:17 by mmm just
[考研] 086502化学工程342求调剂 +6 阿姨复古不过 2026-03-27 6/300 2026-03-28 07:06 by wangy0907
[考研] 材料求调剂 一志愿哈工大总分298分,前三科223分 +5 dongfang59 2026-03-27 5/250 2026-03-28 04:53 by wxiongid
[考研] 求调剂 +4 零八# 2026-03-27 4/200 2026-03-27 18:07 by yu221
[考研] 一志愿北京化工大学 070300 学硕 336分 求调剂 +9 vv迷 2026-03-22 9/450 2026-03-27 15:59 by 不吃魚的貓
[考研] 307求调剂 +8 超级伊昂大王 2026-03-24 9/450 2026-03-27 15:34 by 超级伊昂大王
[考研] 一志愿北京化工大学材料与化工 264分各科过A区国家线 +10 哈哈157349 2026-03-21 10/500 2026-03-27 13:06 by zzll406
[考研] 求调剂 +3 刘柯@ 2026-03-24 4/200 2026-03-27 11:28 by shangxh
[考研] 314求调剂 +3 溪云珂 2026-03-26 3/150 2026-03-27 11:20 by sanrepian
[考研] 考研调剂 +9 小蜡新笔 2026-03-26 9/450 2026-03-27 11:10 by 不吃魚的貓
[考研] 333求调剂 +6 wfh030413@ 2026-03-23 6/300 2026-03-26 22:45 by 学员8dgXkO
[考研] 一志愿北化求调剂 +3 Jsman 2026-03-22 3/150 2026-03-26 21:06 by ajpv风雷
[考研] 生物学 296 求调剂 +4 朵朵- 2026-03-26 6/300 2026-03-26 19:01 by 不吃魚的貓
[考研] 07化学303求调剂 +5 睿08 2026-03-25 5/250 2026-03-25 22:46 by 418490947
[考研] 334分 一志愿武理-080500 材料求调剂 +4 李李不服输 2026-03-25 4/200 2026-03-25 21:26 by 星空星月
[考研] 材料调剂 +3 iwinso 2026-03-23 3/150 2026-03-25 11:29 by greychen00
[考研] 求调剂 +6 研研,接电话 2026-03-24 7/350 2026-03-24 17:01 by barlinike
[考研] 材料考研调剂生 +3 黄粱一梦千年 2026-03-24 3/150 2026-03-24 17:00 by barlinike
[考研] 293求调剂 +3 涛涛Wjt 2026-03-22 5/250 2026-03-22 22:21 by jiangpengfei
信息提示
请填处理意见