CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM (https://www.cfd-online.com/Forums/openfoam/)
-   -   Phi error in implimenting Kepsilon in DPMFoam (https://www.cfd-online.com/Forums/openfoam/190324-phi-error-implimenting-kepsilon-dpmfoam.html)

kinbean July 10, 2017 20:16

Phi error in implimenting Kepsilon in DPMFoam
 
1 Attachment(s)
Hi, I ran into an error I have no idea how to adress. I'm trying to do a eularian-lagrangian analysis of particles in a pipe in crossflow. For this i decided to use DPMFoam, but it's not setup for Kepsilon and i've had to try and impliment it. However I'm getting this error I don't know how to deal with. Any help would be really appreciated.

I've treid to upload to folder without the mesh, If anyone has the time to have a look I'd really appreciate it.

I'm running OpenFOAM 4.1

let me know if you need any more info

Thank you

kinbean July 11, 2017 07:39

2 Attachment(s)
Hi this is no longer a problem, I now am stuck with something else. The simulation runs then gives up after roughly the third timestep.

It didn't let me upload the full case folder. Nothing has changed in constants from last post.

Does anyone have any ideas?
Thank you for your help


Here is the log.DPM file

Code:

Solving 3-D cloud kinematicCloud

Cloud: kinematicCloud injector: model1
    Added 101 new parcels

    389 move-collide subCycles
tail: unrecognised file system type 0x53464846 for ‘log.DPMFoam’. please report this to bug-coreutils@gnu.org. reverting to polling
Cloud: kinematicCloud
    Current number of parcels      = 101
    Current mass in system          = 1.6698e-06
    Linear momentum                = (3.28532e-05 0 -1.26026e-10)
  |Linear momentum|                = 3.28532e-05
    Linear kinetic energy          = 0.000323193
    model1:
        number of parcels added    = 101
        mass introduced            = 1.6698e-06
    Parcel fate (number, mass)      : patch inlet
      - escape                      = 0, 0
      - stick                      = 0, 0
    Parcel fate (number, mass)      : patch outlet
      - escape                      = 0, 0
      - stick                      = 0, 0
    Parcel fate (number, mass)      : patch lowerWall
      - escape                      = 0, 0
      - stick                      = 0, 0
    Parcel fate (number, mass)      : patch upperWall
      - escape                      = 0, 0
      - stick                      = 0, 0
    Parcel fate (number, mass)      : patch frontAndBack
      - escape                      = 0, 0
      - stick                      = 0, 0
    Parcel fate (number, mass)      : patch pipeScale
      - escape                      = 0, 0
      - stick                      = 0, 0
    Rotational kinetic energy      = 0

PIMPLE: iteration 1
smoothSolver:  Solving for U.airx, Initial residual = 1, Final residual = 8.43751e-06, No Iterations 4
smoothSolver:  Solving for U.airy, Initial residual = 1, Final residual = 7.63784e-06, No Iterations 4
smoothSolver:  Solving for U.airz, Initial residual = 1, Final residual = 3.05085e-06, No Iterations 2
GAMG:  Solving for p, Initial residual = 1, Final residual = 0.00525976, No Iterations 6
time step continuity errors : sum local = 7.83773e-09, global = -4.49241e-09, cumulative = -4.49241e-09
GAMG:  Solving for p, Initial residual = 0.337878, Final residual = 9.23512e-07, No Iterations 55
time step continuity errors : sum local = 5.76991e-12, global = 4.81013e-12, cumulative = -4.4876e-09
No finite volume options present

smoothSolver:  Solving for epsilon.air, Initial residual = 1.03029e-05, Final residual = 1.63543e-07, No Iterations 1
bounding epsilon.air, min: -3.34167 max: 1288.46 average: 39.2906
smoothSolver:  Solving for k.air, Initial residual = 1, Final residual = 6.73817e-06, No Iterations 7
ExecutionTime = 12.28 s  ClockTime = 13 s

Courant Number mean: 0.000964355 max: 0.40604
Time = 4e-05

Evolving kinematicCloud

Solving 3-D cloud kinematicCloud
    407 move-collide subCycles
Cloud: kinematicCloud
    Current number of parcels      = 101
    Current mass in system          = 1.6698e-06
    Linear momentum                = (9.46156e-06 -8.20201e-08 -3.47299e-07)
  |Linear momentum|                = 9.46829e-06
    Linear kinetic energy          = 3.65316e-05
    model1:
        number of parcels added    = 101
        mass introduced            = 1.6698e-06
    Parcel fate (number, mass)      : patch inlet
      - escape                      = 0, 0
      - stick                      = 0, 0
    Parcel fate (number, mass)      : patch outlet
      - escape                      = 0, 0
      - stick                      = 0, 0
    Parcel fate (number, mass)      : patch lowerWall
      - escape                      = 0, 0
      - stick                      = 0, 0
    Parcel fate (number, mass)      : patch upperWall
      - escape                      = 0, 0
      - stick                      = 0, 0
    Parcel fate (number, mass)      : patch frontAndBack
      - escape                      = 0, 0
      - stick                      = 0, 0
    Parcel fate (number, mass)      : patch pipeScale
      - escape                      = 0, 0
      - stick                      = 0, 0
    Rotational kinetic energy      = 0

PIMPLE: iteration 1
smoothSolver:  Solving for U.airx, Initial residual = 0.000491792, Final residual = 4.08364e-06, No Iterations 1
smoothSolver:  Solving for U.airy, Initial residual = 0.0346797, Final residual = 6.23164e-06, No Iterations 5
smoothSolver:  Solving for U.airz, Initial residual = 0.195751, Final residual = 9.58755e-06, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.218458, Final residual = 0.0020271, No Iterations 5
time step continuity errors : sum local = 1.47353e-08, global = -1.97603e-09, cumulative = -6.46363e-09
GAMG:  Solving for p, Initial residual = 0.62654, Final residual = 8.9094e-07, No Iterations 60
time step continuity errors : sum local = 2.99806e-12, global = -2.50408e-12, cumulative = -6.46613e-09
smoothSolver:  Solving for epsilon.air, Initial residual = 5.98403e-05, Final residual = 1.47682e-06, No Iterations 1
bounding epsilon.air, min: -27110.7 max: 3.65652e+06 average: 32024.8
smoothSolver:  Solving for k.air, Initial residual = 0.00025891, Final residual = 7.1125e-06, No Iterations 6
ExecutionTime = 17.75 s  ClockTime = 18 s

Courant Number mean: 0.000964269 max: 0.385256
Time = 6e-05

Evolving kinematicCloud

Solving 3-D cloud kinematicCloud
    362 move-collide subCycles
Cloud: kinematicCloud
    Current number of parcels      = 101
    Current mass in system          = 1.6698e-06
    Linear momentum                = (7.43209e-06 -5.82378e-08 -3.71842e-07)
  |Linear momentum|                = 7.44162e-06
    Linear kinetic energy          = 2.72171e-05
    model1:
        number of parcels added    = 101
        mass introduced            = 1.6698e-06
    Parcel fate (number, mass)      : patch inlet
      - escape                      = 0, 0
      - stick                      = 0, 0
    Parcel fate (number, mass)      : patch outlet
      - escape                      = 0, 0
      - stick                      = 0, 0
    Parcel fate (number, mass)      : patch lowerWall
      - escape                      = 0, 0
      - stick                      = 0, 0
    Parcel fate (number, mass)      : patch upperWall
      - escape                      = 0, 0
      - stick                      = 0, 0
    Parcel fate (number, mass)      : patch frontAndBack
      - escape                      = 0, 0
      - stick                      = 0, 0
    Parcel fate (number, mass)      : patch pipeScale
      - escape                      = 0, 0
      - stick                      = 0, 0
    Rotational kinetic energy      = 0

PIMPLE: iteration 1
smoothSolver:  Solving for U.airx, Initial residual = 0.704674, Final residual = 9.70184e-06, No Iterations 24
smoothSolver:  Solving for U.airy, Initial residual = 0.772893, Final residual = 9.88161e-06, No Iterations 73
smoothSolver:  Solving for U.airz, Initial residual = 0.971386, Final residual = 9.89935e-06, No Iterations 110
GAMG:  Solving for p, Initial residual = 0.707649, Final residual = 0.00314431, No Iterations 3
time step continuity errors : sum local = 2.55388e-09, global = -2.52526e-12, cumulative = -6.46866e-09
GAMG:  Solving for p, Initial residual = 0.564502, Final residual = 9.43709e-07, No Iterations 38
time step continuity errors : sum local = 0.000707365, global = 0.000593692, cumulative = 0.000593686
smoothSolver:  Solving for epsilon.air, Initial residual = 1, Final residual = 0.337702, No Iterations 1000
bounding epsilon.air, min: -571.1 max: 2.71286e+26 average: 4.21011e+21
smoothSolver:  Solving for k.air, Initial residual = 0.999799, Final residual = 0.914737, No Iterations 1000
--------------------------------------------------------------------------
mpirun has exited due to process rank 4 with PID 3677 on
node MSI exiting improperly. There are two reasons this could occur:

1. this process did not call "init" before exiting, but others in
the job did. This can cause a job to hang indefinitely while it waits
for all processes to call "init". By rule, if one process calls "init",
then ALL processes must call "init" prior to termination.

2. this process called "init", but exited without calling "finalize".
By rule, all processes that call "init" MUST call "finalize" prior to
exiting or it will be considered an "abnormal termination"

This may have caused other processes in the application to be
terminated by signals sent by mpirun (as reported here).
------------------------------------------------------------------------


chpjz0391 September 12, 2017 09:26

Quote:

Originally Posted by kinbean (Post 656621)
Hi this is no longer a problem, I now am stuck with something else. The simulation runs then gives up after roughly the third timestep.

It didn't let me upload the full case folder. Nothing has changed in constants from last post.

Does anyone have any ideas?
Thank you for your help


Here is the log.DPM file

Code:

Solving 3-D cloud kinematicCloud

Cloud: kinematicCloud injector: model1
    Added 101 new parcels

    389 move-collide subCycles
tail: unrecognised file system type 0x53464846 for ‘log.DPMFoam’. please report this to bug-coreutils@gnu.org. reverting to polling
Cloud: kinematicCloud
    Current number of parcels      = 101
    Current mass in system          = 1.6698e-06
    Linear momentum                = (3.28532e-05 0 -1.26026e-10)
  |Linear momentum|                = 3.28532e-05
    Linear kinetic energy          = 0.000323193
    model1:
        number of parcels added    = 101
        mass introduced            = 1.6698e-06
    Parcel fate (number, mass)      : patch inlet
      - escape                      = 0, 0
      - stick                      = 0, 0
    Parcel fate (number, mass)      : patch outlet
      - escape                      = 0, 0
      - stick                      = 0, 0
    Parcel fate (number, mass)      : patch lowerWall
      - escape                      = 0, 0
      - stick                      = 0, 0
    Parcel fate (number, mass)      : patch upperWall
      - escape                      = 0, 0
      - stick                      = 0, 0
    Parcel fate (number, mass)      : patch frontAndBack
      - escape                      = 0, 0
      - stick                      = 0, 0
    Parcel fate (number, mass)      : patch pipeScale
      - escape                      = 0, 0
      - stick                      = 0, 0
    Rotational kinetic energy      = 0

PIMPLE: iteration 1
smoothSolver:  Solving for U.airx, Initial residual = 1, Final residual = 8.43751e-06, No Iterations 4
smoothSolver:  Solving for U.airy, Initial residual = 1, Final residual = 7.63784e-06, No Iterations 4
smoothSolver:  Solving for U.airz, Initial residual = 1, Final residual = 3.05085e-06, No Iterations 2
GAMG:  Solving for p, Initial residual = 1, Final residual = 0.00525976, No Iterations 6
time step continuity errors : sum local = 7.83773e-09, global = -4.49241e-09, cumulative = -4.49241e-09
GAMG:  Solving for p, Initial residual = 0.337878, Final residual = 9.23512e-07, No Iterations 55
time step continuity errors : sum local = 5.76991e-12, global = 4.81013e-12, cumulative = -4.4876e-09
No finite volume options present

smoothSolver:  Solving for epsilon.air, Initial residual = 1.03029e-05, Final residual = 1.63543e-07, No Iterations 1
bounding epsilon.air, min: -3.34167 max: 1288.46 average: 39.2906
smoothSolver:  Solving for k.air, Initial residual = 1, Final residual = 6.73817e-06, No Iterations 7
ExecutionTime = 12.28 s  ClockTime = 13 s

Courant Number mean: 0.000964355 max: 0.40604
Time = 4e-05

Evolving kinematicCloud

Solving 3-D cloud kinematicCloud
    407 move-collide subCycles
Cloud: kinematicCloud
    Current number of parcels      = 101
    Current mass in system          = 1.6698e-06
    Linear momentum                = (9.46156e-06 -8.20201e-08 -3.47299e-07)
  |Linear momentum|                = 9.46829e-06
    Linear kinetic energy          = 3.65316e-05
    model1:
        number of parcels added    = 101
        mass introduced            = 1.6698e-06
    Parcel fate (number, mass)      : patch inlet
      - escape                      = 0, 0
      - stick                      = 0, 0
    Parcel fate (number, mass)      : patch outlet
      - escape                      = 0, 0
      - stick                      = 0, 0
    Parcel fate (number, mass)      : patch lowerWall
      - escape                      = 0, 0
      - stick                      = 0, 0
    Parcel fate (number, mass)      : patch upperWall
      - escape                      = 0, 0
      - stick                      = 0, 0
    Parcel fate (number, mass)      : patch frontAndBack
      - escape                      = 0, 0
      - stick                      = 0, 0
    Parcel fate (number, mass)      : patch pipeScale
      - escape                      = 0, 0
      - stick                      = 0, 0
    Rotational kinetic energy      = 0

PIMPLE: iteration 1
smoothSolver:  Solving for U.airx, Initial residual = 0.000491792, Final residual = 4.08364e-06, No Iterations 1
smoothSolver:  Solving for U.airy, Initial residual = 0.0346797, Final residual = 6.23164e-06, No Iterations 5
smoothSolver:  Solving for U.airz, Initial residual = 0.195751, Final residual = 9.58755e-06, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.218458, Final residual = 0.0020271, No Iterations 5
time step continuity errors : sum local = 1.47353e-08, global = -1.97603e-09, cumulative = -6.46363e-09
GAMG:  Solving for p, Initial residual = 0.62654, Final residual = 8.9094e-07, No Iterations 60
time step continuity errors : sum local = 2.99806e-12, global = -2.50408e-12, cumulative = -6.46613e-09
smoothSolver:  Solving for epsilon.air, Initial residual = 5.98403e-05, Final residual = 1.47682e-06, No Iterations 1
bounding epsilon.air, min: -27110.7 max: 3.65652e+06 average: 32024.8
smoothSolver:  Solving for k.air, Initial residual = 0.00025891, Final residual = 7.1125e-06, No Iterations 6
ExecutionTime = 17.75 s  ClockTime = 18 s

Courant Number mean: 0.000964269 max: 0.385256
Time = 6e-05

Evolving kinematicCloud

Solving 3-D cloud kinematicCloud
    362 move-collide subCycles
Cloud: kinematicCloud
    Current number of parcels      = 101
    Current mass in system          = 1.6698e-06
    Linear momentum                = (7.43209e-06 -5.82378e-08 -3.71842e-07)
  |Linear momentum|                = 7.44162e-06
    Linear kinetic energy          = 2.72171e-05
    model1:
        number of parcels added    = 101
        mass introduced            = 1.6698e-06
    Parcel fate (number, mass)      : patch inlet
      - escape                      = 0, 0
      - stick                      = 0, 0
    Parcel fate (number, mass)      : patch outlet
      - escape                      = 0, 0
      - stick                      = 0, 0
    Parcel fate (number, mass)      : patch lowerWall
      - escape                      = 0, 0
      - stick                      = 0, 0
    Parcel fate (number, mass)      : patch upperWall
      - escape                      = 0, 0
      - stick                      = 0, 0
    Parcel fate (number, mass)      : patch frontAndBack
      - escape                      = 0, 0
      - stick                      = 0, 0
    Parcel fate (number, mass)      : patch pipeScale
      - escape                      = 0, 0
      - stick                      = 0, 0
    Rotational kinetic energy      = 0

PIMPLE: iteration 1
smoothSolver:  Solving for U.airx, Initial residual = 0.704674, Final residual = 9.70184e-06, No Iterations 24
smoothSolver:  Solving for U.airy, Initial residual = 0.772893, Final residual = 9.88161e-06, No Iterations 73
smoothSolver:  Solving for U.airz, Initial residual = 0.971386, Final residual = 9.89935e-06, No Iterations 110
GAMG:  Solving for p, Initial residual = 0.707649, Final residual = 0.00314431, No Iterations 3
time step continuity errors : sum local = 2.55388e-09, global = -2.52526e-12, cumulative = -6.46866e-09
GAMG:  Solving for p, Initial residual = 0.564502, Final residual = 9.43709e-07, No Iterations 38
time step continuity errors : sum local = 0.000707365, global = 0.000593692, cumulative = 0.000593686
smoothSolver:  Solving for epsilon.air, Initial residual = 1, Final residual = 0.337702, No Iterations 1000
bounding epsilon.air, min: -571.1 max: 2.71286e+26 average: 4.21011e+21
smoothSolver:  Solving for k.air, Initial residual = 0.999799, Final residual = 0.914737, No Iterations 1000
--------------------------------------------------------------------------
mpirun has exited due to process rank 4 with PID 3677 on
node MSI exiting improperly. There are two reasons this could occur:

1. this process did not call "init" before exiting, but others in
the job did. This can cause a job to hang indefinitely while it waits
for all processes to call "init". By rule, if one process calls "init",
then ALL processes must call "init" prior to termination.

2. this process called "init", but exited without calling "finalize".
By rule, all processes that call "init" MUST call "finalize" prior to
exiting or it will be considered an "abnormal termination"

This may have caused other processes in the application to be
terminated by signals sent by mpirun (as reported here).
------------------------------------------------------------------------


Have you solved your problem? I met a similar problem as yours. When I run the case, the case divergence suddenly. The iteration for k and epsilon become 1000 suddenly.

Kind regards,
Ping

dussa December 19, 2017 17:11

Hi kinbean,

In your second post you say you fixed your issue with phi (which is the same problem that I am getting now).

I couldnt figure our what you had changed to get it to work between your different attached files, could you possibly elaborate on how you fixed the phi issue.

Regards,
Robert


All times are GMT -4. The time now is 13:14.