|
[Sponsors] |
![]() |
![]() |
#1 |
New Member
Aydar
Join Date: Dec 2021
Location: Russia, Kazan
Posts: 4
Rep Power: 4 ![]() |
Trying to calculate the velocity profile on the axis of the pipeline. Does not match analytical calculation.
long pipe section, y+ = 1, simpleFoam solver. Please tell me how to set the initial settings. Code:
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: v2106 | | \\ / A nd | Website: www.openfoam.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; object fvSchemes; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // ddtSchemes { default steadyState; } gradSchemes { default Gauss linear; } divSchemes { default none; div(phi,U) bounded Gauss upwind; div(phi,p) bounded Gauss upwind; turbulence bounded Gauss upwind; div(phi,k) $turbulence; div(phi,omega) $turbulence; div((nuEff*dev2(T(grad(U))))) Gauss linear; } laplacianSchemes { default Gauss linear corrected; } interpolationSchemes { default linear; } snGradSchemes { default corrected; } wallDist { method meshWave; } // ************************************************************************* // Code:
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: v2106 | | \\ / A nd | Website: www.openfoam.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; object fvSolution; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // solvers { p { solver GAMG; tolerance 1e-06; relTol 0.1; smoother GaussSeidel; nPreSweeps 0; nPostSweeps 2; cacheAgglomeration on; agglomerator faceAreaPair; nCellsInCoarsestLevel 10; mergeLevels 1; } "(U|k|omega|R|nuTilda)" { solver smoothSolver; smoother symGaussSeidel; tolerance 1e-05; relTol 0.1; } } SIMPLE { nNonOrthogonalCorrectors 0; } relaxationFactors { fields { p 0.3; } equations { U 0.7; k 0.7; omega 0.7; } } // ************************************************************************* //
__________________
Best regards,Mukhametov Aydar |
|
![]() |
![]() |
![]() |
![]() |
#2 |
Member
MNM
Join Date: Aug 2017
Posts: 69
Rep Power: 9 ![]() |
Hello Aydar,
To begin with, you should switch to more accurate div scheme, something like Code:
div(phi,U) bounded Gauss linearUpwind grad(U); Do a checkMesh and see the max non - orthogonality, because now its bit contrdictory, as in fvSolutions you dont have any nonOrthoCorrectors but in fvSchemes you are using "corrector" for snGrad. P.S. I did grep in tutorials and found that, you will never need div(phi, p), unless either you are solving for Compressible multiphase flows or have modified the solver to obtain something else. Last edited by SHUBHAM9595; December 8, 2021 at 15:07. |
|
![]() |
![]() |
![]() |
Tags |
openfoam, simplefoam, yplus |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
[UDF] Relative coordinates in logarithmic velocity profile | cfdworker | Fluent UDF and Scheme Programming | 25 | July 13, 2017 03:44 |
3d velocity profile at velocity inlet | swethaprakash | FLUENT | 3 | September 1, 2015 02:10 |
Logarithmic velocity profile | cfdworker | Fluent UDF and Scheme Programming | 0 | April 23, 2009 19:09 |
[boundary condition] logarithmic velocity profile | cfdworker | FLUENT | 2 | April 17, 2009 23:36 |
Prescribed inflow velocity profile - how to? | Alan | Main CFD Forum | 10 | October 28, 2005 12:14 |