NekMassFluxWeightedSideAverage

Compute mass flux weighted average of a field over a boundary in the NekRS mesh

Description

This postprocessor computes the mass-flux-weighted average of a specified field over a boundary in the NekRS mesh,

$var element = document.getElementById("moose-equation-bc84a574-d2df-4b7c-9c52-89b6912a403a");katex.render("p=\\frac{\\int_{\\Gamma}f\\rho\\vec{V}\\cdot\\hat{n}\\ d\\Gamma}{\\int_{\\Gamma}\\rho\\vec{V}\\cdot\\hat{n}\\ d\\Gamma}", element, {displayMode:true,throwOnError:false});$

where $var element = document.getElementById("moose-equation-40028f05-5d84-4cf8-8e40-dc1256e9c19b");katex.render("p", element, {displayMode:false,throwOnError:false});$ is the value of the postprocessor, $var element = document.getElementById("moose-equation-e0947885-6b8d-497f-ae95-0ff33a90e05d");katex.render("\\Gamma", element, {displayMode:false,throwOnError:false});$ is the boundary of the NekRS mesh, $var element = document.getElementById("moose-equation-7bc09057-28a6-4c87-be4a-c733e47274b6");katex.render("\\rho", element, {displayMode:false,throwOnError:false});$ is the fluid density, $var element = document.getElementById("moose-equation-ba944b97-c2cd-42f0-b9e6-94838f88026f");katex.render("\\vec{V}", element, {displayMode:false,throwOnError:false});$ is the fluid velocity, $var element = document.getElementById("moose-equation-6cd3cf35-dc58-4694-8100-4a0ef1a796d8");katex.render("f", element, {displayMode:false,throwOnError:false});$ is the specified field, and $var element = document.getElementById("moose-equation-e456c07b-3a6c-41b1-9dda-ed76e688ea15");katex.render("\\hat{n}", element, {displayMode:false,throwOnError:false});$ is the surface unit normal.

The boundaries over which to compute this postprocessor in the NekRS mesh are specified with the boundary parameter; these boundaries are the sidesets in NekRS's mesh (i.e. the .re2 file). You can specify more than one boundary, with syntax such as boundary = '1 2'. To be clear, this postprocessor is not evaluated on the NekRSMesh mesh mirror, but instead on the mesh actually used for computation in NekRS.

The field is specified with the field parameter, which may be one of:

pressure
temperature
velocity (magnitude of velocity)
velocity_x ( $var element = document.getElementById("moose-equation-54330620-dc8d-4c74-b663-7acc5dc318e9");katex.render("x", element, {displayMode:false,throwOnError:false});$ -component of velocity)
velocity_y ( $var element = document.getElementById("moose-equation-95d094ce-104d-42c9-9fe6-6f6f5257ce02");katex.render("y", element, {displayMode:false,throwOnError:false});$ -component of velocity)
velocity_z ( $var element = document.getElementById("moose-equation-90b9a7c8-a697-49ea-b15d-670256594240");katex.render("z", element, {displayMode:false,throwOnError:false});$ -component of velocity)
velocity_component (velocity vector projected onto another vector)
scalar01
scalar02
scalar03
unity

Setting field = unity is equivalent to computing 1, since the numerator will be exactly equal to the denominator (unity is of more use for other postprocessors).

If running NekRS in non-dimensional form (and you have indicated the appropriate nondimensional scales by setting nondimensional = true for the [Problem], then the value of this postprocessor is shown in dimensional units. On the otherhand, NekRS scalars (scalar01, scalar02 or scalar03) are never dimensionalized because their dimensions are problem-dependent.

Example Input Syntax

As an example, the outlet_T postprocessor evaluates the mass-flux-weighted average of temperature over the outlet boundary of NekRS's mesh.

[Postprocessors]
  [outlet_T]
    type = NekMassFluxWeightedSideAverage
    field = temperature
    boundary = '3'
  []
  [max_T]
    type = NekVolumeExtremeValue
    field = temperature
  []
[]

Input Parameters

boundaryBoundary ID(s) for which to compute the postprocessor
C++ Type:std::vector<int>
Controllable:No
Description:Boundary ID(s) for which to compute the postprocessor
fieldField to integrate
C++ Type:MooseEnum
Options:velocity_component, velocity_x, velocity_y, velocity_z, velocity, temperature, pressure, scalar01, scalar02, scalar03, unity
Controllable:No
Description:Field to integrate

Required Parameters

execute_onTIMESTEP_ENDThe list of flag(s) indicating when this object should be executed, the available options include NONE, INITIAL, LINEAR, NONLINEAR, POSTCHECK, TIMESTEP_END, TIMESTEP_BEGIN, MULTIAPP_FIXED_POINT_END, MULTIAPP_FIXED_POINT_BEGIN, FINAL, CUSTOM.
Default:TIMESTEP_END
C++ Type:ExecFlagEnum
Options:NONE, INITIAL, LINEAR, NONLINEAR, POSTCHECK, TIMESTEP_END, TIMESTEP_BEGIN, MULTIAPP_FIXED_POINT_END, MULTIAPP_FIXED_POINT_BEGIN, FINAL, CUSTOM, TRANSFER
Controllable:No
Description:The list of flag(s) indicating when this object should be executed, the available options include NONE, INITIAL, LINEAR, NONLINEAR, POSTCHECK, TIMESTEP_END, TIMESTEP_BEGIN, MULTIAPP_FIXED_POINT_END, MULTIAPP_FIXED_POINT_BEGIN, FINAL, CUSTOM.
meshallNekRS mesh to compute postprocessor on
Default:all
C++ Type:MooseEnum
Options:fluid, solid, all
Controllable:No
Description:NekRS mesh to compute postprocessor on
prop_getter_suffixAn optional suffix parameter that can be appended to any attempt to retrieve/get material properties. The suffix will be prepended with a '_' character.
C++ Type:MaterialPropertyName
Controllable:No
Description:An optional suffix parameter that can be appended to any attempt to retrieve/get material properties. The suffix will be prepended with a '_' character.
use_interpolated_stateFalseFor the old and older state use projected material properties interpolated at the quadrature points. To set up projection use the ProjectedStatefulMaterialStorageAction.
Default:False
C++ Type:bool
Controllable:No
Description:For the old and older state use projected material properties interpolated at the quadrature points. To set up projection use the ProjectedStatefulMaterialStorageAction.
velocity_directionDirection in which to evaluate velocity, for 'field = velocity_component'. For example, velocity_direction = '1 0 0' will get the x-component of velocity.
C++ Type:libMesh::Point
Controllable:No
Description:Direction in which to evaluate velocity, for 'field = velocity_component'. For example, velocity_direction = '1 0 0' will get the x-component of velocity.

Optional Parameters

allow_duplicate_execution_on_initialFalseIn the case where this UserObject is depended upon by an initial condition, allow it to be executed twice during the initial setup (once before the IC and again after mesh adaptivity (if applicable).
Default:False
C++ Type:bool
Controllable:No
Description:In the case where this UserObject is depended upon by an initial condition, allow it to be executed twice during the initial setup (once before the IC and again after mesh adaptivity (if applicable).
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:Yes
Description:Set the enabled status of the MooseObject.
execution_order_group0Execution order groups are executed in increasing order (e.g., the lowest number is executed first). Note that negative group numbers may be used to execute groups before the default (0) group. Please refer to the user object documentation for ordering of user object execution within a group.
Default:0
C++ Type:int
Controllable:No
Description:Execution order groups are executed in increasing order (e.g., the lowest number is executed first). Note that negative group numbers may be used to execute groups before the default (0) group. Please refer to the user object documentation for ordering of user object execution within a group.
force_postauxFalseForces the UserObject to be executed in POSTAUX
Default:False
C++ Type:bool
Controllable:No
Description:Forces the UserObject to be executed in POSTAUX
force_preauxFalseForces the UserObject to be executed in PREAUX
Default:False
C++ Type:bool
Controllable:No
Description:Forces the UserObject to be executed in PREAUX
force_preicFalseForces the UserObject to be executed in PREIC during initial setup
Default:False
C++ Type:bool
Controllable:No
Description:Forces the UserObject to be executed in PREIC during initial setup
outputsVector of output names where you would like to restrict the output of variables(s) associated with this object
C++ Type:std::vector<OutputName>
Controllable:No
Description:Vector of output names where you would like to restrict the output of variables(s) associated with this object
use_displaced_meshFalseWhether or not this object should use the displaced mesh for computation. Note that in the case this is true but no displacements are provided in the Mesh block the undisplaced mesh will still be used.
Default:False
C++ Type:bool
Controllable:No
Description:Whether or not this object should use the displaced mesh for computation. Note that in the case this is true but no displacements are provided in the Mesh block the undisplaced mesh will still be used.

Advanced Parameters

Input Files

(test/tests/postprocessors/nek_weighted_side_average/nek.i)
(test/tests/cht/sfr_pincell/nek_vpp.i)
(tutorials/pincell_multiphysics/nek.i)
(test/tests/cht/nondimensional/nek_exact.i)
(test/tests/nek_standalone/conj_ht/nek.i)
(test/tests/cht/nondimensional/nek.i)
(test/tests/cht/sfr_pincell/nek.i)

(tutorials/pincell_multiphysics/nek.i)

inlet_T  = 573.0           # inlet temperature
power = 250                # total power (W)
Re = 500.0                 # Reynolds number

pin_diameter = 0.97e-2     # pin outer diameter
pin_pitch = 1.28e-2        # pin pitch

mu = 8.8e-5                # fluid dynamic viscosity
rho = 723.6                # fluid density
Cp = 5512.0                # fluid isobaric specific heat capacity

flow_area = ${fparse pin_pitch * pin_pitch - pi * pin_diameter * pin_diameter / 4.0}
wetted_perimeter = ${fparse pi * pin_diameter}
hydraulic_diameter = ${fparse 4.0 * flow_area / wetted_perimeter}

U_ref = ${fparse Re * mu / rho / hydraulic_diameter}
mdot = ${fparse rho * U_ref * flow_area}
dT = ${fparse power / mdot / Cp}

[Mesh]
  type = NekRSMesh
  boundary = '1'
  scaling = ${hydraulic_diameter}
  volume = true
[]

[Problem]
  type = NekRSProblem
  casename = 'fluid'

  nondimensional = true
  L_ref = ${hydraulic_diameter}
  T_ref = ${inlet_T}
  U_ref = ${U_ref}
  dT_ref = ${dT}

  rho_0 = ${rho}
  Cp_0 = ${Cp}

  has_heat_source = false
  synchronization_interval = parent_app
[]

[Postprocessors]
  [outlet_T]
    type = NekMassFluxWeightedSideAverage
    field = temperature
    boundary = '3'
  []
  [max_T]
    type = NekVolumeExtremeValue
    field = temperature
  []
[]

[Executioner]
  type = Transient

  [TimeStepper]
    type = NekTimeStepper
  []
[]

[Outputs]
  exodus = true
  interval = 10
[]

(test/tests/postprocessors/nek_weighted_side_average/nek.i)

[Problem]
  type = NekRSProblem
  casename = 'brick'
[]

[Mesh]
  type = NekRSMesh
  boundary = '1 2 3 4 5 6'
[]

[Executioner]
  type = Transient

  [TimeStepper]
    type = NekTimeStepper
  []
[]

[Outputs]
  [out]
    type = CSV
    hide = 'flux_integral'
    execute_on = 'final'
  []
[]

[Postprocessors]
  [mdot_side1]
    type = NekMassFluxWeightedSideIntegral
    field = unity
    boundary = '1'
  []
  [mdot_side2]
    type = NekMassFluxWeightedSideIntegral
    field = unity
    boundary = '2'
  []
  [mdot_side3]
    type = NekMassFluxWeightedSideIntegral
    field = unity
    boundary = '3'
  []
  [mdot_side4]
    type = NekMassFluxWeightedSideIntegral
    field = unity
    boundary = '4'
  []
  [mdot_side5]
    type = NekMassFluxWeightedSideIntegral
    field = unity
    boundary = '5'
  []
  [mdot_side6]
    type = NekMassFluxWeightedSideIntegral
    field = unity
    boundary = '6'
  []
  [weighted_T_side1]
    type = NekMassFluxWeightedSideAverage
    field = temperature
    boundary = '1'
  []
  [weighted_T_side2]
    type = NekMassFluxWeightedSideAverage
    field = temperature
    boundary = '2'
  []
  [weighted_T_side3]
    type = NekMassFluxWeightedSideAverage
    field = temperature
    boundary = '3'
  []
  [weighted_T_side4]
    type = NekMassFluxWeightedSideAverage
    field = temperature
    boundary = '4'
  []
  [weighted_T_side5]
    type = NekMassFluxWeightedSideAverage
    field = temperature
    boundary = '5'
  []
  [weighted_T_side6]
    type = NekMassFluxWeightedSideAverage
    field = temperature
    boundary = '6'
  []
[]

(test/tests/cht/sfr_pincell/nek_vpp.i)

[Problem]
  type = NekRSProblem
  casename = 'sfr_pin'
  conserve_flux_by_sideset = true

  # we only technically need one scratch space slot for this problem,
  # so we can skip allocating extra
  n_usrwrk_slots = 1
[]

[Mesh]
  type = NekRSMesh
  boundary = '1'
[]

[Executioner]
  type = Transient

  [TimeStepper]
    type = NekTimeStepper
  []
[]

[Postprocessors]
  [nek_flux]
    type = NekHeatFluxIntegral
    boundary = '1'
  []
  [average_inlet_T]
    type = NekSideAverage
    field = temperature
    boundary = '3'
    execute_on = initial
  []
  [average_outlet_T]
    type = NekSideAverage
    field = temperature
    boundary = '4'
  []
  [dT]
    type = DifferencePostprocessor
    value1 = average_outlet_T
    value2 = average_inlet_T
  []
  [inlet_mdot]
    type = NekMassFluxWeightedSideIntegral
    field = unity
    boundary = '3'
    execute_on = initial
  []

  # postprocessors for comparing against non-dimensional version in ../nondimensional
  # --> uncomment in order to get the reference values that the nondimensional boundary coupling
  #     was verified against
  #
  # # side integral
  # [area_1]
  #   type = NekSideIntegral
  #   field = unity
  #   boundary = '1'
  # []
  # [pressure_1]
  #   type = NekSideIntegral
  #   field = pressure
  #   boundary = '1'
  # []
  # [temperature_1]
  #   type = NekSideIntegral
  #   field = temperature
  #   boundary = '1'
  # []

  # # side average
  # [avg_area_1]
  #   type = NekSideAverage
  #   field = unity
  #   boundary = '1'
  # []
  # [avg_pressure_1]
  #   type = NekSideAverage
  #   field = pressure
  #   boundary = '1'
  # []
  # [avg_temperature_1]
  #   type = NekSideAverage
  #   field = temperature
  #   boundary = '1'
  # []

  # # volume integral
  # [volume]
  #   type = NekVolumeIntegral
  #   field = unity
  # []
  # [pressure_vol]
  #   type = NekVolumeIntegral
  #   field = pressure
  # []
  # [temperature_vol]
  #   type = NekVolumeIntegral
  #   field = temperature
  # []

  # # volume average
  # [avg_volume]
  #   type = NekVolumeAverage
  #   field = unity
  # []
  # [avg_pressure_vol]
  #   type = NekVolumeAverage
  #   field = pressure
  # []
  # [avg_temperature_vol]
  #   type = NekVolumeAverage
  #   field = temperature
  # []

  # # heat flux integral
  # [nek_flux]
  #   type = NekHeatFluxIntegral
  #   boundary = '1'
  # []

  # # mass flux weighted integral
  # [inlet_mdot]
  #   type = NekMassFluxWeightedSideIntegral
  #   field = unity
  #   boundary = '3'
  #   execute_on = initial
  # []
  # [outlet_T]
  #   type = NekMassFluxWeightedSideIntegral
  #   field = temperature
  #   boundary = '4'
  # []
  # [inlet_P]
  #   type = NekMassFluxWeightedSideIntegral
  #   field = pressure
  #   boundary = '4'
  # []

  # # mass flux weighted integral
  # [inlet_mdot_avg]
  #   type = NekMassFluxWeightedSideAverage
  #   field = unity
  #   boundary = '3'
  #   execute_on = initial
  # []
  # [outlet_T_avg]
  #   type = NekMassFluxWeightedSideAverage
  #   field = temperature
  #   boundary = '4'
  # []
  # [inlet_P_avg]
  #   type = NekMassFluxWeightedSideAverage
  #   field = pressure
  #   boundary = '4'
  # []

  # # extreme value postprocessors - VOLUME
  # [max_T]
  #   type = NekVolumeExtremeValue
  #   field = temperature
  #   value_type = max
  # []
  # [min_T]
  #   type = NekVolumeExtremeValue
  #   field = temperature
  #   value_type = min
  # []
  # [max_p]
  #   type = NekVolumeExtremeValue
  #   field = pressure
  #   value_type = max
  # []
  # [min_p]
  #   type = NekVolumeExtremeValue
  #   field = pressure
  #   value_type = min
  # []
  # [max_1]
  #   type = NekVolumeExtremeValue
  #   field = unity
  #   value_type = max
  # []
  # [min_1]
  #   type = NekVolumeExtremeValue
  #   field = unity
  #   value_type = min
  # []

  # # extreme value postprocessors - SIDE
  # [max_T_out]
  #   type = NekSideExtremeValue
  #   field = temperature
  #   boundary = '4'
  #   value_type = max
  # []
  # [min_T_out]
  #   type = NekSideExtremeValue
  #   field = temperature
  #   boundary = '4'
  #   value_type = min
  # []
  # [max_p_in]
  #   type = NekSideExtremeValue
  #   field = pressure
  #   boundary = '3'
  #   value_type = max
  # []
  # [min_p_in]
  #   type = NekSideExtremeValue
  #   field = pressure
  #   boundary = '3'
  #   value_type = min
  # []
  # [max_1_in]
  #   type = NekSideExtremeValue
  #   field = unity
  #   boundary = '3'
  #   value_type = max
  # []
  # [min_1_in]
  #   type = NekSideExtremeValue
  #   field = unity
  #   boundary = '3'
  #   value_type = min
  # []
[]

[Outputs]
  exodus = true
  execute_on = 'final'
  csv = true

  [screen]
    type = Console
    hide = 'average_inlet_T average_outlet_T'
  []
[]

(tutorials/pincell_multiphysics/nek.i)

inlet_T  = 573.0           # inlet temperature
power = 250                # total power (W)
Re = 500.0                 # Reynolds number

pin_diameter = 0.97e-2     # pin outer diameter
pin_pitch = 1.28e-2        # pin pitch

mu = 8.8e-5                # fluid dynamic viscosity
rho = 723.6                # fluid density
Cp = 5512.0                # fluid isobaric specific heat capacity

flow_area = ${fparse pin_pitch * pin_pitch - pi * pin_diameter * pin_diameter / 4.0}
wetted_perimeter = ${fparse pi * pin_diameter}
hydraulic_diameter = ${fparse 4.0 * flow_area / wetted_perimeter}

U_ref = ${fparse Re * mu / rho / hydraulic_diameter}
mdot = ${fparse rho * U_ref * flow_area}
dT = ${fparse power / mdot / Cp}

[Mesh]
  type = NekRSMesh
  boundary = '1'
  scaling = ${hydraulic_diameter}
  volume = true
[]

[Problem]
  type = NekRSProblem
  casename = 'fluid'

  nondimensional = true
  L_ref = ${hydraulic_diameter}
  T_ref = ${inlet_T}
  U_ref = ${U_ref}
  dT_ref = ${dT}

  rho_0 = ${rho}
  Cp_0 = ${Cp}

  has_heat_source = false
  synchronization_interval = parent_app
[]

[Postprocessors]
  [outlet_T]
    type = NekMassFluxWeightedSideAverage
    field = temperature
    boundary = '3'
  []
  [max_T]
    type = NekVolumeExtremeValue
    field = temperature
  []
[]

[Executioner]
  type = Transient

  [TimeStepper]
    type = NekTimeStepper
  []
[]

[Outputs]
  exodus = true
  interval = 10
[]

(test/tests/cht/nondimensional/nek_exact.i)

[Problem]
  type = NekRSProblem
  casename = 'sfr_pin'

  # This input is run in nondimensional form to verify that all the postprocessors
  # and data transfers in/out of nekRS are properly dimensionalized.
  nondimensional = true
  U_ref = 0.0950466
  T_ref = 628.15
  dT_ref = 50.0
  L_ref = 0.908e-2
  rho_0 = 834.5
  Cp_0 = 1228.0
[]

[Mesh]
  type = NekRSMesh
  boundary = '1'
  exact = true

  # nekRS runs in non-dimensional form, which means that we shrunk the mesh
  # from physical units of meters to our characteristic scale of 0.908e-2 m
  # (the pin pitch, arbitrarily chosen). That means that we must multiply
  # the nekRS mesh by 0.908e-2 to get back in units of meters that BISON is
  # running in.
  scaling = 0.908e-2
[]

[Executioner]
  type = Transient

  [TimeStepper]
    type = NekTimeStepper
  []
[]

[Postprocessors]
# pressure is not a reliable quantity to compare in a short, few-timestep regression test.
# This line can be uncommented for comparing with the dimensional version.
inactive = 'avg_pressure_1 avg_pressure_vol max_p max_p_in min_p_in pressure_1 pressure_vol'

  [synchronization_in]
    type = Receiver
  []

  # side integral
  [area_1]
    type = NekSideIntegral
    field = unity
    boundary = '1'
  []
  [pressure_1]
    type = NekSideIntegral
    field = pressure
    boundary = '1'
  []
  [temperature_1]
    type = NekSideIntegral
    field = temperature
    boundary = '1'
  []

  # side average
  [avg_area_1]
    type = NekSideAverage
    field = unity
    boundary = '1'
  []
  [avg_pressure_1]
    type = NekSideAverage
    field = pressure
    boundary = '1'
  []
  [avg_temperature_1]
    type = NekSideAverage
    field = temperature
    boundary = '1'
  []

  # volume integral
  [volume]
    type = NekVolumeIntegral
    field = unity
  []
  [pressure_vol]
    type = NekVolumeIntegral
    field = pressure
  []
  [temperature_vol]
    type = NekVolumeIntegral
    field = temperature
  []

  # volume average
  [avg_volume]
    type = NekVolumeAverage
    field = unity
  []
  [avg_pressure_vol]
    type = NekVolumeAverage
    field = pressure
  []
  [avg_temperature_vol]
    type = NekVolumeAverage
    field = temperature
  []

  # heat flux integral
  [nek_flux]
    type = NekHeatFluxIntegral
    boundary = '1'
  []

  # mass flux weighted integral
  [inlet_mdot]
    type = NekMassFluxWeightedSideIntegral
    field = unity
    boundary = '3'
    execute_on = initial
  []
  [outlet_T]
    type = NekMassFluxWeightedSideIntegral
    field = temperature
    boundary = '4'
  []
  [inlet_P]
    type = NekMassFluxWeightedSideIntegral
    field = pressure
    boundary = '4'
  []

  # mass flux weighted integral
  [inlet_mdot_avg]
    type = NekMassFluxWeightedSideAverage
    field = unity
    boundary = '3'
    execute_on = initial
  []
  [outlet_T_avg]
    type = NekMassFluxWeightedSideAverage
    field = temperature
    boundary = '4'
  []
  [inlet_P_avg]
    type = NekMassFluxWeightedSideAverage
    field = pressure
    boundary = '4'
  []

  # extreme value postprocessors - VOLUME
  [max_T]
    type = NekVolumeExtremeValue
    field = temperature
    value_type = max
  []
  [min_T]
    type = NekVolumeExtremeValue
    field = temperature
    value_type = min
  []
  [max_p]
    type = NekVolumeExtremeValue
    field = pressure
    value_type = max
  []
  [min_p]
    type = NekVolumeExtremeValue
    field = pressure
    value_type = min
  []
  [max_1]
    type = NekVolumeExtremeValue
    field = unity
    value_type = max
  []
  [min_1]
    type = NekVolumeExtremeValue
    field = unity
    value_type = min
  []

  # extreme value postprocessors - SIDE
  [max_T_out]
    type = NekSideExtremeValue
    field = temperature
    boundary = '4'
    value_type = max
  []
  [min_T_out]
    type = NekSideExtremeValue
    field = temperature
    boundary = '4'
    value_type = min
  []
  [max_p_in]
    type = NekSideExtremeValue
    field = pressure
    boundary = '3'
    value_type = max
  []
  [min_p_in]
    type = NekSideExtremeValue
    field = pressure
    boundary = '3'
    value_type = min
  []
  [max_1_in]
    type = NekSideExtremeValue
    field = unity
    boundary = '3'
    value_type = max
  []
  [min_1_in]
    type = NekSideExtremeValue
    field = unity
    boundary = '3'
    value_type = min
  []
[]

[Outputs]
  exodus = true
  execute_on = 'final'
  csv = true

  [screen]
    type = Console
    hide = 'synchronization_in'
  []
[]

(test/tests/nek_standalone/conj_ht/nek.i)

[Mesh]
  type = NekRSMesh
  volume = true
[]

[Problem]
  type = NekRSStandaloneProblem
  casename = 'conj_ht'
  output = 'temperature'
[]

[Executioner]
  type = Transient

  [TimeStepper]
    type = NekTimeStepper
  []
[]

[Postprocessors]
 [Area_BC3_flow]
   type = NekSideIntegral
   boundary = '3'
   field = unity
   mesh = 'fluid'
 []
 [Area_BC3_all]
   type = NekSideIntegral
   boundary = '3'
   field = unity
   mesh = 'all'
 []
 [SideAverage_T_BC3_flow]
   type = NekSideAverage
   boundary = '3'
   field = temperature
   mesh = 'fluid'
 []
 [SideAverage_T_BC3_all]
   type = NekSideAverage
   boundary = '3'
   field = temperature
   mesh = 'all'
 []
 [HeatFluxIntegral_BC3_flow]
   type = NekHeatFluxIntegral
   boundary = '3'
   mesh = 'fluid'
 []
 [HeatFluxIntegral_BC3_all]
   type = NekHeatFluxIntegral
   boundary = '3'
   mesh = 'all'
 []
 [MassFlowRate_BC1_flow]
   type = NekMassFluxWeightedSideIntegral
   boundary = '1'
   field = unity
   mesh = 'fluid'
 []
 [MassFlowRate_BC1_all]
   type = NekMassFluxWeightedSideIntegral
   boundary = '1'
   field = unity
   mesh = 'all'
 []
 [MflowAvgTemp_BC2_flow]
   type = NekMassFluxWeightedSideAverage
   boundary = '2'
   field = temperature
   mesh = 'fluid'
 []
 [MflowAvgTemp_BC2_all]
   type = NekMassFluxWeightedSideAverage
   boundary = '2'
   field = temperature
   mesh = 'all'
 []
 [Reynolds_BC1_flow]
   type = ReynoldsNumber
   boundary = '1'
   L_ref = 0.5
   mesh = 'fluid'
 []
 [Reynolds_BC1_all]
   type = ReynoldsNumber
   boundary = '1'
   L_ref = 0.5
   mesh = 'all'
 []
 [Peclet_BC1_flow]
   type = PecletNumber
   boundary = '1'
   L_ref = 0.5
   mesh = 'fluid'
 []
 [Peclet_BC1_all]
   type = PecletNumber
   boundary = '1'
   L_ref = 0.5
   mesh = 'all'
 []
#
# Volume post processors
#
 [Vol_flow]
   type = NekVolumeIntegral
   field = unity
   mesh = 'fluid'
 []
 [Vol_all]
   type = NekVolumeIntegral
   field = unity
   mesh = 'all'
 []
 [maxVol_T_flow]
   type = NekVolumeExtremeValue
   field = temperature
   value_type = max
   mesh = 'fluid'
 []
 [maxVol_T_all]
   type = NekVolumeExtremeValue
   field = temperature
   value_type = max
   mesh = 'all'
 []
 [avgVol_T_flow]
   type = NekVolumeAverage
   field = temperature
   mesh = 'fluid'
 []
 [avgVol_T_all]
   type = NekVolumeAverage
   field = temperature
   mesh = 'all'
 []
[]

[Outputs]
  csv = true
  exodus = true
  execute_on = final
[]

(test/tests/cht/nondimensional/nek.i)

[Problem]
  type = NekRSProblem
  casename = 'sfr_pin'

  # This input is run in nondimensional form to verify that all the postprocessors
  # and data transfers in/out of nekRS are properly dimensionalized.
  nondimensional = true
  U_ref = 0.0950466
  T_ref = 628.15
  dT_ref = 50.0
  L_ref = 0.908e-2
  rho_0 = 834.5
  Cp_0 = 1228.0
[]

[Mesh]
  type = NekRSMesh
  boundary = '1'

  # nekRS runs in non-dimensional form, which means that we shrunk the mesh
  # from physical units of meters to our characteristic scale of 0.908e-2 m
  # (the pin pitch, arbitrarily chosen). That means that we must multiply
  # the nekRS mesh by 0.908e-2 to get back in units of meters that BISON is
  # running in.
  scaling = 0.908e-2
[]

[Executioner]
  type = Transient

  [TimeStepper]
    type = NekTimeStepper
  []
[]

[Postprocessors]
  [synchronization_in]
    type = Receiver
  []

  # side integral
  [area_1]
    type = NekSideIntegral
    field = unity
    boundary = '1'
  []
  [pressure_1]
    type = NekSideIntegral
    field = pressure
    boundary = '1'
  []
  [temperature_1]
    type = NekSideIntegral
    field = temperature
    boundary = '1'
  []

  # side average
  [avg_area_1]
    type = NekSideAverage
    field = unity
    boundary = '1'
  []
  [avg_pressure_1]
    type = NekSideAverage
    field = pressure
    boundary = '1'
  []
  [avg_temperature_1]
    type = NekSideAverage
    field = temperature
    boundary = '1'
  []

  # volume integral
  [volume]
    type = NekVolumeIntegral
    field = unity
  []
  [pressure_vol]
    type = NekVolumeIntegral
    field = pressure
  []
  [temperature_vol]
    type = NekVolumeIntegral
    field = temperature
  []

  # volume average
  [avg_volume]
    type = NekVolumeAverage
    field = unity
  []
  [avg_pressure_vol]
    type = NekVolumeAverage
    field = pressure
  []
  [avg_temperature_vol]
    type = NekVolumeAverage
    field = temperature
  []

  # heat flux integral
  [nek_flux]
    type = NekHeatFluxIntegral
    boundary = '1'
  []

  # mass flux weighted integral
  [inlet_mdot]
    type = NekMassFluxWeightedSideIntegral
    field = unity
    boundary = '3'
    execute_on = initial
  []
  [outlet_T]
    type = NekMassFluxWeightedSideIntegral
    field = temperature
    boundary = '4'
  []
  [inlet_P]
    type = NekMassFluxWeightedSideIntegral
    field = pressure
    boundary = '4'
  []

  # mass flux weighted integral
  [inlet_mdot_avg]
    type = NekMassFluxWeightedSideAverage
    field = unity
    boundary = '3'
    execute_on = initial
  []
  [outlet_T_avg]
    type = NekMassFluxWeightedSideAverage
    field = temperature
    boundary = '4'
  []
  [inlet_P_avg]
    type = NekMassFluxWeightedSideAverage
    field = pressure
    boundary = '4'
  []

  # extreme value postprocessors - VOLUME
  [max_T]
    type = NekVolumeExtremeValue
    field = temperature
    value_type = max
  []
  [min_T]
    type = NekVolumeExtremeValue
    field = temperature
    value_type = min
  []
  [max_p]
    type = NekVolumeExtremeValue
    field = pressure
    value_type = max
  []
  [min_p]
    type = NekVolumeExtremeValue
    field = pressure
    value_type = min
  []
  [max_1]
    type = NekVolumeExtremeValue
    field = unity
    value_type = max
  []
  [min_1]
    type = NekVolumeExtremeValue
    field = unity
    value_type = min
  []

  # extreme value postprocessors - SIDE
  [max_T_out]
    type = NekSideExtremeValue
    field = temperature
    boundary = '4'
    value_type = max
  []
  [min_T_out]
    type = NekSideExtremeValue
    field = temperature
    boundary = '4'
    value_type = min
  []
  [max_p_in]
    type = NekSideExtremeValue
    field = pressure
    boundary = '3'
    value_type = max
  []
  [min_p_in]
    type = NekSideExtremeValue
    field = pressure
    boundary = '3'
    value_type = min
  []
  [max_1_in]
    type = NekSideExtremeValue
    field = unity
    boundary = '3'
    value_type = max
  []
  [min_1_in]
    type = NekSideExtremeValue
    field = unity
    boundary = '3'
    value_type = min
  []
[]

[Outputs]
  exodus = true
  execute_on = 'final'

  [screen]
    type = Console
    hide = 'synchronization_in'
  []
[]

(test/tests/cht/sfr_pincell/nek.i)

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

  # we only technically need one scratch space slot for this problem,
  # so we can skip allocating extra
  n_usrwrk_slots = 1
[]

[Mesh]
  type = NekRSMesh
  boundary = '1'
[]

[Executioner]
  type = Transient

  [TimeStepper]
    type = NekTimeStepper
  []
[]

[Postprocessors]
  [nek_flux]
    type = NekHeatFluxIntegral
    boundary = '1'
  []
  [average_inlet_T]
    type = NekSideAverage
    field = temperature
    boundary = '3'
    execute_on = initial
  []
  [average_outlet_T]
    type = NekSideAverage
    field = temperature
    boundary = '4'
  []
  [dT]
    type = DifferencePostprocessor
    value1 = average_outlet_T
    value2 = average_inlet_T
  []
  [inlet_mdot]
    type = NekMassFluxWeightedSideIntegral
    field = unity
    boundary = '3'
    execute_on = initial
  []

  # postprocessors for comparing against non-dimensional version in ../nondimensional
  # --> uncomment in order to get the reference values that the nondimensional boundary coupling
  #     was verified against
  #
  # # side integral
  # [area_1]
  #   type = NekSideIntegral
  #   field = unity
  #   boundary = '1'
  # []
  # [pressure_1]
  #   type = NekSideIntegral
  #   field = pressure
  #   boundary = '1'
  # []
  # [temperature_1]
  #   type = NekSideIntegral
  #   field = temperature
  #   boundary = '1'
  # []

  # # side average
  # [avg_area_1]
  #   type = NekSideAverage
  #   field = unity
  #   boundary = '1'
  # []
  # [avg_pressure_1]
  #   type = NekSideAverage
  #   field = pressure
  #   boundary = '1'
  # []
  # [avg_temperature_1]
  #   type = NekSideAverage
  #   field = temperature
  #   boundary = '1'
  # []

  # # volume integral
  # [volume]
  #   type = NekVolumeIntegral
  #   field = unity
  # []
  # [pressure_vol]
  #   type = NekVolumeIntegral
  #   field = pressure
  # []
  # [temperature_vol]
  #   type = NekVolumeIntegral
  #   field = temperature
  # []

  # # volume average
  # [avg_volume]
  #   type = NekVolumeAverage
  #   field = unity
  # []
  # [avg_pressure_vol]
  #   type = NekVolumeAverage
  #   field = pressure
  # []
  # [avg_temperature_vol]
  #   type = NekVolumeAverage
  #   field = temperature
  # []

  # # heat flux integral
  # [nek_flux]
  #   type = NekHeatFluxIntegral
  #   boundary = '1'
  # []

  # # mass flux weighted integral
  # [inlet_mdot]
  #   type = NekMassFluxWeightedSideIntegral
  #   field = unity
  #   boundary = '3'
  #   execute_on = initial
  # []
  # [outlet_T]
  #   type = NekMassFluxWeightedSideIntegral
  #   field = temperature
  #   boundary = '4'
  # []
  # [inlet_P]
  #   type = NekMassFluxWeightedSideIntegral
  #   field = pressure
  #   boundary = '4'
  # []

  # # mass flux weighted integral
  # [inlet_mdot_avg]
  #   type = NekMassFluxWeightedSideAverage
  #   field = unity
  #   boundary = '3'
  #   execute_on = initial
  # []
  # [outlet_T_avg]
  #   type = NekMassFluxWeightedSideAverage
  #   field = temperature
  #   boundary = '4'
  # []
  # [inlet_P_avg]
  #   type = NekMassFluxWeightedSideAverage
  #   field = pressure
  #   boundary = '4'
  # []

  # # extreme value postprocessors - VOLUME
  # [max_T]
  #   type = NekVolumeExtremeValue
  #   field = temperature
  #   value_type = max
  # []
  # [min_T]
  #   type = NekVolumeExtremeValue
  #   field = temperature
  #   value_type = min
  # []
  # [max_p]
  #   type = NekVolumeExtremeValue
  #   field = pressure
  #   value_type = max
  # []
  # [min_p]
  #   type = NekVolumeExtremeValue
  #   field = pressure
  #   value_type = min
  # []
  # [max_1]
  #   type = NekVolumeExtremeValue
  #   field = unity
  #   value_type = max
  # []
  # [min_1]
  #   type = NekVolumeExtremeValue
  #   field = unity
  #   value_type = min
  # []

  # # extreme value postprocessors - SIDE
  # [max_T_out]
  #   type = NekSideExtremeValue
  #   field = temperature
  #   boundary = '4'
  #   value_type = max
  # []
  # [min_T_out]
  #   type = NekSideExtremeValue
  #   field = temperature
  #   boundary = '4'
  #   value_type = min
  # []
  # [max_p_in]
  #   type = NekSideExtremeValue
  #   field = pressure
  #   boundary = '3'
  #   value_type = max
  # []
  # [min_p_in]
  #   type = NekSideExtremeValue
  #   field = pressure
  #   boundary = '3'
  #   value_type = min
  # []
  # [max_1_in]
  #   type = NekSideExtremeValue
  #   field = unity
  #   boundary = '3'
  #   value_type = max
  # []
  # [min_1_in]
  #   type = NekSideExtremeValue
  #   field = unity
  #   boundary = '3'
  #   value_type = min
  # []
[]

[Outputs]
  exodus = true
  execute_on = 'final'
  csv = true

  [screen]
    type = Console
    hide = 'average_inlet_T average_outlet_T transfer_in'
  []
[]

Description
Example Input Syntax
Input Parameters
Input Files