在某指定网格加入源项，源项为随时间变化的三角函数

[VamosZhi]

上次的问题解决后，面临这个新问题，最近也找了三角函数的定义，但都是关于边界条件的定义，没有明确的类似案例。
总的公式仍然是：

只是公式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];

但是如何在定义三角函数源项以及在网格中添加该源项仍然难。如有了解，希望提示一下。谢谢！

[zym604]

@vamoszhi

我一开始没看明白，是因为你这个公式里面的符号跟一般书上的不一样。

1. 看上去你这里的u'反而是average量，U_i才是湍动量。是么？
2. 公式7的雷诺应力tau是不是少了一项，2/3Kdelta_ij?
   
3. 你的意思是只在某个cell上套用公式9，其他cell继续用公式8么？
4. 我对能量方程不是很熟，公式8里面有P和p'，这就是两个变量了，你是不是还用了一些别的closure law来表达P变量了？

[Cloud1990]

但是建议做个sourceMask，能量方程可以写成
\begin{equation}
LHS=-\frac{\partial p}{\partial t}\times (sourceMask-1)-S_c\times sourceMask
\end{equation}