[Tobi]

Hey all,

at the moment I am checking how FOAM is calculating the source terms for different species. For that purpose I checked the chemFoam. Here we will find the header file where we calculate the species i, named YEqn.H. In that file we see that we call a function from the chemistry object named RR(label).

{                                                                               
    forAll(Y, specieI)                                                          
    {                                                                           
        volScalarField& Yi = Y[specieI];                                        
                                                                                
        solve                                                                   
        (                                                                       
            fvm::ddt(rho, Yi) - chemistry.RR(specieI),                          
            mesh.solver("Yi")                                                   
        );                                                                      
    }                                                                           
}

Howver If I make a source code search for a function of type RR(const label i), I find the function in thermophysicalModels/solidChemistryModel/basicSolidChemistryModel/basicSolidChemistryModel.C

However, it seems that this function just returns zero, right?

const Foam::DimensionedField<Foam::scalar, Foam::volMesh>&                      
Foam::basicSolidChemistryModel::RR(const label i) const                         
{                                                                               
    NotImplemented;                                                             
    return (DimensionedField<scalar, volMesh>::null());                         
}                                                                               
                                                                                
                                                                                
Foam::DimensionedField<Foam::scalar, Foam::volMesh>&                            
Foam::basicSolidChemistryModel::RR(const label i)                               
{                                                                               
    NotImplemented;                                                             
                                                                                
    return dynamic_cast<DimensionedField<scalar, volMesh>&>                     
    (                                                                           
        const_cast<DimensionedField<scalar, volMesh>& >                         
        (                                                                       
            DimensionedField<scalar, volMesh>::null()                           
        )                                                                       
    );                                                                          
}

I am confused :P

[heksel8i]

Quick post:

Apparently the access function is in basicChemistryModel.H and then in the constructor of chemistryModel.C one creates new fields RR_ pointing to the original RR and the RR_ are updated in the end of the solve() routine in chemistryModel.C.

[ebringley]

Hi Tobias,

Were you able figure out how to access the source term? I'm using reactingFoam for a laminar combustion simulation but cannot figure out how to access the RR_ field updated in the reaction->correct(); routine to see the value at runtime.

Thanks,
Eric

[Tobi]

Yes I was able. Check the source and you will get it.

[ebringley]

Hi Tobias,

I think I may be incorrect to have asked my question here. We are trying to do the same thing, but using different methods: You are looking at chemFoam which uses a psiChemistryModel and the source term is called through the chemistry.RR(label specI) method. I am looking at this from reactingFoam which uses a combustionModels:: psiCombustionModel and calls the source term reaction->R(Yi) which returns an fvMatrix. This R(Yi) member function of a combustion model accesses the protected chemistryPtr_ to access the a-fore-mentioned RR function. To my knowledge, I cannot access the chemistry pointer from a combustionModel. I was wondering if there was some C++ knowledge I am missing here to access this pointer to the chemistry.

Some of my questions are as follows:
1. How can I access, from a psiCombustionModel, the chemistry source term? I'd like to print to screen, or ideally, write to file to compare different chemistry methods. It may be that I should switch to chemFoam for this...
2. in the reaction->R(Yi) method, it returns a fvMatrix with units dimMass/dimTime, by adding the RR_[i] of the chemistryModel (kg/m3/s) (accessed through the RR(label i) inline chemistryModel method). What is the reason for this? Why do the units change?
3. Is there any specific references I should look at regarding the way the chemistry models / combustion models / solvers are programmed? I'd like to be able to know whats the best place to find stuff like this on my own.


PS. I am using OpenFOAM v5. There are changes in v6 that may cause some confusion, such as the chemistryModel class is renamed in v6.

[zhangyan]

fvMatrix.C:

template<class Type>
Foam::tmp<Foam::fvMatrix<Type>> Foam::operator+
(
    const fvMatrix<Type>& A,
    const DimensionedField<Type, volMesh>& su
)
{
    checkMethod(A, su, "+");
    tmp<fvMatrix<Type>> tC(new fvMatrix<Type>(A));
    tC.ref().source() -= su.mesh().V()*su.field();
    return tC;
}

[ebringley]

Thank you. I didn't realise the + operator took volume into account.