CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Pre-Processing

Setting up of denseParticleFoam

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 13, 2023, 07:00
Default Setting up of denseParticleFoam
  #1
New Member
 
Jose paul b
Join Date: Aug 2021
Posts: 1
Rep Power: 0
josepaulb is on a distinguished road

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.
josepaulb is offline   Reply With Quote

Old   November 8, 2023, 09:30
Default a possible solution
  #2
New Member
 
Fotis Anagnostopoulos
Join Date: Feb 2023
Location: Athens, Greece
Posts: 10
Rep Power: 3
efalpha is on a distinguished road
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?
efalpha is offline   Reply With Quote

Reply

Tags
openfoam, sediment transport

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


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


All times are GMT -4. The time now is 15:44.