|
[Sponsors] |
fvSolution fvSchemes problem for 3D circular cylinder |
![]() |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
![]() |
![]() |
#1 |
New Member
asalar.14
Join Date: Feb 2024
Posts: 6
Rep Power: 2 ![]() |
Hi everyone
I am not able to find the best fvSchemes and fvSolution for my case. need your advise for best setting The case is a 3D channel (Length= 1.52, Height= 0.5, width- 0.56) running with interFoam RANS kOmegaSST model having a circular cylinder in vicinity of the upstream boundary. Mesh non-orthogonality Max: 63.8068 average: 5.53946 I am checking momentum changes in two zones ( 1) the wake zone =4cm after the cylinder and 2) after the wake zone= 12cm after the cylinder) my optimum mesh size is (0.013 × 0.01 × 0.008) my deltaT is set to be 0.01 with maxCo=0.8 My fvSchemes is: /*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: v2212 | | \\ / A nd | Website: www.openfoam.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; object fvSchemes; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // ddtSchemes { default Euler; } gradSchemes { default cellMDLimited Gauss linear 0.5; grad(U) cellMDLimited Gauss linear 0.5; grad(p) Gauss linear; } divSchemes { div(phi,U) Gauss linearUpwind grad(U); div(phi,alpha) Gauss vanLeer; div(rhoPhi,U) Gauss upwind; div(phi,omega) Gauss linearUpwind default; div(phi,k) Gauss linearUpwind default; div(phirb,alpha) Gauss linear; div((nuEff*dev(T(grad(U))))) Gauss linear; div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear; } laplacianSchemes { default Gauss linear limited 0.777;; } interpolationSchemes { default linear; } snGradSchemes { default limited 0.777; } wallDist { method meshWave; } // ************************************************** *********************** // and my fvSolution is: /*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: v2212 | | \\ / A nd | Website: www.openfoam.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; object fvSolution; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // solvers { "alpha.water.*" { nAlphaCorr 1; nAlphaSubCycles 3; cAlpha 1; MULESCorr yes; nLimiterIter 3; solver smoothSolver; smoother symGaussSeidel; tolerance 1e-8; relTol 0; } p_rgh { solver GAMG; tolerance 1e-04; relTol 0.01; smoother GaussSeidel; nPreSweeps 0; nPostSweeps 2; nFinestSweeps 2; cacheAgglomeration true; nCellsInCoarsestLevel 1000; agglomerator faceAreaPair; mergeLevels 1; } p_rghFinal { $p_rgh; tolerance 1e-06; relTol 0; } "pcorr.*" { $p_rghFinal; tolerance 0.001; } "U.*" { solver PBiCGStab; preconditioner DILU; tolerance 1e-08; relTol 0; nSweeps 1; minIter 3; } "(k|omega|nuTilda).*" { solver smoothSolver; smoother symGaussSeidel; tolerance 1e-08; relTol 0; } } PIMPLE { momentumPredictor yes; nCorrectors 3; nOuterCorrectors 1; nNonOrthogonalCorrectors 1; } relaxationFactors { fields { "p." 0.7; } equations { "p.*" 0.7; "U.*" 0.7; "k.*" 0.7; "omega.*" 0.7; ".*" 1; } } // ************************************************** *********************** // Now what needs to be changed or omitted for better accuracy, faster convergence and stability? |
|
![]() |
![]() |
![]() |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
How to implement time-varying boundary conditions in an enclosed cylinder | Alejandro-FA | OpenFOAM Pre-Processing | 0 | May 14, 2020 06:29 |
SU2-7.0.1 on ubuntu 18.04 | hyunko | SU2 Installation | 7 | March 16, 2020 04:37 |
Drag force coefficient too high for a flow past a cylinder using komega sst | Scabbard | OpenFOAM Running, Solving & CFD | 37 | March 21, 2016 16:16 |
Cylinder motion along a circular path | karamiag | OpenFOAM Running, Solving & CFD | 0 | May 25, 2014 05:05 |
flow around a circular cylinder with velocity inlet and outflow outlet | shuoxue | OpenFOAM Running, Solving & CFD | 0 | November 2, 2013 04:32 |