CFD中文网

Luo & Svendsen的破碎

d00a8e37-5666-4945-8964-2dd0e626a58e-image.png

f也是变量吗？ξmin怎么确定呢？积分好像太复杂了，不知道如何下手？有路过的大佬指点一下不？提前感谢了:chitang:

考虑最简单的动量交换，液相：
\begin{equation}
\frac{\rd \bfU_\rc}{\rd t}=-\bfA
\end{equation}
气相：
\begin{equation}
\frac{\rd \bfU_\rd}{\rd t}=\bfA
\end{equation}
乍一看，如果$\bfA$是负的，那么会导致$\bfU_\rd$向下走（为负），$\bfU_\rc$（为正）。但是这并不符合物理，考虑一个管子的颗粒，如果颗粒向下走，必然会带动空气同时向下走。但是方程缺不是这种体现，为什么呢？

实际上，上述方程并没有写完整，完整形式应该是这样：
液相：
\begin{equation}
\frac{\rd \bfU_\rc}{\rd t}=-\bfA=-\frac{1}{\tau}(\bfU_\rc-\bfU_\rd)
\end{equation}
气相：
\begin{equation}
\frac{\rd \bfU_\rd}{\rd t}=\bfA=\frac{1}{\tau}(\bfU_\rc-\bfU_\rd)
\end{equation}
如果开始的时候$\bfU_\rc=0$，则为

gv只能打开一个eps

wget -O /etc/yum.repos.d/epel.repo http://mirrors.aliyun.com/repo/epel-7.repo && yum -y install ntfs-3g

周期性边界条件主要是这俩个边界的通量一正一负，速度相同。

对称性边界主要是此边界通量为0，边界法向速度为0，只存在切向方向，对称性边界对于标量等于零法向梯度边界。

请教各位老师，外力驱动的流动，常用哪些算例做benchmark，比如Poiseuille流动，非定常的Womersley流动，还有吗？

比如调成2的话，这个时候的时间步长是多少？

你的时间步的2倍。

根据这幅图的说法，隐式耦合计算时，非稳态的courant我也可以设置的比1大吗？

隐式理论上可以使用任意大的时间步，因此没有库狼数限制。所以可以。

我的模型是计算一个泵的外特性以及泵的内部流场，这个泵由一个叶轮加一个导叶组成，叶轮叶片数为3，导叶叶片数为4，现在为了加快计算速度，采用了周期性边界，选取了一片叶轮加一片导叶来计算，所以这个流量的大小应该怎么设置？

看了一下你的模型，之所以你在相贯线处网格画不好，这主要问题在于你的几何模型——相贯线与底面相切，出现了三角区域。从而导致网格在该三角区域网格质量奇差无比。如下图所示

该问题导致出现了如下问题：

解决的办法是将底面拉伸一段距离，避免出现相贯线与底面相切即可，如下图所示
0_1536142189833_4ccfa31a-f8b8-4435-87b6-c322423cd7e1-image.png
如此，便可愉快的进行分块和关联
0_1536142304777_9fb3a5ee-3695-47f4-971f-61fe9c119d23-image.png
接着对其进行O型块划分
0_1536142342178_11fb8b28-9d22-441a-978f-e0f2d4ed1c7c-image.png
最终网格效果图如下
0_1536142436479_986e33bd-8587-4ee0-a010-5f914675e718-image.png

http://www.cfd-china.com/topic/905/几个有意思的研究方向

@下里巴人 流线图知道怎么画了，也能画出整体的矢量图，可是对于同一水平线的矢量图还不知道怎么画，能讲下么，谢谢

放点图有利于讨论..

@bestucan 还是那四个字：适者生存

如题。

看起来很复杂的样子.... :143:

@东岳 老师，我重新画网格后，decomposePar经常会报错，但是不影响后面计算，没有查到原因0_1537492983872_0403f9aa-ee80-4134-a228-ac6289ba58de-image.png

@alex 看到你的网格了，感谢分享，感觉划分的非常不错，好厉害！学习了！

这个看起来需要在你的局部区域做一个网格，然后mapField过去。

\begin{equation}
\zeta =x-ct,\eta=x+ct
\end{equation}
\begin{equation}
\frac{\p\eta}{\p x}=1,\frac{\p\zeta}{\p x}=1
\end{equation}
\begin{equation}
\frac{\p\eta}{\p t}=c,\frac{\p\zeta}{\p t}=-c
\end{equation}
对任意变量关于$\zeta,\eta$的函数
\begin{equation}
\frac{\p}{\p x}=\frac{\p}{\p\eta}\frac{\p\eta}{\p x}+\frac{\p}{\p\zeta}\frac{\p\zeta}{\p x}=\frac{\p}{\p\eta}+\frac{\p}{\p\zeta}
\end{equation}
\begin{equation}
\frac{\p}{\p t}=\frac{\p}{\p\eta}\frac{\p\eta}{\p t}+\frac{\p}{\p\zeta}\frac{\p\zeta}{\p t}=c\frac{\p}{\p\eta}-c\frac{\p}{\p\zeta}
\end{equation}
\begin{equation}
\frac{\p^2}{\p x^2}=\left(\frac{\p}{\p\eta}+\frac{\p}{\p\zeta}\right)\left(\frac{\p}{\p\eta}+\frac{\p}{\p\zeta}\right)=\frac{\p^2}{\p\eta^2}+\frac{\p^2}{\p\zeta^2}+2\frac{\p^2}{\p\eta\p\zeta}
\end{equation}
\begin{equation}
\frac{\p^2}{\p t^2}=\left(c\frac{\p}{\p\eta}-c\frac{\p}{\p\zeta}\right)\left(c\frac{\p}{\p\eta}-c\frac{\p}{\p\zeta}\right)=c^2\frac{\p^2}{\p\eta^2}+c^2\frac{\p^2}{\p\zeta^2}-2c^2\frac{\p^2}{\p\eta\p\zeta}
\end{equation}