CFD中文网

得看人是不是明白人。明白一口口饭都是从哪来的，明白jio底下站的地为什么这么结实，明白前方的路为什么这么清晰。

碰上明白人，给高了会用富裕的资源提升自己（给多少钱干多少活，否则钱拿不稳）；给低了也不担心，毕竟瞎还能活着的概率不是很大，找个不瞎没有特别的困难。

依靠某个人、某个环境下对自己的态度、待遇来判断自己的价值。这显然不是个明白人……给高了骄傲，给低了还怨恨呢。不能自知，怎么都是一件很危险的事情。有核心能力的也有不明白事的：仰天大笑出门去，归来依杖自叹息:136:

而能碰上多少明白人又取决于所在行业对社会精英的吸引能力和自己在行业的层次。总之不努力就整天和糊涂蛋扯吧……:135:

还有这种组织心态，把员工当小孩、追胡萝卜的驴子（有些事情你不需要知道，我在掌着。我要塑造你的认知、影响你的判断）。当然在这个繁复多样的世界，这种方式无法完全避免。但是显然有更好的方式，只是对组织中的每个人要求比较高。

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

看起来挺有意思

微信图片_20200619134942.jpg

我用OpenFOAM计算外流问题，当来流速度设置为0时可以正常计算，但当来流速度大于零时最多计算一步流道就堵塞了，我怀疑是边界条件有问题。我的边界条件设置如下，请教大家这样搭配可以吗？如果不行应该是怎样呢？
U

FoamFile { version 2.0; format ascii; class volVectorField; location "0"; object U; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 1 -1 0 0 0 0]; internalField uniform (0 0 0); boundaryField { //- Set patchGroups for constraint patches #includeEtc "caseDicts/setConstraintTypes" inlet { type fixedValue; value uniform (0 0 50); } outlet { type zeroGradient; } symm { type slip; } propellerW-f { type movingWallVelocity; value uniform (0 0 0); } propellerW-s { type movingWallVelocity; value uniform (0 0 0); } propellerW-h { type fixedValue; value uniform (0 0 0); } propellerW-t { type fixedValue; value uniform (0 0 0); } }

p

dimensions [1 -1 -2 0 0 0 0]; internalField uniform 101325; boundaryField { //- Set patchGroups for constraint patches #includeEtc "caseDicts/setConstraintTypes" inlet { type zeroGradient; } outlet { type fixedValue; value uniform 101325; } symm { type slip; } wall { type zeroGradient; } }

T

dimensions [0 0 0 1 0 0 0]; internalField uniform 293; boundaryField { inlet { type fixedValue; value uniform 293; } outlet { type inletOutlet; inletValue uniform 293; value uniform 293; } symm { type slip; } "inter.*" { type fixedValue; value $internalField; } wall { type zeroGradient; } #includeEtc "caseDicts/setConstraintTypes" }

nuTilda

dimensions [0 2 -1 0 0 0 0]; internalField uniform 1e-5; boundaryField { //- Set patchGroups for constraint patches #includeEtc "caseDicts/setConstraintTypes" inlet { type fixedValue; value uniform 1e-5; } outlet { type inletOutlet; inletValue uniform 1e-5; value uniform 1e-5; } symm { type slip; } wall { type fixedValue; value uniform 0; } }

nut

dimensions [0 2 -1 0 0 0 0]; internalField uniform 8.58e-6; boundaryField { //- Set patchGroups for constraint patches #includeEtc "caseDicts/setConstraintTypes" inlet { type fixedValue; value uniform 8.58e-6; } outlet { type inletOutlet; inletValue uniform 8.58e-06; value uniform 8.58e-06; } symm { type slip; } wall { type nutUSpaldingWallFunction; value uniform 0; } }

alphat

dimensions [1 -1 -1 0 0 0 0]; internalField uniform 0; boundaryField { inlet { type calculated; value uniform 0; } outlet { type calculated; value uniform 0; } symm { type slip; } "inter.*" { type fixedValue; value $internalField; } wall { type compressible::alphatWallFunction; Prt 0.85; value uniform 0; } #includeEtc "caseDicts/setConstraintTypes" }

其中inlet，outlet，分别是进口和出口，几何边界条件为patch；symm为四周的自由边界，为柱状，几何边界条件为patch，propeller是计算的螺旋桨，几何边界条件为Wall；计算所用湍流模型为SA模型，所用求解器为rhoSimpleFoam,设置了两块MRF区域，MRF设置如下：

FoamFile { version 2.0; format ascii; class dictionary; location "constant"; object MRFProperties; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // MRF1 { cellZone sector_part1; active yes; Fixed patches (propellerW-f); nonRotatingPatches (); origin (0 0 0); axis (0 0 1); omega 112.5737; // rad/s } MRF2 { cellZone sector_part2; active yes; Fixed patches (propellerW-s); nonRotatingPatches (); origin (0 0 0); axis (0 0 1); omega -112.5737; // rad/s }

可能北京今天发誓的人有点多，老天看不过去了

嗯 AP 很多AP都很强力 强力AP！

结构化网格非结构化网格按照网格是六面体还是体面体区分好像是不完全正确的。我看一个老师讲解说，结构化和非结构化区分是文件保存格式的区分，结构化是使用索引，非结构化是把每个网格都区分出来了。我看各位大佬上面回答可能是直接讨论的四面体网格，因为我最近也在学openfoam和icem，icem可以通过自动化分成四面体非结构网格，也可以分块化成6面体网格，但是我通常画完6面体网格后，icem有可以转化成非结构网格按钮，将画好的网格转化成非结构体网格（转化后依然是6面体）。
所以如果你要话非结构化网格也可以分块（画了6面体转化一下），如果你说的是画四面体非结构化网格分块的话，这方面知识我就不知道了

毕超的《计算流体力学有限元方法及其编程详解》这本书里面详细介绍了方程各项的离散，不过使用MATLAB编写的，我现在主要就是参考这本书的

你是不是想把曲面展开为一个平面？https://blog.csdn.net/mrbaolong/article/details/106949340

谢谢李老师~

PLIC=piecewise linear interface calculation;
The geometric reconstruction scheme represents the interface between fluids using a piecewise-linear
approach. (Ansys Fluent Theory Guide 2022R1, p579)

@李东岳 好的，谢谢东岳老师

大家好，我最近想做一个这样一个多孔介质：在多孔介质的不同部位，接触角也是不同的。即通过指定两点坐标，来设置两点所包含几何体的接触角，类似于setFields中设置水相分布的功能。这能否通过在openfoam中编译程序来实现？谢谢大家了。

算例为土柱冻结算例，起始土柱内部温度为6.7℃，土柱顶部施以温度为-6℃的循环液流(the circulating fluid)，使土柱从上到下逐步冻结，下图为示意图，
捕获.JPG
温度边界条件：
变热通量上边界条件/a variable heat flux upper boundary condition：
qh=hc(TCoolant-TTop), hc=28 W m-2 K-1
其中：TCoolant是循环液流的温度值(-6 ℃)，而TTop是土柱顶部温度，我理解的是每一次时间步长计算后都会更新，所以顶部通量是一个动态更新的过程。
1、之前有看过关于fixed heat flux boundary的帖子，采用fixedGradient边界条件代替，但没找到关于 variably heat flux boundary condition
2、 采用swak4Foam？
希望能给我一些建议，共同探讨下。非常感谢！

应该还可以用这个办法来判断
if(THREAD_TYPE(t) == THRAD_F_WALL)

实验室有个试验台，里面充满水，水流速度最大60m/s，无叶轮，最近要做单相流实验，想确定增压泵打多少压力水才不会空化。按照我的想法，100℃水温时，水的汽化压力为0.1MPa，要保证不空化，需要线速度最大的地方局部压力达到0.1MPa，按照伯努利方程，需要加上动压，这个时候发现系统压力居然要达到2.8MPa。仔细想想应该是不适用伯努利方程的，所以想问一下，有没有什么好的方法可以估算系统压力呢？

@李东岳 谢谢老师解答，之前一直出现这个问题导致每次画完网格都得要重启，我重新装下系统在试试

$sdsfa$

$sdsfa$

sdsfa

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

不需要