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.

你是不是算的倒虹吸啊

@Zhy2022 好的谢谢！之后有时间看看论文再来填坑:high:

我搜了一下，基本无解，就是用 tecplot 里的命令操作（一般功能比图形界面点点更丰富全面），也只有两个和有效数字相关的命令，还是用于显示的。

它毕竟是个可视化软件，虽然笼统地说做后处理的。可视化是后处理的一部分，但对数据操作的能力毕竟有些欠缺。

@Shane 非常感谢

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

最近在编写程序的时候遇到并行边界上网格面心处数值通量(φ_f)的同步问题，单核运行不存在并行边界，所以没有这个问题。以下图为例：

Domain.png

并行计算时，求解区域分为了4个区域，两个区域之间产生了并行边界面，对于内部网格面元，网格面心处的数值通量(φ_f)的计算跟单核运行时是一样的，但是对于并行边界处的网格面元，数值通量(φ_f)在两个相邻并行边界上应该是一样的，因为并行边界的存在不应该改变原来的解。因此，我用了如下的syncFaceList()函数来进行同步，alpha_f是定义在网格面心上的surfaceScalarField，这个函数在头文件"syncTools.H"中有定义，我用了以下三种同步方式，发现流体穿过并行边界时还是产生了误差，单核运行时误差(数值解与解析解之差)可以达到10E-12次方量级，但是并行运行时流体一穿过并行边界，误差就上升到了10E-4量级，我觉得这是因为并行边界上数据同步没有做好的原因，请问有小伙伴遇到过这个问题嘛？使用syncFaceList()函数的正确方式是什么？如果使用syncFaceList()函数进行同步网格面心处的数值通量(φ_f)还是有问题的话，打算自己写一下同步过程。

// syncTools::syncFaceList(mesh(), alpha_f, maxEqOp<scalar>(), false); syncTools::syncFaceList(mesh(), alpha_f, plusEqOp<scalar>(), false); // syncTools::syncFaceList(mesh(), alpha_f, eqOp<scalar>(), false);

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

advection:可以理解为平流，主要由物质浓度存在梯度引起的变化。convection：可以理解为对流，主要是由热量不均匀引起的垂直方向的对流。好像在流体力学中对这两者区别不是很大:xiabanle:

@李东岳 安装ansys 时就可以同时安装cfx么？按照我经验，目前的ansys学术版里面可以直接安装cfx，不需要任何额外的操作。

https://www.bilibili.com/video/BV1or421x7Ut/?spm_id_from=333.999.0.0&vd_source=5dfca494f85875fa7f031a9f632b1a89
第9分钟介绍这种High Aspect Ratio对于瞬态、稳态的具体影响。

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

烦请各位大佬指点一下，workbench里面mesh画的带角度的翼型的网格，为什么边界层厚度不一样啊？怎么处理一下呢？上面的感觉太窄了

e9d97488-fa0f-453d-9f2b-1884d5136695-image.png

感谢版主哦 :146:

回复: 动网格的旋转运动可以修改旋转轴吗？

你好 我最近也在研究dynamicMeshDict (https://openfoamwiki.net/index.php/DynamicMeshDict). CofG指的是centre of gravity. 在dynamicMeshDict你可以修改 centre of mass, inertia and etc. CofR是centre of rotation,这个我的算例是在controlDict里修改，我用的是interFoam solver，openfoam version7.

@ibelief 在 室内室外空气耦合 中说：

好的，多谢

@李东岳 是的，我手动更改之后可以计算了，谢谢李老师:huahua:

现在想让液滴撞击壁面后完全蒸发或者蒸发一部分，有没有研究过壁面模型的大佬。感觉要将相变模型耦合到壁面模型中，难度有点大啊:135:

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