@李东岳 东岳老师您好,我知道Luo破碎模型的子代分布函数是根据破碎频率计算得出的。但是Liao模型只给出了总体破碎频率无法根据单一气泡破碎频率进行计算。这点使我无法理解。
CFD_xiaobai
帖子
-
-
请教下各位大佬,Baseline closure model for dispersed bubbly flow_ Bubble coalescence__and breakup. 【Yixiang Liao Roland Rzehak, Dirk Lucas, Eckhard Krepper-CES-2015-122-336-349】这篇文章中,作者仅给出破碎速率表达式,请问气泡子尺寸分布函数是如何求出的?
-
@李东岳 原文作者是用CFX做的
-
。请教一下,这篇论文说用LES作为湍流模型分别计算气相和液相的湍流粘度。但是Fluent中LES只能用于单相流,不知道这里是如何实现的。求指点
-
@李东岳 好的,谢谢东岳老师
-
请教下各位,本人的圆管内带有螺旋插件。一开壁面润滑力模型计算几步就发散。不知道是什么原因,恳请赐教。
-
@李东岳 东岳老师,目前本人遇到的问题是湍流模型选取SST-k-omega模型。但是需要调用epsi这个数据,所以只使用k和omega进行计算。现在是认为omega=0.09* k*omega不知道可不可以。希望东岳老师指点下。
-
@李东岳 感谢李老师
-
请教下湍流耗散率(epsi),湍流能(k)和湍流比耗散率(omega)三者之间的关系
-
请问各位有遇到过这个问题吗?谢谢各位答疑
-
这个是气泡破碎频率@cfd_xiaobai
-
请教各位,本人用udf导入的破碎模型,没有考虑气泡的聚并,但是bin方程残差明显有问题。但不知道原因是什么。除了bin方程外,其余方程没有问题。拜托各位大佬解答下。
参考文献 . Theoretical model for drop and bubble breakup in turbulent dispersions. AIChE Journal 42,1996 1225–1233. -
@李东岳 请问东岳老师遇到过,class method中除bin0方程外,其余bin方程残差都一致的情况吗?
-
#include "udf.h" #include "sg_pb.h" #include "sg_mphase.h" #define rhoC 998.2 #define sigma 0.072 #define rho_d 1.2 DEFINE_PB_BREAK_UP_RATE_FREQ(break_up_freq_tav, cell, thread, d_1) { Thread *mixture=THREAD_SUPER_THREAD(thread); real k,epsi,f1,f2,f3,f4,f5,f6,f7,f8,f9,f10,f11,omega,alpha,psi,itea; real pi=3.14159,hf=1.0e-8,h0=1.0e-4; real C1 = 0.675, C2 = 0.1328; int i,s=0; Thread *tm; omega=C_O(cell,mixture); k=C_K(cell, mixture); epsi=0.09*k*omega; alpha=C_VOF(cell,thread); itea=11.4*pow(10,-4.5)/pow(epsi,1/4); psi=pow(epsi,2/3)*pow(d_1,5/3); f1=0.1*pow((1+0.2),2)/pow(0.2,11/3)*(1+exp(-0.26*0.00042/psi/pow(0.2,11/3))+2*exp(-0.1476*0.00042/psi/pow(0.2,11/3))+2*exp(-0.2038*0.00042/psi/pow(0.2,11/3))+2*exp(-0.2365*0.00042/psi/pow(0.2,11/3))+2*exp(-0.2543*0.00042/psi/pow(0.2,11/3)) ); f2=0.1*pow((1+0.3),2)/pow(0.3,11/3)*(1+exp(-0.26*0.00042/psi/pow(0.3,11/3))+2*exp(-0.1476*0.00042/psi/pow(0.3,11/3))+2*exp(-0.2038*0.00042/psi/pow(0.3,11/3))+2*exp(-0.2365*0.00042/psi/pow(0.3,11/3))+2*exp(-0.2543*0.00042/psi/pow(0.3,11/3)) ); f3=0.1*pow((1+0.4),2)/pow(0.4,11/3)*(1+exp(-0.26*0.00042/psi/pow(0.4,11/3))+2*exp(-0.1476*0.00042/psi/pow(0.4,11/3))+2*exp(-0.2038*0.00042/psi/pow(0.4,11/3))+2*exp(-0.2365*0.00042/psi/pow(0.4,11/3))+2*exp(-0.2543*0.00042/psi/pow(0.4,11/3)) ); f4=0.1*pow((1+0.5),2)/pow(0.5,11/3)*(1+exp(-0.26*0.00042/psi/pow(0.5,11/3))+2*exp(-0.1476*0.00042/psi/pow(0.5,11/3))+2*exp(-0.2038*0.00042/psi/pow(0.5,11/3))+2*exp(-0.2365*0.00042/psi/pow(0.5,11/3))+2*exp(-0.2543*0.00042/psi/pow(0.5,11/3)) ); f5=0.1*pow((1+0.6),2)/pow(0.6,11/3)*(1+exp(-0.26*0.00042/psi/pow(0.6,11/3))+2*exp(-0.1476*0.00042/psi/pow(0.6,11/3))+2*exp(-0.2038*0.00042/psi/pow(0.6,11/3))+2*exp(-0.2365*0.00042/psi/pow(0.6,11/3))+2*exp(-0.2543*0.00042/psi/pow(0.6,11/3)) ); f6=0.1*pow((1+0.7),2)/pow(0.7,11/3)*(1+exp(-0.26*0.00042/psi/pow(0.7,11/3))+2*exp(-0.1476*0.00042/psi/pow(0.7,11/3))+2*exp(-0.2038*0.00042/psi/pow(0.7,11/3))+2*exp(-0.2365*0.00042/psi/pow(0.7,11/3))+2*exp(-0.2543*0.00042/psi/pow(0.7,11/3)) ); f7=0.1*pow((1+0.8),2)/pow(0.8,11/3)*(1+exp(-0.26*0.00042/psi/pow(0.8,11/3))+2*exp(-0.1476*0.00042/psi/pow(0.8,11/3))+2*exp(-0.2038*0.00042/psi/pow(0.8,11/3))+2*exp(-0.2365*0.00042/psi/pow(0.8,11/3))+2*exp(-0.2543*0.00042/psi/pow(0.8,11/3)) ); f8=0.1*pow((1+0.9),2)/pow(0.9,11/3)*(1+exp(-0.26*0.00042/psi/pow(0.9,11/3))+2*exp(-0.1476*0.00042/psi/pow(0.9,11/3))+2*exp(-0.2038*0.00042/psi/pow(0.9,11/3))+2*exp(-0.2365*0.00042/psi/pow(0.9,11/3))+2*exp(-0.2543*0.00042/psi/pow(0.9,11/3)) ); f9=0.2*(1+exp(-0.26*0.00042/psi)+2*exp(-0.1476*0.00042/psi)+2*exp(-0.2038*0.00042/psi)+2*exp(-0.2365*0.00042/psi)+2*exp(-0.2543*0.00042/psi) ); f10=(0.05-itea/2)*pow((1+itea),2)/pow(itea,11/3)*(1+exp(-0.26*0.00042/psi/pow(itea,11/3))+2*exp(-0.1476*0.00042/psi/pow(itea,11/3))+2*exp(-0.2038*0.00042/psi/pow(itea,11/3))+2*exp(-0.2365*0.00042/psi/pow(itea,11/3))+2*exp(-0.2543*0.00042/psi/pow(itea,11/3)) ); f11=(0.1-itea/2)*pow((1+0.1),2)/pow(0.1,11/3)*(1+exp(-0.26*0.00042/psi/pow(0.1,11/3))+2*exp(-0.1476*0.00042/psi/pow(0.1,11/3))+2*exp(-0.2038*0.00042/psi/pow(0.1,11/3))+2*exp(-0.2365*0.00042/psi/pow(0.1,11/3))+2*exp(-0.2543*0.00042/psi/pow(0.1,11/3)) ); return 0.923*pow(epsi,1/3)/pow(d_1,2/3)*(1-alpha)*(f1+f2+f3+f4+f5+f6+f7+f8+f9+f10+f11); } 各位大佬,请问我这个根据罗和安教授的破碎核函数写的破碎频率udf,可以正常导入fluent,计算方法使用fluent自带的CM方法,但是除了bin0,剩下的bin方程残差都一样,求解答。下面是子尺寸分布函数 #include "udf.h" #include "sg_pb.h" #include "sg_mphase.h" #define max(a, b) (((a) > (b)) ? (a) : (b)) #define min(a,b) ( ((a)>(b)) ? (b):(a) ) DEFINE_PB_BREAK_UP_RATE_PDF(break_up_pdf_par, cell, thread, d_1, thread_2, d_2) { Thread *mixture=THREAD_SUPER_THREAD(thread); real pdf; real f_v = min(d_2,pow(pow(d_1,3)-pow(d_2,3),1/3))/d_1; real f_bv =pow(f_v,3.) ; real c_f =pow(f_bv,2/3)+pow((1-f_bv),2/3)-1 ; real k,epsi,f1,f2,f3,f4,f5,f6,f7,f8,f9,f10,f11,f12,f13,f14,f15,f16,f17,f18,f19,f20,f21,f22,f23,f24,omega,alpha,psi,itea; Thread *tm; omega=C_O(cell,mixture); k=C_K(cell, mixture); alpha=C_VOF(cell,thread); epsi=0.09*k*omega; itea=11.4*pow(10,-4.5)/pow(epsi,1/4); psi=pow(epsi,2/3)*pow(d_1,5/3); f1=0.1*pow((1+0.2),2)/pow(0.2,11/3)*(1+exp(-0.26*0.00042/psi/pow(0.2,11/3))+2*exp(-0.1476*0.00042/psi/pow(0.2,11/3))+2*exp(-0.2038*0.00042/psi/pow(0.2,11/3))+2*exp(-0.2365*0.00042/psi/pow(0.2,11/3))+2*exp(-0.2543*0.00042/psi/pow(0.2,11/3)) ); f2=0.1*pow((1+0.3),2)/pow(0.3,11/3)*(1+exp(-0.26*0.00042/psi/pow(0.3,11/3))+2*exp(-0.1476*0.00042/psi/pow(0.3,11/3))+2*exp(-0.2038*0.00042/psi/pow(0.3,11/3))+2*exp(-0.2365*0.00042/psi/pow(0.3,11/3))+2*exp(-0.2543*0.00042/psi/pow(0.3,11/3)) ); f3=0.1*pow((1+0.4),2)/pow(0.4,11/3)*(1+exp(-0.26*0.00042/psi/pow(0.4,11/3))+2*exp(-0.1476*0.00042/psi/pow(0.4,11/3))+2*exp(-0.2038*0.00042/psi/pow(0.4,11/3))+2*exp(-0.2365*0.00042/psi/pow(0.4,11/3))+2*exp(-0.2543*0.00042/psi/pow(0.4,11/3)) ); f4=0.1*pow((1+0.5),2)/pow(0.5,11/3)*(1+exp(-0.26*0.00042/psi/pow(0.5,11/3))+2*exp(-0.1476*0.00042/psi/pow(0.5,11/3))+2*exp(-0.2038*0.00042/psi/pow(0.5,11/3))+2*exp(-0.2365*0.00042/psi/pow(0.5,11/3))+2*exp(-0.2543*0.00042/psi/pow(0.5,11/3)) ); f5=0.1*pow((1+0.6),2)/pow(0.6,11/3)*(1+exp(-0.26*0.00042/psi/pow(0.6,11/3))+2*exp(-0.1476*0.00042/psi/pow(0.6,11/3))+2*exp(-0.2038*0.00042/psi/pow(0.6,11/3))+2*exp(-0.2365*0.00042/psi/pow(0.6,11/3))+2*exp(-0.2543*0.00042/psi/pow(0.6,11/3)) ); f6=0.1*pow((1+0.7),2)/pow(0.7,11/3)*(1+exp(-0.26*0.00042/psi/pow(0.7,11/3))+2*exp(-0.1476*0.00042/psi/pow(0.7,11/3))+2*exp(-0.2038*0.00042/psi/pow(0.7,11/3))+2*exp(-0.2365*0.00042/psi/pow(0.7,11/3))+2*exp(-0.2543*0.00042/psi/pow(0.7,11/3)) ); f7=0.1*pow((1+0.8),2)/pow(0.8,11/3)*(1+exp(-0.26*0.00042/psi/pow(0.8,11/3))+2*exp(-0.1476*0.00042/psi/pow(0.8,11/3))+2*exp(-0.2038*0.00042/psi/pow(0.8,11/3))+2*exp(-0.2365*0.00042/psi/pow(0.8,11/3))+2*exp(-0.2543*0.00042/psi/pow(0.8,11/3)) ); f8=0.1*pow((1+0.9),2)/pow(0.9,11/3)*(1+exp(-0.26*0.00042/psi/pow(0.9,11/3))+2*exp(-0.1476*0.00042/psi/pow(0.9,11/3))+2*exp(-0.2038*0.00042/psi/pow(0.9,11/3))+2*exp(-0.2365*0.00042/psi/pow(0.9,11/3))+2*exp(-0.2543*0.00042/psi/pow(0.9,11/3)) ); f9=0.2*(1+exp(-0.26*0.00042/psi)+2*exp(-0.1476*0.00042/psi)+2*exp(-0.2038*0.00042/psi)+2*exp(-0.2365*0.00042/psi)+2*exp(-0.2543*0.00042/psi) ); f10=(0.05-itea/2)*pow((1+itea),2)/pow(itea,11/3)*(1+exp(-0.26*0.00042/psi/pow(itea,11/3))+2*exp(-0.1476*0.00042/psi/pow(itea,11/3))+2*exp(-0.2038*0.00042/psi/pow(itea,11/3))+2*exp(-0.2365*0.00042/psi/pow(itea,11/3))+2*exp(-0.2543*0.00042/psi/pow(itea,11/3)) ); f11=(0.1-itea/2)*pow((1+0.1),2)/pow(0.1,11/3)*(1+exp(-0.26*0.00042/psi/pow(0.1,11/3))+2*exp(-0.1476*0.00042/psi/pow(0.1,11/3))+2*exp(-0.2038*0.00042/psi/pow(0.1,11/3))+2*exp(-0.2365*0.00042/psi/pow(0.1,11/3))+2*exp(-0.2543*0.00042/psi/pow(0.1,11/3)) ); f23=f1+f2+f3+f4+f5+f6+f7+f2+f8+f9+f10+f11; f12=0.1*pow((1+0.2),2)/pow(0.2,11/3)*exp(-0.00042*c_f/psi/pow(0.2,11/3)); f13=0.1*pow((1+0.3),2)/pow(0.3,11/3)*exp(-0.00042*c_f/psi/pow(0.3,11/3)); f14=0.1*pow((1+0.4),2)/pow(0.4,11/3)*exp(-0.00042*c_f/psi/pow(0.4,11/3)); f15=0.1*pow((1+0.5),2)/pow(0.5,11/3)*exp(-0.00042*c_f/psi/pow(0.5,11/3)); f16=0.1*pow((1+0.6),2)/pow(0.6,11/3)*exp(-0.00042*c_f/psi/pow(0.6,11/3)); f17=0.1*pow((1+0.7),2)/pow(0.7,11/3)*exp(-0.00042*c_f/psi/pow(0.7,11/3)); f18=0.1*pow((1+0.8),2)/pow(0.8,11/3)*exp(-0.00042*c_f/psi/pow(0.8,11/3)); f19=0.1*pow((1+0.9),2)/pow(0.9,11/3)*exp(-0.00042*c_f/psi/pow(0.9,11/3)); f20=0.2*exp(-0.00042*c_f/psi); f21=(0.05-itea/2)*pow((1+itea),2)/pow(itea,11/3)*exp(-0.00042*c_f/psi/pow(itea,11/3)); f22=(0.1-itea/2)*pow((1+0.1),2)/pow(0.1,11/3)*exp(-0.00042*c_f/psi/pow(0.1,11/3)); f24=f12+f13+f14+f15+f16+f17+f18+f19+f20+f21+f22; pdf =f24/f23; return pdf; }
PBE中Liao-breakage 模型求助
PBE中Liao-breakage 模型求助
LES如何与Euler-Euler模型耦合,请大佬们指点下
LES如何与Euler-Euler模型耦合,请大佬们指点下
fluent欧拉多相流壁面润滑力选取
fluent欧拉多相流壁面润滑力选取
请教下湍流耗散率(epsi),湍流能(k)和湍流比耗散率(omega)三者之间的关系
请教下湍流耗散率(epsi),湍流能(k)和湍流比耗散率(omega)三者之间的关系
请教下湍流耗散率(epsi),湍流能(k)和湍流比耗散率(omega)三者之间的关系
fluent 出现increasing relaxion sweeps这个报错(仅在迭代中出现一次)
群体平衡模型luo破碎模型求助
群体平衡模型luo破碎模型求助
fluent中pbm破碎核udf求解
fluent中pbm破碎核udf求解