仿真在线:有限元分析咨询与培训中心——Ansys|Fluent|Abaqus|HFSS|SolidWorks|课程培训

fluent中使用非牛頓流體

2017-04-24  by:CAE仿真在線  來源:互聯網



1、非牛頓流體:剪應力與剪切應變率之間滿足線性關系的流體稱為牛頓流體,而把不滿足線性關系的流體稱為非牛頓流體。

2、fluent中使用非牛頓流體

a、層流狀態:直接在材料物性下設置材料的粘度,設置其為非牛頓流體。

b、湍流狀態

fluent在設置湍流模型后,會自動將材料的非牛頓流體性質直接改成了牛頓流體,因此需要做一些修改。最基本的方式有兩種:1、打開隱藏的湍流模型下非牛頓流體功能;2,直接利用UDF宏DEFINE_PROPERTY定義

3、打開隱藏的湍流模型下非牛頓流體功能

方法為:

(1) 在湍流模型中選擇標準的k-e模型;(2) 在Fluent窗口輸入命令:define/models/viscous/turbulence-expert/turb-non-newtonian 然后回車。

(3) 輸入:y 然后回車。


4、利用DEFINE_PROPERTY宏

A:這是一個自定義材料的粘度程序如下,也許對你有幫助。
在記事本中編輯的,另存為“visosity1.c"
#include "udf.h"
DEFINE_PROPERTY(cell_viscosity, cell, thread)
{
real mu_lam;
real trial;
rate=CELL_STRAIN_RATE_MAG(cell, thread);
real temp=C_T(cell, thread);
mu_lam=1.e12;
{
if(rate>1.0e-4 && rate<1.e5)
trial=12830000./rate*log(pow((rate*exp(17440.46/temp)/1.535146e8),0.2817)+pow((1.+pow((rate*exp(17440.46/temp)/1.535146e8),0.5634)),0.5));
else if (rate>=1.e5)
trial=128.3*log(pow((exp(17440.46/temp)/1.535146e8),0.2817)+pow((1.+pow((exp(17440.46/temp)/1.535146e8),0.5634)),0.5));
else
trial=1.283e11*log(pow((exp(17440.46/temp)/1.535146e12),0.2817)+pow((1.+pow((exp(17440.46/temp)/1.535146e12),0.5634)),0.5));
}
else if(temp>=855.&&temp<905.)
{
if(rate>1.0e-4 && rate<1.e5)
trial=12830000./rate*log(pow((rate*4.7063),0.2817)+pow((1.+pow((rate*4.7063),0.5634)),0.5))*pow(10.,-0.06*(temp-855.));
else if (rate>=1.e5)
trial=243.654*pow(10.,-0.06*(temp-855.));
else
trial=1.47897e10*pow(10.,-0.06*(temp-855.));
}
else if(temp>=905.)
{
if(rate>1.0e-4 && rate<1.e5)
trial=12830./rate*log(pow((rate*4.7063),0.2817)+pow((1.+pow((rate*4.7063),0.5634)),0.5));
else if (rate>=1.e5)
trial=0.24365;
else
trial=1.47897e7;
}
if(trial<1.e12&&trial>100.)
mu_lam=trial;
else if(trial<=1.)
mu_lam=1.;
else
mu_lam=1.e12;
return mu_lam;
}

B:在Fluent中使用Herschel-Bulkley粘性模型:

/* UDF for Herschel-Bulkley viscosity */
#include "udf.h”
real T,vis, s_mag, s_mag_c, sigma_y,n,k;
real C_1 = 1.0;
real C_2 = 1.0;
real C_3 = 1.0;
real C_4 = 1.0;
int ia ;
DEFINE_PROPERTY(hb_viscosity,c,t)
{
T=C_T(c, t);
s_mag = CELL_STRAIN_RATE_MAG(c,t);
/* Input parameters for H-B Viscosity */
if (ia==0.0)
{ C_1 = RP_Get_Real("c_1");
C_2 = RP_Get_Real("c_2");
C_3 = RP_Get_Real("c_3");
C_4 = RP_Get_Real("c_4");
ia = 1;}
k= C_1 ;
n= C_2 ;
sigma_y = C_3 ;
s_mag_c = C_4 ;
if (s_mag < s_mag_c)
{vis = sigma_y*(2-s_mag/s_mag_c)/s_mag_c+k*((2-n)+(n-1)*s_mag/s_mag_c)*pow(s_mag_c,(n-1));}
else
{ vis = sigma_y / s_mag + k*pow(s_mag, (n-1));}
return vis;
}

這是一個自定義材料的粘度程序如下,也許對你有幫助。
在記事本中編輯的,另存為“visosity1.c"
#include "udf.h"
DEFINE_PROPERTY(cell_viscosity, cell, thread)
{
real mu_lam;
real trial;
rate=CELL_STRAIN_RATE_MAG(cell, thread);
real temp=C_T(cell, thread);
mu_lam=1.e12;
{
if(rate>1.0e-4 && rate<1.e5)
trial=12830000./rate*log(pow((rate*exp(17440.46/temp)
/1.535146e8),0.2817)+pow((1.+pow((rate*exp(17440.46/
temp)/1.535146e8),0.5634)),0.5));
else if (rate>=1.e5)
trial=128.3*log(pow((exp(17440.46/temp)/
1.535146e8),0.2817)+pow((1.+pow((exp(17440.46/temp)/
1.535146e8),0.5634)),0.5));
else
trial=1.283e11*log(pow((exp(17440.46/temp)/
1.535146e12),0.2817)+pow((1.+pow((exp(17440.46/temp)/
1.535146e12),0.5634)),0.5));
}
else if(temp>=855.&&temp<905.)
{
if(rate>1.0e-4 && rate<1.e5)
trial=12830000./rate*log(pow((rate*4.7063),0.2817)+
pow((1.+pow((rate*4.7063),0.5634)),0.5))*pow(10.,-0.06
*(temp-855.));
else if (rate>=1.e5)
trial=243.654*pow(10.,-0.06*(temp-855.));
else
trial=1.47897e10*pow(10.,-0.06*(temp-855.));
}
else if(temp>=905.)
{
if(rate>1.0e-4 && rate<1.e5)
trial=12830./rate*log(pow((rate*4.7063),0.2817)+
pow((1.+pow((rate*4.7063),0.5634)),0.5));
else if (rate>=1.e5)
trial=0.24365;
else
trial=1.47897e7;
}
if(trial<1.e12&&trial>100.)
mu_lam=trial;
else if(trial<=1.)
mu_lam=1.;
else
mu_lam=1.e12;
return mu_lam;
}

.Fluent中使用Herschel-Bulkley粘性模型:

/* UDF for Herschel-Bulkley viscosity */
#include "udf.h”
real T,vis, s_mag, s_mag_c, sigma_y,n,k;
real C_1 = 1.0;
real C_2 = 1.0;
real C_3 = 1.0;
real C_4 = 1.0;
int ia ;
DEFINE_PROPERTY(hb_viscosity,c,t)
{
T=C_T(c, t);
s_mag = CELL_STRAIN_RATE_MAG(c,t);
/* Input parameters for H-B Viscosity */
if (ia==0.0)
{ C_1 = RP_Get_Real("c_1");
C_2 = RP_Get_Real("c_2");
C_3 = RP_Get_Real("c_3");
C_4 = RP_Get_Real("c_4");
ia = 1;}
k= C_1 ;
n= C_2 ;
sigma_y = C_3 ;
s_mag_c = C_4 ;
if (s_mag < s_mag_c)
{vis = sigma_y*(2-s_mag/s_mag_c)/s_mag_c+k*((2-n)+(n-1)*s_mag/s_mag_c)*pow(s_mag_c,(n-1));}
else
{ vis = sigma_y / s_mag + k*pow(s_mag, (n-1));}
return vis;
}

另外,如果在Fluent中使用了湍流模型,那么Fluent將默認地不允許使用非牛頓流體模型,你的UDF出現CELL_STRAIN_RATE_MAG: undeclared variable這樣的錯誤也可能與此有關,你需要激活在湍流模型條件下使用非牛頓流體模型的選項:

方法為:

(1) 在湍流模型中選擇標準的k-e模型;

(2) 在Fluent窗口輸入命令:define/models/viscous/turbulence-expert/turb-non-newtonian 然后回車。

(3) 輸入:y 然后回車。



開放分享:優質有限元技術文章,助你自學成才

相關標簽搜索:fluent中使用非牛頓流體 Fluent培訓 Fluent流體培訓 Fluent軟件培訓 fluent技術教程 fluent在線視頻教程 fluent資料下載 fluent分析理論 fluent化學反應 fluent軟件下載 UDF編程代做 Fluent、CFX流體分析 HFSS電磁分析 

編輯
在線報名:
  • 客服在線請直接聯系我們的客服,您也可以通過下面的方式進行在線報名,我們會及時給您回復電話,謝謝!
驗證碼

全國服務熱線

1358-032-9919

廣州公司:
廣州市環市中路306號金鷹大廈3800
電話:13580329919
          135-8032-9919
培訓QQ咨詢:點擊咨詢 點擊咨詢
項目QQ咨詢:點擊咨詢
email:kf@1cae.com




主站蜘蛛池模板: 碎料机_拌料机_吸料机_温控箱_烤料斗_模温机_东莞索亿塑料破碎机厂家 | 轻型防化服|重型防化服|全封闭防化服|济南三安安全防护设备有限公司 | 紫外线光疗仪|白癜风光疗仪|牛皮癣治疗仪|308纳米led|SIGMA|上海希格玛高技术有限公司 | 石材雕刻机_墓碑雕刻机_木工雕刻机_雕刻机厂家-合肥沃力数控设备有限责任公司 | 商标注册_北京商标注册费用_申请商标注册代理机构_北京商标注册公司- | 深圳心理咨询-专业心理咨询服务平台「专家免费在线」-深圳从心开始心理O2O | 实木运动地板价格_运动木地板厂家-欧氏实木运动地板网 | 在线BMI计算器 - 免费计算身体质量指数(BMI)| 金刚石线切割-电火花穿孔机-中走丝线切割-泰州锐特 | 气泡膜机_气垫膜机_缠绕膜机_瑞安市优特机械厂 | 在线计算网 - 超多超好用的在线小工具 | 潍坊卓瑞机械有限公司,输送设备,石灰消化设备,餐厨垃圾设备,化机浆设备,污泥脱水 | 质保_天津光栅尺_机床光栅尺数显表-天津方韦仪器有限公司 | 树脂井盖,复合井盖,井盖厂家-山东宝盖新材料 | 无缝钢管生产厂家-20G-L360N无缝钢管价格-天津海诚无缝钢管集团 无尘投料站-真空上料机-旋振筛|超声波振动筛|摇摆筛|筛分机-新乡市欧霖佳机械有限公司 | 艺术网 - 大型艺术类权威门户站 艺考培训-中影人教育 【官网】-中国艺考教育的引航者 | 山东汇河环保科技集团有限公司,水囊水袋,水罐,油囊,预压水袋,吊重水袋_山东汇河环保科技集团有限公司,水囊水袋,水罐,油囊,预压水袋,吊重水袋 | 四方光电(武汉)仪器有限公司_四方仪器首页-烟气分析仪|尾气分析仪|煤气分析仪|沼气分析仪|天然气分析仪|超声波流量计|在线气体分析系统|红外气体传感器 | 上海海外置业展_2024海外置业移民留学展_4月5-7日_企业参展处/免费领门票 | 江西铭鑫冶金设备有限公司-破碎机,铜米机,选矿摇床,电池回收设备 | 日本国际高中学校招生【立学仕日本留学服务】日语中学学校|美术生|艺术生|高中生|本科|研究生|语言学校留学择校申请 | 湖南实验台-防静电工作台-实验设备厂家-长沙实验室设备有限公司-湖南贝塔实验室设备有限公司 | 内蒙古碧云食品有限公司 | 宜昌江峡船用机械有限责任公司欢迎您! | 唐山搬家公司-搬家公司电话价格-四通搬家 | 伸缩机厂家-现货快递物流伸缩皮带机|装车输送机-乐清市华森自动化设备有限公司 | 停车场设计|地下停车场规划设计|智能停车位地下车库设计施工装修|深圳创安顺停车库规划设计有限公司 | 企多网-B2B网站,B2B商务平台免费b2b,推荐注册送积分,免费注册企业商铺, | 久久91精品久久91综合_国产亚洲自拍一区_国产精品第1页_亚洲高清视频一区_91成人午夜在线精品_亚洲国产精品网站在线播放_亚洲国产成人久久综合区_国产精品亚洲专区在线观看_免费视频精品一区二区三区 | 金相切割机-金相磨抛机-显微/维氏/布氏/洛氏硬度计-自准直仪-金相显微镜-万能材料试验机-清洁度检测仪-淋雨试验机-上海中研精密仪器制造有限公司 | 内衬不锈钢复合管,大口径内衬不锈钢复合钢管,双金属复合管,内衬不锈钢复合管厂家-江苏新澎 | 加药装置-排油烟装置-工业滤水器-连云港灵动 | 山西洗煤设备_山西洗煤机_山西离心机-山西东鑫昶机械科技股份有限公司 | 宁波雷豹机电科技有限公司|雷豹冷风机|雷豹工业大风扇|MFC18000|MFC16000|MFC6000|EF3622|EF4222|EF4822|移动工业蒸发式冷风机空气冷却器|大型工业空调扇|雷豹移动式工业大风扇|雷豹大风扇|生产厂家|公司官网 | 长沙变频器维修,变频器维修,ABB变频器维修,西门子变频器维修,施耐德变频器维修,伺服驱动器维修,工业机器人维修,20年专业工控电气维修,长沙文铖电气设备有限公司_长沙文铖电气设备有限公司 | 漫画大全免费看_漫画终身版_羞羞漫画网 | 破碎机设备-锤式颚式反击式圆锥移动冲击式破碎机厂家-成都大宏立机器公司 | 山东洗地机_工业洗地机_驾驶式扫地机_扫地车厂家_鼎洁盛世官网 | 上海律师_上海法律咨询_律师在线咨询网站_上海律师事务所-沪律网 | 无人叉车|智能无人叉车|智能AGV叉车|激光叉车AGV厂家_江西丹巴赫机器人股份有限公司 | 清尼龙滤膜-清洁度检测设备-清洁度分析仪-清洁度萃取机-优昂(百科) |