CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM (https://www.cfd-online.com/Forums/openfoam/)
-   -   Error using DPMFoam (https://www.cfd-online.com/Forums/openfoam/176327-error-using-dpmfoam.html)

raed August 16, 2016 06:00
Error using DPMFoam
 
Hey Foamers,

I am currently working on particle tracking using OpenFOAM and trying to compute my domain with DPMFoam. I am facing a ''couple'' of problems and would need some help.

The error message I am getting is the following:
Code:
--> FOAM FATAL ERROR:
    request for volVectorField U from objectRegistry region0 failed
    available objects of type volVectorField are
3
(
U.air
U.air_0
cloudVolSUSu
)
    From function objectRegistry::lookupObject<Type>(const word&) const
    in file /home/openfoam/OpenFOAM/OpenFOAM-3.0.1/src/OpenFOAM/lnInclude/objectRegistryTemplates.C at line 198.

FOAM aborting
I have read in previous thread that the problem could be a U field not defied as U.air somewhere. However I checked all of my U definitions. I am posting my boundary conditions here:

Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
|  \\    /  O peration    | Version:  3.0.1                                |
|  \\  /    A nd          | Web:      www.OpenFOAM.org                      |
|    \\/    M anipulation  |                                                |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version    2.0;
    format      ascii;
    class      volVectorField;
    location    "0";
    object      U.air;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 1 -1 0 0 0 0];
internalField    uniform (0 0 0);

boundaryField
{
    INLET
    {
        type            pressureInletOutletVelocity;
    phi        phi.air;
    inletValue  uniform (25 0 0);
    value          uniform (25 0 0);
    }
    OUTLET
    {
      //type            zeroGradient;
      type              pressureInletOutletVelocity;
      inletVelocity    uniform (0 0 0);
      value            uniform (0 0 0);
      phi        phi.air;
      alpha        alpha.air;
    }
    AUSSLASSROHR
    {
        type            fixedValue;
        value          uniform (0 0 0);
    }
    AUSSLASSGLOCKE
    {
        type            fixedValue;
        value          uniform (0 0 0);
    }
    EINLASSROHR
    {
        type            fixedValue;
        value          uniform (0 0 0);
    }
    MESSKAMMER
    {
        type            fixedValue;
        value          uniform (0 0 0);
    }
    DUESE_INNEN
    {
        type            fixedValue;
        value          uniform (0 0 0);
    }
    DUESE_STIRNFLAECHE
    {
        type            fixedValue;
        value          uniform (0 0 0);
    }
    DUESE_AUSSEN
    {
        type            fixedValue;
        value          uniform (0 0 0);
    }
}


// ************************************************************************* //
my p file:

Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
|  \\    /  O peration    | Version:  3.0.1                                |
|  \\  /    A nd          | Web:      www.OpenFOAM.org                      |
|    \\/    M anipulation  |                                                |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version    2.0;
    format      ascii;
    class      volScalarField;
    location    "0";
    object      p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 2 -2 0 0 0 0];

internalField uniform 120000;

boundaryField
{
    INLET
    {
        type            totalPressure;
        rho            rho.air;
    phi        phi.air
        psi            none;
        gamma          1.4;
        p0              uniform 147109;
    }
    OUTLET
      {
        type            fixedValue;
        value          uniform 96789.7;
    }
    AUSSLASSROHR
    {
        type            zeroGradient;
    }
    AUSSLASSGLOCKE 
   
    {
        type            zeroGradient;
    }
    EINLASSROHR
    {
        type            zeroGradient;
    }
    MESSKAMMER
    {
        type            zeroGradient;
    }
    DUESE_INNEN
    {
        type            zeroGradient;
    }
    DUESE_STIRNFLAECHE
    {
        type            zeroGradient;
    }
    DUESE_AUSSEN
    {
        type            zeroGradient;
    }
}


// ************************************************************************* //
and here is my kinematicCloudProperties File

Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
|  \\    /  O peration    | Version:  3.0.1                                |
|  \\  /    A nd          | Web:      www.OpenFOAM.org                      |
|    \\/    M anipulation  |                                                |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version    2.0;
    format      ascii;
    class      dictionary;
    location    "constant";
    object      particleProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solution
{
    active          true;
    coupled        true;
    transient      yes;
    cellValueSourceCorrection off;

    interpolationSchemes
    {
        rho.air        cell;
        U.air          cellPoint;
        mu.air          cell;
    }

    integrationSchemes
    {
        U              Euler;
    }

    sourceTerms
    {
        schemes
        {
            U semiImplicit 1;
        }
    }
}

constantProperties
{
    parcelTypeId 1;

    rhoMin          1e-15;
    minParcelMass  1e-15;

    rho0            2526;
    youngsModulus  1e8;
    poissonsRatio  0.35;

    constantVolume  false;

    alphaMax        0.99;
}

subModels
{
    particleForces
    {
        ErgunWenYuDrag
        {
            alphac alpha.air;
        }
        gravity;
    }


    injectionModels
    {
        model1
        {
            type            patchInjection;
            massTotal      0;
            parcelBasisType fixed;
            nParticle      1;
            SOI            0;
            parcelsPerSecond    100;
        duration        1;
        patchName        INLET;
            U0              (0 30 0);
        flowRateProfile constant 1;
            sizeDistribution
            {
                type        general;
                generalDistribution
                {
               
                    distribution
                    (
                        (05e-05        5.19202)
                        (10e-05        33.91157)
                        (15e-05        24.78618)
                        (20e-05        14.25689)
                        (25e-05        8.68968)
                        (30e-05        5.07100)
                        (35e-05        3.130547)
                        (40e-05        1.795223)
                        (45e-05        1.1497498)
                        (50e-05        0.6737130)
                        (55e-05        0.5002420)
                        (60e-05        0.2944973)
                        (65e-05        0.1855736)
                        (70e-05        0.10085525)
                        (75e-05        0.1129578828)
                        (80e-05        0.0726157818)

                    );
               
                }
            }
        }
    }

    dispersionModel none;

    patchInteractionModel localInteraction;

    localInteractionCoeffs
    {
        patches
        (
            INLET
            {
                type rebound;
                e    0.97;
                mu  0.09;
            }
            OUTLET
            {
                type rebound;
                e    0.97;
                mu  0.09;
            }
            AUSSLASSROHR
            {
                type rebound;
                e    0.97;
                mu  0.09;
            }
            AUSSLASSGLOCKE
            {
                type rebound;
                e    0.97;
                mu  0.09;
            }
            EINLASSROHR
            {
                type rebound;
                e    0.97;
                mu  0.09;
            }
            MESSKAMMER
            {
                type rebound;
                e    0.97;
                mu  0.09;
            }
            DUESE_INNEN
            {
                type rebound;
                e    0.97;
                mu  0.09;
            }
            DUESE_STIRNFLAECHE
            {
                type rebound;
                e    0.97;
                mu  0.09;
            }
            DUESE_AUSSEN
            {
                type rebound;
                e    0.97;
                mu  0.09;
            }
        );
    }

    StandardWallInteractionCoeffs
    {
        type rebound;
        e    0.97;
        mu  0.09;
    }

    heatTransferModel none;

    surfaceFilmModel none;

    collisionModel pairCollision;

    pairCollisionCoeffs
    {
        maxInteractionDistance  0.0025;

        writeReferredParticleCloud no;

        pairModel pairSpringSliderDashpot;

        pairSpringSliderDashpotCoeffs
        {
            useEquivalentSize  no;
            alpha              0.02;
            b                  1.5;
            mu                  0.10;
            cohesionEnergyDensity 0;
            collisionResolutionSteps 12;
        };

        wallModel wallSpringSliderDashpot;

        wallSpringSliderDashpotCoeffs
        {
            useEquivalentSize no;
            collisionResolutionSteps 12;
            youngsModulus  1e8;
            poissonsRatio  0.23;
            alpha          0.01;
            b              1.5;
            mu              0.09;
            cohesionEnergyDensity 0;
        };

        UName U.air;
    }

    stochasticCollisionModel none;

    radiation off;
}


cloudFunctions
{}


// ************************************************************************* //
the other files shouldn't have errors since they are almost the same as the ones in the tutorial case, which I could run successfully on the same computer with the same OpenFOAM version.

I would be extremely grateful if someone could help me and maybe report any previous similar errors with the same solver.

Thanks :)
Raed

Sugajen November 16, 2017 17:01
Hi Raed,

I am facing a similar issue. Were you able to figure out the solution for this?
My p_rgh and U.air are below.

p_rgh:
Code:
dimensions      [1 -1 -2 0 0 0 0];

internalField  uniform 1.04e5;

boundaryField
{

    inlet
    {
        type            totalPressure;
        rho            none;
        p0              $internalField;
        value          $internalField;
    }

    outlet
    {
        type                prghPressure;
        p                uniform 0;
    }

    top
    {
        type                prghPressure;
        p                uniform 1.37e5;
    }

    bottom
    {
        type                fixedFluxPressure;
        value            $internalField;
    }

    walls
    {
        type                fixedFluxPressure;
        value                $internalField;
    }

    frontAndBack
    {
        type            empty;
    }
}

U.air:
Code:
dimensions      [0 1 -1 0 0 0 0];

internalField  uniform (0.0 0.0 0.0);

boundaryField
{

    inlet
    {
        type            fixedValue;
        value                $internalField;
    }

    outlet
    {
        type            fixedValue;
        value                $internalField;
    }

    top
    {
        type            fixedValue;
        value                $internalField;
    }

    bottom
    {
        type            noSlip;
    }

    walls
    {
        type            noSlip;
    }

    frontAndBack
    {
        type            empty;
    }

}

Viet-Dung November 19, 2018 13:47
Hi,
Using DPMFoam to study the scour around a cylinder, I started with a classic model without particles in icoFoam. Everything is ok concerning the mesh, fluid behavior...
Then, in DPMFoam, with 1 particle and the same mesh, boundary conditions..., when I set g = 0, the simulation converges well and the results match those of icoFoam. But when g = -9.81, U_z, U_magnitude,Courant Number increase at the outlet and the simulation stops.
Did you solve the problem ?

Leon Nunes November 27, 2021 12:40
DPM Foam Problem
 
I’m trying to simulate a axyssimetric nozzle using wdge and extrudeMesh. The simulation is a multiphasic argon/metal particle. I used a dambreak case with particles as base. I’m getting this error:

"PIMPLE: iteration 1
DILUPBiCG: Solving for U.airx, Initial residual = 1, Final residual = 3.05784e-06, No Iterations 1
DILUPBiCG: Solving for U.airy, Initial residual = 1, Final residual = 1.28577e-06, No Iterations 1
DILUPBiCG: Solving for U.airz, Initial residual = 0.00526608, Final residual = 2.48302e-08, No Iterations 1

--> FOAM FATAL ERROR: (openfoam-2012)

request for volVectorField U from objectRegistry region0 failed
available objects of type volVectorField are
5(U.air U.air_0_0 U.air_0 HbyA cloudVolSUSu)

From const Type& Foam::objectRegistry::lookupObject(const Foam::word&, bool) const [with Type = Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>]
in file /home/pawan/OpenFOAM/OpenFOAM/OpenFOAM-v2012/src/OpenFOAM/lnInclude/objectRegistryTemplates.C at line 463.

FOAM aborting

#0 Foam::error::printStack(Foam::Ostream&)sh: 1: addr2line: not found
addr2line failed
#1 Foam::error::exitOrAbort(int, bool)sh: 1: addr2line: not found
addr2line failed
#2 Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const& Foam::objectRegistry::lookupObject<Foam::Geometric Field<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> >(Foam::word const&, bool) constsh: 1: addr2line: not found
addr2line failed
#3 Foam::totalPressureFvPatchScalarField::updateCoeff s()sh: 1: addr2line: not found
addr2line failed
#4 Foam::fvMatrix<double>::fvMatrix(Foam::GeometricFi eld<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensionSet const&)sh: 1: addr2line: not found
addr2line failed
#5 Foam::fv::gaussLaplacianScheme<double, Foam::SymmTensor<double> >::fvmLaplacianUncorrected(Foam::GeometricField<do uble, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&)sh: 1: addr2line: not found
addr2line failed
#6 Foam::fv::gaussLaplacianScheme<double, double>::fvmLaplacian(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&)sh: 1: addr2line: not found
addr2line failed
#7 Foam::tmp<Foam::fvMatrix<double> > Foam::fvm::laplacian<double, double>(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&)sh: 1: addr2line: not found
addr2line failed
#8 ?sh: 1: addr2line: not found
addr2line failed
#9 __libc_start_mainsh: 1: addr2line: not found
addr2line failed
#10 ?sh: 1: addr2line: not found"[/I][/I]



Thank you in advance,

Leon


All times are GMT -4. The time now is 05:54.