OpenFoam使用rhoReactingFoam过程中出现的问题



  • diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
    smoothSolver:  Solving for Ux, Initial residual = 0.00441557, Final residual = 1.54512e-07, No Iterations 1
    smoothSolver:  Solving for Uy, Initial residual = 0.00624943, Final residual = 1.18395e-07, No Iterations 1
    DILUPBiCG:  Solving for H2, Initial residual = 0.00188285, Final residual = 7.43415e-08, No Iterations 1
    DILUPBiCG:  Solving for O2, Initial residual = 0.00188283, Final residual = 1.01514e-07, No Iterations 1
    DILUPBiCG:  Solving for H, Initial residual = 0, Final residual = 0, No Iterations 0
    DILUPBiCG:  Solving for OH, Initial residual = 0, Final residual = 0, No Iterations 0
    DILUPBiCG:  Solving for H2O, Initial residual = 0.00188283, Final residual = 1.01514e-07, No Iterations 1
    DILUPBiCG:  Solving for O, Initial residual = 0, Final residual = 0, No Iterations 0
    smoothSolver:  Solving for h, Initial residual = 0.00147451, Final residual = 5.51534e-08, No Iterations 1
    --> FOAM Warning : 
        From function Foam::scalar Foam::janafThermo<EquationOfState>::limit(Foam::scalar) const [with EquationOfState = Foam::perfectGas<Foam::specie>; Foam::scalar = double]
        in file /home/ubuntu/OpenFOAM/OpenFOAM-6/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 117
        attempt to use janafThermo<EquationOfState> out of temperature range 200 -> 3500;  T = 30.769
    --> FOAM Warning : 
        From function Foam::scalar Foam::janafThermo<EquationOfState>::limit(Foam::scalar) const [with EquationOfState = Foam::perfectGas<Foam::specie>; Foam::scalar = double]
        in file /home/ubuntu/OpenFOAM/OpenFOAM-6/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 117
        attempt to use janafThermo<EquationOfState> out of temperature range 200 -> 3500;  T = -214.865
    --> FOAM Warning : 
        From function Foam::scalar Foam::janafThermo<EquationOfState>::limit(Foam::scalar) const [with EquationOfState = Foam::perfectGas<Foam::specie>; Foam::scalar = double]
        in file /home/ubuntu/OpenFOAM/OpenFOAM-6/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 117
        attempt to use janafThermo<EquationOfState> out of temperature range 200 -> 3500;  T = 147.685
    --> FOAM Warning : 
        From function Foam::scalar Foam::janafThermo<EquationOfState>::limit(Foam::scalar) const [with EquationOfState = Foam::perfectGas<Foam::specie>; Foam::scalar = double]
        in file /home/ubuntu/OpenFOAM/OpenFOAM-6/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 117
        attempt to use janafThermo<EquationOfState> out of temperature range 200 -> 3500;  T = 147.687
    --> FOAM Warning : 
        From function Foam::scalar Foam::janafThermo<EquationOfState>::limit(Foam::scalar) const [with EquationOfState = Foam::perfectGas<Foam::specie>; Foam::scalar = double]
        in file /home/ubuntu/OpenFOAM/OpenFOAM-6/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 117
        attempt to use janafThermo<EquationOfState> out of temperature range 200 -> 3500;  T = -215.127
    --> FOAM Warning : 
        From function Foam::scalar Foam::janafThermo<EquationOfState>::limit(Foam::scalar) const [with EquationOfState = Foam::perfectGas<Foam::specie>; Foam::scalar = double]
        in file /home/ubuntu/OpenFOAM/OpenFOAM-6/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 117
        attempt to use janafThermo<EquationOfState> out of temperature range 200 -> 3500;  T = 28.5007
    min/max(T) = 200, 534.025
    smoothSolver:  Solving for p, Initial residual = 0.00340835, Final residual = 7.54003e-11, No Iterations 1
    smoothSolver:  Solving for p, Initial residual = 2.66055e-10, Final residual = 2.66055e-10, No Iterations 0
    smoothSolver:  Solving for p, Initial residual = 2.66055e-10, Final residual = 2.66055e-10, No Iterations 0
    diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
    time step continuity errors : sum local = 1.58233e-12, global = 8.75409e-14, cumulative = 6.66904e-12
    smoothSolver:  Solving for p, Initial residual = 2.46585e-05, Final residual = 5.53512e-13, No Iterations 1
    smoothSolver:  Solving for p, Initial residual = 2.76545e-12, Final residual = 2.76545e-12, No Iterations 0
    smoothSolver:  Solving for p, Initial residual = 2.76545e-12, Final residual = 2.76545e-12, No Iterations 0
    diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
    time step continuity errors : sum local = 1.65649e-14, global = 1.10285e-15, cumulative = 6.67015e-12
    #0  Foam::error::printStack(Foam::Ostream&) at ??:?
    #1  Foam::sigFpe::sigHandler(int) at ??:?
    #2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
    #3  log in "/lib/x86_64-linux-gnu/libm.so.6"
    #4  Foam::omegaWallFunctionFvPatchScalarField::calculate(Foam::turbulenceModel const&, Foam::List<double> const&, Foam::fvPatch const&, Foam::Field<double>&, Foam::Field<double>&) at ??:?
    #5  Foam::omegaWallFunctionFvPatchScalarField::calculateTurbulenceFields(Foam::turbulenceModel const&, Foam::Field<double>&, Foam::Field<double>&) at ??:?
    #6  Foam::omegaWallFunctionFvPatchScalarField::updateCoeffs() at ??:?
    #7  Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::Boundary::updateCoeffs() at ??:?
    #8  Foam::kOmegaSST<Foam::eddyViscosity<Foam::RASModel<Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foam::CompressibleTurbulenceModel<Foam::fluidThermo> > > > >, Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foam::CompressibleTurbulenceModel<Foam::fluidThermo> > > >::correct() at ??:?
    #9  ? at ??:?
    #10  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
    #11  ? at ??:?
    浮点数例外 (核心已转储)
    


  • @Rui 我自己尝试过的办法有检查边界,更改离散格式和时间步长等等,换过求解器为reactingFoam.但都未能解决这个问题,希望大佬们能够指点,谢谢~



  • 可压缩那面温度经常降为0K以下,你这个就是,原因不明,只能看看fvOptions能否帮助limit以下

    limitT
    {
        type            limitTemperature;
        active          yes;
        selectionMode   all;
        min             200;
        max             5000;
    }
    

Log in to reply
 

CFD中文网 2016 - 2020 | 京ICP备15017992号-2