如果用DPMFoam求解稀相流会怎么样？误差大么
-
@dzw05
“没错，在离散相的gravity中，表达式是massg_(1.0 - p.rhoc()/p.rho())，即考虑了连续相和离散相的密度比，这个是浮力的%(#c43d3d)[主要来源]。” 浮力还有哪些其它来源呢？
“如果单纯的将连续相动量方程中的重力g删除，而不是设为0，那么其实并没有影响粒子浮力的计算，因为粒子运动方程中根本没有出现连续相的压力。” 粒子运动方程中出现的连续相的压力，不是可以通过压力梯度力来体现出来了吗？这样一来，连续相动量方程中的耦合项或者压力项需做相应调整了。@东岳 您提出存在下面的两套方程，是等效的吗？
-
这俩套方程都是文献中看到的，可以参考本帖24楼。我自己的求解器中采用这种形式
\begin{equation}
m_\mathrm{dpm}\frac{\partial \mathbf{U}_ \mathrm{dpm}}{\partial t}=\mathbf{F}_ \mathrm{drag}+m_\mathrm{dpm}\mathbf{g}+\frac{m_\mathrm{dpm}}{\rho_\mathrm{dpm}}\nabla p_\mathrm{c}+\mathrm{Others}
\end{equation}俩种形式可以互相转换。依据$p=p_0+\rho_\rc \bfg\bfh$即可相互推导
-
\begin{equation}
\label{DPM} m_\dpm\frac{\rd \bfU_\dpm}{\rd t}=-\frac{m_\dpm}{\rho_\dpm}\nabla p+m_\dpm\bfg+\frac{m_\dpm}{\rho_\dpm}\frac{3}{4}\frac{C_\rD\rho_\rc}{d_\rd}\left|\bfU_\rc-\bfU_\dpm\right|\left(\bfU_\rc-\bfU_\dpm\right)+\mathrm{OtherForces}.
\end{equation}
伯努利：
\begin{equation}
\label{bnl} \nabla p=\rho_\rc\bfg
\end{equation}
代入：
\begin{equation}
\label{DPM3} m_\dpm\frac{\rd \bfU_\dpm}{\rd t}=m_\dpm\bfg\left(1-\frac{\rho_\rc}{\rho_\dpm}\right)+\frac{m_\dpm}{\rho_\dpm}\frac{3}{4}\frac{C_\rD\rho_\rc}{d_\rd}\left|\bfU_\rc-\bfU_\dpm\right|\left(\bfU_\rc-\bfU_\dpm\right)+\mathrm{OtherForces}.
\end{equation}
-
-
-
