GrainRigidBodyMotionBase.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 "NonlocalKernel.h"
#include "GrainForceAndTorqueInterface.h"

// Forward Declarations
class GrainTrackerInterface;

class GrainRigidBodyMotionBase : public NonlocalKernel

{
public:
  static InputParameters validParams();

  GrainRigidBodyMotionBase(const InputParameters & parameters);

  virtual void timestepSetup();

protected:
  virtual bool globalDoFEnabled(MooseVariableFEBase & /*var*/, dof_id_type /*dof_index*/);

  virtual void precalculateResidual();
  virtual void precalculateJacobian();
  virtual void precalculateOffDiagJacobian(unsigned int jvar);

  virtual void calculateAdvectionVelocity() {}

  const std::vector<dof_id_type> & _var_dofs;

  unsigned int _c_var;
  const VariableValue & _c;
  const VariableGradient & _grad_c;
  const std::vector<dof_id_type> & _c_dofs;

  const unsigned int _op_num;
  const std::vector<const VariableValue *> _vals;
  const std::vector<unsigned int> _vals_var;
  const std::vector<const VariableGradient *> _grad_vals;

  const std::string _base_name;

  const GrainForceAndTorqueInterface & _grain_force_torque;
  const std::vector<RealGradient> & _grain_forces;
  const std::vector<RealGradient> & _grain_torques;
  const std::vector<Real> & _grain_force_c_jacobians;
  const std::vector<std::vector<Real>> & _grain_force_eta_jacobians;

  const Real _mt;
  const Real _mr;

  const GrainTrackerInterface & _grain_tracker;

  const VectorPostprocessorValue & _grain_volumes;

  unsigned int _total_dofs;

  RealGradient _velocity_advection;
  RealGradient _velocity_advection_jacobian;
  std::vector<unsigned int> _grain_ids;
};