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是出问题的地方。请求各位，这是曳力计算过程中出现问题了吗？
