BicrystalCircleGrainICAction

Bicrystal with a circular grain and an embedding outer grain

Overview

This action creates the initial conditions (ICs) for the two order parameters used to represent a bicrystal. It creates a grain structure with the first grain represented by a circle (2D) or sphere (3D) embedded in the second grain.

The action creates two ICs with the SmoothCircleIC, one for each order parameter. Each are defined with the same center coordinates and radius, but for the first order parameter invalue = 1.0 and outvalue = 0.0 and for the second invalue = 0.0 and outvalue = 1.0.

Example Input File Syntax

The BicrystalCircleGrainICAction is accessed through the ICs/PolycrystalICs block, as shown below.

[ICs]
  [./PolycrystalICs]
    [./BicrystalCircleGrainIC]
      radius = 300
      x = 400
      y = 0
    [../]
  [../]
[]

Input Parameters

op_numNumber of grains, should be 2
C++ Type:unsigned int
Controllable:No
Description:Number of grains, should be 2
radiusVoid radius
C++ Type:double
Controllable:No
Description:Void radius
var_name_basespecifies the base name of the variables
C++ Type:std::string
Controllable:No
Description:specifies the base name of the variables
xThe x coordinate of the circle grain center
C++ Type:double
Controllable:No
Description:The x coordinate of the circle grain center
yThe y coordinate of the circle grain center
C++ Type:double
Controllable:No
Description:The y coordinate of the circle grain center

Required Parameters

3D_sphereTruein 3D, whether the smaller grain is a spheres or columnar grain
Default:True
C++ Type:bool
Controllable:No
Description:in 3D, whether the smaller grain is a spheres or columnar grain
active__all__ If specified only the blocks named will be visited and made active
Default:__all__
C++ Type:std::vector<std::string>
Controllable:No
Description:If specified only the blocks named will be visited and made active
blockBlock restriction for the initial condition
C++ Type:std::vector<SubdomainName>
Controllable:No
Description:Block restriction for the initial condition
control_tagsAdds user-defined labels for accessing object parameters via control logic.
C++ Type:std::vector<std::string>
Controllable:No
Description:Adds user-defined labels for accessing object parameters via control logic.
inactiveIf specified blocks matching these identifiers will be skipped.
C++ Type:std::vector<std::string>
Controllable:No
Description:If specified blocks matching these identifiers will be skipped.
int_width0The interfacial width of the void surface. Defaults to sharp interface
Default:0
C++ Type:double
Controllable:No
Description:The interfacial width of the void surface. Defaults to sharp interface
z0The z coordinate of the circle grain center
Default:0
C++ Type:double
Controllable:No
Description:The z coordinate of the circle grain center

Optional Parameters

(../moose/modules/phase_field/test/tests/grain_growth/test.i)

[Mesh]
  type = GeneratedMesh
  dim = 2
  nx = 10
  ny = 10
  nz = 0
  xmin = 0
  xmax = 400
  ymin = 0
  ymax = 400
  zmin = 0
  zmax = 0
  elem_type = QUAD4
  uniform_refine = 1
[]

[GlobalParams]
  op_num = 2
  var_name_base = gr
[]

[Variables]
  [./PolycrystalVariables]
  [../]
[]

[ICs]
  [./PolycrystalICs]
    [./BicrystalCircleGrainIC]
      radius = 300
      x = 400
      y = 0
    [../]
  [../]
[]

[AuxVariables]
  [./bnds]
    order = FIRST
    family = LAGRANGE
  [../]
[]

[Kernels]
  [./PolycrystalKernel]
  [../]
[]

[AuxKernels]
  [./BndsCalc]
    type = BndsCalcAux
    variable = bnds
  [../]
[]

[Materials]
  [./Copper]
    type = GBEvolution
    T = 500 # K
    wGB = 60 # nm
    GBmob0 = 2.5e-6 #m^4/(Js) from Schoenfelder 1997
    Q = 0.23 #Migration energy in eV
    GBenergy = 0.708 #GB energy in J/m^2
  [../]
[]

[Postprocessors]
  [./gr1area]
    type = ElementIntegralVariablePostprocessor
    variable = gr1
  [../]
[]

[Preconditioning]
  [./SMP]
    type = SMP
    full = true
  [../]
[]

[Executioner]
  type = Transient
  scheme = bdf2
  solve_type = 'NEWTON'

  petsc_options_iname = '-pc_type -pc_hypre_type -ksp_gmres_restart'
  petsc_options_value = 'hypre boomeramg 31'

  l_tol = 1.0e-4
  l_max_its = 30
  nl_max_its = 20
  nl_rel_tol = 1.0e-9
  start_time = 0.0
  num_steps = 5
  dt = 80.0

  [./Adaptivity]
    initial_adaptivity = 2
    refine_fraction = 0.8
    coarsen_fraction = 0.05
    max_h_level = 2
  [../]
[]

[Outputs]
  execute_on = 'timestep_end'
  exodus = true
[]