This material computes thermal conductivity for a PorousMedium - fluid system, by using Thermal conductivity = dry_thermal_conductivity + S^exponent * (wet_thermal_conductivity - dry_thermal_conductivity), where S is the aqueous saturation. More...

#include <PorousFlowThermalConductivityIdeal.h>

Inheritance diagram for PorousFlowThermalConductivityIdeal:

[legend]

Public Member Functions
	PorousFlowThermalConductivityIdeal (const InputParameters &parameters)

Protected Member Functions
virtual void	computeQpProperties () override

Protected Attributes
const RealTensorValue	_la_dry
	Dry thermal conductivity of rock. More...

const bool	_wet_and_dry_differ
	Whether _la_wet has been supplied. More...

const RealTensorValue	_la_wet
	Wet thermal conductivity of rock. More...

const Real	_exponent
	Exponent for saturation. More...

const bool	_aqueous_phase
	Whether this is a fluid simulation. More...

const unsigned	_aqueous_phase_number
	Phase number of the aqueous phase. More...

const MaterialProperty< std::vector< Real > > *const	_saturation_qp
	Saturation of the fluid phases at the quadpoints. More...

const MaterialProperty< std::vector< std::vector< Real > > > *const	_dsaturation_qp_dvar
	d(Saturation)/d(PorousFlow variable) More...

MaterialProperty< RealTensorValue > &	_la_qp
	Thermal conducitivity at the qps. More...

MaterialProperty< std::vector< RealTensorValue > > &	_dla_qp_dvar
	d(thermal conductivity at the qps)/d(PorousFlow variable) More...

const unsigned int	_num_phases
	Number of phases. More...

const unsigned int	_num_components
	Number of fluid components. More...

const unsigned int	_num_var
	Number of PorousFlow variables. More...

Detailed Description

This material computes thermal conductivity for a PorousMedium - fluid system, by using Thermal conductivity = dry_thermal_conductivity + S^exponent * (wet_thermal_conductivity - dry_thermal_conductivity), where S is the aqueous saturation.

Definition at line 26 of file PorousFlowThermalConductivityIdeal.h.

Constructor & Destructor Documentation

◆ PorousFlowThermalConductivityIdeal()

PorousFlowThermalConductivityIdeal::PorousFlowThermalConductivityIdeal ( const InputParameters & parameters )

Definition at line 44 of file PorousFlowThermalConductivityIdeal.C.

   : PorousFlowThermalConductivityBase(parameters),
     _la_dry(getParam<RealTensorValue>("dry_thermal_conductivity")),
     _wet_and_dry_differ(parameters.isParamValid("wet_thermal_conductivity")),
     _la_wet(_wet_and_dry_differ ? getParam<RealTensorValue>("wet_thermal_conductivity")
                                 : getParam<RealTensorValue>("dry_thermal_conductivity")),
     _exponent(getParam<Real>("exponent")),
     _aqueous_phase(_num_phases > 0),
     _aqueous_phase_number(getParam<unsigned>("aqueous_phase_number")),
     _saturation_qp(_aqueous_phase
                        ? &getMaterialProperty<std::vector<Real>>("PorousFlow_saturation_qp")
                        : nullptr),
     _dsaturation_qp_dvar(_aqueous_phase ? &getMaterialProperty<std::vector<std::vector<Real>>>(
                                               "dPorousFlow_saturation_qp_dvar")
                                         : nullptr)
 {
   if (_aqueous_phase && (_aqueous_phase_number >= _num_phases))
     mooseError("PorousFlowThermalConductivityIdeal: Your aqueous phase number, ",
                _aqueous_phase_number,
                " must not exceed the number of fluid phases in the system, which is ",
                _num_phases,
                "\n");
 }

Member Function Documentation

◆ computeQpProperties()

void PorousFlowThermalConductivityIdeal::computeQpProperties ( )

overrideprotectedvirtual

Definition at line 70 of file PorousFlowThermalConductivityIdeal.C.

 {
   _la_qp[_qp] = _la_dry;
   if (_aqueous_phase && _wet_and_dry_differ)
     _la_qp[_qp] +=
         std::pow((*_saturation_qp)[_qp][_aqueous_phase_number], _exponent) * (_la_wet - _la_dry);
  
   _dla_qp_dvar[_qp].assign(_num_var, RealTensorValue());
   if (_aqueous_phase && _wet_and_dry_differ)
     for (unsigned v = 0; v < _num_var; ++v)
       _dla_qp_dvar[_qp][v] =
           _exponent * std::pow((*_saturation_qp)[_qp][_aqueous_phase_number], _exponent - 1.0) *
           (*_dsaturation_qp_dvar)[_qp][_aqueous_phase_number][v] * (_la_wet - _la_dry);
 }