进口马赫数0.7左右，用rhosimpleFoam不收敛的问题
-
介质就是空气
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 3.0.1 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; location "constant"; object thermophysicalProperties; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // thermoType { type hePsiThermo; mixture pureMixture; transport sutherland; thermo hConst; equationOfState perfectGas; specie specie; energy sensibleInternalEnergy; } mixture { specie { nMoles 1; molWeight 28.9; } thermodynamics { Cp 1007; Hf 0; } transport { As 1.4792e-06; Ts 116; } } // ************************************************************************* // -
| ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 3.0.1 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; location "constant"; object turbulenceProperties; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // simulationType RAS; RAS { // Tested with kEpsilon, realizableKE, kOmega, kOmegaSST, v2f, // ShihQuadraticKE, LienCubicKE. RASModel kOmegaSST; turbulence on; printCoeffs on; } // ************************************************************************* // -
网格checkmesh OK
| ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 3.0.1 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class polyBoundaryMesh; location "constant/polyMesh"; object boundary; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 8 ( hub { type wall; nFaces 16560; startFace 2399630; } inlet { type patch; nFaces 3430; startFace 2416190; } outlet { type patch; nFaces 2842; startFace 2419620; } per1 { type cyclicAMI; nFaces 9163; startFace 2422462; neighbourPatch per2; matchTolerance 0.01; transform translational; separationVector (0 0.022 0); } per2 { type cyclicAMI; nFaces 9163; startFace 2431625; neighbourPatch per1; matchTolerance 0.01; transform translational; separationVector (0 -0.022 0); } ps { type wall; nFaces 5831; startFace 2440788; } ss { type wall; nFaces 5831; startFace 2446619; } sym { type symmetryPlane; nFaces 16560; startFace 2452450; } ) // ************************************************************************* // -
| ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 3.0.1 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; object controlDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // application rhoSimpleFoam; startFrom latestTime; startTime 0; stopAt endTime; endTime 100; deltaT 1e-4 writeControl timeStep; writeInterval 1000; purgeWrite 0; writeFormat ascii; writePrecision 6; writeCompression off; timeFormat general; timePrecision 6; graphFormat raw; runTimeModifiable true; adjustTimeStep on; // ************************************************************************* // -
| ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 3.0.1 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ 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(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear; div(phi,e) bounded Gauss upwind; div(phi,epsilon) bounded Gauss upwind; div(phi,k) bounded Gauss upwind; div(phi,Ekp) bounded Gauss upwind; div(phi,K) bounded Gauss upwind; div(phi,h) bounded Gauss upwind; div(phi,omega) bounded Gauss upwind; div(phid,p) bounded Gauss upwind; } laplacianSchemes { default Gauss linear corrected; } interpolationSchemes { default linear; } snGradSchemes { default corrected; } wallDist { method meshWave; } fluxRequired { default no; p ; } // ************************************************************************* // -
| ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 3.0.1 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; object fvSolution; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // solvers { p { solver GAMG; tolerance 1e-05; relTol 0.1; smoother GaussSeidel; nPreSweeps 0; nPostSweeps 2; nFinestSweeps 2; cacheAgglomeration true; nCellsInCoarsestLevel 20; agglomerator faceAreaPair; mergeLevels 1; } "(U|omega|e|k|epsilon)" { solver smoothSolver; smoother GaussSeidel; tolerance 1e-8; relTol 0.1; nSweeps 1; } } SIMPLE { nNonOrthogonalCorrectors 1; rhoMin 0.1; rhoMax 2; // consistent yes; transonic yes; residualControl { p 1e-5; U 1e-5; e 1e-5; // possibly check turbulence fields "(k|epsilon|omega)" 1e-5; } } relaxationFactors { fields { p 0.3; rho 0.7; } equations { U 0.7; "(k|epsilon)" 0.7; e 0.7; } } // ************************************************************************* // -
omega飙升
Time = 0.0275 Duration of pickling 0.0523500442505 too long. Extending frequency from 1.0 to 2.61750221252 smoothSolver: Solving for Ux, Initial residual = 0.00152996, Final residual = 0.000123741, No Iterations 5 smoothSolver: Solving for Uy, Initial residual = 0.00222516, Final residual = 0.000169915, No Iterations 5 smoothSolver: Solving for Uz, Initial residual = 0.00285202, Final residual = 0.000244214, No Iterations 5 smoothSolver: Solving for e, Initial residual = 0.00798671, Final residual = 0.000494355, No Iterations 3 GAMG: Solving for p, Initial residual = 0.0160753, Final residual = 0.0013806, No Iterations 4 GAMG: Solving for p, Initial residual = 0.00143423, Final residual = 0.000143128, No Iterations 94 time step continuity errors : sum local = 0.32827, global = 0.0097095, cumulative = 1.35785 rho max/min : 1.34304 1.20262 smoothSolver: Solving for omega, Initial residual = 1, Final residual = 1.08444e-49, No Iterations 1 bounding omega, min: -5.71564e+59 max: 7.78487e+77 average: 1.01456e+72 smoothSolver: Solving for k, Initial residual = 8.8899e-07, Final residual = 1.0766e-40, No Iterations 1 bounding k, min: -8.90538e-08 max: 713.841 average: 2.87127 ExecutionTime = 3087.11 s ClockTime = 3100 s Time = 0.0276 Duration of pickling 0.0545029640198 too long. Extending frequency from 1.0 to 2.72514820099 smoothSolver: Solving for Ux, Initial residual = 0.0015129, Final residual = 0.000121723, No Iterations 5 smoothSolver: Solving for Uy, Initial residual = 0.00220558, Final residual = 0.000167395, No Iterations 5 smoothSolver: Solving for Uz, Initial residual = 0.00284466, Final residual = 0.000244908, No Iterations 5 smoothSolver: Solving for e, Initial residual = 0.00802377, Final residual = 0.000497091, No Iterations 3 GAMG: Solving for p, Initial residual = 0.0172986, Final residual = 0.00170545, No Iterations 3 GAMG: Solving for p, Initial residual = 0.00174044, Final residual = 0.000172569, No Iterations 68 time step continuity errors : sum local = 0.326715, global = 0.0096409, cumulative = 1.36749 rho max/min : 1.34304 1.20279 smoothSolver: Solving for omega, Initial residual = 1, Final residual = 2.9831e-72, No Iterations 1 bounding omega, min: -1.96379e+59 max: 1.91134e+73 average: 2.35575e+67 smoothSolver: Solving for k, Initial residual = 1.89765e-06, Final residual = 1.06009e-29, No Iterations 1 bounding k, min: -8.15996e-06 max: 718.357 average: 2.79046 ExecutionTime = 3096.76 s ClockTime = 3109 s Time = 0.0277 Duration of pickling 0.0812048912048 too long. Extending frequency from 1.0 to 4.06024456024 smoothSolver: Solving for Ux, Initial residual = 0.0014993, Final residual = 0.000119976, No Iterations 5 smoothSolver: Solving for Uy, Initial residual = 0.00218975, Final residual = 0.000165293, No Iterations 5 smoothSolver: Solving for Uz, Initial residual = 0.00283678, Final residual = 0.000245435, No Iterations 5 smoothSolver: Solving for e, Initial residual = 0.0080669, Final residual = 0.00049988, No Iterations 3 GAMG: Solving for p, Initial residual = 0.0190055, Final residual = 0.00162944, No Iterations 3 GAMG: Solving for p, Initial residual = 0.00165444, Final residual = 0.000165306, No Iterations 58 time step continuity errors : sum local = 0.325202, global = 0.0095477, cumulative = 1.37704 rho max/min : 1.34304 1.20296 smoothSolver: Solving for omega, Initial residual = 1, Final residual = 1.65172e-63, No Iterations 1 bounding omega, min: -3.89243e+67 max: 4.18007e+77 average: 6.1041e+71 smoothSolver: Solving for k, Initial residual = 4.19584e-08, Final residual = 1.05679e-23, No Iterations 1 bounding k, min: -2.76439e-08 max: 720.659 average: 2.71377 ExecutionTime = 3105.59 s ClockTime = 3118 s Time = 0.0278 Duration of pickling 0.0585401058197 too long. Extending frequency from 1.0 to 2.92700529099 smoothSolver: Solving for Ux, Initial residual = 0.00149028, Final residual = 0.000118581, No Iterations 5 smoothSolver: Solving for Uy, Initial residual = 0.00218168, Final residual = 0.000163814, No Iterations 5 smoothSolver: Solving for Uz, Initial residual = 0.00282467, Final residual = 0.000245514, No Iterations 5 smoothSolver: Solving for e, Initial residual = 0.00812588, Final residual = 0.000503407, No Iterations 3 GAMG: Solving for p, Initial residual = 0.0211776, Final residual = 0.00167735, No Iterations 3 GAMG: Solving for p, Initial residual = 0.00169416, Final residual = 0.000167178, No Iterations 52 time step continuity errors : sum local = 0.323742, global = 0.0094258, cumulative = 1.38646 rho max/min : 1.34304 1.18013 smoothSolver: Solving for omega, Initial residual = 1, Final residual = 5.93357e-72, No Iterations 1 bounding omega, min: -1.81139e+71 max: 4.20153e+75 average: 5.17821e+69 smoothSolver: Solving for k, Initial residual = 1.6527e-05, Final residual = 5.61923e-22, No Iterations 1 bounding k, min: -6.39732e-08 max: 711.177 average: 2.6393 ExecutionTime = 3113.97 s ClockTime = 3127 s Time = 0.0279 Duration of pickling 0.0739150047302 too long. Extending frequency from 1.0 to 3.69575023651 smoothSolver: Solving for Ux, Initial residual = 0.00149312, Final residual = 0.000117778, No Iterations 5 smoothSolver: Solving for Uy, Initial residual = 0.0021903, Final residual = 0.000163573, No Iterations 5 smoothSolver: Solving for Uz, Initial residual = 0.00283009, Final residual = 0.00024665, No Iterations 5 smoothSolver: Solving for e, Initial residual = 0.00820162, Final residual = 0.000508201, No Iterations 3 GAMG: Solving for p, Initial residual = 0.0246477, Final residual = 0.00198208, No Iterations 3 GAMG: Solving for p, Initial residual = 0.00200497, Final residual = 0.000198396, No Iterations 55 time step continuity errors : sum local = 0.322345, global = 0.00924313, cumulative = 1.39571 rho max/min : 1.34304 1.17683 smoothSolver: Solving for omega, Initial residual = 1, Final residual = 1.13229e-69, No Iterations 1 bounding omega, min: -4.25636e+69 max: 8.12744e+67 average: -5.13238e+63 smoothSolver: Solving for k, Initial residual = 2.60154e-07, Final residual = 7.34386e-26, No Iterations 1 bounding k, min: -5.22454e-09 max: 713.695 average: 2.56462 ExecutionTime = 3122.56 s ClockTime = 3135 s Time = 0.028 Duration of pickling 0.067948102951 too long. Extending frequency from 1.0 to 3.39740514755 smoothSolver: Solving for Ux, Initial residual = 0.00150484, Final residual = 0.000117816, No Iterations 5 smoothSolver: Solving for Uy, Initial residual = 0.00220887, Final residual = 0.000164596, No Iterations 5 smoothSolver: Solving for Uz, Initial residual = 0.00285433, Final residual = 0.000249085, No Iterations 5 smoothSolver: Solving for e, Initial residual = 0.00832458, Final residual = 0.00051524, No Iterations 3 GAMG: Solving for p, Initial residual = 0.0284494, Final residual = 0.00250348, No Iterations 4 GAMG: Solving for p, Initial residual = 0.00262877, Final residual = 0.000259279, No Iterations 92 time step continuity errors : sum local = 0.321011, global = 0.00897452, cumulative = 1.40468 rho max/min : 1.34305 1.09793 smoothSolver: Solving for omega, Initial residual = 1, Final residual = 3.60599e-62, No Iterations 1 bounding omega, min: -1.30817e+62 max: 5.18859e+60 average: -1.70377e+56 smoothSolver: Solving for k, Initial residual = 0.650092, Final residual = 5.24055e-24, No Iterations 1 bounding k, min: -5.18956e-11 max: 710.995 average: 2.49681 ExecutionTime = 3134.49 s ClockTime = 3147 s Time = 0.0281 Duration of pickling 0.0619690418243 too long. Extending frequency from 1.0 to 3.09845209122 smoothSolver: Solving for Ux, Initial residual = 0.00152832, Final residual = 0.000117872, No Iterations 5 smoothSolver: Solving for Uy, Initial residual = 0.00223761, Final residual = 0.000166456, No Iterations 5 smoothSolver: Solving for Uz, Initial residual = 0.00292184, Final residual = 0.000253548, No Iterations 5 smoothSolver: Solving for e, Initial residual = 0.00846883, Final residual = 0.000523429, No Iterations 3 GAMG: Solving for p, Initial residual = 0.0354473, Final residual = 0.00331063, No Iterations 6 GAMG: Solving for p, Initial residual = 0.00364502, Final residual = 0.000362752, No Iterations 125 time step continuity errors : sum local = 0.319631, global = 0.00861118, cumulative = 1.41329 rho max/min : 1.34305 1.06718 smoothSolver: Solving for omega, Initial residual = 1, Final residual = 9.78348e-48, No Iterations 1 bounding omega, min: -4.71803e+50 max: 1.27516e+67 average: 1.57147e+61 smoothSolver: Solving for k, Initial residual = 0.993077, Final residual = 6.53872e-31, No Iterations 1 bounding k, min: -1.47294e-10 max: 799.109 average: 2.43261 ExecutionTime = 3149.37 s ClockTime = 3162 s Time = 0.0282 Duration of pickling 0.0754199028015 too long. Extending frequency from 1.0 to 3.77099514008 smoothSolver: Solving for Ux, Initial residual = 0.00160754, Final residual = 0.000122563, No Iterations 5 smoothSolver: Solving for Uy, Initial residual = 0.00234167, Final residual = 0.000174383, No Iterations 5 smoothSolver: Solving for Uz, Initial residual = 0.00305302, Final residual = 0.000264702, No Iterations 5 smoothSolver: Solving for e, Initial residual = 0.00890675, Final residual = 0.000544697, No Iterations 3 GAMG: Solving for p, Initial residual = 0.0456857, Final residual = 0.00450887, No Iterations 11 GAMG: Solving for p, Initial residual = 0.00532448, Final residual = 0.000530223, No Iterations 148 time step continuity errors : sum local = 0.318517, global = 0.00804085, cumulative = 1.42133 rho max/min : 1.34304 0.944974 smoothSolver: Solving for omega, Initial residual = 1, Final residual = 5.03085e-61, No Iterations 1 bounding omega, min: -6.12926e+56 max: 1.89257e+63 average: 2.58001e+57 smoothSolver: Solving for k, Initial residual = 0.0986294, Final residual = 2.34926e-17, No Iterations 1 bounding k, min: -3.64697e-08 max: 1561.98 average: 2.37097 ExecutionTime = 3166.71 s ClockTime = 3180 s Time = 0.0283 Duration of pickling 0.0689640045166 too long. Extending frequency from 1.0 to 3.44820022583 smoothSolver: Solving for Ux, Initial residual = 0.00188211, Final residual = 0.00013498, No Iterations 5 smoothSolver: Solving for Uy, Initial residual = 0.0026032, Final residual = 0.000188799, No Iterations 5 smoothSolver: Solving for Uz, Initial residual = 0.0033437, Final residual = 0.000287587, No Iterations 5 smoothSolver: Solving for e, Initial residual = 0.00997634, Final residual = 0.000591898, No Iterations 3 GAMG: Solving for p, Initial residual = 0.0732098, Final residual = 0.00729724, No Iterations 26 GAMG: Solving for p, Initial residual = 0.0102857, Final residual = 0.00102317, No Iterations 160 time step continuity errors : sum local = 0.316687, global = 0.00719629, cumulative = 1.42853 rho max/min : 1.368 0.265827 smoothSolver: Solving for omega, Initial residual = 1, Final residual = 2.76757e-55, No Iterations 1 bounding omega, min: -2.79393e+51 max: 1.89919e+65 average: 2.72917e+59 smoothSolver: Solving for k, Initial residual = 0.00338243, Final residual = 8.55618e-21, No Iterations 1 bounding k, min: -2.59751e-09 max: 2425.35 average: 2.31181 ExecutionTime = 3186.64 s ClockTime = 3200 s Time = 0.0284 Duration of pickling 0.0688588619232 too long. Extending frequency from 1.0 to 3.44294309616 smoothSolver: Solving for Ux, Initial residual = 0.00844988, Final residual = 0.000751562, No Iterations 3 smoothSolver: Solving for Uy, Initial residual = 0.00963742, Final residual = 0.000907189, No Iterations 3 smoothSolver: Solving for Uz, Initial residual = 0.00423285, Final residual = 0.000343988, No Iterations 5 smoothSolver: Solving for e, Initial residual = 0.24105, Final residual = 0.013761, No Iterations 2 #0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 ? in "/lib/x86_64-linux-gnu/libc.so.6" #3 Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate() at ??:? #4 Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct() at ??:? #5 ? at ??:? #6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #7 ? at ??:? Floating point exception (core dumped) cc@ubuntu:~/OpenFOAM/cc-3.0.1/run/tutorials/compressible/rhoSimpleFoam$ -
要不算个初始场？
在fvSolution中加
solvers { ... Phi { solver PCG; preconditioner PIC; tolerance 1e-12; relTol 0; } } potentialFlow { nCorrectors 10; }然后运行potentialFlow搞个速度初始场试试？
或者在controlDict中加一些debugSwitches
DebugSwitches { lduMatrix 2; ...// 去$FOAM_ETC/controlDict下找和湍流有关的。 }还有就是加functionObject : fieldMinMax看看是哪个位置的场出了问题
fieldMinMax1 { type fieldMinMax; libs ("libfieldFunctionObjects.so"); write yes; log yes; location yes; // 关键是这个！ mode magnitude; fields ( U p omega ); } -
-
@RebelYoung http://dyfluid.com/rhoSimpleFoam.html 这面有一些rhoSimpleFoam的算例，可以先研究研究
