|
[Sponsors] |
August 18, 2010, 13:51 |
Please confirm, if this is a bug or a typo
|
#1 |
Member
N. A.
Join Date: May 2010
Posts: 64
Rep Power: 16 |
Dear Foamers,
Can someone confirm the model of surface reaction using lagrangian particles. I am simulating cola-particle consumption. The consumpotio of coal is as follows: C + O2 --> CO2 dmC, cMO2 and dmCO2 are being calculated to find the amount of species consumed. Later the particle mass is updated and so is the gas phase species mass. I am wondering if dMassSolid[CslocalId_] +- dmC; (instead of -dmC) The part of file, COxidationDiffusionLimitedRate, is as below: ... ... // Local mass fraction of O2 in the carrier phase const scalar YO2 = this->owner().mcCarrierThermo().Y(O2GlobalId_)[cellI]; // Change in C mass [kg] scalar dmC = 4.0*mathematicalConstant:i*d*D_*YO2*Tc*rhoc/(Sb_*(T + Tc))*dt; // Limit mass transfer by availability of C dmC = min(mass*fComb, dmC); // Change in O2 mass [kg] const scalar dmO2 = dmC/WC_*Sb_*WO2_; // Mass of newly created CO2 [kg] const scalar dmCO2 = dmC + dmO2; // Update local particle C mass dMassSolid[CsLocalId_] += dmC; // Update carrier O2 and CO2 mass dMassSRCarrier[O2GlobalId_] -= dmO2; dMassSRCarrier[CO2GlobalId_] += dmCO2; // Heat of reaction [J] return -HcCO2_*dmCO2; } ------------- Should I change dMassSolid[CsLocalId_] += dmC; to dMassSolid[CsLocalId_] -= dmC; ? Thanks, N |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Serious bug in LES interface | fs82 | OpenFOAM Bugs | 21 | November 16, 2009 09:15 |
Maybe a bug in CFX12.0 BETA? | frank | CFX | 4 | November 26, 2008 19:08 |
Bug reports | Mattijs Janssens (Mattijs) | OpenFOAM | 0 | January 10, 2005 11:05 |
Win XP hot fix bug!!! | Emilio | Siemens | 1 | March 27, 2003 13:06 |
Forum y2k Bug | Jonas Larsson | Main CFD Forum | 1 | January 5, 2000 11:22 |