CFD中文网

@李东岳 感谢东岳老师的回复！

这个潘工的联系方式给我呀，500块买不了上当啊

上传一下我的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; // ************************************************************************* //

@HITSC30 icem不是很熟练，我去试一下

各位大佬，我要模拟一个圆柱体，上表面逐步加载，上面表面没有流入和流出。水从下表面流出，算出来的结果很奇怪，要求计算10S，前3S都算正常，水向下流，但算到5S就停止了，最后一步水反向流出了。以下是流速的计算结果，帮忙看看，边界条件设置是否有问题，为什么结果怪怪的，谢谢。

circle.jpg

circle2.jpg

压力边界条件
inlet
{
type fixedValue;
value uniform 0.01;
}

pipe
{
type zeroGradient;
}
outlet
{
type zeroGradient;

}

流速边界条件
inlet
{
type zeroGradient;

}

pipe
{
type noSlip;
}
outlet
{
type inletOutlet;
inletValue uniform (0 0 0.1);
value uniform (0 0 0.1);
}

thinkpad t480s i7 8550u 16G内存
平时写论文、做ppt用，也可以开虚拟机改求解器算算小规模问题，真要数值计算还得远程服务器:140:

@李东岳 好的。以前在CFD-DEM的时候要几千个颗粒，DNS的话应该得缩短一下管道长度，算少一点颗粒量。

兄弟，这个是树木简化为多孔介质吧，这个比较简单，不要用UDF就能实现，直接采用多孔介质模型就可以了，就是需要转化一下，多孔介质模型所填的阻力系数=2*LAD*Cd。另外，兄弟，我也是在做这个，我想在多孔介质区域加一个附加项，考虑树木对例子的吸附，会吗？欢迎加我讨论QQ1225220118

fluent初始化显示mass diffusivity：invalid diffusivity please check material。用的组分输运涡耗散模型

image.png

@李东岳 感谢东岳老师解答

As an application of this geometric formulation we prove that the solution map of the incompressible Euler equation, associating intial data in Hs(Rn) to the corresponding solution at time t > 0, is nowhere locally uniformly continuous and nowhere differentiable. by Hasan Inci (2013)

这个是自己写的求解器，所以有可能差别确实很大

icem网格划分导入fluent时检查失败，显示Error: interface zone 3 has two adjacent cell zones.
4f68eb65-132e-4314-b0f6-6964dc466673-image.png

@Amadeus 好啊好啊！加个微信吧，我微信号：wx-200115qsw:146:

自问自答一下
C-AMG和SA-AMG的区别在于粗化方案不同以及插值矩阵的构造是不一样的，C-AMG基于细网格生成粗网格，可以理解为更精细，因此计算收敛性会更好；SA-AMG根据聚类的思想来的，所以在setup阶段会大大缩减时间，但是收敛性不如C-AMG；而两者都比较鲁棒性，因此并不是只针对某种方法有效。
而后一个问题目前调研到的是，AMG在商软中均可以独自使用或者作为预条件使用；而Fluent、COMSOL、STAR-CCM+这三个软件的官网都有对应的介绍，基本上作为预条件都是与CG/BiCGSTAB联用，其他的还有GMRES。
找了好久也算是把资料找到了

高雷诺数低雷诺数模型不是看流场的雷诺数，而是看壁面附近。

你2000用kEpsilon模型+壁面函数就可以。 如果壁面网格不好弄，那就需要高解析壁面网格+低雷诺数湍流模型。

@学流体的小明 好的，谢谢大佬，准备用starccm+试试画多面体网格，感觉六面体这么搞没得玩

@李东岳 感谢回复。我是想先写个简单的特定算例学习一下算法，格式边界条件啥的都直接写死，然后再慢慢做的通用一点，能搞出来肯定开源。