Temperature-dependent diffusion

Usability
1

Hello everyone,

I would like to use the ComponentTransport process to perform HTC simulations and reproduce the work of Galsa et al. (2025), but I have encountered an issue with the diffusion coefficient.

I would like to prescribe a temperature-dependent diffusion coefficient with a quadratic dependence on temperature, that is, a form involving (T^2).

I first tried to use the TemperatureDependentDiffusion property type, but from what I understand, this implementation only supports a predefined temperature dependence and not a general quadratic expression.

So my question is:

Is there currently a way in OGS to define a quadratic temperature-dependent diffusion coefficient for pore_diffusion in ComponentTransport?

More specifically, is there a way, as for some other properties such as density or viscosity, to use a customized Function type for pore_diffusion?

Thank you in advance for your help.

2

Hello rose-nelly,
yes there is also the Function tag for properties. I think something like this should work.

<property>
   <name>pore_diffusion</name>
   <type>Function</type>
      <value>
         <expression>(temperature * temperature) *1e-10</expression>
      </value>
</property>

Hope i could help!
Best
Max