Esempio n. 1
0
PresetVelocity::PresetVelocity(const std::string & name, InputParameters parameters) :
  PresetNodalBC(name, parameters),
  _u_old(valueOld()),
  _velocity(parameters.get<Real>("velocity")),
  _function(parameters.get<FunctionName>("function") != "" ? &getFunction("function") : NULL)
{
}
Esempio n. 2
0
PresetAcceleration::PresetAcceleration(const InputParameters & parameters)
  : PresetNodalBC(parameters),
    _u_old(valueOld()),
    _scale_factor(parameters.get<Real>("scale_factor")),
    _function(getFunction("function")),
    _vel_old(coupledValueOld("velocity")),
    _accel_old(coupledValueOld("acceleration")),
    _beta(getParam<Real>("beta"))
{
}
Esempio n. 3
0
CHBulk::CHBulk(const std::string & name, InputParameters parameters)
  :KernelGrad(name, parameters),
   _mob_name(getParam<std::string>("mob_name")),
   _Dmob_name(getParam<std::string>("Dmob_name")),
   _M(getMaterialProperty<Real>(_mob_name)),
   _implicit(getParam<bool>("implicit")),
   _u_old(valueOld()),
   _grad_u_old(_implicit ? _grad_zero : gradientOld()),
   _has_MJac(getParam<bool>("has_MJac")),
   _DM(_has_MJac ? &getMaterialProperty<Real>(_Dmob_name) : NULL)
{
}
Esempio n. 4
0
PoroFullSatTimeDerivative::PoroFullSatTimeDerivative(const InputParameters & parameters) :
    DerivativeMaterialInterface<TimeDerivative>(parameters),
    _u_old(valueOld()),
    _volstrain(getMaterialProperty<Real>("volumetric_strain")),
    _volstrain_old(getMaterialPropertyOld<Real>("volumetric_strain")),

    _disp_x_var(coupled("disp_x")),
    _disp_y_var(coupled("disp_y")),
    _disp_z_var(coupled("disp_z")),

    _alpha(getMaterialProperty<Real>("biot_coefficient")),

    _one_over_biot_modulus(getMaterialProperty<Real>("one_over_biot_modulus")),
    _done_over_biot_modulus_dP(getMaterialPropertyDerivative<Real>("one_over_biot_modulus", _var.name())),
    _done_over_biot_modulus_dep(getMaterialPropertyDerivative<Real>("one_over_biot_modulus", "volumetric_strain"))

{
}
Esempio n. 5
0
PoroFullSatTimeDerivative::PoroFullSatTimeDerivative(const InputParameters & parameters)
  : DerivativeMaterialInterface<TimeDerivative>(parameters),
    _u_old(valueOld()),
    _volstrain(getMaterialProperty<Real>("volumetric_strain")),
    _volstrain_old(getMaterialPropertyOld<Real>("volumetric_strain")),

    _ndisp(coupledComponents("displacements")),
    _disp_var_num(_ndisp),

    _alpha(getMaterialProperty<Real>("biot_coefficient")),

    _one_over_biot_modulus(getMaterialProperty<Real>("one_over_biot_modulus")),
    _done_over_biot_modulus_dP(
        getMaterialPropertyDerivative<Real>("one_over_biot_modulus", _var.name())),
    _done_over_biot_modulus_dep(
        getMaterialPropertyDerivative<Real>("one_over_biot_modulus", "volumetric_strain"))
{
  for (unsigned i = 0; i < _ndisp; ++i)
    _disp_var_num[i] = coupled("displacements", i);
}
Esempio n. 6
0
CoupledBEEquilibriumSub::CoupledBEEquilibriumSub(const InputParameters & parameters) :
    Kernel(parameters),
    _weight(getParam<Real>("weight")),
    _log_k(getParam<Real>("log_k")),
    _sto_u(getParam<Real>("sto_u")),
    _sto_v(getParam<std::vector<Real> >("sto_v")),
    _porosity(getMaterialProperty<Real>("porosity")),
    _u_old(valueOld())
{
  const unsigned int n = coupledComponents("v");
  _vars.resize(n);
  _v_vals.resize(n);
  _v_vals_old.resize(n);

  for (unsigned int i=0; i < n; ++i)
  {
    _vars[i] = coupled("v", i);
    _v_vals[i] = &coupledValue("v", i);
    _v_vals_old[i] = & coupledValueOld("v", i);
  }
}
Esempio n. 7
0
EnthalpyTimeDerivative::EnthalpyTimeDerivative(const std::string & name,
                                                     InputParameters parameters):TimeDerivative(name, parameters),
_density(getMaterialProperty<Real>("density")),
_time_old_density(getMaterialProperty<Real>("time_old_density")),

_temperature(getMaterialProperty<Real>("material_temperature")),
_time_old_temperature(getMaterialProperty<Real>("time_old_material_temperature")),

_dTdH_P(getMaterialProperty<Real>("dTdH_P")),
_dTdP_H(getMaterialProperty<Real>("dTdP_H")),
_ddensitydH_P(getMaterialProperty<Real>("ddensitydH_P")),
_ddensitydp_H(getMaterialProperty<Real>("ddensitydp_H")),

_p_var(coupled("pressure")),

_porosity(getMaterialProperty<Real>("porosity")),
_specific_heat_rock(getMaterialProperty<Real>("specific_heat_rock")),
_density_rock(getMaterialProperty<Real>("density_rock")),
_u_old(valueOld())

{}
Esempio n. 8
0
SecondDerivative::SecondDerivative(const std::string & name, InputParameters parameters) :
    TimeKernel(name, parameters),
    _u_old(valueOld()),
    _u_older(valueOlder())
{}
// DEPRECATED CONSTRUCTOR
SecondDerivativeImplicitEuler::SecondDerivativeImplicitEuler(const std::string & deprecated_name, InputParameters parameters) :
    TimeKernel(deprecated_name, parameters),
    _u_old(valueOld()),
    _u_older(valueOlder())
{}
SecondDerivativeImplicitEuler::SecondDerivativeImplicitEuler(const InputParameters & parameters) :
    TimeKernel(parameters),
    _u_old(valueOld()),
    _u_older(valueOlder())
{}