CFD中文网

在breaking wave领域，这种海浪结构，会将空气卷吸并形成很多小气泡。最近提交一个稿件，审稿人推荐我们讨论一下这方面的研究。我详细看了一下，确实有点意思。

看下图，一些实验研究已经证明，这种海浪拍打过程卷吸进的气泡具有尺寸分布。见过很多文献，都是用的DNS直接模拟研究的，结果也都能证实。

但目前我想的是，能否用PBE模型，或者其他类似模型做相关的研究。

捕获.JPG

气泡直径和分布的关系
捕获.JPG

一些参考文献

Deike, L., Melville, W. K. & Popinet, S. 2016. Air entrainment and bubble statistics in breaking waves. J. Fluid Mech. 801, 91–129.
Wang, Z., Yang, J., Stern, F., 2016. High-fidelity simulations of bubble, droplet and spray formation in breaking waves. J. Fluid Mech. 792, 307–327.
G. Soligo, A. Roccon, and A. Soldati, Breakage, coalescence and size distribution of surfactant-laden droplets in turbulent flow, J. Fluid Mech. 881, 244 (2019).
Ahmed, Z., Izbassarov, D., Costa, P., Muradoglu, M., Tammisola, O. 2020. Turbulent bubbly channel flows: Effects of soluble surfactant and viscoelasticity, Computers & Fluids 212, 104717.

个人建议，看这个之前，看看icoFoam，那个是绝对的经典~:146:

@chengan-wang 或许你可以搞一个局部细加密。

@qhwj1990 这几个数值就是单位的差别把，kj=1000j, kgmol=1000mol

@东岳 光看流场的矢量图可能确实看着还凑合对称，但是假如我要在出口得到RTD曲线（停留时间分布曲线），就会得到对称出口明显不同的RTD曲线？

东岳大神、各路大神好：
最近在模拟一个上下振动混合器中的单相三维流场，如图。 考虑到计算时间和成本问题，想采用moving frame of reference的方法来代替动网格进行振动的模拟。目前设置inlet速度为振动速度，壁面和出口都设置的是wall，速度也是振动速度，而板和轴的速度为0.但得到的流场似乎不太对，想请教下边界条件的设置存在哪些问题？WeChat Image_20190929112939.jpg

是的，理想气体，明天我更新下。

另外，温度本身不是守恒变量，因此，CFD一般求解守恒变量，但是温度方程通常比较简单，一些不太成熟的代码会使用温度防尘

谢谢李老师~

你需要系统看一下流体力学的相关知识，可能你没有这方面的背景所以不太理解。对物理空间中的一块区域而言，其所包含的物理量有几种变化方式，要么它自己随时间在变化（变化率，比如密度因为温度升高降低了），要么因为有流动带它进来或者出去（对流，比如能量），要么它自己在向外扩散或者外部在向这块区域扩散（扩散，比如这一块是高浓度区域而周围是低浓度），要么这块区域中这个物理量自己在生成或者消灭（源项，比如有化学反应）。所有的这些变化，都需要满足物理学的基础即守恒定律，其数学表达就是所谓的“广义”NS方程。这里“广义”是指这个方程描述了物理量的一般过程。

phi等于1，这个方程代表质量守恒；phi等于速度，这个方程代表动量守恒；phi等于内能，这个方程代表能量守恒；phi等于某个物质的浓度，那么这个方程就代表该物质的质量守恒，等等。

@东岳 谢谢东岳老师哈

各位前辈好！
我对fluent的export ASCII file 文件有些疑惑如下：
该输出数据可以选择节点或者是网格，我输出两个点在同一个网格内，按网格输出是两点的数据是一样的，但是按节点输出是不一样的。（我现在对节点的认识将其认为是网格的端点，不知这种认识是否正确。）按我现在的认识的话，两个点的位置是不在任一节点上，但又因为两个点的输出数据不一样，所以应该是有数值差分做了处理，那么fluent是如何进行差分？fluent的理论指南和用户指南，只粗略介绍是按节点输出，不是很能理解。还望大佬们指点一二，感谢不尽！！

我用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 }

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

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

填坑，参考官方t13教程（没记错的话，反正题目有个stl，似乎只有这种方法）

@bart 慢慢增大进出口压差计算，在看看结果

运动物体计算，浸入边界。

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

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

各位大佬，请问CFD-DEM coupling时，流体的控制方程组求解，用耦合式解法还是分离式解法？
非常感谢！