CFD中文网

x3d
用这个导出，导入到blender

@vbcwl 这里调g相当于在方程后面添加了一个源项吧

@李东岳 老师您好，我经常往返上海和江苏南通，经过苏通大桥，感觉现在城市的公路建设里，有个非常草率的问题就是，一条路堵了，那就加条路（并联流动，或者说是加大管径，增大“车”流量），但是实际有个关键问题，就是堵路的原因经常是一条公路源源不断的增加路口，以及多车道的增增减减，那就是相当于管道流动中的增加了沿程阻力产生了压损，即理解为渐缩渐扩，三通管道等模式，这些沿程阻力点就是人驾驶的车辆的不确定性产生的，（即车辆的行驶会影响到其他车辆的行驶），人的思考考虑超车减速，能否用涡粘模型来建立，最终这些交通的建设可以通过流体力学的优化来进行设计？
虽然现在增加了铁路桥，但是苏通还是堵。。。

@李东岳 之前算过一个case，是通过分段画网格的方法，把关注区域的体速度场按10e-5s的时间存了1000个样本，一个cas的数据就存了1T。然后求u_i(t)-1/nΣu_i得到脉动速度，算湍动能的。
后来发现有个好东西，在fluent里面有个这样的子模型直接求到亚格子应力湍动能的
对应openfoam里是kEqu的方法，可以直接得到湍动能的分布趋势图。

@小张同学 主要是你需要存储的数据比较多，这就和piv测量的数据一样的道理，看你计算多少物理时间和时间步，如果存不下的话，可以设置只存需要的点线面区域的速度。

@小张同学 要么用单方程TKE模型，或者保存1000+的速度场数据，用湍动能的公式计算啊

@如果随风 先转数据格式为矢量； 然后计算涡量；然后绘制iso-surface， ，最后在当前涡量值等值面iso-surface上设置在这个面上显示其他变量，帖子最上的图应该显示的是速度脉动吧，带正负的，可能被调过颜色的，蓝红跨度这么大。

@ztchen https://www.twblogs.net/a/5e50f79fbd9eee2117c018c1可以参考这个搞一波 lambda2

@如果随风 Analyze 然后Calculate Turbulent Functions，里面可以直接计算一些流体的参量。

@aimer 就是你上面写的这个式子，原文是哪里的 。

您好，我想问个题外话，这个时变边界条件的文章是哪里看到的啊，能不能发一下链接

不过我碰到一个问题，这样的计算的时候单核计算可以进行，但是并行会碰到“mpirun noticed that process rank 0 with pid on node exited on signal 8”这个问题。谷歌也没找到

我也是在研究这个的，我的想法是首先研究湍流，对单相流计算或者PIV测量得到一个速度场，拟合出脉动速度的方程f(x,y,z,t)，随后在设置的方腔流域采用LES计算，使用上面的方程来设置边界条件,但是现在算出来的流场我用湍动能谱来对比原来的流场的，不知道有没有别的流场检验方法。
这样的话可以仅仅研究湍流作用下的气泡。下面是我测试的一个湍流边界进口。

![5e0120de-e12c-47b9-b79f-e79e95220c42-image.png](/assets/uploads/files/1565439493774-5e0120de-e12c-47b9-b79f-e79e95220c42-image.png) IN    
{	
type     codedFixedValue;	
value     $internalField;	
redirectType 	velocityProfile;	
code	
#{	
fixedValueFvPatchVectorFieldmyPatc(*this);	
forAll(this->patch().Cf(),i)	
{	
double x= this->patch().Cf()[i].x();	
double y= this->patch().Cf()[i].y();	
double z= this->patch().Cf()[i].z();
double t= this->db().time().value;
//Info<<“current time step ==” << t << endl;	
double r= sqrt(y*y+z*z);	
double SIN = y/r;	
double CON = x/r;	
double ux = -50*r*SIN;	
double uy = 50*r*CON;	
double uz = (1-cos(r/0.0265*3.14/2))*10;	
myPatch[i]=vector(ux,uy,uz);	
}	
operator==(myPatch);	
#};