Skip to content

OpenFOAM

OpenFOAM交流区

4.8k 主题 28.5k 帖子
  • 关于postProcess probe的问题

    5
    5 帖子
    5k 浏览
    hurricane007H

    @范春永 成了!非常感谢!!

  • 3 帖子
    3k 浏览
    李东岳

    最近看了几个extend的问题,暂时没用过extend的求解器...

  • 9 帖子
    8k 浏览
    R

    @xpqiu 感谢回复!
    现在把inlet,outlet和本体分开转换成stl文件 然后snappyHexMesh了一下
    0_1499148821716_屏幕快照 2017-07-04 下午3.07.55.png

    圆柱是r=18mm h=20mm inlet尺寸是r=1.5mm,h=3mm这样够不够大?现在想按您说的把空气区域和喷口区域网格如何变的密集?该添加哪些代码。。 求指导一下
    这个是我现在的snappyHexMeshDic

    FoamFile
    {
    version 2.0;
    format ascii;
    class dictionary;
    object snappyHexMeshDict;
    }

    castellatedMesh true; // make basic mesh ?
    snap true; // decide to snap back to surface ?
    addLayers true; // decide to add viscous layers ?

    geometry // Load in STL files here
    {
    inlet.stl {type triSurfaceMesh; name inlet;}
    outlet.stl {type triSurfaceMesh; name outlet;}
    wall.stl {type triSurfaceMesh; name wall;}
    volume.stl {type triSurfaceMesh; name volume;}
    refinementBox {type searchableBox; min (-0.03 -0.03 -0.03); max ( 0.03 0.3 0.3);}
    };

    castellatedMeshControls
    {
    maxLocalCells 1000000; //max cells per CPU core
    maxGlobalCells 2000000; //max cells to use before mesh deletion step
    minRefinementCells 10; //was 0 - zero means no bad cells are allowed during refinement stages
    maxLoadUnbalance 0.10;
    nCellsBetweenLevels 1; // expansion factor between each high & low refinement zone

    // Explicit feature edge refinement // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ features // taken from STL from each .eMesh file created by "SurfaceFeatureExtract" command ( {file "inlet.eMesh"; level 2;} {file "outlet.eMesh"; level 2;} {file "wall.eMesh"; level 2;} ); // Surface based refinement // ~~~~~~~~~~~~~~~~~~~~~~~~ refinementSurfaces // Surface-wise min and max refinement level { inlet {level (0 0);}

    outlet {level (0 0);}
    wall {level (3 3);}
    }

    resolveFeatureAngle 80; // Resolve sharp angles // Default 30 refinementRegions // In descending levels of fine-ness {volume {mode distance; levels ((0.0006 3) (0.002 2) (0.01 1));}} // was ((0.001 4) (0.003 3) (0.01 2)) locationInMesh (0.01 0 0); //to decide which side of mesh to keep ** allowFreeStandingZoneFaces true;

    }

    // Settings for the snapping.
    snapControls
    {
    nSmoothPatch 3;
    tolerance 4.0;
    nSolveIter 30;
    nRelaxIter 5;
    nFeatureSnapIter 15; // default is 10

    // New settings from openfoam 2.2 onwards for SHMesh

    implicitFeatureSnap false; // default is false - detects without doing surfaceFeatureExtract
    explicitFeatureSnap true; // default is true
    multiRegionFeatureSnap false; // deafault is false - detects features between multiple surfaces

    }

    // Settings for the layer addition.
    addLayersControls //add the PATCH names from inside the STL file so STLpatchName_insideSTLName
    {
    relativeSizes false; // was true
    layers
    {
    wall
    {nSurfaceLayers 3;} // was 3
    }

    expansionRatio 1.3; finalLayerThickness 0.00016; //was 0.00016 minThickness 0.00008; //was 0.00008 nGrow 0; // was 1 // Advanced settings featureAngle 80; // was 70 //- When not to extrude surface. 0 is flat, 90 is right angle. nRelaxIter 3; //- Max# of snapping relaxation iter. Should stop before upon reaching a correct mesh. nSmoothSurfaceNormals 1; // Number of smoothing iterations of surface normals nSmoothNormals 3; // Number of smoothing iterations of interior mesh movement direction nSmoothThickness 10; // Smooth layer thickness over surface patches maxFaceThicknessRatio 0.5; // Stop layer growth on highly warped cells maxThicknessToMedialRatio 0.3; // Reduce layer growth where ratio thickness to medial distance is large minMedianAxisAngle 130; // Angle used to pick up medial axis points nBufferCellsNoExtrude 0; // Create buffer region for new layer terminations nLayerIter 50; // Overall max number of layer addition iterations

    }

    // Generic mesh quality settings. At any undoable phase these determine
    // where to undo.
    meshQualityControls
    {
    maxNonOrtho 65;
    maxBoundarySkewness 20;
    maxInternalSkewness 4;
    maxConcave 80;
    minFlatness 0.5;
    minVol 1e-13;
    minTetQuality 1e-9;
    minArea -1;
    minTwist 0.02;
    minDeterminant 0.001;
    minFaceWeight 0.02;
    minVolRatio 0.01;
    minTriangleTwist -1;

    // Advanced nSmoothScale 4; errorReduction 0.75;

    }

    // Advanced

    debug 0;

    // Merge tolerance. Is fraction of overall bounding box of initial mesh.
    // Note: the write tolerance needs to be higher than this.
    mergeTolerance 1E-6;

    // ************************************************************************* //

  • 判断alpha==零点几的语句怎么编写

    9
    9 帖子
    7k 浏览
    X

    ForAll 是一个宏,定义如下

    #define forAll (list, i ) for (Foam::label i=0; i<(list).size(); ++i)
  • 如何获得cell信息?

    10
    10 帖子
    10k 浏览
    Y

    @程迪 good stuff, thanks

  • 关于数值模拟离心泵过程中发散的一些问题

    15
    15 帖子
    11k 浏览
    C

    明白了,谢谢前辈

  • 3 帖子
    3k 浏览
    bfamqB

    @赵一铭 画完网格之后 输入paraFoam命令,点击apply闪退。现在已经没问题了,之前p和U边界没有改。感谢!

  • DPM和流体耦合的时候是否考虑了压力?

    4
    4 帖子
    4k 浏览
    李东岳

    0_1498740892208_oie_291452346Q7QP69q.png

    谢谢,今天还真遇到一个无压力的NS方程,有意思
    http://epubs.siam.org/doi/pdf/10.1137/S0036142901398040

  • C++中if语句效率问题

    2
    2 帖子
    2k 浏览
    Y

    @cfd-china 这个不太好说,可以测试一下。但是gcc编译器的话应该优化的很好,这可能并不会是主要瓶颈,毕竟矩阵求解才是最耗时的

  • 这种特别长的公式一般怎么处理?

    2
    2 帖子
    2k 浏览
    Y

    以latex的经验是可以手动拆分的

  • 关于MRFSimpleFoam求解离心泵的一些问题

    6
    6 帖子
    5k 浏览
    C

    多谢前辈,我正在用probes监测

  • 通量分裂在多相流中的应用小测验

    2
    2 帖子
    2k 浏览

    比MULES更好的算法?看起来很尖锐啊,目前有没有完全尖锐的算法?好像教材里面说都有耗散性?
    发JFM啊岳哥,带我第9作者。

  • 请教各位大神关于动网格的问题

    2
    2 帖子
    3k 浏览
    M

    @lv1995 不知道是否正确,innerdistance的意思是随着动边界一起运动的网格厚度,不变形,outerdistance是网格变形域的范围。仅供参考。

  • 8 帖子
    6k 浏览
    李东岳

    0_1498032046222_2017-06-21 15-59-41屏幕截图.png

    我的一生,都在和震荡作斗争。

  • OpenFOAM后处理问题?

    4
    4 帖子
    4k 浏览

    @小龙 openfoam后处理forcesIncompressible可以用吗?因为我没有找到这个的详细解释,可以给我说一下这个后处理计算指定patch上压力的原理是什么吗?

  • 9 帖子
    10k 浏览
    W

    @yhdthu

    调解后的时间步进的长度delta_t都要小于设置的delta_T

    调节后的时间步大小取决于controlDict中的maxCo值,maxCo的默认值为1。代码见src/finiteVolume/cfdTools/general/include/createTimeControls.H 以及
    src/finiteVolume/cfdTools/general/include/readTimeControls.H。

    在每个delta_t中都要跑一遍nOuterCorrectors(50)和nCorrectors (2)对不?

    没错。

    干吗不用piso算法自己手动调节时间步长

    可以通过增大maxCo来增大自动调节的时间步长,据说PIMPLE算法的Co数可以达到10以上,不过对于用DES或LES计算的湍流问题另当别论。

    如何给定nOuterCorrctors 的值呢?

    这个只能靠经验了,不过一般肯定用不着50。

    这句话怎么理解?什么叫high turbulent flow??

    当simulationType设置为RAS或LES时,求解的是考虑湍流粘度($\nu_t$)影响的RANS/LES方程,而湍流粘度取决于求解湍流输运方程后得到的各湍流系数。

    turbOnFinalIterOnly为false代表每个PIMPLE循环内求解多次湍流输运方程,每次PIMPLE循环求解RANS/LES方程时的湍流系数都用的是上个PIMPLE循环内得到的湍流系数。

    turbOnFinalIterOnly为true则代表每个PIMPLE循环内只求解一次湍流输运方程,每次PIMPLE循环求解RANS/LES方程时的湍流系数都用的是上个时间步得到的湍流系数。

    高雷诺数流动为了精确得到速度压力场,需要每个PIMPLE循环内都求解多次湍流输运方程。

  • interior

    已移动
    7
    7 帖子
    6k 浏览
    李东岳

    http://www.cfd-china.com/topic/862/多个网格合并问题/9

  • DPMFoam 能不能处理颗粒的沉积

    已移动
    2
    2 帖子
    2k 浏览
    hurricane007H

    似乎可以,用stick壁面?

  • 基于直接数值模拟的新曳力系数

    已移动
    3
    3 帖子
    3k 浏览
    hurricane007H

    话说,这个作者我在CFB11的时候还见过,和她还有张合影。

  • 湍流模型是否需要考虑空隙率

    已移动
    8
    8 帖子
    8k 浏览
    hurricane007H

    OpenFOAM 里面,比如DPMFoam 需要重新编译湍流模型,是这个意思么?