CFD中文网

网上实在搜不到
作者: J.O. Hinze (Author)
出版社: McGraw-Hill Inc.,US; 2nd Revised edition edition

ISBN: 9780070290372 出版时间: 1975-05

邮箱 yuan_neu@163.com

@Junren-Hou 在 ICEM划分网格 中说：

可以试试check block，有些时候其实没有问题，就它软件自己的毛病，然后block里最后有一个check block，可以让他自己检查一下。里面有好多选项，都试试，什么check/fixed block

回忆了一下，我的旋转区域用ICEM画的应该没有问题。出问题的应该是我画的外流场区域，一方面追求地面边界层，致使网格尺度过小，一方面为了平衡网格数量，过渡不均匀。猜测应该是外流区域的地面边界层设置导致了相关报错的产生。

此问题同类于湍流动能耗散率的定义。湍流动能耗散率的定义为：
\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}

@东岳 好的，麻烦您啦！我在一点一点仔细研究下。如果能研究明白，万一以后也有人问，我再做个帖子出来♪(^∇^*)

想问下，有人做过这些网格划分和几何建模工具的对比吗？也想了解了解

小弟现在在模拟往一个封闭罐体里充气的过程，然后已知当地是三级风（风速在3.4到5.4 之间），哪位知道有大概的公式（最好能备注文献来源）根据风速算对流换热系数，与空气对流换热的物理结构是中间一段直线加两端四分之一椭圆弧。

页数按照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页，高数中明确的对速度散度的定义

在Drag Report Definition里，create Output Parameter是什么意思，一旦勾选，不能自动保存了，帮助里解释太少了，谢谢!!!

@zhe 谢谢，我按照方法排查尝试了一下，发现并不是粒子的加入影响了液体速度场，二而是MPPICInterFoam对于setFields对于输入速度场就有问题。下图一个是MPPICInterFoam计算的，一个是interFoam计算的，同样设置的5m/s液体下落速度，但结果截然不同。@李东岳 也请麻烦李老师能否帮忙看一下，谢谢！
2a82db74-6e61-43e1-a687-5a6ca08d9055-image.png
7751b300-7c97-4692-97a7-48f19f694107-image.png
96ee3830-8f39-4ffe-8986-3988a9d011df-image.png

@李东岳 果然跑起来了！

如果气含率比较大呢，比如0.3甚至更大呢，压力30 kPa和300 kPa感觉如果考虑压缩性体积差别很大，气含率这些差别应该也很大？

我觉得你更应该从实验或者文献上获取一下信息，因为有的需要考虑，有的不需要考虑。不太普适性。

我那个考虑了PBE，比较麻烦，你可以不考虑PBE，直接植入一个跟压力相关的线性粒径函数，比如$d=f(p)$，非常简单

不过多谢引用我那个文章，现在引用数低了点拉低了h因子，帮我引用我很高兴 :papa:

湍流普朗特数.png
请问有同志接触过S-A模型中将湍流普朗特数改到0.4以下的情况或者有提及此操作的文献吗？烦请赐教，谢谢！

@两月三年 嗯 有点了解了

@coolhhh 应该是平均值

大家好，我用 icoUncoupledKinematicParcelFoam求解器模拟静电除尘器中颗粒荷电运动轨迹，颗粒考虑曳力、重力和电场力，计算电场力需要电荷量q（通过方程计算，是动态变化的值）和场强E(矢量场，给定）两个参数，我在submodels/Kinematic/particleforces路径下添加了一个电场力子模型，怎么能做到q的实时读取和电场力的实时更新呢，求大神指点。
0.png

大都忽略重力源项，哪位大神给解释下为什么忽略重力源项

重力对流体的作用就是引起水动力学压力，如果是密闭的流体，其不会引起流动变化。因此可以忽略重力项，这样压力边界条件也可以好给出。

一般考虑重力的主要是重力为一种重要的体积力，如密度不均引起的浮力，其中密度不均进一步由温度导致。

还有一种情况要考虑重力是模拟流体的自由落体。

其他情况可能跟具体的应用有关，如旋转系、海洋热风那种普遍都添加重力。

如果不忽略，怎样恰当的处理？

如果考虑重力，通常求解的压力不是先前的压力p，而是压力p减去水力学压力$\rho g h$。从算法层面理解，重力引致下落，压力对这种效果抵消，采用这种压力计算的时候，可以省去这种来回下落-抵消的求解过程。求解出$p-\rho gh$之后，真实的压力即为$p-\rho gh+\rho gh$

@ibelief 摇头电风扇的话就涉及到三个域了，一个是流场静止域，两外两个旋转域，其中一个旋转域是支撑轴的旋转，另外一个是风扇叶片的定轴旋转，这就涉及到风扇转轴跟随着支撑轴旋转（风扇摇头）了，其中第二个应该涉及到欧拉角的定义了，这个该怎么设置没搞明白。
另外一个问题就是风扇叶片旋转区域网格需要跟随支撑轴摇头，也就是自身旋转的同时要绕另外一个轴旋转。
还有就是，如果扇叶定轴旋转可以用MRF实现，摇头运动通过滑移网格实现，这就涉及到流场中一个是非定常运动区域，另外一个是设置为准定常区域，这种不知道在物理逻辑上能否说通。如果不能说通的话，就需两个旋转域均设置为滑移网格计算，而当叶片转速较高的话，就会限制时间步长过小，带来计算代价过高的问题。
万分感谢！

@ann 在 chemkin 中说：

在网上找到的C12H26化学反应机理文件，不管是直接运行还是通过chemkinToFoam转换都显示类似格式错误的问题

因为OpenFOAM里chemisrtyReader和chemkinToFoam只会读transportProperties这个文件里的数据，类似下面这种。用在sutherland公式里。

".*" { transport { As 1.512e-06; Ts 120.; } } "H2" { transport { As 6.362e-07; Ts 72.; } } "CO2" { transport { As 1.572e-06; Ts 240.; } }

@ann 在 chemkin 中说：

我在网上找到的trans文件都长类似的样子

chemkin软件的三个文件之一就是这个六列分子特性数据的trans.dat，这才是正确的。只是OpenFOAM因为用sutherland公式就不需要这个文件了。sutherland公式需要的是上面的As和Ts。
使用OpenFOAM自带的chemkinToFoam命令：

chemkinToFoam chem.inp therm.dat transportProperties reactions.dat thermo.dat