|
[Sponsors] |
[snappyHexMesh] defining layers in snappyhexmesh |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
February 9, 2011, 06:55 |
|
#21 |
New Member
Leon Eilders
Join Date: Jan 2011
Location: Enschede, The Netherlands
Posts: 13
Rep Power: 15 |
in reality there would be an outlet. that is what i want to simulate.
|
|
February 15, 2011, 09:44 |
|
#22 |
New Member
Leon Eilders
Join Date: Jan 2011
Location: Enschede, The Netherlands
Posts: 13
Rep Power: 15 |
Here we are again... Adapted my snappyhexmesh untill it passed the checkmesh. I used a higher refinement and set the maxskewness to 0.5. Giving quite good results.
Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create polyMesh for time = 0 Time = 0 Mesh stats points: 677206 faces: 1798140 internal faces: 1710919 cells: 561545 boundary patches: 8 point zones: 0 face zones: 0 cell zones: 0 Overall number of cells of each type: hexahedra: 460753 prisms: 30962 wedges: 0 pyramids: 0 tet wedges: 124 tetrahedra: 0 polyhedra: 69706 Checking topology... Boundary definition OK. Point usage OK. Upper triangular ordering OK. Face vertices OK. Number of regions: 1 (OK). Checking patch topology for multiply connected surfaces ... Patch Faces Points Surface topology maxY 210 249 ok (non-closed singly connected) minX 210 249 ok (non-closed singly connected) maxX 210 249 ok (non-closed singly connected) minY 210 249 ok (non-closed singly connected) minZ 799 974 ok (non-closed singly connected) maxZ 797 972 ok (non-closed singly connected) stlSurface_Sail 81905 85179 ok (non-closed singly connected) stlSurface_Suck 2880 3844 ok (non-closed singly connected) Checking geometry... Overall domain bounding box (-10.75 -10.3595 -5.52462e-06) (10 10 30) Mesh (non-empty, non-wedge) directions (1 1 1) Mesh (non-empty) directions (1 1 1) Boundary openness (2.53646e-16 1.24211e-16 -1.78792e-16) OK. Max cell openness = 9.50189e-16 OK. Max aspect ratio = 13.2989 OK. Minumum face area = 6.84828e-05. Maximum face area = 16.8984. Face area magnitudes OK. Min volume = 2.70464e-06. Max volume = 33.7968. Total volume = 12544.9. Cell volumes OK. Mesh non-orthogonality Max: 62.8324 average: 8.96609 Non-orthogonality check OK. Face pyramids OK. Max skewness = 1.2706 OK. Mesh OK. End Code:
Time = 70 smoothSolver: Solving for Ux, Initial residual = 0.00795705, Final residual = 7.13943e-05, No Iterations 4 smoothSolver: Solving for Uy, Initial residual = 0.0110534, Final residual = 0.000101015, No Iterations 4 smoothSolver: Solving for Uz, Initial residual = 0.101537, Final residual = 0.000387196, No Iterations 5 GAMG: Solving for p, Initial residual = 0.0691879, Final residual = 0.000240485, No Iterations 4 time step continuity errors : sum local = 0.000322106, global = -2.5441e-06, cumulative = -0.000247151 smoothSolver: Solving for epsilon, Initial residual = 0.384302, Final residual = 0.000275313, No Iterations 2 bounding epsilon, min: -577.335 max: 4.34232e+07 average: 272.56 smoothSolver: Solving for k, Initial residual = 0.0316469, Final residual = 0.000295912, No Iterations 3 bounding k, min: -1.23476 max: 3342.86 average: 0.483003 ExecutionTime = 260.45 s ClockTime = 261 s forceCoeffs output: Cd = 1.04372 Cl = 1.91814 Cm = -0.564935 forceCoeffs output: Cd = 0.966914 Cl = 1.87556 Cm = -0.564033 forceCoeffs output: Cd = 0.0768036 Cl = 0.0425825 Cm = -0.000902006 Time = 71 smoothSolver: Solving for Ux, Initial residual = 0.308738, Final residual = 0.0023828, No Iterations 5 smoothSolver: Solving for Uy, Initial residual = 0.167585, Final residual = 0.000735947, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.23801, Final residual = 0.00192168, No Iterations 4 GAMG: Solving for p, Initial residual = 0.0718548, Final residual = 0.000221095, No Iterations 4 time step continuity errors : sum local = 0.000302917, global = 1.43559e-05, cumulative = -0.000232795 smoothSolver: Solving for epsilon, Initial residual = 0.108878, Final residual = 0.00106995, No Iterations 4 bounding epsilon, min: -27638.7 max: 2.86248e+07 average: 249.215 smoothSolver: Solving for k, Initial residual = 0.271598, Final residual = 0.00251725, No Iterations 5 bounding k, min: -67.7914 max: 4138.36 average: 0.486237 ExecutionTime = 264.09 s ClockTime = 265 s forceCoeffs output: Cd = 1.24648 Cl = 1.95098 Cm = -0.468876 forceCoeffs output: Cd = 1.16514 Cl = 1.90592 Cm = -0.467914 forceCoeffs output: Cd = 0.0813401 Cl = 0.0450539 Cm = -0.00096219 Time = 72 smoothSolver: Solving for Ux, Initial residual = 0.999998, Final residual = 0.00776452, No Iterations 5 smoothSolver: Solving for Uy, Initial residual = 0.999999, Final residual = 0.00811549, No Iterations 4 smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 0.00648876, No Iterations 4 GAMG: Solving for p, Initial residual = 0.674664, Final residual = 0.00246658, No Iterations 4 time step continuity errors : sum local = 74714.8, global = 2712.55, cumulative = 2712.55 smoothSolver: Solving for epsilon, Initial residual = 1, Final residual = 0.00265908, No Iterations 5 bounding epsilon, min: -1.75341e+13 max: 3.03906e+14 average: 8.17846e+09 smoothSolver: Solving for k, Initial residual = 1, Final residual = 0.00157535, No Iterations 2 bounding k, min: -2.28545e+11 max: 9.35617e+12 average: 1.09552e+08 ExecutionTime = 267.75 s ClockTime = 269 s forceCoeffs output: Cd = -1.93271e+07 Cl = -3.8579e+07 Cm = -6.25219e+06 forceCoeffs output: Cd = -1.93341e+07 Cl = -3.85828e+07 Cm = -6.25209e+06 forceCoeffs output: Cd = 6971.54 Cl = 3785.83 Cm = -98.9127 Time = 73 smoothSolver: Solving for Ux, Initial residual = 0.765646, Final residual = 0.00289439, No Iterations 4 smoothSolver: Solving for Uy, Initial residual = 0.797186, Final residual = 0.00462204, No Iterations 3 smoothSolver: Solving for Uz, Initial residual = 0.656947, Final residual = 0.00441734, No Iterations 4 GAMG: Solving for p, Initial residual = 0.999727, Final residual = 0.00550768, No Iterations 5 time step continuity errors : sum local = 41235.4, global = 17613.6, cumulative = 20326.2 #0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam171/lib/linux64GccDPOpt/libOpenFOAM.so" #1 Foam::sigFpe::sigFpeHandler(int) in "/opt/openfoam171/lib/linux64GccDPOpt/libOpenFOAM.so" #2 in "/lib/libc.so.6" #3 Foam::LimitedScheme<double, Foam::limitedLinearLimiter<Foam::NVDTVD>, Foam::limitFuncs::magSqr>::limiter(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/opt/openfoam171/lib/linux64GccDPOpt/libfiniteVolume.so" #4 Foam::limitedSurfaceInterpolationScheme<double>::weights(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/opt/openfoam171/lib/linux64GccDPOpt/libfiniteVolume.so" #5 Foam::fv::gaussConvectionScheme<double>::fvmDiv(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&) const in "/opt/openfoam171/lib/linux64GccDPOpt/libfiniteVolume.so" #6 Foam::tmp<Foam::fvMatrix<double> > Foam::fvm::div<double>(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::word const&) in "/opt/openfoam171/lib/linux64GccDPOpt/libincompressibleRASModels.so" #7 Foam::tmp<Foam::fvMatrix<double> > Foam::fvm::div<double>(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&) in "/opt/openfoam171/lib/linux64GccDPOpt/libincompressibleRASModels.so" #8 Foam::incompressible::RASModels::RNGkEpsilon::correct() in "/opt/openfoam171/lib/linux64GccDPOpt/libincompressibleRASModels.so" #9 in "/opt/openfoam171/applications/bin/linux64GccDPOpt/simpleFoam" #10 __libc_start_main in "/lib/libc.so.6" #11 in "/opt/openfoam171/applications/bin/linux64GccDPOpt/simpleFoam" created temporary '35degr2.OpenFOAM' Code:
FoamFile { version 2.0; format ascii; class dictionary; object fvSolution; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // solvers { p { solver GAMG; tolerance 1e-8; relTol 0.01; //smoother GaussSeidel; smoother DIC; nPreSweeps 0; nPostSweeps 1; cacheAgglomeration on; agglomerator faceAreaPair; nCellsInCoarsestLevel 10; mergeLevels 1; }; U { solver smoothSolver; smoother GaussSeidel; tolerance 1e-7; relTol 0.01; nSweeps 1; }; k { solver smoothSolver; smoother GaussSeidel; tolerance 1e-6; relTol 0.01; nSweeps 1; }; epsilon { solver smoothSolver; smoother GaussSeidel; tolerance 1e-6; relTol 0.01; nSweeps 1; }; } SIMPLE { nNonOrthogonalCorrectors 0; //0; // 0, 1 or 2 pRefCell 0; pRefValue 0; } relaxationFactors { p 0.3; //0.3; U 0.7; //0.7; k 0.7; //0.7; epsilon 0.7; //0.7; } Code:
FoamFile { version 2.0; format ascii; class dictionary; object fvSchemes; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // ddtSchemes { default steadyState; } gradSchemes { default Gauss linear; //was 0.5 } divSchemes { default none; div(phi,U) Gauss GammaV 0.5; //was 0.5 1.0 div(phi,k) Gauss limitedLinear 1; div(phi,epsilon) Gauss limitedLinear 1; div((nuEff*dev(grad(U).T()))) Gauss linear; //div((nuEff*dev(grad(U).T()))) Gauss linear faceMDLimited Gauss linear 0.85; //div(phi,U) Gauss linearUpwindV faceMDLimited linear 0.68; //was 0.5 1.0 //div(phi,k) Gauss linearUpwind faceMDLimited Gauss linear 0.2; //div(phi,epsilon) Gauss linearUpwind faceMDLimited Gauss linear 0.91; //div((nuEff*dev(grad(U).T()))) Gauss linearUpwind faceMDLimited Gauss linear 0.85; } laplacianSchemes { default Gauss linear corrected; } interpolationSchemes { default linear; } snGradSchemes { default corrected; } fluxRequired { default no; p; } Do you have any suggestions ? Again many thanks for helping me so far! really appreciate it! Leon |
|
March 10, 2011, 05:15 |
|
#23 |
New Member
Leon Eilders
Join Date: Jan 2011
Location: Enschede, The Netherlands
Posts: 13
Rep Power: 15 |
nNonOrthogonalCorrectors 0;
had to be nNonOrthogonalCorrectors 2; Did the trick! Thanks for your help! |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[snappyHexMesh] Error defining boundary layer around cube snappyHexMesh | crizpi21 | OpenFOAM Meshing & Mesh Conversion | 5 | October 16, 2021 10:56 |
[snappyHexMesh] Triangular elements in boundary layers generated by snappyHexMesh | mike.franky | OpenFOAM Meshing & Mesh Conversion | 0 | June 3, 2018 09:46 |
[snappyHexMesh] SnappyHexMesh no layers and no decent mesh for complex geometry | pizzaspinate | OpenFOAM Meshing & Mesh Conversion | 1 | February 25, 2015 07:05 |
[snappyHexMesh] 2D snappyHexMesh with layers | grjmell | OpenFOAM Meshing & Mesh Conversion | 2 | May 27, 2013 12:00 |
[snappyHexMesh] SnappyHexMesh doesn't succeed in adding several layers | eti0123456789 | OpenFOAM Meshing & Mesh Conversion | 5 | April 9, 2013 02:36 |