Importing a constant number (maybe via ConstantFunction) to damage dependent perm (K) calculation

[Traiwit]

Hi guys,

so here is my current version of PorousFlowPermeabilityConstCSV is a bit messy and have constant numbers (1e-13 1e-15 etc) which I have to update every time when I work on different projects.

is there a way that I can define those number in the input file and call those numbers in my PorousFlowPermeabilityConstCSV?

I see that ConstantFunction allows me to define a constant number, but how to I take the value from there to my perm calculation?

Please help if you have any idea

Thank you!
Traiwit

@GiudGiud @cpgr tagging you guys, coz you guys are very good at this.

PorousFlowPermeabilityConstCSV::PorousFlowPermeabilityConstCSV(const InputParameters & parameters)
  : PorousFlowPermeabilityBase(parameters), _damage(getMaterialProperty<Real>("damage")),
  _read_prop_user_object(isParamValid("read_prop_user_object")
                             ? &getUserObject<ElementPropertyReadFile>("read_prop_user_object")
                             : nullptr),
  _mat_prop(declareProperty<RealVectorValue>("perm")),
  _RMD(declareProperty<RealVectorValue>("RMD")),
  _kmax(declareProperty<RealVectorValue>("kmax")),
  _perm(declareProperty<RealVectorValue>("perm_tensor")),
  _flag(getMaterialProperty<Real>("flag"))
{


PorousFlowPermeabilityConstCSV::computeQpProperties()
{

 
  //  read values from .CSV file
    _mat_prop[_qp](0) = _read_prop_user_object->getData(_current_elem, 0);
    _mat_prop[_qp](1) = _read_prop_user_object->getData(_current_elem, 1);
    _mat_prop[_qp](2) = _read_prop_user_object->getData(_current_elem, 2);
    _RMD[_qp](0) = _read_prop_user_object->getData(_current_elem, 6);
    _kmax[_qp](0)    = _read_prop_user_object->getData(_current_elem, 7);
    // // calculate perm tensor

        if (_damage[_qp] >= 3 ) { // backfill  
        _perm[_qp](0) =  1e-15; 
        _perm[_qp](1) =  1e-15; 
        _perm[_qp](2) =  1e-15; 

       }

        else if (_damage[_qp] >= 2 && _damage[_qp] < 3 ) {  // CAVE 

        _perm[_qp](0) =  1e-13; 
        _perm[_qp](1) =  1e-13; 
        _perm[_qp](2) =  1e-13; 

       }

    else {

    _perm[_qp](0) = std::min(exp(_damage[_qp]*(1/(_RMD[_qp](0))*log(_kmax[_qp](0)/_mat_prop[_qp](0))))*_mat_prop[_qp](0)*1e-7, _kmax[_qp](0)*1e-7); 
    _perm[_qp](1) = std::min(exp(_damage[_qp]*(1/(_RMD[_qp](0))*log(_kmax[_qp](0)/_mat_prop[_qp](1))))*_mat_prop[_qp](1)*1e-7, _kmax[_qp](0)*1e-7); 
    _perm[_qp](2) = std::min(exp(_damage[_qp]*(1/(_RMD[_qp](0))*log(_kmax[_qp](0)/_mat_prop[_qp](2))))*_mat_prop[_qp](2)*1e-7, _kmax[_qp](0)*1e-7); 

       }

[cpgr]

Have you tried adding these as parameters? Something like (I'm just guessing names for these, you probably have better ones you could use)

params.addParam<Real>("perm_backfill", 1e-15, "Permeability for backfill (default 1e-15)");
params.addParam<Real>("perm_cave", 1e-13, "Permeability for cave (default 1e-13)");

in validParams(), then

_perm_backfill(getParam<Real>("perm_backfill")

etc in the constructor, then

if (_damage[_qp] >= 3 ) { // backfill  
        _perm[_qp](0) =  _perm_backfill; 
        _perm[_qp](1) =   _perm_backfill; 
        _perm[_qp](2) =   _perm_backfill; 

       }

etc

[Traiwit]

omg I feel so stupid after reading your comment

Thank you very much @cpgr