Learning Forum

[Vijay Narayan]

Hi Ghislian,

There is no known issue with mass conservation.

Reaction Kinetics only dictates the rate of formation of species.

In the plot, the mass corresponds to the total mass within the closed domain ?

I'm assuming you have only 1 cell zone and the mass is reported as Volume intergral report. 

Also, I'm assuming you are running it as Transient and the domain is a premixed system of your reactant species ( I see a mention of CH4, but there is an oxidizer also rt ?). Plus patching the entire domain to high temp of 2400K could trigger the reaction in the complete domain and could pose challenges in the simulation run.

In Summary, I would require details based on my above comments to assist your further.

[Ghislain Madiot]

 

Hi Vijay,

Thank you for your reply.

I can confirm your assumptions are correct.

The plotrepresents the total mass within the enclosed domain, and this mass is calculated as a volume integral.

This simulation is transient.

The initial conditions are set with methane (CH4) at 100%, a temperature of 2400K, and a pressure of 39 atm.

The reaction modeled is the one-step thermal decomposition of methane: CH4 ==> 2H2 + C.

Additionally, I have utilized this method for more complex reactions and compared the results with Chemkin, finding them to be similar, without any issues during the simulation runs.

Ghislain

 

 

[Vijay Narayan]

Hi Ghislian,

Thanks for sharing the details and sorry for the delay in response. Could you share the below details as well.

a) FR/No TCI or ED/FR model and which chemistry solver you are using.

b) Also what is the rate exponent you are using ? For global reactions the rate exponent need not be unity.

c) what is the time step size and solver (SIMPLEC or PISO) you are using ?

[Vijay Narayan]

One  more thing, what is the density formulation being used. By default it is incompressible ideal gas. Please change it to ideal gas formulation as your closed system pressure is also changing.

Try using a smaller time step of 1e-8 sec - you can increase it later and share your observation

 

[Ghislain Madiot]

Hi Vijay,

Here the information you asked:

• I use Finit Rate / No TCI with the Explicit source for chemistry solver.
• Density formulation in the initial case is ideal gas 
• Rate exposure is equal to 1 
• The initial time step in the plot were 10e-5 which corresponds to the Chemkin methodology, but I noticed the same behavior with lower time step like 10e-7. I quick run a case with 10e-8 same result.
• Scheme used is Coupled, however I ran the same case with SIMPLEC and PISO and obtained even higher value of mass

 

I put the screenshots for the differents menu

Perhaps, it might be usefull if I share the model or the setting file with you.

Regards,

Ghislain