|
[Sponsors] |
May 7, 2018, 05:03 |
non physical value of alpha
|
#1 | |
Member
L S
Join Date: Apr 2016
Posts: 63
Rep Power: 9 |
Hi,
I am trying to simulate condensation problem using phaseChangeHeatFoam using Dynamic Mesh. But, I am getting following error. Where Value of alpha is going -751. WHich is non physical. Quote:
Code: Code:
ddtSchemes { default Euler; } gradSchemes { default Gauss linear; } divSchemes { div(rho*phi,U) Gauss upwind; div(rhoPhi,U) Gauss vanLeerV; div(phi,T) Gauss vanLeer; div(rhoPhi,T) Gauss vanLeer; div(phi,alpha) Gauss vanLeer01; div(phirb,alpha) Gauss interfaceCompression; div((muEff*dev(T(grad(U))))) Gauss linear; } laplacianSchemes { default Gauss linear corrected; } interpolationSchemes { default linear; } snGradSchemes { default corrected; } fluxRequired { default no; p_rgh; pcorr; alpha; } Code:
solvers { alpha1 { maxUnboundedness 1e-5; CoCoeff 2; maxIter 5; nLimiterIter 2; solver PBiCG; preconditioner DILU; tolerance 1e-12; relTol 0.01; }; "(U|UFinal|T|TFinal)" { solver PBiCG; preconditioner DILU; tolerance 1e-06; relTol 0; }; "(p_rgh|p)" { solver PCG; preconditioner { preconditioner GAMG; tolerance 1e-08; relTol 0; nVcycles 2; smoother GaussSeidel; nPreSweeps 0; nPostSweeps 2; nFinestSweeps 2; cacheAgglomeration true; nCellsInCoarsestLevel 10; agglomerator faceAreaPair; mergeLevels 1; } tolerance 1e-08; relTol 0.001; maxIter 20; }; pcorr { solver PCG; preconditioner { preconditioner GAMG; tolerance 0.001; relTol 0; smoother GaussSeidel; nPreSweeps 0; nPostSweeps 2; nFinestSweeps 2; cacheAgglomeration false; nCellsInCoarsestLevel 10; agglomerator faceAreaPair; mergeLevels 1; } tolerance 0.0001; relTol 0; maxIter 100; }; "(p_rghFinal|pFinal)" { solver PCG; preconditioner { preconditioner GAMG; tolerance 1e-08; relTol 0; nVcycles 2; smoother GaussSeidel; nPreSweeps 2; nPostSweeps 0; nFinestSweeps 2; cacheAgglomeration true; nCellsInCoarsestLevel 10; agglomerator faceAreaPair; mergeLevels 1; } tolerance 1e-08; relTol 0; }; p_rgh { solver GAMG; tolerance 1e-08; relTol 0.05; smoother GaussSeidel; nPreSweeps 0; nPostSweeps 2; nFinestSweeps 2; cacheAgglomeration true; nCellsInCoarsestLevel 10; agglomerator faceAreaPair; mergeLevels 1; } p_rghFinal { solver PCG; preconditioner { preconditioner GAMG; tolerance 1e-08; relTol 0; nVcycles 2; smoother GaussSeidel; nPreSweeps 0; nPostSweeps 2; nFinestSweeps 2; cacheAgglomeration true; nCellsInCoarsestLevel 10; agglomerator faceAreaPair; mergeLevels 1; } tolerance 1e-08; relTol 0; maxIter 20; } U { solver smoothSolver; smoother GaussSeidel; tolerance 1e-06; relTol 0; nSweeps 1; } "(k|B|nuTilda)" { solver PBiCG; preconditioner DILU; tolerance 1e-08; relTol 0; } /* T //add { solver PBiCG; preconditioner DILU; tolerance 1e-10; relTol 0; }; /*TFinal //added { solver PBiCG; preconditioner DILU; tolerance 1e-10; relTol 0; }; */ "(T|TFinal)" { solver PBiCG; preconditioner DILU; tolerance 1e-10; relTol 0; }; rho //add { solver Diagonal; tolerance 1e-8; relTol 0; }; } relaxationFactors { T 0.3 U 0.3 p 0.3; } "(PISO|PIMPLE)" { momentumPredictor yes; nCorrectors 3; nNonOrthogonalCorrectors 0; nAlphaCorr 1; nAlphaSubCycles 3; cAlpha 1; pRefPoint (0.51 0.51 0.51); pRefValue 0; smoothItr 2; kSmoothItr 0; } |
||
May 18, 2018, 09:43 |
|
#2 |
Member
L S
Join Date: Apr 2016
Posts: 63
Rep Power: 9 |
I got the solution. I just reduced the value of deltaT from 0.001 to 0.000001.
|
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[Gmsh] Nasa Nozzle: Plot 3D -->gmsh -->fluent - no physical groups | Nico89 | OpenFOAM Meshing & Mesh Conversion | 1 | June 14, 2016 06:40 |
alphaEqn.H in twoPhaseEulerFoam | cheng1988sjtu | OpenFOAM Bugs | 15 | May 1, 2016 17:12 |
Superlinear speedup in OpenFOAM 13 | msrinath80 | OpenFOAM Running, Solving & CFD | 18 | March 3, 2015 06:36 |
Partition does not start on physical sector boundary!!?? | immortality | Lounge | 1 | December 8, 2013 16:04 |
Determining alpha and beta for porous baffle | Liaqat Khan | Siemens | 1 | October 27, 2000 05:44 |