#include <InertialForce.h>

Inheritance diagram for InertialForce:

Public Member Functions
	InertialForce (const InputParameters &parameters)

Protected Member Functions
virtual Real	computeQpResidual ()

virtual Real	computeQpJacobian ()

Private Attributes
const MaterialProperty< Real > &	_density

const VariableValue *	_u_old

const VariableValue *	_vel_old

const VariableValue *	_accel_old

const bool	_has_beta

const bool	_has_gamma

const Real	_beta

const Real	_gamma

const bool	_has_velocity

const bool	_has_acceleration

const MaterialProperty< Real > &	_eta

const Real	_alpha

const VariableValue *	_u_dot

const VariableValue *	_u_dotdot

const VariableValue *	_u_dot_old

const VariableValue *	_du_dot_du

const VariableValue *	_du_dotdot_du

Detailed Description

Definition at line 22 of file InertialForce.h.

Constructor & Destructor Documentation

◆ InertialForce()

InertialForce::InertialForce ( const InputParameters & parameters )

Definition at line 44 of file InertialForce.C.

   : TimeKernel(parameters),
     _density(getMaterialProperty<Real>("density")),
     _has_beta(isParamValid("beta")),
     _has_gamma(isParamValid("gamma")),
     _beta(_has_beta ? getParam<Real>("beta") : 0.1),
     _gamma(_has_gamma ? getParam<Real>("gamma") : 0.1),
     _has_velocity(isParamValid("velocity")),
     _has_acceleration(isParamValid("acceleration")),
     _eta(getMaterialProperty<Real>("eta")),
     _alpha(getParam<Real>("alpha"))
 {
   if (_has_beta && _has_gamma && _has_velocity && _has_acceleration)
   {
     _vel_old = &coupledValueOld("velocity");
     _accel_old = &coupledValueOld("acceleration");
     _u_old = &valueOld();
   }
   else if (!_has_beta && !_has_gamma && !_has_velocity && !_has_acceleration)
   {
     _u_dot = &(_var.uDot());
     _u_dotdot = &(_var.uDotDot());
     _u_dot_old = &(_var.uDotOld());
     _du_dot_du = &(_var.duDotDu());
     _du_dotdot_du = &(_var.duDotDotDu());
   }
   else
     mooseError("InertialForce: Either all or none of `beta`, `gamma`, `velocity`and `acceleration` "
                "should be provided as input.");
 }

Member Function Documentation

◆ computeQpJacobian()

Real InertialForce::computeQpJacobian ( )

protectedvirtual

Definition at line 96 of file InertialForce.C.

 {
   if (_dt == 0)
     return 0;
   else if (_has_beta)
     return _test[_i][_qp] * _density[_qp] / (_beta * _dt * _dt) * _phi[_j][_qp] +
            _eta[_qp] * (1 + _alpha) * _test[_i][_qp] * _density[_qp] * _gamma / _beta / _dt *
                _phi[_j][_qp];
   else
     return _test[_i][_qp] * _density[_qp] * (*_du_dotdot_du)[_qp] * _phi[_j][_qp] +
            _eta[_qp] * (1 + _alpha) * _test[_i][_qp] * _density[_qp] * (*_du_dot_du)[_qp] *
                _phi[_j][_qp];
 }

◆ computeQpResidual()

Real InertialForce::computeQpResidual ( )

protectedvirtual

Definition at line 76 of file InertialForce.C.

 {
   if (_dt == 0)
     return 0;
   else if (_has_beta)
   {
     Real accel = 1. / _beta *
                  (((_u[_qp] - (*_u_old)[_qp]) / (_dt * _dt)) - (*_vel_old)[_qp] / _dt -
                   (*_accel_old)[_qp] * (0.5 - _beta));
     Real vel = (*_vel_old)[_qp] + (_dt * (1 - _gamma)) * (*_accel_old)[_qp] + _gamma * _dt * accel;
     return _test[_i][_qp] * _density[_qp] *
            (accel + vel * _eta[_qp] * (1 + _alpha) - _alpha * _eta[_qp] * (*_vel_old)[_qp]);
   }
   else
     return _test[_i][_qp] * _density[_qp] *
            ((*_u_dotdot)[_qp] + (*_u_dot)[_qp] * _eta[_qp] * (1.0 + _alpha) -
             _alpha * _eta[_qp] * (*_u_dot_old)[_qp]);
 }