CFD中文网

zhanghan

fvSolution中修改每一项的tolerance，我的case是bubbleColumnEvaporatingReacting，自带的tolerance是1e-20，我修改为1e-10后计算的非常快，而且计算结果几乎相同。

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
    "alpha.*"
    {
        nAlphaCorr      1;
        nAlphaSubCycles 3;
    }

    p_rgh
    {
        solver          GAMG;
        smoother        DIC;
        tolerance       1e-10;
        relTol          0.001;
    }

    p_rghFinal
    {
        $p_rgh;
        relTol          0;
    }

    "U.*"
    {
        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance       1e-10;
        relTol          0;
        minIter         1;
    }

    "(e|h).*"
    {
        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance       1e-10;
        relTol          0;
        minIter         1;
    }

    "(k|epsilon|Theta).*"
    {
        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance       1e-10;
        relTol          0;
        minIter         1;
    }

    Yi
    {
        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance       1e-10;
        relTol          0;
        minIter         1;
        residualAlpha   1e-8;
    }
}

PIMPLE
{
    nOuterCorrectors 1;
    nCorrectors      1;
    nNonOrthogonalCorrectors 0;
}

relaxationFactors
{
    equations
    {
        ".*"            0.3;
    }
}


// ************************************************************************* //

zhanghan

@李东岳我直接在左边的框里打字，提交之后就是现在帖子。没有其他操作了。东岳老师，您知道我的关于reactingTwoPhaseEulerFoam中速度的残差residual问题问题该怎么解决么？

zhanghan

@bestucan 特别感谢您的回答，这个帖子我也看过但是不知道该怎么办，怎么操作？您能指点指点么？我用的是OpenFOAM6.0reactingTwoPhaseEulerFoam
求解器。以下是UEqn.H的代码

Info<< "Constructing momentum equations" << endl;

fvVectorMatrix U1Eqn(U1, rho1.dimensions()*U1.dimensions()*dimVol/dimTime);
fvVectorMatrix U2Eqn(U2, rho2.dimensions()*U2.dimensions()*dimVol/dimTime);

{
    autoPtr<phaseSystem::momentumTransferTable>
        momentumTransferPtr(fluid.momentumTransfer());

    phaseSystem::momentumTransferTable&
        momentumTransfer(momentumTransferPtr());

    {
        U1Eqn =
        (
            phase1.UEqn()
         ==
           *momentumTransfer[phase1.name()]
          + fvOptions(alpha1, rho1, U1)
        );
        U1Eqn.relax();
        fvOptions.constrain(U1Eqn);
        fvOptions.correct(U1);
    }

    {
        U2Eqn =
        (
            phase2.UEqn()
         ==
           *momentumTransfer[phase2.name()]
          + fvOptions(alpha2, rho2, U2)
        );
        U2Eqn.relax();
        fvOptions.constrain(U2Eqn);
        fvOptions.correct(U2);
    }
}

zhanghan

有哪位大神知道为什么用reactingTwoPhaseEulerFoam求解器计算时不输出气相和液相速度的残差吗？
该求解器使用PIMPLE算法，我查看了pimpleFoam求解器里的U.Eqn，在速度矩后面有以下代码

if (pimple.momentumPredictor())
{
    solve(UEqn == -fvc::grad(p));

    fvOptions.correct(U);
}

reactingTwoPhaseEulerFoam求解器的U.Eqn中无该项，是这个原因吗？如果没solve，为什么会有速度场呢？
我该怎么做才能让reactingTwoPhaseEulerFoam求解器计算的算例显示速度残差？
在线求指点

zhanghan

OpenFOAM中有pbm群体平衡模型的求解器吗？？？求大神指点

zhanghan

OpenFOAM中有pbm群体平衡模型的求解器吗？？？求指点

zhanghan

@李东岳
为什么呢？

zhanghan

感谢大师

zhanghan

我用的是六面体，求解器是buoyantBoussinesSimpleFoam

zhanghan

Time = 150

DILUPBiCG:  Solving for Ux, Initial residual = 0.00826367, Final residual = 0.000199357, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.00351916, Final residual = 7.27402e-05, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.00281649, Final residual = 5.96235e-05, No Iterations 1
DILUPBiCG:  Solving for T, Initial residual = 0.00370137, Final residual = 0.000199429, No Iterations 1
DICPCG:  Solving for p_rgh, Initial residual = 0.0258986, Final residual = 0.000255685, No Iterations 25
time step continuity errors : sum local = 9.32347e-06, global = 1.95768e-07, cumulative = -1.43517e-07
DILUPBiCG:  Solving for epsilon, Initial residual = 0.00141499, Final residual = 2.88278e-05, No Iterations 2
bounding epsilon, min: -1.3845e-05 max: 0.0318741 average: 0.00019363
DILUPBiCG:  Solving for k, Initial residual = 0.0035897, Final residual = 8.74542e-05, No Iterations 2
ExecutionTime = 49.64 s  ClockTime = 50 s

Time = 151

DILUPBiCG:  Solving for Ux, Initial residual = 0.00806578, Final residual = 0.00019513, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.00348652, Final residual = 7.20562e-05, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.00276889, Final residual = 5.93574e-05, No Iterations 1
DILUPBiCG:  Solving for T, Initial residual = 0.00372472, Final residual = 0.000206805, No Iterations 1
DICPCG:  Solving for p_rgh, Initial residual = 0.0256067, Final residual = 0.000252871, No Iterations 209
time step continuity errors : sum local = 9.21069e-06, global = -1.68826e-08, cumulative = -1.60399e-07
DILUPBiCG:  Solving for epsilon, Initial residual = 0.00138175, Final residual = 2.92088e-05, No Iterations 2
bounding epsilon, min: -9.34359e-06 max: 0.0318873 average: 0.000193974
DILUPBiCG:  Solving for k, Initial residual = 0.00354032, Final residual = 8.86824e-05, No Iterations 2
ExecutionTime = 50.04 s  ClockTime = 50 s

Time = 152

DILUPBiCG:  Solving for Ux, Initial residual = 0.007934, Final residual = 0.00019362, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.00340541, Final residual = 7.07251e-05, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.00275302, Final residual = 5.91439e-05, No Iterations 1
DILUPBiCG:  Solving for T, Initial residual = 0.0037123, Final residual = 0.000213689, No Iterations 1
DICPCG:  Solving for p_rgh, Initial residual = 0.0250098, Final residual = 0.000247275, No Iterations 23
time step continuity errors : sum local = 8.99532e-06, global = -3.24882e-08, cumulative = -1.92888e-07
DILUPBiCG:  Solving for epsilon, Initial residual = 0.00136729, Final residual = 3.04028e-05, No Iterations 2
bounding epsilon, min: -2.09042e-06 max: 0.0319003 average: 0.000194322
DILUPBiCG:  Solving for k, Initial residual = 0.00356288, Final residual = 9.68293e-05, No Iterations 2
ExecutionTime = 50.29 s  ClockTime = 50 s

Time = 153

DILUPBiCG:  Solving for Ux, Initial residual = 0.0078555, Final residual = 0.00019842, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.0034158, Final residual = 7.32936e-05, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.00274224, Final residual = 6.04341e-05, No Iterations 1
DILUPBiCG:  Solving for T, Initial residual = 0.00389528, Final residual = 0.000253005, No Iterations 1
DICPCG:  Solving for p_rgh, Initial residual = 0.024297, Final residual = 0.000241194, No Iterations 55
time step continuity errors : sum local = 8.74315e-06, global = -2.1097e-07, cumulative = -4.03858e-07
DILUPBiCG:  Solving for epsilon, Initial residual = 0.00133526, Final residual = 3.0388e-05, No Iterations 2
bounding epsilon, min: -4.67459e-07 max: 0.0319129 average: 0.000194669
DILUPBiCG:  Solving for k, Initial residual = 0.00350007, Final residual = 9.68988e-05, No Iterations 2
ExecutionTime = 50.57 s  ClockTime = 51 s

Time = 154

DILUPBiCG:  Solving for Ux, Initial residual = 0.00766995, Final residual = 0.000189005, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.00333889, Final residual = 7.02533e-05, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.00268644, Final residual = 5.87275e-05, No Iterations 1
DILUPBiCG:  Solving for T, Initial residual = 0.00367952, Final residual = 0.000212829, No Iterations 1
DICPCG:  Solving for p_rgh, Initial residual = 0.024061, Final residual = 0.00023743, No Iterations 24
time step continuity errors : sum local = 8.59228e-06, global = -3.679e-07, cumulative = -7.71758e-07
DILUPBiCG:  Solving for epsilon, Initial residual = 0.00131178, Final residual = 2.96102e-05, No Iterations 2
bounding epsilon, min: -1.24139e-05 max: 0.0319253 average: 0.000195017
DILUPBiCG:  Solving for k, Initial residual = 0.00343625, Final residual = 9.3087e-05, No Iterations 2
ExecutionTime = 50.81 s  ClockTime = 51 s

Time = 155

DILUPBiCG:  Solving for Ux, Initial residual = 0.0075257, Final residual = 0.000185498, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.00330689, Final residual = 6.95719e-05, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.00263843, Final residual = 5.68287e-05, No Iterations 1
DILUPBiCG:  Solving for T, Initial residual = 0.00362838, Final residual = 0.000209773, No Iterations 1
DICPCG:  Solving for p_rgh, Initial residual = 0.0236067, Final residual = 0.000219917, No Iterations 206
time step continuity errors : sum local = 7.9418e-06, global = 3.02227e-08, cumulative = -7.41536e-07
DILUPBiCG:  Solving for epsilon, Initial residual = 0.00130737, Final residual = 2.83841e-05, No Iterations 2
bounding epsilon, min: -4.76185e-06 max: 0.0319373 average: 0.000195346
DILUPBiCG:  Solving for k, Initial residual = 0.00339847, Final residual = 9.30199e-05, No Iterations 2
ExecutionTime = 51.22 s  ClockTime = 51 s

Time = 156

DILUPBiCG:  Solving for Ux, Initial residual = 0.00760435, Final residual = 0.000185143, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.00334821, Final residual = 6.95993e-05, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.00265672, Final residual = 5.67016e-05, No Iterations 1
DILUPBiCG:  Solving for T, Initial residual = 0.00358951, Final residual = 0.000190184, No Iterations 1
DICPCG:  Solving for p_rgh, Initial residual = 0.023508, Final residual = 0.000234265, No Iterations 26
time step continuity errors : sum local = 8.4461e-06, global = 2.18412e-07, cumulative = -5.23124e-07
DILUPBiCG:  Solving for epsilon, Initial residual = 0.00131108, Final residual = 2.78022e-05, No Iterations 2
bounding epsilon, min: -1.64757e-06 max: 0.0319493 average: 0.000195658
DILUPBiCG:  Solving for k, Initial residual = 0.0033769, Final residual = 8.83642e-05, No Iterations 2
ExecutionTime = 51.46 s  ClockTime = 51 s

Time = 157

DILUPBiCG:  Solving for Ux, Initial residual = 0.00748129, Final residual = 0.00018526, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.00333077, Final residual = 6.95565e-05, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.00262409, Final residual = 5.65079e-05, No Iterations 1
DILUPBiCG:  Solving for T, Initial residual = 0.00361847, Final residual = 0.000204245, No Iterations 1
DICPCG:  Solving for p_rgh, Initial residual = 0.0230253, Final residual = 0.000228085, No Iterations 226
time step continuity errors : sum local = 8.19885e-06, global = 5.21378e-09, cumulative = -5.1791e-07
DILUPBiCG:  Solving for epsilon, Initial residual = 0.00127826, Final residual = 2.81699e-05, No Iterations 2
bounding epsilon, min: -1.27848e-05 max: 0.0319612 average: 0.000195978
DILUPBiCG:  Solving for k, Initial residual = 0.00335192, Final residual = 9.93353e-05, No Iterations 2
ExecutionTime = 51.88 s  ClockTime = 52 s

Time = 158

DILUPBiCG:  Solving for Ux, Initial residual = 0.00751949, Final residual = 0.000183825, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.003262, Final residual = 6.87015e-05, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.00261507, Final residual = 5.61683e-05, No Iterations 1
DILUPBiCG:  Solving for T, Initial residual = 0.00363588, Final residual = 0.000197698, No Iterations 1
DICPCG:  Solving for p_rgh, Initial residual = 0.0224584, Final residual = 0.000220381, No Iterations 22
time step continuity errors : sum local = 7.90328e-06, global = -2.23323e-08, cumulative = -5.40243e-07
DILUPBiCG:  Solving for epsilon, Initial residual = 0.00125523, Final residual = 2.86238e-05, No Iterations 2
bounding epsilon, min: -3.54997e-06 max: 0.0319729 average: 0.000196303
DILUPBiCG:  Solving for k, Initial residual = 0.00330904, Final residual = 9.39092e-05, No Iterations 2
ExecutionTime = 52.12 s  ClockTime = 52 s

Time = 159

DILUPBiCG:  Solving for Ux, Initial residual = 0.00739305, Final residual = 0.000181325, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.00325567, Final residual = 6.93665e-05, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.0026006, Final residual = 5.6702e-05, No Iterations 1
DILUPBiCG:  Solving for T, Initial residual = 0.00373395, Final residual = 0.000222137, No Iterations 1
DICPCG:  Solving for p_rgh, Initial residual = 0.0220611, Final residual = 0.000209639, No Iterations 54
time step continuity errors : sum local = 7.4913e-06, global = -1.98319e-07, cumulative = -7.38562e-07
DILUPBiCG:  Solving for epsilon, Initial residual = 0.00123809, Final residual = 2.97651e-05, No Iterations 2
bounding epsilon, min: -8.03126e-07 max: 0.0319843 average: 0.00019663
DILUPBiCG:  Solving for k, Initial residual = 0.00326329, Final residual = 9.44291e-05, No Iterations 2
ExecutionTime = 52.39 s  ClockTime = 52 s

Time = 160

DILUPBiCG:  Solving for Ux, Initial residual = 0.00707382, Final residual = 0.000174981, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.00322515, Final residual = 7.02473e-05, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.00256808, Final residual = 5.83962e-05, No Iterations 1
DILUPBiCG:  Solving for T, Initial residual = 0.00378956, Final residual = 0.000240966, No Iterations 1
DICPCG:  Solving for p_rgh, Initial residual = 0.0220337, Final residual = 0.000218272, No Iterations 25
time step continuity errors : sum local = 7.77203e-06, global = -3.57213e-07, cumulative = -1.09577e-06
DILUPBiCG:  Solving for epsilon, Initial residual = 0.00120289, Final residual = 2.9362e-05, No Iterations 2
bounding epsilon, min: -1.94654e-07 max: 0.0319956 average: 0.000196958
DILUPBiCG:  Solving for k, Initial residual = 0.00326015, Final residual = 0.000107391, No Iterations 2
ExecutionTime = 52.64 s  ClockTime = 53 s

Time = 161

DILUPBiCG:  Solving for Ux, Initial residual = 0.00694157, Final residual = 0.000173456, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.00322068, Final residual = 6.96394e-05, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.00252925, Final residual = 5.73262e-05, No Iterations 1
DILUPBiCG:  Solving for T, Initial residual = 0.0037208, Final residual = 0.000230451, No Iterations 1
DICPCG:  Solving for p_rgh, Initial residual = 0.022425, Final residual = 0.000224014, No Iterations 226
time step continuity errors : sum local = 7.94352e-06, global = 1.60691e-08, cumulative = -1.07971e-06
DILUPBiCG:  Solving for epsilon, Initial residual = 0.00119974, Final residual = 2.7392e-05, No Iterations 2
bounding epsilon, min: -1.14276e-05 max: 0.0320066 average: 0.000197268
DILUPBiCG:  Solving for k, Initial residual = 0.00318392, Final residual = 0.00010371, No Iterations 2
ExecutionTime = 53.05 s  ClockTime = 53 s

Time = 162

DILUPBiCG:  Solving for Ux, Initial residual = 0.00383098, Final residual = 9.45489e-05, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.0025646, Final residual = 5.49044e-05, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.00173677, Final residual = 3.95809e-05, No Iterations 1
DILUPBiCG:  Solving for T, Initial residual = 0.00210807, Final residual = 0.000135087, No Iterations 1
DICPCG:  Solving for p_rgh, Initial residual = 0.0229007, Final residual = 0.000227679, No Iterations 23
time step continuity errors : sum local = 8.03969e-06, global = 2.00147e-07, cumulative = -8.79559e-07
DILUPBiCG:  Solving for epsilon, Initial residual = 0.000991875, Final residual = 2.18963e-05, No Iterations 2
bounding epsilon, min: -9.09326e-06 max: 0.0320174 average: 0.000197563
DILUPBiCG:  Solving for k, Initial residual = 0.00285986, Final residual = 0.000101956, No Iterations 2
ExecutionTime = 53.3 s  ClockTime = 53 s

Time = 163

DILUPBiCG:  Solving for Ux, Initial residual = 0.00731843, Final residual = 0.000205183, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.00325374, Final residual = 7.02807e-05, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.00257374, Final residual = 5.79608e-05, No Iterations 1
DILUPBiCG:  Solving for T, Initial residual = 0.0036288, Final residual = 0.000195479, No Iterations 1
DICPCG:  Solving for p_rgh, Initial residual = 0.0231494, Final residual = 0.00021041, No Iterations 230
time step continuity errors : sum local = 7.39353e-06, global = -1.54056e-08, cumulative = -8.94964e-07
DILUPBiCG:  Solving for epsilon, Initial residual = 0.00118534, Final residual = 2.67827e-05, No Iterations 2
bounding epsilon, min: -7.90048e-07 max: 0.032028 average: 0.000197865
DILUPBiCG:  Solving for k, Initial residual = 0.00314674, Final residual = 9.33877e-05, No Iterations 2
ExecutionTime = 53.71 s  ClockTime = 54 s

Time = 164

DILUPBiCG:  Solving for Ux, Initial residual = 0.00109829, Final residual = 5.09955e-05, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.00196787, Final residual = 0.000119689, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.00115028, Final residual = 2.72874e-05, No Iterations 1
DILUPBiCG:  Solving for T, Initial residual = 0.00324789, Final residual = 0.000112044, No Iterations 1
DICPCG:  Solving for p_rgh, Initial residual = 0.021766, Final residual = 0.00020783, No Iterations 23
time step continuity errors : sum local = 7.25895e-06, global = -4.02801e-08, cumulative = -9.35244e-07
DILUPBiCG:  Solving for epsilon, Initial residual = 0.0010497, Final residual = 5.82234e-05, No Iterations 2
bounding epsilon, min: -1.17141e-06 max: 0.0320386 average: 0.000198149
DILUPBiCG:  Solving for k, Initial residual = 0.00313566, Final residual = 0.000111555, No Iterations 2
ExecutionTime = 53.95 s  ClockTime = 54 s

Time = 165

DILUPBiCG:  Solving for Ux, Initial residual = 0.916116, Final residual = 0.0687375, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.628033, Final residual = 0.0490923, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.684389, Final residual = 0.0540322, No Iterations 1
DILUPBiCG:  Solving for T, Initial residual = 0.00363733, Final residual = 0.000197744, No Iterations 1
DICPCG:  Solving for p_rgh, Initial residual = 0.93311, Final residual = 0.00907299, No Iterations 26
time step continuity errors : sum local = 0.00554995, global = -2.15723e-07, cumulative = -1.15097e-06
DILUPBiCG:  Solving for epsilon, Initial residual = 0.975191, Final residual = 0.0108415, No Iterations 3
bounding epsilon, min: -5.60805e-09 max: 0.652271 average: 0.00047027
DILUPBiCG:  Solving for k, Initial residual = 0.999823, Final residual = 0.0110168, No Iterations 2
ExecutionTime = 54.21 s  ClockTime = 54 s

Time = 166

DILUPBiCG:  Solving for Ux, Initial residual = 0.254913, Final residual = 0.0107531, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.221955, Final residual = 0.00933945, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.255228, Final residual = 0.00779565, No Iterations 1
DILUPBiCG:  Solving for T, Initial residual = 0.00559863, Final residual = 8.41868e-05, No Iterations 2
DICPCG:  Solving for p_rgh, Initial residual = 0.797795, Final residual = 0.00767407, No Iterations 114
time step continuity errors : sum local = 0.00145442, global = -9.6189e-07, cumulative = -2.11286e-06
DILUPBiCG:  Solving for epsilon, Initial residual = 0.992882, Final residual = 0.0350133, No Iterations 3
bounding epsilon, min: -1.37942e-07 max: 980.413 average: 0.189155
DILUPBiCG:  Solving for k, Initial residual = 0.959965, Final residual = 0.0355944, No Iterations 3
ExecutionTime = 54.53 s  ClockTime = 55 s

Time = 167

DILUPBiCG:  Solving for Ux, Initial residual = 0.129662, Final residual = 0.00406521, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.133136, Final residual = 0.00674105, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.142888, Final residual = 0.00531125, No Iterations 1
DILUPBiCG:  Solving for T, Initial residual = 0.137353, Final residual = 0.00850811, No Iterations 1
DICPCG:  Solving for p_rgh, Initial residual = 0.686287, Final residual = 0.00662569, No Iterations 105
time step continuity errors : sum local = 0.00117898, global = -2.06758e-06, cumulative = -4.18044e-06
DILUPBiCG:  Solving for epsilon, Initial residual = 0.0396736, Final residual = 0.000498224, No Iterations 1
bounding epsilon, min: -15.0105 max: 10176.7 average: 0.511179
DILUPBiCG:  Solving for k, Initial residual = 0.29001, Final residual = 0.0215973, No Iterations 2
ExecutionTime = 54.84 s  ClockTime = 55 s

Time = 168

DILUPBiCG:  Solving for Ux, Initial residual = 0.109495, Final residual = 0.00390908, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.109165, Final residual = 0.00482594, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.106908, Final residual = 0.0038473, No Iterations 1
DILUPBiCG:  Solving for T, Initial residual = 0.159025, Final residual = 0.00980602, No Iterations 1
DICPCG:  Solving for p_rgh, Initial residual = 0.782107, Final residual = 0.0045842, No Iterations 4
time step continuity errors : sum local = 4.06227, global = -3.19262e-06, cumulative = -7.37306e-06
DILUPBiCG:  Solving for epsilon, Initial residual = 0.889542, Final residual = 0.046541, No Iterations 2
bounding epsilon, min: -21.7381 max: 3.1549e+08 average: 3290.74
DILUPBiCG:  Solving for k, Initial residual = 0.213535, Final residual = 0.0095986, No Iterations 2
ExecutionTime = 55.08 s  ClockTime = 55 s

Time = 169

DILUPBiCG:  Solving for Ux, Initial residual = 0.174216, Final residual = 0.0154308, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.247336, Final residual = 0.0161903, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.287834, Final residual = 0.0191366, No Iterations 1
DILUPBiCG:  Solving for T, Initial residual = 0.124848, Final residual = 0.00426001, No Iterations 1
DICPCG:  Solving for p_rgh, Initial residual = 0.876163, Final residual = 0.00848489, No Iterations 134
time step continuity errors : sum local = 2.73843, global = -0.00923847, cumulative = -0.00924584
DILUPBiCG:  Solving for epsilon, Initial residual = 0.782446, Final residual = 0.0339288, No Iterations 2
DILUPBiCG:  Solving for k, Initial residual = 0.796286, Final residual = 0.0700807, No Iterations 2
ExecutionTime = 55.41 s  ClockTime = 55 s

Time = 170

DILUPBiCG:  Solving for Ux, Initial residual = 0.144657, Final residual = 0.00724012, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.207885, Final residual = 0.00705765, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.231172, Final residual = 0.0081708, No Iterations 1
DILUPBiCG:  Solving for T, Initial residual = 0.0825186, Final residual = 0.00407425, No Iterations 1
DICPCG:  Solving for p_rgh, Initial residual = 0.878503, Final residual = 0.00878104, No Iterations 124
time step continuity errors : sum local = 2.82593, global = -0.0534922, cumulative = -0.062738
DILUPBiCG:  Solving for epsilon, Initial residual = 0.0625576, Final residual = 0.000684989, No Iterations 1
bounding epsilon, min: -41.3991 max: 3.85702e+11 average: 1.66281e+07
DILUPBiCG:  Solving for k, Initial residual = 0.366084, Final residual = 0.00800287, No Iterations 2
ExecutionTime = 55.74 s  ClockTime = 56 s

Time = 171

DILUPBiCG:  Solving for Ux, Initial residual = 0.398345, Final residual = 0.015244, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.507564, Final residual = 0.0411822, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.586626, Final residual = 0.0394722, No Iterations 1
DILUPBiCG:  Solving for T, Initial residual = 0.0499476, Final residual = 0.00218177, No Iterations 1
DICPCG:  Solving for p_rgh, Initial residual = 0.825913, Final residual = 0.00759051, No Iterations 15
time step continuity errors : sum local = 6504.32, global = -0.0451193, cumulative = -0.107857
DILUPBiCG:  Solving for epsilon, Initial residual = 0.995574, Final residual = 0.029691, No Iterations 2
bounding epsilon, min: -14151 max: 1.48707e+15 average: 7.1395e+10
DILUPBiCG:  Solving for k, Initial residual = 0.143321, Final residual = 0.0120984, No Iterations 1
bounding k, min: -1.27293 max: 2.7193e+11 average: 5.90916e+06
ExecutionTime = 55.97 s  ClockTime = 56 s

Time = 172

DILUPBiCG:  Solving for Ux, Initial residual = 0.21569, Final residual = 0.00362517, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.0930196, Final residual = 0.00287307, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.264727, Final residual = 0.00398071, No Iterations 1
DILUPBiCG:  Solving for T, Initial residual = 0.0401774, Final residual = 0.00127523, No Iterations 1
DICPCG:  Solving for p_rgh, Initial residual = 0.918266, Final residual = 0.00896664, No Iterations 106
time step continuity errors : sum local = 4248.06, global = -52.9898, cumulative = -53.0977
DILUPBiCG:  Solving for epsilon, Initial residual = 0.0469099, Final residual = 0.0019286, No Iterations 1
bounding epsilon, min: -1.05659e+12 max: 2.19523e+18 average: 4.89037e+13
DILUPBiCG:  Solving for k, Initial residual = 0.382871, Final residual = 0.0198646, No Iterations 2
ExecutionTime = 56.29 s  ClockTime = 56 s

Time = 173

DILUPBiCG:  Solving for Ux, Initial residual = 0.157742, Final residual = 0.00172639, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.0865542, Final residual = 0.00219012, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.0649434, Final residual = 0.00138597, No Iterations 1
DILUPBiCG:  Solving for T, Initial residual = 0.0488159, Final residual = 0.00152357, No Iterations 1
DICPCG:  Solving for p_rgh, Initial residual = 0.959374, Final residual = 0.00776998, No Iterations 16
time step continuity errors : sum local = 4.79311e+06, global = -54.0686, cumulative = -107.166
DILUPBiCG:  Solving for epsilon, Initial residual = 0.0341313, Final residual = 1.46307e-05, No Iterations 1
bounding epsilon, min: -1.24625e+18 max: 1.15866e+20 average: 2.15355e+15
DILUPBiCG:  Solving for k, Initial residual = 0.299962, Final residual = 0.00352435, No Iterations 2
bounding k, min: -2.1925e+06 max: 7.58498e+14 average: 2.16185e+10
ExecutionTime = 56.52 s  ClockTime = 57 s

Time = 174

DILUPBiCG:  Solving for Ux, Initial residual = 0.156213, Final residual = 0.00156379, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.233462, Final residual = 0.00243138, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.175802, Final residual = 0.00172358, No Iterations 1
DILUPBiCG:  Solving for T, Initial residual = 0.0239022, Final residual = 0.000850046, No Iterations 1
DICPCG:  Solving for p_rgh, Initial residual = 0.921057, Final residual = 0.00779221, No Iterations 8
time step continuity errors : sum local = 2.53692e+09, global = -43.9838, cumulative = -151.15
DILUPBiCG:  Solving for epsilon, Initial residual = 0.0246204, Final residual = 0.000719125, No Iterations 2
bounding epsilon, min: -2.64601e+16 max: 4.38454e+26 average: 5.73001e+21
DILUPBiCG:  Solving for k, Initial residual = 0.0219395, Final residual = 0.000897031, No Iterations 2
bounding k, min: -4.36139e+06 max: 2.2271e+19 average: 5.04632e+14
ExecutionTime = 56.75 s  ClockTime = 57 s

Time = 175

DILUPBiCG:  Solving for Ux, Initial residual = 0.0928123, Final residual = 0.000826536, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.267698, Final residual = 0.00191604, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.209699, Final residual = 0.00159421, No Iterations 1
DILUPBiCG:  Solving for T, Initial residual = 0.0338199, Final residual = 0.00127989, No Iterations 1
DICPCG:  Solving for p_rgh, Initial residual = 0.911825, Final residual = 0.00892103, No Iterations 13
time step continuity errors : sum local = 5.71972e+12, global = 6.16787e+08, cumulative = 6.16786e+08
DILUPBiCG:  Solving for epsilon, Initial residual = 0.0906286, Final residual = 3.67525e-05, No Iterations 6
bounding epsilon, min: -1.90126e+21 max: 2.49982e+34 average: 2.11881e+30
DILUPBiCG:  Solving for k, Initial residual = 3.03287e-06, Final residual = 3.03287e-06, No Iterations 0
ExecutionTime = 56.99 s  ClockTime = 57 s

Time = 176

DILUPBiCG:  Solving for Ux, Initial residual = 0.181862, Final residual = 0.00187564, No Iterations 6
DILUPBiCG:  Solving for Uy, Initial residual = 0.192574, Final residual = 0.00113445, No Iterations 5
DILUPBiCG:  Solving for Uz, Initial residual = 0.160795, Final residual = 0.01166, No Iterations 4
DILUPBiCG:  Solving for T, Initial residual = 0.0229856, Final residual = 0.00073597, No Iterations 3
DICPCG:  Solving for p_rgh, Initial residual = 5.6219e-05, Final residual = 9.57783e-05, No Iterations 1001
time step continuity errors : sum local = 2.16016e+20, global = 1.04027e+17, cumulative = 1.04027e+17
DILUPBiCG:  Solving for epsilon, Initial residual = 0.54281, Final residual = 0.0254158, No Iterations 6
DILUPBiCG:  Solving for k, Initial residual = 1.50684e-07, Final residual = 1.50684e-07, No Iterations 0
ExecutionTime = 58.08 s  ClockTime = 58 s

Time = 177

DILUPBiCG:  Solving for Ux, Initial residual = 0.224522, Final residual = 0.00350717, No Iterations 23
DILUPBiCG:  Solving for Uy, Initial residual = 0.269353, Final residual = 0.00369183, No Iterations 25
DILUPBiCG:  Solving for Uz, Initial residual = 0.229731, Final residual = 0.00893461, No Iterations 18
DILUPBiCG:  Solving for T, Initial residual = 0.103618, Final residual = 0.00722955, No Iterations 23
--------------------------------------------------------------------------
mpirun noticed that process rank 1 with PID 43896 on node IP-219-216-120-163 exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------

感谢大师的回复！！！！！

zhanghan

请问这是怎么回事呢？？？？

Time = 176

DILUPBiCG:  Solving for Ux, Initial residual = 0.181862, Final residual = 0.00187564, No Iterations 6
DILUPBiCG:  Solving for Uy, Initial residual = 0.192574, Final residual = 0.00113445, No Iterations 5
DILUPBiCG:  Solving for Uz, Initial residual = 0.160795, Final residual = 0.01166, No Iterations 4
DILUPBiCG:  Solving for T, Initial residual = 0.0229856, Final residual = 0.00073597, No Iterations 3
DICPCG:  Solving for p_rgh, Initial residual = 5.6219e-05, Final residual = 9.57783e-05, No Iterations 1001
**time step continuity errors : sum local = 2.16016e+20, global = 1.04027e+17, cumulative = 1.04027e+17******
DILUPBiCG:  Solving for epsilon, Initial residual = 0.54281, Final residual = 0.0254158, No Iterations 6
DILUPBiCG:  Solving for k, Initial residual = 1.50684e-07, Final residual = 1.50684e-07, No Iterations 0
ExecutionTime = 58.08 s  ClockTime = 58 s

Time = 177

DILUPBiCG:  Solving for Ux, Initial residual = 0.224522, Final residual = 0.00350717, No Iterations 23
DILUPBiCG:  Solving for Uy, Initial residual = 0.269353, Final residual = 0.00369183, No Iterations 25
DILUPBiCG:  Solving for Uz, Initial residual = 0.229731, Final residual = 0.00893461, No Iterations 18
DILUPBiCG:  Solving for T, Initial residual = 0.103618, Final residual = 0.00722955, No Iterations 23
[5] #0  Foam::error::printStack(Foam::Ostream&)[8] #0  Foam::error::printStack(Foam::Ostream&)[9] #0  Foam::error::printStack(Foam::Ostream&)[7] #0  Foam::error::printStack(Foam::Ostream&)[0] #0  Foam::error::printStack(Foam::Ostream&)[1] #0  Foam::error::printStack(Foam::Ostream&)[2] #0  Foam::error::printStack(Foam::Ostream&)[3] #0  Foam::error::printStack(Foam::Ostream&)[4] #0  Foam::error::printStack(Foam::Ostream&) at ??:?
[4] #1  Foam::sigFpe::sigHandler(int) at ??:?
[2] #1  Foam::sigFpe::sigHandler(int) at ??:?

zhanghan

@队长别开枪
您好！我试过您的meshInfoFoam了，非常棒的求解器。
我有一个疑问，请问我怎样才能将导出的faceArea与我的cellCenter一一对应上呢？？？

zhanghan

@赵一铭

能指点我一下吗，实在不知道该怎么做

zhanghan

请问有哪位大神提取过网格体积吗？？？？

zhanghan

@chixt
我在这里@了您。
http://www.cfd-china.com/topic/872/传热边界/6

zhanghan

@zhanghan 在传热边界中说：

@李东岳

是的!!!!!
但是计算过程中出现了一个问题如下,我不知道该怎么办了,麻烦您帮我看一下,万分感谢!!!

/---------------------------------------------------------------------------
| ========= | |
| \ / 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-18e358711b4e
Exec : tgbuoyantw
Date : Jun 07 2017
Time : 09:17:04
Host : "zh-QiTianM4500-N000"
PID : 4374
Case : /home/zh/Downloads/Desktop/gT/wallfg
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_rgh tolerance 1e-05
field U tolerance 1e-06
field T tolerance 1e-05
field "(k|epsilon)" tolerance 0.001

Reading thermophysical properties

Reading field T

--> FOAM FATAL IO ERROR:
none is not in enumeration:
4
(
directionalSolidThermo
fluidThermo
lookup
solidThermo
)

file: /home/zh/Downloads/Desktop/gT/wallfg/0/T.boundaryField.WALL.kappa at line 37.
From function NamedEnum<Enum, nEnum>::read(Istream&) const
in file /home/openfoam/OpenFOAM/OpenFOAM-3.0.1/src/OpenFOAM/lnInclude/NamedEnum.C at line 72.
FOAM exiting

@chixt

以上是我设置这种边界条件后出现的问题。

zhanghan

@chixt

亲问您是怎么应用的这种边界条件，能指点我一下吗？万分感谢！

zhanghan

@李东岳

是的!!!!!
但是计算过程中出现了一个问题如下,我不知道该怎么办了,麻烦您帮我看一下,万分感谢!!!

/---------------------------------------------------------------------------
| ========= | |
| \ / 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-18e358711b4e
Exec : tgbuoyantw
Date : Jun 07 2017
Time : 09:17:04
Host : "zh-QiTianM4500-N000"
PID : 4374
Case : /home/zh/Downloads/Desktop/gT/wallfg
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_rgh tolerance 1e-05
field U tolerance 1e-06
field T tolerance 1e-05
field "(k|epsilon)" tolerance 0.001

Reading thermophysical properties

Reading field T

--> FOAM FATAL IO ERROR:
none is not in enumeration:
4
(
directionalSolidThermo
fluidThermo
lookup
solidThermo
)

file: /home/zh/Downloads/Desktop/gT/wallfg/0/T.boundaryField.WALL.kappa at line 37.

From function NamedEnum<Enum, nEnum>::read(Istream&) const
in file /home/openfoam/OpenFOAM/OpenFOAM-3.0.1/src/OpenFOAM/lnInclude/NamedEnum.C at line 72.

FOAM exiting

zhanghan

我没有用过fluent，在CFX中是热通量，好像也可以说成热流密度，单位是W/m2

zhanghan

我在buoyantBoussinesqSimpleFoam下计算不可压缩流体温度场和流场，涉及到壁面散热边界条件。
请问各位我将怎样设置壁面散热这个边界条件呢。