联系我们 Fluent案例｜挥发冷凝UDF
联系我们
本案例演示在Fluent中哄骗UDF添加质料源项及能量源项的形势模拟水沸腾联系我们。
注：本算例来自Fluent官方老案例，仅为演示源项法模拟水沸腾。在早期版块的Fluent中莫得现成的挥发冷凝模子，需要借助UDF指定质料源与能量源来模拟沸腾气候。在较新的Fluent版块中依然内置了挥发冷凝模子，本算例中的UDF等同于Fluent中内置的Lee模子。本案例的UDF也不错历程修改用于更复杂的沸腾气候模拟。
”1 问题形色盘算模子如下图所示。
图片
盘算区域内介质为液态水，开动温度372 K，两侧面为绝热壁面。容器底部中心区域为高温T=573 K。
2 Fluent诞生以2D、Double Precision形势启动Fluent哄骗菜单File → Read → Mesh…读取网格文献boil.msh2.1 General诞生弃取Transient，指定重力加快度为Y场所-9.81 m/s2图片
2.2 Models诞生激活能量方程图片
选拔层流盘算图片
选拔Mixture多相流模子，激活选项Implicit Body Force图片
上升选手之中李淑瑛在女子韩巡乐天公开赛中并列位于39位，获得0.62分，世界排名从368位上升到354位。
亚军玛丽娅-埃尔南德斯（Maria Hernandez）来自西班牙，获得10.8分，世界排名从251位上升到188位。
2.3 Materials诞生添加水与水蒸气图片
按下图所示修改水蒸气材料参数图片
注：本算例哄骗UDF指定传质与传热，因此圭臬情景焓与参考温度不错率性诞生，不参与盘算。
”添加液态水，按下图所示指定参数图片
2.4 诞生多相流模子指定多相流模子主相为water-liquid图片
指定次相为vapor，诞生Diameter为0.0002 m图片
相间作用如下图所示，不需要诞生传质与传热图片
2.5 诠释注解UDF右键弃取模子树节点User Defined Functions，点击弹出菜单项Interpreted...掀开UDF诠释注解对话框图片
在掀开的对话框中添加源文献source.c，点击按钮Interpret诠释注解源代码图片
2.6 诞生盘算区域诞生流体盘算区域fluid，弃取Phase为mixture，激活选项Source Terms，软件定制开发点击Energy后的按钮Edit...裁剪能量源图片
指定能量源为udf enrg_src，点击OK按钮关闭对话框图片
小程序开发诞生Phase为liquid，如下图所示诞生Mass后方的Edit... 按钮添加质料源图片
诞生液相质料源为udf liq_src图片
增多气相质料源图片
诞不满相质料源为udf vap_src图片
2.7 边界条款诞生诞生出口边界poutlet的出口静压为0 Pa图片
指定出口回流温度为372 K图片
指定出语气相回流体积分数为0图片
指定边界wall-hot的温度为573 K图片
2.8 Methods诞生盘算秩序诞生图片
2.9 秀丽区域右键弃取模子树节点Cell Registers，点击弹出菜单项New → Boundary…掀开诞生对话框图片
如下图所示秀丽边界区域wall-hot相邻的第一层网格图片
2.10 开动化全局开动化图片
将秀丽的网格区域温度开动化为373.15 K图片
2.11 最先盘算指定手艺步长0.005 s，手艺步数为1000，进行盘算图片
2.12 盘算恶果盘算恶果如下图所示图片
3 UDF要领简解这个UDF照旧相比浅薄的，主要包含两个质料源与一个能量源。盘算公式选拔的是Lee模子。
UDF要领如下所示。
#include "udf.h"#include "sg_mphase.h"#define T_SAT 373.15 //富余温度373.15 K#define LAT_HT 1.e3 //汽化潜热1000 J/kg // 液相质料源项DEFINE_SOURCE(liq_src, cell, pri_th, dS, eqn){ Thread *mix_th, *sec_th; real m_dot_l; // mix_th存储主相的Thread指针；sec_th存储次相指针，注视ID mix_th = THREAD_SUPER_THREAD(pri_th); sec_th = THREAD_SUB_THREAD(mix_th, 1); // 温度大于富余温度，示意为挥发 if (C_T(cell, mix_th) >= T_SAT) { //盘算挥发量，注视液相质料流量为负值，0.1为coeff值 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) { // 获得挥发量，Lee模子 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;}文献下载贯串：https://pan.baidu.com/s/189oVaShGNLP2tmYAFu83uw 索取码：n77p
”本站仅提供存储管事，系数实质均由用户发布，如发现存害或侵权实质，请点击举报。