OpenFOAM chemistry solver

alicharchi · February 27, 2017, 14:15

I am working on reactive flow and I do not understand the way reactingFOAM solves the Yi equations.

To be more specific, the question is how the source terms (reaction->R(Yi)) are discretized? I am quite sure that the Ode solver is used to calculate that term in reaction->correct() member function while the temporal term (fvm::ddt(rho, Yi)) is discretized by Euler ddt scheme.

I have noticed the R member function returns a pointer of type fvScalarMatrix.
Can anyone give some insight into the procedure of discretizing the source term?

Note that I am using OF 2.3.1.

Thanks

kbeat · September 2, 2019, 03:44

step 1 The reaction-> correct() actually solves the stiff chemical ode to get R[i]
(order of integration is 10e-9 integrated upto ~ deltaT)

step 2 the species equation is solved based normal implicit scheme

Since the implementation allows interchangeable use of fvScalarMatric and volScalarField( difference being the boundary)

February 27, 2017, 14:15	OpenFOAM chemistry solver	#1
alicharchi New Member Ali Charchi Join Date: Oct 2013 Location: South Carolina, USA Posts: 8 Rep Power: 12	I am working on reactive flow and I do not understand the way reactingFOAM solves the Yi equations. To be more specific, the question is how the source terms (reaction->R(Yi)) are discretized? I am quite sure that the Ode solver is used to calculate that term in reaction->correct() member function while the temporal term (fvm::ddt(rho, Yi)) is discretized by Euler ddt scheme. I have noticed the R member function returns a pointer of type fvScalarMatrix. Can anyone give some insight into the procedure of discretizing the source term? Note that I am using OF 2.3.1. Thanks

September 2, 2019, 03:44	Reacting Flow Chemistry	#2
kbeat New Member Karthik Thyagarjan Join Date: Jul 2019 Posts: 4 Rep Power: 6	step 1 The reaction-> correct() actually solves the stiff chemical ode to get R[i] (order of integration is 10e-9 integrated upto ~ deltaT) step 2 the species equation is solved based normal implicit scheme Since the implementation allows interchangeable use of fvScalarMatric and volScalarField( difference being the boundary) alicharchi, randolph and h_kuriakose like this.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Frequently Asked Questions about Installing OpenFOAM	wyldckat	OpenFOAM Installation	3	November 14, 2023 11:58
How to rewrite a standard OpenFOAM solver as a C++ class	cfbaptista	OpenFOAM Programming & Development	7	March 23, 2016 04:50
Can't get data from OpenFoam to external solver using externalCoupled	perry	OpenFOAM Running, Solving & CFD	4	May 26, 2014 08:09
Working directory via command line	Luiz	CFX	4	March 6, 2011 20:02
CH4 gas-phase reaction and stiff chemistry solver	sonosun	FLUENT	0	April 15, 2010 21:01