CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM

Additional Parameter "C" for combustionProperties

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 26, 2022, 05:44
Default Additional Parameter "C" for combustionProperties
  #1
New Member
 
Chris
Join Date: Oct 2022
Posts: 2
Rep Power: 0
dwarfgang is on a distinguished road
Hello everyone, I am a student who is quite new to using OpenFOAM.
I recently installed a custom solver called rhoReactingCentralFoam (based on OpenFOAM v7) in order to simulate one-step global combustion in a supersonic environment (H2-Air reaction).

From setting up the combustionProperties file in the constant folder (in which I chose the diffusion combustion model), the program asked me to input a scalar value for variable "C." After searching up on the internet, this was the most detailed description I could find for the parameter.
"Additional parameter C is used to distribute the heat release rate in time."
Does it mean that I have to specify the reaction's heat release rate after manual calculation?
I was able to find a correlation for the reaction rate, which was found on the thesis at the bottom of https://repository.tudelft.nl/island...8-9eda89bb04e7 (page 64, under diffusion), but still cannot exactly understand what the variable indicates.
It would be great if someone could offer a thorough explanation on what this is.
Thank you in advance!


-Chris
dwarfgang is offline   Reply With Quote

Old   December 2, 2022, 06:54
Default
  #2
New Member
 
Rocco Tortorolo
Join Date: Dec 2022
Location: Rome
Posts: 4
Rep Power: 3
rocco_Tororolo is on a distinguished road
Hi I am doing simulations in reactingFoam and i encountred the same issues.
did you find any answers?
rocco_Tororolo is offline   Reply With Quote

Old   December 15, 2022, 01:52
Default
  #3
New Member
 
Chris
Join Date: Oct 2022
Posts: 2
Rep Power: 0
dwarfgang is on a distinguished road
Quote:
Originally Posted by rocco_Tororolo View Post
Hi I am doing simulations in reactingFoam and i encountred the same issues.
did you find any answers?
Unfortunately no I implemented an arbitrary value of C = 10, but it does not seem to suit supersonic combustion accurately as the concentration of H2 diminishes much slower compared to similar cases.
Will be looking further into this parameter once I am done writing my Undergrad research paper.
dwarfgang is offline   Reply With Quote

Old   December 16, 2022, 18:13
Default
  #4
Member
 
Ali B.
Join Date: Mar 2020
Location: abzrg.github.io
Posts: 34
Rep Power: 6
reverseila is on a distinguished road
I don't know anything about combustion...

but There are two combustion models (in $FOAM_SRC/combustionModels) that have this C parameter in them, infinitelyFastChemistry and diffusion

According to the source code, C is a private member variable of type scalar which they commented it as "Model constant"

It is read from a dictionary (not sure from where exactly but I could only find C in infinitelyFastChemistryCoeffs inside combustionProperties files of the $FOAM_TUTORIALS/combustion.)

In the source code of the above combustion models, we have something like these

Code:
// file: infinitelyFastChemistry.C
// Constructor definition

template<class ReactionThermo, class ThermoType>
void infinitelyFastChemistry<ReactionThermo, ThermoType>::correct()
{
    this->wFuel_ == dimensionedScalar(dimMass/dimVolume/dimTime, Zero);

    if (this->active())
    {
        this->singleMixturePtr_->fresCorrect();

        const label fuelI = this->singleMixturePtr_->fuelIndex();

        const volScalarField& YFuel =
            this->thermo().composition().Y()[fuelI];

        const dimensionedScalar s = this->singleMixturePtr_->s();

        if (this->thermo().composition().contains("O2"))
        {
            const volScalarField& YO2 = this->thermo().composition().Y("O2");

             // wFuel: Fuel consumption rate
            this->wFuel_ ==
                this->rho()/(this->mesh().time().deltaT()*C_)
               *min(YFuel, YO2/s.value());
        }
    }
}
...
Code:
// file: diffusion.C
...
// Constructor definition
template<class ReactionThermo, class ThermoType>
void diffusion<ReactionThermo, ThermoType>::correct()
{
    this->wFuel_ == dimensionedScalar(dimMass/dimVolume/dimTime, Zero);

    if (this->active())
    {
        this->singleMixturePtr_->fresCorrect();

        const label fuelI = this->singleMixturePtr_->fuelIndex();

        const volScalarField& YFuel =
            this->thermo().composition().Y()[fuelI];

        if (this->thermo().composition().contains(oxidantName_))
        {
            const volScalarField& YO2 =
                this->thermo().composition().Y(oxidantName_);

            this->wFuel_ ==
                C_*this->turbulence().muEff()
               *mag(fvc::grad(YFuel) & fvc::grad(YO2))
               *pos0(YFuel)*pos0(YO2);
        }
    }
}
...
Hope this helps!
reverseila is offline   Reply With Quote

Reply

Tags
c parameter, combustionproperties

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
[PyFoam] running pyFoam(Plot)Runner.py in parallel Studi OpenFOAM Community Contributions 10 September 6, 2015 16:38
Force can not converge colopolo CFX 13 October 4, 2011 22:03
emag beta feature: charge density charlotte CFX 4 March 22, 2011 09:14
plot additional parameter Jane Siemens 2 March 15, 2004 19:40
Replace periodic by inlet-outlet pair lego CFX 3 November 5, 2002 20:09


All times are GMT -4. The time now is 19:55.