CFD中文网

非常感谢您的回复与帮助，我尝试一下。

想开个贴讨论下多相流界面追踪类方法目前的前沿水平，例如能算到的最高参数，各类方法都有什么限制等等。目前多相流算法里，计算界面追踪的多相动力学，比如droplet collision, droplet splashing, droplet impact on solid wall，这几个典型案例，目前有用VOF，level-set，phase field，LBM, SPH, 欧拉拉格朗日耦合其它模型的。这几个案例在不同Re和We下会产生许多种演化形态，实验都做出来了，但感觉计算还没跟上。

调研过不下百篇的文献，基本上关键词相关的都翻了一遍，在文献中看到droplet collision大概最好的也就是Re=6500,We=1500，密度比DR=1000左右的，大致是这个级别，像droplet impact有篇08年的用TransAT算到了RE=11000，而最近16年左右的用levelset这些也算到Re=8000，We=1500左右。大多数都是采用了一些改进的数值技巧来稳定模拟，而且网格也都比较密，计算资源消耗也都比较大，估计都得上超算，当然和实验也都对的上，基本都是到了prompt splattering这样的形态。

想讨论下大家都做过或者见过的这类模拟最高参数是多少? 或者往上算各类方法的限制在哪？ 因为多相流在算到高RE和We的时候会数值发散，稳定性应该是主要限制，还有例如VOF和level set会在薄液层处由于数值格式产生人工破碎的非物理情况，需要单独使用人工识别薄液层再稳定的方式，对不规则碰撞并不实用。

测试一下能否随机显示主题

您好，我在使用comsol网格自适应算两相流会出现网格自适应速度跟不上两相运动速度，导致，这方面您是如何解决的呢

本人目前在用rhoPorousSimpleFoam求解器进行多孔介质模型的计算，对angleDuctImplicit算例进行修改。计算条件：固相提供恒定热流密度边界条件。目的：获得固相与液相的温度场。
现在遇到两个问题：
1 无法设置恒定的热流密度边界条件。算例中给的是恒定温度，如果要设置成热流密度需要定义换热区域（固相&液相）。
2 rhoPorousSimpleFoam求解器中的thermophysicalProperties只能设置流体的参数，无法设置固体参数。

请问各位大神这个问题有什么好建议么

请问像这种狭窄处要划分边界层，怎样处理能保证网格高质量呢

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

@ibelief 在 求教fluent后处理-导入曲面（线）-输出曲面（线）参数的实现方法 中说：

或者是分割或者什么，探索一下。 祝

好的 谢谢老师

1.我猜是symmetry这个边界条件，试一下改成wall，然后在0文件夹里自己定义上下的slip边界条件。openfoam会对生成的网格编号进行检查，外部软件的网格节点编号和blockMesh的逻辑不一样。
2.可能是前后的empty边界问题，你可以在starccm里生成3D网格后，重新在openfoam里extrude来生成前后empty边界
类似这样
constructFrom patch;
sourceCase "../w3-d_hc1-3";
sourcePatches (symFront);

// If construct from patch: patch to use for back (can be same as sourcePatch)
exposedPatchName symBack;

// Flip surface normals before usage. Valid only for extrude from surface or
// patch.
flipNormals false;

//- Linear extrusion in point-normal direction
extrudeModel linearNormal;

nLayers 1;

expansionRatio 1.0;

linearNormalCoeffs
{
thickness 1;
}

// Do front and back need to be merged? Usually only makes sense for 360
// degree wedges.
mergeFaces false; //true;

// Merge small edges. Fraction of bounding box.
mergeTol 0;

@yhdthu

感谢，没想到知乎上也有大神研究这个！

你好，请问现在有basilisk的交流群吗，

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

我看了下，确实不是连续的。在ManualInjection.C文件中：

template<class CloudType> Foam::scalar Foam::ManualInjection<CloudType>::timeEnd() const { // Injection is instantaneous - but allow for a finite interval to // avoid numerical issues when interval is zero return ROOTVSMALL; }

我看了一下，可以用cellZoneInject：

Description Injection positions specified by a particle number density within a cell set. User specifies: - Number density of particles in cell set (effective) - Total mass to inject - Initial parcel velocity Properties: - Parcel diameters obtained by PDF model - All parcels introduced at SOI

你试试这个？

PS. 你需要用topoSet做一个cellSet出来toturial里面有

@promisegwj 啊不对 是用拉梅系数···

@流动魅力 我比较好奇你的相间作用力的直径怎么给的

aurora，可以在word里用latex

现在很多AI号，翻译sci基本就是10秒钟一篇，就比如那种jcp昨天见刊的文章，今天国内的马上出来一篇中文爆火全网。这些都是ai写的。感兴趣的可以把类似这些文章的文字复制到这里然后识别一下。

我测试了一半ai写的，一半人工写的，准确率非常高，太厉害了。

https://smodin.io/zh-cn/ai内容检测器

捕1获.JPG

请问有谁对 bodyPoint（refAttachmentPt_）这个函数了解吗？
这是函数的具体实现：

inline Foam::point Foam::RBD::restraint::bodyPoint ( const point& p ) const { return (model_.X0(bodyID_).inv() && spatialVector(Zero, p)).l(); // 双内积 }

@lyc 你好，我很少遇到。即便遇到都是重启就能解决的那种。

导出时候选择了如下这些选项
image.png