CFD中文网

虚拟机上openfoam10和11没有问题

各位老师w71h11131532_1601213641_696.png ，我想模拟一下油滴在超临界水里下落过程中的对流和扩散问题：
初步打算用VOF模型来捕捉相界面，但有一个问题VOF并不能解决扩散的问题。请问有没有什么比较好的解决方案？
如果用UDF/UDS来解决扩散项的问题，请问能不能分享下类似的代码？

欧拉欧拉那面由于已经在统计学方面已经处理过了，没有噪音。欧拉拉格朗日模拟这面目前有什么文章在处理么？找了找，还没看到。

As an alternative, Lagrangian methods ‘‘discretize” the density function into ‘‘parcels” that are simulated using Monte-Carlo methods. While quite accurate, as in any statistical approach, Lagrangian methods require a relatively large number of parcels to control statistical noise, and thus are computationally expensive.

最近看到了一个新的并行分解算法，KaHIP，依然存在疑问。不清楚这种不同的挂载方式的数学本质是什么。

http://algo2.iti.kit.edu/kahip/

替代文字

@youhaoyu 我是使用compressibleVoF中的throttle算例更改了模型和边界条件

楼上说的对。版本差的太大，移植相对麻烦。

0_1516841984913_捕获.PNG
动图1：河床沙石冲刷运动

河流的流动会改变河床的形状。底部的沙子以及大石头通常会分为两个区：

顶层的bed-load zone（有研究河海的么，这个怎么翻译？），其中的颗粒会随之移动； 底层：沙子以及石头会被困住，但是会缓慢的爬行；

在碎石较多的河床中，较大的沙石倾向于堆积在顶层，这种现象被称之为armoring。个人感觉之所以称之为armoring是因为较大的沙石趋向于在表面，类似一种武装。

实验表明，在这个区域，大沙石有一个净得向上的运动速度，小颗粒具有一个小的向下的移动速度。 目前对于这种现象的原因尚不清楚。不过现存了一些潜在的理论。例如支持程度比较高的巴西果效应（Brazil Nut Effect）。

0_1516841984913_捕获.PNG
动图2：有人晃荡大米里面的密度差不多的小石头

替代文字
动图3：进行的计算模拟

巴西果效应是指如果把两种颗粒的混合物置于容器中，然后施加外加的振荡，体积比较大的颗粒会上升到表层，而较小的颗粒会沉降到底部。关于这种古老的效应的动力学机制至今仍众说纷纭。

DEM模拟：文中作者采用LIGGGHTS对这个现象进行了模拟，下图可以看出，随着时间的推移，大颗粒漂浮出来。

0_1516841984913_捕获.PNG
图4：Nature Communication中的LIGGGHTS模拟

看起来很有意思 :cheeky:

本文来自Nature Communication

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

博士生和博士后要求：

中华人民共和国公民，拥护中国共产党的领导，品德良好，遵纪守法； 取得工程热物理、能源工程、航空航天、或流体力学等学科的本/硕/博学位； 具有良好的团队合作和科研热情，严谨的科学态度，吃苦耐劳，积极向上，具有良好的英语写作能力； 毕业于双一流、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类高校首个在深圳市举办的本硕博培养体系齐全的高水平校区。
校区在办学过程中始终坚持扎根深圳，面向国家重大需求，立足国际学术前沿，紧密结合广东省、深圳市经济社会发展战略目标，坚持“高等教育改革的试验田、汇聚高端人才的桥头堡、培育创新成果的策源地、国际合作办学的示范区”的办学定位，贯彻落实哈工大“一校三区”融合发展要求，开展党建思政、人才培养、师资队伍建设、学术研究、国际合作交流等各项工作。

@东岳 谢谢教授

@肖艺 STAR CCM+
例子： https://www.cfd-schuck.de/en/cfd/anwendungsbeispiele/simulation-of-gear.html

@冰冰玩 好的，谢谢！我之后尝试一下。

@blgy 非常的感谢，你这个是在哪看到的啊？你这个图片也看不出来用什么软件，用什么方法划分的啊，仅仅是欣赏一下:ok:

@东岳 是我搞错了，向量混合积跟微分算子不一样

@李东岳 http://dx.doi.org/10.1016/j.compchemeng.2015.09.008
这篇文章对三种重构方法做了对比。

之前有点问题，最终如下
\par $C_p$采用4次多项式分段拟合
\begin{equation}\label{equ:NS_Cp}
C_p(T)=a_1+a_2T+a_3T^2+a_4T^3+a_5T^4
\end{equation}
\par静焓
\begin{equation}\label{equ:NS_H}
H(T)=\int_{T_{0}}^{T_{x}}C_p(T)dT=a_1T+\frac{a_2}{2}T^2+\frac{a_3}{3}T^3+\frac{a_4}{4}T^4+\frac{a_5}{5}T^5+a_6
\end{equation}
\par熵
\begin{equation}\label{equ:NS_S}
S(T)=\int_{T_{0}}^{T_{x}}C_p(T)\frac{dT}{T}=a_1\ln{T}+a_2T+\frac{a_3}{2}T^2+\frac{a_4}{3}T^3+\frac{a_5}{4}T^4+a_7
\end{equation}
\par$\bullet$求解静温(已知总温和马赫数)
\begin{equation}\label{equNSUs}
U_s^2=2\left(H(T_{tot})-H(T_{sta})\right)
\end{equation}
\begin{equation}\label{equNSmach}
Mach^2=\frac{U^2}{\gamma(T)R_gT}=\frac{U^2}{\frac{C_p(T)}{C_p(T)-R_g}R_gT}
\end{equation}
\par由(\ref{equ:NSmach})和(\ref{equ:NSUs})得
\begin{equation}\label{equNSTsta}
T_{sta}=\frac{2\left(H(T_{tot})-H(T_{sta})\right)}{Mach^2\frac{C_p(T_{sta})R_g}{C_p(T_{sta})-R_g}}
\end{equation}
$\bullet$求解静压(已知总温、总压和静温)
\par由p等熵过程
\begin{equation}
ds = C_p(T)\frac{dT}{T} -R_g\frac{d p}{p}=0
\end{equation}
\par两边同时积分有
\begin{equation}
\int_{T_{tot}}^{T_{sta}}C_p(T)\frac{dT}{T} =\int_{p_{tot}}^{p_{sta}} R_g\frac{d p}{p}
\end{equation}
\par记
\begin{equation}
S(T_{x})=\int_{T_{0}}^{T_{x}}C_p(T)\frac{dT}{T}
\end{equation}
\par则
\begin{equation}
S(T_{sta}) - S(T_{tot}) = R_g\ln\frac{p_{sta}}{p_{tot}}
\end{equation}
\par那么
\begin{equation}\label{equNS_psta}
p_{sta}=p_{tot}e^{\left(\frac{S(T_{sta})-S(T_{tot})}{R_g}\right)}
\end{equation}

@两月三年 是的，这样显示的时候是表示选上了。注意选的位置就好。

空气龄这个概念在通风领域应用比较广泛，我对此不大了解，但觉得这个物理概念上和停留时间相似。
因此扩散系数本质上于空气或者示踪气体扩散系数相似。
相关理论可以参考一下 https://age-of-air.forumotion.com/

openfoam 里面的空气龄有个类似的solver， 可以参考一下 （19楼）
https://www.cfd-online.com/Forums/openfoam-post-processing/68300-mean-age-air.html

只有一个图。那就修网格吧。

在方柱计算中，对升力系数进行FFT分析，取样数据点在10000左右，得到如下结果，主频与实验较为吻合，在高频处出现谐振，是由于方柱前缘分离剪切层中K-H不稳定性导致的升力系数高频振荡吗？谐振又是怎么产生的呢？
诚心向各位前辈请教，求指导！
FFT Power.jpg