jaguar_stripheatsource.i

[Mesh]
  file = 2d_moving.msh
[]

# [MeshModifiers]
#   [./layer_1_sideset]
#     type = SideSetsBetweenSubdomains
#     master_block = layer_1
#     paired_block = layer_2
#     new_boundary = layer_1_top
#   [../]
#   [./layer_2_sideset]
#     type = SideSetsBetweenSubdomains
#     master_block = layer_2
#     paired_block = layer_3
#     new_boundary = layer_2_top
#   [../]
#   [./layer_3_sideset]
#     type = SideSetsBetweenSubdomains
#     master_block = layer_3
#     paired_block = layer_4
#     new_boundary = layer_3_top
#   [../]
  # Define ID for subdomain
  # [./layer_1_id]
  #   type = AssignSubdomainID
  #   subdomain_id = 11
  # [../]
  # [./layer_2_id]
  #   type = AssignSubdomainID
  #   subdomain_id = 10
  # [../]
  # [./layer_3_id]
  #   type = AssignSubdomainID
  #   subdomain_id = 9
  # [../]
  # [./layer_4_id]
  #   type = AssignSubdomainID
  #   subdomain_id = 8
  # [../]
# []

[Variables]
  [./temperature]
    initial_condition = 296 # Initial temperature
    scaling = 1e-1
  [../]
[]

[Functions]
# Assign different thermal conductivity along time to materials in different layers
  # [./material_layer2]
  #   type = ParsedFunction
  #   value = 'if(t>12,6.7,0.0)'
  # [../]
  # [./material_layer3]
  #   type = ParsedFunction
  #   value = 'if(t>8,6.7,0.0)'
  # [../]
  # [./material_layer4]
  #   type = ParsedFunction
  #   value = 'if(t>12,6.7,0.0)'
  # [../]
# Parsed function determines the position of heat point_heat_source
  [./heatsource_layer1]
    type = ParsedFunction
    value = 'if(t=8,1900,296)'
  [../]
[]

[Kernels]
  [./heat_conduction_layer1]
    type = HeatConduction
    # diffusion_coefficient = 6.7 # Thermal conductivity W/m*K * 1e3
    variable = temperature
  [../]
  [./heat_conduction_time_derivative]
    type = HeatConductionTimeDerivative
    # density_name = 0.00443 # g/mm^3
    # specific_heat = 526 # J/g-K * 1e3
    variable = temperature
  [../]
  [./heatsource]
    type = HeatSource
    block = layer_2
    variable = temperature
    function = heatsource_layer1
    # value = 400
  [../]
[]

# [DiracKernels]
#  [./point_heat_source]
#    type = MovingDirac2d
#    variable = temperature
#    value = 200
#    point = '0 0 0'
#  [../]
# []
 # Interface kernel for potential use of heat source
 # [InterfaceKernels]
 #   [./interface_layer1]
 #     type = InterfaceDiffusion
 #     variable = temperature1
 #     neighbor_var = temperature2
 #     boundary = layer_1_top
 #   [../]
 # []

[BCs]
  [./outlet_temperature]
    type = DirichletBC
    variable = temperature
    boundary = bottom_outlet
    value = 296 # (K)
  [../]
#  [./middle]
#   type = MatchedValueBC
#   variable = temperature1
#   boundary = 'layer_1_top'
#   temperature1 = temperature2
# [../]
  # [./layer_1_top_bc]
  #   type = z
  #   variable = temperature
  #   boundary = layer_1_top
  #   function = heatsource_layer1
  # [../]
#   [./layer_2_top_bc]
#     type = DirichletBC
#     variable = temperature
#     boundary = layer_2_top
#     value = 296
#   [../]
#   [./layer_3_top_bc]
#     type = DirichletBC
#     variable = temperature
#     boundary = layer_3_top
#    value = 296
#   [../]
# [./heat_influx]
#    type = NeumannBC
#    variable = temperature
#    boundary =
#    value = 1000
#  [../]
[]

 # [Controls]
 # [./period_0]
 #   type = TimePeriod
 #   disable_objects = 'Kernels::heatsource'
 #   start_time = 8.2
 #   end_time = 16
 #   execute_on = 'initial timestep_begin'
 # [../]
 # [./period_1]
 #   type = TimePeriod
 #   enable_objects = 'Kernels::heatsource'
 #   start_time = 8
 #   end_time = 13
 #   execute_on = 'initial timestep_begin'
 # [../]

  # Disable the layers
  # [./period_0]
  #   type = TimePeriod
  #   disable_objects = 'layer_2'
  #   start_time = 0
  #   end_time = 4
  #   execute_on = 'initial timestep_begin'
  # [../]
  # [./period_1]
  #   type = TimePeriod
  #   enable_objects = 'MeshModifiers::layer_2_id'
  #   disable_objects = 'MeshModifiers::layer_3_id MeshModifiers::layer_4_id'
  #   start_time = 4
  #   end_time = 8
  # [../]
  # [./period_2]
  #   type = TimePeriod
  #   enable_objects = 'MeshModifiers::layer_2_id MeshModifiers::layer_3_id'
  #   disable_objects = 'MeshModifiers::layer_4_id'
  #   start_time = 8
  #   end_time = 12
  # [../]
  # [./period_3]
  #   type = TimePeriod
  #   enable_objects = 'MeshModifiers::layer_2_id MeshModifiers::layer_3_id MeshModifiers::layer_4_id'
  #   start_time = 12
  #   end_time = 16
  # [../]
 # []

[Materials]
  # Update conductivity of materials along simulation time
  [./Ti64_layer1]
    type = GenericConstantMaterial
    # block = layer_1
    prop_names = 'thermal_conductivity specific_heat density'
    prop_values = '0.0067 0.526 0.0043' # W/m*K * 1e3, J/kg-K * 1e3, kg/m^3 @ 296K
  [../]
  # [./Ti64_layer2_conduct]
  #   type = GenericFunctionMaterial
  #   block = layer_2
  #   prop_names = 'thermal_conductivity'
  #   prop_values = material_layer2
  # [../]
  # [./Ti64_layer2]
  #   type = GenericConstantMaterial
  #   block = layer_2
  #   prop_names = 'specific_heat density'
  #   prop_values = '0.526 4430' # W/m*K * 1e3, J/kg-K * 1e3, kg/m^3 @ 296K
  # [../]
  # [./Ti64_layer3_conduct]
  #   type = GenericFunctionMaterial
  #   block = layer_3
  #   prop_names = 'thermal_conductivity'
  #   prop_values = material_layer3
  # [../]
  # [./Ti64_layer3]
  #   type = GenericConstantMaterial
  #   block = layer_3
  #   prop_names = 'specific_heat density'
  #   prop_values = '0.526 4430' # W/m*K * 1e3, J/kg-K * 1e3, kg/m^3 @ 296K
  # [../]
  # [./Ti64_layer4_conduct]
  #   type = GenericFunctionMaterial
  #   block = layer_4
  #   prop_names = 'thermal_conductivity'
  #   prop_values = material_layer4
  # [../]
  # [./Ti64_layer4]
  #   type = GenericConstantMaterial
  #   block = layer_4
  #   prop_names = 'specific_heat density'
  #   prop_values = '0.526 4430' # W/m*K * 1e3, J/kg-K * 1e3, kg/m^3 @ 296K
  # [../]
[]

[Problem]
  type = FEProblem
  coord_type = RZ
  rz_coord_axis = X
[]

[Executioner]
  type = Transient
  nl_rel_tol = 1e-10
  l_tol = 1e-08
  num_steps = 80
  end_time = 16
  solve_type = 'PJFNK'
  petsc_options_iname = '-pc_type -pc_hypre_type'
  petsc_options_value = 'hypre boomeramg'
[]

[Outputs]
  exodus = true
[]