应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » C/C++ » 【求助】请教差分得到的代数方程组如何C语言实现
61/1返回列表
查看: 727  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

夏之篮

木虫 (正式写手)

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

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

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

谢谢！

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

» 猜你喜欢
1楼 2009-12-29 13:54:15
已阅   回复此楼   关注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 ]
赞一下(12人) 回复此楼
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 ]
赞一下(12人) 回复此楼
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 ]
赞一下(10人) 回复此楼
6楼2009-12-29 21:31:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 夏之篮 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 3/150 2026-05-17 00:55 by ue3ir18jc3
[论文投稿] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 x0mp7owy2b 2026-05-15 5/250 2026-05-17 00:29 by xiangfeng
[考研] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 x0mp7owy2b 2026-05-15 4/200 2026-05-17 00:25 by ue3ir18jc3
[公派出国] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 v9tggjlwd0 2026-05-15 4/200 2026-05-17 00:15 by ue3ir18jc3
[考研] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 xx7gd5zq4e 2026-05-15 5/250 2026-05-16 23:57 by ue3ir18jc3
[找工作] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 l7k6xnh0yc 2026-05-14 4/200 2026-05-16 23:10 by ue3ir18jc3
[考博] 现在不知道怎么办，感觉很痛苦 +5 qweww 2026-05-11 6/300 2026-05-16 19:57 by Equinoxhua
[基金申请] 请问大佬b0816评完了吗 +4 市民华南虎 2026-05-12 8/400 2026-05-16 19:54 by Equinoxhua
[高分子] 本人最近太闲了，谁有问题可以提，每天会统一回复 +9 一切都是空工 2026-05-12 20/1000 2026-05-16 19:52 by Equinoxhua
[有机交流] 求助2，4-二氯-5-嘧啶甲醛的合成方法 20+3 光吃不拉 2026-05-14 6/300 2026-05-16 19:46 by Equinoxhua
[有机交流] 求有机合成大神指点三硫酸乙烯酯（CAS：2793408-99-6）的合成路线 30+3 Leekmid 2026-05-13 10/500 2026-05-16 16:37 by czyzsu
[公派出国] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 l7k6xnh0yc 2026-05-14 5/250 2026-05-16 16:35 by x28q7dxf75
[有机交流] 如何实现卤原子转化 10+3 BT20230424 2026-05-15 5/250 2026-05-16 16:20 by czyzsu
[硕博家园] 申请博士 +3 呃?呃 2026-05-15 3/150 2026-05-16 11:01 by a4742549
[基金申请] 这年头没有找到涵评专家，还有中面上的可能吗 +11 dd921ww 2026-05-12 13/650 2026-05-16 09:16 by Howard28
[教师之家] 上海大学实验技术岗位非升即走 +5 嘻嘻哈哈乐呵呵 2026-05-15 5/250 2026-05-16 00:17 by caiyun
[文学芳草园] 风把牡丹吹跑了 +5 myrtle 2026-05-12 9/450 2026-05-15 15:27 by myrtle
[基金申请] 青C资助名额大幅增加！ +11 西葫芦炒鸡蛋 2026-05-13 15/750 2026-05-15 14:36 by ambravo
[教师之家] 教学课件你会给同学吗 +8 硕士研究生吗 2026-05-13 8/400 2026-05-14 22:23 by 常规沥青
[考博] 材料类只有一篇综述能申博么 +4 乐逍遥谷 2026-05-13 4/200 2026-05-14 12:05 by zhyzzh
信息提示
关闭
请填处理意见
关闭