geometry
{
  import_mesh false ; Use built-in mesh generator
  dim 3 ; dimension of the domain
  length 36.0e-3 ; [m]
  height 6.0e-3 ; [m] In 3D, the third parameters is width
  width  6.0e-3
  length_divisions 24 ; Number of cell layers in the length direction
  height_divisions 3 ; Number of cell layers in the height direction
  width_divisions  3
  material_deposition true
  material_deposition_method scan_paths
  deposition_length 0.3e-3
  deposition_height 2.0e-3
  deposition_width 2.0e-3
  deposition_lead_time 1.0e0
  material_height 6.0e-3
}

physics
{
  thermal true
  mechanical true
}

discretization
{
  thermal 
  {
    fe_degree 3
    quadrature gauss ; Optional parameter. Possibilities: gauss or lobatto
  }
  mechanical 
  {
    fe_degree 1
  }
}

boundary
{
  type adiabatic ; convective,radiative
  boundary_4
  {
    type adiabatic,clamped
  }
}

refinement
{
  n_refinements 1 ; Number of time the cells on the paths of the beams are
                  ; refined
  time_steps_between_refinement 100 ; number of time steps after which
                                    ; the refinement process is performed
}

materials
{
  n_materials 1

  property_format polynomial
  material_0
  {
    solid
    {
      density 7904; [kg/m^3] For now all the states needs to have the same
                    ; density.
      specific_heat 714; [J/kg K]
      thermal_conductivity_x 31.4 ; [W/m K]
      thermal_conductivity_y 31.4 ; [W/m K]
      thermal_conductivity_z 31.4 ; [W/m K]
      lame_first_parameter 92.0e9 ; [Pa]
      lame_second_parameter 79.0e9 ; [Pa]
      thermal_expansion_coef 17.2e-6 ; [1/K]
    }

    powder
    {
      specific_heat 714; [J/kg K]
      density 7904; [kg/m^3]
      thermal_conductivity_x 0.314 ; [W/m K]
      thermal_conductivity_y 0.314 ; [W/m K]
      thermal_conductivity_z 0.314 ; [W/m K]
    }

    liquid
    {
      specific_heat 847; [J/kg K]
      density 7904; [kg/m^3]
      thermal_conductivity_x 37.3 ; [W/m k]
      thermal_conductivity_y 37.3 ; [W/m k]
      thermal_conductivity_z 37.3 ; [W/m k]

      ; Not all three states need to define the same properties or to exist
    }

    solidus 1675; [K]
    liquidus 1708; [K]
    latent_heat 290000 ; [J/kg]
    
  }
}

sources
{
  n_beams 1

  beam_0
  {
    type goldak ; goldak (laser) or electron_beam
    depth 3.0e-3 ; [m] maximum depth reached by the laser
    diameter 3.0e-3 ; [m]
    scan_path_file thermoelastic_bare_plate_scan_path.txt
    scan_path_file_format segment
    absorption_efficiency 0.6 ; number between 0 and 1 equivalent to
                              ; energy_conversion_efficiency * control_efficiency
                              ; for an electron beam
    max_power 9000.0 ; [W], current * voltage for an electron beam
  }
}

time_stepping
{
  method forward_euler ; Possibilities: forward_euler, rk_third_order,
                        ; rk_fourth_order
  duration 4.0e0 ; [s]
  time_step 1.0e-3 ; [s]
}

post_processor
{
  filename_prefix output
  time_steps_between_output 500
}


profiling
{
  timer false
  caliper "spot(profile.mpi),loop-report,runtime-report"
}

memory_space host ; Always run on the host