当前位置:首页 >> 能源/化工 >>

FLUENT水沸腾相变算例C语言编辑代码


FLUENT 水沸腾相变算例 C 语言编辑代码,想要就过来

#include "udf.h" /*包括常规宏*/ #include "sg_mphase.h" /*包括体积分数宏*/ #define T_SAT 373 /*定义饱和温度*/ #define LAT_HT 1.e3 /*定义水蒸汽潜热*/ DEFINE_SOURCE(liq_src, cell, pri_th, dS, eqn) /*定义液相源项*/ { Thread *mix_th, *sec_th; /*混合相、第二相 定义计算区域指针*/ real m_dot_l; /*蒸发冷凝速率 定义液相质量转移*/ mix_th = THREAD_SUPER_THREAD(pri_th); /*指向混合区的主相即液相的指针*/ sec_th = THREAD_SUB_THREAD(mix_th, 1); /* 指向单相控制区的气相的指针*/ if(C_T(cell, mix_th)>=T_SAT){ m_dot_l = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)* fabs(C_T(cell, pri_th) - T_SAT)/T_SAT; /*如果液相单元的温度高于蒸发温度, 液 相向气相的质量转移*/ dS[eqn] = -0.1*C_R(cell, pri_th)* fabs(C_T(cell, pri_th) - T_SAT)/T_SAT; /*定义源项对质量转移偏导*/ } else { m_dot_l = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)* fabs(T_SAT-C_T(cell,mix_th))/T_SAT; dS[eqn] = 0.; /*于是气相向液相转移,所以液相的质量源项对质量转移的偏导为零*/ } return m_dot_l; } DEFINE_SOURCE(vap_src, cell, sec_th, dS, eqn) /*定义气相源项*/ { Thread * mix_th, *pri_th; /*混合相、第一相*/ real m_dot_v; mix_th = THREAD_SUPER_THREAD(sec_th); /*指向混合区的第二相即气相的指针*/ pri_th = THREAD_SUB_THREAD(mix_th, 0); /*指向单相控制区的液相的指针,液相为主 相*/ if(C_T(cell, mix_th)>=T_SAT){ m_dot_v = 0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)* fabs(C_T(cell, mix_th) - T_SAT)/T_SAT; /*如果混合区单元的温度高于蒸发温度, 液相向气相的质量质量转移*/ dS[eqn] = 0.; /*由于是液相向气相转移,所以气相的质量源项对来自液相的质量转移的 偏导为零*/ } else { m_dot_v = -0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*

fabs(T_SAT-C_T(cell,mix_th))/T_SAT; dS[eqn] = -0.1*C_R(cell, sec_th)* fabs(C_T(cell, sec_th) - T_SAT)/T_SAT; /*由于是气相向液相转移,所以气相的质量源 项对自身的质量转移的偏导不为零*/ } return m_dot_v; } DEFINE_SOURCE(enrg_src, cell, mix_th, dS, eqn) /*混合模型能量源项*/ { Thread *pri_th, *sec_th; real m_dot; pri_th = THREAD_SUB_THREAD(mix_th, 0); /*指向混合区的液相的指针*/ sec_th = THREAD_SUB_THREAD(mix_th, 1); /*指向混合区的气相的指针*/ if(C_T(cell, mix_th)>=T_SAT){ m_dot = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)* fabs(C_T(cell, pri_th) - T_SAT)/T_SAT; /*如果混合区的单元温度高于蒸发温度。 质量转移由液相向气相转移,吸热,质量转移量前有负号*/ dS[eqn] = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)/T_SAT; } else { m_dot = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)* fabs(T_SAT-C_T(cell,mix_th))/T_SAT; dS[eqn] = -0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)/T_SAT;} /*相反, 气相向液相转移 则放热*/ return LAT_HT*m_dot; /*能量大小气化潜热与质量转移率相乘得能量源项*/



相关文章:
fluent 简单算例
传热与相变、化学反应与燃烧、多相流、旋转机械、动/变 形网格、噪声、材料...(UDF) ; FLUENT 软件采用 C/C++语言编写,从而大大提高了对计算机内存的利用率...
FLUENT算例 (5)搅拌桨底部十字挡板的流场分析
将材料设成水, motion type 设成 moving reference...相流模型,混合多相流模型,颗粒相模型,空穴两 相流...(UDF); FLUENT 软件采用 C/C++语言编写,从而大大...
FLUENT算例
一:前处理——利用 GAMBIT 建立计算模型,其网格画法 如下图 换热器的网格图 设置边界类型,如下图所示 定义边界类型设置对话框 二:利用 FLUENT 进行混合器内流动...
ansys fluent中的气穴算例之中文版
c、关闭 FLUENT Database Materials 窗口. B、Density 中输入 中输入0.02558 ...2、Specify water vapor as the secondary phase.指定水蒸气为第二相 Phases》...
FLUENT算例 (13)
FLUENT算例 (3) 14页 1下载券 fluent实例 22页 ...(c) (d) (a)间隙内充入水 (b) 间隙内充入...差, 这与表 1 中计算得到 的导热热阻的大小相...
FLUENT算例 (4)
FLUENT算例 (4)_工程科技_专业资料。Fluent算例课程:计算流体力学 题目:FLUENT...容器内部 4 3 的介质密度为 1000 kg / m 的流体介质,其定压比热 C P ?...
fluent多项流算例
fluent多项流算例_IT/计算机_专业资料。初学者骤...(a)指定水为基本相 选择 Phase-1,并点击 Set ...在 Interpreted UDFs 面板中,编辑 UDF (*.c)" ...
Fluent_操作手册
网格代码FLUENT 4 生成的网格,那么在 FLUENT 中...FLUENT 处理这类网格的技巧和滑移网格模型的技巧相 ...在 Figure 3 的例子中,界面区域由面 A-B、B-C...
fluent算例-船舶航行
fluent算例-船舶航行_电力/水利_工程科技_专业资料...项和第一相 基本相是水(water) ,第一相是空气(...压力分布云图 3.16 图 3.16 c) 速度分布矢量图...
fluent中多孔介质设置问题和算例
fluent中多孔介质设置问题和算例_机械/仪表_工程科技_专业资料。多孔介质流体仿真...(c) Close the Iterate panel. 7. Save the case and data files (...
更多相关标签: