@wwzhao 好没关系 谢谢您

怎样在主程序中添加correct( ); 来在更新类中的　rho1()rho2() p()

暂时还没研究过这个类:sad: ....

@gzh3369

如果公式是这样的话：$x=\frac{\sum V_i x_i}{\sum V_i}$。可以这样：

forAll(x, celli)//对所有的体x进行计算 { x[celli] = x[celli] *mesh.V()[celli];//此处即为每个网格的x*网格体积 } scalar sumx = gSum(x);//对所有的x*网格体积进行加和 scalar V = gSum(mesh.V());//计算总网格体积 scalar Final = sumx/V;//计算体均湍流动能耗散率 Info << runTime.timeName() <<" volume average x is "<<x<<endl;//输出结果

楼主好人 w-y10@foxmail.com

据我所知，在面向对象由于添加了RTTI机制、构造函数机制会使得程序相对于使用一个函数变得稍微慢。当然这个慢我们基本看不到。

在实际使用中，在我改写一个非常复杂的函数为一个类的时候，反而程序变得更快了。这是什么原因造成的？

比如我执行这个外挂函数：

... A(a, b, c, d, e, f, g, h, e, j, k, m, n, z); ...

反而不如把他改写为类快：

... A B(a, b, c, d, e, f, g, h, e, j, k, m, n, z); B.calculateA(); ...

目前的认知存在于类要比公式慢一点点。不过这个认知被自身的测试推翻了。

Any ideas?

:sunglasses:

@金石为开 本质上是 hashTable 的使用。在基类里调用宏函数 declareRunTimeSelectionTable 创建了一个 hashTable，其 key 的类型为 “word”，value 类型为函数指针，这个函数指针指向一个返回派生类临时对象的智能指针（autoPtr）。

在派生类里，每次调用 addToRunTimeSelectionTable 的时候，相当于往 hashTable 里添加一组元素。

在求解器里，创建的是基类的智能指针，这个指针指向基类中定义的 New 函数的返回值。New 的作用是选择具体的派生类，这里涉及到C++的一个特性，即基类指针可以指向派生类。New 函数根据你在算例里提供的模型的名字（typeName），去 hashTable 里匹配，然后将匹配到的派生类的对象返回。这样，就实现了模型的选择。

具体参考：http://xiaopingqiu.github.io/2016/03/12/RTS1/

@cfd-china 多谢啦，我试试！

@Aeronastro

在tutorial中，通常如此：

outlet { type inletOutlet; inletValue uniform (0 0 0); value uniform (20 0 0); }

这可以有效地防止回流
如果直接给定zeroGradient，可以出现回流;

可以拿一个回流特别明显的算例来试试，比如一个涡旋伸展后方的流动。进口速度(1,0,0)，出口的inletValue是(-1,0,0)。看是否不可压缩算例能否守恒。

:kiss:

@mark 信息量太少，我无法判断问题原因。不过你可以根据截图中的back trace去找错误。

nparticles 是不是就是那个颗粒数量？

0_1462776631144_QQ截图20160509145537.png

dimensions [0 1 -1 0 0 0 0];

internalField uniform (0 0 0);

boundaryField
{

INLET
{
　　　 type fixedValue;
　
　　　value uniform (0 0 0);
}

WALL { type fixedValue; value uniform (0 0 0); } OUTLET { type pressureInletOutletVelocity; inletValue uniform (0 0 0); value uniform (0 0 0); }

}

// ************************************************************************* //
｀｀｀
这样写对不对？

知道了。内部场和边界场不是分开的，所以用等号赋值的时候也就更新了边界场。

@mark

那就可以不读取密度:sunglasses:

谢谢回答。现在知道把模板参数 typedef 是为了可以输出模板参数的类型。参见C++ traits

@mark 初始条件不同？用setFields可以

@wwzhao 既然大小圆柱要分开画，那么它们之间重叠的边界面应该怎么处理，不是同一个边界面吗？

@mark 能否提供测试代码和算例？

@mark 这个是个构造函数，alpha10进去注册的场，其值为alpha1。:expressionless:

@中国轰炸机飞机斯曼

Try this:

scalar xxx=readScalar(RASProperties.kOmegaSSTCoeffs.lookup("Prt"));

向x轴正方向运动，看来我设置对了:sunglasses: thanks

能否留下联系方式，我也用到interphasechangefoam

@cfd-china

42万的网格。

我觉得应该就是全部由三角形面单元组成的surface。

@cfd-china 嗯嗯，好像是没什么影响，谢谢哈:cheeky:

@cfd-china 没有必要手动进行分块，我只要保证有一个processor的INLET边界是完整的就好了，其他的怎么分块比较随意，实际手动分块的话也可以，我为了简单，就用了simple了。

记得之前好像有关于手动分块的方法的帖子，我找找看:happy:

之前那个问题解决后又出现0_1461572346439_upload-5c98c499-d439-44d1-96ce-31e1b96bbd8a 这个错误，求教submap是什么文件

@cfd-china 十分感谢！:happy:

@wwzhao 嗯嗯是的，我因为要在主函数中调用一个函数，而这个函数里要用到一些量，所以我把这些要用到的量定义成了全局变量，这样的话在执行完这个函数的时候，那些变量还在，供主函数继续操作。而别的量我已经在createField.H中定义了。

@lhzhu

在计算大规模问题时，采用多个GPU进行计算，在不同的GPU之间进行数据传输是个问题。不过现在好像已经有这个技术了。

另外还有GPU和CPU数据传输速率的限制，不过这个在Maxwell架构中似乎已经不成问题了。

不管怎么说，GPU在高性能计算领域普及开来还有很长一段路需要走。

@散漫守望2016 非常感谢，看到了，正在尝试联系帖子主人。

@mark

这三个参数是在system/controlDict中打开adjustTimeStep（即自适应时间步）的情况下才起作用的。

求解器通过这三个参数来控制时间步长。这三个参数根据不同问题可能需要设置不同的值，我一般设maxCo 25; maxAlphaCo 15; maxDeltaT 0.01;，你可以参考一下。

@mark alpha1_ 代表VOF模型中的体积分数，是一个volScalarField类型的对象。

@Wayne 多谢多谢，研究一下。有啥不懂的再向大神请教吧。

@张晓宇 在 关于toalpressure边界条件的一点疑问 中说：

在C++ documentation里面看的toalpressure的使用，有一个参数gama是气体热力学的性质参数，而且是必须要给的数值。这个意思是不是就是说toalpressure只能用于气体呢？

总压方程，例如$p_0 = p + \frac{1}{2}u^2$来自伯努利方程。对于马赫数比较小的情况，适应于上述方程。当马赫数较大的时候，接近音速或者超音速的时候。则变为可压伯努利方程。这些在空气动力学书里面一般都可以找到。一般气体可以看错可压的，大部分情况下液体一般看做不可压缩的。并且，载推到可压缩伯努利方程的时候用到了理想气体状态方程。所以，totalPressure的可压缩部分用于气体比较合适，但是不能说不能用于液体。

这个p0是相对压力还是绝对压力呢？

公式中的p0为总压

离心泵进口条件有全压的设置，这个全压能用toalpressure么？

什么是全压？总压？可以把进口设置为总压。

@Wayne 在openfoam里面，　访问未定义的内存是什么意思呢？

@cfd-china 非常感谢！

@xpqiu 又学到了 非常感谢！

@xpqiu 十分感谢！我会仔细根据您提的意见学习的，谢谢！

谢谢...

感谢管理员，

这个minimal example 不写文件确实是因为我忘了加上runTime.write():cold_sweat: 。它并不能重现我原来solver里面的问题。我再找找原因。

@李东岳 谢谢版主:laughing:

感谢，我看了下，UserGuide里面U28说的是把scale关掉，这会导致所有的箭头长度和箭头头部大小都一样。目前从ParaView的控制面板上还找不到保持箭头头部大小不变，而箭头长度随速度大小scale的选项。

一般的工具（Matlab, matplotlib, Paraview）里面， quiver图 里面的箭头大小都跟着箭头长度一起变，所以在矢量很小时候，效果很糟糕。如下图就是Matlab里面quiver画出来的，
0_1461126738148_upload-25983765-853b-4a0a-b72b-a0ffce9f781b

不过后来找到一个比较好的办法，在latex里面使用pgfplot宏包可以解决。
比如如下代码(来自这里)

\documentclass{standalone} \usepackage{tikz} \usepackage{pgfplots} \begin{document} \begin{tikzpicture} \begin{axis}[domain=-3:3, view={0}{90}] \addplot3[blue, quiver={u={1}, v={(x-y)}, scale arrows=0.15}, -stealth,samples=20] {0}; \end{axis} \end{tikzpicture} \end{document}

可以产生画出：
alt text
可以看到箭头大小不随箭头长度改变。

@李东岳 这一项和interpolationSchemes改为linear就可以运行，用MUSCL和venLeer就不可以

@聆听whatever

Hi,

可以看看这个：链接文本。欢迎补充交流。:baby:

我的并行出问题了，不并行倒是可以，我想具体学一学怎么用Debug模式来调试，能否举个例子教一下我呢？谢谢各位了~

@lhzhu 在 如何拷贝/创建 volVectorField 在每一个时间步文件夹中 中说：

只要0目录里面有一个场文件，decomposePar 都会在processorX里面生成子区域的场文件。是不是你写的0/HAG_BC 的file header 里面的object 不是 HAG_BC ?

decomposePar 后 在“processorX” 中的 "0" 目录里有这个HAG_BC, 问题是下一时间步的时候没法创建这个HAG_BC啊，
所以会有那个error.
不知道怎么解决

@李东岳 在 有关侵入边界法。 中说：

Hi，

OpenFOAM-extend里面有IBM方法，Extend和官方的OpenFOAM侧重不同，几个我了解到的不同主要是IBM，粘弹流体，动网格表面跟踪。

还有可压缩流动高分辨率格式。https://github.com/Unofficial-Extend-Project-Mirror/Unofficial-Extend-Project-Mirror.github.io/wiki

http://www.tfd.chalmers.se/~hani/kurser/OS_CFD/HrvojeJasak/DiscretisationBestPractice.pdf