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

Phi error in implimenting Kepsilon in DPMFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 10, 2017, 20:16
Default Phi error in implimenting Kepsilon in DPMFoam
  #1
New Member
 
Join Date: Jun 2017
Posts: 8
Rep Power: 2
kinbean is on a distinguished road
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
Attached Files
File Type: zip partlamRe40 - Copy.zip (132.1 KB, 1 views)

Last edited by kinbean; July 11, 2017 at 07:43.
kinbean is offline   Reply With Quote

Old   July 11, 2017, 07:39
Default
  #2
New Member
 
Join Date: Jun 2017
Posts: 8
Rep Power: 2
kinbean is on a distinguished road
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).
------------------------------------------------------------------------
Attached Files
File Type: zip 0.orig.zip (6.8 KB, 1 views)
File Type: zip system.zip (8.8 KB, 1 views)
kinbean is offline   Reply With Quote

Old   September 12, 2017, 09:26
Default
  #3
Member
 
Ping Chang
Join Date: Feb 2016
Location: Perth
Posts: 92
Rep Power: 3
chpjz0391 is on a distinguished road
Quote:
Originally Posted by kinbean View Post
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
chpjz0391 is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Question about phi jeanpinto24| OpenFOAM Pre-Processing 0 April 15, 2017 04:13
Correction procedure on a non-orthogonal skewed Mesh me.ouda OpenFOAM Programming & Development 0 January 18, 2017 10:05
DPMFoam - Serious Error --particle-laden flow in simple geometric config benz25 OpenFOAM Running, Solving & CFD 26 July 21, 2016 11:44
mixerVesselAMI2D's mass is not balancing sharonyue OpenFOAM Running, Solving & CFD 6 June 10, 2013 09:34
Turbulence Model phi vs phi_ doug OpenFOAM Running, Solving & CFD 4 November 10, 2009 05:33


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