implementation of combustion in fireFoam 1.6-1.0

hz283 · January 2, 2013, 16:51

Hi All,

I read the source codes of fireFoam 1.6-1.0, and have some questions about how the combustion is implemented in this code. The source is downloaded from the website: http://code.google.com/p/firefoam-dev/.

In this code, the equations about density, momentum, mixture fraction and mixture fraction variance, pressure, total enthalpy as well as the corresponding turbulence quantities are solved. However, in all these equations, there is no combustion source terms solved. For example, in the enthalpy equation:
fvm::ddt(rho, h)
+ mvConvection->fvmDiv(phi,h)
- fvm::laplacian(turbulence->alphaEff(), h)
==
dpdt
+ fvc::div(phi/fvc::interpolate(rho)*fvc::interpolate(p))
- p*fvc::div(phi/fvc::interpolate(rho))
+ radiation->Sh(thermo)

The terms in the RHS of the above expression are Dp/Dt+radiation source terms.

Does anybody have experience in using this code? Thank you in advance.

best regards,
H

adhiraj · January 3, 2013, 16:04

The total enthalpy equation should not have a chemical reaction term,
unless I am missing something.

hz283 · January 3, 2013, 16:07

yes, you are correct. I did not question the correctness of this equation, I am just ask the question about how the combustion is implemented. I read the fireFoam related papers and have already had some understanding now.

January 2, 2013, 16:51	implementation of combustion in fireFoam 1.6-1.0	#1
hz283 Senior Member Join Date: Nov 2012 Posts: 171 Rep Power: 14	Hi All, I read the source codes of fireFoam 1.6-1.0, and have some questions about how the combustion is implemented in this code. The source is downloaded from the website: http://code.google.com/p/firefoam-dev/. In this code, the equations about density, momentum, mixture fraction and mixture fraction variance, pressure, total enthalpy as well as the corresponding turbulence quantities are solved. However, in all these equations, there is no combustion source terms solved. For example, in the enthalpy equation: fvm::ddt(rho, h) + mvConvection->fvmDiv(phi,h) - fvm::laplacian(turbulence->alphaEff(), h) == dpdt + fvc::div(phi/fvc::interpolate(rho)fvc::interpolate(p)) - pfvc::div(phi/fvc::interpolate(rho)) + radiation->Sh(thermo) The terms in the RHS of the above expression are Dp/Dt+radiation source terms. Does anybody have experience in using this code? Thank you in advance. best regards, H

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
[mesh manipulation] extrudeMesh from STL Surface	TRT	OpenFOAM Meshing & Mesh Conversion	1	March 26, 2018 08:15
Modelling Combustion in Porous Zone	tanjinjack	FLUENT	2	September 26, 2016 05:10
DieselFoam spray	thumthae	OpenFOAM Running, Solving & CFD	98	December 24, 2014 16:55
[snappyHexMesh] Giving patch names within stl files	Nik_	OpenFOAM Meshing & Mesh Conversion	2	October 27, 2010 14:41
UDF : Chemical Reaction in Catalytic Converter	adhimac	FLUENT	1	March 27, 2001 15:04

January 3, 2013, 16:04		#2
adhiraj Senior Member Adhiraj Join Date: Sep 2010 Location: Karnataka, India Posts: 187 Rep Power: 16	The total enthalpy equation should not have a chemical reaction term, unless I am missing something.

January 3, 2013, 16:07		#3
hz283 Senior Member Join Date: Nov 2012 Posts: 171 Rep Power: 14	yes, you are correct. I did not question the correctness of this equation, I am just ask the question about how the combustion is implemented. I read the fireFoam related papers and have already had some understanding now.