CFD中文网

在一篇博文里看到，对于LES模型，过于细密的网格可能会导致大的误差和非物理解，请问各位大佬有人知道这是什么原因导致的吗？

想问一下nutroughwallfunction下设置的几个参数具有什么样的物理意义啊？其中roughnessHeight应该是常用的粗糙度吧？另外设置的value，roughnessFactor，roughnessConstant分别代表什么含义呢？

等待赐教，谢谢

嗯嗯，矢量分量等值面有了，谢谢！

利用paraview后处理计算结果绘制等值面时，等值面contour by后面的选项只有标量可以选择，想求助一下有没有哪位大佬知道怎么绘制矢量分量的等值面，比如说vorticity_z的等值面，或者某一速度分量的等值面？在此先行谢过了

@random_ran 再次打扰了，我想做一下尾部三维的涡结构，可是为什么创建contour等值面时，contour by后面的选项只有压力P和nut，并没有后来vorticity命令计算出来的变量vorticity，vorticity只能以云图的形式出现，我需要怎样创建wake的结构呢？

@random_ran 嗯呢，我开始也算过最大库朗数1.3左右的，后来担心结果的稳定性，最近调小了，可是这计算效率实在太熬人了，在您的建议下我决定再将时间步长迈大一点，哈哈，等时间充裕可以研究一下时间步长的收敛性。
感谢回复！

@random_ran 您好，最近看了一部分您的毕业论文，也在应用OpenFOAM的LES模型，看您在论文中提到计算过程中保证局部库朗数小于1 ，是指计算过程中的最大库朗数吗？
目前我的计算将最大库朗数调节到了0.3左右，可是这样的时间步长计算起来实在太慢了，请问LES模型的平均库朗数和最大库朗数大概应控制在什么范围内计算就相对比较稳定呢？

@东岳 好的，明白了。谢谢李老师！

@东岳 老师，我对LES对3D的要求有些疑问，2D一般是只在某一平面上进行网格划分，OF在2D模拟中进行了纵向拉伸了一层网格，那LES理论下所要求的3D是必须要完整的需要模拟的几何模型，还是只需要在2D网格的基础上在纵向延伸两层及以上的网格就可以满足3D模拟了？

@东岳 李老师的意思是运动粘性系数是用动力粘性系数除以密度得到的，但设置过程中我没有看到动力粘性系数的存在，那我怎么通过设置运动粘性系数保证求解过程中流体的密度是我计算工况条件下需要的密度大小呢，这和我设的运动粘性系数=1.139×10^-6Pa*s有什么关系吗？麻烦李老师能解答一下

@ryo 不太明白为什么动力粘度就不需要单独设置密度了呢，那如果我需要设置成水的属性，密度为1000kg/m^3，运动粘性系数=1.139×10^-6Pa*s。除了修改属性设置文件transportProperties中的粘度，我怎么保证密度也是我想要的1000kg/m^3呢？希望能指导一二，跪谢

在simpleFoam，pisoFoam等不可压缩求解器中的流体属性设置文件transportProperties中发现只有流体粘性的设置，那流体密度怎么给定呢，比如说我需要设置成水的属性，密度为1000kg/m^3，运动粘性系数=1.139×10^-6Pa*s。我该如何设置密度呢？

@赵一铭 还是不太明白，那动力粘度的值在哪呢，怎么根据动力粘度和运动粘度设置流体的密度值呢？希望能指导一二

想讲提取的自由液面波动的曲线数据输出，利用plot date工具，其中只涉及自由表面波动的问题能得到理想的曲线，如下面图1和2，可是当自由液面出现气泡时，曲线变得波动起伏，不能按理想曲线排列，如下面图3和4，想知道该用哪个filter怎样改变设置条件能得到云图中那样理想的曲线数据，希望懂的大神能指点一二，谢谢！

@东岳 奥奥，我明白李老师的意思了，用您的方法也得到了实心的圆截面rotationalExtrusion2，谢谢！可是我还有个问题想请教一下，就是我通过slice或者由slice1变换之后的transform2得到的rotationalExtrusion1就应该是个空心的圆柱体吗？还是我哪里设置有问题出现了错误？

@东岳 完整的圆柱计算域有了，可是为什么只有这个圆柱体的外表面，做这个rotationalExtrusion的截面只获得了一条圆形曲线，这是什么原因呢，麻烦李东岳老师帮忙解答一下，谢谢

@wwzhao 嗯，镜像成功了，谢谢，还想问一下，计算轴对称模型后处理能实现显示整个轴对称物体吗？

@东岳 找到了，谢谢李东岳老师！

@wwzhao 链接不存在了，能再重新分享一下吗，谢谢

kEqnCoeffs
{
    Ce              1.048;
    Ck              0.094;
}


Reading g

Reading hRef
Calculating field g.h

No MRF models present

No finite volume options present

GAMGPCG:  Solving for pcorr, Initial residual = 1, Final residual = 8.8753743e-06, No Iterations 8
GAMGPCG:  Solving for pcorr, Initial residual = 0.024255991, Final residual = 2.1939681e-06, No Iterations 3
time step continuity errors : sum local = 4.7008361e-15, global = -6.386216e-17, cumulative = -6.386216e-17
Courant Number mean: 0.001130128 max: 63.68082

Starting time loop

Courant Number mean: 0.001130128 max: 63.68082
Interface Courant Number mean: 7.3912381e-05 max: 9.8596523
Time = 0.0021

PIMPLE: iteration 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.50002593  Min(alpha.water) = -2.5864875e-05  Max(alpha.water) = 1.0000001
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.50002492  Min(alpha.water) = -0.010606339  Max(alpha.water) = 1.0000001
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.50002392  Min(alpha.water) = -0.0074242149  Max(alpha.water) = 1.0000001
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.50002292  Min(alpha.water) = -0.0078825052  Max(alpha.water) = 1.0000001
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::DICPreconditioner::calcReciprocalD(Foam::Field<double>&, Foam::lduMatrix const&) at ??:?
#4  Foam::DICSmoother::DICSmoother(Foam::word const&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&) at ??:?
#5  Foam::lduMatrix::smoother::addsymMatrixConstructorToTable<Foam::DICSmoother>::New(Foam::word const&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&) at ??:?
#6  Foam::lduMatrix::smoother::New(Foam::word const&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::dictionary const&) at ??:?
#7  Foam::GAMGSolver::initVcycle(Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::lduMatrix::smoother>&, Foam::Field<double>&, Foam::Field<double>&) const at ??:?
#8  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#9  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
#10  Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
#11  ? at ??:?
#12  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#13  ? at ??:?
浮点数例外 (核心已转储)

求助，跪谢！