CFD中文网

@低碳生活 请问帮助文档中有关随机数的教程在哪里？我也没有找到/:xinxin2:

这个问题也可以参考我之前的回答：

http://120.25.245.130:4567/topic/94/les和rans中-离散的对象是相同的吗/3

@李东岳 在 请问UprimeMean是如何计算的 中说：

所以你要计算得是$\overline{{\bfU '}^2}=\overline{\bfU^2}-{\overline{\bfU}}^2$？$\bfU '$是脉动速度，$\bfU$是速度，$\overline{\bfU}$是平均速度

老师，请问openFoam中可压缩求解器（rhoreactingFoam）使用Les模型，使用uprimemean计算的是favre平均的值吗还是雷诺平均的值，如果想要分别计算favre平均和雷诺平均可以做到吗

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

这个看起来是边界条件设置的问题，正常用的话不会这样，网格少的话发在论坛我给你看看

https://www.cfd-online.com/Forums/fluent/40645-non-newtonian-fluid-fluent.html

@李东岳 在 OpenFOAM里面的dnsFoam为什么这么菜？ 中说：

@ZifeiYin

Hi,

请参考下述内容，把第一行第二句话自行删去：

鉴于CFD圈子的小众性，“CFD中国”执行最严厉的友善管理制度：所有用户请以专业的态度对待所有问题。不得在所有版块发表轻蔑性言论；
例如，下面的语句被认为是具有轻蔑性的：“楼主回去看看书吧”，“这种问题你也问？”，“不会用百度？”
“CFD中国”建议用户对每个问题严谨友善处理。例如上述内容可以这样表述：“建议楼主参考《数值传热学》第48页，有对xxx问题的详细表述”，“个人认为这个问题是一个非常普适性的问题，在大部分的教材中都有涉及，楼主可以参考xxx”，“这个问题和CFD的关系不大，建议楼主去百度。google去搜索相关内容”

哈哈哈哈，东岳老师干得漂亮

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

libevent_pthreads-2.1.so.7

libevent-2.1.so.7

libhwloc.so.15

libmpi.so.40

libopen-pal.so.40

libopen-rte.so.40

这几个包是在一个刚装完ubuntu系统的机器上 1）如果能从别的地方拷贝一个装好的openfoam 2）把这些包复制到lib下，就可以跑openfoam

求助模拟大神给看下，这个case中结晶设置是否有问题？想模拟液固流化床氯化铵冷却结晶过程，不知道传质那部分怎么设置？没找到帮助文件中mass mechanism 下拉列表中的population balance机制

2016年autodesk公司把CNN用于不同外形的流场预测。但是CNN这面最爆炸性的成果当属2021年google那面发表在PNAS的工作。目前被引700多次。

在Kochkov等的工作中，作者们在4倍糙化的网格上，实现了传统CFD做直接模拟的精度。这一篇工作并不是首创。google在2021年的另外一篇文章中，将他们的思想最开始应用于一个简单的对流传输方程中。

2021年稍早，哈佛大学与google公司合作，在Zhuang等的工作中提出了Learned Intepolation。即通过大量的样本来训练出来的CFD格式。

这对于咱们做CFD的很好理解，传统的面插值格式，可以直接从网格体心插值过来，也可以从2个网格体心或者多个网格体心取平均过来。对于CNN的插值格式，是通过针对样本进行多次采样，实现一个类似可以理解的在糙网格下，实现细网格上同样结果的插值格式。更好的理解可以这样：

在细网格下，传统CFD预测的面插值是1，糙网格下，CFD预测的是0.5。CNN经过多次特征与标签的采样之后，在糙网格下预测的插值可能是0.9。通过这样来实现糙网格的高精度模拟。

在google公司的工作中，除此之外，全部用的是CFD传统的解法。这对于使用深度学习从头到尾的计算，有很多优势。最大的优势就是，可以保证传统CFD的守恒特性。

在上图中，第1行是CFD计算的细网格的结果。第4行是CFD糙网格计算的结果。第2行，是CFD细网格的结果，通过平均技术投影到糙网格的结果。第3行，是CNN计算的结果。

可以看出，第3行与第1行结果相当！也就是CNN的“学习来的数值格式”可以在糙网格上进行直接模拟。

问题还是有的，这部分数据样本还是来自于CFD。如果CFD本身样本不对。那学习的样本本身就是错误的。

需要保持警惕的是：从错误中来学习，必然是大错特错。

google的工作是通过他们自家的机器学习代码jax来实现的。google的工程师为了弄这个，在jax上面写了个错位网格的CFD代码。。

CFD自家亲戚都懂的，现在很少有错位网格的CFD代码了。

因为google那面CFD层面比较初级，所以需要处理成同位网格，时间步进还没有植入高阶格式，多重网格求解器也还没有，边界条件只有周期边界条件。。。

目前，google团队还在公开欢迎各类CFD大佬去参与合作。岳子就是没有闲工夫，要不然岳子也要躺一躺这个浑水。

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

🏷️【职位信息】
研究方向为计算流体力学、多相流、界面追踪、空天发动机喷雾燃烧、壁面冷却、燃烧室优化等。

博士生和博士后要求：

中华人民共和国公民，拥护中国共产党的领导，品德良好，遵纪守法； 取得工程热物理、能源工程、航空航天、或流体力学等学科的本/硕/博学位； 具有良好的团队合作和科研热情，严谨的科学态度，吃苦耐劳，积极向上，具有良好的英语写作能力； 毕业于双一流、985院校优先。

博士生福利待遇：

将为博士生提供有竞争力的补贴和待遇，月薪不低于5000元，科研成绩突出者可申请诸多奖学金，并获得课题组额外奖励； 博士期间有多次参加国际、国内会议和交流的机会；与众多知名大学（如斯坦福大学、北京大学、浙江大学等）建立了良好的合作关系，可为科研能力优秀的同学提供国际合作、访学以及海外高校深造机会等； 优秀的学生博士毕业后可留在课题组继续工作或推荐到海外著名高校做博士后。

博士后福利待遇：

年薪不低于30万（含18万元深圳市政府补贴），优秀者可以面议； 鼓励申请“博士后创新人才支持计划”，国家可给予每人两年60万元的资助，其中40万元为博士后日常经费，20万元为博士后科学基金； 优秀博士后出站前可申请学校优秀博士后奖励（比例约30%），一次性奖励12万； 福利：享受五险一金； 住房：可租住深圳市人才公寓，或享受深圳市政府租房补贴； 海外全球前200高校毕业生待遇更优，详见深圳市政策。

申请流程：

申请博士和博士后的学生，请以“【应聘博士/博士后+姓名】”为主题发送个人简历至邮箱：shaochangxiao@hit.edu.cn。申请截止日期为2024年9月22日。 博士研究生采用申请考核制，不需要进行博士入学考试。对于优秀学子，推荐通过《哈尔滨工业大学博士研究生招生优秀生源快响行动》来申请。

🏷️【课题组简介】
哈尔滨工业大学（深圳）流体和燃烧模拟课题组的负责人是邵长孝，目前担任教授、博士生导师。加入哈工大深圳之前，曾在浙江大学能源工程学院获得博士学位，导师为罗坤教授、樊建人教授和岑可法院士，之后在美国斯坦福大学湍流研究中心（Center for Turbulence Research）从事两年博士后研究。
课题组主要从事空天发动机复杂多相流动和燃烧的模拟和工业软件开发等研究，近年来在算法、机理、模型方面取得了多项重要成果，并发表在Prog. Energ. Combust. Sci., J. Comput. Phys., J. Fluid Mech., P. Combust. Inst.等期刊。主持国家自然科学基金面上项目、青年项目，以及多项省市项目，荣获广东省珠江青年人才、深圳市海外高层次人才、深圳市南山区领航人才等荣誉。详细信息见个人主页：https://faculty.hitsz.edu.cn/shaochangxiao。

🏷【学校简介】
哈尔滨工业大学深圳校区由哈工大与深圳市政府合作共建，以全日制本科生与研究生教育为主、非全日制教育为辅，是国家“985工程”建设高校和“双一流”建设A类高校首个在深圳市举办的本硕博培养体系齐全的高水平校区。
校区在办学过程中始终坚持扎根深圳，面向国家重大需求，立足国际学术前沿，紧密结合广东省、深圳市经济社会发展战略目标，坚持“高等教育改革的试验田、汇聚高端人才的桥头堡、培育创新成果的策源地、国际合作办学的示范区”的办学定位，贯彻落实哈工大“一校三区”融合发展要求，开展党建思政、人才培养、师资队伍建设、学术研究、国际合作交流等各项工作。

天天扯淡第一名

Jacobi是并行性最好的算法，LU-SGS并行性较差。

管内气体流动，管壁定温加热，我在计算的时候设置了report：热流密度，对流换热系数，Nu等参数。发现用热流密度反推的对流换热系数和Nu和报告的对不上，相差很大，请教各位老师这是什么原因？谢谢

工科末流985硕士毕业，想在武汉找CFD相关的工作，家人们有没有推荐的公司。

@搬运工

高斯迭代法处理20阶矩阵就已经非常慢了。通常和多重网格+高斯赛德尔求解器一起处理。并且OpenFOAM这种成熟的商软里面矩阵求解器算法非常精炼。通常由主求解器+预处理+光顺+重组构成。多重网格更是复杂。另外，在编程的时候使用动态、寻址的方法寄存系数。

所以如果你要写一个快速的求解器：

使用外挂包； 潜心研究稀疏线性系统求解理论；我建议花大量时间研读：YousefSaad, 2003。看看大型稀疏矩阵怎么求解的。而不是3/4阶矩阵的高斯塞德尔求解器。:cheeky:

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