Material designed to form a std::vector of property and derivatives of these wrt the nonlinear variables from the individual phase properties. More...

#include <PorousFlowJoiner.h>

Inheritance diagram for PorousFlowJoiner:

[legend]

Public Member Functions
	PorousFlowJoiner (const InputParameters &parameters)

Protected Member Functions
virtual void	initQpStatefulProperties () override

virtual void	computeQpProperties () override

Protected Attributes
const std::string	_pf_prop
	Name of material property to be joined. More...

const MaterialProperty< std::vector< std::vector< Real > > > &	_dporepressure_dvar
	Derivatives of porepressure variable wrt PorousFlow variables at the qps or nodes. More...

const MaterialProperty< std::vector< std::vector< Real > > > &	_dsaturation_dvar
	Derivatives of saturation variable wrt PorousFlow variables at the qps or nodes. More...

const MaterialProperty< std::vector< Real > > &	_dtemperature_dvar
	Derivatives of temperature variable wrt PorousFlow variables at the qps or nodes. More...

MaterialProperty< std::vector< Real > > &	_property
	Computed property of the phase. More...

MaterialProperty< std::vector< std::vector< Real > > > &	_dproperty_dvar
	d(property)/d(PorousFlow variable) More...

std::vector< const MaterialProperty< Real > * >	_phase_property
	Property of each phase. More...

std::vector< const MaterialProperty< Real > * >	_dphase_property_dp
	d(property of each phase)/d(pressure) More...

std::vector< const MaterialProperty< Real > * >	_dphase_property_ds
	d(property of each phase)/d(saturation) More...

std::vector< const MaterialProperty< Real > * >	_dphase_property_dt
	d(property of each phase)/d(temperature) 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

Material designed to form a std::vector of property and derivatives of these wrt the nonlinear variables from the individual phase properties.

Values at the quadpoint or the nodes are formed depending on _at_qps

Properties can be viscosities, densities, thermal conductivities , etc and the user specifies the property they are interested in using the pf_prop string.

Also, using d(property)/dP, d(property)/dS, etc, and dP/dvar, dS/dvar, etc, the matrix of derivatives of property with respect to the nonlinear Variables, var, are computed.

Only values at the nodes are used - not at the quadpoints

Definition at line 36 of file PorousFlowJoiner.h.

Constructor & Destructor Documentation

◆ PorousFlowJoiner()

PorousFlowJoiner::PorousFlowJoiner ( const InputParameters & parameters )

Definition at line 29 of file PorousFlowJoiner.C.

   : PorousFlowMaterialVectorBase(parameters),
     _pf_prop(getParam<std::string>("material_property")),
     _dporepressure_dvar(!_nodal_material ? getMaterialProperty<std::vector<std::vector<Real>>>(
                                                "dPorousFlow_porepressure_qp_dvar")
                                          : getMaterialProperty<std::vector<std::vector<Real>>>(
                                                "dPorousFlow_porepressure_nodal_dvar")),
     _dsaturation_dvar(!_nodal_material ? getMaterialProperty<std::vector<std::vector<Real>>>(
                                              "dPorousFlow_saturation_qp_dvar")
                                        : getMaterialProperty<std::vector<std::vector<Real>>>(
                                              "dPorousFlow_saturation_nodal_dvar")),
     _dtemperature_dvar(
         !_nodal_material
             ? getMaterialProperty<std::vector<Real>>("dPorousFlow_temperature_qp_dvar")
             : getMaterialProperty<std::vector<Real>>("dPorousFlow_temperature_nodal_dvar")),
     _property(declareProperty<std::vector<Real>>(_pf_prop)),
     _dproperty_dvar(declareProperty<std::vector<std::vector<Real>>>("d" + _pf_prop + "_dvar"))
 {
   _phase_property.resize(_num_phases);
   _dphase_property_dp.resize(_num_phases);
   _dphase_property_ds.resize(_num_phases);
   _dphase_property_dt.resize(_num_phases);
  
   for (unsigned int ph = 0; ph < _num_phases; ++ph)
   {
     std::string phase = Moose::stringify(ph);
     _phase_property[ph] = &getMaterialProperty<Real>(_pf_prop + phase);
     _dphase_property_dp[ph] =
         &getMaterialPropertyDerivative<Real>(_pf_prop + phase, _pressure_variable_name);
     _dphase_property_ds[ph] =
         &getMaterialPropertyDerivative<Real>(_pf_prop + phase, _saturation_variable_name);
     _dphase_property_dt[ph] =
         &getMaterialPropertyDerivative<Real>(_pf_prop + phase, _temperature_variable_name);
   }
 }

Member Function Documentation

◆ computeQpProperties()

void PorousFlowJoiner::computeQpProperties ( )

overrideprotectedvirtual

Definition at line 75 of file PorousFlowJoiner.C.

 {
   initQpStatefulProperties();
  
   _dproperty_dvar[_qp].resize(_num_phases);
   for (unsigned int ph = 0; ph < _num_phases; ++ph)
   {
     _dproperty_dvar[_qp][ph].resize(_num_var);
     for (unsigned v = 0; v < _num_var; ++v)
     {
       // the "if" conditions in the following are because a nodal_material's derivatives might
       // not have been defined.  If that is the case, then DerivativeMaterial passes back a
       // MaterialProperty with zeroes (for the derivatives), but that property will be sized
       // by the number of quadpoints in the element, which may be smaller than the number of
       // nodes!
       _dproperty_dvar[_qp][ph][v] = 0.0;
       if ((*_dphase_property_dp[ph]).size() > _qp)
         _dproperty_dvar[_qp][ph][v] +=
             (*_dphase_property_dp[ph])[_qp] * _dporepressure_dvar[_qp][ph][v];
       if ((*_dphase_property_ds[ph]).size() > _qp)
         _dproperty_dvar[_qp][ph][v] +=
             (*_dphase_property_ds[ph])[_qp] * _dsaturation_dvar[_qp][ph][v];
       if ((*_dphase_property_dt[ph]).size() > _qp)
         _dproperty_dvar[_qp][ph][v] += (*_dphase_property_dt[ph])[_qp] * _dtemperature_dvar[_qp][v];
     }
   }
 }

◆ initQpStatefulProperties()

void PorousFlowJoiner::initQpStatefulProperties ( )

overrideprotectedvirtual

Definition at line 66 of file PorousFlowJoiner.C.

 {
   _property[_qp].resize(_num_phases);
  
   for (unsigned int ph = 0; ph < _num_phases; ++ph)
     _property[_qp][ph] = (*_phase_property[ph])[_qp];
 }