CFD中文网

psb.jpg

春风催碧叶
夜灯盼文章
横批：诸事如意

改了一下
春风催碧叶
夜灯盼收敛
横批：哈哈哈哈

碧叶双关
收敛双关

@东岳 嗯嗯，确实是尺度问题，谢谢东岳老师！

确实和概率论有关，谢谢，记录一下

\begin{equation}
\left( \begin{matrix}
k\\
x
\end{matrix}\right )=\frac{k!}{(k-x)!x!}
\end{equation}

@东岳 李老师，我没推出来

@shuchang769 一样的思路啊

@李东岳 李老师您好，感谢您提出的意见，现在只是算的有些慢但是似乎问题已经解决了，谢谢。33.png

生诞快乐 :chitang:

@李东岳 李老师您好，D()是一个volVectorField。主要是发现建立两个fvMatrix时，前面一个fvMatrix会对后面的fvMatrix产生影响，主要体现在源项上。按理说这是两个独立的fvMatrix对象，不应该互相影响啊

@bestucan 谢谢老师，我试了下，还是不行~

他是由于none量 他把它删了，看来我的问题就是那个omega，只要用

const volScalarField& omega_ = ( this->mesh_.objectRegistry::template lookupObject<volScalarField>("omega") );

就会找不到，估计是结构问题，难道只能重写了。

有些站友提出问题之后，其他站友们给了追问或者讨论之后，就没有消息了.....
方案能不能用，倒是回复一下啊
完全不能理解这种人提出问题就再也不管的人，不靠谱

湍流普朗特数.png
请问有同志接触过S-A模型中将湍流普朗特数改到0.4以下的情况或者有提及此操作的文献吗？烦请赐教，谢谢！

感谢大神！后期操作是是指什么？我的运动的对象壁面的编号是23，但是如果我运动部分有多个部分该怎么办，就是运动几何体有多个避免边界？

替代文字

甲烷与空气燃烧反应，使用10种组分+10个基元反应，首先使用部分预混燃烧模型计算稳定后，更改为EDC模型，计算过程出现dasac failure at temperature，如下图所示：
捕获.PNG
请问这是什么原因

想开个贴讨论下多相流界面追踪类方法目前的前沿水平，例如能算到的最高参数，各类方法都有什么限制等等。目前多相流算法里，计算界面追踪的多相动力学，比如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会在薄液层处由于数值格式产生人工破碎的非物理情况，需要单独使用人工识别薄液层再稳定的方式，对不规则碰撞并不实用。

Data-Create Zone-Rectangular提取出想截取的位置
d2469623-64c0-45c3-952d-a6339a3b560b-image.png
但是提取出来的Zone没有数据，Data-Interpolate根据原数据对创建的Zone进行插值就可以
abd5cf18-caaf-4a49-a528-3b9afd696f35-image.png
这样就得到一个名为Rectangular zone的Zone
758ce629-1262-40bf-83b2-1477d15a474d-image.png

@李东岳 在 VOF方法中如何准确计算自由表面曲率？ 中说：

这个方法在函数calculateK()里面，可直接调用。计算结果在3232网格上能得到光滑的结果，但不对。在6464网格上出现了数值振荡。

VOF结果在网格比较细的时候结果不对？这个结论很有挑战性

李老师，我用interIsoFoam算了一下，还是不对。在32$\times$32的网格上可以得到光滑的结果，但是和解析解不一样；在64$\times$64的网格上结果出现了振荡。interfoam也跑了一下，结果和interisofoam不一样，也不对。。。（这个算例的设置参照和解析解在主贴的参考文献[1]中）

算例在这里，大家有兴趣可以跑一下。要是算例设置有问题也可以指出来。另外，我是在v2006环境下跑的(这个链接可能需要注册坚果云，我没找到其他的上传方法:xinlei: )
capillaryWave

如果不考虑时间项
\begin{equation}
\int^{f+1}_ {f} \frac{\partial}{\partial x}\left(\alpha \frac{\partial \phi}{\partial x}\right)\mathrm{d}x=\left(\alpha \frac{\partial \phi}{\partial x}\right)_ {f+1}-\left(\alpha \frac{\partial \phi}{\partial x}\right)_ {f}
\end{equation}
把后面两项$\left(\alpha \frac{\partial \phi}{\partial x}\right)_ {f+1}-\left(\alpha \frac{\partial \phi}{\partial x}\right)_{f}$把所有的面连起来，看起来是守恒的，剩下的只有边界条件

c运行库太老了，新的接口，软件用了，系统还不认识。更新C库可以解决。
按说这库装哪都可以，家目录也行，只要配置好环境让软件找到就行。
实际操作比较复杂，这种操作个性化太强，没有手把手式的教程，自己配方案，只能正着学，不能倒着学。倒着一堆坑，磨死人