@李东岳 李老师,yang2009这篇sci论文里面我没看到有壁面剪应力的相关设置。由于这部分内容是我师兄做的,我先依据硕士论文解释您提出的几个问题,具体的我还需要再查证之后再给您回复。
硕士论文中引用的入流条件也是出自这篇sci文献,有出入的是壁面函数的粗糙高度取值,硕士论文中的粗糙高度是依据文献(Blocken et al. 2007, CFD simulation of the atmospheric boundary layer: wall function problems)计算所得:
壁面剪应力的取值是依据硕士论文中的一些假定提出的,而非yang2009中相关参数。
H
Hope
@Hope
帖子
-
-
@李东岳 李老师您好,我们这边之前在fluent里面有验证过该入流条件,具体结果可参见刘安彬的硕士论文《平衡大气边界层的CFD数值模拟》。但是在模拟时上部流场的湍动能并无明显衰减,只是近壁面有衰减,并在施加壁面剪应力之后可有效改善。您可以考虑下顶面边界条件、壁面剪应力这两方面的设置。
下图模拟结果截取自上述硕士论文,未加任何修正
下图是施加壁面剪应力后的模拟结果
-
你是否比较Upwind格式下的残差曲线和linearUpwind 格式的差别呢?这两方法在OpenFOAM中的收敛性是不一样的。或者设置Porbes探针观察下数据采样点最后的波动情况。我觉的是收敛性影响的。有复杂地形或者建筑群的边界层计算时,OpenFOAM高阶阶方案的收敛性并不好,会影响计算精度@李东岳
-
请教各位老师两个问题(1)在采用有限体积法对方程进行离散之后,动量控制方程中左侧的惯性力项如何计算?
(2)在高湍流度的流动中,方程右侧的粘性力项很小,流体的受力状态如何理解?(左侧惯性力与右侧压力梯度的平衡吗?广义源项的加入如何影响流动的发展呢?) 谢谢~~~
-
@东岳 好的,我试一下。谢谢老师~
-
@东岳 嗯嗯,是在同一个文件夹0/下的
-
@东岳 李老师,我边界条件代码修改后如下:
INLET { type fixedProfile; profile tableFile; profileCoeffs { nHeaderLine 0; // Number of header lines refColumn 0; // Reference column index componentColumns (1 2 3); // Component column indices separator ","; // Optional (defaults to ",") mergeSeparators no; // Merge multiple separators file "inVelocity.txt";//"inletProfiles/inVelocity.csv"; outOfBounds clamp; // Optional out-of-bounds handling interpolationScheme linear; // Optional interpolation scheme } direction (0 1 0); origin 0; }
同时inVelocity.txt数据按照下面格式给出,但还是有同样的报错信息
2.31474,0,0 3.07832,0,0 3.44144,0,0 3.64768,0,0 3.78464,0,0 3.88646,0,0 3.96812,0,0 4.03721,0,0 4.09802,0,0 4.15319,0,0 4.20449,0,0 4.25314,0,0 4.30008,0,0 4.34601,0,0 4.39149,0,0 4.43702,0,0 4.483,0,0 4.5298,0,0
-
模拟二维大气边界层流动,入口添加ESDU的对数律剖面时,采用如下方式定义,
INLET { type fixedProfile; profile tableFile; profileCoeffs { nHeaderLine 0; // Number of header lines refColumn 1; // Reference column index componentColumns (1); // Component column indices separator " "; // Optional (defaults to ",") mergeSeparators no; // Merge multiple separators file "inVelocity.txt";//"inletProfiles/inVelocity.csv"; outOfBounds clamp; // Optional out-of-bounds handling interpolationScheme linear; // Optional interpolation scheme } direction (0 1 0); origin 0; }
运行
simpleFoam
时,提示如下错误
inVelocity.txt文件为0.899251 3.78464 3.88646 3.96812 4.03721 4.09802 4.15319 4.20449 4.25314 4.30008 4.34601 4.39149 4.43702 4.483 4.5298 4.57773 4.62706 4.67807 4.73099 4.78608 4.84357 4.90368 4.96665 5.03268 5.10199 5.17476 5.25119 5.33143 5.41563 5.50394 5.59645 5.69327 5.79447 5.90012 6.01023
不知道问题在哪里?还请各位指教
中性大气环境湍流动能的自保持 | 附有算例下载
中性大气环境湍流动能的自保持 | 附有算例下载
算例比较upwind与linearUpwind对结果的影响
请教各位老师关于动量控制方程在数值求解过程中的问题
边界条件tableFile读取错误问题
边界条件tableFile读取错误问题
边界条件tableFile读取错误问题
边界条件tableFile读取错误问题