在某指定网格加入源项,源项为随时间变化的三角函数
-
上次的问题解决后,面临这个新问题,最近也找了三角函数的定义,但都是关于边界条件的定义,没有明确的类似案例。
总的公式仍然是:只是公式8中的源项换了,换成了:
我已经写好程序去找到这个网格,label cellI; label cellIob; DynamicList<label> p1I; scalar value(1000); vector p1loc(1,0.5,0.5); forAll(mesh.cells(),cellI) { if(mag(U.mesh().C()[cellI]-p1loc)<=value) { cellIob=cellI; value=mag(U.mesh().C()[cellI]-p1loc); } } p1I.append(cellIob); Info<<"location"<<U.mesh().C()[cellIob];
但是如何在定义三角函数源项以及在网格中添加该源项仍然难。如有了解,希望提示一下。谢谢!
-
@vamoszhi
我一开始没看明白,是因为你这个公式里面的符号跟一般书上的不一样。- 看上去你这里的u'反而是average量,U_i才是湍动量。是么?
- 公式7的雷诺应力tau是不是少了一项,2/3Kdelta_ij?
- 你的意思是只在某个cell上套用公式9,其他cell继续用公式8么?
- 我对能量方程不是很熟,公式8里面有P和p',这就是两个变量了,你是不是还用了一些别的closure law来表达P变量了?
-
但是建议做个sourceMask,能量方程可以写成
\begin{equation}
LHS=-\frac{\partial p}{\partial t}\times (sourceMask-1)-S_c\times sourceMask
\end{equation}
论坛登录问题反馈可联系 li.dy@dyfluid.com