Boundary conditions and ignition hyfoam

[PranaykumarSingeetham]

What I understood was that this solver calculates the reaction rates using arrhenius kinetic expression and turbulence chemistry coupling is ignored.I will be working on this turbulence chemistry coupling as a part of my phd work.

Currently, I am simulating supersonic combustion with your hyfoam solver. To know the boundary conditions, I read createFields header file.I have few doubts.

I could not able find the terminology of "kappatr" and "kappave" and how to ignite the mixture ? Do I need to patch the region with high temperature as similar to the reactingFoam cases ? Do I have to give energy released per unit time (dQ) as B.C's/I.C' condition or is this field created during the calculations or what exactly this dQ means here?

Thanks

[vincentcasseau]

Hi Pranay,

1. The terminology is given in my PhD thesis that you can find in the doc/ folder.
   kappa -> thermal conductivity
   subscript tr refers to the trans-rotational energy mode
   subscript ve refers to the vibro-electronic energy mode

2. Make sure chemistry is enabled in constant/hTCProperties and constant/chemistryProperties and that your set of reactions is present in the relevant constant/hTCReactionsXXX dict. Please refer to the Wiki on this topic: https://github.com/vincentcasseau/hyStrath/wiki/How-to-::-Chemistry

3. dQ works the same way as in reactingFoam. It's just a field being printed and you gave the right meaning. If you're seeing zeros in the entire domain, it means you did not enable chemistry.

Thanks

[jimmyjohnh]

Hi Pranay,

Out of curiosity, what type of TCI are you thinking about? Note that it would be good to first run the code without TCI on a test case (2D). I did work on that at some point but gave up due to the unavailability of computational power.

cheers JJ

On Mon, 28 Oct 2019 at 13:36, Vince notifications@github.com wrote:

Hi Pranay,

1.

The terminology is given in my PhD thesis that you can find in the doc/ folder. kappa -> thermal conductivity subscript tr refers to the trans-rotational energy mode subscript ve refers to the vibro-electronic energy mode 2.

Make sure chemistry is enabled in constant/hTCProperties and constant/chemistryProperties and that your set of reactions is present in the relevant constant/hTCReactionsXXX dict. Please refer to the Wiki on this topic: https://github.com/vincentcasseau/hyStrath/wiki/How-to-::-Chemistry 3.

dQ works the same way as in reactingFoam. It's just a field being printed and you gave the right meaning. If you're seeing zeros in the entire domain, it means you did not enable chemistry.

Thanks

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/vincentcasseau/hyStrath/issues/28?email_source=notifications&email_token=ADALDFS2KUNE7F4TS4LICCTQQ3MERA5CNFSM4JFZHZT2YY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOECMXCWI#issuecomment-546926937, or unsubscribe https://github.com/notifications/unsubscribe-auth/ADALDFRQZJXBHTXAJWCXWOLQQ3MERANCNFSM4JFZHZTQ .

-- Dr. ir. Jimmy-John Hoste Research Fellow at German Aerospace Center (DLR) Institute of Aerodynamics and Flow Technology ---- Spacecraft Bunsenstraße 10 | D-37073 Göttingen | Germany

[PranaykumarSingeetham]

Hi Pranay,

1. The terminology is given in my PhD thesis that you can find in the doc/ folder.
   _kappa_ -> thermal conductivity
   subscript _tr_ refers to the trans-rotational energy mode
   subscript _ve_ refers to the vibro-electronic energy mode

2. Make sure chemistry is enabled in _constant/hTCProperties_ and _constant/chemistryProperties_ and that your set of reactions is present in the relevant _constant/hTCReactionsXXX_ dict. Please refer to the Wiki on this topic: https://github.com/vincentcasseau/hyStrath/wiki/How-to-::-Chemistry

3. dQ works the same way as in reactingFoam. It's just a field being printed and you gave the right meaning. If you're seeing zeros in the entire domain, it means you did not enable chemistry.

Thanks

Thanks. I will go through your PhD thesis and other documents. I think I made some mistakes while converting twotemperature to a singletemperature solver. I will try to rectify them.

[PranaykumarSingeetham]

Hi Pranay, Out of curiosity, what type of TCI are you thinking about? Note that it would be good to first run the code without TCI on a test case (2D). I did work on that at some point but gave up due to the unavailability of computational power. cheers JJ … On Mon, 28 Oct 2019 at 13:36, Vince @.**> wrote: Hi Pranay, 1. The terminology is given in my PhD thesis that you can find in the doc/ folder. kappa -> thermal conductivity subscript tr refers to the trans-rotational energy mode subscript ve refers to the vibro-electronic energy mode 2. Make sure chemistry is enabled in constant/hTCProperties and constant/chemistryProperties and that your set of reactions is present in the relevant constant/hTCReactionsXXX* dict. Please refer to the Wiki on this topic: https://github.com/vincentcasseau/hyStrath/wiki/How-to-::-Chemistry 3. dQ works the same way as in reactingFoam. It's just a field being printed and you gave the right meaning. If you're seeing zeros in the entire domain, it means you did not enable chemistry. Thanks — You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub <#28?email_source=notifications&email_token=ADALDFS2KUNE7F4TS4LICCTQQ3MERA5CNFSM4JFZHZT2YY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOECMXCWI#issuecomment-546926937>, or unsubscribe https://github.com/notifications/unsubscribe-auth/ADALDFRQZJXBHTXAJWCXWOLQQ3MERANCNFSM4JFZHZTQ . -- Dr. ir. Jimmy-John Hoste Research Fellow at German Aerospace Center (DLR) Institute of Aerodynamics and Flow Technology ---- Spacecraft Bunsenstraße 10 | D-37073 Göttingen | Germany

Thank you very much for your suggestion. I am trying to model scramjet without turbuelnce chemistry coupling as you already said that would be a good point to start with it. I just want to see some simulations running :-). I am thinking about EDC model for a validation study. I also wanted to study the non premixed combustion with flamelet modeling in OpenFOAM. So, I am currently reading computational combustion theory and its implementation. With some interest, I would like to know what type of simulations that you had worked ? Were you trying to do LES in Combustion ?

[vincentcasseau]

Jimmy-John's publications (PhD, papers) are in our Wiki: Publications and you will find answers to your questions there.

As a side note, you can NOT run a LES simulation in hyStrath yet for the same reason that there is no chemistry-turbulence interaction model yet. It was not a priority, a bit of work is needed to reuse what is available for all single-temperature OpenFOAM solvers (edits to be made: alpha -> alphatr, kappa -> kappatr, T -> Tt).

[PranaykumarSingeetham]

Jimmy-John's publications (PhD, papers) are in our Wiki: Publications and you will find answers to your questions there.

As a side note, you can NOT run a LES simulation in hyStrath yet for the same reason that there is no chemistry-turbulence interaction model yet. It was not a priority, a bit of work is needed to reuse what is available for all single-temperature OpenFOAM solvers (edits to be made: alpha -> alphatr, kappa -> kappatr, T -> Tt).

Okay. Thanks for the information. :-)

[jimmyjohnh]

Hi Pranay,

If I recall there is an EDC solver somewhere in openfoam, perhaps an older version. I did "laminar chemistry" and EDM for TCI as it runs fast but of course a huge assumption. EDC is extremely costly, would advice PaSR then or maybe the LEM. LES I wouldn't advice in general, the cost is huge, of course the advantage is that the TCI should have a smaller impact as it involves the modeled scale if not mistaken.

cheers JJ

On Tue, 29 Oct 2019 at 08:39, PranaykumarSingeetham < notifications@github.com> wrote:

Hi Pranay, Out of curiosity, what type of TCI are you thinking about? Note that it would be good to first run the code without TCI on a test case (2D). I did work on that at some point but gave up due to the unavailability of computational power. cheers JJ … <#m-5148405090621158539> On Mon, 28 Oct 2019 at 13:36, Vince @.**> wrote: Hi Pranay, 1. The terminology is given in my PhD thesis that you can find in the doc/ folder. kappa -> thermal conductivity subscript tr refers to the trans-rotational energy mode subscript ve refers to the vibro-electronic energy mode 2. Make sure chemistry is enabled in constant/hTCProperties and constant/chemistryProperties and that your set of reactions is present in the relevant constant/hTCReactionsXXX* dict. Please refer to the Wiki on this topic: https://github.com/vincentcasseau/hyStrath/wiki/How-to-::-Chemistry 3. dQ works the same way as in reactingFoam. It's just a field being printed and you gave the right meaning. If you're seeing zeros in the entire domain, it means you did not enable chemistry. Thanks — You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub <#28 https://github.com/vincentcasseau/hyStrath/issues/28?email_source=notifications&email_token=ADALDFS2KUNE7F4TS4LICCTQQ3MERA5CNFSM4JFZHZT2YY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOECMXCWI#issuecomment-546926937>, or unsubscribe https://github.com/notifications/unsubscribe-auth/ADALDFRQZJXBHTXAJWCXWOLQQ3MERANCNFSM4JFZHZTQ . -- Dr. ir. Jimmy-John Hoste Research Fellow at German Aerospace Center (DLR) Institute of Aerodynamics and Flow Technology ---- Spacecraft Bunsenstraße 10 | D-37073 Göttingen | Germany

Thank you very much for your suggestion. I am trying to model scramjet without turbuelnce chemistry coupling as you already said that would be a good point to start with it. I just want to see some simulations running :-). I am thinking about EDC model for a validation study. I also wanted to simulate non premixed combustion with flamelet modeling in OpenFOAM. I am currently reading computational combustion theory and its implementation. With some interest, I would like to know what type of simulations that you had worked ? Were you trying to do LES in Combustion ?

— You are receiving this because you commented. Reply to this email directly, view it on GitHub https://github.com/vincentcasseau/hyStrath/issues/28?email_source=notifications&email_token=ADALDFRLZXJU5TWH75U7K7LQQ7SDDA5CNFSM4JFZHZT2YY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOECPQ3MQ#issuecomment-547294642, or unsubscribe https://github.com/notifications/unsubscribe-auth/ADALDFVBTNQSN7NVNP2RAP3QQ7SDDANCNFSM4JFZHZTQ .

-- Dr. ir. Jimmy-John Hoste Research Fellow at German Aerospace Center (DLR) Institute of Aerodynamics and Flow Technology ---- Spacecraft Bunsenstraße 10 | D-37073 Göttingen | Germany

[PranaykumarSingeetham]

Hi Pranay, If I recall there is an EDC solver somewhere in openfoam, perhaps an older version. I did "laminar chemistry" and EDM for TCI as it runs fast but of course a huge assumption. EDC is extremely costly, would advice PaSR then or maybe the LEM. LES I wouldn't advice in general, the cost is huge, of course the advantage is that the TCI should have a smaller impact as it involves the modeled scale if not mistaken. cheers JJ On Tue, 29 Oct 2019 at 08:39, PranaykumarSingeetham < notifications@github.com> wrote: … Hi Pranay, Out of curiosity, what type of TCI are you thinking about? Note that it would be good to first run the code without TCI on a test case (2D). I did work on that at some point but gave up due to the unavailability of computational power. cheers JJ … <#m-5148405090621158539> On Mon, 28 Oct 2019 at 13:36, Vince @.**> wrote: Hi Pranay, 1. The terminology is given in my PhD thesis that you can find in the doc/ folder. kappa -> thermal conductivity subscript tr refers to the trans-rotational energy mode subscript ve refers to the vibro-electronic energy mode 2. Make sure chemistry is enabled in constant/hTCProperties and constant/chemistryProperties and that your set of reactions is present in the relevant constant/hTCReactionsXXX* dict. Please refer to the Wiki on this topic: https://github.com/vincentcasseau/hyStrath/wiki/How-to-::-Chemistry 3. dQ works the same way as in reactingFoam. It's just a field being printed and you gave the right meaning. If you're seeing zeros in the entire domain, it means you did not enable chemistry. Thanks — You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub <#28 <#28>?email_source=notifications&email_token=ADALDFS2KUNE7F4TS4LICCTQQ3MERA5CNFSM4JFZHZT2YY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOECMXCWI#issuecomment-546926937>, or unsubscribe https://github.com/notifications/unsubscribe-auth/ADALDFRQZJXBHTXAJWCXWOLQQ3MERANCNFSM4JFZHZTQ . -- Dr. ir. Jimmy-John Hoste Research Fellow at German Aerospace Center (DLR) Institute of Aerodynamics and Flow Technology ---- Spacecraft Bunsenstraße 10 | D-37073 Göttingen | Germany Thank you very much for your suggestion. I am trying to model scramjet without turbuelnce chemistry coupling as you already said that would be a good point to start with it. I just want to see some simulations running :-). I am thinking about EDC model for a validation study. I also wanted to simulate non premixed combustion with flamelet modeling in OpenFOAM. I am currently reading computational combustion theory and its implementation. With some interest, I would like to know what type of simulations that you had worked ? Were you trying to do LES in Combustion ? — You are receiving this because you commented. Reply to this email directly, view it on GitHub <#28?email_source=notifications&email_token=ADALDFRLZXJU5TWH75U7K7LQQ7SDDA5CNFSM4JFZHZT2YY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOECPQ3MQ#issuecomment-547294642>, or unsubscribe https://github.com/notifications/unsubscribe-auth/ADALDFVBTNQSN7NVNP2RAP3QQ7SDDANCNFSM4JFZHZTQ . -- Dr. ir. Jimmy-John Hoste Research Fellow at German Aerospace Center (DLR) Institute of Aerodynamics and Flow Technology ---- Spacecraft Bunsenstraße 10 | D-37073 Göttingen | Germany

Thank you very much for your advice and giving me an overview on computational aspects of different TCI models. Once after completing the basic test case, then I will implement the TCI.