CFD中文网

@vbcwl 原来是因为硬盘里权限问题，不能出现带：的文件:136:

这不是用来推导离散Boltzmann方程的方法吗？Hermite expansion and Gauss-Hermite quadrature????

@Junren-Hou 在 ICEM划分网格 中说：

可以试试check block，有些时候其实没有问题，就它软件自己的毛病，然后block里最后有一个check block，可以让他自己检查一下。里面有好多选项，都试试，什么check/fixed block

回忆了一下，我的旋转区域用ICEM画的应该没有问题。出问题的应该是我画的外流场区域，一方面追求地面边界层，致使网格尺度过小，一方面为了平衡网格数量，过渡不均匀。猜测应该是外流区域的地面边界层设置导致了相关报错的产生。

我记得有些鼓泡塔的文献提过升力会影响计算收敛，尤其是气泡聚集时（云图出现大面积红色）残差会变得很大。这里建议先检查网格和时间步长是否匹配，适当减小时间步长，然后在UDF中添加
if (vof > 0.99)
cl = 0;
祝收敛！

@Wayne

我已联系你:sunglasses:

退Fortran用c++保平安，虽然没用过Fortran，但看知乎上专门论述Fortran的缺陷，发现用c++一下就好了。但是c++的学习周期长，东西多，也是问题。

@Wayne :confused:

各位大佬们好，我想在fsifoam流固耦合求解器中添加两种固体材料属性，请问有人做过相关的算例嘛？是否可以直接添加？还是需要修改求解器呢？

@xpqiu
.................
25T硬盘确实花了不少钱。:sunglasses:

icem也有这个功能，只是我没用过，

https://cfd-china.com/topic/6069/搅拌反应器icem结构网格划分方式/7

Screenshot from 2024-09-28 14-40-10.png

用的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 其他的都没动

c7a27198-c865-41aa-bc5f-ce28c85f32d3-1728486466827.jpg

请教各位大佬专家，立方体里边，挖掉一个圆柱体，对圆柱顶面进行切割，上半部分为流体入口边界，请问又什么好的网格划分思路吗？

单纯立方体里挖一个圆柱能够实现，但中间多切了一条线，不知道怎么操作比较好。

另外，有没有好的网格划分方面的教程资料推荐？

恳请大佬们指点。

@bestucan 感谢，我先试试。

@李东岳 还是受到李老师的启发想到这种操作:chouchou:

@李东岳 在 reconstructPar 一边组一边删除数据 中说：

最近在算几个LES，我勒个去数据太大了。大算例一般算完了只有几百个G，比如500G，reconstructPar重组的时候并不会删processor的数据，这将导致占据双倍的硬盘，直接干到1T。利用下面这个命令，可以重组一个时间步，删除一个时间步，对于硬盘吃紧的时候非常有用：

#!/bin/bash for i in $(foamListTimes -case processor0); do reconstructPar -time ${i} && rm -r processo*/${i} done

Pretty Simple :high: :high: :high:

@李东岳 都设置沉积边界条件的话，颗粒物算的就大很多了。它文献中是1微米的颗粒物稳态情况下充满整个房间，但我算的10微米，也是充满整个房间，跟文献中的数据差别很大。

advection:可以理解为平流，主要由物质浓度存在梯度引起的变化。convection：可以理解为对流，主要是由热量不均匀引起的垂直方向的对流。好像在流体力学中对这两者区别不是很大:xiabanle:

@bestucan 非常感谢老师指导:xinxin: 我回去试一试

上传一下我的snappyhexMeshDict文件，请各位老师帮忙看看```
code_text

========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | Website: https://openfoam.org \\ / A nd | Version: 10 \\/ M anipulation | \*---------------------------------------------------------------------------*/ FoamFile { format ascii; class dictionary; object snappyHexMeshDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // #includeEtc "caseDicts/mesh/generation/snappyHexMeshDict.cfg" castellatedMesh true; snap true; addLayers true; geometry { cylinder { type triSurfaceMesh; file "cylinder.stl"; } refinementBox { type searchableBox; min (0.0 0.0 0.08); max (1.2 0.2 0.22); } }; castellatedMeshControls { maxLocalCells 200000; maxGlobalCells 5000000; minRefinementCells 10; maxLoadUnbalance 0.10; nCellsBetweenLevels 6; features ( { file "cylinder.eMesh"; level 3; } ); refinementSurfaces { cylinder { level (3 3); patchInfo { type wall; } } } refinementRegions { refinementBox { mode inside; level 1; } } insidePoint (0.4 0.1 0.15); } snapControls { nSmoothPatch 3; tolerance 2.0; nSolveIter 30; nRelaxIter 5; nFeatureSnapIter 10; explicitFeatureSnap true; implicitFeatureSnap false; multiRegionFeatureSnap false; } addLayersControls { layers { "(cylinder_stl_face<stlunit=MM>|cylinder_stl_top<stlunit=MM>)" { nSurfaceLayers 10; } } relativeSizes true; expansionRatio 1.05; finalLayerThickness 0.8; minThickness 0.5; } meshQualityControls {} writeFlags ( // scalarLevels // layerSets // layerFields ); mergeTolerance 1e-6; // ************************************************************************* //

我碰巧最近在看这方面的工作，打算在OKS课进行植入，中文还对不上，但是应该就是楼上说的 转捩模型。我参考的是Menter的文章A One-Equation Local Correlation-Based Transition Model，他们的参考案例也很有意思，T3A系列，单独用层流、湍流都不行，只能用附加转变的模型。下图：

捕获.JPG