CFD中文网

是的，@hangsz 说的第二种和行业关系很大。
就是在计算机行业也是如此，化工行业的CS说是待遇20万+也是普遍。

:sunglasses:

最近遇到了类似的问题。在做高炉风口回旋区模拟，就是热风和焦炭两相流。热风从炉侧鼓入，炉顶排出；焦炭从炉顶加入，下部排出。所以，炉顶，也就是计算域的上边界，对于热风来说是出口边界，对于焦炭相来说是入口边界。在fluent中设置边界条件的时候遇到了问题。看楼上李老师说，速度可以设置，意思是不是在上部边界选择速度入口，然后热风相给负速度，焦炭相给正速度？

补充一下，不知道这些图中的这些数据点是如何提取的，谢谢大家

页数按照pdf的页数，不是书上的页数

上册：

https://www.jianguoyun.com/p/DVcobNoQ9s3ZBhjejV8

73页【直线运动的速度】：导数的直观描述，CFD方程这面全是导数，大部分为偏导数，后面会介绍 75页、76页【公式1-4】 138页【公式3-3】：泰勒公式，CFD很多方程都是基于泰勒公式推导出来的

下册

https://www.jianguoyun.com/p/DQnanmcQ9s3ZBhjijV8

23页：流量计算公式 71-75页：偏导数，CFD这面都是偏导，重点看一下，理解一下 77页提到了个拉普拉斯方程，随便扫一眼 78页全微分，跟CFD的物质导数关联 79页3-3 81页有个例子 112页介绍了梯度，但是感觉描述的不是很简单，可以看看 141-144页 重积分，为之后理解有限体积法做准备 224-225页 曲面积分的定义，为之后理解有限体积法做准备。其中225页有个图，比较形象 230-231页，看高数中对通量的定义 243-245页，高数中明确的对速度散度的定义

这个主要看cfd之道公众号流沙的教程就可以学会了（哈哈哈哈，你是我拿到注册码之后回复的第一个人）

@random_ran

我给他代码高亮的时候显示成斜体了，现在改好了

此问题同类于湍流动能耗散率的定义。湍流动能耗散率的定义为：
\begin{equation}
\varepsilon=2\nu_t\overline{S_{ij}S_{ij}}
\end{equation}
经过各向同性假定后有：
\begin{equation}
\varepsilon=\nu_t\overline{\frac{\p u_i'}{\p x_k}\frac{\p u_i'}{\p x_k}}
\end{equation}

@bestucan 在 关于fluentmeshing画网格 中说：

@hitsc30 我没搞过边界层 y+ 什么的。
但是，对于库郎数大于1的状况，fluent 确实比 OpenFOAM 能算，但是能算是怎么算的就不知道了。

如果是高速流动，边界处很难出现垂直界面的速度，大纵横比就没问题，建筑上吹的风，有个回流什么的就……

当然要直接解决就是减小时间步，特别小的那种，保证在垂直界面方向一个时间步也跨不过一个网格。但是估计会很慢。

snappyhexmesh确实玩的可行性比较高，我画了一个1300万的网格，运行了4h，质量还是不太行。目前对网格有点麻木了，谢谢指导！

0_1460680809776_捕获.JPG

我在昨天发现了这个错误，怀疑也是内存出的问题。仔细回想了一下，昨天在程序里面写了几行：

scalar* zeta = new double [2*n];

后面忘了删除。昨晚上计算的时候三个算例全部killed：

diagonal: Solving for M11, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for M12, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for M13, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for M14, Initial residual = 0, Final residual = 0, No Iterations 0 PIMPLE: iteration 1 GAMG: Solving for p, Initial residual = 0.0004130904583986163, Final residual = 2.405751520795531e-06, No Iterations 4 time step continuity errors : sum local = 8.330982304837603e-10, global = 1.267524541545882e-12, cumulative = -2.029919160042075e-06 Killed

无任何消息。然后我在源代码中将这个new出来的数组delete之后，就不会出错了。

另外，我测试了一下，不删除数组计算的时候内存会一直上升。因此，动态数组记得释放！

@Wayne :confused:

@bestucan 在别人电脑上是正常的 ，很奇怪，我最后重装电脑解决了

@him 谢谢解答，另外一般商业软件的定常求解器是用的时间推进来求解的吗？如果是时间推进的定常求解器，这种不稳定的振荡可否理解为流场随时间的变化？

@d-benjamin 这个是正则表达式，是用来匹配字符串的。"."代表任意字符，"*"代表出现多次（包括0次），合起来".*"代表匹配任意字符。

"epsilon.*"可以匹配以epsilon开头的任意字符（包括本身），如epsilon，epsilonFinal等。其中epsilonFinal是用于PIMPLE算法最后一个迭代循环的控制参数。

以上规则同样适用于"U.*"和"k.*"。

如果使用standard KE，计算二维流场，在某个高雷诺数下把Cd, Cl, St都算对，那么说明碰对了，意思就是好运。这并不能说明standard KE可以算对漩涡脱离。

这个是理论大坑啊。

如果压力基的求解器，求解的是速度和压力场，那么得用连续性方程来求压力。
如果是密度基的，连续性可以求密度，然后用状态方程求压力。

@徐小双 应该是没有删除干净的原因。

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

0_1537523541221_IMG_20180921_173609.jpg

@东岳 好的，谢谢李东岳老师，在您的建议下我继续尝试一下，再次表示感谢！