CFD中文网

个人建议去看一下北大吴望一老师的流体力学第一章，看了之后推导这类方程简直小菜一碟

@jacob-z

mpirun --oversubscribe -np 20 xxFoam -parallel

应该是进程超限锁死，需要添加申请允许超限，就可以了

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

求问各位，以往的OpenFOAM Workshop 的论文集在哪可以找到？在这个网站里没有呢？https://openfoamworkshop.org/previous-workshops/
我想找一个OpenFOAM Workshop 2016年的一些报告，请问有没有知晓如何获取的？:quwan:

@夜阑烟寒 从昨天开始domecross用不了了，你那是不是也是？我用了得一年半了吧，第一次开始这么长时间用不了

TEMOM只展开$(v+v_1)^{1/2}$，任意阶矩展开任意阶都是泰勒的截断误差，是个原式的高阶无穷小，原式也就是$1/2$次方。
但DEMM把$(v+v_1)^k-v^k-v_1^k$也带上展开了，这个截断误差是变的，k越大，误差越大。越高阶矩，误差的量级越大，虽然是都是无穷小。

所以我觉得越高阶这个方法的优势越不明显。当然都是我猜的:135:
在提出DEMM的那篇文章里的数据也有这样的苗头，低阶矩符合很好，高阶的出点问题又修复了。
我把相关资料发给老师。

@李东岳 并行文件夹里面的2e-08文件夹是有pointProcAddressing文件的，但是后面算出来的2e-02，4e-02就没有这个文件了。所以reconstructPar说找不到这个文件，报错了。

fluent里重力方向应该是可以指定的

@杨英狄 就是说想以这样的方式注入拉格朗日粒子：每个时间步注入流场的粒子的数量、位置都与当前时间步的流场有关，这样的注入方式怎么实现呢？

好像是分内部点（所有面都不是边界）和边界点，内部点根据cell插值（距离倒数分之一），边界点根据face插值。

在某个时期怀疑人生很正常

找工作被拒很正常

找工作跟找对象一样

不只是我，你，也会看到很多人很流弊，工作没怎么找好。或者很多人大学谈恋爱过去了，毕业找了个年薪40万的外企。同样，隔壁老五老实巴交家穷，被富婆看上了，从此人生巅峰。你勤奋读书考研考博，现在还是单身无人问津。

以我自身举例，我总感觉自己很流弊，找工作都是没人要的，这个具体就不细说了。

所以很多时候，努力只是一部分，你需要机遇。机遇只能靠等。

做好定位再打算

如果确定要出国、读博。你英语水平还很好。可以一直找下去。要么是“实验室已满”，要么是“我没有钱如果你自费我们再考虑你”，这种太正常了。我经常遇到这种情况。这说明不了什么。换别的。硕士期间还没有定位研究方向，可以满世界找。

在这里我强调一下，就是不一定要找顶级高校顶级课题组。以我自身举例，2014年出国，我去的是一个世界排名400的学校，导师那个时候还是副教授。但人非常好，我学到了很多东西。很多，包括国内，最强产出者很多来自于讲师与博士后。相反，大教授往往集中精力与拿项目组团队，具体的技术都是博后来做。

所以扩大选择面。

人生三年一变

你回头看你3年前在做什么，6年前在做什么。所以，你才26岁，机会多的很。我2018年1月回国后，还失业了半年，也没影响我做CFD。

给你个2018年我申请北京某高校的拒信截图

捕获.JPG

@bestucan 谢谢老师，我试了下，还是不行~

他是由于none量 他把它删了，看来我的问题就是那个omega，只要用

const volScalarField& omega_ = ( this->mesh_.objectRegistry::template lookupObject<volScalarField>("omega") );

就会找不到，估计是结构问题，难道只能重写了。

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

@yfclark 你好，你这个功能实现了吗，可以请教一下如何实现的吗？

用的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 其他的都没动

@李东岳 哦，明白，谢谢李老师

@李东岳 好滴，感谢李老师，我编一下

对对对，这样是对的，我之前把应变率和湍动能搞混了，谢谢李老师了！

@李东岳 李老师您好，可能是网格质量的问题。我采用简单的单管模型六面体结构网格模拟相同工况（fluent参数设置一致），模拟收敛也没有超温。但是我的实际的结构有点复杂，自己能力有限且时间紧画六面体不太现实。我在ICEM中绘制四面体网格（QUICK(delaunay)）出现一个很奇怪的问题，就是要么体网格就是一个方块（将圆柱计算域包围的方块）且面网格不太好，要么就是有较好的面网格但没有生成体网格，二者的设置区别就是图中的这个选择不一样!另外，采用简单模型则没有出现这种问题

替代文字

reFineMesh可以加密一个方向，这个需要指定区域，你试试topoSet+refineMesh呢》