Non-isothermal component transport

Dear OGS community,
I noticed that non-isothermal component transport (HTC) was recently included in OGS.

  1. I was wondering, if variable density and viscosity can depend on both temperature and component concentration?
  2. Can fractures be used in the HTC modeling?
  3. Are there any test or benchmark file for an HTC model?

I would be glad about any help. Thank you.
Jonas

Hi Jonas,

the non-isothermal “feature” was added in merge request 4837.

regarding your answers:

  1. as far as I can see in the code, temperature is passed on to the material models and should work. A concentration is set, too, but the data structure passed to the material models only allows for a single concentration, so far. So, if you only have one chemical species, it should work, if you have more you should double-check the results. Maybe somebody else can help in that case.
  2. (and 3.) Said merge request contains an example with “fracture” in its name. Maybe that can serve as a starting point for you.

Best regards,
Christoph

Hi Christoph,
thanks you a lot for your answer and the given hints. One species would be fine for me and I will check the fracture benchmark.
Best, Jonas

Hi Jonas,

Thanks for your interest in HTC model with fractures.

Currently, the variable density cannot depend on the temperature and component concentration at the same time, this is limited by the material properties. If you have specific expressions, we can consider to further implement it.

For your reference, please check the recently published paper using HTC model with fracture networks: Modeling of heat and solute transport in a fracture-matrix mine thermal energy storage system

Best regards,
Chaofan