Line data    Source code

       1             : //* This file is part of the MOOSE framework
       2             : //* https://mooseframework.inl.gov
       3             : //*
       4             : //* All rights reserved, see COPYRIGHT for full restrictions
       5             : //* https://github.com/idaholab/moose/blob/master/COPYRIGHT
       6             : //*
       7             : //* Licensed under LGPL 2.1, please see LICENSE for details
       8             : //* https://www.gnu.org/licenses/lgpl-2.1.html
       9             : 
      10             : #include "PolycrystalColoringICAction.h"
      11             : #include "Factory.h"
      12             : #include "FEProblem.h"
      13             : #include "Conversion.h"
      14             : #include "PolycrystalICTools.h"
      15             : 
      16             : registerMooseAction("PhaseFieldApp", PolycrystalColoringICAction, "add_ic");
      17             : 
      18             : InputParameters
      19         446 : PolycrystalColoringICAction::validParams()
      20             : {
      21         446 :   InputParameters params = Action::validParams();
      22         446 :   params.addClassDescription("Action to create ICs for polycrystal variables from a UserObject");
      23         892 :   params.addRequiredParam<unsigned int>("op_num", "number of order parameters to create");
      24         892 :   params.addRequiredParam<std::string>("var_name_base", "specifies the base name of the variables");
      25         892 :   params.addRequiredParam<UserObjectName>("polycrystal_ic_uo", "Optional: TODO");
      26         892 :   params.addParam<std::vector<SubdomainName>>(
      27             :       "block", {}, "Block restriction for the initial condition");
      28         892 :   params.addParam<bool>(
      29         892 :       "linearized_interface", false, "Whether to use linearized interface or the standard model");
      30         892 :   params.addParam<Real>("bound_value",
      31         892 :                         5.0,
      32             :                         "Bound value used to keep variable "
      33             :                         "between +/-bound. Must be positive.");
      34         892 :   params.addParamNamesToGroup("linearized_interface", "LinearizedInterface");
      35         892 :   params.addParamNamesToGroup("bound_value", "LinearizedInterface");
      36             : 
      37         446 :   return params;
      38           0 : }
      39             : 
      40         446 : PolycrystalColoringICAction::PolycrystalColoringICAction(const InputParameters & params)
      41             :   : Action(params),
      42         446 :     _op_num(getParam<unsigned int>("op_num")),
      43         892 :     _var_name_base(getParam<std::string>("var_name_base")),
      44        1338 :     _linearized_interface(getParam<bool>("linearized_interface"))
      45             : {
      46         446 : }
      47             : 
      48             : void
      49         446 : PolycrystalColoringICAction::act()
      50             : {
      51             :   // Loop through the number of order parameters
      52        3354 :   for (unsigned int op = 0; op < _op_num; op++)
      53             :   {
      54        2908 :     std::string IC_type = "PolycrystalColoringIC";
      55        2908 :     if (_linearized_interface)
      56             :       IC_type = "PolycrystalColoringICLinearizedInterface";
      57             : 
      58             :     // Set parameters for IC
      59        2908 :     InputParameters poly_params = _factory.getValidParams(IC_type);
      60        8724 :     poly_params.set<VariableName>("variable") = _var_name_base + Moose::stringify(op);
      61        2908 :     poly_params.set<unsigned int>("op_index") = op;
      62        2908 :     poly_params.applySpecificParameters(parameters(), {"polycrystal_ic_uo", "block"});
      63        5816 :     poly_params.set<UserObjectName>("polycrystal_ic_uo") =
      64        2908 :         getParam<UserObjectName>("polycrystal_ic_uo");
      65        2908 :     if (_linearized_interface)
      66          91 :       poly_params.applySpecificParameters(parameters(),
      67             :                                           {"polycrystal_ic_uo", "block", "bound_value"});
      68             :     else
      69        2817 :       poly_params.applySpecificParameters(parameters(), {"polycrystal_ic_uo", "block"});
      70             : 
      71             :     // Add initial condition
      72        2908 :     _problem->addInitialCondition(IC_type, IC_type + Moose::stringify(op), poly_params);
      73        2908 :   }
      74         446 : }