请教关于进口条件的设置的问题
-
@cfd-china 用的不可压，只是想尝试算个结果，曾经在fluent中算过，用的komega-sst，湍流模型，看曲线k和omega特别大，达到了10e35量级，2_1490340702077_omega 1_1490340702077_k 0_1490340702077_epsilon 错误信息0_1490340438293_PyFoamRunner.simpleFoam.logfile
-
- 你这个网格哪里调用了交界面边界？抱歉看到了。
- 结果从第一步就没有保证连续性
- 将网格贴图（而非显示面）
- 贴
chechMesh结果 - 贴
fvSolution - 贴
fvScheme
注意是贴哦不是上传附近如：
Time = 2001 GAMG: Solving for Ux, Initial residual = 0.00638252, Final residual = 0.000400813, No Iterations 1 GAMG: Solving for Uy, Initial residual = 0.00114759, Final residual = 5.63775e-05, No Iterations 1 GAMG: Solving for Uz, Initial residual = 0.0283682, Final residual = 0.000925109, No Iterations 1 GAMG: Solving for p, Initial residual = 0.867854, Final residual = 0.0670769, No Iterations 4 time step continuity errors : sum local = 43.0388, global = -0.914839, cumulative = -0.914839 GAMG: Solving for omega, Initial residual = 0.0563256, Final residual = 0.00245155, No Iterations 1 GAMG: Solving for k, Initial residual = 1, Final residual = 0.0431817, No Iterations 1 ExecutionTime = 32.24 s ClockTime = 33 s Time = 2002 GAMG: Solving for Ux, Initial residual = 0.00830621, Final residual = 0.000239452, No Iterations 1 GAMG: Solving for Uy, Initial residual = 0.00165236, Final residual = 4.55896e-05, No Iterations 1 GAMG: Solving for Uz, Initial residual = 0.0188184, Final residual = 0.000526957, No Iterations 1 GAMG: Solving for p, Initial residual = 0.837434, Final residual = 0.0711894, No Iterations 2 time step continuity errors : sum local = 148.659, global = -0.294784, cumulative = -1.20962 GAMG: Solving for omega, Initial residual = 0.0143057, Final residual = 8.19713e-06, No Iterations 1 bounding omega, min: -2.58031e+07 max: 3.30734e+11 average: 4.70801e+09 GAMG: Solving for k, Initial residual = 0.11034, Final residual = 0.000207143, No Iterations 1 ExecutionTime = 39.79 s ClockTime = 40 s Time = 2003 GAMG: Solving for Ux, Initial residual = 0.00392685, Final residual = 7.6632e-05, No Iterations 1 GAMG: Solving for Uy, Initial residual = 0.00110286, Final residual = 1.82308e-05, No Iterations 1 GAMG: Solving for Uz, Initial residual = 0.0143816, Final residual = 0.000210414, No Iterations 1 GAMG: Solving for p, Initial residual = 0.791161, Final residual = 0.0685184, No Iterations 2 time step continuity errors : sum local = 111.167, global = -0.185086, cumulative = -1.39471 GAMG: Solving for omega, Initial residual = 1, Final residual = 0.0313884, No Iterations 3 bounding omega, min: -1.11841e+23 max: 5.89039e+26 average: 3.21879e+22 GAMG: Solving for k, Initial residual = 0.999599, Final residual = 4.49893e-17, No Iterations 1 ExecutionTime = 47.54 s ClockTime = 48 s Time = 2004 GAMG: Solving for Ux, Initial residual = 0.0028955, Final residual = 7.92784e-05, No Iterations 1 GAMG: Solving for Uy, Initial residual = 0.0009558, Final residual = 2.10947e-05, No Iterations 1 GAMG: Solving for Uz, Initial residual = 0.0218884, Final residual = 0.000373211, No Iterations 1 GAMG: Solving for p, Initial residual = 0.777876, Final residual = 0.065162, No Iterations 3 time step continuity errors : sum local = 78.0342, global = 0.432795, cumulative = -0.961914 GAMG: Solving for omega, Initial residual = 0.538686, Final residual = 3.00162e-09, No Iterations 1 bounding omega, min: -9.64586e+19 max: 1.32289e+42 average: 7.89521e+35 GAMG: Solving for k, Initial residual = 0.178897, Final residual = 1.16811e-20, No Iterations 1 ExecutionTime = 54.91 s ClockTime = 56 s Time = 2005另外，这是一个极为简单的几何，建议采用全六面体。
-
| ========= | | | \\ / 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. -
| ========= | | | \\ / 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 | | \*---------------------------------------------------------------------------*/ 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 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 您看有什么问题吗
| ========= | | | \\ / 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; 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; } } // ************************************************************************* // -
| ========= | | | \\ / 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; } } // ************************************************************************* // -
我注意到你的log 文件里面从第二个时间步长开始， specific dissipation rate (omega)就已经发散了(10^9)。
omega 文件里
FoamFile { version 2.0; format ascii; class volScalarField; object omega; } 你的 internalField uniform $omegaInlet = 630 这个神奇的数字是怎么算的？ 从你的log file 看出，你是从2000 s 开始算的，贴出2000s的初始条件或许能帮助解决问题？ 对于你的的进，出口边界条件不是特别理解。用进口速度边界条件，出口用压力边界条件不能满足你的需求吗？气动问题大家通常都是用这个方法。 -
@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 | | \*---------------------------------------------------------------------------*/ 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 -
@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; } ) // ************************************************************************* //
