你好@random_ran :
我只针对我算的压气机的问题解释一下,我所指的转捩就是您所说的转捩,但是对于k-omega SST这些经典的湍流模型,对于叶轮机械,飞行器等流动问题,算出来的结果都是一个全湍流的状态。无法准确的预测转捩发生的位置和大小,而转捩对于气动性能,分离的影响不容忽视,因此有学者利用经验关系式,耦合流场当地变量来预测转捩发生的位置和大小。
小葱小虫
帖子
-
-
高马赫数,可压缩的
-
我想咨询一下openfoam有基于k-omega SST的转捩模型吗,类似于cfx和fluent中的。有的话如何设置那
-
@程迪 感谢!正在尝试
-
-
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$
-
我给的进口总温总压,出口背压,速度让它自己计算,速度是有方向的,不知道如何给稳定些
-
| ========= | | | \\ / 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; } } // ************************************************************************* //
-
| ========= | | | \\ / 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 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; // ************************************************************************* //
-
网格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; location "constant"; object turbulenceProperties; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // simulationType RAS; RAS { // Tested with kEpsilon, realizableKE, kOmega, kOmegaSST, v2f, // ShihQuadraticKE, LienCubicKE. RASModel kOmegaSST; turbulence on; printCoeffs on; } // ************************************************************************* //
-
介质就是空气
/*--------------------------------*- 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 volVectorField; object U; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 1 -1 0 0 0 0]; internalField uniform (126 113 0); boundaryField { hub { type fixedValue; value uniform (0 0 0); } inlet { type pressureDirectedInletVelocity; inletDirection uniform (0.743 0.669 0); value uniform (126 113 0); } outlet { type inletOutlet; inletValue $internalField; value uniform (63 55.5 0); } per1 { type cyclicAMI; } per2 { type cyclicAMI; } ps { type fixedValue; value uniform (0 0 0); } ss { type fixedValue; value uniform (0 0 0); } sym { type symmetryPlane; } } // ************************************************************************* //
-
| ========= | | | \\ / 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 volScalarField; object T; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 0 0 1 0 0 0]; internalField uniform 298; boundaryField { hub { type zeroGradient; } inlet { type totalTemperature; T0 uniform 300; gamma 1.4; } outlet { type zeroGradient; } per1 { type cyclicAMI; } per2 { type cyclicAMI; } ps { type zeroGradient; } ss { type zeroGradient; } sym { type symmetryPlane; } } // ************************************************************************* //
-
与上一个帖子一样,http://www.cfd-china.com/topic/744/请教关于进口条件的设置的问题/17
还是算的平面叶栅,只是改成可压缩的,但是总是挂掉,我贴出所有边界条件以及网格文件,希望大家给些建议。![替代文字]
弦长是40mm,不知道这是不是有问题Turbulence length scale,要不要成个0.07。| ========= | | | \\ / 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 volScalarField; location "0"; object k; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 2 -2 0 0 0 0]; internalField uniform 194; boundaryField { hub { type symmetryPlane; } inlet { type turbulentIntensityKineticEnergyInlet; intensity 0.05; value uniform 194; } outlet { type zeroGradient; } per1 { type cyclicAMI; value uniform 190; } per2 { type cyclicAMI; value uniform 190; } ps { type kqRWallFunction; value uniform 190; } ss { type kqRWallFunction; value uniform 190; } sym { type symmetryPlane; } } // ************************************************************************* // ```/*--------------------------------*- 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 volScalarField; location "0"; object epsilon; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 2 -3 0 0 0 0]; internalField uniform 6043; boundaryField { hub { type epsilonWallFunction; value uniform 6043; } inlet { type turbulentMixingLengthDissipationRateInlet; mixingLength 0.05; value uniform 6043; } outlet { type zeroGradient; } per1 { type cyclicAMI; } per2 { type cyclicAMI; } ps { type epsilonWallFunction; value uniform 6043; } ss { type epsilonWallFunction; value uniform 6043; ; } sym { type symmetryPlane; } } // ************************************************************************* // ```/*--------------------------------*- 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 volScalarField; object omega; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // omegaInlet 500; dimensions [0 0 -1 0 0 0 0]; internalField uniform $omegaInlet; boundaryField { hub { type omegaWallFunction; value $internalField; } inlet { type turbulentMixingLengthFrequencyInlet; mixingLength 0.05; k k; value $internalField; } outlet { type zeroGradient; } per1 { type cyclicAMI; value $internalField; } per2 { type cyclicAMI; value $internalField; } ps { type omegaWallFunction; value $internalField; } ss { type omegaWallFunction; value $internalField; } sym { type symmetryPlane; } } // ************************************************************************* // ```/*--------------------------------*- 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 volScalarField; object p; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [1 -1 -2 0 0 0 0]; internalField uniform 103325; boundaryField { hub { type zeroGradient; } inlet { type totalPressure; p0 uniform 119000; gamma 1.4; value uniform 119000; } outlet { type fixedValue; value uniform 101325; } per1 { type cyclicAMI; } per2 { type cyclicAMI; } ps { type zeroGradient; } ss { type zeroGradient; } sym { type symmetryPlane; } } // ************************************************************************* // ```/*--------------------------------*- 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 volScalarField; object T; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 0 0 1 0 0 0]; internalField uniform 298; boundaryField { hub { type zeroGradient; } inlet { type totalTemperature; T0 uniform 300; gamma 1.4; } outlet { type zeroGradient; } per1 { type cyclicAMI; } per2 { type cyclicAMI; } ps { type zeroGradient; } ss { type zeroGradient; } sym { type symmetryPlane; } } // ************************************************************************* //
-
@赵一铭 你好,我这个是从群里下载的,从来没有编译过,应该是直接可以用的那,这个怎么弄
-
cc@ubuntu:~/Downloads/rhoCentralLTS-OF-2.1$ wmake Making dependency list for source file rhoCentralLTS.C SOURCE=rhoCentralLTS.C ; g++ -m32 -Dlinux -DWM_DP -Wall -Wextra -Wno-unused-parameter -Wold-style-cast -O3 -DNoRepository -ftemplate-depth-100 -IBCs/lnInclude -I/opt/openfoam211/src/finiteVolume/lnInclude -I/opt/openfoam211/src/thermophysicalModels/basic/lnInclude -I/opt/openfoam211/src/thermophysicalModels/specie/lnInclude -I/opt/openfoam211/src/turbulenceModels/compressible/turbulenceModel -I/opt/openfoam211/src/dynamicMesh/lnInclude -I/opt/openfoam211/src/meshTools/lnInclude -IlnInclude -I. -I/opt/openfoam211/src/OpenFOAM/lnInclude -I/opt/openfoam211/src/OSspecific/POSIX/lnInclude -fPIC -c $SOURCE -o Make/linuxGccDPOpt/rhoCentralLTS.o In file included from createFields.H:35:0, from rhoCentralLTS.C:48: rhoBoundaryTypes.H: In function ‘int main(int, char**)’: rhoBoundaryTypes.H:12:36: error: ‘fixedRhoFvPatchScalarField’ has not been declared /opt/openfoam211/src/finiteVolume/lnInclude/readTimeControls.H:32:12: warning: unused variable ‘adjustTimeStep’ [-Wunused-variable] /opt/openfoam211/src/finiteVolume/lnInclude/readTimeControls.H:32:12: warning: unused variable ‘adjustTimeStep’ [-Wunused-variable] /opt/openfoam211/src/finiteVolume/lnInclude/readTimeControls.H:35:8: warning: unused variable ‘maxCo’ [-Wunused-variable] /opt/openfoam211/src/finiteVolume/lnInclude/readTimeControls.H:38:8: warning: unused variable ‘maxDeltaT’ [-Wunused-variable] make: *** [Make/linuxGccDPOpt/rhoCentralLTS.o] Error 1 cc@ubuntu:~/Downloads/rhoCentralLTS-OF-2.1$ cd .. cc@ubuntu:~/Downloads$ wmake wmake error: 'Make' directory does not exist in /home/cc/Downloads cc@ubuntu:~/Downloads$ cd rhoCentralLTS-OF-2.1/ cc@ubuntu:~/Downloads/rhoCentralLTS-OF-2.1$ wmake SOURCE=rhoCentralLTS.C ; g++ -m32 -Dlinux -DWM_DP -Wall -Wextra -Wno-unused-parameter -Wold-style-cast -O3 -DNoRepository -ftemplate-depth-100 -IBCs/lnInclude -I/opt/openfoam211/src/finiteVolume/lnInclude -I/opt/openfoam211/src/thermophysicalModels/basic/lnInclude -I/opt/openfoam211/src/thermophysicalModels/specie/lnInclude -I/opt/openfoam211/src/turbulenceModels/compressible/turbulenceModel -I/opt/openfoam211/src/dynamicMesh/lnInclude -I/opt/openfoam211/src/meshTools/lnInclude -IlnInclude -I. -I/opt/openfoam211/src/OpenFOAM/lnInclude -I/opt/openfoam211/src/OSspecific/POSIX/lnInclude -fPIC -c $SOURCE -o Make/linuxGccDPOpt/rhoCentralLTS.o In file included from createFields.H:35:0, from rhoCentralLTS.C:48: rhoBoundaryTypes.H: In function ‘int main(int, char**)’: rhoBoundaryTypes.H:12:36: error: ‘fixedRhoFvPatchScalarField’ has not been declared /opt/openfoam211/src/finiteVolume/lnInclude/readTimeControls.H:32:12: warning: unused variable ‘adjustTimeStep’ [-Wunused-variable] /opt/openfoam211/src/finiteVolume/lnInclude/readTimeControls.H:32:12: warning: unused variable ‘adjustTimeStep’ [-Wunused-variable] /opt/openfoam211/src/finiteVolume/lnInclude/readTimeControls.H:35:8: warning: unused variable ‘maxCo’ [-Wunused-variable] /opt/openfoam211/src/finiteVolume/lnInclude/readTimeControls.H:38:8: warning: unused variable ‘maxDeltaT’ [-Wunused-variable] make: *** [Make/linuxGccDPOpt/rhoCentralLTS.o] Error 1 cc@ubuntu:~/Downloads/rhoCentralLTS-OF-2.1$ wmake SOURCE=rhoCentralLTS.C ; g++ -m32 -Dlinux -DWM_DP -Wall -Wextra -Wno-unused-parameter -Wold-style-cast -O3 -DNoRepository -ftemplate-depth-100 -IBCs/lnInclude -I/opt/openfoam211/src/finiteVolume/lnInclude -I/opt/openfoam211/src/thermophysicalModels/basic/lnInclude -I/opt/openfoam211/src/thermophysicalModels/specie/lnInclude -I/opt/openfoam211/src/turbulenceModels/compressible/turbulenceModel -I/opt/openfoam211/src/dynamicMesh/lnInclude -I/opt/openfoam211/src/meshTools/lnInclude -IlnInclude -I. -I/opt/openfoam211/src/OpenFOAM/lnInclude -I/opt/openfoam211/src/OSspecific/POSIX/lnInclude -fPIC -c $SOURCE -o Make/linuxGccDPOpt/rhoCentralLTS.o In file included from createFields.H:35:0, from rhoCentralLTS.C:48: rhoBoundaryTypes.H: In function ‘int main(int, char**)’: rhoBoundaryTypes.H:12:36: error: ‘fixedRhoFvPatchScalarField’ has not been declared /opt/openfoam211/src/finiteVolume/lnInclude/readTimeControls.H:32:12: warning: unused variable ‘adjustTimeStep’ [-Wunused-variable] /opt/openfoam211/src/finiteVolume/lnInclude/readTimeControls.H:32:12: warning: unused variable ‘adjustTimeStep’ [-Wunused-variable] /opt/openfoam211/src/finiteVolume/lnInclude/readTimeControls.H:35:8: warning: unused variable ‘maxCo’ [-Wunused-variable] /opt/openfoam211/src/finiteVolume/lnInclude/readTimeControls.H:38:8: warning: unused variable ‘maxDeltaT’ [-Wunused-variable] make: *** [Make/linuxGccDPOpt/rhoCentralLTS.o] Error 1 cc@ubuntu:~/Downloads/rhoCentralLTS-OF-2.1$
-
@Izumi 是增加nNonOrthogonalCorrector这个数吗,谢谢我试试
-
@cfd-china 老师您好,由于我的算例进口马赫数为0.7,所以要考虑可压缩性,但是算着算着就发散了,用的rhosimplefoam,请问,能否帮我调试一下,由于我的上面是对称面,下面是壁面,所以网格10万很难满足壁面要求,能否帮忙调试,谢谢
-
@qjh888 请问这些求解器是在extend中的吗,还有可以计算稳态吗
-
交流群李又说rhocentralfoam可以,但为发现稳态的算例供参考
-
@qjh888 请问你算高马赫数的叶轮机械用的什么求解器,rhocentralfoam有稳态的吗
-
请问,openfoam有可以像fluent和cfx那样的插值功能吗,即用原来的算例(网格疏密有变化)做初场,算新的算例。是否也可以将fluent的结果作为初场那。
-
@cfd-china 我计算的是50%高度,另一半按照对称面设置的,所以给出的是一半叶片的高度。平面叶栅算出来分离形状都差不多,但of算出来的分离范围太小了(这是同比例的,不能局部放大)
-
@cfd-china :crying: 我出的壁面极限流线,实验的分离范围占了2/3,但是openfoam算的1/3都不到
-
我采用sst湍流模型计算平面叶栅,采用可压缩,rhosimpleFoam计算,进口速度0.22Ma,算出来的结果吸力面分离过小,不知是什么原因。采用同样模型cfx计算与实验吻合较好
实验参数设置见 http://www.cfd-china.com/topic/744/请教关于进口条件的设置的问题/24 -
| ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 3.0.1 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ Build : 3.0.1-119cac7e8750 Exec : rhoSimpleFoam Date : Mar 30 2017 Time : 04:56:26 Host : "ubuntu" PID : 6792 Case : /home/cc/OpenFOAM/cc-3.0.1/run/tutorials/compressible/rhoSimpleFoam/b-sysm nProcs : 1 sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). fileModificationChecking : Monitoring run-time modified files using timeStampMaster allowSystemOperations : Allowing user-supplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create mesh for time = 0 SIMPLE: convergence criteria field p tolerance 1e-05 field U tolerance 1e-05 field e tolerance 1e-05 field "(k|epsilon|omega)" tolerance 1e-05 Reading thermophysical properties Selecting thermodynamics package { type hePsiThermo; mixture pureMixture; transport sutherland; thermo hConst; equationOfState perfectGas; specie specie; energy sensibleInternalEnergy; } AMI: Creating addressing and weights between 6006 source faces and 6006 target faces AMI: Patch source sum(weights) min/max/average = 0.999619, 1.00005, 0.999982 AMI: Patch target sum(weights) min/max/average = 0.999924, 1.00042, 1.00003 Reading field U Reading/calculating face flux field phi Creating turbulence model Selecting RAS turbulence model kEpsilon kEpsilonCoeffs { Cmu 0.09; C1 1.44; C2 1.92; C3 -0.33; sigmak 1; sigmaEps 1.3; } No MRF models present No finite volume options present Starting time loop Time = 1 smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 0.069783, No Iterations 4 smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 0.0513316, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 0.0997659, No Iterations 4 smoothSolver: Solving for e, Initial residual = 1, Final residual = 0.0104707, No Iterations 4 GAMG: Solving for p, Initial residual = 1, Final residual = 0.0959957, No Iterations 5 time step continuity errors : sum local = 4.20761, global = 1.18992, cumulative = 1.18992 rho max/min : 1.19542 1.17035 smoothSolver: Solving for epsilon, Initial residual = 0.0778234, Final residual = 0.0010543, No Iterations 2 smoothSolver: Solving for k, Initial residual = 1, Final residual = 0.0337686, No Iterations 2 ExecutionTime = 6.36 s ClockTime = 7 s Time = 2 smoothSolver: Solving for Ux, Initial residual = 0.428754, Final residual = 0.032496, No Iterations 4 smoothSolver: Solving for Uy, Initial residual = 0.638559, Final residual = 0.0467907, No Iterations 4 smoothSolver: Solving for Uz, Initial residual = 0.635954, Final residual = 0.0516729, No Iterations 4 smoothSolver: Solving for e, Initial residual = 0.0653834, Final residual = 0.000793882, No Iterations 4 GAMG: Solving for p, Initial residual = 0.878525, Final residual = 0.0461803, No Iterations 2 time step continuity errors : sum local = 10.1308, global = -3.28196, cumulative = -2.09205 rho max/min : 1.18866 1.16877 smoothSolver: Solving for epsilon, Initial residual = 0.00247931, Final residual = 3.84594e-05, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.532281, Final residual = 0.0125941, No Iterations 2 ExecutionTime = 8.76 s ClockTime = 9 s Time = 3 smoothSolver: Solving for Ux, Initial residual = 0.772756, Final residual = 4.26156e-07, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.586361, Final residual = 6.45035e-06, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.999934, Final residual = 3.78056e-06, No Iterations 2 smoothSolver: Solving for e, Initial residual = 0.895471, Final residual = 1.25174e-08, No Iterations 2 GAMG: Solving for p, Initial residual = 0.753867, Final residual = 0.0690092, No Iterations 2 time step continuity errors : sum local = 13.8686, global = -3.34011, cumulative = -5.43216 rho max/min : 1.18714 1.17017 smoothSolver: Solving for epsilon, Initial residual = 0.0100085, Final residual = 2.52959e-08, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.999995, Final residual = 4.04155e-08, No Iterations 2 ExecutionTime = 10.99 s ClockTime = 11 s Time = 4 smoothSolver: Solving for Ux, Initial residual = 0.293869, Final residual = 0.0235438, No Iterations 4 smoothSolver: Solving for Uy, Initial residual = 0.102553, Final residual = 0.00811364, No Iterations 4 smoothSolver: Solving for Uz, Initial residual = 0.233767, Final residual = 0.0189625, No Iterations 4 smoothSolver: Solving for e, Initial residual = 0.0705048, Final residual = 0.00567789, No Iterations 2 GAMG: Solving for p, Initial residual = 0.658214, Final residual = 0.0621858, No Iterations 3 time step continuity errors : sum local = 9.7774, global = 2.79685, cumulative = -2.63531 rho max/min : 1.188 1.17357 smoothSolver: Solving for epsilon, Initial residual = 0.00186519, Final residual = 5.42206e-05, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.134148, Final residual = 0.00401518, No Iterations 2 ExecutionTime = 13.37 s ClockTime = 14 s Time = 5 smoothSolver: Solving for Ux, Initial residual = 0.205107, Final residual = 0.0165675, No Iterations 4 smoothSolver: Solving for Uy, Initial residual = 0.0770089, Final residual = 0.0055433, No Iterations 4 smoothSolver: Solving for Uz, Initial residual = 0.191079, Final residual = 0.0154472, No Iterations 4 smoothSolver: Solving for e, Initial residual = 0.0379826, Final residual = 0.0032221, No Iterations 2 GAMG: Solving for p, Initial residual = 0.538864, Final residual = 0.0453438, No Iterations 3 time step continuity errors : sum local = 5.75337, global = 1.2244, cumulative = -1.41091 rho max/min : 1.18967 1.17214 smoothSolver: Solving for epsilon, Initial residual = 0.00158084, Final residual = 4.74005e-05, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.101466, Final residual = 0.00294749, No Iterations 2 ExecutionTime = 15.82 s ClockTime = 16 s Time = 6 smoothSolver: Solving for Ux, Initial residual = 0.16938, Final residual = 0.0129302, No Iterations 4 smoothSolver: Solving for Uy, Initial residual = 0.275953, Final residual = 0.0202823, No Iterations 4 smoothSolver: Solving for Uz, Initial residual = 0.25257, Final residual = 0.020254, No Iterations 4 smoothSolver: Solving for e, Initial residual = 0.353213, Final residual = 0.00210383, No Iterations 2 GAMG: Solving for p, Initial residual = 0.384399, Final residual = 0.0268333, No Iterations 4 time step continuity errors : sum local = 2.92741, global = 0.637186, cumulative = -0.773722 rho max/min : 1.19094 1.1677 smoothSolver: Solving for epsilon, Initial residual = 0.00161393, Final residual = 4.37915e-05, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.0916599, Final residual = 0.00156534, No Iterations 2 ExecutionTime = 18.34 s ClockTime = 19 s Time = 7 smoothSolver: Solving for Ux, Initial residual = 0.138107, Final residual = 0.0104691, No Iterations 4 smoothSolver: Solving for Uy, Initial residual = 0.386117, Final residual = 0.0304146, No Iterations 4 smoothSolver: Solving for Uz, Initial residual = 0.376104, Final residual = 0.021135, No Iterations 4 smoothSolver: Solving for e, Initial residual = 0.282986, Final residual = 0.00235044, No Iterations 2 GAMG: Solving for p, Initial residual = 0.377236, Final residual = 0.0276733, No Iterations 6 time step continuity errors : sum local = 2.70399, global = 0.639034, cumulative = -0.134688 rho max/min : 1.19184 1.1637 smoothSolver: Solving for epsilon, Initial residual = 0.00161466, Final residual = 4.03465e-05, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.0392622, Final residual = 0.00112142, No Iterations 2 ExecutionTime = 21.03 s ClockTime = 21 s Time = 8 smoothSolver: Solving for Ux, Initial residual = 0.10629, Final residual = 0.00759966, No Iterations 4 smoothSolver: Solving for Uy, Initial residual = 0.172133, Final residual = 0.0105368, No Iterations 4 smoothSolver: Solving for Uz, Initial residual = 0.345993, Final residual = 0.0228648, No Iterations 4 smoothSolver: Solving for e, Initial residual = 0.300534, Final residual = 0.00107786, No Iterations 2 GAMG: Solving for p, Initial residual = 0.317855, Final residual = 0.0222564, No Iterations 7 time step continuity errors : sum local = 1.98611, global = -0.444252, cumulative = -0.578941 rho max/min : 1.1925 1.1601 smoothSolver: Solving for epsilon, Initial residual = 0.00146714, Final residual = 3.61448e-05, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.0664401, Final residual = 0.00082753, No Iterations 2 ExecutionTime = 23.84 s ClockTime = 24 s Time = 9 smoothSolver: Solving for Ux, Initial residual = 0.0989662, Final residual = 0.0078023, No Iterations 4 smoothSolver: Solving for Uy, Initial residual = 0.222443, Final residual = 0.0175467, No Iterations 4 smoothSolver: Solving for Uz, Initial residual = 0.275044, Final residual = 0.0222035, No Iterations 4 smoothSolver: Solving for e, Initial residual = 0.280869, Final residual = 0.00125374, No Iterations 2 GAMG: Solving for p, Initial residual = 0.358093, Final residual = 0.0228239, No Iterations 5 time step continuity errors : sum local = 1.78002, global = -0.281345, cumulative = -0.860286 rho max/min : 1.19322 1.15686 smoothSolver: Solving for epsilon, Initial residual = 0.00137547, Final residual = 3.39391e-05, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.050327, Final residual = 0.000665553, No Iterations 2 ExecutionTime = 26.45 s ClockTime = 27 s Time = 10 smoothSolver: Solving for Ux, Initial residual = 0.0931371, Final residual = 0.00748747, No Iterations 4 smoothSolver: Solving for Uy, Initial residual = 0.151795, Final residual = 0.0123852, No Iterations 4 smoothSolver: Solving for Uz, Initial residual = 0.279421, Final residual = 0.0209241, No Iterations 4 smoothSolver: Solving for e, Initial residual = 0.0364338, Final residual = 0.00354662, No Iterations 2 GAMG: Solving for p, Initial residual = 0.31275, Final residual = 0.0300344, No Iterations 4 time step continuity errors : sum local = 2.28464, global = -0.443375, cumulative = -1.30366 rho max/min : 1.19397 1.15395 smoothSolver: Solving for epsilon, Initial residual = 0.0027927, Final residual = 3.47056e-05, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.0227749, Final residual = 0.000590414, No Iterations 2 ExecutionTime = 28.98 s ClockTime = 29 s Time = 11 smoothSolver: Solving for Ux, Initial residual = 0.0797475, Final residual = 0.00640376, No Iterations 4 smoothSolver: Solving for Uy, Initial residual = 0.0571989, Final residual = 0.00457257, No Iterations 4 smoothSolver: Solving for Uz, Initial residual = 0.42095, Final residual = 0.0153136, No Iterations 4 smoothSolver: Solving for e, Initial residual = 0.0825759, Final residual = 0.00390404, No Iterations 2 GAMG: Solving for p, Initial residual = 0.319658, Final residual = 0.0280046, No Iterations 5 time step continuity errors : sum local = 2.15133, global = 0.43661, cumulative = -0.867051 rho max/min : 1.20824 1.15132 smoothSolver: Solving for epsilon, Initial residual = 0.0788348, Final residual = 0.000295584, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.132764, Final residual = 0.00065471, No Iterations 2 ExecutionTime = 31.58 s ClockTime = 32 s Time = 12 smoothSolver: Solving for Ux, Initial residual = 0.788096, Final residual = 0.00198053, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.933599, Final residual = 0.00158327, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.997667, Final residual = 0.000460968, No Iterations 2 smoothSolver: Solving for e, Initial residual = 0.205344, Final residual = 0.000880161, No Iterations 2 GAMG: Solving for p, Initial residual = 0.374334, Final residual = 0.0310307, No Iterations 7 time step continuity errors : sum local = 2.63811, global = 0.0465569, cumulative = -0.820494 rho max/min : 1.26961 1.05412 smoothSolver: Solving for epsilon, Initial residual = 0.390939, Final residual = 0.00159594, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.783255, Final residual = 0.000409466, No Iterations 2 ExecutionTime = 34.25 s ClockTime = 34 s Time = 13 smoothSolver: Solving for Ux, Initial residual = 0.604962, Final residual = 0.0183882, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.721848, Final residual = 0.00919717, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.826181, Final residual = 0.0156051, No Iterations 2 smoothSolver: Solving for e, Initial residual = 0.988811, Final residual = 0.0208175, 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)
-
@小葱小虫 我只是换了一套网格,其他没动
-
@cfd-china 我的算例可以算了,但是不知道为什么,壁面的分离算的特别小
-
@cfd-china 使用cyclic会有这个问题,因为两侧网格并不完全匹配
/*---------------------------------------------------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 3.0.1 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ Build : 3.0.1-119cac7e8750 Exec : simpleFoam Date : Mar 25 2017 Time : 05:33:37 Host : "ubuntu" PID : 48984 Case : /home/cc/OpenFOAM/cc-3.0.1/run/tutorials/incompressible/simpleFoam/cynasa nProcs : 1 sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). fileModificationChecking : Monitoring run-time modified files using timeStampMaster allowSystemOperations : Allowing user-supplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create mesh for time = 0 --> FOAM FATAL ERROR: face 47 area does not match neighbour by 10.4554% -- possible face ordering problem. patch:per1 my area:1.63789e-10 neighbour area:1.47515e-10 matching tolerance:0.1 Mesh face:5421439 fc:(0.0411928 0.00446869 5e-07) Neighbour fc:(0.040347 0.0264296 5e-07) If you are certain your matching is correct you can increase the 'matchTolerance' setting in the patch dictionary in the boundary file. Rerun with cyclic debug flag set for more information. From function cyclicPolyPatch::calcTransforms() in file meshes/polyMesh/polyPatches/constraint/cyclic/cyclicPolyPatch.C at line 221. FOAM exiting
-
@cfd-china 我感觉这个我理解的可能有问题,是应该给平移距离的大小(0.022m),还是(0,1,0)这么给。当时从这里看见的
separationVector (0 -0.022 0);
-
@cfd-china 好的,我这么设置周期边界没有问题吧,两个面的距离是22mm
| ========= | | | \\ / 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 30748; startFace 5380555; } inlet { type patch; nFaces 5546; startFace 5411303; } outlet { type patch; nFaces 4543; startFace 5416849; } per1 { type cyclicAMI; nFaces 15340; startFace 5421392; neighbourPatch per2; matchTolerance 0.01; transform translational; separationVector (0 0.022 0); } per2 { type cyclicAMI; nFaces 15340; startFace 5436732; neighbourPatch per1; matchTolerance 0.01; transform translational; separationVector (0 -0.022 0); } ps { type wall; nFaces 10266; startFace 5452072; } ss { type wall; nFaces 11151; startFace 5462338; } sym { type symmetryPlane; nFaces 30748; startFace 5473489; } ) // ************************************************************************* //
-
这是最新算的几部,麻烦老师们给看一下
| ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 3.0.1 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ Build : 3.0.1-119cac7e8750 Exec : simpleFoam Date : Mar 25 2017 Time : 02:33:53 Host : "ubuntu" PID : 44090 Case : /home/cc/OpenFOAM/cc-3.0.1/run/tutorials/incompressible/simpleFoam/nasa65 nProcs : 1 sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). fileModificationChecking : Monitoring run-time modified files using timeStampMaster allowSystemOperations : Allowing user-supplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create mesh for time = 0 SIMPLE: convergence criteria field U tolerance 1e-05 field p tolerance 1e-05 field "(k|epsilon|omega)" tolerance 1e-05 Reading field p AMI: Creating addressing and weights between 15340 source faces and 15340 target faces AMI: Patch source sum(weights) min/max/average = 0.999911, 1.00059, 1.00001 AMI: Patch target sum(weights) min/max/average = 0.999839, 1.00646, 1.00006 Reading field U Reading/calculating face flux field phi Selecting incompressible transport model Newtonian Selecting turbulence model type RAS Selecting RAS turbulence model kOmegaSST Selecting patchDistMethod meshWave kOmegaSSTCoeffs { alphaK1 0.85; alphaK2 1; alphaOmega1 0.5; alphaOmega2 0.856; gamma1 0.555556; gamma2 0.44; beta1 0.075; beta2 0.0828; betaStar 0.09; a1 0.31; b1 1; c1 10; F3 false; } No MRF models present No finite volume options present Starting time loop Time = 1 smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 0.0813711, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 0.0714855, No Iterations 5 smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 0.0771402, No Iterations 4 GAMG: Solving for p, Initial residual = 1, Final residual = 0.0876464, No Iterations 5 time step continuity errors : sum local = 10.8537, global = 4.01707, cumulative = 4.01707 smoothSolver: Solving for omega, Initial residual = 0.259035, Final residual = 0.0221475, No Iterations 3 smoothSolver: Solving for k, Initial residual = 1, Final residual = 0.0890898, No Iterations 5 ExecutionTime = 34.31 s ClockTime = 36 s Time = 2 smoothSolver: Solving for Ux, Initial residual = 0.886786, Final residual = 0.058124, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.979345, Final residual = 0.0631154, No Iterations 6 smoothSolver: Solving for Uz, Initial residual = 0.810242, Final residual = 0.0670691, No Iterations 6 GAMG: Solving for p, Initial residual = 0.986904, Final residual = 0.0565546, No Iterations 4 time step continuity errors : sum local = 1345.89, global = -10.0278, cumulative = -6.0107 smoothSolver: Solving for omega, Initial residual = 0.0850778, Final residual = 0.00220061, No Iterations 1 smoothSolver: Solving for k, Initial residual = 0.814193, Final residual = 0.0238861, No Iterations 2 ExecutionTime = 45.36 s ClockTime = 47 s Time = 3 smoothSolver: Solving for Ux, Initial residual = 0.0193303, Final residual = 0.00152468, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.011814, Final residual = 0.00102405, No Iterations 6 smoothSolver: Solving for Uz, Initial residual = 0.422504, Final residual = 0.0419746, No Iterations 6 GAMG: Solving for p, Initial residual = 0.791662, Final residual = 0.0642592, No Iterations 4 time step continuity errors : sum local = 957.279, global = -143.421, cumulative = -149.432 smoothSolver: Solving for omega, Initial residual = 0.00274129, Final residual = 4.02132e-05, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.169936, Final residual = 0.0157838, No Iterations 1 ExecutionTime = 56.05 s ClockTime = 58 s Time = 4 smoothSolver: Solving for Ux, Initial residual = 0.0983965, Final residual = 0.00679746, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.0657851, Final residual = 0.00457164, No Iterations 6 smoothSolver: Solving for Uz, Initial residual = 0.206983, Final residual = 0.0189299, No Iterations 6 GAMG: Solving for p, Initial residual = 0.980449, Final residual = 0.0634521, No Iterations 4 time step continuity errors : sum local = 1261.26, global = -92.2404, cumulative = -241.672 smoothSolver: Solving for omega, Initial residual = 0.000609973, Final residual = 3.54731e-05, No Iterations 1 smoothSolver: Solving for k, Initial residual = 0.10451, Final residual = 0.0098634, No Iterations 1 ExecutionTime = 66.79 s ClockTime = 69 s Time = 5 smoothSolver: Solving for Ux, Initial residual = 0.0261376, Final residual = 0.00214716, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.0187723, Final residual = 0.00145964, No Iterations 6 smoothSolver: Solving for Uz, Initial residual = 0.0626224, Final residual = 0.00595206, No Iterations 6 GAMG: Solving for p, Initial residual = 0.921146, Final residual = 0.0735434, No Iterations 4 time step continuity errors : sum local = 1305.03, global = -174.17, cumulative = -415.842 smoothSolver: Solving for omega, Initial residual = 0.000273909, Final residual = 9.03935e-06, No Iterations 1 smoothSolver: Solving for k, Initial residual = 0.0349887, Final residual = 0.00218001, No Iterations 3 ExecutionTime = 78.4 s ClockTime = 81 s Time = 6 smoothSolver: Solving for Ux, Initial residual = 0.0566235, Final residual = 0.00458913, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.0400099, Final residual = 0.00302616, No Iterations 6 smoothSolver: Solving for Uz, Initial residual = 0.0384566, Final residual = 0.00324406, No Iterations 6 GAMG: Solving for p, Initial residual = 0.916166, Final residual = 0.0603851, No Iterations 4 time step continuity errors : sum local = 668.506, global = -94.6009, cumulative = -510.443 smoothSolver: Solving for omega, Initial residual = 0.000292328, Final residual = 1.62378e-05, No Iterations 1 smoothSolver: Solving for k, Initial residual = 0.046623, Final residual = 0.00262596, No Iterations 2 ExecutionTime = 89.53 s ClockTime = 92 s Time = 7 smoothSolver: Solving for Ux, Initial residual = 0.233958, Final residual = 0.0176323, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.128911, Final residual = 0.00936957, No Iterations 6 smoothSolver: Solving for Uz, Initial residual = 0.0442786, Final residual = 0.00376373, No Iterations 6 GAMG: Solving for p, Initial residual = 0.952315, Final residual = 0.0728294, No Iterations 4 time step continuity errors : sum local = 1263.9, global = -195.294, cumulative = -705.737 smoothSolver: Solving for omega, Initial residual = 0.000206228, Final residual = 1.22213e-05, No Iterations 1 smoothSolver: Solving for k, Initial residual = 0.0167436, Final residual = 0.00109875, No Iterations 4 ExecutionTime = 100.52 s ClockTime = 103 s Time = 8 smoothSolver: Solving for Ux, Initial residual = 0.064942, Final residual = 0.00538914, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.0316677, Final residual = 0.00284742, No Iterations 6 smoothSolver: Solving for Uz, Initial residual = 0.0769102, Final residual = 0.00692843, No Iterations 6 GAMG: Solving for p, Initial residual = 0.791177, Final residual = 0.0625584, No Iterations 6 time step continuity errors : sum local = 645.649, global = -319.458, cumulative = -1025.19 smoothSolver: Solving for omega, Initial residual = 0.000258477, Final residual = 1.59141e-05, No Iterations 1 smoothSolver: Solving for k, Initial residual = 0.0181236, Final residual = 0.00171421, No Iterations 3 ExecutionTime = 113.31 s ClockTime = 116 s Time = 9 smoothSolver: Solving for Ux, Initial residual = 0.219725, Final residual = 0.0183762, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.158121, Final residual = 0.0111089, No Iterations 7 smoothSolver: Solving for Uz, Initial residual = 0.303528, Final residual = 0.0257045, No Iterations 6 GAMG: Solving for p, Initial residual = 0.94797, Final residual = 0.064793, No Iterations 4 time step continuity errors : sum local = 1113.34, global = -309.45, cumulative = -1334.64 smoothSolver: Solving for omega, Initial residual = 0.000122954, Final residual = 1.00576e-05, No Iterations 1 smoothSolver: Solving for k, Initial residual = 0.0197785, Final residual = 0.00176744, No Iterations 3 ExecutionTime = 126.36 s ClockTime = 129 s Time = 10 smoothSolver: Solving for Ux, Initial residual = 0.0821976, Final residual = 0.00721542, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.03459, Final residual = 0.00247223, No Iterations 7 smoothSolver: Solving for Uz, Initial residual = 0.1429, Final residual = 0.0115741, No Iterations 6 GAMG: Solving for p, Initial residual = 0.861997, Final residual = 0.084008, No Iterations 5 time step continuity errors : sum local = 1083.47, global = 389.776, cumulative = -944.868 smoothSolver: Solving for omega, Initial residual = 0.000136441, Final residual = 2.43071e-06, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.0235868, Final residual = 0.00172782, No Iterations 3 ExecutionTime = 138.69 s ClockTime = 141 s Time = 11 smoothSolver: Solving for Ux, Initial residual = 0.163168, Final residual = 0.0136899, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.0793402, Final residual = 0.00752288, No Iterations 5 smoothSolver: Solving for Uz, Initial residual = 0.386931, Final residual = 0.0294783, No Iterations 6 GAMG: Solving for p, Initial residual = 0.870129, Final residual = 0.0711586, No Iterations 6 time step continuity errors : sum local = 919.618, global = 516.016, cumulative = -428.851 smoothSolver: Solving for omega, Initial residual = 0.000106145, Final residual = 1.01946e-05, No Iterations 1 smoothSolver: Solving for k, Initial residual = 0.0196961, Final residual = 0.0018465, No Iterations 3 ExecutionTime = 151.28 s ClockTime = 154 s Time = 12 smoothSolver: Solving for Ux, Initial residual = 0.274151, Final residual = 0.0245354, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.186447, Final residual = 0.0165338, No Iterations 5 smoothSolver: Solving for Uz, Initial residual = 0.156144, Final residual = 0.0136401, No Iterations 6 GAMG: Solving for p, Initial residual = 0.887478, Final residual = 0.070432, No Iterations 6 time step continuity errors : sum local = 1119.32, global = 661.637, cumulative = 232.786 smoothSolver: Solving for omega, Initial residual = 0.000100293, Final residual = 2.24457e-06, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.0248731, Final residual = 0.00212404, No Iterations 3 ExecutionTime = 164.01 s ClockTime = 167 s Time = 13 smoothSolver: Solving for Ux, Initial residual = 0.205525, Final residual = 0.0187076, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.0923904, Final residual = 0.00894338, No Iterations 5 smoothSolver: Solving for Uz, Initial residual = 0.196898, Final residual = 0.0175523, No Iterations 6 GAMG: Solving for p, Initial residual = 0.817651, Final residual = 0.0790382, No Iterations 8 time step continuity errors : sum local = 1128.57, global = 728.044, cumulative = 960.83 smoothSolver: Solving for omega, Initial residual = 0.00010952, Final residual = 2.80076e-06, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.0218195, Final residual = 0.00158256, No Iterations 4 ExecutionTime = 178.57 s ClockTime = 181 s Time = 14 smoothSolver: Solving for Ux, Initial residual = 0.209248, Final residual = 0.0174402, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.234617, Final residual = 0.0220479, No Iterations 6 smoothSolver: Solving for Uz, Initial residual = 0.330718, Final residual = 0.0281098, No Iterations 6 GAMG: Solving for p, Initial residual = 0.871195, Final residual = 0.0759581, No Iterations 10 time step continuity errors : sum local = 1172.79, global = 858.035, cumulative = 1818.86 smoothSolver: Solving for omega, Initial residual = 0.000111802, Final residual = 3.26755e-06, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.0256925, Final residual = 0.00171392, No Iterations 4 ExecutionTime = 193.33 s ClockTime = 196 s Time = 15 smoothSolver: Solving for Ux, Initial residual = 0.173376, Final residual = 0.0142151, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.141625, Final residual = 0.0112524, No Iterations 7 smoothSolver: Solving for Uz, Initial residual = 0.203415, Final residual = 0.0168199, No Iterations 6 GAMG: Solving for p, Initial residual = 0.850352, Final residual = 0.0820699, No Iterations 8 time step continuity errors : sum local = 1261.82, global = 852.091, cumulative = 2670.96 smoothSolver: Solving for omega, Initial residual = 0.000119946, Final residual = 3.04546e-06, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.0523867, Final residual = 0.00451596, No Iterations 2 ExecutionTime = 206.27 s ClockTime = 209 s Time = 16 smoothSolver: Solving for Ux, Initial residual = 0.189402, Final residual = 0.0147132, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.172538, Final residual = 0.0148313, No Iterations 6 smoothSolver: Solving for Uz, Initial residual = 0.216101, Final residual = 0.0175322, No Iterations 6 GAMG: Solving for p, Initial residual = 0.831111, Final residual = 0.0800892, No Iterations 13 time step continuity errors : sum local = 1124.03, global = -888.775, cumulative = 1782.18 smoothSolver: Solving for omega, Initial residual = 0.000132821, Final residual = 2.71341e-06, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.0429148, Final residual = 0.00279213, No Iterations 3 ExecutionTime = 220.4 s ClockTime = 223 s Time = 17 smoothSolver: Solving for Ux, Initial residual = 0.240906, Final residual = 0.0185245, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.481738, Final residual = 0.0406032, No Iterations 5 smoothSolver: Solving for Uz, Initial residual = 0.306464, Final residual = 0.0227411, No Iterations 6 GAMG: Solving for p, Initial residual = 0.829137, Final residual = 0.0797512, No Iterations 5 time step continuity errors : sum local = 1176.95, global = 495.085, cumulative = 2277.27 smoothSolver: Solving for omega, Initial residual = 0.000148683, Final residual = 1.27859e-05, No Iterations 1 smoothSolver: Solving for k, Initial residual = 0.0542324, Final residual = 0.0033722, No Iterations 3 ExecutionTime = 231.68 s ClockTime = 235 s Time = 18 smoothSolver: Solving for Ux, Initial residual = 0.206106, Final residual = 0.0161844, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.255402, Final residual = 0.0223374, No Iterations 5 smoothSolver: Solving for Uz, Initial residual = 0.341787, Final residual = 0.0275786, No Iterations 6 GAMG: Solving for p, Initial residual = 0.827099, Final residual = 0.082525, No Iterations 7 time step continuity errors : sum local = 1134.39, global = 819.278, cumulative = 3096.54 smoothSolver: Solving for omega, Initial residual = 9.46476e-05, Final residual = 2.25188e-06, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.0505738, Final residual = 0.00348548, No Iterations 3 ExecutionTime = 242.31 s ClockTime = 246 s Time = 19 smoothSolver: Solving for Ux, Initial residual = 0.199931, Final residual = 0.0167285, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.361021, Final residual = 0.0275666, No Iterations 6 smoothSolver: Solving for Uz, Initial residual = 0.234612, Final residual = 0.0186392, No Iterations 6 GAMG: Solving for p, Initial residual = 0.829252, Final residual = 0.0781392, No Iterations 18 time step continuity errors : sum local = 1182.87, global = 1048.32, cumulative = 4144.86 smoothSolver: Solving for omega, Initial residual = 0.000103443, Final residual = 9.28303e-06, No Iterations 1 smoothSolver: Solving for k, Initial residual = 0.160832, Final residual = 0.0122526, No Iterations 1 ExecutionTime = 255.21 s ClockTime = 259 s Time = 20 smoothSolver: Solving for Ux, Initial residual = 0.203207, Final residual = 0.0177404, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.328249, Final residual = 0.0269088, No Iterations 6 smoothSolver: Solving for Uz, Initial residual = 0.224519, Final residual = 0.0164899, No Iterations 6 GAMG: Solving for p, Initial residual = 0.838448, Final residual = 0.0774968, No Iterations 16 time step continuity errors : sum local = 1282.19, global = 1084.48, cumulative = 5229.34 smoothSolver: Solving for omega, Initial residual = 0.000137159, Final residual = 1.18779e-05, No Iterations 1 smoothSolver: Solving for k, Initial residual = 0.105212, Final residual = 0.00733032, No Iterations 2 ExecutionTime = 268.59 s ClockTime = 273 s Time = 21 smoothSolver: Solving for Ux, Initial residual = 0.195891, Final residual = 0.0166656, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.295498, Final residual = 0.0264212, No Iterations 6 smoothSolver: Solving for Uz, Initial residual = 0.256556, Final residual = 0.019235, No Iterations 6 GAMG: Solving for p, Initial residual = 0.820159, Final residual = 0.0809549, No Iterations 27 time step continuity errors : sum local = 1210.13, global = -1110.62, cumulative = 4118.72 smoothSolver: Solving for omega, Initial residual = 0.000209472, Final residual = 1.33081e-05, No Iterations 1 smoothSolver: Solving for k, Initial residual = 0.0963909, Final residual = 0.00779254, No Iterations 2 ExecutionTime = 283.95 s ClockTime = 288 s Time = 22 smoothSolver: Solving for Ux, Initial residual = 0.210533, Final residual = 0.0180219, No Iterations 5 smoothSolver: Solving for Uy, Initial residual = 0.451871, Final residual = 0.0368704, No Iterations 6 smoothSolver: Solving for Uz, Initial residual = 0.267514, Final residual = 0.0242343, No Iterations 5 GAMG: Solving for p, Initial residual = 0.828406, Final residual = 0.0758993, No Iterations 30 time step continuity errors : sum local = 1183.8, global = 1080.35, cumulative = 5199.06 smoothSolver: Solving for omega, Initial residual = 0.00040392, Final residual = 1.58342e-05, No Iterations 1 smoothSolver: Solving for k, Initial residual = 0.116199, Final residual = 0.00839545, No Iterations 2 ExecutionTime = 299.9 s ClockTime = 304 s Time = 23 smoothSolver: Solving for Ux, Initial residual = 0.223249, Final residual = 0.0207695, No Iterations 5 smoothSolver: Solving for Uy, Initial residual = 0.479522, Final residual = 0.0438371, No Iterations 5 smoothSolver: Solving for Uz, Initial residual = 0.287509, Final residual = 0.0233311, No Iterations 5 GAMG: Solving for p, Initial residual = 0.819938, Final residual = 0.0816971, No Iterations 21 time step continuity errors : sum local = 1364, global = -1206.8, cumulative = 3992.26 smoothSolver: Solving for omega, Initial residual = 0.000308337, Final residual = 2.31198e-05, No Iterations 1 smoothSolver: Solving for k, Initial residual = 0.183739, Final residual = 0.0103415, No Iterations 2 ExecutionTime = 314.31 s ClockTime = 319 s
-
@random_ran 我开始没有开湍流模型,想算个初场,2000部之前都是没有湍流模型的,omega用这个公式算的
omega = k^0.5/(C^0.25 * l) c=0.09
我的是进口条件是进口已知总压,和气流方向出口知道静压,算出来速度大约是0.67Ma -
| ========= | | | \\ / 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 volScalarField; object omega; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // omegaInlet 630.0; dimensions [0 0 -1 0 0 0 0]; internalField uniform $omegaInlet; boundaryField { hub { type omegaWallFunction; value $internalField; } inlet { type turbulentMixingLengthFrequencyInlet; mixingLength 0.025; k k; value $internalField; } outlet { type zeroGradient; } per1 { type cyclicAMI; } per2 { type cyclicAMI; } ps { type omegaWallFunction; value $internalField; } ss { type omegaWallFunction; value $internalField; } sym { type symmetryPlane; } } // ************************************************************************* //
-
| ========= | | | \\ / 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 volScalarField; location "0"; object k; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 2 -2 0 0 0 0]; internalField uniform 194.6; boundaryField { hub { type kqRWallFunction; value uniform 194.6; } inlet { type turbulentIntensityKineticEnergyInlet; intensity 0.05; U U; phi phi; value uniform 194.6; } outlet { type zeroGradient; } per1 { type cyclicAMI; } per2 { type cyclicAMI; } ps { type kqRWallFunction; value uniform 194.6; } ss { type kqRWallFunction; value uniform 194.6; } sym { type symmetryPlane; } } // ************************************************************************* //
-
@cfd-china 您看有什么问题吗
| ========= | | | \\ / 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 volVectorField; object U; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 1 -1 0 0 0 0]; internalField uniform (126 113 0); boundaryField { hub { type fixedValue; value uniform (0 0 0); } inlet { type pressureDirectedInletVelocity; inletDirection uniform (0.743 0.669 0); value uniform (0 0 0); } outlet { type zeroGradient; } per1 { type cyclicAMI; } per2 { type cyclicAMI; } ps { type fixedValue; value uniform (0 0 0); } ss { type fixedValue; value uniform (0 0 0); } sym { type symmetryPlane; } } // ************************************************************************* //
-
| ========= | | | \\ / 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 volScalarField; object p; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 2 -2 0 0 0 0]; internalField uniform 122825; boundaryField { hub { type zeroGradient; } inlet { type totalPressure; p0 uniform 132700; value uniform 132700; gamma 1.4; } outlet { type fixedValue; value uniform 111500; } per1 { type cyclicAMI; } per2 { type cyclicAMI; } ps { type zeroGradient; } ss { type zeroGradient; } sym { type symmetryPlane; } }
-
@cfd-china
Failed 1 mesh checks. 这个是不是网格不好的意思
-
-
| ========= | | | \\ / 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 { efault none; div(phi,U) bounded Gauss limitedLinearV 1; div(phi,k) bounded Gauss limitedLinear 1; div(phi,epsilon) bounded Gauss limitedLinear 1; div(phi,R) bounded Gauss limitedLinear 1; div(R) Gauss linear; div(phi,nuTilda) bounded Gauss limitedLinear 1; div((nuEff*dev2(T(grad(U))))) Gauss linear; } laplacianSchemes { default Gauss linear corrected; } interpolationSchemes { default linear; } snGradSchemes { default corrected; } wallDist { method meshWave; } // ************************************************************************* //
-
| ========= | | | \\ / 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; smoother GaussSeidel; cacheAgglomeration true; nCellsInCoarsestLevel 10; agglomerator faceAreaPair; mergeLevels 1; tolerance 1e-06; relTol 0.05; } pFinal { $p; tolerance 1e-06; relTol 0; } "(U|k|epsilon)" { solver smoothSolver; smoother symGaussSeidel; tolerance 1e-05; relTol 0.1; } "(U|k|epsilon)Final" { $U; tolerance 1e-05; relTol 0; } } PIMPLE { nOuterCorrectors 4; nCorrectors 1; nNonOrthogonalCorrectors 0; pRefCell 0; pRefValue 0; } SIMPLE { nNonOrthogonalCorrectors 0; residualControl { p 1e-2; U 1e-3; "(k|epsilon)" 1e-3; } } relaxationFactors { fields { p 0.3; } equations { U 0.7; k 0.7; "epsilon.*" 0.7; } } cache { grad(U); } // ************************************************************************* //
-
| ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 3.0.1 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ Build : 3.0.1-119cac7e8750 Exec : checkMesh Date : Mar 23 2017 Time : 23:51:13 Host : "ubuntu" PID : 41427 Case : /home/cc/OpenFOAM/cc-3.0.1/run/tutorials/incompressible/simpleFoam/nasa65 nProcs : 1 sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). fileModificationChecking : Monitoring run-time modified files using timeStampMaster allowSystemOperations : Allowing user-supplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create polyMesh for time = 0 Time = 0 Mesh stats points: 1876500 faces: 5504237 internal faces: 5380555 cells: 1814132 faces per cell: 6 boundary patches: 8 point zones: 0 face zones: 0 cell zones: 0 Overall number of cells of each type: hexahedra: 1814132 prisms: 0 wedges: 0 pyramids: 0 tet wedges: 0 tetrahedra: 0 polyhedra: 0 Checking topology... Boundary definition OK. Cell to face addressing OK. Point usage OK. Upper triangular ordering OK. Face vertices OK. Number of regions: 1 (OK). Checking patch topology for multiply connected surfaces... Patch Faces Points Surface topology hub 30748 31275 ok (non-closed singly connected) inlet 5546 5700 ok (non-closed singly connected) outlet 4543 4680 ok (non-closed singly connected) per1 15340 15660 ok (non-closed singly connected) per2 15340 15660 ok (non-closed singly connected) ps 10266 10500 ok (non-closed singly connected) ss 11151 11400 ok (non-closed singly connected) sym 30748 31275 ok (non-closed singly connected) Checking geometry... Overall domain bounding box (-0.0925 -0.0130406 -1.48492e-15) (0.13 0.0278965 0.02) Mesh has 3 geometric (non-empty/wedge) directions (1 1 1) Mesh has 3 solution (non-empty) directions (1 1 1) Boundary openness (-1.04226e-16 4.59207e-16 4.83911e-16) OK. ***High aspect ratio cells found, Max aspect ratio: 10084.9, number of cells 119355 <<Writing 119355 cells with high aspect ratio to set highAspectRatioCells Minimum face area = 9.85732e-14. Maximum face area = 3.49678e-06. Face area magnitudes OK. Min volume = 9.85915e-19. Max volume = 2.44324e-09. Total volume = 9.66559e-05. Cell volumes OK. Mesh non-orthogonality Max: 54.3818 average: 6.13782 Non-orthogonality check OK. Face pyramids OK. Max skewness = 0.752069 OK. Coupled point location match (average 0) OK. Failed 1 mesh checks.
-
@cfd-china 用的不可压,只是想尝试算个结果,曾经在fluent中算过,用的komega-sst,湍流模型,看曲线k和omega特别大,达到了10e35量级,2_1490340702077_omega 1_1490340702077_k 0_1490340702077_epsilon 错误信息0_1490340438293_PyFoamRunner.simpleFoam.logfile
-
@cfd-china 您好,我算的马赫数大约0.6,是一个平面的压气机叶栅,目前用的是simpleFoam,网格用pointwise画的,检查质量ok,但总是挂掉,K,Omega总是超级大。
-
我需要给定进口总温总压,进口气流速度角度,以及出口静压请问,这么给是否有问题,总是挂。1_1490324845490_U 0_1490324845490_p
请问openfoam中有转捩模型吗(transition model)
请问openfoam中有转捩模型吗(transition model)
请问openfoam中有转捩模型吗(transition model)
进口马赫数0.7左右,用rhosimpleFoam不收敛的问题
Commutation error一直再升是什么错误
进口马赫数0.7左右,用rhosimpleFoam不收敛的问题
进口马赫数0.7左右,用rhosimpleFoam不收敛的问题
进口马赫数0.7左右,用rhosimpleFoam不收敛的问题
进口马赫数0.7左右,用rhosimpleFoam不收敛的问题
进口马赫数0.7左右,用rhosimpleFoam不收敛的问题
进口马赫数0.7左右,用rhosimpleFoam不收敛的问题
进口马赫数0.7左右,用rhosimpleFoam不收敛的问题
进口马赫数0.7左右,用rhosimpleFoam不收敛的问题
进口马赫数0.7左右,用rhosimpleFoam不收敛的问题
进口马赫数0.7左右,用rhosimpleFoam不收敛的问题
进口马赫数0.7左右,用rhosimpleFoam不收敛的问题
进口马赫数0.7左右,用rhosimpleFoam不收敛的问题
我想用rhocentralfoamlts,稳态的求解器,但是编译的时候有错误提示,不知道什么原因。
我想用rhocentralfoamlts,稳态的求解器,但是编译的时候有错误提示,不知道什么原因。
请教关于进口条件的设置的问题
请教关于进口条件的设置的问题
如何得到processor 的编号?
请问openfoam中有算高亚音Ma=0.67的求解器吗,稳态的,用rhoSimpleFoam总是算发散。
如何得到processor 的编号?
openfoam初场插值的问题
平面叶栅计算,采用sst湍流模型,吸力面分离过小
平面叶栅计算,采用sst湍流模型,吸力面分离过小
平面叶栅计算,采用sst湍流模型,吸力面分离过小
计算平面叶栅,可压缩出现问题
请教关于进口条件的设置的问题
请教关于进口条件的设置的问题
请教关于进口条件的设置的问题
请教关于进口条件的设置的问题
请教关于进口条件的设置的问题
请教关于进口条件的设置的问题
请教关于进口条件的设置的问题
请教关于进口条件的设置的问题
请教关于进口条件的设置的问题
请教关于进口条件的设置的问题
请教关于进口条件的设置的问题
请教关于进口条件的设置的问题
请教关于进口条件的设置的问题
请教关于进口条件的设置的问题
请教关于进口条件的设置的问题
请教关于进口条件的设置的问题
请教关于进口条件的设置的问题
请教关于进口条件的设置的问题
请教关于进口条件的设置的问题