 |
 |
 |
|
|
|
[Sponsors] | ||||
Phi error in implimenting Kepsilon in DPMFoam |
 |
|
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
July 10, 2017, 21:16
|
Phi error in implimenting Kepsilon in DPMFoam
|
#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 08:43. |
|
|
|
|
|
July 11, 2017, 08:39
|
|
#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).
------------------------------------------------------------------------
|
|
|
|
|
|
|
September 12, 2017, 10:26
|
|
#3 | |
|
Member
Ping Chang
Join Date: Feb 2016
Location: Perth
Posts: 93
Rep Power: 10 |
Quote:
Kind regards, Ping |
||
|
|
|
||
|
December 19, 2017, 18:11
|
|
#4 |
|
Member
Robert
Join Date: Sep 2016
Posts: 32
Rep Power: 10 |
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 |
|
|
|
|
|
|
|
|
Similar Threads
|
||||
| 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 21:47 |
| Question about phi | jeanpinto24| | OpenFOAM Pre-Processing | 0 | April 15, 2017 05:13 |
| Correction procedure on a non-orthogonal skewed Mesh | me.ouda | OpenFOAM Programming & Development | 0 | January 18, 2017 10:05 |
| mixerVesselAMI2D's mass is not balancing | sharonyue | OpenFOAM Running, Solving & CFD | 6 | June 10, 2013 10:34 |
| Turbulence Model phi vs phi_ | doug | OpenFOAM Running, Solving & CFD | 4 | November 10, 2009 05:33 |
Linear Mode
