Missing dt in OxidationKineticDiffusionLimitedRate

omarzouk · May 14, 2010, 06:00

Hello,

There is a bug in the kinetic- and diffusion-limited carbon oxidation model (to be used with coalChemistryFoam) solver.
It is in the source file:
OpenFOAM-1.6/src/lagrangian/coalCombustion/submodels/surfaceReactionModel/COxidationKineticDiffusionLimitedRate/COxidationKineticDiffusionLimitedRate.C
The following formula is missing (dt). Thus, it actually calculates dmC/dt (mass-rate of carbon consumption) instead of dmC (increment of carbon consumed over a time step):
scalar dmC = Ap*rhoc*specie::RR*Tc*YO2/WO2_*D0*Rk/(D0 + Rk);

henry · May 17, 2010, 11:07

Thanks for the bug report, this is now fixed in OpenFOAM-1.6.x.

H

May 14, 2010, 06:00	Missing dt in OxidationKineticDiffusionLimitedRate	#1
omarzouk New Member Join Date: Apr 2009 Posts: 6 Rep Power: 17	Hello, There is a bug in the kinetic- and diffusion-limited carbon oxidation model (to be used with coalChemistryFoam) solver. It is in the source file: OpenFOAM-1.6/src/lagrangian/coalCombustion/submodels/surfaceReactionModel/COxidationKineticDiffusionLimitedRate/COxidationKineticDiffusionLimitedRate.C The following formula is missing (dt). Thus, it actually calculates dmC/dt (mass-rate of carbon consumption) instead of dmC (increment of carbon consumed over a time step): scalar dmC = Aprhocspecie::RRTcYO2/WO2_D0Rk/(D0 + Rk);

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
CircumferentialAveragingInterpolationName.C file is missing	tian	OpenFOAM Bugs	1	November 29, 2009 17:43
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

May 17, 2010, 11:07		#2
henry Senior Member Join Date: Mar 2009 Posts: 854 Rep Power: 22	Thanks for the bug report, this is now fixed in OpenFOAM-1.6.x. H