请教各位大佬关于shm画网格的问题
-
各位大佬好!我想请教一下为什么下面这个suboff在尖锐处和尾翼会丢失边界层,另外想请教一下可以优化控制文件的哪些部分来进行改善,感觉snap这一步尾翼附近的网格就变得不太好了,感谢各位大佬,祝各位科研顺利
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.3.0 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; location system; object snappyHexMeshDict; } castellatedMesh true; snap true; addLayers true; geometry { model.stl { type triSurfaceMesh; name model; patchInfo { type wall; } } sail { type searchableBox; min (0.8 -0.1 0.19); max (1.4 0.1 0.47); } wake { type searchableBox; min (-0.5 -0.5 -0.5); max ( 6.0 0.5 0.5); } propeller { type searchableBox; min (3.6 -0.3 -0.3); max (4.4 0.3 0.3); } } castellatedMeshControls { features ( { file "model.eMesh"; level 1; } ); refinementSurfaces { model { level (3 4); } } refinementRegions { sail { mode inside; levels ((3 4)); } wake { mode inside; levels ((2 2)); } propeller { mode inside; levels ((3 4)); } } locationInMesh ( 6.0 1.0 1.0); maxLocalCells 10000000; maxGlobalCells 50000000; minRefinementCells 0; nCellsBetweenLevels 3; resolveFeatureAngle 30; allowFreeStandingZoneFaces false; planarAngle 30; maxLoadUnbalance 0.10; } snapControls { nSolveIter 300; nSmoothPatch 10; tolerance 1.5; nRelaxIter 10; nFeatureSnapIter 10; implicitFeatureSnap false; explicitFeatureSnap true; multiRegionFeatureSnap false; } addLayersControls { layers { "model" {nSurfaceLayers 8;} } relativeSizes true; expansionRatio 1.3; finalLayerThickness 0.4; minThickness 0.0001; nGrow 0; featureAngle 180; slipFeatureAngle 80; nRelaxIter 5; nSmoothSurfaceNormals 1; nSmoothNormals 3; nSmoothThickness 10; maxFaceThicknessRatio 0.9; maxThicknessToMedialRatio 0.3; minMedialAxisAngle 90; nBufferCellsNoExtrude 0; nLayerIter 50; nRelaxedIter 20; writeVTK false; noErrors false; layerRecovery 1; growZoneLayers false; projectGrownUp 0.0; } meshQualityControls { maxNonOrtho 65; maxBoundarySkewness 4; maxInternalSkewness 4; maxConcave 180; minFlatness 0.5; minVol 1.00E-13; minArea -1; minTetQuality 1e-30; minTwist 0.01; minDeterminant 0.01; minFaceWeight 0.01; minVolRatio 0.01; minTriangleTwist -1; nSmoothScale 4; errorReduction 0.75; } debug 0; mergeTolerance 1E-6; autoBlockMesh false;
-
如果你们看openfoam自带的,也这样。不过自带的应该是只是加了一个面。SHM加边界层都是老大难了。这部分工作应该没有人赞助。Mattijs也不在基金会这么。前几年ESI那面snappyHexMesh改动的挺多,不过这也是老大难。
商软跟开源区别挺大。商软卖个几千万,团队人员咔咔咔有活就干不闲着,招过来本身就是做网格的,那就得一直开发网格。开源这面核心没几个人,不存在固定做某个分支的人,有赞助就全干这个。没赞助就靠别人贡献。有贡献就有bug提升。没贡献就没有。
OpenFOAM那个AMI之前也一直没有赞助,也没贡献,也属于老大难。后来Will过去之后,把这个给弄了。