请教关于进口条件的设置的问题
-
| ========= | | | \\ / 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; } ) // ************************************************************************* // -
@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
