CFD中文网

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

sed -i "s/oldString/newString/g" `grep oldString -rl /path`

@李东岳 谢谢您

@J 我们课题还是偏这个软件，谢谢您的回复

但是RC插值没有离散，因为它是显式计算源项 -- “速度插值+高阶压力导数项”。

为什么我这个推导的反了？应该是大于0 。

0_1529644803725_微信图片_20180622131807.jpg

可压缩PBE方程形式的比较

咦、之前写的这个研究方向已经发出去了

SuperBee好像就会锐化界面。

Yes, SUPERBEE加了下风格式，是叫下风么？Down-wind

哈哈，东岳大神，感谢你做了这个平台 ... 感觉学习CFD就是要不断地重复才能领悟其中的真谛，当然有通俗易懂的资料会更为便捷。

@mohui 多谢支持，CFD任重道远比较耐玩！只要有兴趣，是能研究一辈子的，不像王者荣耀，玩几年到头了 哈哈 :mihu: 资料以后会持续更新，精力有限，稍微慢点

请问像这种狭窄处要划分边界层，怎样处理能保证网格高质量呢

搜索功能可用了

@sgggggg 硕士论文我也想拜读一下，可以加个联系方式吗2132306535

张量标识预习测验

下面的方程取自各个不同的CFD教材。方程并没有采用代码的形式而是直接截图，这样更加直观的可以感受到文献里面方程这种多样性的书写方式。各位老师和同学需要能够拆分为分量形式。

0_1529193840362_1.JPG

0_1529193880752_2.JPG

0_1529194030583_3.JPG

0_1529194061311_4.JPG

0_1529194127887_5.JPG

0_1529195669196_6.JPG

0_1529195734307_7.JPG

0_1529195783721_8.JPG

0_1529196217099_9.JPG

0_1529196245498_10.JPG

0_1529196323392_11.JPG

0_1529196403565_12.JPG

0_1529196500098_13.JPG

0_1529197325295_14.JPG

VIV.zip 这是我配置的openFoam动网格算VIV的算例，从目前的结果来看，程序能正常运行，但程序运行结束以后，圆柱不动，想请各位大佬帮我看看，哪儿出问题了，谢谢大家。

@东岳 谢谢:baobao:

:141: :142: :143: :144: :145: :baobao:

额，自问自答？

1.我猜是symmetry这个边界条件，试一下改成wall，然后在0文件夹里自己定义上下的slip边界条件。openfoam会对生成的网格编号进行检查，外部软件的网格节点编号和blockMesh的逻辑不一样。
2.可能是前后的empty边界问题，你可以在starccm里生成3D网格后，重新在openfoam里extrude来生成前后empty边界
类似这样
constructFrom patch;
sourceCase "../w3-d_hc1-3";
sourcePatches (symFront);

// If construct from patch: patch to use for back (can be same as sourcePatch)
exposedPatchName symBack;

// Flip surface normals before usage. Valid only for extrude from surface or
// patch.
flipNormals false;

//- Linear extrusion in point-normal direction
extrudeModel linearNormal;

nLayers 1;

expansionRatio 1.0;

linearNormalCoeffs
{
thickness 1;
}

// Do front and back need to be merged? Usually only makes sense for 360
// degree wedges.
mergeFaces false; //true;

// Merge small edges. Fraction of bounding box.
mergeTol 0;

各位大佬好，我是CFDEM新手，最近尝试参照twoSpheresGlowinskiMPI构建自己的算例，采用dynamicRefine，在计算过程中断构建的算例基本与参照算例一致，但试跑算例出现中断，不知道是否是哪里没有注意到，求大佬指导
报错描述
当不使用自适应网格时，算例正常计算；当采用自适应网格功能interDyMFoam，网格未加密前（未达到设置加密条件），正常计算;网格开始加密计算中断。
报错代码如下：

Time = 0.0001 Selected 0 cells for refinement out of 1613716. Selected 0 split points out of a possible 0. Courant Number mean: 0.000199582 max: 0.0804916 - evolve() timeStepFraction() = 2 Starting up LIGGGHTS Executing command: 'run 10 ' Setting up run at Fri Jul 26 12:59:28 2024 Memory usage per processor = 6.75634 Mbytes Step Atoms KinEng rke Volume dragtota dragtota dragtota 1 1 2.5261073e-08 0 0.0260508 0 0 0 CFD Coupling established at step 10 11 1 3.5918126e-08 0 0.0260508 0 0 0 Loop time of 0.000962037 on 16 procs for 10 steps with 1 atoms, finish time Fri Jul 26 12:59:28 2024 Pair time (%) = 4.70437e-06 (0.489001) Neigh time (%) = 0 (0) Comm time (%) = 4.65044e-06 (0.483395) Outpt time (%) = 6.52643e-05 (6.78397) Other time (%) = 0.000887418 (92.2436) Nlocal: 0.0625 ave 1 max 0 min Histogram: 15 0 0 0 0 0 0 0 0 1 Nghost: 0 ave 0 max 0 min Histogram: 16 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 16 0 0 0 0 0 0 0 0 0 Total # of neighbors = 0 Ave neighs/atom = 0 Neighbor list builds = 0 Dangerous builds = 0 LIGGGHTS finished Foam::cfdemCloudIB::reAllocArrays() nr particles = 1 evolve done. DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 1.06984e-07, No Iterations 1 DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 4.99282e-07, No Iterations 1 DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 6.95219e-10, No Iterations 2 DICPCG: Solving for p, Initial residual = 1, Final residual = 9.26695e-07, No Iterations 573 time step continuity errors : sum local = 3.69904e-10, global = -1.28591e-13, cumulative = -1.28591e-13 DICPCG: Solving for p, Initial residual = 0.0476412, Final residual = 9.26405e-07, No Iterations 513 time step continuity errors : sum local = 3.16922e-08, global = -3.49546e-12, cumulative = -3.62406e-12 DICPCG: Solving for p, Initial residual = 0.0055667, Final residual = 9.51869e-07, No Iterations 500 time step continuity errors : sum local = 3.3706e-08, global = -6.10801e-12, cumulative = -9.73206e-12 DICPCG: Solving for p, Initial residual = 0.00178549, Final residual = 9.94251e-07, No Iterations 473 time step continuity errors : sum local = 3.53268e-08, global = -5.71629e-12, cumulative = -1.54484e-11 No finite volume options present particleCloud.calcVelocityCorrection() DICPCG: Solving for phiIB, Initial residual = 1, Final residual = 9.63924e-07, No Iterations 501 ExecutionTime = 30.7 s ClockTime = 32 s Time = 0.0002 Selected 784 cells for refinement out of 1613716. Refined from 1613716 to 1619204 cells. Selected 0 split points out of a possible 784. Courant Number mean: 0.126921 max: 9.46381 - evolve() timeStepFraction() = 2 Starting up LIGGGHTS Executing command: 'run 10 ' Setting up run at Fri Jul 26 12:59:56 2024 Memory usage per processor = 6.75634 Mbytes Step Atoms KinEng rke Volume dragtota dragtota dragtota 11 1 3.5918126e-08 0 0.0260508 0 0 0.017347056 CFD Coupling established at step 20 21 1 3.6972707e-08 0 0.0260508 0 0 0.017347056 Loop time of 0.00149142 on 16 procs for 10 steps with 1 atoms, finish time Fri Jul 26 12:59:56 2024 Pair time (%) = 3.75399e-06 (0.251706) Neigh time (%) = 0 (0) Comm time (%) = 3.65471e-06 (0.24505) Outpt time (%) = 8.7527e-05 (5.86871) Other time (%) = 0.00139648 (93.6345) Nlocal: 0.0625 ave 1 max 0 min Histogram: 15 0 0 0 0 0 0 0 0 1 Nghost: 0 ave 0 max 0 min Histogram: 16 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 16 0 0 0 0 0 0 0 0 0 Total # of neighbors = 0 Ave neighs/atom = 0 Neighbor list builds = 0 Dangerous builds = 0 LIGGGHTS finished nr particles = 1 -------------------------------------------------------------------------- MPI_ABORT was invoked on rank 15 in communicator MPI_COMM_WORLD with errorcode 1. NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes. You may or may not see output from other processes, depending on exactly when Open MPI kills them. --------------------------------------------------------------------------

导出时候选择了如下这些选项
image.png

@李东岳 好的李老师，感谢感谢

您好，我在使用comsol网格自适应算两相流会出现网格自适应速度跟不上两相运动速度，导致，这方面您是如何解决的呢