有关超临界流的仿真 水 物性参数什么的

cfkristen

都在小木虫里面了 就不细说了
http://emuch.net/bbs/viewthread.php?tid=5939030&pid=1#pid1
最近在做超临界的东西 参考一篇文献 选了一个水的先试试 就是物性参数（比热 密度 黏度 导热系数）用公式取的点用Origin拟合成多项式用UDF输入进去的 和物性参数曲线实际曲线对比了下基本没问题 相关的设置全部依照文献做的  就是一个圆形管道 算的是流固耦合 进口流量进口 出口压力出口 出口一直有回流 一直不能收敛

1.UDF代码：

/***********************************************************************/
/* viscosity.c */
/* UDF for specifying a temperature-dependent viscosity property */
/***********************************************************************/
#include "udf.h"
#include "math.h"
DEFINE_PROPERTY(cell_viscosity, cell, thread)
{
real mu;
real temp;
temp= C_T(cell, thread);

if (temp >=540 && temp<670)
mu=-6.26542E-13*pow(temp,4)+1.37963E-9*pow(temp,3)-1.12934E-6*pow(temp,2)+4.0331E-4*temp-0.05139;
else if (temp >=670 && temp<=720)
mu=-8.9484E-11*pow(temp,4)+2.49246E-7*pow(temp,3)-2.60191E-4*pow(temp,2)+0.12064*temp-20.96297;
else
        mu=3.05888E-14*pow(temp,4)-1.11035E-10*pow(temp,3)+1.50831E-7*pow(temp,2)-9.0531E-5*temp+0.02054;
return mu;
}

DEFINE_PROPERTY(cell_ktc, cell, thread)
{
real ktc;
real temp;
temp= C_T(cell, thread);

if (temp >=540 && temp<650)
ktc=-2.27855E-10*pow(temp,4)+5.12332E-7*pow(temp,3)-4.37236E-4*pow(temp,2)+0.16643*temp-23.04113;
else if (temp >=650 && temp<700)
ktc=5.60049E-8*pow(temp,4)-1.51738E-4*pow(temp,3)+0.15408*pow(temp,2)-69.50274*temp+11751.89412;
else if (temp >=700 && temp<850)
ktc=1.01777E-9*pow(temp,4)-3.24317E-6*pow(temp,3)+0.00387*pow(temp,2)-2.05748*temp+409.82069;
else
        ktc=3.53259E-12*pow(temp,4)-1.42586E-8*pow(temp,3)+2.17375E-5*pow(temp,2)-0.01474*temp+3.85046;

return ktc;
}

DEFINE_PROPERTY(cell_density, cell, thread)
{
real rho;
real temp;
temp= C_T(cell, thread);

if (temp >=540 && temp<650)
rho=-5.23274E-7*pow(temp,4)+0.00118*pow(temp,3)-1.00937*pow(temp,2)+382.40622*temp-53339.08108;
else if (temp >=650 && temp<700)
rho=8.71409E-5*pow(temp,4)-0.23348*pow(temp,3)+234.48023*pow(temp,2)-104613.64289*temp+1.74966E7;
else
        rho=9.81174E-8*pow(temp,4)-3.50877E-4*pow(temp,3)+0.4696*pow(temp,2)-278.99576*temp+62232.92532;

return rho;
}

DEFINE_SPECIFIC_HEAT (cp, T, Tref, h, yi)
{
  real cp;

if (T >=540 && T<650)
cp=4.12126E-5*pow(T,4)-0.0947*pow(T,3)+81.65593*pow(T,2)-31302.17859*T+4.5041E6;
else if (T >=650 && T<688)
cp=-0.02116*pow(T,4)+56.4407*pow(T,3)-56449.66828*pow(T,2)+2.50871E7*T-4.18E9;
else if (T >=688 && T<750)
cp=-0.00139*pow(T,4)+3.98394*pow(T,3)-4266.71572*pow(T,2)+2.0283E6*T-3.61059E8;
else
        cp=1.97284E-6*pow(T,4)-0.0073*pow(T,3)+10.13277*pow(T,2)-6260.7141*T+1.45669E6;
return cp;
}

看了网上的很多帖子 基本问题也不大

2.结果
相关结果如图
3.边界条件
如图

刚看了下 好像水力直径给错了 估计这个问题不大的吧 不影响结果的

[ 本帖最后由 cfkristen 于 2013-5-26 11:41 编辑 ]

通流

流体只有在临界点附近的物性比较难搞。只要远离临界点，不管是超临界，还是亚临界，都还是比较容易处理的。

我只能写这些。那个木虫，我进不了。

cfkristen

原帖由 通流 于 2013-5-26 08:47 发表

                               
登录/注册后可看大图

流体只有在临界点附近的物性比较难搞。只要远离临界点，不管是超临界，还是亚临界，都还是比较容易处理的。

我只能写这些。那个木虫，我进不了。

通版有时间再看看吧 超临界的确不难 我都大概怎么搞

通流

你想拟合整个区域的物性？有这个必要吗？

cfkristen

原帖由 通流 于 2013-5-26 19:15 发表

                               
登录/注册后可看大图

你想拟合整个区域的物性？有这个必要吗？

嗯 有这个必要的 温度在540K到1020K左右的四个物性参数 输入进去 要不版主你说怎么搞啊 物性参数在整个区域是随温度变化的啊

浪子轻狂

从读书时起，一直对科研没什么兴趣，小木虫一次也没上过，所以就不过去看了，你这个直接在物性参数库里把物性参数设置成多项式不就完事了吗，写了这么长一个udf岂不麻烦？回流可能由如下几种情况造成：
1、初始化不合理，从一开始就回流，魔鬼缠身一样一直存在；
2、离散格式、边界上的湍流参数设置不合理，一直回流，计算甚至会发散；
3、计算区域不够长；
4、边界层网格不够细。

通流

在临界点附近，物性是没有办法用多项式表达的（其实也有可能，只是表达式不能是温度的函数）。因为，临界点是个奇点。在临界点，比热是无穷，音速为零。

飞翔鸟

超临界没有什么难搞的，说白了就是在临界点前后的剧烈物性变化，自己写进去了就好。

浪子轻狂

原帖由 通流 于 2013-5-27 01:18 发表

                               
登录/注册后可看大图

在临界点附近，物性是没有办法用多项式表达的（其实也有可能，只是表达式不能是温度的函数）。因为，临界点是个奇点。在临界点，比热是无穷，音速为零。

不是表达式，那udf是没辙。可以写一个数据点组成的profile文本导进去用。

5超7

比热的宏为什么要用？我在论坛听到人反映过，比热这个宏不要轻易用，计算的时候可能有问题。

超超

楼主，我也遇到同样问题，求解答

没晓晓

楼主最后解决了吗？我现在也做这方面的

不挑食的小瘦子

楼主问题解决了吗？