- 积分
- 36
- 注册时间
- 2004-2-6
- 仿真币
-
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2008-10-25 16:54:22
|
显示全部楼层
来自 大连理工大学
这里在fepg公司的论坛上面找到了两个非线性抛物方程耦合求解的算例,大家容易做错,参考参考。
一个是关于压强的方程,一个是关于温度的方程.
大概意思就是 一端注入甲烷,但注入的甲烷会发生反应,生成热,随着温度升高,
管内压强会变大.
FEPG有限元文件见附件里面的proto文件.虽然是一维的程序,难点并不在维数.
计算程序见附件.就是收敛比较慢.
需要注意的是第二个场的nfe文件,比较一下和普通的nparb的区别.
主要是read,和write语句. read和write的最后一个vect和本问题相关.可以不用管.
注意:
两个非线性抛物方程耦合求解第二个场的nfe文件不能直接用nparb.
因为在第一个场的非线性循环内,第二个场的初值每次都得用u1.
而如果直接用nparb,则会变成,在第一个场的非线性循环内,第二个场已经在进行以后几个时间步的计算了. |
本帖子中包含更多资源
您需要 登录 才可以下载或查看,没有账号?注册
×
|