Skip to content
  • OpenFOAM 长波辐射模拟

    OpenFOAM
    1
    0 赞同
    1 帖子
    1k 浏览
    S

    想请教一下大家openfoam中有涉及长波辐射模拟吗,如何实现呢?

  • 压缩?不可压缩?

    Fluent
    3
    0 赞同
    3 帖子
    6k 浏览
    X

    但是你这个速度不高,马赫数很小,应该是可以按照不可压缩的

  • 0 赞同
    11 帖子
    17k 浏览
    ZhangZhiPengZ

    @Vortex
    感谢学长^_^

  • 一个关于K方程SGS模型的问题:

    Algorithm
    5
    0 赞同
    5 帖子
    10k 浏览

    谢谢,那个我也看了(但没仔细看,只是特意找了扩散项)。应该也是和Menon一样的结论。因为Menon的文章截图比较方便,才引用他的

  • 后向台阶流动的入口速度给定问题

    OpenFOAM
    2
    0 赞同
    2 帖子
    5k 浏览
    A

    我算358

  • 搅拌桨叶附近压力异常

    OpenFOAM
    21
    0 赞同
    21 帖子
    25k 浏览
    1

    @青山 好的谢谢!

  • 双层壁面函数

    Fluent
    4
    0 赞同
    4 帖子
    7k 浏览

    @CFD中文网 壁面函数是否只针对物面的第一层网格?

  • snappyHexMesh的网格加密问题。

    Meshy
    17
    0 赞同
    17 帖子
    17k 浏览
    G

    感谢大佬,我解决了问题,但似乎又出现了新的问题。起因是我希望用snap来加密多孔介质,但是我不知道为什么突然多了stl中不存在的部分。如下图,第一张图是cellZone,第二个是用来生成cellZone的stl文件。我在最下面贴出了我的snap设置。这让我十分困惑,我不理解为什么会出现这种情况。

    image.png

    image.png
    /--------------------------------- C++ -----------------------------------
    | ========= | |
    | \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
    | \ / O peration | Version: 2.2.0 |
    | \ / A nd | Web: www.OpenFOAM.org |
    | \/ M anipulation | |
    *---------------------------------------------------------------------------*/
    FoamFile
    {
    version 2.0;
    format ascii;
    class dictionary;
    object snappyHexMeshDict;
    }
    // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

    // Which of the steps to run
    castellatedMesh true; // make basic mesh ?
    snap true; // decide to snap back to surface ?
    addLayers false; // decide to add viscous layers ?

    geometry // 录入STL文件
    {

    side_mesh_0_45.stl {type triSurfaceMesh; name side_0;} side_mesh_45_90.stl {type triSurfaceMesh; name side_1;} side_mesh_90_135.stl {type triSurfaceMesh; name side_2;} side_mesh_135_180.stl {type triSurfaceMesh; name side_3;} side_mesh_180_225.stl {type triSurfaceMesh; name side_4;} side_mesh_225_270.stl {type triSurfaceMesh; name side_5;} side_mesh_270_315.stl {type triSurfaceMesh; name side_6;} side_mesh_315_360.stl {type triSurfaceMesh; name side_7;} top_mesh.stl {type triSurfaceMesh; name top;} building.stl {type triSurfaceMesh; name building;} grassland.stl {type triSurfaceMesh; name grassland;} ground.stl {type triSurfaceMesh; name ground;} water.stl {type triSurfaceMesh; name water;} tree.stl {type triSurfaceMesh; name tree;} refinementBox_building { type searchableBox; min ( 720 -600 15); max ( 2250 860 280); }

    };

    castellatedMeshControls
    {
    maxLocalCells 2000000; // 每个 CPU 核心允许使用的最大网格单元数,避免某个核心负载过重
    maxGlobalCells 20000000; // 全局允许的最大网格单元数,在细化阶段如果超过这个数量,系统将触发删除操作,防止内存溢出
    maxRefinementIterations 5; // 有效地限制整体迭代次数
    minRefinementCells 10; // 提高触发终止的阈值 maxLoadUnbalance 0.10; // 允许的最大负载不平衡度,表示各个 CPU 核心之间的工作负载差异最多为 10%,确保负载均衡
    nCellsBetweenLevels 3; // 不同细化等级之间的网格扩展因子,设置为 1 表示高细化区域与低细化区域之间的过渡最平滑
    maxShellRefinementIter 10;
    // 边缘特征细化程度
    features
    (
    {file "building.eMesh"; level 2;}
    {file "tree.eMesh"; level 1;}
    );

    refinementSurfaces { building { level (2 3); patchInfo { type wall; inGroups (building); } } tree { // level (minLevel maxLevel) level (2 3); // 例:最少细 3 级、必要时再细 1 级 faceZone treeFZ; // 可省;若想留出面区可保留 cellZone treeZone; // 关键:让 snappy 直接生成 treeZone cellZoneInside inside; // inside = STL 内部全进 cellZone patchInfo{ type patch; inGroups (treeSurf); } } side_0 { level (0 0); patchInfo{ type patch; inGroups (side); } } side_1 { level (0 0); patchInfo{ type patch; inGroups (side); } } side_2 { level (0 0); patchInfo{ type patch; inGroups (side); } } side_3 { level (0 0); patchInfo{ type patch; inGroups (side); } } side_4 { level (0 0); patchInfo{ type patch; inGroups (side); } } side_5 { level (0 0); patchInfo{ type patch; inGroups (side); } } side_6 { level (0 0); patchInfo{ type patch; inGroups (side); } } side_7 { level (0 0); patchInfo{ type patch; inGroups (side); } } top { level (0 0); patchInfo { type patch; inGroups (top); } } grassland { level (0 1); patchInfo { type wall; inGroups (ground); } } ground { level (0 1); patchInfo { type wall; inGroups (ground); } } water { level (0 1); patchInfo { type wall; inGroups (ground); } } } resolveFeatureAngle 60; // 特征边缘角度,小于这个角度的边缘将被忽略 refinementRegions { refinementBox_building { mode inside; levels ((1 2)); } tree { mode inside; // STL 内部 levels ((2 3)); // 在 level 3–4 的基础上再加 1 级 } } // {wholeDomain {mode inside; levels ((0 0));}} // 细化区域,mode 为细化模式,levels 为细化等级 // {building {mode distance; levels ((1 2) (3 1));}} // 细化区域,mode 为细化模式,levels 为细化等级mode 为细化模式(包括inside对内部网格进行细化 outside distance根据距表面网格距离进行细化),levels 为细化等级 locationInMesh (100 100 100); //规定哪边网格是流域 allowFreeStandingZoneFaces true; // 允许自由区域面

    }

    snapControls
    {
    nSmoothPatch 3; // 平滑次数
    tolerance 0.5; // 两个网格点之间的最大距离
    nSolveIter 50; // 迭代次数
    nRelaxIter 8; // 松弛次数
    nFeatureSnapIter 10; // 特征边缘迭代次数
    implicitFeatureSnap true; // 隐式特征边缘
    explicitFeatureSnap true; // 显式特征边缘
    multiRegionFeatureSnap false; // 多区域特征边缘
    }

    addLayersControls
    {
    relativeSizes false; // 相对大小/绝对大小
    layers
    {
    building {nSurfaceLayers 1;}
    }

    expansionRatio 1; // 扩展比例 finalLayerThickness 0.3; // 最外层厚度 minThickness 0.001; // 最小厚度 nGrow 0; // 增长层数 // 高级设置: featureAngle 80; // 表面层生成的特征角度。0 表示平面,90 表示直角。如果几何表面的角度超过 80 度,则停止挤出层。 nRelaxIter 5; // 层生成过程中,最大允许的松弛迭代次数。用于平滑网格,防止扭曲或变形。 nSmoothSurfaceNormals 1; // 对表面法线的平滑处理次数,用于减少表面网格的尖锐突变。 nSmoothNormals 5; // 对网格内部运动方向的平滑处理迭代次数,确保网格平滑地向外扩展。 nSmoothThickness 10; // 平滑表面层厚度的迭代次数,数值越大,厚度分布越均匀。 maxFaceThicknessRatio 0.5; // 如果网格面太过扭曲(厚度超过此比值),则停止层生长。防止生成变形严重的网格。 maxThicknessToMedialRatio 0.3; // 层厚度与中线距离的比值,防止厚度过大时停止层的生长。 minMedianAxisAngle 130; // 中线轴的最小角度。如果角度过小,可能会跳过层生成。通常用于复杂几何体的细化。 nBufferCellsNoExtrude 0; // 用于新层终止时创建的缓冲区单元数,防止层过度扩展。 nLayerIter 50; // 添加层的最大迭代次数,限制网格生成过程中层生长的迭代次数

    }
    // 控制网格质量的设置。任何无法处理的阶段将根据这些设置回退操作。
    meshQualityControls
    {
    maxNonOrtho 65; // 允许的最大非正交角(度数),控制网格的正交性。
    maxBoundarySkewness 4; // 网格边界倾斜度的最大允许值。
    maxInternalSkewness 4; // 内部网格单元的最大倾斜度。
    maxConcave 80; // 允许的最大凹角(度数),控制网格单元的形状。
    maxFaceDiff 2;
    minFlatness 0.5; // 允许的最小平面度,确保网格面不过于扭曲。
    minVol 1e-13; // 允许的最小单元体积,避免生成过小的网格单元。
    minTetQuality 1e-8; // 四面体网格单元的最小质量。
    minArea -1; // 最小网格面面积,-1 表示忽略此限制。
    minTwist 0.02; // 网格单元的最小扭曲度,确保网格质量。
    minDeterminant 0.001; // 网格最小行列式值,衡量网格形状的质量。
    minFaceWeight 0.02; // 网格面权重的最小值,用于衡量网格面质量。
    minVolRatio 0.01; // 最小网格体积比,衡量相邻单元的体积差异。
    minTriangleTwist -1; // 最小三角形扭曲度,-1 表示忽略此限制。

    // 高级设置: nSmoothScale 5; // 控制网格平滑的比例,数值越大,网格越平整。 errorReduction 0.75; // 错误减少系数,控制每次迭代中如何减少网格问题。

    }

    // 高级设置:是否启用调试信息
    debug 0; // 设置为 0 表示关闭调试信息。

    // 网格合并容差。相邻网格单元如果在此容差范围内,将被合并。
    // 容差值是初始网格边界框尺寸的一个分数。
    mergeTolerance 1E-5; // 合并容差,值越大,网格越容易被合并。

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

  • Fluent中带有凝固的压力场问题

    Fluent
    6
    0 赞同
    6 帖子
    8k 浏览
    bestucanB

    @wangchangjun_neu 流场影响温度,温度影响凝固,凝固影响应力,应力分布不均会导致成孔,结构不稳固。仿真凝固就是为了这个,不能考虑这一过程,应该没啥存在的意义。……还是算没有吧。我从没接触过这种仿真,我说的都是我瞎猜现搜的
    :mihu:

  • 0 赞同
    5 帖子
    9k 浏览
    Q

    我就是这个情况,笔记本会出错,win10系统win7都会出错,至今还没在台式机上发现

  • 0 赞同
    3 帖子
    3k 浏览
    H

    @hurricane007 老铁,发个联系方式

  • interPlicFoam

    OpenFOAM
    91
    2 赞同
    91 帖子
    324k 浏览
    队长别开枪

    interPlicFoam求解器已经并入扩展包geometricVofExt,仓库地址:https://github.com/daidezhi/geometricVofExt

  • blockMesh楔形体网格划分

    OpenFOAM
    5
    0 赞同
    5 帖子
    5k 浏览
    J

    @李东岳 ```
    /--------------------------------- C++ -----------------------------------\

    \ / 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 extrudeProperties;
    }
    // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

    constructFrom patch;
    sourceCase "$FOAM_CASE";

    sourcePatches (front);
    exposedPatchName back;

    extrudeModel wedge;

    sectorCoeffs
    {
    axisPt (0 0 0);
    axis (0 -1 0);
    angle 5;
    }

    flipNormals false;
    mergeFaces false;

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

    老师,这是我的extrudeMeshDict文件,我用这个文件进行拉伸,它提示我的wedge不是一个平面,为什么会这样呢?

    /--------------------------------- C++ -----------------------------------\

    \ / 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 blockMeshDict;
    }
    // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

    convertToMeters 0.001;

    vertices
    (
    (0 0 0)//0
    (7 0 0)//1
    (10 0 0)//2
    (1010 0 0)//3
    (2010 0 0)//4

    (0 5 0)//5 (5 5 0)//6 (10 5 0)//7 (1010 5 0)//8 (2010 5 0)//9 (0 40 0)//10 (7 40 0)//11 (10 40 0)//12 (1010 40 0)//13 (2010 40 0)//14 (0 1640 0)//15 (7 1640 0)//16 (10 1640 0)//17 (0 0 10)//0 (7 0 10)//1 (10 0 10)//2 (1010 0 10)//3 (2010 0 10)//4 (0 5 10)//5 (7 5 10)//6 (10 5 10)//7 (1010 5 10)//8 (2010 5 10)//9 (0 40 10)//10 (7 40 10)//11 (10 40 10)//12 (1010 40 10)//13 (2010 40 10)//14 (0 1640 10)//15 (7 1640 10)//16 (10 1640 10)//17

    );

    blocks
    (
    hex (0 1 6 5 18 19 24 23) (20 50 1) simpleGrading (1 ((1 1 50)(0 0 0)) 1)//1

    hex (1 2 7 6 19 20 25 24) (30 50 1) simpleGrading (((1 1 0.05)(0 0 0)) ((1 1 50)(0 0 0)) 1)//2 hex (2 3 8 7 20 21 26 25) (300 50 1) simpleGrading (((1 1 1500)(0 0 0)) ((1 1 50)(0 0 0)) 1)//3 hex (3 4 9 8 21 22 27 26) (40 50 1) simpleGrading (1 ((1 1 50)(0 0 0)) 1)//4 hex (5 6 11 10 23 24 29 28) (20 100 1) simpleGrading (1 1 1)//5 hex (6 7 12 11 24 25 30 29) (30 100 1) simpleGrading (((1 1 0.05)(0 0 0)) 1 1)//6 hex (7 8 13 12 25 26 31 30) (300 100 1) simpleGrading (((1 1 1500)(0 0 0)) 1 1)//7 hex (8 9 14 13 26 27 32 31) (40 100 1) simpleGrading (1 1 1)//8 hex (10 11 16 15 28 29 34 33) (20 400 1) simpleGrading (1 ((1 1 50)(0 0 0)) 1)//9 hex (11 12 17 16 29 30 35 34) (30 400 1) simpleGrading (((1 1 0.05)(0 0 0)) ((1 1 50)(0 0 0)) 1)//10

    );

    boundary
    (
    inlet
    {
    type patch;
    faces
    (
    (15 33 34 16)
    (16 34 35 17)
    );
    }

    outlet { type patch; faces ( (4 9 27 22) (9 14 32 27) ); } upper { type symmetry; faces ( (12 30 31 13) (13 31 32 14) ); } axis { type symmetry; faces ( (0 18 23 5) (5 23 28 10) (10 28 33 15) ); } bottom { type wall; faces ( (0 1 19 18) (1 2 20 19) (2 3 21 20) (3 4 22 21) ); } wall { type wall; faces ( (12 17 35 30) ); } back { type symmetry; faces ( (0 5 6 1) (1 6 7 2) (2 7 8 3) (3 8 9 4) (5 10 11 6) (6 11 12 7) (7 12 13 8) (8 13 14 9) (10 15 16 11) (11 16 17 12) ); } front { type symmetry; faces ( (18 19 24 23) (19 20 25 24) (20 21 26 25) (21 22 27 26) (23 24 29 28) (24 25 30 29) (25 26 31 30) (26 27 32 31) (28 29 34 33) (29 30 35 34) ); }

    );

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

    这是我的blockMesh文件,根据你之前的提示,我看那个算例里面是先定义了一个矩形出来,然后进行拉伸的,想仿照人家的来做
  • 关于snGrad格式里面的corrected格式

    Algorithm
    3
    0 赞同
    3 帖子
    7k 浏览
    金石为开

    @cfd-china 我是在看红皮书里面介绍的,这样的话我看看论文

  • fluent中使用的是同位网格吗?

    Fluent
    2
    0 赞同
    2 帖子
    3k 浏览
    L

    fluent使用的是交错网格加重叠网格技术

  • 自定义函数

    Fluent
    2
    0 赞同
    2 帖子
    3k 浏览
    R

    你看看是不是有U=0的点,specify时强制要求选中全部zone,可能是固体边界也被计入,导致无法计算。

  • 0 赞同
    4 帖子
    5k 浏览
    李东岳

    @zhaoym 老铁,后来这个做的怎么样了,有进展么

  • 快速多极子展开

    Algorithm
    5
    0 赞同
    5 帖子
    11k 浏览
    Elibathe康E

    @Wayne 对,电磁场里面会用到

  • 0 赞同
    2 帖子
    3k 浏览
    学流体的小明

    并行计算时候,每个进程掌握的网格就固定那些。
    举例b点只存在于processor x1,不在processor x2当中,processor x2自然无法获取b点的信息。
    要想processor x2获取b点信息,可以用reduce函数,也就是进程之间通讯一下。
    首先初始化intU为零向量,然后赋值,再reduce。

    point b (1,1,1); label cellID = mesh.findCell(b); vector intU(0,0,0); intU = U[cellID]; reduce(U,sumOp<vector>());

    具体代码你再看一下,我这个是思路,代码可能哪里不对。

  • 论坛搜索功能大幅度改善

    CFD彩虹条
    1
    0 赞同
    1 帖子
    2k 浏览
    李东岳

    不知道大家感受到没有。换了一个开源的搜索引擎,叫meilisearch。简直屌爆了。。。。。

    :146: :146: :146: