CFD中文网

@lgw https://github.com/gerlero/openfoam-app/tree/v1.8.4

这个用interFoam+动网格不就可以么

@bestucan 运行第一行，物理CPU直接就是32

@李东岳 非常感谢。
又读了几篇类似的文献，真是用时均处理的仿真模型找不到，高雷诺数下的这种流动都是在用实验测试了，但是我需要计算一个比较准确的结果，合作单位的同志希望能找到统一的规律，形成经验公式类似的结论，然后用作产品开发:crying: 任重道远啊，有个师兄推荐用DES试试，目前在测试，希望别掉进坑了。

@李东岳 好的~谢谢东岳老师:xinxin:

@下里巴人 最近稍微试了一下均匀流下粗糙网格的3组试验，一组时间步长给定为旋转2度的时间，一组为0.2度，一组Fluent自调整。因为时间原因仅做了一小段时间内的对比（未达到稳定周期性变化），关注的两个物理量（一个轴向力，一个转矩）在三组试验中趋势几乎一致，数值误差小于3%.另外对一种非均匀流下的精细网格进行了同样的试验（同样时间较短），趋势几乎一致，但不如均匀流下趋势符合得那么好，误差在9%以内。
受时间限制，对比不严谨，仅作为前期考察

@bestucan 感谢，follow path 和follow data 都很好用

@李东岳 感谢李老师。我去学习下:140:

请问各位大佬，对于可压缩流大涡模拟的总湍动能如何定义呢？

对于不可压缩，大涡模型采用的亚格子模型采用的雷诺时均的形式，因此可以用雷诺时均计算大涡湍动能和亚格子湍动能，两者的和就是总湍动能。
$K=0.5*(\overline{uu}-\bar{u}\bar{u}+\overline{vv}-\bar{v}\bar{v}+\overline{ww}-\bar{w}\bar{w})+\bar{k}_{sgs}$ 但是对于可压缩流的大涡模拟，亚格子模型采用的为法夫尔时均形式，此时如何定义总湍动能呢？我的理解是
$\widetilde{K}=0.5*(\overline{\rho uu}-\bar{\rho}\widetilde{u}\widetilde{u}+\overline{\rho vv}-\bar{\rho}\widetilde{v}\widetilde{v}+\overline{\rho ww}-\bar{\rho}\widetilde{w}\widetilde{w})+\widetilde{k}_{sgs}$
请问各位大佬，我的理解对么

@李东岳 除了phase的图像，其他压力图，密度图等都没有问题。0_1504141697200_ea5fe613-9544-43ae-acce-f76b5a6f726e-image.png 这是密度图。我在网上有人提到会不会是由于相与相之间的界限不明显，而产生这个问题。老师，请问您觉得这种说法可靠么？谢谢

@Samuel-Tu :
谢兄之玉言，弟牢记

然兄亦暗示：
有的参数需调整，以保证物理意义一致。

另外，如兄及各位大学士不弃，望再赐一教：
OF能自动算雷诺数吗？

感激不尽

有的时候icem就是会报错，只能没做一些操作之后多保存，然后再进行。。

想问下这个职位还招人吗，谢谢！

@李东岳

李老师您好。

目前出口是采用如下方程
3.png

这里Uc尝试过：

使用上下层速度的平均值1.5。 在X方向上采集最靠近边界的三个点的速度值进行平均，作为Uc。

如果采用2的话，Uc是一个随着Y坐标以及时间变换的函数Uc(y,t)

Phi是压力泊松方程构造的中间量，这个查了好多文章都说的很模糊，我也就直接用了零梯度。

我用openfoam也出现这个问题，有种说法是不可压假设导致的，我还没试过

OpenFOAM v2006中，adiabaticPerfectFluid状态方程如下：

$\rho (p)=\rho_{ref}(\frac{p+B}{p_{ref}+B})^{\frac{1}{n_{T}}}$

根据熵的计算公式：

$dS=c_{p}\frac{dT}{T}-(\frac{\partial v}{\partial T})_{p}dp$

求解器中将这一公式分为两部分进行积分计算，其中第一部分$c_{p}\frac{dT}{T}$ 的积分计算在thermophysicalModels/specie/thermo/hConst/hConstThermoI.H中实现，即$c_{p}*log(\frac{T}{T_{std}})$
但第二部分的积分，由于公式是绝热状态的，没有温度作为参数，想请教一下官方求解器中为什么会进行以下处理：

通过源代码公式的求导反推，官方求解器在计算第二部分积分时，在状态方程的分母上直接添加了温度T，如下式：
CodeCogsEqn (1).gif
公式源代码详见adiabaticPerfectFluidI.H文件的S函数：

template<class Specie> inline Foam::scalar Foam::adiabaticperfectFluid<specie>::S ( scalar p, scalar T ) const { scalar n=1-1.0/gamma_; return -pow(p0_+B_,1.0/gamma_)*(pow((p+B_),n)-pow((Pstd+B_),n)) /(rho0_*T*n); }

除上述绝热状态方程外，线性变化的状态方程linear也采用了同样的处理方式，想知道为什么可以这样处理？
另外一个问题是，如果采用上述处理之后，对于气体的$c_{p}-c_{v}$，根据公式:
CodeCogsEqn (2).gif 该如何处理？

ICEM画二维网格导出时报错：
WARNING: Mesh has uncovered edges. ANSYS Fluent needs a complete boundary (lines in 2D) or it will give a variety of errors and not read in the mesh! If this was 2D Hexa, perhaps your edges are not associated with perimeter curves

模型描述：二维管道流体的测试网格，使用ICEM导入模型后：
1、划分块
2、在preMesh 里 使用Edge Paras进行边界层的绘制，
3、preMesh
4、在preMesh里Convert to Unstruct Mesh
导出网格选二维，其他未改变，出现上述报错
新手从B站学了一点划分网格操作就尝试了这一个，不知道哪里出错了，还希望前辈们解答

:xinxin2:

@程迪 666啊老铁，听君一席话，胜读百篇论文:146: :146: :146: :146: :146:

@李东岳 李老师有没有兴趣看下风工程这块的有个商业软件WT，里面是速度-压力耦合求解器，真的非常快非常快。我看了下，软件底层是从PHOENICS来的，然后搞了一套叫MIGAL的算法。

现在CFD用的SIMPLE这种segregated solver确实太慢了，城市这种尺度下的流动本身大概率又比较简单，不可压，也没有化学反应之类的，可以暂时不考虑温度？可能会有颗粒流的问题，但是不是也可以先只考虑流动对颗粒的影响。这种情况下感觉流场本身的求解其实就看求解器的效率了。

我搜了下，OpenFOAM这边其实已经有一些coupled solver 的工作，不知道成不成熟。