#include <HHPFCRFFSplitKernelAction.h>

Inheritance diagram for HHPFCRFFSplitKernelAction:

Public Member Functions
	HHPFCRFFSplitKernelAction (const InputParameters &params)

virtual void	act ()

Protected Attributes
const unsigned int	_num_L

const std::string	_L_name_base

const VariableName	_n_name

Detailed Description

Definition at line 20 of file HHPFCRFFSplitKernelAction.h.

Constructor & Destructor Documentation

◆ HHPFCRFFSplitKernelAction()

HHPFCRFFSplitKernelAction::HHPFCRFFSplitKernelAction ( const InputParameters & params )

Definition at line 40 of file HHPFCRFFSplitKernelAction.C.

   : Action(params),
     _num_L(getParam<unsigned int>("num_L")),
     _L_name_base(getParam<std::string>("L_name_base")),
     _n_name(getParam<VariableName>("n_name"))
 {
 }

Member Function Documentation

◆ act()

void HHPFCRFFSplitKernelAction::act ( )

virtual

Reimplemented in PFCRFFKernelAction.

Definition at line 49 of file HHPFCRFFSplitKernelAction.C.

 {
   // Loop over the L_variables
   for (unsigned int l = 0; l < _num_L; ++l)
   {
     // Create L base name
     std::string L_name = _L_name_base + Moose::stringify(l);
  
     // Create real  and imaginary L variable names
     std::string real_name = L_name + "_real";
     std::string imag_name = L_name + "_imag";
  
     //
     // Create the kernels for the real L variable
     //
  
     // Create the diffusion kernel for L_real_l
     InputParameters poly_params = _factory.getValidParams("Diffusion");
     poly_params.set<NonlinearVariableName>("variable") = real_name;
     poly_params.set<bool>("use_displaced_mesh") = getParam<bool>("use_displaced_mesh");
     _problem->addKernel("Diffusion", "diff_" + real_name, poly_params);
  
     // Create the (alpha^R_m L^R_m) term
     poly_params = _factory.getValidParams("HHPFCRFF");
     poly_params.set<NonlinearVariableName>("variable") = real_name;
     poly_params.set<bool>("positive") = true;
     poly_params.set<MaterialPropertyName>("prop_name") = "alpha_R_" + Moose::stringify(l);
     poly_params.set<bool>("use_displaced_mesh") = getParam<bool>("use_displaced_mesh");
     _problem->addKernel("HHPFCRFF", "HH1_" + real_name, poly_params);
  
     // **Create the -(alpha^I_m L^I_m) term
     if (l > 0)
     {
       poly_params = _factory.getValidParams("HHPFCRFF");
       poly_params.set<NonlinearVariableName>("variable") = real_name;
       poly_params.set<bool>("positive") = false;
       poly_params.set<std::vector<VariableName>>("coupled_var").push_back(imag_name);
       poly_params.set<bool>("use_displaced_mesh") = getParam<bool>("use_displaced_mesh");
       poly_params.set<MaterialPropertyName>("prop_name") = "alpha_I_" + Moose::stringify(l);
       _problem->addKernel("HHPFCRFF", "HH2_" + real_name, poly_params);
     }
  
     // **Create the -(A^R_m n) term
     poly_params = _factory.getValidParams("HHPFCRFF");
     poly_params.set<NonlinearVariableName>("variable") = real_name;
     poly_params.set<bool>("positive") = false;
     poly_params.set<std::vector<VariableName>>("coupled_var").push_back(_n_name);
     poly_params.set<bool>("use_displaced_mesh") = getParam<bool>("use_displaced_mesh");
     poly_params.set<MaterialPropertyName>("prop_name") = "A_R_" + Moose::stringify(l);
     _problem->addKernel("HHPFCRFF", "HH3_" + real_name, poly_params);
  
     //
     // Create the kernels for the imaginary L variable, l > 0
     //
     if (l > 0)
     {
       // Create the diffusion kernel for L_imag_l
       InputParameters poly_params = _factory.getValidParams("Diffusion");
       poly_params.set<NonlinearVariableName>("variable") = imag_name;
       poly_params.set<bool>("use_displaced_mesh") = getParam<bool>("use_displaced_mesh");
       _problem->addKernel("Diffusion", "diff_" + imag_name, poly_params);
  
       // **Create the (alpha^R_m L^I_m) term
       poly_params = _factory.getValidParams("HHPFCRFF");
       poly_params.set<NonlinearVariableName>("variable") = imag_name;
       poly_params.set<bool>("positive") = true;
       poly_params.set<bool>("use_displaced_mesh") = getParam<bool>("use_displaced_mesh");
       poly_params.set<MaterialPropertyName>("prop_name") = "alpha_R_" + Moose::stringify(l);
       _problem->addKernel("HHPFCRFF", "HH1_" + imag_name, poly_params);
  
       // **Create the (alpha^I_m L^R_m) term
       poly_params = _factory.getValidParams("HHPFCRFF");
       poly_params.set<NonlinearVariableName>("variable") = imag_name;
       poly_params.set<bool>("positive") = true;
       poly_params.set<std::vector<VariableName>>("coupled_var").push_back(real_name);
       poly_params.set<bool>("use_displaced_mesh") = getParam<bool>("use_displaced_mesh");
       poly_params.set<MaterialPropertyName>("prop_name") = "alpha_I_" + Moose::stringify(l);
       _problem->addKernel("HHPFCRFF", "HH2_" + imag_name, poly_params);
  
       // **Create the -(A^I_m n) term
       poly_params = _factory.getValidParams("HHPFCRFF");
       poly_params.set<NonlinearVariableName>("variable") = imag_name;
       poly_params.set<bool>("positive") = false;
       poly_params.set<std::vector<VariableName>>("coupled_var").push_back(_n_name);
       poly_params.set<bool>("use_displaced_mesh") = getParam<bool>("use_displaced_mesh");
       poly_params.set<MaterialPropertyName>("prop_name") = "A_I_" + Moose::stringify(l);
       _problem->addKernel("HHPFCRFF", "HH3_" + imag_name, poly_params);
     }
   }
 }