大家好,问一个很弱智的问题。。。用了openfoam的sampling surface,输出格式为raw。请问用什么软件处理呢?
CYW
帖子
-
-
@Cav Charin, Amanda HLM, et al. "On the dynamic behavior of rising droplets." International Journal of Multiphase Flow 110 (2019): 165-178.
你好, 请问你知道怎么计算气泡的升阻力吗 -
用interFoam模拟气泡上升,如何计算气泡的升阻力?文献里是用施加在气泡上的pressure和viscous stress来计算的。请问大神有没有知道如何操作的?用paraview可以直接计算吗?我看到了pressure,但是viscous stress如何计算呢?
-
最近在模拟有传热情况下额气泡上升问题,求解器为compressibleInterFoam。当连续计算时,出现如下错误,有没有大神指导一下?不胜感激!!
Lookup solver for TFinal
smoothSolver: Solving for T, Initial residual = 0.0076371, Final residual = 8.40008e-09, No Iterations 5
[60] #0 Foam::error::printStack(Foam::Ostream&)[36] #0 Foam::error::printStack(Foam::Ostream&)[15] #0 Foam::error::printStack(Foam::Ostream&) at ??:?
at ??:?
[15] #1 Foam::sigFpe::sigHandler(int)[36] #1 Foam::sigFpe::sigHandler(int) at ??:?
[60] #1 Foam::sigFpe::sigHandler(int) at ??:?
[15] #2 ? at ??:?
[36] #2 ? at ??:?
[60] #2 ?Lookup gradScheme for nHat
Cache: Calculating nHat, originating from alpha.water event No. 624
in /lib64/libc.so.6
[15] #3 in /lib64/libc.so.6
[36] #3 in /lib64/libc.so.6
[60] #3 ???Lookup interpolationScheme for interpolate(nHat)
at ??:?
[15] #4 Foam::H2O::rho(double, double) const at ??:?
[36] #4 Foam::H2O::rho(double, double) const at ??:?
[60] #4 Foam::H2O::rho(double, double) const at ??:?
[15] #5 Foam::species::thermo<Foam::thermophysicalPropertiesSelectorFoam::liquidProperties, Foam::sensibleInternalEnergy>::Es(double, double) const at ??:?
[36] #5 Foam::heThermo<Foam::rhoThermo, Foam::pureMixture<Foam::species::thermo<Foam::thermophysicalPropertiesSelectorFoam::liquidProperties, Foam::sensibleInternalEnergy> > >::he(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const at ??:?
[60] #5 Foam::heThermo<Foam::rhoThermo, Foam::pureMixture<Foam::species::thermo<Foam::thermophysicalPropertiesSelectorFoam::liquidProperties, Foam::sensibleInternalEnergy> > >::he(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const at ??:?
[15] #6 Foam::species::thermo<Foam::thermophysicalPropertiesSelectorFoam::liquidProperties, Foam::sensibleInternalEnergy>::T(double, double, double, double (Foam::species::thermo<Foam::thermophysicalPropertiesSelectorFoam::liquidProperties, Foam::sensibleInternalEnergy>::)(double, double) const, double (Foam::species::thermo<Foam::thermophysicalPropertiesSelectorFoam::liquidProperties, Foam::sensibleInternalEnergy>::)(double, double) const, double (Foam::species::thermo<Foam::thermophysicalPropertiesSelectorFoam::liquidProperties, Foam::sensibleInternalEnergy>::*)(double) const) const at ??:?
[36] #6 Foam::twoPhaseMixtureThermo::correctThermo() at ??:?
[60] #6 Foam::twoPhaseMixtureThermo::correctThermo() at ??:?
[36] #7 at ??:?
[60] #7 at ??:?
[15] #7 Foam::heRhoThermo<Foam::rhoThermo, Foam::pureMixture<Foam::species::thermo<Foam::thermophysicalPropertiesSelectorFoam::liquidProperties, Foam::sensibleInternalEnergy> > >::calculate(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, bool)?? at ??:?
[15] #8 Foam::heRhoThermo<Foam::rhoThermo, Foam::pureMixture<Foam::species::thermo<Foam::thermophysicalPropertiesSelectorFoam::liquidProperties, Foam::sensibleInternalEnergy> > >::correct() at ??:?
[36] #8 __libc_start_main at ??:?
[60] #8 __libc_start_main at ??:?
[15] #9 Foam::twoPhaseMixtureThermo::correctThermo() in /lib64/libc.so.6
[36] #9 in /lib64/libc.so.6
[60] #9 at ??:?
[15] #10 ??? at /home/abuild/rpmbuild/BUILD/glibc-2.22/csu/../sysdeps/x86_64/start.S:121
at /home/abuild/rpmbuild/BUILD/glibc-2.22/csu/../sysdeps/x86_64/start.S:121
at ??:?
[15] #11 __libc_start_main in /lib64/libc.so.6
[15] #12 ? at /home/abuild/rpmbuild/BUILD/glibc-2.22/csu/../sysdeps/x86_64/start.S:121 -
@桎梏 也可以用不可压缩的算例,不过气泡体积就不会变化了。你可以在interFoam加入能量方程,就可以算啦。
-
@桎梏 你好,目前还没有相变。求解器的话可以用compressibleInterFoam
-
@mohui 谢谢!
-
大家好,单个气泡处于无重力静水中,气泡跟水之间有温差如下图,请问气泡内的温度分布有没有理论解?
谢谢大家! -
大家好,下图模型,用blockMeshDict可以建立模型,然后topoSet 和createPatch可以建立速度入口,喷嘴怎么搞呢?非常感谢!
-
最近在interFoam中添加了传热模型,但是没找到相关的验证算例,有没有推荐的,谢谢!
-
@东岳
谢谢!这下明白了! -
@东岳 谢谢!然后这样并不能切出一个球。。。我试过用contour可以搞出一个球。那第二个图,气泡下边的流线呢?也有sphere?
-
@buaa_hgj 谢谢!这个是二维截面的,但也很有帮助。
-
谢谢!如果设置sphere的话,球心代表什么呢?我试了一下,如果放在圆心的话,流线还是很乱,并不能保证在气泡内。第二幅图也是用sphere做出来的吗?
-
大家好,想请教一下各位大神,下图中的流线图是怎么用paraview画出来的?试了很多次,paraview里的stream trace一直掌握不好,或者有没有大神推荐个教程之类的。谢谢各位!!
-
-
@MrZZQi 你好,目前这个问题我还没有解决。在忙别的事情了。。。
-
@东岳
在interfaceProperties.C中编程如下
fvc::interpolate(fvc::grad(sigma) - fvc::grad(alpha)*(fvc::grad(alpha)&fvc::grad (sigma))*mag(fvc::grad (sigma))) & mesh.Sf()
将此项添加到速度和压力方程。其中当计算压力项时发散了。 -
@东岳 是的。在VOF动量方程添加。n 是 后两项,我的想法是fvc::snGrad(sigma) - fvc::snGrad(alpha)*(fvc::snGrad(alpha)&fvc::snGrad (sigma)),但是编译一直不成功。
-
想研究一下surface tension gradient
-
@东岳 后两项
-
有人用openfoam模拟thermocapillary motion of deformable drops? Marangoni stress 需要考虑吗?直接在表面张力项加入
-
interfoam模拟俩气泡,想输出一下俩气泡重心位置,请问大家怎么计算?在controlDict中写function
-
starFrom:latestTime
-
我在interFoam加入了传热模型,在模拟单个无重力气泡加热的时候,气泡体积不变,请问可能的原因是什么?谢谢
-
@东岳 谢谢!我试了一下,效果如下。
但是怎么样可以产生一个气泡一个气泡的效果呢?调节速度进口值吗?
如果没有速度,可以由于进口压力大于水压从而产生气泡吗? -
各位大神,我想模拟一下气泡喷射到水里然后上升的过程,但是喷嘴如何设置呢?用什么边界条件?
-
求解相方程时
MULES::explicitSolve ( geometricOneField(), alpha1, phi, talphaPhiCorr.ref(), vDotvmcAlphal, (divU*alpha1 + vDotcAlphal)(), 1, 0 );
第二项phi为什么不是AlphaPhi?
-
@飞火流星jyj 我的理解是
vDotvmcAlphal = vDotvAlpha - vDotcAlpha (凝结与汽化的差)
而相方程的mDot = vDotvAlpha,故又加上了vDotcAlpha。
不知道这样理解对不对 -
@飞火流星jyj 啊,不好意思,我说的是另外一个solver,interCondensatingEvaporatingFoam.
-
@飞火流星jyj 谢谢,是差了个负号。mDotAlpha*Alpha为mDot,这个相变模型是modified Lee model,方程如下.
mDot即为Sg。质量流与体积分数有关。
如果我换一个相变模型,比如Tanasawa简化的Schrage模型,方程如下
这里的mDotAlpha应该是什么呢?难道是mDot再除以体积分数吗? -
各位大神,最近在研究interCondensatingEvaporatingFoam,有点疑问。
mDot
和mDotAlpha
差一个Alpha,mDotAlpha
为return Pair<tmp<volScalarField>> ( coeffC_*mixture_.rho2()*max(TSat - T.oldTime(), T0), -coeffE_*mixture_.rho1()*max(T.oldTime() - TSat, T0) );
mDot为
return Pair<tmp<volScalarField>> ( coeffC_*mixture_.rho2()*limitedAlpha2*max(TSat - T.oldTime(), T0), coeffE_*mixture_.rho1()*limitedAlpha1*max(T.oldTime() - TSat, T0) ); }
相方程源项为
mDotAlpha
乘以
1.0/mixture_.rho1() - mixture_.alpha1()*(1.0/mixture_.rho1() - 1.0/mixture_.rho2())
。mDot
与mDotAlpha
量纲相同,为什么相方程中用mDotAlpha
? -
@飞火流星jyj 你知道mDot 和mDotAlpha 的区别吗
-
@东岳 Thanks!
-
编译以下语句时
volVectorField gradAlpha(fvc::grad(alpha1_)): volScalarField mDot(this->mDot()); volScalarField vmDot(mDot*gradAlpha);
错误提示如下
no matching function for call to Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> >::tmp(Foam::tmp<Foam::Field<Foam::Vector<double> > >)
各位大神有知道为什么的吗?
-
@东岳 我对这个期刊也不熟。只是在CFDOnline上看到了这个solver。拿来用一下,结果老是发散。
Time = 0.0001 MULES: Solving for alpha1 Liquid phase volume fraction = 0.999782 Min(alpha1) = 0 Max(alpha1) = 1 GAMG: Solving for p_rgh, Initial residual = 1, Final residual = 9.17326e-07, No Iterations 32 GAMG: Solving for p_rgh, Initial residual = 4.56128e-07, Final residual = 4.56128e-07, No Iterations 0 time step continuity errors : sum local = 1.15564e-10, global = 4.67937e-12, cumulative = 4.67937e-12 GAMG: Solving for p_rgh, Initial residual = 0.000258322, Final residual = 7.79975e-07, No Iterations 4 GAMG: Solving for p_rgh, Initial residual = 7.79882e-07, Final residual = 7.79882e-07, No Iterations 0 time step continuity errors : sum local = 1.96234e-10, global = 4.41825e-12, cumulative = 9.09762e-12 GAMG: Solving for p_rgh, Initial residual = 3.261e-06, Final residual = 6.33275e-07, No Iterations 1 GAMGPCG: Solving for p_rgh, Initial residual = 6.33275e-07, Final residual = 9.32048e-09, No Iterations 5 time step continuity errors : sum local = 3.247e-11, global = -1.80463e-12, cumulative = 7.29299e-12 **** ****Pressure range: 313.455 Pa ****Max velocity: 0.168515 m/s ****Phase change energy: 0 W ****Volume change: 0 m^3/s DILUPBiCG: Solving for H, Initial residual = 0.00100367, Final residual = 7.26908e-09, No Iterations 9 DILUPBiCG: Solving for H, Initial residual = 0.000932397, Final residual = 7.34408e-09, No Iterations 9 DILUPBiCG: Solving for H, Initial residual = 0.000911756, Final residual = 7.5305e-09, No Iterations 9 ExecutionTime = 15.77 s ClockTime = 17 s Courant Number mean: 2.65292e-05 max: 0.283422 Interface Courant Number mean: 0 max: 0 Fourier number mean: 0.00271598 max: 1.81326 Time = 0.0002 MULES: Solving for alpha1 Liquid phase volume fraction = 0.999782 Min(alpha1) = -1.24343e-09 Max(alpha1) = 1 GAMG: Solving for p_rgh, Initial residual = 0.158707, Final residual = 9.59991e-07, No Iterations 33 GAMG: Solving for p_rgh, Initial residual = 9.90688e-07, Final residual = 9.90688e-07, No Iterations 0 time step continuity errors : sum local = 2.63908e-08, global = -1.91596e-08, cumulative = -1.91524e-08 GAMG: Solving for p_rgh, Initial residual = 0.00214601, Final residual = 9.69321e-07, No Iterations 12 GAMG: Solving for p_rgh, Initial residual = 9.70692e-07, Final residual = 9.70692e-07, No Iterations 0 time step continuity errors : sum local = 2.63875e-08, global = -1.9196e-08, cumulative = -3.83484e-08 GAMG: Solving for p_rgh, Initial residual = 0.000147321, Final residual = 6.99269e-07, No Iterations 4 GAMGPCG: Solving for p_rgh, Initial residual = 6.99262e-07, Final residual = 4.28891e-09, No Iterations 6 time step continuity errors : sum local = 2.62568e-08, global = -1.91768e-08, cumulative = -5.75252e-08 **** ****Pressure range: 314.046 Pa ****Max velocity: 0.158209 m/s ****Phase change energy: -0.0533324 W ****Volume change: -4.60434e-10 m^3/s DILUPBiCG: Solving for H, Initial residual = 0.0112104, Final residual = 7.17193e-09, No Iterations 11 DILUPBiCG: Solving for H, Initial residual = 0.0204525, Final residual = 6.41804e-09, No Iterations 11 DILUPBiCG: Solving for H, Initial residual = 0.0172126, Final residual = 4.83197e-09, No Iterations 11 ExecutionTime = 30.71 s ClockTime = 32 s Courant Number mean: 3.50588e-05 max: 0.283714 Interface Courant Number mean: 1.07527e-05 max: 0.225319 Fourier number mean: 0.00259948 max: 1.81296 Time = 0.0003 MULES: Solving for alpha1 Liquid phase volume fraction = 0.999788 Min(alpha1) = -4.86607e-10 Max(alpha1) = 1 GAMG: Solving for p_rgh, Initial residual = 0.248811, Final residual = 8.889e-07, No Iterations 56 GAMG: Solving for p_rgh, Initial residual = 4.8953e-08, Final residual = 4.8953e-08, No Iterations 0 time step continuity errors : sum local = 7.55439e-06, global = -5.52413e-06, cumulative = -5.58166e-06 GAMG: Solving for p_rgh, Initial residual = 0.000592748, Final residual = 8.74594e-07, No Iterations 21 GAMG: Solving for p_rgh, Initial residual = 8.86839e-07, Final residual = 8.86839e-07, No Iterations 0 time step continuity errors : sum local = 7.55687e-06, global = -5.52428e-06, cumulative = -1.11059e-05 GAMG: Solving for p_rgh, Initial residual = 2.6548e-05, Final residual = 6.06255e-07, No Iterations 2 GAMGPCG: Solving for p_rgh, Initial residual = 6.06256e-07, Final residual = 9.35864e-09, No Iterations 5 time step continuity errors : sum local = 7.55429e-06, global = -5.52414e-06, cumulative = -1.66301e-05 **** ****Pressure range: 488.866 Pa ****Max velocity: 0.275259 m/s ****Phase change energy: -15.3631 W ****Volume change: -1.32635e-07 m^3/s DILUPBiCG: Solving for H, Initial residual = 0.048899, Final residual = 9.83152e-09, No Iterations 11 DILUPBiCG: Solving for H, Initial residual = 0.0988781, Final residual = 9.59266e-09, No Iterations 11 DILUPBiCG: Solving for H, Initial residual = 0.0302533, Final residual = 4.60914e-09, No Iterations 11 ExecutionTime = 49.86 s ClockTime = 51 s Courant Number mean: 0.00111238 max: 0.451657 Interface Courant Number mean: 2.23494e-05 max: 0.451657 Fourier number mean: 0.00249786 max: 1.81736 Time = 0.0004 MULES: Solving for alpha1 Liquid phase volume fraction = 0.999796 Min(alpha1) = -0.0745499 Max(alpha1) = 1.0008 GAMG: Solving for p_rgh, Initial residual = 0.0207549, Final residual = 9.63291e-07, No Iterations 41 GAMG: Solving for p_rgh, Initial residual = 3.93396e-06, Final residual = 8.35811e-07, No Iterations 8 time step continuity errors : sum local = 9.83598e-06, global = -6.79663e-06, cumulative = -2.34267e-05 GAMG: Solving for p_rgh, Initial residual = 0.00157222, Final residual = 8.93124e-07, No Iterations 24 GAMG: Solving for p_rgh, Initial residual = 9.05354e-07, Final residual = 9.05354e-07, No Iterations 0 time step continuity errors : sum local = 9.83602e-06, global = -6.79665e-06, cumulative = -3.02234e-05 GAMG: Solving for p_rgh, Initial residual = 0.000103908, Final residual = 5.83835e-07, No Iterations 3 GAMGPCG: Solving for p_rgh, Initial residual = 5.83854e-07, Final residual = 7.0116e-09, No Iterations 5 time step continuity errors : sum local = 9.83541e-06, global = -6.79658e-06, cumulative = -3.70199e-05 **** ****Pressure range: 350.415 Pa ****Max velocity: 0.218803 m/s ****Phase change energy: -18.9019 W ****Volume change: -1.63186e-07 m^3/s DILUPBiCG: Solving for H, Initial residual = 0.1639, Final residual = 9.47499e-09, No Iterations 12 DILUPBiCG: Solving for H, Initial residual = 0.272875, Final residual = 5.27952e-09, No Iterations 12 DILUPBiCG: Solving for H, Initial residual = 0.214969, Final residual = 6.45167e-09, No Iterations 12 ExecutionTime = 68.7 s ClockTime = 70 s Courant Number mean: 0.00134918 max: 0.480577 Interface Courant Number mean: 3.19172e-05 max: 0.404985 Fourier number mean: 0.00247607 max: 1.8178 Time = 0.0005 MULES: Solving for alpha1 Liquid phase volume fraction = 0.999804 Min(alpha1) = -0.251971 Max(alpha1) = 1.04115 [4] #0 Foam::error::printStack(Foam::Ostream&) at ??:? [4] #1 Foam::sigFpe::sigHandler(int) at ??:? [4] #2 in "/lib64/libc.so.6" [4] #3 Foam::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const at ??:? [4] #4 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:? [4] #5 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:? [4] #6 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:? [4] #7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:? [4] #8 [4] at ??:? [4] #9 __libc_start_main in "/lib64/libc.so.6" [4] #10 [4] at /usr/src/packages/BUILD/glibc-2.11.3/csu/../sysdeps/x86_64/elf/start.S:116 APPLICATION TERMINATED WITH THE EXIT STRING: Floating point exception (signal 8)
-
@东岳 https://www.sciencedirect.com/science/article/pii/S2352711016300309 这个不是官方的,个人开发并开源的
-
最近在用interThermalPhaseChangeFoam算气泡冷凝的验证算例,利用1.3bar水和水蒸气的参数进行模拟,不过总是发散,个人认为是密度比太大了?有没有人做过这方面的研究啊?
-
@东岳 非常感谢。dynamicCode对于后处理有用吗?
-
各位好,我用从网上下载的solver运算算例会产生dynamicCode,请问有人了解这个吗?作用是什么?谢谢!
-
@东岳 谢谢。阻力和浮力不是平衡的。我想计算一下单个气泡在重力条件下上升过程的升阻力。
-
各位大神,有没有计算气液两相流中气泡受到的升阻力的?openfoam 的Utility中有个Force,但好像不是针对多相流的。twophaseEulerFoam好像是双流体模型。像interFoam这种只解一个mixture的通过什么实现吗?在controlDict中添加function?然后重写force这个Utility?各位有什么建议吗?
-
@东岳 谢谢。这个好像是针对2.2版本的。改一下还是可以用在2.3版本。
-
@bestucan 谢谢!
-
@东岳 谢谢!低版本如2.3好像不可以?
-
如题, controlDict中的function可以在模拟结束之后调用吗?命令是什么?
-
各位大神,有人用interFoam 模拟single bubble rising in quiescent liquids问题吗?跟文章对比过吗?小弟根据参考文献,设置了算例,但不能重复所有的region。文献作者是采用Gerris软件模拟的,这个软件模拟结果好像比openfoam要准确吗?Anyway, 我附上了我的算例,恳求大神检查计算文件!
reference
[1] Sharaf, D. M., et al. "Shapes and paths of an air bubble rising in quiescent liquids." Physics of Fluids 29.12 (2017): 122104.
[0_1522156837938_pofcase.tar.gz](正在上传 100%) -
之前在朋友圈有人分享compressibleInterPhaseTransportFoam的git链接,但现在打开连接之后提示404,有人知道为什么吗?或者提供新的网址链接?如果有人下载到源文件可以发我一下吗?
-
各位foamer, 气液两相流,比如气泡在水中。但两相温度不同,如何计算两相之间heat flux 或者 heat transfer coefficient?
-
@李东岳 非常感谢。写好之后,重新编译solver会提示,在UEqn.H和pEqn.H, "no matching for operator -",请问你知道是什么原因吗
.raw格式文件如何处理
气泡升阻力的计算
气泡升阻力的计算
compressibleInterFoam error
单个气泡传热理论解
单个气泡传热理论解
喷嘴建模问题
单个气泡传热理论解
喷嘴建模问题
气泡传热benchmark
怎么用paraview画出streamlines inside bubble?
怎么用paraview画出streamlines inside bubble?
怎么用paraview画出streamlines inside bubble?
怎么用paraview画出streamlines inside bubble?
怎么用paraview画出streamlines inside bubble?
mDot 与 mDotAlpha
mDot 与 mDotAlpha
tangential surface tension
tangential surface tension
tangential surface tension
tangential surface tension
tangential surface tension
interfoam 两个气泡重心位置
openfoam如何热启动
静止气泡加热不扩张
interFoam 气泡喷射
interFoam 气泡喷射
MULES::explicitSolve求解相方程
interPhaseChangeFoam相方程
interPhaseChangeFoam相方程
mDot 与 mDotAlpha
mDot 与 mDotAlpha
interPhaseChangeFoam相方程
编译错误,求指正
编译错误,求指正
大密度比两相流发散
大密度比两相流发散
大密度比两相流发散
dynamicCode
dynamicCode
interFoam 气泡升阻力计算
interFoam 气泡升阻力计算
模拟完调用controlDict function
compressibleInterPhase
模拟完调用controlDict function
模拟完调用controlDict function
interFoam 模拟单个气泡在静水中上升问题
compressibleInterPhase
气液两相流热传递系数计算
表面张力与温度的函数