|
[Sponsors] |
Missing division by dt in COxidationMurphyShaddix |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
May 17, 2010, 15:43 |
Missing division by dt in COxidationMurphyShaddix
|
#1 |
New Member
Join Date: Apr 2009
Posts: 6
Rep Power: 17 |
Hello,
In the Murphy-Shaddix submodel for char oxidation: src/lagrangian/coalCombustion/submodels/surfaceReactionModel/COxidationMurphyShaddix/COxidationMurphyShaddix.C In the statement: const scalar qVol = sum(dMassVolatile)/(WVol_*Ap); This is supposed to be a molar rate [kmol/m^2/s], so it should be divided by dt. |
|
May 17, 2010, 15:59 |
Calculating heat-of-reaction in submodel: COxidationMurphyShaddix
|
#2 |
New Member
Join Date: Apr 2009
Posts: 6
Rep Power: 17 |
In the char-oxidation submodels under:
OpenFOAM-1.6/src/lagrangian/coalCombustion/submodels/surfaceReactionModel/ for: C(s) + O2(g) => CO2(g) [COxidationDiffusionLimitedRate] and [COxidationKineticDiffusionLimitedRate] calculate the heat of reaction based on the enthalpies of formation (i.e., chemical specific enthalpy), whereas [COxidationMurphyShaddix] calculates it based on the total (i.e., chemical plus sensible specific enthalpy) specific enthalpies. |
|
May 18, 2010, 12:22 |
|
#3 |
Senior Member
Join Date: Mar 2009
Posts: 854
Rep Power: 22 |
Are you running OpenFOAM-1.6.x? It appears that the code snippet you posted is out of date and the latest version is in our git repository.
H |
|
May 18, 2010, 15:19 |
|
#4 |
New Member
Join Date: Apr 2009
Posts: 6
Rep Power: 17 |
Thanks for the reply.
I am running 1.6. However, I see that you advise me to switch to 1.6-x which seems to have addressed several implementation mistakes in 1.6. Thanks. |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
UDF: DEFINE_CG_MOTION for vertical jump motion of an electrode! | alban | Fluent UDF and Scheme Programming | 2 | June 8, 2010 18:54 |
Problem with paraFoam | ata | OpenFOAM | 14 | November 30, 2009 04:23 |
Missing output files | stefanke | OpenFOAM Running, Solving & CFD | 4 | September 8, 2006 04:51 |
Missing tets along line when meshing with ICEM CFD | Georges P. Côté | CFX | 6 | March 23, 2006 00:34 |
about the missing component by setting BC. | Kim | FLUENT | 1 | May 28, 2004 11:33 |