| 查看: 636 | 回复: 1 | ||
[求助]
急求UDF中用PR方程求解超临界流体co2密度的UDF程序
|
|
我在编写用PR方程求解超临界co2密度的UDF程序时采用的是盛金公式，但我只考虑了盛金公式中得Δ=B^2－4AC>0时的情况，所以在初始化流场时通不过，求哪位大侠可以帮我分析分析，急求有类似的程序，或也写过此程序的能够帮助我一下，谢谢了。或者可以提供一个C语言求三次方程根的程序，其中a,b,c,d都是变量的，附上我写的程序： #include "udf.h" #define TC 305.25 #define PC 7380000 #define W 0.225 #define R 8.314 #define tc 305.25 #define pc 7375000 #define w 0.225 #define r 8.314 DEFINE_PROPERTY(cell_density, cell, thread) { real dd; real temp = C_T(cell, thread); real pres = C_P(cell, thread); real tr,k,a1,a2,a,b,a3,a4,a5,a6,b1,b2,b3,y1,y2,y3,x; tr = temp/TC; k = 0.37464+1.54226*W-0.26992*W*W; a1 = (1+k*(1-sqrt(tr)))*(1+k*(1-sqrt(tr))); a2 = 0.45724*R*R*TC*TC/PC; a = a1*a2; b=0.0778*R*TC/PC; a3=pres; a4=pres*b-R*temp; a5=a-3*pres*b*b-2*R*temp*b; a6=pres*b*b*b+R*temp*b*b-a*b; b1 = a4*a4-3*a3*a5; b2 = a4*a5-9*a3*a6; b3 = a5*a5-3*a4*a6; y1 = b1*a4+3*a3*(-b2+sqrt(b2*b2-4*b1*b3))/2; y2 = b1*a4+3*a3*(-b2-sqrt(b2*b2-4*b1*b3))/2; if (y1<0) {y3=-y1;} else {y3=y1;} x = (-a4-(pow(y3,0.33333)-pow(-y2,0.33333)))/(3*a3); dd =44/(1000*x); return dd; } UID99953 帖子5 精华0 积分10 阅读权限10 来自湖南 在线时间11 小时 注册时间2010-8-31 最后登录2012-7-18 查看详细资料 |
回复此楼» 猜你喜欢
实验室接单子
已经有4人回复
-
全日制（定向）博士
已经有4人回复
-
假如你的研究生提出不合理要求
已经有6人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
不自信的我
已经有12人回复
-
所感
已经有4人回复
-
论文终于录用啦！满足毕业条件了
已经有28人回复
-
要不要辞职读博？
已经有7人回复
-
北核录用
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 等密度下超临界CO2对聚合物的溶解性
已经有10人回复
- 好心人帮帮忙吧，为什么这个UDF导入fluent以后，在边界条件处没有显示~~~
已经有17人回复
- fluent中用udf编写三维速度边界条件出现错误，大家来指导一下啊
已经有9人回复
- 多相流中用udf定义混合物的材料属性
已经有10人回复
- 这个求解浓度方程的udf怎么不对呢？急啊急啊~SOS
已经有7人回复
- 能量方程源项udf
已经有4人回复
- 曳力函数的的UDF
已经有15人回复
- udf profile update interval 是什么意思？
已经有5人回复
- udf 重力方向变化
已经有17人回复
- 【求助】udf读取混合物密度失败、及平均流速读取？
已经有8人回复
- 【求助】单机多核fluent并行udf
已经有13人回复
- 【求助】Warning: incorrect cg motion UDF value on zone 35
已经有5人回复
- 【求助】fluent计算中，udf编译通过，初始化出错
已经有11人回复
- 【求助】物料属性可以用UDF修改吗
已经有8人回复
Edy159357
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 67.4
- 帖子: 29
- 在线: 21.5小时
- 虫号: 10675975
- 注册: 2018-10-14
- 专业: 工程热物理相关交叉领域
2楼2019-04-23 21:48:52