- 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
Controllable:No
Description:Field to integrate
NekSideAverage
Compute the average of a field over a boundary of the NekRS mesh
Description
This postprocessor computes the average of a specified field over a boundary in the NekRS mesh,
where is the value of the postprocessor, is the boundary of the NekRS mesh, and is the specified field.
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
(-component of velocity)velocity_y
(-component of velocity)velocity_z
(-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 pressure_in
postprocessor will evaluate the average pressure on the NekRS inlet boundary.
[Postprocessors]
[boundary_flux]
type = NekHeatFluxIntegral
boundary = ${fluid_solid_interface}
[]
[max_nek_T]
type = NekVolumeExtremeValue
field = temperature
value_type = max
[]
[min_nek_T]
type = NekVolumeExtremeValue
field = temperature
value_type = min
[]
[pressure_in]
type = NekSideAverage
field = pressure
boundary = '5'
[]
[mdot_in]
type = NekMassFluxWeightedSideIntegral
field = unity
boundary = '5'
[]
[mdot_out]
type = NekMassFluxWeightedSideIntegral
field = unity
boundary = '6'
[]
[]
(tutorials/fhr_reflector/cht/nek.i)Input 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
- (tutorials/fhr_reflector/cht/nek.i)
- (test/tests/userobjects/nek_scalar_value/nek.i)
- (test/tests/postprocessors/dimensionless_numbers/nondimensional/nek.i)
- (test/tests/nek_standalone/conj_ht/nek.i)
- (test/tests/cht/nondimensional/nek.i)
- (test/tests/nek_standalone/channel/nek.i)
- (tutorials/gas_compact_multiphysics/nek.i)
- (tutorials/gas_compact_cht/nek.i)
- (test/tests/nek_standalone/lowMach/nek.i)
- (test/tests/cht/nondimensional/nek_exact.i)
- (test/tests/postprocessors/nek_side_average/nek.i)
- (tutorials/pebble_67/nek.i)
- (test/tests/cht/sfr_pincell/nek_vpp.i)
- (tutorials/standalone/nek.i)
- (test/tests/nek_standalone/lowMach/nek_boundary.i)
- (test/tests/cht/sfr_pincell/nek.i)
- (test/tests/postprocessors/dimensionless_numbers/dimensional/nek.i)
- (test/tests/conduction/identical_volume/cube/nek.i)
- (test/tests/userobjects/nek_scalar_value/nek_standalone.i)
- (test/tests/multiple_nek_apps/subdirectory/nek.i)
(tutorials/fhr_reflector/cht/nek.i)
fluid_solid_interface = '1 2 7'
[Mesh]
type = NekRSMesh
boundary = ${fluid_solid_interface}
scaling = 0.006
[]
[Problem]
type = NekRSProblem
casename = 'fluid'
nondimensional = true
U_ref = 0.0575
T_ref = 923.15
dT_ref = 10.0
L_ref = 0.006
rho_0 = 1962.13
Cp_0 = 2416.0
[]
[Executioner]
type = Transient
timestep_tolerance = 1e-9
[./TimeStepper]
type = NekTimeStepper
[../]
[]
[Outputs]
exodus = true
[]
[Postprocessors]
[boundary_flux]
type = NekHeatFluxIntegral
boundary = ${fluid_solid_interface}
[]
[max_nek_T]
type = NekVolumeExtremeValue
field = temperature
value_type = max
[]
[min_nek_T]
type = NekVolumeExtremeValue
field = temperature
value_type = min
[]
[pressure_in]
type = NekSideAverage
field = pressure
boundary = '5'
[]
[mdot_in]
type = NekMassFluxWeightedSideIntegral
field = unity
boundary = '5'
[]
[mdot_out]
type = NekMassFluxWeightedSideIntegral
field = unity
boundary = '6'
[]
[]
(tutorials/fhr_reflector/cht/nek.i)
fluid_solid_interface = '1 2 7'
[Mesh]
type = NekRSMesh
boundary = ${fluid_solid_interface}
scaling = 0.006
[]
[Problem]
type = NekRSProblem
casename = 'fluid'
nondimensional = true
U_ref = 0.0575
T_ref = 923.15
dT_ref = 10.0
L_ref = 0.006
rho_0 = 1962.13
Cp_0 = 2416.0
[]
[Executioner]
type = Transient
timestep_tolerance = 1e-9
[./TimeStepper]
type = NekTimeStepper
[../]
[]
[Outputs]
exodus = true
[]
[Postprocessors]
[boundary_flux]
type = NekHeatFluxIntegral
boundary = ${fluid_solid_interface}
[]
[max_nek_T]
type = NekVolumeExtremeValue
field = temperature
value_type = max
[]
[min_nek_T]
type = NekVolumeExtremeValue
field = temperature
value_type = min
[]
[pressure_in]
type = NekSideAverage
field = pressure
boundary = '5'
[]
[mdot_in]
type = NekMassFluxWeightedSideIntegral
field = unity
boundary = '5'
[]
[mdot_out]
type = NekMassFluxWeightedSideIntegral
field = unity
boundary = '6'
[]
[]
(test/tests/userobjects/nek_scalar_value/nek.i)
[Problem]
type = NekRSProblem
casename = 'pyramid'
has_heat_source = false
output = 'temperature'
[]
[Mesh]
type = NekRSMesh
volume = true
[]
[Executioner]
type = Transient
[TimeStepper]
type = NekTimeStepper
[]
[]
[UserObjects]
[scalar1]
type = NekScalarValue
[]
[scalar2]
type = NekScalarValue
[]
[]
[Controls]
[func_control]
type = RealFunctionControl
parameter = 'UserObjects/scalar1/value'
function = 'val'
execute_on = 'timestep_begin'
[]
[func_control2]
type = RealFunctionControl
parameter = 'UserObjects/scalar2/value'
function = 'val2'
execute_on = 'timestep_begin'
[]
[]
[Functions]
[val]
type = ParsedFunction
expression = 'if (t > 1, 200.0, 100.0)'
[]
[val2]
type = ParsedFunction
expression = 'if (t > 1, 400.0, 300.0)'
[]
[]
[Postprocessors]
[avg_t_5]
type = NekSideAverage
field = temperature
boundary = '5'
[]
[avg_t_6]
type = NekSideAverage
field = temperature
boundary = '6'
[]
[]
[Outputs]
exodus = true
csv = true
[]
(test/tests/postprocessors/dimensionless_numbers/nondimensional/nek.i)
[Problem]
type = NekRSStandaloneProblem
casename = 'brick'
nondimensional = true
L_ref = 0.25
U_ref = 0.001
rho_0 = 834.5
Cp_0 = 1228.0
T_ref = 573.0
dT_ref = 10.0
[]
[Mesh]
type = NekRSMesh
boundary = '3'
scaling = 0.25
[]
[Executioner]
type = Transient
[TimeStepper]
type = NekTimeStepper
[]
[]
[Outputs]
[out]
type = CSV
execute_on = 'final'
[]
[]
[Postprocessors]
[Pe]
type = PecletNumber
boundary = '1'
[]
[Re]
type = ReynoldsNumber
boundary = '1'
[]
[area]
type = NekSideIntegral
field = unity
boundary = '1'
[]
[mdot]
type = NekMassFluxWeightedSideIntegral
field = unity
boundary = '1'
[]
[inlet_v]
type = NekSideAverage
field = velocity
boundary = '1'
[]
[]
(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/nek_standalone/channel/nek.i)
[Mesh]
type = NekRSMesh
volume = true
[]
[Problem]
type = NekRSStandaloneProblem
casename = 'channel'
output = 'pressure velocity'
[]
[Executioner]
type = Transient
[TimeStepper]
type = NekTimeStepper
[]
[]
[Postprocessors]
[max_Vx]
type = NekVolumeExtremeValue
field = velocity_x
value_type = max
[]
[max_Vx_output]
type = NodalExtremeValue
variable = vel_x
value_type = max
[]
[max_Vx_diff]
type = DifferencePostprocessor
value1 = max_Vx
value2 = max_Vx_output
[]
[min_Vx]
type = NekVolumeExtremeValue
field = velocity_x
value_type = min
[]
[min_Vx_output]
type = NodalExtremeValue
variable = vel_x
value_type = min
[]
[min_Vx_diff]
type = DifferencePostprocessor
value1 = min_Vx
value2 = min_Vx_output
[]
[max_Vy]
type = NekVolumeExtremeValue
field = velocity_y
value_type = max
[]
[max_Vy_output]
type = NodalExtremeValue
variable = vel_y
value_type = max
[]
[max_Vy_diff]
type = DifferencePostprocessor
value1 = max_Vy
value2 = max_Vy_output
[]
[min_Vy]
type = NekVolumeExtremeValue
field = velocity_y
value_type = min
[]
[min_Vy_output]
type = NodalExtremeValue
variable = vel_y
value_type = min
[]
[max_p]
type = NekVolumeExtremeValue
field = pressure
value_type = max
[]
[max_p_output]
type = NodalExtremeValue
variable = P
value_type = max
[]
[max_p_diff]
type = DifferencePostprocessor
value1 = max_p
value2 = max_p_output
[]
[min_p]
type = NekVolumeExtremeValue
field = pressure
value_type = min
[]
[min_p_output]
type = NodalExtremeValue
variable = P
value_type = min
[]
[min_p_diff]
type = DifferencePostprocessor
value1 = min_p
value2 = min_p_output
[]
[area]
type = NekSideIntegral
field = unity
boundary = '1'
[]
[area_output]
type = AreaPostprocessor
boundary = '1'
[]
[area_diff]
type = DifferencePostprocessor
value1 = area
value2 = area_output
[]
[volume]
type = NekVolumeIntegral
field = unity
[]
[volume_output]
type = VolumePostprocessor
[]
[volume_diff]
type = DifferencePostprocessor
value1 = volume
value2 = volume_output
[]
[max_Vx_side]
type = NekSideExtremeValue
field = velocity_x
value_type = max
boundary = '1'
[]
[max_Vx_side_output]
type = NodalExtremeValue
variable = vel_x
value_type = max
boundary = '1'
[]
[max_Vx_side_diff]
type = DifferencePostprocessor
value1 = max_Vx_side
value2 = max_Vx_side_output
[]
[max_Vy_side]
type = NekSideExtremeValue
field = velocity_y
value_type = max
boundary = '1'
[]
[max_Vy_side_output]
type = NodalExtremeValue
variable = vel_y
value_type = max
boundary = '1'
[]
[max_Vy_side_diff]
type = DifferencePostprocessor
value1 = max_Vy_side
value2 = max_Vy_side_output
[]
[min_Vx_side]
type = NekSideExtremeValue
field = velocity_x
value_type = min
boundary = '1'
[]
[min_Vx_side_output]
type = NodalExtremeValue
variable = vel_x
value_type = min
boundary = '1'
[]
[min_Vx_side_diff]
type = DifferencePostprocessor
value1 = min_Vx_side
value2 = min_Vx_side_output
[]
[min_Vy_side]
type = NekSideExtremeValue
field = velocity_y
value_type = min
boundary = '1'
[]
[min_Vy_side_output]
type = NodalExtremeValue
variable = vel_y
value_type = min
boundary = '1'
[]
[min_Vy_side_diff]
type = DifferencePostprocessor
value1 = min_Vy_side
value2 = min_Vy_side_output
[]
[max_p_side]
type = NekSideExtremeValue
field = pressure
value_type = max
boundary = '1'
[]
[max_p_side_output]
type = NodalExtremeValue
variable = P
value_type = max
boundary = '1'
[]
[min_p_side_diff]
type = DifferencePostprocessor
value1 = min_p_side
value2 = min_p_side_output
[]
[min_p_side]
type = NekSideExtremeValue
field = pressure
value_type = min
boundary = '1'
[]
[min_p_side_output]
type = NodalExtremeValue
variable = P
value_type = min
boundary = '1'
[]
[max_p_side_diff]
type = DifferencePostprocessor
value1 = max_p_side
value2 = max_p_side_output
[]
[avg_p]
type = NekVolumeAverage
field = pressure
[]
[avg_p_output]
type = ElementAverageValue
variable = P
[]
[avg_p_diff]
type = DifferencePostprocessor
value1 = avg_p
value2 = avg_p_output
[]
[avg_Vx]
type = NekVolumeAverage
field = velocity_x
[]
[avg_Vx_output]
type = ElementAverageValue
variable = vel_x
[]
[avg_Vx_diff]
type = DifferencePostprocessor
value1 = avg_Vx
value2 = avg_Vx_output
[]
[avg_Vy]
type = NekVolumeAverage
field = velocity_y
[]
[avg_Vy_output]
type = ElementAverageValue
variable = vel_y
[]
[avg_Vy_diff]
type = DifferencePostprocessor
value1 = avg_Vy
value2 = avg_Vy_output
[]
[avg_Vx_side]
type = NekSideAverage
field = velocity_x
boundary = '1'
[]
[avg_Vx_side_output]
type = SideAverageValue
variable = vel_x
boundary = '1'
[]
[avg_Vx_side_diff]
type = DifferencePostprocessor
value1 = avg_Vx_side
value2 = avg_Vx_side_output
[]
[avg_Vy_side]
type = NekSideAverage
field = velocity_y
boundary = '1'
[]
[avg_Vy_side_output]
type = SideAverageValue
variable = vel_y
boundary = '1'
[]
[avg_Vy_side_diff]
type = DifferencePostprocessor
value1 = avg_Vy_side
value2 = avg_Vy_side_output
[]
[avg_p_side]
type = NekSideAverage
field = pressure
boundary = '1'
[]
[avg_p_side_output]
type = SideAverageValue
variable = P
boundary = '1'
[]
[avg_p_side_diff]
type = DifferencePostprocessor
value1 = avg_p_side
value2 = avg_p_side_output
[]
[]
[Outputs]
csv = true
exodus = true
execute_on = 'final'
hide = 'max_Vx_output min_Vx_output max_Vy_output min_Vy_output max_p_output min_p_output area_output volume_output max_Vx_side max_Vx_side_output max_Vy_side max_Vy_side_output max_p_side max_p_side_output min_Vx_side min_Vx_side_output min_Vy_side min_Vy_side_output min_p_side min_p_side_output avg_p avg_p_output avg_Vx avg_Vx_output avg_Vy avg_Vy_output avg_Vx_side avg_Vx_side_output avg_Vy_side avg_Vy_side_output avg_p_side avg_p_side_output max_Vx max_Vy max_p min_Vx min_Vy min_p area'
[]
(tutorials/gas_compact_multiphysics/nek.i)
# copy-pasta from common_input.i
channel_diameter = 0.016 # diameter of the coolant channels (m)
height = 6.343 # height of the full core (m)
inlet_T = 598.0 # inlet fluid temperature (K)
power = 200e6 # full core power (W)
mdot = 117.3 # fluid mass flowrate (kg/s)
fluid_density = 5.5508 # fluid density (kg/m3)
fluid_Cp = 5189.0 # fluid isobaric specific heat (J/kg/K)
n_bundles = 12 # number of bundles in the full core
n_coolant_channels_per_block = 108 # number of coolant channels per assembly
unit_cell_height = 1.6 # unit cell height - arbitrarily selected
[Mesh]
type = NekRSMesh
boundary = '3'
volume = true
scaling = ${channel_diameter}
[]
[Problem]
type = NekRSProblem
casename = 'ranstube'
has_heat_source = false
n_usrwrk_slots = 2
nondimensional = true
U_ref = ${fparse mdot / (n_bundles * n_coolant_channels_per_block) / fluid_density / (pi * channel_diameter * channel_diameter / 4.0)}
T_ref = ${inlet_T}
dT_ref = ${fparse power / mdot / fluid_Cp * unit_cell_height / height}
L_ref = ${channel_diameter}
rho_0 = ${fluid_density}
Cp_0 = ${fluid_Cp}
synchronization_interval = parent_app
[]
[Executioner]
type = Transient
[TimeStepper]
type = NekTimeStepper
[]
[]
[Outputs]
exodus = true
[screen]
type = Console
hide = 'flux_integral transfer_in'
[]
[csv]
file_base = 'csv/nek'
type = CSV
[]
[]
[Postprocessors]
[inlet_T]
type = NekSideAverage
field = temperature
boundary = '1'
[]
[outlet_T]
type = NekSideAverage
field = temperature
boundary = '2'
[]
[max_T]
type = NekVolumeExtremeValue
field = temperature
value_type = max
[]
[]
(tutorials/gas_compact_cht/nek.i)
# copy-pasta from common_input.i
channel_diameter = 0.016 # diameter of the coolant channels (m)
height = 6.343 # height of the full core (m)
inlet_T = 598.0 # inlet fluid temperature (K)
power = 200e6 # full core power (W)
mdot = 117.3 # fluid mass flowrate (kg/s)
fluid_density = 5.5508 # fluid density (kg/m3)
fluid_Cp = 5189.0 # fluid isobaric specific heat (J/kg/K)
n_bundles = 12 # number of bundles in the full core
n_coolant_channels_per_block = 108 # number of coolant channels per assembly
unit_cell_height = 1.6 # unit cell height - arbitrarily selected
num_layers_for_plots = 50 # number of layers to average fields over for plotting
[Mesh]
type = NekRSMesh
boundary = '3'
volume = true
scaling = ${channel_diameter}
[]
[Problem]
type = NekRSProblem
casename = 'ranstube'
has_heat_source = false
nondimensional = true
U_ref = ${fparse mdot / (n_bundles * n_coolant_channels_per_block) / fluid_density / (pi * channel_diameter * channel_diameter / 4.0)}
T_ref = ${inlet_T}
dT_ref = ${fparse power / mdot / fluid_Cp * unit_cell_height / height}
L_ref = ${channel_diameter}
rho_0 = ${fluid_density}
Cp_0 = ${fluid_Cp}
synchronization_interval = parent_app
[]
[Executioner]
type = Transient
[TimeStepper]
type = NekTimeStepper
min_dt = 1e-10
[]
[]
[Outputs]
exodus = true
[screen]
type = Console
hide = 'boundary_flux inlet_T outlet_T max_T flux_integral transfer_in'
[]
[csv]
file_base = 'csv/nek'
type = CSV
[]
[]
[Postprocessors]
[boundary_flux]
type = NekHeatFluxIntegral
boundary = '3'
[]
[inlet_T]
type = NekSideAverage
field = temperature
boundary = '1'
[]
[outlet_T]
type = NekSideAverage
field = temperature
boundary = '2'
[]
[max_T]
type = NekVolumeExtremeValue
field = temperature
value_type = max
[]
[]
[UserObjects]
active = '' # simply comment out this line to evaluate these user objects
[layered_bin]
type = LayeredBin
num_layers = ${num_layers_for_plots}
direction = z
[]
[wall_temp]
type = NekBinnedSideAverage
bins = 'layered_bin'
boundary = '3'
field = temperature
map_space_by_qp = true
[]
[bulk_temp]
type = NekBinnedVolumeAverage
bins = 'layered_bin'
field = temperature
map_space_by_qp = true
[]
[]
[VectorPostprocessors]
active = '' # simply comment out this line to evaluate these user objects
[wall]
type = SpatialUserObjectVectorPostprocessor
userobject = wall_temp
[]
[bulk]
type = SpatialUserObjectVectorPostprocessor
userobject = bulk_temp
[]
[]
(test/tests/nek_standalone/lowMach/nek.i)
[Mesh]
type = NekRSMesh
volume = true
order = SECOND
[]
[Problem]
type = NekRSStandaloneProblem
casename = 'lowMach'
output = 'pressure velocity temperature'
# We omit the non-dimensional settings here in order to just extract the
# non-dimensional solution as-is, without dimensionalizing it.
[]
[Executioner]
type = Transient
[TimeStepper]
type = NekTimeStepper
[]
[]
[Postprocessors]
# All the following postprocessors are applying operations both (a) directly to the NekRS
# solution arrays, and (b) to the variables extracted with the 'outputs = ...' syntax.
# Rather than check the actual values of these postprocessors (which might change if the
# NekRS development team changes the nature of their CI tests), we can just check that
# the difference between the Nek-style postprocessors from the MOOSE-style postprocessors
# (acting on the extract solution) are nearly zero. We only check the absolute value of
# the min/max volume values for Vx, temperature, and pressure because those values are printed to
# the screen and offer quick confirmation of any changes that are due to changes in NekRS itself.
[max_Vx]
type = NekVolumeExtremeValue
field = velocity_x
value_type = max
[]
[max_Vx_output]
type = NodalExtremeValue
variable = vel_x
value_type = max
[]
[max_Vx_diff]
type = DifferencePostprocessor
value1 = max_Vx
value2 = max_Vx_output
[]
[min_Vx]
type = NekVolumeExtremeValue
field = velocity_x
value_type = min
[]
[min_Vx_output]
type = NodalExtremeValue
variable = vel_x
value_type = min
[]
[min_Vx_diff]
type = DifferencePostprocessor
value1 = min_Vx
value2 = min_Vx_output
[]
[max_p]
type = NekVolumeExtremeValue
field = pressure
value_type = max
[]
[max_p_output]
type = NodalExtremeValue
variable = P
value_type = max
[]
[max_p_diff]
type = DifferencePostprocessor
value1 = max_p
value2 = max_p_output
[]
[min_p]
type = NekVolumeExtremeValue
field = pressure
value_type = min
[]
[min_p_output]
type = NodalExtremeValue
variable = P
value_type = min
[]
[min_p_diff]
type = DifferencePostprocessor
value1 = min_p
value2 = min_p_output
[]
[max_T]
type = NekVolumeExtremeValue
field = temperature
value_type = max
[]
[max_T_output]
type = NodalExtremeValue
variable = temp
value_type = max
[]
[max_T_diff]
type = DifferencePostprocessor
value1 = max_T
value2 = max_T_output
[]
[min_T]
type = NekVolumeExtremeValue
field = temperature
value_type = min
[]
[min_T_output]
type = NodalExtremeValue
variable = temp
value_type = min
[]
[min_T_diff]
type = DifferencePostprocessor
value1 = min_T
value2 = min_T_output
[]
[area]
type = NekSideIntegral
field = unity
boundary = '1'
[]
[area_output]
type = AreaPostprocessor
boundary = '1'
[]
[area_diff]
type = DifferencePostprocessor
value1 = area
value2 = area_output
[]
[volume]
type = NekVolumeIntegral
field = unity
[]
[volume_output]
type = VolumePostprocessor
[]
[volume_diff]
type = DifferencePostprocessor
value1 = volume
value2 = volume_output
[]
[max_T_side]
type = NekSideExtremeValue
field = temperature
value_type = max
boundary = '1'
[]
[max_T_side_output]
type = NodalExtremeValue
variable = temp
value_type = max
boundary = '1'
[]
[max_T_side_diff]
type = DifferencePostprocessor
value1 = max_T_side
value2 = max_T_side_output
[]
[min_T_side]
type = NekSideExtremeValue
field = temperature
value_type = min
boundary = '1'
[]
[min_T_side_output]
type = NodalExtremeValue
variable = temp
value_type = min
boundary = '1'
[]
[min_T_side_diff]
type = DifferencePostprocessor
value1 = min_T_side
value2 = min_T_side_output
[]
[max_Vx_side]
type = NekSideExtremeValue
field = velocity_x
value_type = max
boundary = '1'
[]
[max_Vx_side_output]
type = NodalExtremeValue
variable = vel_x
value_type = max
boundary = '1'
[]
[max_Vx_side_diff]
type = DifferencePostprocessor
value1 = max_Vx_side
value2 = max_Vx_side_output
[]
[min_Vx_side]
type = NekSideExtremeValue
field = velocity_x
value_type = min
boundary = '1'
[]
[min_Vx_side_output]
type = NodalExtremeValue
variable = vel_x
value_type = min
boundary = '1'
[]
[min_Vx_side_diff]
type = DifferencePostprocessor
value1 = min_Vx_side
value2 = min_Vx_side_output
[]
[max_p_side]
type = NekSideExtremeValue
field = pressure
value_type = max
boundary = '1'
[]
[max_p_side_output]
type = NodalExtremeValue
variable = P
value_type = max
boundary = '1'
[]
[min_p_side_diff]
type = DifferencePostprocessor
value1 = min_p_side
value2 = min_p_side_output
[]
[min_p_side]
type = NekSideExtremeValue
field = pressure
value_type = min
boundary = '1'
[]
[min_p_side_output]
type = NodalExtremeValue
variable = P
value_type = min
boundary = '1'
[]
[max_p_side_diff]
type = DifferencePostprocessor
value1 = max_p_side
value2 = max_p_side_output
[]
[avg_T]
type = NekVolumeAverage
field = temperature
[]
[avg_T_output]
type = ElementAverageValue
variable = temp
[]
[avg_T_diff]
type = DifferencePostprocessor
value1 = avg_T
value2 = avg_T_output
[]
[avg_Vx]
type = NekVolumeAverage
field = velocity_x
[]
[avg_Vx_output]
type = ElementAverageValue
variable = vel_x
[]
[avg_Vx_diff]
type = DifferencePostprocessor
value1 = avg_Vx
value2 = avg_Vx_output
[]
[avg_T_side]
type = NekSideAverage
field = temperature
boundary = '1'
[]
[avg_T_side_output]
type = SideAverageValue
variable = temp
boundary = '1'
[]
[avg_T_side_diff]
type = DifferencePostprocessor
value1 = avg_T_side
value2 = avg_T_side_output
[]
[avg_Vx_side]
type = NekSideAverage
field = velocity_x
boundary = '1'
[]
[avg_Vx_side_output]
type = SideAverageValue
variable = vel_x
boundary = '1'
[]
[avg_Vx_side_diff]
type = DifferencePostprocessor
value1 = avg_Vx_side
value2 = avg_Vx_side_output
[]
[avg_p_side]
type = NekSideAverage
field = pressure
boundary = '1'
[]
[avg_p_side_output]
type = SideAverageValue
variable = P
boundary = '1'
[]
[avg_p_side_diff]
type = DifferencePostprocessor
value1 = avg_p_side
value2 = avg_p_side_output
[]
[]
[Outputs]
csv = true
exodus = true
execute_on = 'final'
hide = 'max_Vx_output min_Vx_output max_p_output min_p_output area_output volume_output max_Vx_side max_Vx_side_output max_p_side max_p_side_output min_Vx_side min_Vx_side_output min_p_side min_p_side_output avg_Vx avg_Vx_output avg_Vx_side avg_Vx_side_output avg_p_side avg_p_side_output max_T_output min_T_output max_T_side max_T_side_output min_T_side min_T_side_output avg_T avg_T_output avg_T_side avg_T_side_output'
[]
(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/postprocessors/nek_side_average/nek.i)
[Problem]
type = NekRSProblem
casename = 'pyramid'
[]
[Mesh]
type = NekRSMesh
boundary = '1 2 3 4 5 6 7 8'
[]
[Executioner]
type = Transient
[TimeStepper]
type = NekTimeStepper
[]
[]
[Outputs]
[out]
type = CSV
hide = 'flux_integral'
execute_on = 'final'
[]
[]
[Postprocessors]
[temp_avg1]
type = NekSideAverage
field = temperature
boundary = '1'
[]
[temp_avg2]
type = NekSideAverage
field = temperature
boundary = '2'
[]
[temp_avg3]
type = NekSideAverage
field = temperature
boundary = '3'
[]
[temp_avg4]
type = NekSideAverage
field = temperature
boundary = '4'
[]
[temp_avg5]
type = NekSideAverage
field = temperature
boundary = '5'
[]
[temp_avg6]
type = NekSideAverage
field = temperature
boundary = '6'
[]
[temp_avg7]
type = NekSideAverage
field = temperature
boundary = '7'
[]
[temp_avg8]
type = NekSideAverage
field = temperature
boundary = '8'
[]
[pressure_avg1]
type = NekSideAverage
field = pressure
boundary = '1'
[]
[pressure_avg2]
type = NekSideAverage
field = pressure
boundary = '2'
[]
[pressure_avg3]
type = NekSideAverage
field = pressure
boundary = '3'
[]
[pressure_avg4]
type = NekSideAverage
field = pressure
boundary = '4'
[]
[pressure_avg5]
type = NekSideAverage
field = pressure
boundary = '5'
[]
[pressure_avg6]
type = NekSideAverage
field = pressure
boundary = '6'
[]
[pressure_avg7]
type = NekSideAverage
field = pressure
boundary = '7'
[]
[pressure_avg8]
type = NekSideAverage
field = pressure
boundary = '8'
[]
[velocity_avg1]
type = NekSideAverage
field = velocity
boundary = '1'
[]
[velocity_avg2]
type = NekSideAverage
field = velocity
boundary = '2'
[]
[velocity_avg3]
type = NekSideAverage
field = velocity
boundary = '3'
[]
[velocity_avg4]
type = NekSideAverage
field = velocity
boundary = '4'
[]
[velocity_avg5]
type = NekSideAverage
field = velocity
boundary = '5'
[]
[velocity_avg6]
type = NekSideAverage
field = velocity
boundary = '6'
[]
[velocity_avg7]
type = NekSideAverage
field = velocity
boundary = '7'
[]
[velocity_avg8]
type = NekSideAverage
field = velocity
boundary = '8'
[]
[x_velocity_avg1]
type = NekSideAverage
field = velocity_x
boundary = '1'
[]
[x_velocity_avg2]
type = NekSideAverage
field = velocity_x
boundary = '2'
[]
[x_velocity_avg3]
type = NekSideAverage
field = velocity_x
boundary = '3'
[]
[x_velocity_avg4]
type = NekSideAverage
field = velocity_x
boundary = '4'
[]
[x_velocity_avg5]
type = NekSideAverage
field = velocity_x
boundary = '5'
[]
[x_velocity_avg6]
type = NekSideAverage
field = velocity_x
boundary = '6'
[]
[x_velocity_avg7]
type = NekSideAverage
field = velocity_x
boundary = '7'
[]
[x_velocity_avg8]
type = NekSideAverage
field = velocity_x
boundary = '8'
[]
[y_velocity_avg1]
type = NekSideAverage
field = velocity_y
boundary = '1'
[]
[y_velocity_avg2]
type = NekSideAverage
field = velocity_y
boundary = '2'
[]
[y_velocity_avg3]
type = NekSideAverage
field = velocity_y
boundary = '3'
[]
[y_velocity_avg4]
type = NekSideAverage
field = velocity_y
boundary = '4'
[]
[y_velocity_avg5]
type = NekSideAverage
field = velocity_y
boundary = '5'
[]
[y_velocity_avg6]
type = NekSideAverage
field = velocity_y
boundary = '6'
[]
[y_velocity_avg7]
type = NekSideAverage
field = velocity_y
boundary = '7'
[]
[y_velocity_avg8]
type = NekSideAverage
field = velocity_y
boundary = '8'
[]
[z_velocity_avg1]
type = NekSideAverage
field = velocity_z
boundary = '1'
[]
[z_velocity_avg2]
type = NekSideAverage
field = velocity_z
boundary = '2'
[]
[z_velocity_avg3]
type = NekSideAverage
field = velocity_z
boundary = '3'
[]
[z_velocity_avg4]
type = NekSideAverage
field = velocity_z
boundary = '4'
[]
[z_velocity_avg5]
type = NekSideAverage
field = velocity_z
boundary = '5'
[]
[z_velocity_avg6]
type = NekSideAverage
field = velocity_z
boundary = '6'
[]
[z_velocity_avg7]
type = NekSideAverage
field = velocity_z
boundary = '7'
[]
[z_velocity_avg8]
type = NekSideAverage
field = velocity_z
boundary = '8'
[]
[velocity_comp1]
type = NekSideIntegral
field = velocity_component
velocity_direction = '0.1 -0.2 0.3'
boundary = '1'
[]
[velocity_comp2]
type = NekSideIntegral
field = velocity_component
velocity_direction = '0.1 -0.2 0.3'
boundary = '2'
[]
[velocity_comp3]
type = NekSideIntegral
field = velocity_component
velocity_direction = '0.1 -0.2 0.3'
boundary = '3'
[]
[velocity_comp4]
type = NekSideIntegral
field = velocity_component
velocity_direction = '0.1 -0.2 0.3'
boundary = '4'
[]
[velocity_comp5]
type = NekSideIntegral
field = velocity_component
velocity_direction = '0.1 -0.2 0.3'
boundary = '5'
[]
[velocity_comp6]
type = NekSideIntegral
field = velocity_component
velocity_direction = '0.1 -0.2 0.3'
boundary = '6'
[]
[velocity_comp7]
type = NekSideIntegral
field = velocity_component
velocity_direction = '0.1 -0.2 0.3'
boundary = '7'
[]
[velocity_comp8]
type = NekSideIntegral
field = velocity_component
velocity_direction = '0.1 -0.2 0.3'
boundary = '8'
[]
[]
(tutorials/pebble_67/nek.i)
Re = 1460.0
dp = 0.06
cylinder_diameter = ${fparse 4.4 * dp}
rho = 3.645
Cp = 3121.0
mu = 2.93e-5
power = 2400.0
inlet_area = ${fparse pi * cylinder_diameter^2 / 4.0}
[Mesh]
type = NekRSMesh
boundary = 4
scaling = ${dp}
[]
[Problem]
type = NekRSProblem
casename = 'pb67'
output = 'velocity'
has_heat_source = false
n_usrwrk_slots = 2
nondimensional = true
L_ref = ${dp}
U_ref = ${fparse Re * mu / rho / dp}
T_ref = 523.0
dT_ref = ${fparse power * dp / Re / inlet_area / mu / Cp}
rho_0 = ${rho}
Cp_0 = ${Cp}
synchronization_interval = parent_app
[]
[Executioner]
type = Transient
[TimeStepper]
type = NekTimeStepper
[]
[]
[Outputs]
exodus = true
interval = 100
[]
[Postprocessors]
[max_nek_T]
type = NekVolumeExtremeValue
field = temperature
value_type = max
[]
[min_nek_T]
type = NekVolumeExtremeValue
field = temperature
value_type = min
[]
[average_nek_pebble_T]
type = NekSideAverage
boundary = '4'
field = temperature
[]
[]
(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/standalone/nek.i)
[Mesh]
type = NekRSMesh
volume = true
order = SECOND
[]
[Problem]
type = NekRSStandaloneProblem
casename = 'turbPipe'
output = 'pressure velocity'
[]
[Postprocessors]
[outlet_p]
type = NekSideAverage
boundary = '2'
field = pressure
[]
[inlet_p]
type = NekSideAverage
boundary = '1'
field = pressure
[]
[mdot]
type = NekMassFluxWeightedSideIntegral
boundary = '1'
field = unity
[]
# subtracts the two pressure postprocessors
[dP]
type = DifferencePostprocessor
value1 = outlet_p
value2 = inlet_p
[]
[]
[UserObjects]
[axial_bins]
type = LayeredBin
direction = z
num_layers = 20
[]
[radial_bins]
type = RadialBin
vertical_axis = z
rmax = 0.5
nr = 12
growth_r = 0.9
[]
[volume_averages]
type = NekBinnedVolumeAverage
bins = 'radial_bins axial_bins'
field = velocity_z
map_space_by_qp = true
[]
[]
[AuxVariables]
[volume_averages]
family = MONOMIAL
order = CONSTANT
[]
[]
[AuxKernels]
[volume_averages]
type = SpatialUserObjectAux
variable = volume_averages
user_object = volume_averages
[]
[]
[MultiApps]
[sub]
type = TransientMultiApp
input_files = sub.i
execute_on = timestep_end
[]
[]
[Transfers]
[uo_to_sub]
type = MultiAppGeneralFieldUserObjectTransfer
to_multi_app = sub
source_user_object = volume_averages
variable = avg_velocity
[]
[]
[Executioner]
type = Transient
[TimeStepper]
type = NekTimeStepper
[]
[]
[Outputs]
exodus = true
csv = true
# this hides these values from the screen for neater output
hide = 'outlet_p inlet_p'
[]
(test/tests/nek_standalone/lowMach/nek_boundary.i)
[Mesh]
type = NekRSMesh
boundary = '1'
order = SECOND
[]
[Problem]
type = NekRSStandaloneProblem
casename = 'lowMach'
output = 'pressure velocity temperature'
# We omit the non-dimensional settings here in order to just extract the
# non-dimensional solution as-is, without dimensionalizing it.
[]
[Executioner]
type = Transient
[TimeStepper]
type = NekTimeStepper
[]
[]
[Postprocessors]
[max_T_side]
type = NekSideExtremeValue
field = temperature
value_type = max
boundary = '1'
[]
[max_T_side_output]
type = NodalExtremeValue
variable = temp
value_type = max
[]
[max_T_side_diff]
type = DifferencePostprocessor
value1 = max_T_side
value2 = max_T_side_output
[]
[min_T_side]
type = NekSideExtremeValue
field = temperature
value_type = min
boundary = '1'
[]
[min_T_side_output]
type = NodalExtremeValue
variable = temp
value_type = min
[]
[min_T_side_diff]
type = DifferencePostprocessor
value1 = min_T_side
value2 = min_T_side_output
[]
[max_Vx_side]
type = NekSideExtremeValue
field = velocity_x
value_type = max
boundary = '1'
[]
[max_Vx_side_output]
type = NodalExtremeValue
variable = vel_x
value_type = max
[]
[max_Vx_side_diff]
type = DifferencePostprocessor
value1 = max_Vx_side
value2 = max_Vx_side_output
[]
[min_Vx_side]
type = NekSideExtremeValue
field = velocity_x
value_type = min
boundary = '1'
[]
[min_Vx_side_output]
type = NodalExtremeValue
variable = vel_x
value_type = min
[]
[min_Vx_side_diff]
type = DifferencePostprocessor
value1 = min_Vx_side
value2 = min_Vx_side_output
[]
[max_p_side]
type = NekSideExtremeValue
field = pressure
value_type = max
boundary = '1'
[]
[max_p_side_output]
type = NodalExtremeValue
variable = P
value_type = max
[]
[min_p_side_diff]
type = DifferencePostprocessor
value1 = min_p_side
value2 = min_p_side_output
[]
[min_p_side]
type = NekSideExtremeValue
field = pressure
value_type = min
boundary = '1'
[]
[min_p_side_output]
type = NodalExtremeValue
variable = P
value_type = min
[]
[max_p_side_diff]
type = DifferencePostprocessor
value1 = max_p_side
value2 = max_p_side_output
[]
[avg_T_side]
type = NekSideAverage
field = temperature
boundary = '1'
[]
[avg_T_side_output]
type = ElementAverageValue
variable = temp
[]
[avg_T_side_diff]
type = DifferencePostprocessor
value1 = avg_T_side
value2 = avg_T_side_output
[]
[avg_Vx_side]
type = NekSideAverage
field = velocity_x
boundary = '1'
[]
[avg_Vx_side_output]
type = ElementAverageValue
variable = vel_x
[]
[avg_Vx_side_diff]
type = DifferencePostprocessor
value1 = avg_Vx_side
value2 = avg_Vx_side_output
[]
[avg_p_side]
type = NekSideAverage
field = pressure
boundary = '1'
[]
[avg_p_side_output]
type = ElementAverageValue
variable = P
[]
[avg_p_side_diff]
type = DifferencePostprocessor
value1 = avg_p_side
value2 = avg_p_side_output
[]
[]
[Outputs]
csv = true
exodus = true
execute_on = 'final'
hide = 'max_Vx_side max_Vx_side_output max_p_side max_p_side_output min_Vx_side min_Vx_side_output min_p_side min_p_side_output avg_Vx_side avg_Vx_side_output avg_p_side avg_p_side_output max_T_side max_T_side_output min_T_side min_T_side_output avg_T_side avg_T_side_output'
[]
(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'
[]
[]
(test/tests/postprocessors/dimensionless_numbers/dimensional/nek.i)
[Problem]
type = NekRSProblem
casename = 'brick'
[]
[Mesh]
type = NekRSMesh
boundary = '3'
[]
[Executioner]
type = Transient
[TimeStepper]
type = NekTimeStepper
[]
[]
[Outputs]
[out]
type = CSV
hide = 'flux_integral'
execute_on = 'final'
[]
[]
[Postprocessors]
[Re]
type = ReynoldsNumber
L_ref = 0.25
boundary = '1'
[]
[Pe]
type = PecletNumber
L_ref = 0.25
boundary = '1'
[]
[area]
type = NekSideIntegral
field = unity
boundary = '1'
[]
[mdot]
type = NekMassFluxWeightedSideIntegral
field = unity
boundary = '1'
[]
[inlet_v]
type = NekSideAverage
field = velocity
boundary = '1'
[]
[]
(test/tests/conduction/identical_volume/cube/nek.i)
[Problem]
type = NekRSProblem
casename = 'cube'
# we only technically need two scratch space slots for this problem,
# (the first is reserved for heat flux, but ultimately not used in this
# volume-only problem, while the second is reserved for the volumetric heat
# source, which is actually used),
# so we can skip allocating extra
n_usrwrk_slots = 2
[]
[Mesh]
type = NekRSMesh
order = SECOND
volume = true
[]
[Executioner]
type = Transient
[TimeStepper]
type = NekTimeStepper
[]
[]
[Postprocessors]
[flux_out]
type = NekHeatFluxIntegral
boundary = '1 2 3 4 5 6'
[]
[max_T]
type = NekVolumeExtremeValue
field = temperature
value_type = max
[]
[avg_T_volume]
type = NekVolumeIntegral
field = temperature
[]
[avg_T_back]
type = NekSideAverage
field = temperature
boundary = '6'
[]
[]
[Outputs]
exodus = true
execute_on = 'final'
# for this tests purposes, we only want to check temperature. This keeps the gold file smaller
hide = 'heat_source'
[]
(test/tests/userobjects/nek_scalar_value/nek_standalone.i)
[Problem]
type = NekRSStandaloneProblem
casename = 'pyramid'
n_usrwrk_slots = 1
[]
[Mesh]
type = NekRSMesh
volume = true
[]
[Executioner]
type = Transient
[TimeStepper]
type = NekTimeStepper
[]
[]
[UserObjects]
[scalar1]
type = NekScalarValue
[]
[scalar2]
type = NekScalarValue
[]
[]
[Controls]
[func_control]
type = RealFunctionControl
parameter = 'UserObjects/scalar1/value'
function = 'val'
execute_on = 'timestep_begin'
[]
[func_control2]
type = RealFunctionControl
parameter = 'UserObjects/scalar2/value'
function = 'val2'
execute_on = 'timestep_begin'
[]
[]
[Functions]
[val]
type = ParsedFunction
expression = 'if (t > 1, 200.0, 100.0)'
[]
[val2]
type = ParsedFunction
expression = 'if (t > 1, 400.0, 300.0)'
[]
[]
[Postprocessors]
[avg_t_5]
type = NekSideAverage
field = temperature
boundary = '5'
[]
[avg_t_6]
type = NekSideAverage
field = temperature
boundary = '6'
[]
[]
[Outputs]
exodus = true
csv = true
[]
(test/tests/multiple_nek_apps/subdirectory/nek.i)
[Problem]
type = NekRSProblem
casename = 'pyramid/pyramid'
[]
[Mesh]
type = NekRSMesh
boundary = '1 2 3 4 5 6 7 8'
[]
[Executioner]
type = Transient
[TimeStepper]
type = NekTimeStepper
[]
[]
[Outputs]
[out]
type = CSV
hide = 'flux_integral'
execute_on = 'final'
[]
[]
[Postprocessors]
[temp_avg1]
type = NekSideAverage
field = temperature
boundary = '1'
[]
[temp_avg2]
type = NekSideAverage
field = temperature
boundary = '2'
[]
[temp_avg3]
type = NekSideAverage
field = temperature
boundary = '3'
[]
[temp_avg4]
type = NekSideAverage
field = temperature
boundary = '4'
[]
[temp_avg5]
type = NekSideAverage
field = temperature
boundary = '5'
[]
[temp_avg6]
type = NekSideAverage
field = temperature
boundary = '6'
[]
[temp_avg7]
type = NekSideAverage
field = temperature
boundary = '7'
[]
[temp_avg8]
type = NekSideAverage
field = temperature
boundary = '8'
[]
[pressure_avg1]
type = NekSideAverage
field = pressure
boundary = '1'
[]
[pressure_avg2]
type = NekSideAverage
field = pressure
boundary = '2'
[]
[pressure_avg3]
type = NekSideAverage
field = pressure
boundary = '3'
[]
[pressure_avg4]
type = NekSideAverage
field = pressure
boundary = '4'
[]
[pressure_avg5]
type = NekSideAverage
field = pressure
boundary = '5'
[]
[pressure_avg6]
type = NekSideAverage
field = pressure
boundary = '6'
[]
[pressure_avg7]
type = NekSideAverage
field = pressure
boundary = '7'
[]
[pressure_avg8]
type = NekSideAverage
field = pressure
boundary = '8'
[]
[velocity_avg1]
type = NekSideAverage
field = velocity
boundary = '1'
[]
[velocity_avg2]
type = NekSideAverage
field = velocity
boundary = '2'
[]
[velocity_avg3]
type = NekSideAverage
field = velocity
boundary = '3'
[]
[velocity_avg4]
type = NekSideAverage
field = velocity
boundary = '4'
[]
[velocity_avg5]
type = NekSideAverage
field = velocity
boundary = '5'
[]
[velocity_avg6]
type = NekSideAverage
field = velocity
boundary = '6'
[]
[velocity_avg7]
type = NekSideAverage
field = velocity
boundary = '7'
[]
[velocity_avg8]
type = NekSideAverage
field = velocity
boundary = '8'
[]
[x_velocity_avg1]
type = NekSideAverage
field = velocity_x
boundary = '1'
[]
[x_velocity_avg2]
type = NekSideAverage
field = velocity_x
boundary = '2'
[]
[x_velocity_avg3]
type = NekSideAverage
field = velocity_x
boundary = '3'
[]
[x_velocity_avg4]
type = NekSideAverage
field = velocity_x
boundary = '4'
[]
[x_velocity_avg5]
type = NekSideAverage
field = velocity_x
boundary = '5'
[]
[x_velocity_avg6]
type = NekSideAverage
field = velocity_x
boundary = '6'
[]
[x_velocity_avg7]
type = NekSideAverage
field = velocity_x
boundary = '7'
[]
[x_velocity_avg8]
type = NekSideAverage
field = velocity_x
boundary = '8'
[]
[y_velocity_avg1]
type = NekSideAverage
field = velocity_y
boundary = '1'
[]
[y_velocity_avg2]
type = NekSideAverage
field = velocity_y
boundary = '2'
[]
[y_velocity_avg3]
type = NekSideAverage
field = velocity_y
boundary = '3'
[]
[y_velocity_avg4]
type = NekSideAverage
field = velocity_y
boundary = '4'
[]
[y_velocity_avg5]
type = NekSideAverage
field = velocity_y
boundary = '5'
[]
[y_velocity_avg6]
type = NekSideAverage
field = velocity_y
boundary = '6'
[]
[y_velocity_avg7]
type = NekSideAverage
field = velocity_y
boundary = '7'
[]
[y_velocity_avg8]
type = NekSideAverage
field = velocity_y
boundary = '8'
[]
[z_velocity_avg1]
type = NekSideAverage
field = velocity_z
boundary = '1'
[]
[z_velocity_avg2]
type = NekSideAverage
field = velocity_z
boundary = '2'
[]
[z_velocity_avg3]
type = NekSideAverage
field = velocity_z
boundary = '3'
[]
[z_velocity_avg4]
type = NekSideAverage
field = velocity_z
boundary = '4'
[]
[z_velocity_avg5]
type = NekSideAverage
field = velocity_z
boundary = '5'
[]
[z_velocity_avg6]
type = NekSideAverage
field = velocity_z
boundary = '6'
[]
[z_velocity_avg7]
type = NekSideAverage
field = velocity_z
boundary = '7'
[]
[z_velocity_avg8]
type = NekSideAverage
field = velocity_z
boundary = '8'
[]
[velocity_comp1]
type = NekSideIntegral
field = velocity_component
velocity_direction = '0.1 -0.2 0.3'
boundary = '1'
[]
[velocity_comp2]
type = NekSideIntegral
field = velocity_component
velocity_direction = '0.1 -0.2 0.3'
boundary = '2'
[]
[velocity_comp3]
type = NekSideIntegral
field = velocity_component
velocity_direction = '0.1 -0.2 0.3'
boundary = '3'
[]
[velocity_comp4]
type = NekSideIntegral
field = velocity_component
velocity_direction = '0.1 -0.2 0.3'
boundary = '4'
[]
[velocity_comp5]
type = NekSideIntegral
field = velocity_component
velocity_direction = '0.1 -0.2 0.3'
boundary = '5'
[]
[velocity_comp6]
type = NekSideIntegral
field = velocity_component
velocity_direction = '0.1 -0.2 0.3'
boundary = '6'
[]
[velocity_comp7]
type = NekSideIntegral
field = velocity_component
velocity_direction = '0.1 -0.2 0.3'
boundary = '7'
[]
[velocity_comp8]
type = NekSideIntegral
field = velocity_component
velocity_direction = '0.1 -0.2 0.3'
boundary = '8'
[]
[]