CFD中文网

用的windAroundBuildings算例，840万网格的二维模型，网格质量检查过了，计算域是14000*4000米。epsilonInlet 0.0003; // Cmu^0.75 * k^1.5 / L ; L =10

dimensions [0 2 -3 0 0 0 0];

internalField uniform $epsilonInlet;

boundaryField
{
INLET
{
type fixedValue;
value uniform $epsilonInlet;
}

OUTLET { type zeroGradient; } WALL { type epsilonWallFunction; value uniform $epsilonInlet; } frontAndBackPlanes { type empty; } #includeEtc "caseDicts/setConstraintTypes"

}
kInlet 1.5; // approx k = 1.5*(I*U)^2 ; I = 0.1

dimensions [0 2 -2 0 0 0 0];

internalField uniform $kInlet;

boundaryField
{
INLET
{
type fixedValue;
value uniform $kInlet;
}
OUTLET
{
type zeroGradient;
}

WALL { type kqRWallFunction; value uniform $kInlet; } frontAndBackPlanes { type empty; } #includeEtc "caseDicts/setConstraintTypes"

}
nut
internalField uniform 0;

boundaryField
{
INLET
{
type calculated;
value uniform 0;
}

OUTLET { type calculated; value uniform 0; } WALL { type nutUSpaldingWallFunction; value uniform 0; } frontAndBackPlanes { type empty; } #includeEtc "caseDicts/setConstraintTypes"

}
p
internalField uniform 0;

boundaryField
{
INLET
{
type zeroGradient;
}
OUTLET
{
type fixedValue;
value uniform 0;
p0 uniform 0;
}
WALL
{
type zeroGradient;
}
frontAndBackPlanes
{
type empty;
}

#includeEtc "caseDicts/setConstraintTypes"

}
Uinlet (10 0 0);

dimensions [0 1 -1 0 0 0 0];

internalField uniform (10 0 0);

boundaryField
{
INLET
{
type fixedValue;
value uniform $Uinlet;
}
OUTLET
{
type zeroGradient;
}

WALL { type noSlip; } frontAndBackPlanes { type empty; } #includeEtc "caseDicts/setConstraintTypes"

}
时间步长试过1 0.1 0.01 其他的都没动

Direct numerical simulation of reactor two-phase flows enabled by high performance computing

捕获.JPG

$\alpha$

@FluidGao 我不是啥前辈:chouchou: 高校，博士后。

@程迪 请问linux如何并行计算，在ubuntu下没有问题，但是在linux下并行，显示这样的错误，希望您指教

OpenFOAM-6(23) mpirun -np 2 interFoam parallel -------------------------------------------------------------------------- The value of the MCA parameter "plm_rsh_agent" was set to a path that could not be found: plm_rsh_agent: ssh : rsh Please either unset the parameter, or check that the path is correct

@李东岳 谢谢李老师。

感谢东岳老师，我尝试着改小了，问题得到了部分解决。

@东岳 谢谢！

图一梯度算子的分量$\frac{\p }{\p x}$，不是某个物理量的梯度分量

$\frac{\p }{\p x}$不是一个函数么，你要写一个函数？

@King_RVM 你好，使用starCCM绘制网格，然后导入到openFoam里面跑算例，计算结果会出现质量不好的情况么？请教一下~

@李东岳 哈哈哈，物理隔绝了:mihu:

@yfleung @荷兰后腰阿贾克斯 @波流力
大佬朋友，请问你们处理400万以上网格，速度很慢，无法出结果的问题，解决了吗？
感谢

DPM限制.PNG
如上图，帮助文档中在dpm模型下使用mrf、动网格、sliding mesh会导致颗粒追踪失去意义，请问在dpm模型中使用srf会不会导致颗粒追踪失去意义？或者说在fluent里是不是把srf归为mrf？

@yxg2018 谢谢，我将网格按照单元类型分离并将边界层网格也设定为刚体运动区域，计算依旧会报错，提示Divergence detected in AMG solver: x-momentum，我将迭代次数增加并改小了动量的松弛因子，计算还是发散。
微信截图_20200713155718.png

你说的把Q的插值变成比值，应该是Ω方法吧？？？
Liutex-Ω的话好像比较麻烦

@西湖冷月 你具体编程过吗？

能用钱解决的问题都不是问题 :haqi: :haqi: :haqi:

@李东岳 在 HyperMesh为STARCCM+生成面网格 中说：

感谢分享！有益于后来人

共同进步！

@东岳 那截面云图总是有小的变动，这个影响结果分析吗？谢谢老师～

DEM加速比很高
CFD好像一般

@李东岳 感谢李老师，我试试看能不能画出来