CFD中文网

大家好，问一个很弱智的问题。。。用了openfoam的sampling surface，输出格式为raw。请问用什么软件处理呢？

@cav 在 气泡升阻力的计算 中说：

@东岳 谢谢李老师回复！在询问原作者。

您好，请问原作者有说怎么计算气泡表面的粘性阻力吗？我最近也碰到了这个问题:zoule:

最近在模拟有传热情况下额气泡上升问题，求解器为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

@CYW 谢谢。在interFoam里加入过能量方程，但是在加物理量比如普康特数啥的出现了点问题。。。我对照着compressibleInterFoam研究研究，，有问题再向你请教。不好意思没帮到你什么忙

@mohui 谢谢！

可以算一下compressibleInterFoam下的传热:duang:

@东岳
谢谢！这下明白了！

@东岳
在interfaceProperties.C中编程如下
fvc::interpolate(fvc::grad(sigma) - fvc::grad(alpha)*(fvc::grad(alpha)&fvc::grad (sigma))*mag(fvc::grad (sigma))) & mesh.Sf()
将此项添加到速度和压力方程。其中当计算压力项时发散了。

同问，重心的位置怎么求，，，

@cyw 你好，请问你是在哪个版本中添加的传热模型呢？
网上只找到了of230的方法。
我看过模拟加热气泡的文献，晚几天把名字发给你可以参考下。
话说，方便细聊么？我最近也在弄这方面的东西
vx：zhigu4610

你的网格尺寸太大了，缩小1/10或者1/100试试

phi是用来调用判断速度方向的： surfaceScalarField phiBD(upwind<scalar>(psi.mesh(), phi).flux(psi));

@MrZZQi 感谢分享 ！:xiexie: 前几天刚跟海东说完更新一下这个内容，顺便参考一下你的内容

@东岳 Thanks!

@东岳 我对这个期刊也不熟。只是在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)

@东岳 您好，我想知道#calc 可以使用的标准函数包括哪些呢？哪里能查得到？

interFoam是直接模拟，不能模化升力。
双流体模型是模型化后的模型，可以模化升力。
因此

升力的方程表达式中出现气相与液相速度，这两者速度与interFoam控制方程中的速度有什么关联呢？

二者毫无关系

@东岳 谢谢。这个好像是针对2.2版本的。改一下还是可以用在2.3版本。

@AppleKiller 我知道一个开放数据集，但是里边单气泡这部分没开放，我现在密码过期了看不到是u否开放了，可以问问他们http://141.44.132.124/spp1740/index.html要是开放了也@我一下

@yuan_neu https://github.com/OpenFOAM/OpenFOAM-dev/tree/master/tutorials/multiphase/compressibleInterFoam/laminar/climbingRod 最新版本代码已经把climbingRod算例换位置了。所以曾经的连接失效了。5版本就有着个求解器compressibleInterFoam，不过没这个算例climbingRod