CFD中文网

0_1501229627502_figure1.PNG

学OpenFOAM至今一直水在某大佬的OpenFOAM千人群中，个人认为还是比较活跃的之一(毕竟也水到50多级了)。
但是今天早上一起床，发现我被群内某管理员踢掉了：
qq_pic_merged_1629032541800.jpg

我寻思我也没发什么违禁言论，也没做什么违法群规的事情啊，且前天我还在群里解答疑惑:134:
关于这个问题我私聊了一下管理员，结果发现我问完他直接给我删了？:136:
30bfc521-cae5-4276-8a6b-89c6da167351-Screenshot_20210815_210758_com.tencent.mobileqq.jpg
145f5c52-5a8e-47e9-9341-302bf79191c7-Screenshot_20210815_210832_com.tencent.mobileqq_e.jpg

我不是很懂这波操作，按理说踢人不应该有个规范吗，还是按管理员想踢谁就踢谁吗？

PS：不知道这里是否允许对这类问题进行讨论，侵权麻烦管理员大大直接删帖吧

@西湖冷月 我理解了。我是做的水锤问题。
我现实中的问题是，从关闭阀门到阀门完全关闭耗时4s。开始流速2m/s，管子内径25mm，管子长度150mm，一共画了1.5w个网格。
还有下面两个问题：
1、由于我画的是个长方形流动区域，我看bilibili上它设置的出口瞬间关闭(0.000001s)
可是我想让出口均匀的关闭，以加速度为-0.5m/s2关闭。不知道能不能这样设置。
2、关于计算参数设置：计算10000步，每步的时间把原来的0.000001s设置成0.0004s，即(10000*0.0004=4s)，这样是不是就代表算出来的是实际物理时间的4s的过程？但是不知道这样是否对，是否可行？
望指点

感谢分享 楼主好人 :xiexie: :xiexie:

前几天OpenFOAM基金会在群里面发了个消息，说OpenFOAM-9发布啦！然后一片死寂，一个人都没有，跟咱们中国的夸夸群差远了，于是基金会补上了一张图，就是这面这个

tumbleweed.gif

给我笑爆了..

@东岳 多谢提示，我再研究研究。

奇怪，我在tutorials里面allrun了一个例子，paraview的时候就打开不了，会有error，说是polymesh/faces里面符号错误。。
按理说，tutorials里面的应该没问题啊。。我找了很久的原因，发现可能是二进制的问题，所以我把controlDict里面的输出从Binary改成了ASCII，结果allrun之后就可以paraview打开看见了。。
好奇怪，这是什么原因

@东岳 如果能够通过代码让边界层网格按湍流模型和y+自行调整就好了

\begin{equation}
m_0=w_1+w_2
\end{equation}
\begin{equation}
m_1=z_1w_1\xi_1+z_2w_2\xi_2
\end{equation}
\begin{equation}
m_2=z_1^4w_1\xi_1^2+z_2^4w_2\xi_2^2
\end{equation}
\begin{equation}
m_3=z_1^9w_1\xi_1^3+z_2^9w_2\xi_2^3
\end{equation}
\begin{equation}
m_4=z_1^{16}w_1\xi_1^4+z_2^{16}w_2\xi_2^4
\end{equation}
\begin{equation}
m_5=z_1^{25}w_1\xi_1^5+z_2^{25}w_2\xi_2^5
\end{equation}
左边的$m$都已知，右边的6个未知数$z_1$, $z_2$, $w_1$, $w_2$, $\xi_1$, $\xi_2$需要求。

除了Newton-Raphson迭代方法？

interFoam是直接模拟，不能模化升力。
双流体模型是模型化后的模型，可以模化升力。
因此

升力的方程表达式中出现气相与液相速度，这两者速度与interFoam控制方程中的速度有什么关联呢？

二者毫无关系

@bestucan 好的，谢谢

“Error: Invalid result name”

就是保存的太频繁了，电脑忙懵了。也可能是你把所有的核全用上了。按说应该留两个给电脑喘喘气。
或者你的c盘满了？

请教大家一个VOF自适应问题，问题如下：
利用fluent(Ansys 2021 R1)中的VOF+自适应网格模拟双液滴碰撞（三维）。液滴在计算几步之后，液滴表面变的凹凸不平。在尝试本站同学“加密后再patch”的方法后，初始时刻液滴表面变的十分光滑，但计算几步之后，液滴表面依然凹凸不平，如图1所示。
a1431f41-b0c8-4365-8228-6349b1ea8bdd-image.png
自适应设置和自适应网格效果如图2所示
3fd32bf9-f414-425d-a40d-141af0a77e5d-image.png
其他详细设置如图3所示
5addc093-d9b4-4a41-9c0c-44e24ef3bb45-image.png
谢谢大家，恳请大家指点帮组！

实验室有个试验台，里面充满水，水流速度最大60m/s，无叶轮，最近要做单相流实验，想确定增压泵打多少压力水才不会空化。按照我的想法，100℃水温时，水的汽化压力为0.1MPa，要保证不空化，需要线速度最大的地方局部压力达到0.1MPa，按照伯努利方程，需要加上动压，这个时候发现系统压力居然要达到2.8MPa。仔细想想应该是不适用伯努利方程的，所以想问一下，有没有什么好的方法可以估算系统压力呢？

欧拉拉格朗日模型中离散型用ODE方程描述，追踪每个粒子，这和PBM模型是冲突的，PBM是数量密度函数的模型。

目前还没见到过文献讨论这个 :cheeky:

用fsiFoam计算湍流，将流体求解器由constentIcoFlow改为pisoFlow，再添加湍流相关文件后开始计算。计算时间设置为6秒，时间步为0.00005，刚开始计算时很容易收敛，当计算到2秒的时候，每一个时间步都要迭代到最大迭代步数才会进行下一个时间步的计算，有什么方法可以改进吗？以下是计算过程迭代截图、流体求解器设置、耦合求解器设置以及固体属性截图。迭代1.png 迭代2.png 流体求解器.png 耦合求解器.png 固体属性.png

Von Neumann Stability Analysis 在普通的CFD教材里面讲的太简单了
但是到复杂非线性问题（比如Euler方程）中分析时间推进格式的阻尼特性 (damping characteritics) 的时候就比较复杂
而且一般期刊文章中只给应用结果而且缺乏细节
它背后的逻辑是什么，我一直没有搞清楚。。。
希望有大神能推荐一些详细一点的教材，能把思路讲的比较透彻的

我好像发现什么了，这个结构是一样的，不过tutorial case里obj生成polyMesh中，有分组。。。

( aerofoil { type wall; inGroups List<word> 1(wall); nFaces 120; startFace 31760; } inlet { type patch; inGroups List<word> 1(freestream); nFaces 200; startFace 31880; } outlet { type patch; inGroups List<word> 1(freestream); nFaces 160; startFace 32080; } back { type empty; inGroups List<word> 1(empty); nFaces 16000; startFace 32240; } front { type empty; inGroups List<word> 1(empty); nFaces 16000; startFace 48240; } )

所以，/0下文件头上freestream不是类型名，是边界（组）名。。。
感觉这个地方有误导性，导致我以为文件结构不同。。。
建议改个组名！

终于注册成功了，感谢李老师！

@东岳 感谢李老师！我去找一下相关内容！