CFD中文网

我导收了一个湘潭大学的搅拌器的英文sci。实验+cfd，用openfoam搞的。我针对算例设置这一段提出了1个不痛不痒的问题，不能一个不提啊。目前已返回我导。估计明后天就返回编辑吧？哈哈。我导没有抵制，不过提了2页问题...anyway，问题不大。:sunglasses:

@bestucan 你是行家呀，厉害厉害，终于知道原因了

@东岳 用interDyMFoam算潮流能水轮机，没有仔细看，计算速度大概有不到10%的提高

http://dyfluid.com/icoFoam.html 另外 这个例子你看过么 最后那个扰流

@TINGHAIK 谢谢您，问题已解决

$\epsilon \rightarrow 0$ 并不是$\epsilon = 0$，否则不会有...if $u$ is smooth at $x_1$ and $x_2$...。在粘度趋向于很小的时候，不连续变成具备一定厚度的光滑解，同样承认有厚度的激波。所以
$$
\epsilon\int_\Omega\frac{\p}{\p u}\left(\frac{\p \eta}{\p u}\right)\left(\frac{\p u}{\p x}\right)^2\rd x\rd t \geq 0
$$
另外，
\begin{equation}
\int_{x_2}^{x_1}\left(\epsilon(\eta_q q_x)_x\rd x -\epsilon\eta(\eta_q)_q q_x^2\right)\rd x=\epsilon\left(\eta_q q_x|_{x=x_1}-\eta_q q_x|_{x=x_2}\right)-\epsilon\eta(\eta_q)_q q_x^2\Delta x
\end{equation}
考虑一个非常小的$\epsilon=1e-10$，在控制体内$\epsilon\left(\eta_q q_x|_{x=x_1}-\eta_q q_x|_{x=x_2}\right)\rightarrow 0$，$\epsilon\eta(\eta_q)_q q_x^2\Delta x$还是大于0.

直接的宏应该是没有的 你可以自己用UDF写个类似于level set的函数计算每个位置到界面的距离

牛顿内摩擦定律与stokes定律间有什么关系？二者计算壁面应力有何区别？

@搬运工 好的，谢谢回复:yes:

@东岳 Thin-Layer Navier—Stokes Approximation

苹果retina屏那个本也不错。就是能不能装win系统？

http://item.jd.com/2819597.html

@xpqiu
.................
25T硬盘确实花了不少钱。:sunglasses:

@肖艺 STAR CCM+
例子： https://www.cfd-schuck.de/en/cfd/anwendungsbeispiele/simulation-of-gear.html

@东岳 因为时间成本问题吧，现阶段还没有这方面的打算，之后应该会学。。。其实只是想问下这个求解器，如果简单方便入手的话想要了解一下。。。

@东岳 总压减小是因为速度变小，还是静压变小，或二者都有？无粘不可压的时候就是动能与压力势能的转化，是守恒的。

谢谢李老师的回复！因为雷诺应力不仅出现在动量方程里面，也出现在turbulent kinetic energy equation 与 turbulent frequency equation中，我的想法是如果我分别修改湍流模型的文件（kOmegaSSTBase）与动量方程的文件（linearEddyViscosity）就会要定义两次非线性项，所以我选择就在kOmegaSSTBase中更新turbulence->divDevRhoReff(U),在其中加入非线性项，参考的资料是OpenFoam里面的nonlineareEddyViscosity。由于这样修改需要在kOmegaSSTBase.C文件的constructors中加入eddyViscosity 也就是上图中蓝色后面的部分，我直接在原来的代码Turbulence后面加这一项会报错，想请问这怎么解决？非常感谢李老师耐心看我的问题，不知道我表述清楚没有？

@blgy 非常的感谢，你这个是在哪看到的啊？你这个图片也看不出来用什么软件，用什么方法划分的啊，仅仅是欣赏一下:ok:

@百善孝为先0626 可以的话还是想麻烦您把我的联系方式告知他，这边看能不能合作一下

@东岳 李博，我准备试试，之前写过fluent多面体的网格转换程序，不过这个发不了论文:mihu: