CFD中文网

好吧，既然没有人回答，我就把我的笨办法贴出来
foamToVTK-->paraview-->打开算例文件和粒子Cloud文件-->File/Save Data-->得到每个时间步每个粒子信息.csv-->用python提取有用信息

@OF九月九 nParticel 是每个parcel里面包含的颗粒数。。 主要还是你的颗粒总数。关键看你要模拟什么了，这个数字要符合实际不是说瞎改的。。比如你实际系统中就5w个颗粒 。你怎么改体积分数和质量分数都很小。。

@cfd-china 谢谢！我找到了，好像是主文件中的parcel.evolve()。

@zhanghan 客气:sunglasses:

@cfd-china :happy: THANK YOU!

0_1469972545332_foamToTecplot360编译.JPG 下载OF-2.3.0对应的tecplot360版本进行编译即可，我在OF-2.4.0版本安装没问题，你可以试试，我的理解是，OF-2.3.0后tecplot没有变

@李东岳 李老师，您好。volScalarField rAU(1.0/UEqn.A());这个重理论上分析应该是矢量才对，但是细节上的操作确实做成了体标量场。

template<class Type> Foam::tmp<Foam::volScalarField> Foam::fvMatrix<Type>::A() const { tmp<volScalarField> tAphi ( volScalarField::New ( "A("+psi_.name()+')', psi_.mesh(), dimensions_/psi_.dimensions()/dimVol, extrapolatedCalculatedFvPatchScalarField::typeName ) ); tAphi.ref().primitiveFieldRef() = D()/psi_.mesh().V(); tAphi.ref().correctBoundaryConditions(); return tAphi; }

$ A=\frac{A_p}{\Delta V}$
在计算系数Diag的过程中，将边界对主对角线的系数进行了平均更新：

template<class Type> Foam::tmp<Foam::scalarField> Foam::fvMatrix<Type>::D() const { tmp<scalarField> tdiag(new scalarField(diag())); //这里进行了平均处理并修正了主对角线系数。 addCmptAvBoundaryDiag(tdiag.ref()); return tdiag; }

为什么这样操作呢？这样会使得系数矩阵不一样了!

更新科氏力源项：$\mathbf{f}\times\bfU=2\Omega\sin \lambda \times\bfU$，其中$\mathbf{f}$表示科氏力参数，一些文献会直接给出$\mathbf{f}$的值。如果绕z旋转，那么$\mathbf{f}$在$z$方向存在值。在这种情况下，$\mathbf{f}\times\bfU$的值为：
\begin{equation}
\mathbf{f}\times\bfU=\left[
\begin{matrix}
f_2u_3-f_3u_2\\
f_3u_1-f_1u_3\\
f_1u_2-f_2u_1\
\end{matrix}
\right]=\left[
\begin{matrix}
-f_3u_2\\
f_3u_1\\
0\
\end{matrix}
\right]
\end{equation}
如果进一步考虑速度方向是x方向的，那么变为：
\begin{equation}
\mathbf{f}\times\bfU=\left[
\begin{matrix}
0\\
f_3u_1\\
0\
\end{matrix}
\right]
\end{equation}
可以看出如果地球绕z旋转，对于x方向的速度，科氏力导致y方向的偏离。

$\Omega$表示地球旋转速率，单位为[弧度$/\mathrm{s}$]=$1/\mathrm{s}$，$\lambda$表示维度[无单位]。因此$f$的单位为[$1/\mathrm{s}$]