fvOptions中添加的源项需要考虑网格的体积吗?
-
大家好!
最近我在OpenFOAM中使用fvOptions设置了一个动量方程的源项,定义为
\begin{equation}
S=-C_z \times |U| \times U
\end{equation}
$C_z$是一个常数,$|U|$是速度的大小,$U$是速度,我在fvOptions中是直接写了eqn += -Cz*mag(U)
我看到源项的表达式应该是写成,是乘以网格体积了的
\begin{equation}
S=S_U \times V_P + S_P \times\ V_P \times T_P
\end{equation}
请问我在fvOptions中写的方式是需要乘以网格的体积吗 -
@allanZhong 2年过去了。。。是需要乘以体积。例如我刚写了个科氏力源项:
nergySource { type coded; selectionMode all; field U; codeInclude #{ #}; codeAddSup #{ Pout<< "**codeAddSup**" << endl; const volVectorField& U = eqn.psi(); vectorField& USource = eqn.source(); vector omega(0, 0, 5e-5); const DimensionedField<scalar, volMesh>& V = mesh().V(); USource += (2.0*omega)^U*V; #}; codeAddRhoSup #{ Pout<< "**codeAddRhoSup**" << endl; #}; codeAddAlphaRhoSup #{ Pout<< "**codeAddAlphaRhoSup**" << endl; #}; }