Automatically generates all variables, Kernels, and Materials to ensure the correct derivatives of the elastic free energy in a non-split Cahn-Hilliard simulation are assembled. More...

#include <DisplacementGradientsAction.h>

Inheritance diagram for DisplacementGradientsAction:

[legend]

Public Member Functions
	DisplacementGradientsAction (const InputParameters &params)

virtual void	act ()

Private Attributes
std::vector< VariableName >	_displacements

std::vector< VariableName >	_displacement_gradients

Detailed Description

Automatically generates all variables, Kernels, and Materials to ensure the correct derivatives of the elastic free energy in a non-split Cahn-Hilliard simulation are assembled.

Definition at line 20 of file DisplacementGradientsAction.h.

Constructor & Destructor Documentation

◆ DisplacementGradientsAction()

DisplacementGradientsAction::DisplacementGradientsAction ( const InputParameters & params )

Definition at line 39 of file DisplacementGradientsAction.C.

   : Action(params),
     _displacements(getParam<std::vector<VariableName>>("displacements")),
     _displacement_gradients(getParam<std::vector<VariableName>>("displacement_gradients"))
 {
 }

Member Function Documentation

◆ act()

void DisplacementGradientsAction::act ( )

virtual

Definition at line 47 of file DisplacementGradientsAction.C.

 {
   unsigned int ngrad = _displacement_gradients.size();
  
   if (_current_task == "add_variable")
   {
     // Loop through the gij variables
     Real scaling = getParam<Real>("scaling");
     for (unsigned int i = 0; i < ngrad; ++i)
     {
       auto var_params = _factory.getValidParams("MooseVariable");
       var_params.set<MooseEnum>("family") = "LAGRANGE";
       var_params.set<MooseEnum>("order") = "FIRST";
       var_params.set<std::vector<Real>>("scaling") = {scaling};
  
       // Create displacement gradient variables
       _problem->addVariable("MooseVariable", _displacement_gradients[i], var_params);
     }
   }
   else if (_current_task == "add_material")
   {
     InputParameters params = _factory.getValidParams("StrainGradDispDerivatives");
     params.set<std::vector<VariableName>>("displacement_gradients") = _displacement_gradients;
     params.set<std::vector<SubdomainName>>("block") = {"0"}; // TODO: add parameter for this
     _problem->addMaterial("StrainGradDispDerivatives", "strain_grad_disp_derivatives", params);
   }
   else if (_current_task == "add_kernel")
   {
     unsigned int ndisp = _displacements.size();
     if (ndisp * ndisp != ngrad)
       paramError("displacement_gradients",
                  "Number of displacement gradient variables must be the square of the number of "
                  "displacement variables.");
  
     // Loop through the displacements
     unsigned int i = 0;
     for (unsigned int j = 0; j < ndisp; ++j)
       for (unsigned int k = 0; k < ndisp; ++k)
       {
         InputParameters params = _factory.getValidParams("GradientComponent");
         params.set<NonlinearVariableName>("variable") = _displacement_gradients[i];
         params.set<std::vector<VariableName>>("v") = {_displacements[j]};
         params.set<unsigned int>("component") = k;
         _problem->addKernel(
             "GradientComponent", _displacement_gradients[i] + "_grad_kernel", params);
         ++i;
       }
   }
   else
     mooseError("Internal error.");
 }