CFD中文网

测试编辑

SV_locate_point (real [], CX_Cell_Id *);
你可以用这个函数来获取这个点的所在网格的信息 要自己去看下头文件
如果这是你的要求的话
而如果不是的话 那么只是一个纯解析几何的问题

@李东岳 感谢李老师的回复，一下子就通透了

@东岳 请问东岳老师有关于双欧拉解决相变问题的例子吗，附带算例内没有找到有关相变的设置

@bestucan 谢谢大神的解答，受教了

@李东岳 多谢李老师鼓励，还得感谢李老师的讨论答疑和提供的交流平台:ok: 。

@allvic 是不是可以在边界层内设监测点，得到频率分布，与T-S波的频率进行相关性验证呢？

@King_RVM 你好，使用starCCM绘制网格，然后导入到openFoam里面跑算例，计算结果会出现质量不好的情况么？请教一下~

在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.

@李东岳 感谢李老师，我试试看能不能画出来

见过一些，目前也不好解释。:141:

如题。

说一点我的看法，不一定对，仅供参考：

non-orthogonality

dca809a8-9cd5-4f74-9329-6d8459966264-image.png

non-orthogonality由snappyHexMeshDict.meshQualityControls.maxNonOrtho直接控制，你的snappyHexMeshDict里这个值是45，这个值一般不需要控制的这么小，65是更常见的配置值。

另外，你的snappyHexMeshDict里，relaxed.maxNonOrtho给的是75，这个值通常是和maxNonOrtho = 65相配合的。如果你真的要把maxNonOrtho设为45，按逻辑来说relaxed.maxNonOrtho也应该相应的减少一点。

网格纵横比

这里截取你的背景网格配置：

vertices ( (-340 -302 -3) // Slightly smaller than the STL model bounds (332 -302 -3) (332 341 -3) (-340 341 -3) (-340 -302 95) // Slightly larger than the STL model bounds (332 -302 95) (332 341 95) (-340 341 95) ); blocks ( hex (0 1 2 3 4 5 6 7) (100 100 100) simpleGrading (1 1 1) // Adjust mesh density as needed );

如上面B老师所说，你网格的大纵横比就是这里导致的。如果不是有意要在Z方向上做加密的话没必要这么分块，我口算一下，大概(70 65 10)就可以。

castellate

从上面可以看到，你背景网格的量级是100 * 100 * 100 = 1,000,000；而你SHM的maxGlobalCells给了2,000,000。这看起来不太对，这几乎没给SHM进一步细分的空间，从结果来看也是这样，网格划分的不够细。

另外，截取一段level的配置：

features ( { file "building.eMesh"; level 3; } ...... ); refinementSurfaces { building { level (3 4); } ...... }

我理解features的level应该不低于Surface的。

胡乱写了一些个人看法，希望能有所帮助。

可能北京今天发誓的人有点多，老天看不过去了

Really? 那得谢谢王博士！

@Mina-Lee 在 precice tourorials中的受力问题 中说：

@RolandLent precice中没有计算OF力的能力，力都是在OF里计算的。OF中的力也是通过force函数计算的。
pecice只干三件事:
1、引用of中的力变量，提取流固耦合交界面的力变量；引用位移变量。2、precice内部运算。3、覆盖引用的位移变量。

我阅读了第一位回答者的源代码，力的计算是precice的of适配器实现的，我在适配器的forcebase类里边找到了计算力的方法，这里附上计算压力和黏性力的部分

if (solverType_.compare("incompressible") == 0) { forceField.boundaryFieldRef()[patchID] = surface * pb[patchID] * rhob[patchID]; } else if (solverType_.compare("compressible") == 0) { forceField.boundaryFieldRef()[patchID] = surface * pb[patchID]; } forceField.boundaryFieldRef()[patchID] += surface & devRhoReffb[patchID];

如果是用of的solid求解器就是向您这样的执行方式了

// Check if a force field with the requested name exists. // If yes (e.g., solids4Foam), bind Force_ to that field. // If not (e.g., pimpleFoam without the Forces function object), create it. if (mesh_.foundObject<volVectorField>(nameForce)) { Force_ = &const_cast<volVectorField&>( mesh_.lookupObject<volVectorField>(nameForce)); } else { ForceOwning_.reset(new volVectorField( IOobject( nameForce, mesh_.time().timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), mesh, dimensionedVector( "fdim", dimensionSet(1, 1, -2, 0, 0, 0, 0), Foam::vector::zero))); Force_ = ForceOwning_.get(); }

能用钱解决的问题都不是问题 :haqi: :haqi: :haqi:

@walten 你的回答很有帮助

你算的是什么工业过程有没有参考的SCI可以看看？

计算结果的不合理之处在于，壁面处沉降的固相体积分数仅为5%。

跟实验对不上？这不好说，影响因素太多了。你用的什么模型？