|
[Sponsors] |
PisoFoam and PimpleFoam for my problem can not converge |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
August 1, 2013, 06:56 |
PisoFoam and PimpleFoam for my problem can not converge
|
#1 |
Senior Member
Join Date: Jun 2011
Posts: 163
Rep Power: 15 |
I want to solve flow around a cylinder and thin plate
the thin plate is placed horizontally after cylinder I have test many different settings in fvSchemes and fvSolutions but all of them give uncorrected results I will be so grateful if anybody can help me Code:
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.1.1 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; object fvSchemes; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // ddtSchemes { default Euler; } gradSchemes { default leastSquares ; // grad(p) leastSquares ; // grad(U) leastSquares ; } divSchemes { default none; div(phi,U) Gauss limitedLinearV 1;//Gauss linearUpwind grad(U);// div(phi,k) Gauss limitedLinear 1; div(phi,omega) Gauss limitedLinear 1; div((nuEff*dev(T(grad(U))))) Gauss linear; div((nuEff*dev(grad(U).T()))) Gauss linear; } laplacianSchemes { default Gauss linear corrected;//Gauss linear limited 0.5;// } interpolationSchemes { default linear; } snGradSchemes { default corrected; } fluxRequired { default no; // pcorr ; p; } // ************************************************************************* // Code:
/*---------------------------------------------------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 1.3 | | \\ / A nd | Web: http://www.openfoam.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; object fvSolution; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // solvers { p { solver GAMG; tolerance 1e-5; relTol 0; smoother GaussSeidel; nPreSweeps 1; nPostSweeps 2; cacheAgglomeration true; directSolveCoarsest true; agglomerator faceAreaPair; nCellsInCoarsestLevel 400; mergeLevels 1; minIter 1; } pFinal { solver GAMG; tolerance 1e-9; relTol 0; smoother GaussSeidel; nPreSweeps 1; nPostSweeps 2; cacheAgglomeration true; agglomerator faceAreaPair; nCellsInCoarsestLevel 400; mergeLevels 1; minIter 1; } "(U|k|omega)" { // solver PBiCG; // preconditioner DILU; solver smoothSolver; smoother GaussSeidel; tolerance 1e-5; relTol 0; nSweeps 1; minIter 1; } "(UFinal|kFinal|omegaFinal)" { // solver PBiCG; // preconditioner DILU; solver smoothSolver; smoother GaussSeidel; tolerance 1e-6; relTol 0.0; nSweeps 1; minIter 1; } cellMotionU { solver GAMG; tolerance 1e-6; relTol 1e-3; smoother GaussSeidel; cacheAgglomeration true; directSolveCoarsest true; nCellsInCoarsestLevel 20; agglomerator faceAreaPair; mergeLevels 1; } } PIMPLE { nOuterCorrectors 3; nCorrectors 2; nNonOrthogonalCorrectors 3; pRefCell 0; pRefValue 0; correctPhi no; // momentumPredictor yes; } PISO { nOuterCorrectors 2; nCorrectors 2; nNonOrthogonalCorrectors 3; pRefCell 0; pRefValue 0; // ddtPhiCorr no; } relaxationFactors { fields { p 0.3; } equations { U 0.7; k 0.7; omega 0.7; } } potentialFlow { nNonOrthogonalCorrectors 0; } cache { grad(U); } // ************************************************************************* // |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
pimpleFoam vs simpleFoam vs pisoFoam vs icoFoam? | phsieh2005 | OpenFOAM Running, Solving & CFD | 45 | March 22, 2021 10:14 |
I got complete different results for pisofoam and pimplefoam? | lakeat | OpenFOAM Running, Solving & CFD | 6 | December 27, 2012 05:25 |