View Issue Details
| ID | Project | Category | View Status | Date Submitted | Last Update |
|---|---|---|---|---|---|
| 0000395 | OpenFOAM | Bug | public | 2012-01-25 17:16 | 2012-01-30 14:00 |
| Reporter | Assigned To | ||||
| Priority | normal | Severity | minor | Reproducibility | have not tried |
| Status | resolved | Resolution | no change required | ||
| Summary | 0000395: concentrations updated by molar stoichiometric coefficients | ||||
| Description | In openFoam2.0.x in the ODEChemistryModel.C-file in the updateConcsInReactionI-function it updates the concentrations by Code: scalar sl = R.lhs()[s].stoichCoeff; c[si] -= dt*sl*omeg; but since omega is in [kg/(m^3*s] the mass stoichiometric coefficient (normalised) instead of the molar stoichiometric coefficient should be taken, right ?? gregor | ||||
| Tags | No tags attached. | ||||
|
|
You are looking at the wrong place. The internal chemistry integration is using molar concentration, but only to obtain the new concentrations. Scroll down to the bottom of the file, where you will see that the reaction rate is linearised and multiplied with the mol. weight. dc = c - c0; for (label i=0; i<nSpecie_; i++) { RR_[i][celli] = dc[i]*specieThermo_[i].W()/deltaT; } |
|
|
Ok I somehow missed to see that molar concentrations are used Thanks |
| Date Modified | Username | Field | Change |
|---|---|---|---|
| 2012-01-25 17:16 |
|
New Issue | |
| 2012-01-26 13:04 |
|
Note Added: 0000958 | |
| 2012-01-26 14:32 |
|
Note Added: 0000960 | |
| 2012-01-30 14:00 |
|
Status | new => resolved |
| 2012-01-30 14:00 |
|
Fixed in Version | => 2.1.x |
| 2012-01-30 14:00 |
|
Resolution | open => no change required |
| 2012-01-30 14:00 |
|
Assigned To | => user2 |
