Setting up of denseParticleFoam

josepaulb · June 13, 2023, 07:00

Hello

I am currently validating a study on sediment transport problem in openfoam 10.I have used MPPICCloud in cloud properties.
This is the error i am receiving while running in denseParticleFoam as a trial.The geometry is cuboid with inlet,outlet,bed(wall) and symmetry boundaries.

Code:

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


PIMPLE: No convergence criteria found


PIMPLE: No corrector convergence criteria found
        Calculations will do 2 corrections


PIMPLE: Operating solver in transient mode with 2 outer correctors



Reading g
Reading field U

Reading field p

Reading/calculating continuous-phase face flux field phic

Creating turbulence model

Selecting viscosity model Newtonian
Creating field alphac

Constructing clouds
Selecting parcelCloud MPPICCloud
Constructing particle forces
    Selecting particle force sphereDrag
    Selecting particle force gravity
Constructing cloud functions
    none
Constructing particle injection models
Creating injector: model1
Selecting injection model patchFlowRateInjection
    Constructing 3-D injection
    Choosing nParticle to be a fixed value, massTotal variable now does not determine anything.
Selecting distribution model fixedValue
    Distribution min: 0.00017 max: 0.00017 mean: 0.00017
Selecting dispersion model none
Selecting patch interaction model localInteraction
    Interaction fields will not be written
Selecting stochastic collision model none
Selecting surface film model none
Selecting U integration scheme Euler
Selecting packing model implicit
Selecting particle stress model HarrisCrighton
Selecting damping model none
Selecting isotropy model stochastic
Selecting time scale model isotropic
Selecting turbulence model type RAS
Selecting RAS turbulence model kOmegaSST
Selecting patchDistMethod meshWave
RAS
{
    RASModel        kOmegaSST;
    turbulence      on;
    printCoeffs     on;
    alphaK1         0.85;
    alphaK2         1;
    alphaOmega1     0.5;
    alphaOmega2     0.856;
    gamma1          0.555555555556;
    gamma2          0.44;
    beta1           0.075;
    beta2           0.0828;
    betaStar        0.09;
    a1              0.31;
    b1              1;
    c1              10;
    F3              false;
}

No fvModels present
No fvConstraints present

Starting time loop

Courant Number mean: 2.58125000001e-05 max: 0.00413
deltaT = 0.000119904076739
Time = 0.000119904s


Solving 3-D cloud cloud


--> FOAM FATAL ERROR: 

    request for surfaceScalarField phi from objectRegistry region0 failed
    available objects of type surfaceScalarField are

3
(
phi.water
alphaPhi.water
alphacf
)


    From function const Type& Foam::objectRegistry::lookupObject(const Foam::word&) const [with Type = Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>]
    in file /home/ubuntu/OpenFOAM/OpenFOAM-10/src/OpenFOAM/lnInclude/objectRegistryTemplates.C at line 211.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::error::abort() at ??:?
#2  Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const& Foam::objectRegistry::lookupObject<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> >(Foam::word const&) const in "/opt/openfoam10/platforms/linux64GccDPInt32Opt/bin/denseParticleFoam"
#3  Foam::PatchFlowRateInjection<Foam::MomentumCloud<Foam::ParcelCloudBase<Foam::MPPICParcel<Foam::MomentumParcel<Foam::particle> > > > >::flowRate() const at ??:?
#4  Foam::PatchFlowRateInjection<Foam::MomentumCloud<Foam::ParcelCloudBase<Foam::MPPICParcel<Foam::MomentumParcel<Foam::particle> > > > >::parcelsToInject(double, double) at ??:?
#5  void Foam::InjectionModel<Foam::MomentumCloud<Foam::ParcelCloudBase<Foam::MPPICParcel<Foam::MomentumParcel<Foam::particle> > > > >::inject<Foam::MPPICCloud<Foam::MomentumCloud<Foam::ParcelCloudBase<Foam::MPPICParcel<Foam::MomentumParcel<Foam::particle> > > > > >(Foam::MPPICCloud<Foam::MomentumCloud<Foam::ParcelCloudBase<Foam::MPPICParcel<Foam::MomentumParcel<Foam::particle> > > > >&, Foam::MPPICParcel<Foam::MomentumParcel<Foam::particle> >::trackingData&) at ??:?
#6  void Foam::MomentumCloud<Foam::ParcelCloudBase<Foam::MPPICParcel<Foam::MomentumParcel<Foam::particle> > > >::evolveCloud<Foam::MPPICCloud<Foam::MomentumCloud<Foam::ParcelCloudBase<Foam::MPPICParcel<Foam::MomentumParcel<Foam::particle> > > > > >(Foam::MPPICCloud<Foam::MomentumCloud<Foam::ParcelCloudBase<Foam::MPPICParcel<Foam::MomentumParcel<Foam::particle> > > > >&, Foam::MPPICParcel<Foam::MomentumParcel<Foam::particle> >::trackingData&) at ??:?
#7  void Foam::MomentumCloud<Foam::ParcelCloudBase<Foam::MPPICParcel<Foam::MomentumParcel<Foam::particle> > > >::solve<Foam::MPPICCloud<Foam::MomentumCloud<Foam::ParcelCloudBase<Foam::MPPICParcel<Foam::MomentumParcel<Foam::particle> > > > > >(Foam::MPPICCloud<Foam::MomentumCloud<Foam::ParcelCloudBase<Foam::MPPICParcel<Foam::MomentumParcel<Foam::particle> > > > >&, Foam::MPPICParcel<Foam::MomentumParcel<Foam::particle> >::trackingData&) at ??:?
#8  virtual thunk to Foam::ParcelCloud<Foam::MPPICCloud<Foam::MomentumCloud<Foam::ParcelCloudBase<Foam::MPPICParcel<Foam::MomentumParcel<Foam::particle> > > > > >::evolve() at ??:?
#9  Foam::parcelCloudList::evolve() at ??:?
#10  ? in "/opt/openfoam10/platforms/linux64GccDPInt32Opt/bin/denseParticleFoam"
#11  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#12  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#13  ? in "/opt/openfoam10/platforms/linux64GccDPInt32Opt/bin/denseParticleFoam"
Aborted (core dumped)

I searched for previous posts regarding same problem and tried to change according to the recommendations.
I don't where else should I define phi.But still I am getting the error.

These are the boundary conditions used.I am combining into a single code block

.

Code:

/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |                
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  10                                       
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    format      ascii;
    class       volScalarField;
    object      k.water;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

// value computed for a turbulence intensity = 5 %    -    0.24
internalField   uniform 0.00063963375;

boundaryField
{
    side_1
    {
     type symmetry;
   //  phi     phi.water;
    }
    side_2
    {
      type symmetry;
     // phi     phi.water;
    }
    Top
    {
      type symmetry;
    //  phi     phi.water;
    }
    Bottom
    {
      type kqRWallFunction;
      
    }
    Inlet
    {
      type fixedValue;
      value uniform 0.00063963375; 
    //        phi phi.water;
    }
    Outlet
    {
       type            inletOutlet;
    phi             phi.water;
        inletValue      $internalField;
        value           $internalField;

    
    }
    
}

// ************************************************************************* 
/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |                
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  10                                       
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    format      ascii;
    class       volScalarField;
    object      nut.water;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform 0;

boundaryField
{
side_1
    {
     type symmetry;
    }
    side_2
    {
     type symmetry;    
    }
    Top
    {
     type symmetry;    
    }
    Bottom
    {
     type nutkRoughWallFunction;
     value uniform 0;
     Ks uniform 425e-6;
     Cs  uniform 0.5;
    }
    Inlet
    {
     type   fixedValue;
     value  uniform 0; 
    }
    Outlet
    {
      type zeroGradient;
    //  phi phi.water;
    
    }
}

// ************************************************************************* 
/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |                
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  10                                      
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    format      ascii;
    class       volScalarField;
    object      omega.water;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

//value computed for a turbulence intensity = 5%     -    345

internalField   uniform 18.46992149415;

boundaryField
{
   side_1
    {
     type symmetry;
    }
    side_2
    {
      type symmetry;
    }
    Top
    {
      type symmetry;
    }
    Bottom
    {
      type omegaWallFunction;
      value $internalField;
    }
    Inlet
    {
      type fixedValue;
      value uniform 18.46992149415;
    }
    Outlet
    {
     type            inletOutlet;
    phi             phi.water;
        inletValue      $internalField;
        value           $internalField;

    
    }
}

// ************************************************************************* 
/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |                
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  10                                       
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    format      ascii;
    class       volScalarField;
    object      p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform 0;

boundaryField
{
    side_1
    {
         type symmetry;
    //  phi     phi.water;
          }
    side_2
    {
         type symmetry;
   //   phi     phi.water;    
      }
    Top
    {
     type symmetry;
     // phi     phi.water;    
      }
    Bottom
    {
    type fixedFluxPressure;
    value $internalField;
//      phi     phi.water;    
      }
    Inlet
    {
    type totalPressure;
    p0 uniform 88.76115825;
    value uniform 88.76115825;
  //  phi phi.water;
    }
    Outlet
    {
     type            fixedFluxPressure;
        gradient        uniform 0;
        value           uniform 0;
        phi phi.water;
    }
    
}

// ************************************************************************* 
/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |                
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  10                                      
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    format      ascii;
    class       volVectorField;
    object      U.water;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
    side_1
    {
         type symmetry;
   //   phi     phi.water;
    }
    side_2
    {
         type symmetry;
   //            phi     phi.water;
    }
    Top
    {
     type symmetry;
  //    phi     phi.water;
    }
    Bottom
    {
     type noSlip;
//      phi     phi.water;
    }
    Inlet
    {
      type            fixedValue;
        value           uniform (0 0.413 0);
    }
    Outlet
    {
     type            pressureInletOutletVelocity;
        inletValue      uniform (0 0 0);
        value           uniform (0 0 0);
     
     phi phi.water;
     
    }
    
}

 // ************************************************************************* //

Could someone give any suggestions for the change I should make?Any suggestions are fine.

efalpha · November 8, 2023, 09:30

Hi,
It is a bit late, however i post my solution just in case. I had the same problem and i solved it via setting "phi phi.x;" where x is a placeholder for the name of your phase.

However, I see that you have a similar line per patch commented out. Did you try it?

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Fluent Parallelization Problem After AC Power Dropped	pawl	Hardware	5	November 13, 2016 06:08
using chemkin	JMDag2004	OpenFOAM Pre-Processing	2	March 8, 2016 22:38
[snappyHexMesh] determining displacement for added points	CFDnewbie147	OpenFOAM Meshing & Mesh Conversion	1	October 22, 2013 09:53
Cells with t below lower limit	Purushothama	Siemens	2	May 31, 2010 21:58
Warning 097-	AB	Siemens	6	November 15, 2004 04:41

November 8, 2023, 09:30	a possible solution	#2
efalpha New Member Fotis Anagnostopoulos Join Date: Feb 2023 Location: Athens, Greece Posts: 10 Rep Power: 3	Hi, It is a bit late, however i post my solution just in case. I had the same problem and i solved it via setting "phi phi.x;" where x is a placeholder for the name of your phase. However, I see that you have a similar line per patch commented out. Did you try it?