CFD中文网

天天净整这乱七八糟的

sad....

0_1460680809776_捕获.JPG

@李东岳 老师你好，我想请问下，设置了重力加速度g后，dynamicDict下不设置初速度，物体不能下落，这是哪里没设置好吗？

@LuoS 在 snappyHexMesh对复杂地形添加边界层报“浮点数例外 (核心已转储)”错误 中说：

教程中似乎有个建筑风的案例：
/Tutorials/inco

好的谢谢 我试一下

@cloud1990 在 VOF求解自由水面问题，求解空气区域吗？ 中说：

Image

Left: Correct jet shape predicted by TruVOF technique used in FLOW-3D. Right: Incorrect jet shape predicted by pseudo-VOF technique used by other CFD codes.

Image
Left: FLOW-3D‘s TruVOF technique predicts jet impingement on wall and some outflow. Right: Pseudo-VOF methods don’t predict realistic jetting of fluid on side walls.

从结果上看，flow3D 所谓的TruVOF是有问题的，air velocity 和　air resistance = 0?　这个可能比较适合粗网格计算，作为一个tricky的存在．但理论上，这个TruVOF是无法解决wind wave的问题．

此人在2005年左右淡出CFD-online.com，个人和他有一些私交，谈到了这个话题：为什么不活跃与论坛了？我们从您的帖子中学到了很多东西。他如下回答（翻译）：

活跃于论坛是非常耗费时间，许多问题都是重复的，越来越多的人想不劳而获的获取帮助，他们也不愿意为帮助他们的人提供一点好处。CFD不是我的业余爱好，CFD是我的职业，我需要靠CFD谋生。提供免费的服务不可持续。

是的，@hangsz 说的第二种和行业关系很大。
就是在计算机行业也是如此，化工行业的CS说是待遇20万+也是普遍。

:sunglasses:

@刘豪 可以的，可以设置Fluent 缓存位置

我上传到夸克网盘中了，请自行下载。我用夸克网盘给你分享了「ex15-2.zip」，点击链接或复制整段内容，打开「夸克APP」即可获取。
/~de573769eK~:/
链接：https://pan.quark.cn/s/e371aa252cc8
提取码：uaz2

@东岳 是我搞错了，向量混合积跟微分算子不一样

@李东岳 最近在潜水

@南旅 好的，感谢指导

都柏林大学发来贺电！

可以试试改为速度入口

黑圈那些edge，都不能关联，要取消关联

0_1519692297906_捕获.JPG
是否是这样？应力张量我可以理解。我的疑问是速度梯度在笛卡尔坐标系下变换之后直接就变成在极坐标系中了速度梯度了？

@李东岳 老师我想请问一下，设置了重力加速度g后，我将dynamicDict文件下将 velcity设为（0 0 0），物体就停着不动了，请问需要怎么设置一下呢。

欧拉欧拉那面由于已经在统计学方面已经处理过了，没有噪音。欧拉拉格朗日模拟这面目前有什么文章在处理么？找了找，还没看到。

As an alternative, Lagrangian methods ‘‘discretize” the density function into ‘‘parcels” that are simulated using Monte-Carlo methods. While quite accurate, as in any statistical approach, Lagrangian methods require a relatively large number of parcels to control statistical noise, and thus are computationally expensive.

之前有点问题，最终如下
\par $C_p$采用4次多项式分段拟合
\begin{equation}\label{equ:NS_Cp}
C_p(T)=a_1+a_2T+a_3T^2+a_4T^3+a_5T^4
\end{equation}
\par静焓
\begin{equation}\label{equ:NS_H}
H(T)=\int_{T_{0}}^{T_{x}}C_p(T)dT=a_1T+\frac{a_2}{2}T^2+\frac{a_3}{3}T^3+\frac{a_4}{4}T^4+\frac{a_5}{5}T^5+a_6
\end{equation}
\par熵
\begin{equation}\label{equ:NS_S}
S(T)=\int_{T_{0}}^{T_{x}}C_p(T)\frac{dT}{T}=a_1\ln{T}+a_2T+\frac{a_3}{2}T^2+\frac{a_4}{3}T^3+\frac{a_5}{4}T^4+a_7
\end{equation}
\par$\bullet$求解静温(已知总温和马赫数)
\begin{equation}\label{equNSUs}
U_s^2=2\left(H(T_{tot})-H(T_{sta})\right)
\end{equation}
\begin{equation}\label{equNSmach}
Mach^2=\frac{U^2}{\gamma(T)R_gT}=\frac{U^2}{\frac{C_p(T)}{C_p(T)-R_g}R_gT}
\end{equation}
\par由(\ref{equ:NSmach})和(\ref{equ:NSUs})得
\begin{equation}\label{equNSTsta}
T_{sta}=\frac{2\left(H(T_{tot})-H(T_{sta})\right)}{Mach^2\frac{C_p(T_{sta})R_g}{C_p(T_{sta})-R_g}}
\end{equation}
$\bullet$求解静压(已知总温、总压和静温)
\par由p等熵过程
\begin{equation}
ds = C_p(T)\frac{dT}{T} -R_g\frac{d p}{p}=0
\end{equation}
\par两边同时积分有
\begin{equation}
\int_{T_{tot}}^{T_{sta}}C_p(T)\frac{dT}{T} =\int_{p_{tot}}^{p_{sta}} R_g\frac{d p}{p}
\end{equation}
\par记
\begin{equation}
S(T_{x})=\int_{T_{0}}^{T_{x}}C_p(T)\frac{dT}{T}
\end{equation}
\par则
\begin{equation}
S(T_{sta}) - S(T_{tot}) = R_g\ln\frac{p_{sta}}{p_{tot}}
\end{equation}
\par那么
\begin{equation}\label{equNS_psta}
p_{sta}=p_{tot}e^{\left(\frac{S(T_{sta})-S(T_{tot})}{R_g}\right)}
\end{equation}

@jdzhang 你的颗粒会动么？不会动，为什么颗粒间还有距离？会动，那对着“多相流”查比对着“多孔介质”查更容易找到解决方案些。

这有一个画网格的：
https://www.bilibili.com/video/BV1eS4y1F7Wm/?vd_source=c01e8d4fff0db460c75b32e0d0191b24

这有些相关的：
https://www.researchgate.net/publication/308996206_Pore-scale_numerical_simulation_of_mud_erosion_in_the_subsea_sand-mud_alternate_layer_using_lattice_Boltzmann_method

https://ascelibrary.org/doi/10.1061/(ASCE)HY.1943-7900.0001489