[Swagga5aur]

Hello Foamers,

I am trying to make a chemical reactor solver, with both catalyst temperature and the thermal properties of the reacting fluid.

The issue is with the reaction.R(Yi) in YiEqn, as its solved as a function of the fluid temperature, however, it should be based upon the catalyst temperature and I'm unable to alter the reaction.R in the solver. Therefore, I'm trying to define a new fvMatrix based upon the reaction.R(Yi) which I then may alter to depend on catalyst temperature and remove its dependence on the reacting fluid temperature.

I tried to do this be adding the following code to the following directories.

Createfields:

PtrList<fvScalarMatrix> RrateReactor(reactorRegions.size());

    Info<< "    Adding to RrateReactor\n" << endl;
    RrateReactor.set
    (
        i,
        new fvScalarMatrix
        (
            IOobject
            (
                "Rrate",
                runTime.timeName(),
                reactorRegions[i],
                IOobject::NO_READ,
                IOobject::AUTO_WRITE
            ),
            reactorRegions[i],
	    fvMatrix<>("Rrate", dimEnergy/dimVolume/dimTime,0)
        )
    );

SetRegionFields

    	    fvMatrix& Rrate = reaction.R(Y[i]);

YiEqn;

    fvScalarMatrix& Rrate=RrateReactor[i];

However, this results in fvMatrix typing issues with the fvMatrix<>, and was hoping anyone could help me with the fvMatrix typing or have another solution approach for the above issue.

Here is a link to the solver which I'm trying to add the above in the reactor domain for OpenFOAM v4: https://www.dropbox.com/s/affha5r5ht...am.tar.gz?dl=0

Regards Lasse.

[faeze.d]

Hi,
have a look at this thread:
Species Reaction Rate RR()
I think this would be helpful

[Swagga5aur]

Hi Faezeh,
Thank you for the link that helped a lot
I have attached a reactingRateFoam solver to this post for others to look at (for ofv6).

https://www.dropbox.com/s/6si5ecxjzm...am.tar.xz?dl=0

Going to implement it in a multi region solver as well.