NodalRotationalInertia.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html

#pragma once

#include "TimeNodalKernel.h"
#include "RankTwoTensor.h"

// Forward Declarations
class TimeIntegrator;
class AuxiliarySystem;

class NodalRotationalInertia : public TimeNodalKernel
{
public:
  static InputParameters validParams();

  NodalRotationalInertia(const InputParameters & parameters);

protected:
  virtual Real computeQpResidual() override;

  virtual Real computeQpJacobian() override;

  virtual Real computeQpOffDiagJacobian(unsigned int jvar) override;

  const bool _has_beta;
  const bool _has_gamma;
  const bool _has_rot_velocities;
  const bool _has_rot_accelerations;
  const bool _has_x_orientation;
  const bool _has_y_orientation;

  AuxiliarySystem * _aux_sys;

  unsigned int _nrot;

  std::vector<const VariableValue *> _rot;

  std::vector<const VariableValue *> _rot_old;

  std::vector<unsigned int> _rot_vel_num;

  std::vector<unsigned int> _rot_accel_num;

  std::vector<unsigned int> _rot_variables;

  std::vector<Real> _rot_accel;

  std::vector<Real> _rot_vel;

  std::vector<Real> _rot_vel_old;

  const Real _beta;

  const Real _gamma;

  const Real & _eta;

  const Real & _alpha;

  const unsigned int _component;

  RankTwoTensor _inertia;

  std::vector<const VariableValue *> _rot_dot_residual;

  std::vector<const VariableValue *> _rot_vel_old_value;

  std::vector<const VariableValue *> _rot_dotdot_residual;

  const VariableValue * _du_dot_du;

  const VariableValue * _du_dotdot_du;

  const TimeIntegrator & _time_integrator;
};