View Issue Details

IDProjectCategoryView StatusLast Update
0001325OpenFOAMBugpublic2014-06-11 11:36
Reportergwierink Assigned Towill  
PrioritynormalSeveritycrashReproducibilityalways
Status resolvedResolutionfixed 
Platformx86_64OSUbuntuOS Version12.04
Summary0001325: JohnsonJacksonParticleSlip cannot find gs0.particles
DescriptionJohnsonJacksonParticleSlip for $FOAM_TUTORIALS/multiphase/twoPhaseEulerFoam/RAS/fluidisedBed/0/U.particles cannot find gs0:

--> FOAM FATAL ERROR:

    request for volScalarField gs0.particles from objectRegistry region0 failed
    available objects of type volScalarField are

8
(
alpha.particles
thermo:psi.particles
thermo:rho.particles
thermo:mu.particles
thermo:alpha.particles
p
T.particles
e.particles
)
Steps To ReproduceApply JohnsonJacksonParticleSlip to $FOAM_TUTORIALS/multiphase/twoPhaseEulerFoam/RAS/fluidisedBed/0/U.particles and run twoPhaseEulerFoam.

    walls
    {
        // type fixedValue;
        // value uniform (0 0 0);
        //
        type JohnsonJacksonParticleSlip;
        specularityCoefficient 0.01;
    }
Additional InformationtwoPhaseEulerFoam also does not run properly on $FOAM_TUTORIALS/multiphase/twoPhaseEulerFoam/RAS/fluidisedBed without changes.
TagsNo tags attached.

Activities

gwierink

2014-06-11 11:21

reporter   ~0003128

Error when running twoPhaseEulerFoam on fluidisedBed:

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.x |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.3.x-fe4674e7ca68
Exec : twoPhaseEulerFoam
Date : Jun 11 2014
Time : 12:18:10
Host :
PID : 10080
Case : /home/gijsbert/OpenFOAM/OpenFOAM-2.3.x/tutorials/multiphase/twoPhaseEulerFoam/RAS/fluidisedBed
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
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;
    }
}

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 heatTransferModel for (particles in air): RanzMarshall
Calculating field DDtU1 and DDtU2

Creating field dpdt

Creating field kinetic energy K

No MRF models present

No finite volume options present

Courant Number mean: 0.00725 max: 0.01
Max Ur Courant Number = 0.01

PIMPLE: no residual control data found. Calculations will employ 3 corrector loops


Starting time loop

fieldAverage fieldAverage1:
    Starting averaging at time 0

Courant Number mean: 0.00725 max: 0.01
Max Ur Courant Number = 0.01
Time = 0.0002

PIMPLE: iteration 1
MULES: Solving for alpha.particles
MULES: Solving for alpha.particles
smoothSolver: Solving for alpha.particles, Initial residual = 4.69513e-17, Final residual = 1.54587e-17, No Iterations 1
alpha.particles volume fraction = 0.275008 Min(alpha1) = 0 Max(alpha1) = 0.55
smoothSolver: Solving for e.particles, Initial residual = 1, Final residual = 5.33467e-08, No Iterations 2
smoothSolver: Solving for e.air, Initial residual = 1, Final residual = 2.54982e-07, No Iterations 3
GAMG: Solving for p, Initial residual = 1, Final residual = 0.00687942, No Iterations 5
PIMPLE: iteration 2
MULES: Solving for alpha.particles
MULES: Solving for alpha.particles
smoothSolver: Solving for alpha.particles, Initial residual = 1.02469e-05, Final residual = 2.4109e-10, No Iterations 4
alpha.particles volume fraction = 0.275 Min(alpha1) = 0 Max(alpha1) = 0.550058
smoothSolver: Solving for e.particles, Initial residual = 0.981449, Final residual = 1.05203e-07, No Iterations 1
smoothSolver: Solving for e.air, Initial residual = 0.067867, Final residual = 4.94025e-07, No Iterations 3
GAMG: Solving for p, Initial residual = 0.199092, Final residual = 0.000992419, No Iterations 4
PIMPLE: iteration 3
MULES: Solving for alpha.particles
MULES: Solving for alpha.particles
smoothSolver: Solving for alpha.particles, Initial residual = 1.08562e-05, Final residual = 3.98266e-10, No Iterations 1
alpha.particles volume fraction = 0.275 Min(alpha1) = 0 Max(alpha1) = 0.55007
smoothSolver: Solving for e.particles, Initial residual = 0.0133306, Final residual = 5.25778e-07, No Iterations 1
smoothSolver: Solving for e.air, Initial residual = 0.0803781, Final residual = 5.18721e-07, No Iterations 3
GAMG: Solving for p, Initial residual = 0.103434, Final residual = 4.11953e-09, No Iterations 15
word::stripInvalid() called for word laplacian(kappa,Theta)
    For debug level (= 2) > 1 this is considered fatal
Aborted (core dumped)

will

2014-06-11 11:36

manager   ~0003129

Thanks for the report(s).

Fixed in the following commits:
c70020edaebc0d66ebea82b55ca4786d4c05f49c
ddafa9836ad9fdbb12c9ee18916f55d419596a24

You will need to specify a "value" entry for the JohnsonJackson BC-s now.

Issue History

Date Modified Username Field Change
2014-06-11 10:55 gwierink New Issue
2014-06-11 11:21 gwierink Note Added: 0003128
2014-06-11 11:36 will Note Added: 0003129
2014-06-11 11:36 will Status new => resolved
2014-06-11 11:36 will Resolution open => fixed
2014-06-11 11:36 will Assigned To => will