[hiuluom]

Hi guys,

I also problem with kEpsilon model in DPMFoam OF2.2.x. I tried to add terms in scheme but it stills errors.

--> FOAM FATAL ERROR:

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

6
(
phic_0
((interpolate((((1|A(Uc))*cloudVolSUSu)|rhoc))&S)+(Dp*(g&S)))
alphaPhic
phic
Dp
alphacf
)

Does anybody give me a suggest to fix it?

Thanks

[hiuluom]

Hi all,

I solved the above problems by adding phic at inlet and outlet BC for k and epsilon

k

 inlet_pipe
    {
        type                    turbulentIntensityKineticEnergyInlet; 
        value                   $initialField; 
        phi             phic;
        alpha           alphac;
        intensity               0.03; 
    }

outlet
    {
        type               inletOutlet;
        phi             phic;
        alpha           alphac;
        inletValue         $internalField;
        value              $internalField;
    }

epsilon

inlet_pipe
    {
        type                    turbulentMixingLengthDissipationRateInlet; 
        value                   $initialField; 
        phi             phic;
        alpha           alphac;
        mixingLength            0.042; 
    }

outlet
    {
        type               inletOutlet;
        phi             phic;
        alpha           alphac;
        inletValue         $internalField;
        value              $internalField;
    }

But I get a new error with that:

Courant Number mean: 5.18008e-007 max: 0.0244965
Time = 2e-005

Evolving kinematicCloud

Solving 3-D cloud kinematicCloud
Cloud: kinematicCloud
    Current number of parcels       = 0
    Current mass in system          = 0
    Linear momentum                 = (0 0 0)
   |Linear momentum|                = 0
    Linear kinetic energy           = 0
    Rotational kinetic energy       = 0
    model1:
        number of parcels added     = 0
        mass introduced             = 0
    Parcels absorbed into film      = 0
    New film detached parcels       = 0
    Parcel fate (number, mass)      : patch inlet_pipe
      - escape                      = 0, 0
      - stick                       = 0, 0
    Parcel fate (number, mass)      : patch outside_wall_pipe
      - escape                      = 0, 0
      - stick                       = 0, 0
    Parcel fate (number, mass)      : patch inside_wall_pipe_1
      - escape                      = 0, 0
      - stick                       = 0, 0
    Parcel fate (number, mass)      : patch symmetry_pipe
      - escape                      = 0, 0
      - stick                       = 0, 0
    Parcel fate (number, mass)      : patch symmetry_2
      - escape                      = 0, 0
      - stick                       = 0, 0
    Parcel fate (number, mass)      : patch outside_wall
      - escape                      = 0, 0
      - stick                       = 0, 0
    Parcel fate (number, mass)      : patch outlet_ambient
      - escape                      = 0, 0
      - stick                       = 0, 0
    Parcel fate (number, mass)      : patch inside_wall_pipe_2
      - escape                      = 0, 0
      - stick                       = 0, 0

DILUPBiCG:  Solving for Ucx, Initial residual = 1, Final residual = 1.25761e-009, No Iterations 2
DILUPBiCG:  Solving for Ucy, Initial residual = 1, Final residual = 9.96601e-010, No Iterations 2
DILUPBiCG:  Solving for Ucz, Initial residual = 1, Final residual = 1.65935e-006, No Iterations 1
GAMG:  Solving for p, Initial residual = 1, Final residual = 0.00719945, No Iterations 18
GAMG:  Solving for p, Initial residual = 0.157115, Final residual = 0.00130879, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.0174056, Final residual = 0.000154434, No Iterations 3
time step continuity errors : sum local = 1.30207e-007, global = -1.40036e-008, cumulative = -1.40036e-008
GAMG:  Solving for p, Initial residual = 0.00302915, Final residual = 1.90974e-005, No Iterations 5
GAMG:  Solving for p, Initial residual = 0.00119684, Final residual = 6.68533e-006, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.000211614, Final residual = 7.64506e-007, No Iterations 5
time step continuity errors : sum local = 6.32305e-010, global = -7.29831e-011, cumulative = -1.40765e-008
DILUPBiCG:  Solving for epsilon, Initial residual = 0.000435878, Final residual = 1.6674e-008, No Iterations 1


--> FOAM FATAL ERROR:

    request for volVectorField U from objectRegistry region0 failed
    available objects of type volVectorField are

3
(
Uc_0
Uc
cloudVolSUSu
)


    From function objectRegistry::lookupObject<Type>(const word&) const

From a message OF, I do not know "request for volVectorField U" where BCs are or solution?

The velocity Uc as the same Goldschmidt OF tutorial. I hope foamers give me a recommendation to solve problem.

Thanks

[hiuluom]

Hi again,

Finally, I solved the problems that some BCs require phi and U but these conditions do not get parameter for phi, U because DPM solver is using Uc and phic, rhoc for continuous or fluid phase. Therefore, it should add phic and Uc in BCs.

Here is the case

https://www.cfd-online.com/Forums/at...1&d=1496850233

[minzhang]

Hello,

I am wondering whether you could have a look at my post about MPPICFoam.
Simulate particle/liquid flow in a converging geometry with MPPICFoam

Thank you so much!