twoPhaseEulerFoam求解器求解固液两相流问题出错,请指教。



  • 终端提示出错信息如下:

    /*---------------------------------------------------------------------------*\
    | =========                 |                                                 |
    | \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
    |  \\    /   O peration     | Version:  4.1                                   |
    |   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
    |    \\/     M anipulation  |                                                 |
    \*---------------------------------------------------------------------------*/
    Build  : 4.1
    Exec   : twoPhaseEulerFoamHydrate
    Date   : Aug 29 2017
    Time   : 10:50:03
    Host   : "JK"
    PID    : 5392
    Case   : /home/jk/OpenFOAM/jk-4.1/run/paperCase/fluidisedBedModify
    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
    
    
    PIMPLE: no residual control data found. Calculations will employ 3 corrector loops
    
    
    Reading g
    
    Reading hRef
    Creating twoPhaseSystem
    
    Selecting thermodynamics package 
    {
        type            heRhoThermo;
        mixture         pureMixture;
        transport       const;
        thermo          hConst;
        equationOfState rhoConst;
        specie          specie;
        energy          sensibleInternalEnergy;
    }
    
    Calculating face flux field phi.particles
    Selecting diameterModel for phase particles: constant
    Selecting turbulence model type RAS
    Selecting RAS turbulence model kineticTheory
    Selecting viscosityModel Gidaspow
    Selecting conductivityModel Gidaspow
    Selecting radialModel SinclairJackson
    Selecting granularPressureModel Lun
    Selecting frictionalStressModel JohnsonJackson
    kineticTheoryCoeffs
    {
        equilibrium     off;
        e               0.8;
        alphaMax        0.62;
        alphaMinFriction 0.5;
        residualAlpha   0.0001;
        viscosityModel  Gidaspow;
        conductivityModel Gidaspow;
        granularPressureModel Lun;
        frictionalStressModel JohnsonJackson;
        radialModel     SinclairJackson;
        JohnsonJacksonCoeffs
        {
            Fr              0.05;
            eta             2;
            p               5;
            phi             28.5;
            alphaDeltaMin   0.05;
        }
    }
    
    Selecting thermodynamics package 
    {
        type            heRhoThermo;
        mixture         pureMixture;
        transport       const;
        thermo          hConst;
        equationOfState perfectGas;
        specie          specie;
        energy          sensibleInternalEnergy;
    }
    
    Calculating face flux field phi.air
    Selecting diameterModel for phase air: constant
    Selecting turbulence model type RAS
    Selecting RAS turbulence model kEpsilon
    kEpsilonCoeffs
    {
        Cmu             0.09;
        C1              1.44;
        C2              1.92;
        C3              -0.33;
        sigmak          1;
        sigmaEps        1.3;
    }
    
    Selecting default blending method: none
    Selecting dragModel for (particles in air): GidaspowErgunWenYu
    Selecting swarmCorrection for (particles in air): none
    Selecting swarmCorrection for (particles in air): none
    Selecting swarmCorrection for (particles in air): none
    Selecting heatTransferModel for (particles in air): RanzMarshall
    Calculating field g.h
    
    Reading field p_rgh
    
    Creating field dpdt
    
    Creating field kinetic energy K
    
    No MRF models present
    
    No finite volume options present
    
    Courant Number mean: 0.000180506 max: 0.227071
    Max Ur Courant Number = 0
    Calculating field DDtU1 and DDtU2
    
    
    Starting time loop
    
    fieldAverage fieldAverage1:
        Starting averaging at time 0
    
    Courant Number mean: 0.000180506 max: 0.227071
    Max Ur Courant Number = 0
    Time = 0.0001
    
    PIMPLE: iteration 1
    MULES: Solving for alpha.particles
    MULES: Solving for alpha.particles
    alpha.particles volume fraction = 0.3  Min(alpha.particles) = 0.3  Max(alpha.particles) = 0.3
    Constructing momentum equations
    GAMG:  Solving for p_rgh, Initial residual = 0.99745, Final residual = 9.5914e-09, No Iterations 53
    PIMPLE: iteration 2
    MULES: Solving for alpha.particles
    MULES: Solving for alpha.particles
    smoothSolver:  Solving for alpha.particles, Initial residual = 1.95039e-05, Final residual = 2.69854e-14, No Iterations 1
    alpha.particles volume fraction = 0.300103  Min(alpha.particles) = 0.276552  Max(alpha.particles) = 0.389144
    Constructing momentum equations
    GAMG:  Solving for p_rgh, Initial residual = 0.606125, Final residual = 9.92728e-09, No Iterations 37
    PIMPLE: iteration 3
    MULES: Solving for alpha.particles
    MULES: Solving for alpha.particles
    smoothSolver:  Solving for alpha.particles, Initial residual = 7.59494e-06, Final residual = 1.31385e-14, No Iterations 1
    alpha.particles volume fraction = 0.299443  Min(alpha.particles) = 0.00196107  Max(alpha.particles) = 0.357417
    Constructing momentum equations
    GAMG:  Solving for p_rgh, Initial residual = 0.993696, Final residual = 8.98953e-09, No Iterations 53
    smoothSolver:  Solving for Theta.particles, Initial residual = 1, Final residual = 3.78152e-17, No Iterations 1
    smoothSolver:  Solving for epsilon.air, Initial residual = 0.0282379, Final residual = 1.65533e-06, No Iterations 1
    smoothSolver:  Solving for k.air, Initial residual = 1, Final residual = 5.48024e-08, No Iterations 2
    ExecutionTime = 23.69 s
    
    
        Reading/initialising field U.particlesMean
        Reading/initialising field U.airMean
        Reading/initialising field alpha.particlesMean
        Reading/initialising field pMean
    fieldAverage fieldAverage1 write:
        Calculating averages
    
    Courant Number mean: 4.33164 max: 1456.71
    Max Ur Courant Number = 1666.61
    Time = 0.0002
    
    PIMPLE: iteration 1
    MULES: Solving for alpha.particles
    MULES: Solving for alpha.particles
    smoothSolver:  Solving for alpha.particles, Initial residual = 6.32659e-06, Final residual = 7.27769e-15, No Iterations 1
    alpha.particles volume fraction = 0.302681  Min(alpha.particles) = 0.0538449  Max(alpha.particles) = 3.36489
    Constructing momentum equations
    #0  Foam::error::printStack(Foam::Ostream&) at ??:?
    #1  Foam::sigFpe::sigHandler(int) at ??:?
    #2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
    #3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
    #4  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
    #5  Foam::dragModels::Ergun::CdRe() const at ??:?
    #6  Foam::dragModels::GidaspowErgunWenYu::CdRe() const at ??:?
    #7  Foam::dragModel::Ki() const at ??:?
    #8  Foam::dragModel::K() const at ??:?
    #9  Foam::BlendedInterfacialModel<Foam::dragModel>::K() const at ??:?
    #10  Foam::twoPhaseSystem::Kd() const at ??:?
    #11  ? at ??:?
    #12  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
    #13  ? at ??:?
    浮点数例外 (核心已转储)
    

    Ps:用of时间较短,架构理解不是很深。网上查找相关问题时解答时,有前人说#3和#5给出了此问题的根源,#3是由于分母为0导致的,#5是出问题的地方。请求各位,这是曳力计算过程中出现问题了吗?



  • 你这个divide表示除掉了0无法 计算,openfoam里面cdre已经进行过稳定处理了,你看看你的Wenyu模型里面的CdRe函数除的变量是不是在你计算的时候变成了0.



  • @李东岳 好的 谢谢东岳老师


Log in to reply
 

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