|
[Sponsors] |
March 11, 2015, 10:40 |
twoPhaseEulerFoam execution error
|
#1 |
New Member
Tomas Trewhela
Join Date: Mar 2015
Posts: 1
Rep Power: 0 |
Hi,
I'm new in OpenFoam and I'm trying to model a snowdrift in a rectangular grid, with an inlet velocity from the left face and outlet at the right face of the grid. I'm getting this error and there are a few things I can't explain, such as the negative courant number and a negative volume fraction for alpha.particles. ¿Could this be a BC error? Thanks, /*---------------------------------------------------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.3.0 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ Build : 2.3.0-f5222ca19ce6 Exec : twoPhaseEulerFoam Date : Mar 11 2015 Time : 11:27:51 Host : "TT-est-rhma" PID : 32230 Case : /home/tomas/OpenFOAM/tomas-2.3.0/run/snowdrift nProcs : 1 sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). fileModificationChecking : Monitoring run-time modified files using timeStampMaster allowSystemOperations : Disallowing user-supplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create mesh for time = 0 Reading g 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 on; e 0.8; alphaMax 0.62; alphaMinFriction 0.5; residualAlpha 1e-06; viscosityModel Gidaspow; conductivityModel Gidaspow; granularPressureModel Lun; frictionalStressModel JohnsonJackson; radialModel SinclairJackson; JohnsonJacksonCoeffs { Fr 0.05; eta 2; p 5; phi 28.5; } } 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; 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 virtualMassModel for (particles in air): constantCoefficient Selecting heatTransferModel for (particles in air): RanzMarshall Calculating field DDtU1 and DDtU2 Creating field dpdt Creating field kinetic energy K No MRF models present Courant Number mean: -0.0797 max: -0.04 Max Ur Courant Number = -0.04 PIMPLE: Operating solver in PISO mode Starting time loop fieldAverage fieldAverage1: Starting averaging at time 0 Courant Number mean: -0.0797 max: -0.04 Max Ur Courant Number = -0.04 Time = 0.0002 MULES: Solving for alpha.particles MULES: Solving for alpha.particles smoothSolver: Solving for alpha.particles, Initial residual = 0, Final residual = 0, No Iterations 1 alpha.particles volume fraction = -0 Min(alpha1) = 0 Max(alpha1) = 0.5 smoothSolver: Solving for e.particles, Initial residual = 1, Final residual = 1.12571e-16, No Iterations 1 smoothSolver: Solving for e.air, Initial residual = 1, Final residual = 0.994753, No Iterations 1000 #0 Foam::error:rintStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 in "/lib/x86_64-linux-gnu/libc.so.6" #3 Foam::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const at ??:? #4 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMa trix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:? #5 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:? #6 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:? #7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:? #8 at ??:? #9 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #10 at ??:? |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[OpenFOAM.org] compile error in dynamicMesh and thermophysicalModels libraries | NickG | OpenFOAM Installation | 3 | December 30, 2019 01:21 |
[swak4Foam] GroovyBC the dynamic cousin of funkySetFields that lives on the suburb of the mesh | gschaider | OpenFOAM Community Contributions | 300 | October 29, 2014 19:00 |
checking the system setup and Qt version | vivek070176 | OpenFOAM Installation | 22 | June 1, 2010 13:34 |
POSDAT problem | piotka | STAR-CD | 4 | June 12, 2009 09:43 |
Compiling problems with hello worldC | fw407 | OpenFOAM Installation | 21 | January 6, 2008 18:38 |