CFD中文网

我也来贴一个好东西吧：https://github.com/barbagroup/CFDPython
叫做12 Steps to Navier-Stokes
语言是python，IDE推荐anaconda
内嵌视频是youtube的，翻墙自己想办法吧。

不匿名，中航工业624所搞航空发动机的，工作9年，本科生，明年评高工，普通职工，不是领导干部。
税后年收入十万上下，在成都市郊区，房子单价五千左右。住140㎡的房子，开十万左右的车。自己觉得还行，就是每隔一段时间工作压力略大。

9年年收入10万年。
看来还是得CS

SA是Spalart-Allmaras，Reynolds-Averaged Navier-Stokes Spalart-Allmaras turbulence和Spalart-Allmaras Detached Eddy
Spalart-Allmaras采用的是标准的SA，不是高雷诺数的。DES也是标准的DES模型

在气液欧拉多相流的大涡模拟中，有没有什么模型是把liquid或者particle相的SGS viscosity和主相的SGS viscosity联系起来的？
就像RANS中的一种模型：
$\begin{equation}
\frac{\nu_{t,particle}}{\nu_{t,gas}}=(\frac{k_{particle}}{k_{gas}})^2=\frac{1}{1+\frac{t_p}{t_{fl}}}
\end{equation}$
其中$\nu_t$是湍流粘度，k是湍动能，$t_p$是particle relaxation time，$t_{fl}$是Lagarangian fluid time scale。
现在像找一个类似的模型用在LES上，但是翻了好久也没翻到 :135:

@leilei ls试了下和在外面是一样的，pwd看了下是我进singularity之前的路径。不方便贴图上来，还请见谅哈。

参数波动跟监测点位置有关，放在分离区跟放在直管中间是不一样的，另外，感觉稳态计算是监测面上的通量吧，通量守恒就可以认为收敛了吧。

solidworks导入icem时候有些地方确实有点问题，我现在先用sw把大框架画出来，自己在icem里面再编辑一下小框架，最后再切分一下，有点麻烦。你现在已经解决把，有啥好办法吗

@xpqiu
.................
25T硬盘确实花了不少钱。:sunglasses:

loop的是面网格上的所有网格。你的需求应该能实现，最简单的办法是，每个面网格中心去读一遍文档，把散点距离自己最近的那个点的速度赋给所在网格。但是会有其他问题，散点疏密是不是和面网格疏密比较匹配，也许一个面网格配多个散点或者反之。

x3d
用这个导出，导入到blender

参考文献
Demirdzic, I.; Peric, M., Space conservation law in finite volume calculations of fluid flow.

李博的桌面？

@Wayne :confused:

@wsy11 谢谢您~我现在已经可以成功构建网格啦^^

请教一下你是每一层加密都对三个方向进行加密吗，shm中能能不能只对z方向进行加密呢。@一朵苍凉

我模拟的对象是反应堆内的稳态流动传热，由于堆芯内组件数量庞大，精细建模复杂，因此采用多孔介质模型对堆芯区域进行简化。

为了验证这种等效的合理性和误差，我先对单个组件栅元进行多孔介质简化，并将多孔介质模拟结果与未简化的模型结果对比验证多孔介质模型的合理性。

单个组件栅元的截面图如下：
图片1.png

其二维轴对称模型如下图：
f510f5de-0d5e-4389-b577-ad2d3defe29c-image.png

模拟结果
(1)未采用多孔介质模型
我用了二维轴对称模型来模拟，设置固体流体内热源，入口质量流量和出口压力边界条件，流体和固体分别加载体积热源P1和P2（W/m3），流体和固体的导热系数均为常数，得到的稳态温度场和压损分布如下：
bca065f1-7cc4-4385-8113-d7e0acd6d23e-image.png
流体出口温度700℃，固体由于带内热源温度高于流体，固体最高温度为799℃，固体平均温度为756℃，组件压损为2400Pa

(2)采用多孔介质模型模拟
由于固体带内热源，流固存在较大温差，采用非热平衡模型
根据栅元几何结构指定：孔隙率γ=0.1 ，界面换热系数hfs=5000W/M2K（与未采用多孔介质模型算得的换热系数一致），界面面积密度Afs=6.71m-1

通过流速-压力实验关系确定粘性损失项系数和惯性损失项系数，只考虑x方向的流动损失，y方向系数放大100倍

别给多孔介质流体和多孔介质固体设置了内热源，分别为P1×孔隙率和P2×（1-孔隙率），初始化后检查内热源加载总量无误；流体固体的热导率输入值和上个模型一致

得到的模拟结果：
af1efa10-a95f-4e49-a69c-63b4b4357ca0-image.png
采用多孔介质非热平衡模型模拟得到的流体温度是正确的，压损也是正确的，但固体温度分布明显不对，固体域温度热点仅为725℃，平均温度为711℃，均远小于实际模型算得的固体温度值。

我自己的猜测：
多孔介质模型将单组件变为流固均匀介质后，不能模拟固体区域沿导热路径带来的固体温升？刚好我这个栅元模型中固体由于导热产生的温升非常高，如果是这样，说明多孔介质非热平衡模型不适用于我所研究的问题，因为简化后会显著低估固体温度热点。

请问我的理解是正确的么

此人在2005年左右淡出CFD-online.com，个人和他有一些私交，谈到了这个话题：为什么不活跃与论坛了？我们从您的帖子中学到了很多东西。他如下回答（翻译）：

活跃于论坛是非常耗费时间，许多问题都是重复的，越来越多的人想不劳而获的获取帮助，他们也不愿意为帮助他们的人提供一点好处。CFD不是我的业余爱好，CFD是我的职业，我需要靠CFD谋生。提供免费的服务不可持续。

@卡洛 请问你是用的哪种方法？

我错了，是$S_{ii}=0$。因为满足不可压缩流体的连续性方程。尴尬。。。:shangxue:

topo.png
我需要用ICEM画multiblock网格，几何类似于上图的模型，我知道三角形Y-block和圆形O-block，那请问这两个如何结合起来呢？在这种结构中它的edge并不共用相同的vertices
第二个问题是icem的multiblock的网格如何光顺？我尝试使用了pre-mesh smooth，但是这个功能很慢，容易卡死。我只知道可以用edit edge调整spline来改变edge形状，如果一个一个edge调是不是有点笨？