一个2D的算例更好。可以调节进口速度,来测试是否是正常的。
1D算例由于通量守恒,不好测试
administrators
负责CFD中文网站务
私有
帖子
-
activePressureForceBaffleVelocity/activeBaffleVelocity边界条件如何使用 -
关于求解器编译报错估计你下载的不是v11的版本。你要卸载互相适配的版本。
-
activePressureForceBaffleVelocity/activeBaffleVelocity边界条件如何使用感谢分享!!
-
关于求解器编译报错这个求解器不适用于openfoam11
-
移植sprayFoam液相功能屌爆了老铁,聊天记录已收藏,明天就发JFM
-
activePressureForceBaffleVelocity/activeBaffleVelocity边界条件如何使用测试这个边界条件,最好使用一个基础的求解器比如pisoFoam,使用自带的2D pitzdaily算例来单核运行测试,测试通过后将其衍生到你的算例。我在近期可以顺着我的路子弄一个。不过如果你可以弄一个上传的话,可以加快这个进程。
-
activePressureForceBaffleVelocity/activeBaffleVelocity边界条件如何使用你有测试pisoFoam这种单相流么
-
activePressureForceBaffleVelocity/activeBaffleVelocity边界条件如何使用https://cfd-china.com/topic/8018
算例发上来我看看
-
bug in constrainPressureconstrainPressure主要是更新压力第二类边界条件,公式如下:
\begin{equation}
\left( \nabla p_{rgh} \right)_f \cdot\bfn_f=
\frac{\left(\mathbf{HbyA}_f^{*}- \frac{1}{{{A^n_{\mathrm{P},f}}}}(\bfg\cdot\bfh\nabla\rho)_f - \mathbf{U}_f \right)\cdot\bfS_f}
{
|\bfS_f|
\frac{1}{{{A^n_{\mathrm{P},f}}}}
}
\end{equation}
上述公式与代码并不一致。在OpenFOAM中,constrainPressure为forAll(pBf, patchi) { if (isA<fixedFluxPressureFvPatchScalarField>(pBf[patchi])) { refCast<fixedFluxPressureFvPatchScalarField> ( pBf[patchi] ).updateCoeffs ( ( phiHbyABf[patchi] - rho.boundaryField()[patchi] *MRF.relative(SfBf[patchi] & UBf[patchi], patchi) ) /(magSfBf[patchi]*rhorAUBf[patchi]) ); } }
多乘了一个密度。应该改为:
forAll(pBf, patchi) { if (isA<fixedFluxPressureFvPatchScalarField>(pBf[patchi])) { refCast<fixedFluxPressureFvPatchScalarField> ( pBf[patchi] ).updateCoeffs ( ( phiHbyABf[patchi] - MRF.relative(SfBf[patchi] & UBf[patchi], patchi) ) /(magSfBf[patchi]*rhorAUBf[patchi]/rho.boundaryField()[patchi]) ); } }
在非常老的OpenFOAM版本中,看起来是正确的,与公式一致
setSnGrad<fixedFluxPressureFvPatchScalarField> ( p_rgh.boundaryField(), ( phiHbyA.boundaryField() - fvOptions.relative(mesh.Sf().boundaryField() & U.boundaryField()) )/(mesh.magSf().boundaryField()*rAUf.boundaryField()) );
-
内部面及其边界条件修改activeBaffleVelocity
在我那个笔记里面有写