CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM

Please confirm, if this is a bug or a typo

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Display Modes
Old   August 18, 2010, 12:51
Default Please confirm, if this is a bug or a typo
N. A.
Join Date: May 2010
Posts: 64
Rep Power: 9
N. A. is on a distinguished road
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; ?

N. A. is offline   Reply With Quote


Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On

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

All times are GMT -4. The time now is 00:31.