关于cyclic boundary condition的一些问题

[明小米]

@wangfei9088 谢谢大佬，可是在跑的过程中还是有一些报错耶。

[wangfei9088]

@明小米 方便把你的U文件贴出来吗？

[明小米]

@wangfei9088 是这样的。

[wangfei9088]

@明小米
你的pi那一行应该是写错了，100和pi写反了。
正确的应该是：

	            field[faceI] = vector(U_0*sin(constant::mathematical::pi*100*t),0,0);

这句constant::mathematical::pi可以这么理解：
常数::数学::圆周率

[明小米]

@wangfei9088 请问大佬这是指u=sin100pi t吗？如果只用u=sint该如何表达呢？ 改了之后还是有报错耶

[wangfei9088]

@明小米
对。就是u=sin(100pit)的意思，是我对pitzDaily那个case自己改的。

1. u=sint可以这么写：

	        scalar U_0 = 1;
		scalar t = this->db().time().value();
	 
	        forAll(Cf, faceI)
	        {
	            field[faceI] = vector(U_0*sin(t),0,0);
	        }

2. 最后那个报错：浮点数例外（核心已转储），可能的原因是哪个位置把0作为分母了。

[明小米]

@wangfei9088 好的 报错的问题解决了 是由于初始速度设定为0了，修改为1就不报错了，非常感谢大佬！

[明小米]

@wangfei9088 大佬 如果要改称sinx的绝对值应该如何修改代码呢？这里的速度代码似乎与C++的表示有些不同

[wangfei9088]

@明小米
两种方法吧，第一种平方再开方，第二种就是用mag(sin(t))表示求模：

	        scalar U_0 = 1;
		scalar t = this->db().time().value();
	 
	        forAll(Cf, faceI)
	        {
	            field[faceI] = vector(U_0*mag(sin(t)),0,0);
	        }

[明小米]

@wangfei9088 好的，非常感谢大佬！