NekMeshDeformation

Description

NekMeshDeformation is a FieldTransfer that couples NekRS and MOOSE via fluid-structure interaction. This is an upcoming capability and is not yet documented

Input Parameters

directionDirection in which to send data
C++ Type:MooseEnum
Options:to_nek, from_nek
Controllable:No
Description:Direction in which to send data

Required Parameters

usrwrk_slotWhen 'direction = to_nek', the slot(s) in the usrwrk array to write the incoming data; provide one entry for each quantity being passed
C++ Type:std::vector<unsigned int>
Controllable:No
Description:When 'direction = to_nek', the slot(s) in the usrwrk array to write the incoming data; provide one entry for each quantity being passed

Optional Parameters

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.
enableTrueSet the enabled status of the MooseObject.
Default:True
C++ Type:bool
Controllable:No
Description:Set the enabled status of the MooseObject.

Advanced Parameters

Input Files

(test/tests/deformation/vol-areas/nek.i)
(test/tests/nek_errors/deformation/nek.i)
(test/tests/nek_errors/deformation/user/nek.i)
(test/tests/deformation/simple-cube/nek.i)
(test/tests/deformation/mesh-velocity-areas/nek.i)

(test/tests/deformation/vol-areas/nek.i)

[Mesh]
  type = NekRSMesh
  order = SECOND
  volume = true
  parallel_type = replicated
  displacements = 'disp_x disp_y disp_z'
[]

[Problem]
  type = NekRSProblem
  casename = 'nekbox'

  [FieldTransfers]
    [heat_source]
      type = NekVolumetricSource
      direction = to_nek
      usrwrk_slot = 0
      postprocessor_to_conserve = source_integral
    []
    [disp]
      type = NekMeshDeformation
      usrwrk_slot = '1 2 3'
      direction = to_nek
    []
    [temp]
      type = NekFieldVariable
      field = temperature
      direction = from_nek
    []
  []
[]

[Executioner]
  type = Transient
  [TimeStepper]
    type = NekTimeStepper
  []
  [Quadrature]
    type = GAUSS_LOBATTO
    order = SECOND
  []
[]

[Postprocessors]
  [nekbdry_icar1]
    type = AreaPostprocessor
    boundary = '1'
    use_displaced_mesh = true
    execute_on = INITIAL
  []
  [nekbdry_ar1]
    type = AreaPostprocessor
    boundary = '1'
    use_displaced_mesh = true
  []
  [nekbdry_ar2]
    type = AreaPostprocessor
    boundary = '2'
    use_displaced_mesh = true
  []
  [nekbdry_ar3]
    type = AreaPostprocessor
    boundary = '3'
    use_displaced_mesh = true
  []
  [nekbdry_ar4]
    type = AreaPostprocessor
    boundary = '4'
    use_displaced_mesh = true
  []
  [nekbdry_ar5]
    type = AreaPostprocessor
    boundary = '5'
    use_displaced_mesh = true
  []
  [nekbdry_ar6]
    type = AreaPostprocessor
    boundary = '6'
    use_displaced_mesh = true
  []
[]

[Outputs]
  csv = true
  execute_on = 'final'

  # uncomment the temp_ansol to see that the solution matches very well
  hide = 'source_integral heat_source'
[]

(test/tests/nek_errors/deformation/nek.i)

[Mesh]
  type = NekRSMesh
  order = SECOND
  volume = true
  parallel_type = replicated
[]

[Problem]
  type = NekRSProblem
  casename = 'elast_nomv'

  [FieldTransfers]
    [disp]
      type = NekMeshDeformation
      direction = to_nek
      usrwrk_slot = '1 2 3'
    []
  []
[]

[Executioner]
  type = Transient

  [TimeStepper]
    type = NekTimeStepper
  []
[]

(test/tests/nek_errors/deformation/user/nek.i)

[Mesh]
  type = NekRSMesh
  volume = true
  boundary = '2'
  parallel_type = replicated
[]

[Problem]
  type = NekRSProblem
  casename = 'user'

  [FieldTransfers]
    [disp]
      type = NekMeshDeformation
      direction = to_nek
      usrwrk_slot = '0 1 2'
    []
    [temp]
      type = NekFieldVariable
      field = temperature
      direction = from_nek
    []
  []
[]

[Executioner]
  type = Transient
  [TimeStepper]
    type = NekTimeStepper
  []
[]

(test/tests/deformation/simple-cube/nek.i)

[Mesh]
  type = NekRSMesh
  order = SECOND
  volume = true
  parallel_type = replicated
  displacements = 'disp_x disp_y disp_z'
[]

[Problem]
  type = NekRSProblem
  casename = 'nekbox'
  synchronization_interval = parent_app

  [FieldTransfers]
    [heat_source]
      type = NekVolumetricSource
      direction = to_nek
      usrwrk_slot = 0
      postprocessor_to_conserve = source_integral
    []
    [disp]
      type = NekMeshDeformation
      usrwrk_slot = '1 2 3'
      direction = to_nek
    []
    [temp]
      type = NekFieldVariable
      field = temperature
      direction = from_nek
    []
  []
[]

[AuxVariables]
  [temp_ansol]
    order = SECOND
  []
[]

[Executioner]
  type = Transient
  [TimeStepper]
    type = NekTimeStepper
  []
[]

[Postprocessors]
  [integral]
    type = ElementL2Difference
    variable = temp
    other_variable = temp_ansol
    execute_on = timestep_end
  []
  [max_T]
    type = NekVolumeExtremeValue
    field = temperature
    value_type = max
    execute_on = timestep_end
  []
  [min_T]
    type = NekVolumeExtremeValue
    field = temperature
    value_type = min
    execute_on = timestep_end
  []
[]

[Outputs]
  exodus = true
  execute_on = 'final'

  # uncomment the temp_ansol to see that the solution matches very well
  hide = 'heat_source transfer_in source_integral temp_ansol'
[]

(test/tests/deformation/mesh-velocity-areas/nek.i)

[Mesh]
  type = NekRSMesh
  order = SECOND
  boundary = '2'
  parallel_type = replicated
  displacements = 'disp_x disp_y disp_z'
  use_displaced_mesh =true
[]

[Problem]
  type = NekRSProblem
  casename = 'pipe'

  [FieldTransfers]
    [disp]
      type = NekMeshDeformation
      direction = to_nek
      usrwrk_slot = '1 2 3'
    []
    [temp]
      type = NekFieldVariable
      field = temperature
      direction = from_nek
    []
  []
[]

[Executioner]
  type = Transient

  [TimeStepper]
    type = NekTimeStepper
  []

  [Quadrature]
    type = GAUSS_LOBATTO
    order = SECOND
  []
[]

[Postprocessors]
  [nek_area_initial]
    type = NekSideIntegral
    field = unity
    boundary = '2'
    execute_on = initial
  []
  [nek_area]
    type = NekSideIntegral
    field = unity
    boundary = '2'
  []
  [nek_diff]
    type = DifferencePostprocessor
    value1 = nek_area
    value2 = nek_area_initial
  []
  [moose_area_initial]
    type = VolumePostprocessor # for a side-only mesh, this is area
    execute_on = initial
    use_displaced_mesh = true
  []
  [moose_area]
    type = VolumePostprocessor
    use_displaced_mesh = true
  []
  [moose_diff]
    type = DifferencePostprocessor
    value1 = moose_area
    value2 = moose_area_initial
  []
[]

[Outputs]
  csv = true
  execute_on = 'final'
  hide = 'nek_area_initial nek_area moose_area_initial moose_area'
[]

Description
Input Parameters
Input Files