CFD中文网

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

@东岳 让孩子赢在起跑线上:high:

@红豆沙 并没有。。。我们这边很神奇，几个学校的cluster是共享的，有各种CPU，5118，6142，EPYC 7501什么乱七八糟的，还有V100加速卡，然而我自己实验室给的一台破E5450，自己买了台Y540 打游戏。。。

原实验论文里面把这个case称作‘subturbulent’，雷诺数是6020，T0 = 297.15 K,
P0 = 0.0927 MPa, D = 0.0274 m，V0=3.73 m/s，文中好像有讲到流动有趋于层流的趋势。

@李东岳 好的~谢谢东岳老师:xinxin:

对简单的长方体进行网格划分，再检查网格质量时发现在长方体的周围有很多有问题的face产生，也不知原因，很疑惑，求大神指点，万分感谢！

0T5OC4QVS4OWJO998JSUAU7.png

@东岳 谢谢东岳老师的关心，后来还是有点进展，不过进度不是很大，现在所走的路就是修改求解器，文章还没有出来。等出来之后我会详细说下这个问题。

据小弟了解：湍流数值模拟方法包括DNS，RANS，DES和LES等，其基本控制方程都是N-S方程耦合自身封闭方程。 请问：
1.对于数值模拟中SPH方法，和上述方法的归类情况如何？
2.还有那些常见的数值方法？

谢谢各位！

icoFoam twoPhaseEulerFoam

@东岳 好的谢谢！
我想问下圆周沿线要怎么画？:chouchou: ，另外流向的平均就是流向沿线每个点的u+相加起来然后除以点数吗？
那y方向上的位置要选多少个？谢谢！

Navier-Stokes 方程中的对流项写作 $\nabla \cdot (\textbf{u}\phi)$，扩散项写作 $\nabla \cdot (\nu \nabla \phi)$。

对流表示流体微团从空间某点运动到另一点的过程，而扩散则表示流体物理性质浓度变化的过程。

对动量方程而言 $\phi=\textbf{u}$，因此对流项变为 $\nabla \cdot (\textbf{u} \textbf{u})$，为非线性项。一些算法将其中一个速度作为已知量，另一个作为未知量离散求解。可参考icoFoam解析。

@samuel-tu
前辈好，请问交界面网格畸形的问题您是否解决了？我也遇到了类似的问题:135:

@黄昏晓为谁破 好的感谢前辈回复

不清楚商软那面怎么回事，本身我也没太关注这个，我查了一下，OpenFOAM本身这面的streamline看起来是你说的patiline：

This function object generates streamline data by sampling a set of user-specified fields along a particle track, transported by a user-specified velocity field.

后处理软件paraview的streamline是如何定义的还没找到。

vim保存窗口：

:mksession ~/.session.vim vim -S ~/.session.vim

@benf 不客气不客气，欢迎引用我们组的论文~这个模型有外国组已经验证过了，所以如果你有兴趣编应该可以原样复现

格式不对噻

@李东岳 ansys18版本 linux
链接：https://pan.baidu.com/s/1QpXN_GhW5tCv1rNJrlBUYA
提取码：1234