CFD中文网

@李东岳 :big_mouth: ，下半年要找工作了，不知能不能看上:crying:

0_1462072825790_14.jpg
调大此值，如果还出现警告，建议改为标准初始化

@李东岳 好的，谢谢李老师耐心解答

@wwzhao

简单的感受了一下两种赋值方式，效率确实很大提高。

#include <stdio.h> #include <string> #include <chrono> #include <iostream> using namespace std ; class Timer { private: // Type aliases to make accessing nested type easier using clock_t = std::chrono::high_resolution_clock; using second_t = std::chrono::duration<double, std::ratio<1> >; std::chrono::time_point<clock_t> m_beg; public: Timer() : m_beg(clock_t::now()) { } void reset() { m_beg = clock_t::now(); } double elapsed() const { return std::chrono::duration_cast<second_t>(clock_t::now() - m_beg).count(); } }; class SomethingBefore11 { private: int m_array[5]; public: SomethingBefore11() // zero the member array { // If we want the array to have values, we'll have to use assignment here m_array[0] = 1; m_array[1] = 2; m_array[2] = 3; m_array[3] = 4; m_array[4] = 5; } }; class SomethingAfter11 { private: int m_array[5]; public: SomethingAfter11(): m_array { 1, 2, 3, 4, 5 } //zero the member array { } }; int main(){ Timer tBefore11; SomethingBefore11 m_array_before11; std::cout << "Time elapsed: " << tBefore11.elapsed() << ‘n’; Timer tAfter11; SomethingAfter11 m_array_after11; std::cout << "Time elapsed: " << tAfter11.elapsed() << ‘n’; return 0; } [xx OFtutorial0_helloWorld]$ whatAboutThisGuy Time elapsed: 7.506e-06 Time elapsed: 1.47e-07 [xx OFtutorial0_helloWorld]$ whatAboutThisGuy Time elapsed: 8.664e-06 Time elapsed: 1.9e-07 [xx OFtutorial0_helloWorld]$ whatAboutThisGuy Time elapsed: 7.646e-06 Time elapsed: 1.89e-07

博主您好，我用v1812版本得wedge网格旋转一周生成360°的3D网格似乎很顺利，没有遇到博主您说的问题，可能差别在于axis这里：

nLayers 36; expansionRatio 1.0; sectorCoeffs { axisPt (0 0 0.0); axis (0 0 -1); angle 360; // For nLayers=1 assume symmetry so angle/2 on each side }

但是仿真的时候就有点问题了。请教一下博主，我后续的计算是用wedge网格生成的3D圆管状网格进行两相流动计算，在对称轴位置处，界面波动非常剧烈，（举例来说就像包子的褶皱），请问博主您是否有解决类似问题的经验？

2c3fb71d-127a-479f-b43f-1a41360f7552-image.png
9d8b3854-1345-4cb7-bf4c-5159ab3e4ea2-image.png

感谢。

@anubis 集合重构不求解方程，而是通过通量变化计算体积分数，最后再进行界面重构。得到的体积分数必然是间断（指界面处发生跃变）的一个场。

@cfd-china
最近还是沿袭了
http://cfd-china.com/topic/336/越多的网格得到越接近实验的解
的思路

采用符合网格无关性的网格进行模拟，模拟介质是水，fluent计算后yplus均小于200.
实验值均在0.72左右， 模拟值在0.80左右，误差在11%左右浮动，应该也满足了要求（这个还需核实）
至于为什么不采用网格少但不符合网格无关性的网格进行模拟，在原理上我还没有找到合适的解释。

但是模拟Re数范围调整了一下，以前均小于30万，现在Re达到了200万，如果继续让其满足Y+值要求，网格数会增加很多，正在想办法取舍。

@bestucan 嗯嗯，感谢大佬。检查过，check和Quallity其实都没有问题，再次感谢~

毕超的《计算流体力学有限元方法及其编程详解》这本书里面详细介绍了方程各项的离散，不过使用MATLAB编写的，我现在主要就是参考这本书的

各位大佬们，本人现在研究方向是slamming，能推荐一些书籍看看么，并且最近才接触openfoam，感觉要学的好多！

我导收了一个湘潭大学的搅拌器的英文sci。实验+cfd，用openfoam搞的。我针对算例设置这一段提出了1个不痛不痒的问题，不能一个不提啊。目前已返回我导。估计明后天就返回编辑吧？哈哈。我导没有抵制，不过提了2页问题...anyway，问题不大。:sunglasses:

@望山跑 多谢，已经解决了。。

周老师，这里我测试发现好像这种结构网格的加密方式对于三维模型可以，对于二维网格，在经过局部加密后再导入到OF里面是会报错的，错误是：--> FOAM FATAL ERROR:
face 11528151 does not belong to cell 210100. Face owner: -1 neighbour: 210099

From Foam::cellShape Foam::extrudedQuadCellShape(Foam::label, const labelList&, const faceList&, const labelList&, const labelList&, Foam::label, Foam::faceList&) in file extrudedQuadCellShape.C at line 105.

FOAM aborting
感觉好像是在z方向拉伸出来的面没有识别到具体的cell中

$\mathrm{TV}=\sum_{i}^{}|\phi_{i+1}-\phi_{i}|$

coalChemisryFoam变密度的求解器无法计算常密度的情况吗？

我确定这个可以。http://dyfluid.com/class.html#id4 我在这个算例里面做的就是一个基本常密度的，附加密度的求解器，压力方程跟coalChemistryFoam压力方程类似，。

具体coalChemisryFoam这个我不太清楚了。这求解器没玩过。

基类fvPatchField<Type>的snGrad()定义就是返回face normal gradient... 也就是$\nabla_{\perp}\phi=\nabla\phi\cdot\vec n $

@hustzzw 微信图片_20221101093137.png
就像这个图中里面圆柱区域和外面圆环柱区域

@李东岳 最后还是用了后处理的方法进行计算，又要做实验又要做模拟有点赶，用的是老板之前研究的方法几何平均粒径、几何标准差与浓度的关系，算了浓度再来得到颗粒数关系。不过布朗运动要加进去我感觉得用群体平衡方法。

这个画出来了吗