用compressibleinterfoam好像能量方程有问题,经常出现低速和超音速都出现,求解决?



  • /*---------------------------------------------------------------------------*\
    | =========                 |                                                 |
    | \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
    |  \\    /   O peration     | Version:  2.3.1                                 |
    |   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
    |    \\/     M anipulation  |                                                 |
    \*---------------------------------------------------------------------------*/
    Build  : 2.3.1-262087cdf8db
    Exec   : PBCavitatingFoam
    Date   : Jan 04 2021
    Time   : 21:20:13
    Host   : "mgt"
    PID    : 292963
    Case   : /home/test/liubo/liubo/PBC
    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
    
    
    Reading g
    
    PIMPLE: max iterations = 20
        field "(U|p_rgh|p)" : relTol 0.01, tolerance 1e-06
    
    Reading field p_rgh
    
    Reading field U
    
    Reading/calculating face flux field phi
    
    Constructing phaseChangeTwoPhaseMixture
    
    Selecting phaseChange model SchnerrSauer
    Selecting thermodynamics package 
    {
        type            heRhoThermo;
        mixture         pureMixture;
        transport       const;
        thermo          hConst;
        equationOfState rhoConst;
        specie          specie;
        energy          sensibleInternalEnergy;
    }
    
    Selecting thermodynamics package 
    {
        type            heRhoThermo;
        mixture         pureMixture;
        transport       const;
        thermo          hConst;
        equationOfState rhoConst;
        specie          specie;
        energy          sensibleInternalEnergy;
    }
    
    Selecting thermodynamics package 
    {
        type            heRhoThermo;
        mixture         pureMixture;
        transport       const;
        thermo          hConst;
        equationOfState perfectGas;
        specie          specie;
        energy          sensibleInternalEnergy;
    }
    
    leave from phaseChangeTwoPhaseMixture Constructor
    
    Reading thermophysical properties
    
    Calculating field g.h
    
    Selecting turbulence model type RASModel
    Selecting RAS turbulence model kEpsilon
    kEpsilonCoeffs
    {
        Cmu             0.09;
        C1              1.44;
        C2              1.92;
        C3              -0.33;
        sigmak          1;
        sigmaEps        1.3;
        Prt             1;
    }
    
    Creating field kinetic energy K
    
    Courant Number mean: 0.02677338955 max: 0.2119866949
    No finite volume options present
    
    
    Starting time loop
    
    Courant Number mean: 0.02677338955 max: 0.2119866949
    Interface Courant Number mean: 0 max: 0
    ^o^ deltaT = 1e-05
    Time = 1e-05
    
    PIMPLE: iteration 1
    MULES: Solving for alpha.air
    MULES: Solving for alpha.vapour
    Liquid phase volume fraction = 0.9999999964
      Min(alpha.water) = 0  Max(alpha.water) = 1
      Min(alpha.vapour) = 0  Max(alpha.vapour) = 0
      Min(alpha.air) = 0  Max(alpha.air) = 1
    MULES: Solving for alpha.air
    MULES: Solving for alpha.vapour
    Liquid phase volume fraction = 0.9999999927
      Min(alpha.water) = 0  Max(alpha.water) = 1
      Min(alpha.vapour) = 0  Max(alpha.vapour) = 0
      Min(alpha.air) = 0  Max(alpha.air) = 1
    MULES: Solving for alpha.air
    MULES: Solving for alpha.vapour
    Liquid phase volume fraction = 0.999999989
      Min(alpha.water) = 0  Max(alpha.water) = 1
      Min(alpha.vapour) = 0  Max(alpha.vapour) = 0
      Min(alpha.air) = 0  Max(alpha.air) = 1
    MULES: Solving for alpha.air
    MULES: Solving for alpha.vapour
    Liquid phase volume fraction = 0.9999999853
      Min(alpha.water) = 0  Max(alpha.water) = 1
      Min(alpha.vapour) = 0  Max(alpha.vapour) = 0
      Min(alpha.air) = 0  Max(alpha.air) = 1
    diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
    smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 1.565205693e-09, No Iterations 5
    smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 1.423403593e-09, No Iterations 5
    smoothSolver:  Solving for Uz, Initial residual = 1, Final residual = 1.232277888e-09, No Iterations 5
    smoothSolver:  Solving for T, Initial residual = 0.9999963229, Final residual = 0.09144259942, No Iterations 6
    min(T) 292.9992852
    GAMG:  Solving for p_rgh, Initial residual = 1, Final residual = 6085871.068, No Iterations 50
    GAMG:  Solving for p_rgh, Initial residual = 2.170075185e-15, Final residual = 2.170075185e-15, No Iterations 0
    max(U) 4889963.244
    min(p_rgh) -849012.5866
    GAMG:  Solving for p_rgh, Initial residual = 2.661173357e-10, Final residual = 2.661173357e-10, No Iterations 0
    GAMG:  Solving for p_rgh, Initial residual = 2.661173357e-10, Final residual = 2.661173357e-10, No Iterations 0
    max(U) 1.142439092e+14
    min(p_rgh) -849012.5866
    GAMG:  Solving for p_rgh, Initial residual = 2.676702599e-10, Final residual = 2.676702599e-10, No Iterations 0
    GAMG:  Solving for p_rgh, Initial residual = 2.676702599e-10, Final residual = 2.676702599e-10, No Iterations 0
    max(U) 1.425091029e+14
    min(p_rgh) -849012.5866
    GAMG:  Solving for p_rgh, Initial residual = 2.68040637e-10, Final residual = 2.68040637e-10, No Iterations 0
    GAMG:  Solving for p_rgh, Initial residual = 2.68040637e-10, Final residual = 2.68040637e-10, No Iterations 0
    max(U) 1.489703193e+14
    min(p_rgh) -849012.5866
    GAMG:  Solving for p_rgh, Initial residual = 2.681313485e-10, Final residual = 2.681313485e-10, No Iterations 0
    GAMG:  Solving for p_rgh, Initial residual = 2.681313485e-10, Final residual = 2.681313485e-10, No Iterations 0
    max(U) 1.504634378e+14
    min(p_rgh) -849012.5866
    PIMPLE: iteration 2
    MULES: Solving for alpha.air
    MULES: Solving for alpha.vapour
    Liquid phase volume fraction = 100322.0152
      Min(alpha.water) = 0  Max(alpha.water) = 2.810378902e+11
      Min(alpha.vapour) = 0  Max(alpha.vapour) = 0
      Min(alpha.air) = -2.810378902e+11  Max(alpha.air) = 2.858846245e+11
    MULES: Solving for alpha.air
    MULES: Solving for alpha.vapour
    Liquid phase volume fraction = 3.059042308e+16
      Min(alpha.water) = 0  Max(alpha.water) = 8.173458136e+22
      Min(alpha.vapour) = 0  Max(alpha.vapour) = 0
      Min(alpha.air) = -8.173458136e+22  Max(alpha.air) = 6.422063805e+21
    MULES: Solving for alpha.air
    MULES: Solving for alpha.vapour
    Liquid phase volume fraction = 2.539802609e+27
      Min(alpha.water) = 0  Max(alpha.water) = 1.817814046e+34
      Min(alpha.vapour) = 0  Max(alpha.vapour) = 0
      Min(alpha.air) = -1.817814046e+34  Max(alpha.air) = 2.285753964e+34
    MULES: Solving for alpha.air
    MULES: Solving for alpha.vapour
    Liquid phase volume fraction = 1.050092687e+39
      Min(alpha.water) = 0  Max(alpha.water) = 6.532567778e+45
      Min(alpha.vapour) = 0  Max(alpha.vapour) = 0
      Min(alpha.air) = -6.532567778e+45  Max(alpha.air) = 5.600607821e+44
    diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
    smoothSolver:  Solving for Ux, Initial residual = 0.07551363127, Final residual = 4.375766351e-09, No Iterations 5
    smoothSolver:  Solving for Uy, Initial residual = 0.08313923201, Final residual = 8.360047038e-09, No Iterations 5
    smoothSolver:  Solving for Uz, Initial residual = 0.08809270175, Final residual = 8.908948216e-10, No Iterations 6
    smoothSolver:  Solving for T, Initial residual = 1, Final residual = 6.382659004e+119, No Iterations 1000
    
    
    --> FOAM FATAL ERROR: 
    Maximum number of iterations exceeded
    
        From function thermo<Thermo, Type>::T(scalar f, scalar T0, scalar (thermo<Thermo, Type>::*F)(const scalar) const, scalar (thermo<Thermo, Type>::*dFdT)(const scalar) const, scalar (thermo<Thermo, Type>::*limit)(const scalar) const) const
        in file /opt/ohpc/pub/apps/OpenFOAM/OpenFOAM-2.3.1/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 76.
    
    FOAM aborting
    
    #0  Foam::error::printStack(Foam::Ostream&) at ??:?
    #1  Foam::error::abort() at ??:?
    #2  Foam::species::thermo<Foam::hConstThermo<Foam::rhoConst<Foam::specie> >, Foam::sensibleInternalEnergy>::TEs(double, double, double) const at ??:?
    #3  Foam::heRhoThermo<Foam::rhoThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::rhoConst<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate() at ??:?
    #4  Foam::heRhoThermo<Foam::rhoThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::rhoConst<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct() at ??:?
    #5  Foam::threePhaseMixtureThermo::correct() at ??:?
    #6  
     at ??:?
    #7  __libc_start_main in "/lib64/libc.so.6"
    #8  
     at ??:?

  • 版主 副教授

    我猜是网格质量问题
    https://cfd-china.com/topic/4357/checkmesh中发现的几个网格问题如何检查修改/4?_=1610591734371

    参考这个帖子描述的,四两拨千斤,千斤拨四两:chigua:


Log in to reply
 


CFD中文网 | 东岳流体学术 | 东岳流体商业 | 吉ICP备20003622号-1