inletOutlet边界条件
-
inletOutlet主要过程是这个,这个边界条件主要调用mix边界条件,mix是
$$
\psi_f=w \psi_{\text {refValue }}+(1-w)\left(\psi_P+\nabla_f \psi y\right)
$$
其种的$\nabla_f \psi=0$(因为this->refGrad() = pTraits<Type>::zero;
),$w$通过inletOutlet计算出来:this->valueFraction() = 1.0 - pos(phip);
this->valueFraction()
就是$w$,phip
取决于往外出流量还是往里进,$w$的值要么是0,要么是1。所以$\psi_f$要么是固定值
$$
\psi_f= \psi_{\text {refValue }}
$$
要么是零法向梯度:
$$
\psi_f= \psi_P
$$
从代码上看倒是没问题。我一会找个算例debug一下