CFD中文网

@turbulence 02525037-654f-4a52-bc51-dcef3cc03688-image.png

描述热扩散的无量纲物理量热扩散因子究竟是什么？有无针对不同情形下的具体公式？
1.png 2.png 3.png

用的windAroundBuildings算例，840万网格的二维模型，网格质量检查过了，计算域是14000*4000米。epsilonInlet 0.0003; // Cmu^0.75 * k^1.5 / L ; L =10

dimensions [0 2 -3 0 0 0 0];

internalField uniform $epsilonInlet;

boundaryField
{
INLET
{
type fixedValue;
value uniform $epsilonInlet;
}

OUTLET { type zeroGradient; } WALL { type epsilonWallFunction; value uniform $epsilonInlet; } frontAndBackPlanes { type empty; } #includeEtc "caseDicts/setConstraintTypes"

}
kInlet 1.5; // approx k = 1.5*(I*U)^2 ; I = 0.1

dimensions [0 2 -2 0 0 0 0];

internalField uniform $kInlet;

boundaryField
{
INLET
{
type fixedValue;
value uniform $kInlet;
}
OUTLET
{
type zeroGradient;
}

WALL { type kqRWallFunction; value uniform $kInlet; } frontAndBackPlanes { type empty; } #includeEtc "caseDicts/setConstraintTypes"

}
nut
internalField uniform 0;

boundaryField
{
INLET
{
type calculated;
value uniform 0;
}

OUTLET { type calculated; value uniform 0; } WALL { type nutUSpaldingWallFunction; value uniform 0; } frontAndBackPlanes { type empty; } #includeEtc "caseDicts/setConstraintTypes"

}
p
internalField uniform 0;

boundaryField
{
INLET
{
type zeroGradient;
}
OUTLET
{
type fixedValue;
value uniform 0;
p0 uniform 0;
}
WALL
{
type zeroGradient;
}
frontAndBackPlanes
{
type empty;
}

#includeEtc "caseDicts/setConstraintTypes"

}
Uinlet (10 0 0);

dimensions [0 1 -1 0 0 0 0];

internalField uniform (10 0 0);

boundaryField
{
INLET
{
type fixedValue;
value uniform $Uinlet;
}
OUTLET
{
type zeroGradient;
}

WALL { type noSlip; } frontAndBackPlanes { type empty; } #includeEtc "caseDicts/setConstraintTypes"

}
时间步长试过1 0.1 0.01 其他的都没动

@chaoscfd 我用了 但还是mesh.findcell最快

按说，那个 pdf 的大部分勘误记录都会跳转到这里
https://cfd-china.com/topic/2338/无痛苦n-s方程笔记-勘误/38

“错”和“误”不一样，（所以某些理论上无法消除、工具条件限制带来的与事实不符的差距叫“误差”不叫“错差”）（我对你措辞没意见，只是看你挺有互助热情，也“互助”你一下，说给你听听，参考参考。）

可以直接说出笔误位置，细微的笔误难以发现（容易发现也不会误了）。比如：

"接下来定义温度 T 在 (x1, y1, z1) 点的物质导数为移动的无穷小微团通过 (x1, y1, z1) 点的时候"
这句第二个坐标（或者第二个（x1, y1, z1））标号应该是(x2, y2, z2)

我看你之前也发过类似的帖子，可以看完了一遍书，把发现的笔误统一发一下，也方便修改。

欢迎继续:146:

libs ("libmyfvOptions.so"); 这句不能像一些官方提供的案例那样放在fvoptions文件里，得放在controldict文件下
Screenshot from 2021-12-03 17-31-56.png

@caijinjin 是计算停在这里还是终端会停止，方便的话可以交流一下吗，我微信15098932699

topoSet有一些比如寻找点最近的面，寻找距离stl最近的面之类的，你试试

各位安好，
想请教一下，对于离心叶轮来说，他不同叶高的截面不是一个圆柱，而是一个曲面。这种情况下，怎么提取不同叶高面。我们知道，在CFD-Post里面有turbo模板，可以生成对应的turboSurface，但是换了后处理软件，就没有办法生成了。
因此，想请教一下，不同叶高曲面表达是什么（图片中蓝色曲面）？应该是在笛卡尔坐标系和柱坐标的转换。查阅了CFD-Post的帮助文档和Fluent的帮助文档里面，都没有这一块内容。所以想发个帖咨询下大家
各位如果有什么线索也可以一起讨论，谢谢！
14da7003-9457-43c8-a206-c5bde261a26f-image.png
8672775e-2671-4021-b7ed-35fe6e11b1fa-image.png

[相关资料]
https://www.youtube.com/watch?v=C6xUFZD6wik - 创建柱坐标截面

实际工程上说的压力是什么压，压力表测到的是什么压，还有就是数值模拟计算得到的压力结果我们应该关注什么压呢？或者说从动压静压总压分布图上我们能看出什么东西。谢谢

业内基本是公认的单通道多通道对CFD影响很大

不过我还没自己实测过，所以不敢说100%，只敢说98%

多谢老师解答，问题已经解决了

@李东岳
感谢李老师的回复，前段时间做了其他的事情，就搁置了这个。哈哈
这个问题解决了，是因为边界条件的设置。
另外请教李老师另外的问题：

OpenFOAM中有求解瞬态不可压多孔介质的求解器吗？ 或者关于二相流中造单纯的流，关于相界面处的边界条件设置问题，李老师有更好的建议吗？
谢谢。:146:

@李东岳

李老师您好。

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

这里Uc尝试过：

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

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

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

可以设计一个量纲的类，里面包含一个整形数组来记录量纲的分量，如[L, M, T]还要定义量纲的运算。然后设计一个标量的类，里面除了浮点数据之外，再加一个量纲，定义标量的动作时也要注意其中的量纲计算。同理还可以定义向量和张量。比如距离的量纲是[1, 0, 0]，速度的量纲是[1, 0, -1]。

参考这个标题Ten years of industrial experience with the SST turbulence model，我想我可以写一个Ten years of academic experience with OpenFOAM CFD code :mihu: :mihu:

重新画了一下，之前那个确实太丑了点。屏幕截图 2024-06-03 221931.png

@桑葚熟了 请问下为什么EDEM软件不行呀？DEM模拟不是会更精确吗？

这个放在他们搞涡那面，估计能讨论一上午 :duang:

我个人也不太熟悉。在我的文章里面，也不太关心涡的情况。

$sdsfa$

$sdsfa$

sdsfa