|
[Sponsors] |
Phi error in implimenting Kepsilon in DPMFoam |
![]() |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
![]() |
![]() |
#1 |
New Member
Join Date: Jun 2017
Posts: 11
Rep Power: 9 ![]() |
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 Last edited by kinbean; July 11, 2017 at 07:43. |
|
![]() |
![]() |
![]() |
![]() |
#2 |
New Member
Join Date: Jun 2017
Posts: 11
Rep Power: 9 ![]() |
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). ------------------------------------------------------------------------ |
|
![]() |
![]() |
![]() |
![]() |
#3 | |
Member
Ping Chang
Join Date: Feb 2016
Location: Perth
Posts: 93
Rep Power: 10 ![]() |
Quote:
Kind regards, Ping |
||
![]() |
![]() |
![]() |
![]() |
#4 |
Member
Robert
Join Date: Sep 2016
Posts: 32
Rep Power: 9 ![]() |
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 |
|
![]() |
![]() |
![]() |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
DPMFoam - Serious Error --particle-laden flow in simple geometric config | benz25 | OpenFOAM Running, Solving & CFD | 27 | December 19, 2017 20:47 |
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 09:05 |
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 04:33 |