Skip to content
  • 0 赞同
    5 帖子
    9k 浏览
    李东岳

    我没太跟随数值传热学中压力修正方程的推导,目前还不能给出推导:tishizi:

    我自己的推导 http://dyfluid.com/docs/RhieChow.html

  • 将文件中的空格替换为逗号

    CFD彩虹条
    2
    0 赞同
    2 帖子
    3k 浏览
    bestucanB

    如果只是查看,emacs 有个 whitespace mode。
    空格、tab、换行之类的特殊字符都可以显示
    vim按说应该也有,没搜到

    Screenshot from 2021-05-22 15-16-48.png

  • Word里面你们怎么打公式

    Algorithm
    21
    0 赞同
    21 帖子
    37k 浏览
    S

    aurora,可以在word里用latex

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

    OpenFOAM
    6
    0 赞同
    6 帖子
    10k 浏览
    C

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

  • 0 赞同
    2 帖子
    3k 浏览
    zousiyuZ

    用 DEFINE_SOURCE 宏,根据需求挂载在 k 或者 epsilon 的方程后面,试试

  • Icem混合网格设置

    Meshy
    3
    0 赞同
    3 帖子
    6k 浏览
    李东岳

    你的图就是节点对齐。

    有文献里是中心六面体网格,边界层为四面体问题网格。

    你可以试试hex-core,在做体网格的时候,可以成圣hex-core

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

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

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

  • 0 赞同
    6 帖子
    9k 浏览

    @cccrrryyy 感谢大佬答疑解惑,也谢谢您的分享:140:

  • 0 赞同
    9 帖子
    16k 浏览
    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;

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

  • 0 赞同
    1 帖子
    3k 浏览

    想开个贴讨论下多相流界面追踪类方法目前的前沿水平,例如能算到的最高参数,各类方法都有什么限制等等。目前多相流算法里,计算界面追踪的多相动力学,比如droplet collision, droplet splashing, droplet impact on solid wall,这几个典型案例,目前有用VOF,level-set,phase field,LBM, SPH, 欧拉拉格朗日耦合其它模型的。这几个案例在不同Re和We下会产生许多种演化形态,实验都做出来了,但感觉计算还没跟上。

    调研过不下百篇的文献,基本上关键词相关的都翻了一遍,在文献中看到droplet collision大概最好的也就是Re=6500,We=1500,密度比DR=1000左右的,大致是这个级别,像droplet impact有篇08年的用TransAT算到了RE=11000,而最近16年左右的用levelset这些也算到Re=8000,We=1500左右。大多数都是采用了一些改进的数值技巧来稳定模拟,而且网格也都比较密,计算资源消耗也都比较大,估计都得上超算,当然和实验也都对的上,基本都是到了prompt splattering这样的形态。

    想讨论下大家都做过或者见过的这类模拟最高参数是多少? 或者往上算各类方法的限制在哪? 因为多相流在算到高RE和We的时候会数值发散,稳定性应该是主要限制,还有例如VOF和level set会在薄液层处由于数值格式产生人工破碎的非物理情况,需要单独使用人工识别薄液层再稳定的方式,对不规则碰撞并不实用。

  • Fluent meshing totorials

    Meshy
    6
    0 赞同
    6 帖子
    11k 浏览

    这个主要看cfd之道公众号流沙的教程就可以学会了(哈哈哈哈,你是我拿到注册码之后回复的第一个人)

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

    压力没有一个绝对值的话,必须需要给一个参考值。但是你是自己写的代码,不好debug

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

    目前想通过编写“.js脚本文件”在DesignModeler中创建三维几何模型,之后再划分网格用Fluent模拟。求教“.js脚本文件”的编写方法,教程资源。感谢!
    1.png

  • PANS model 在CFX/FLUENT中的应用问题

    Fluent
    1
    0 赞同
    1 帖子
    3k 浏览

    参考文献:Turbulent Flows Over a Backward Facing Step Simulated Using a Modified Partially Averaged Navier–Stokes Model

    文献中利用CFX研究PANS方法对二维后台阶流动模拟的影响,由于刚开始研究此方面内容,不知如何下手,希望各位大佬进行指导!十分感谢!!figure2.jpg figure1.jpg

  • fluent封闭形成流体区域

    Meshy
    1
    0 赞同
    1 帖子
    4k 浏览

    QQ浏览器截图20220113163626.png
    问题:2021r2版本,在封闭形成流体区域的过程中,无法选择面。有哪位前辈指点一下吗

  • 0 赞同
    3 帖子
    6k 浏览
    Prometheus10P

    @李东岳 好的~谢谢东岳老师:xinxin:

  • 0 赞同
    1 帖子
    3k 浏览
    B

    对简单的长方体进行网格划分,再检查网格质量时发现在长方体的周围有很多有问题的face产生,也不知原因,很疑惑,求大神指点,万分感谢!

    0T5OC4QVS4OWJO998JSUAU7.png

  • 0 赞同
    4 帖子
    7k 浏览
    A

    @bestucan
    感谢您的回复!“不同方向的流动计算差异就大”这个解释很直接明了,我突然就有点明白这个网格线出现的原因了:xinxin:

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

    @流动魅力 我比较好奇你的相间作用力的直径怎么给的

  • ipmitool 调节风扇阈值

    CFD彩虹条
    1
    0 赞同
    1 帖子
    2k 浏览
    李东岳
    ipmitool -I lan -U ADMIN -H 192.168.1.18 sensor thresh FANB lower 100 150 200 ipmitool -I lan -U ADMIN -H 192.168.1.18 sensor thresh FANB lower 0 0 0