NekPointValue

Uses NekRS's pointInterpolation to query the NekRS solution at a point (does not need to be a grid point).

Description

This postprocessor interpolates the NekRS solution onto a given point (using the actual NekRS solution). The provided point does not need to be a Gauss-Lobatto-Legendre (GLL) point.

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

pressure: pressure
temperature: temperature, a.k.a. zeroth passive scalar
velocity: magnitude of velocity
velocity_x: $var element = document.getElementById("moose-equation-b505e312-ffd0-43e2-a075-26a04059a51d");katex.render("x", element, {displayMode:false,throwOnError:false});$ -component of velocity
velocity_y: $var element = document.getElementById("moose-equation-fef71143-d126-4fe0-9faf-02844c0dc263");katex.render("y", element, {displayMode:false,throwOnError:false});$ -component of velocity
velocity_z: $var element = document.getElementById("moose-equation-1988b174-91a1-4011-a550-e9aaa6d03017");katex.render("z", element, {displayMode:false,throwOnError:false});$ -component of velocity
velocity_x_squared: $var element = document.getElementById("moose-equation-32cde448-30cc-40ce-89c9-61400a28fc4f");katex.render("x", element, {displayMode:false,throwOnError:false});$ -component of velocity, squared
velocity_y_squared: $var element = document.getElementById("moose-equation-26f0db05-a977-4f3b-a087-7af592078218");katex.render("y", element, {displayMode:false,throwOnError:false});$ -component of velocity, squared
velocity_z_squared: $var element = document.getElementById("moose-equation-987015c4-2506-41be-bb23-ef7b8959eb72");katex.render("z", element, {displayMode:false,throwOnError:false});$ -component of velocity, squared
velocity_component: velocity vector dotted against another vector
scalar01: first passive scalar
scalar02: second passive scalar
scalar03: third passive scalar
unity: the value 1.0

If running NekRS in non-dimensional form (and you have indicated the appropriate nondimensional scales with the Dimensionalize sub-block for the [Problem]) then the value of this object is shown in dimensional units.

Example Input Syntax

As an example, the following code snippet will interpolate the spectral NekRS solution onto a provided point.

[Postprocessors<<<{"href": "../../syntax/Postprocessors/index.html"}>>>]
  [vx]
    type = NekPointValue<<<{"description": "Uses NekRS's pointInterpolation to query the NekRS solution at a point (does not need to be a grid point).", "href": "NekPointValue.html"}>>>
    field<<<{"description": "Field to apply this object to"}>>> = velocity_x
    point<<<{"description": "The physical point where the field will be evaluated"}>>> = '0.25 0.3 0.27'
  []
  [vy]
    type = NekPointValue<<<{"description": "Uses NekRS's pointInterpolation to query the NekRS solution at a point (does not need to be a grid point).", "href": "NekPointValue.html"}>>>
    field<<<{"description": "Field to apply this object to"}>>> = velocity_y
    point<<<{"description": "The physical point where the field will be evaluated"}>>> = '0.25 0.3 0.27'
  []
  [vz]
    type = NekPointValue<<<{"description": "Uses NekRS's pointInterpolation to query the NekRS solution at a point (does not need to be a grid point).", "href": "NekPointValue.html"}>>>
    field<<<{"description": "Field to apply this object to"}>>> = velocity_z
    point<<<{"description": "The physical point where the field will be evaluated"}>>> = '0.25 0.3 0.27'
  []
  [comp]
    type = NekPointValue<<<{"description": "Uses NekRS's pointInterpolation to query the NekRS solution at a point (does not need to be a grid point).", "href": "NekPointValue.html"}>>>
    field<<<{"description": "Field to apply this object to"}>>> = velocity_component
    velocity_direction<<<{"description": "Unit vector to dot with velocity, for 'field = velocity_component'. For example, velocity_direction = '1 1 0' will get the velocity dotted with (1/sqrt(2), 1/sqrt(2), 0)."}>>> = '0.5 0.5 0.5'
    point<<<{"description": "The physical point where the field will be evaluated"}>>> = '0.25 0.3 0.27'
  []
  [vx2]
    type = NekPointValue<<<{"description": "Uses NekRS's pointInterpolation to query the NekRS solution at a point (does not need to be a grid point).", "href": "NekPointValue.html"}>>>
    field<<<{"description": "Field to apply this object to"}>>> = velocity_x_squared
    point<<<{"description": "The physical point where the field will be evaluated"}>>> = '0.25 0.3 0.27'
  []
  [vy2]
    type = NekPointValue<<<{"description": "Uses NekRS's pointInterpolation to query the NekRS solution at a point (does not need to be a grid point).", "href": "NekPointValue.html"}>>>
    field<<<{"description": "Field to apply this object to"}>>> = velocity_y_squared
    point<<<{"description": "The physical point where the field will be evaluated"}>>> = '0.25 0.3 0.27'
  []
  [vz2]
    type = NekPointValue<<<{"description": "Uses NekRS's pointInterpolation to query the NekRS solution at a point (does not need to be a grid point).", "href": "NekPointValue.html"}>>>
    field<<<{"description": "Field to apply this object to"}>>> = velocity_z_squared
    point<<<{"description": "The physical point where the field will be evaluated"}>>> = '0.25 0.3 0.27'
  []
  [vel]
    type = NekPointValue<<<{"description": "Uses NekRS's pointInterpolation to query the NekRS solution at a point (does not need to be a grid point).", "href": "NekPointValue.html"}>>>
    field<<<{"description": "Field to apply this object to"}>>> = velocity
    point<<<{"description": "The physical point where the field will be evaluated"}>>> = '0.25 0.3 0.27'
  []
  [temp]
    type = NekPointValue<<<{"description": "Uses NekRS's pointInterpolation to query the NekRS solution at a point (does not need to be a grid point).", "href": "NekPointValue.html"}>>>
    field<<<{"description": "Field to apply this object to"}>>> = temperature
    point<<<{"description": "The physical point where the field will be evaluated"}>>> = '0.25 0.3 0.27'
  []
  [p]
    type = NekPointValue<<<{"description": "Uses NekRS's pointInterpolation to query the NekRS solution at a point (does not need to be a grid point).", "href": "NekPointValue.html"}>>>
    field<<<{"description": "Field to apply this object to"}>>> = pressure
    point<<<{"description": "The physical point where the field will be evaluated"}>>> = '0.25 0.3 0.27'
  []
  [scalar01]
    type = NekPointValue<<<{"description": "Uses NekRS's pointInterpolation to query the NekRS solution at a point (does not need to be a grid point).", "href": "NekPointValue.html"}>>>
    field<<<{"description": "Field to apply this object to"}>>> = scalar01
    point<<<{"description": "The physical point where the field will be evaluated"}>>> = '0.25 0.3 0.27'
  []
  [scalar02]
    type = NekPointValue<<<{"description": "Uses NekRS's pointInterpolation to query the NekRS solution at a point (does not need to be a grid point).", "href": "NekPointValue.html"}>>>
    field<<<{"description": "Field to apply this object to"}>>> = scalar02
    point<<<{"description": "The physical point where the field will be evaluated"}>>> = '0.25 0.3 0.27'
  []
  [scalar03]
    type = NekPointValue<<<{"description": "Uses NekRS's pointInterpolation to query the NekRS solution at a point (does not need to be a grid point).", "href": "NekPointValue.html"}>>>
    field<<<{"description": "Field to apply this object to"}>>> = scalar03
    point<<<{"description": "The physical point where the field will be evaluated"}>>> = '0.25 0.3 0.27'
  []
  [unity]
    type = NekPointValue<<<{"description": "Uses NekRS's pointInterpolation to query the NekRS solution at a point (does not need to be a grid point).", "href": "NekPointValue.html"}>>>
    field<<<{"description": "Field to apply this object to"}>>> = unity
    point<<<{"description": "The physical point where the field will be evaluated"}>>> = '0.25 0.3 0.27'
  []
  [usrwrk00]
    type = NekPointValue<<<{"description": "Uses NekRS's pointInterpolation to query the NekRS solution at a point (does not need to be a grid point).", "href": "NekPointValue.html"}>>>
    field<<<{"description": "Field to apply this object to"}>>> = usrwrk00
    point<<<{"description": "The physical point where the field will be evaluated"}>>> = '0.25 0.3 0.27'
  []
  [usrwrk01]
    type = NekPointValue<<<{"description": "Uses NekRS's pointInterpolation to query the NekRS solution at a point (does not need to be a grid point).", "href": "NekPointValue.html"}>>>
    field<<<{"description": "Field to apply this object to"}>>> = usrwrk01
    point<<<{"description": "The physical point where the field will be evaluated"}>>> = '0.25 0.3 0.27'
  []
  [usrwrk02]
    type = NekPointValue<<<{"description": "Uses NekRS's pointInterpolation to query the NekRS solution at a point (does not need to be a grid point).", "href": "NekPointValue.html"}>>>
    field<<<{"description": "Field to apply this object to"}>>> = usrwrk02
    point<<<{"description": "The physical point where the field will be evaluated"}>>> = '0.25 0.3 0.27'
  []
[]

Input Parameters

fieldField to apply this object to
C++ Type:MooseEnum
Options:velocity_x, velocity_y, velocity_z, velocity, velocity_component, velocity_x_squared, velocity_y_squared, velocity_z_squared, temperature, pressure, scalar01, scalar02, scalar03, unity, usrwrk00, usrwrk01, usrwrk02
Controllable:No
Description:Field to apply this object to
pointThe physical point where the field will be evaluated
C++ Type:libMesh::Point
Controllable:No
Description:The physical point where the field will be evaluated

Required Parameters

meshallNekRS mesh to compute postprocessor on
Default:all
C++ Type:MooseEnum
Options:fluid, solid, all
Controllable:No
Description:NekRS mesh to compute postprocessor on
velocity_componentuserDirection in which to evaluate velocity when 'field = velocity_component.' Options: user (you then need to specify a direction with 'velocity_direction'); normal
Default:user
C++ Type:MooseEnum
Options:normal, user
Controllable:No
Description:Direction in which to evaluate velocity when 'field = velocity_component.' Options: user (you then need to specify a direction with 'velocity_direction'); normal
velocity_directionUnit vector to dot with velocity, for 'field = velocity_component'. For example, velocity_direction = '1 1 0' will get the velocity dotted with (1/sqrt(2), 1/sqrt(2), 0).
C++ Type:libMesh::Point
Controllable:No
Description:Unit vector to dot with velocity, for 'field = velocity_component'. For example, velocity_direction = '1 1 0' will get the velocity dotted with (1/sqrt(2), 1/sqrt(2), 0).

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).
execute_onTIMESTEP_ENDThe list of flag(s) indicating when this object should be executed. For a description of each flag, see https://mooseframework.inl.gov/source/interfaces/SetupInterface.html.
Default:TIMESTEP_END
C++ Type:ExecFlagEnum
Options:NONE, INITIAL, LINEAR, LINEAR_CONVERGENCE, NONLINEAR, NONLINEAR_CONVERGENCE, 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. For a description of each flag, see https://mooseframework.inl.gov/source/interfaces/SetupInterface.html.
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

Execution Scheduling Parameters

control_tagsAdds user-defined labels for accessing object parameters via control logic.
C++ Type:std::vector<std::string>
Controllable:No
Description:Adds user-defined labels for accessing object parameters via control logic.
enableTrueSet the enabled status of the MooseObject.
Default:True
C++ Type:bool
Controllable:Yes
Description:Set the enabled status of the MooseObject.
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

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
Unit:(no unit assumed)
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.

Material Property Retrieval Parameters

Input Files

(test/tests/postprocessors/nek_point_value/usrwrk_units.i)
(test/tests/postprocessors/nek_point_value/points_nondimensional.i)
(test/tests/postprocessors/nek_point_value/points.i)
(test/tests/postprocessors/nek_point_value/unknown.i)

(test/tests/postprocessors/nek_point_value/points.i)

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

[Mesh]
  type = NekRSMesh
  volume = true
[]

[Executioner]
  type = Transient

  [TimeStepper]
    type = NekTimeStepper
  []
[]

[Postprocessors]
  [vx]
    type = NekPointValue
    field = velocity_x
    point = '0.25 0.3 0.27'
  []
  [vy]
    type = NekPointValue
    field = velocity_y
    point = '0.25 0.3 0.27'
  []
  [vz]
    type = NekPointValue
    field = velocity_z
    point = '0.25 0.3 0.27'
  []
  [comp]
    type = NekPointValue
    field = velocity_component
    velocity_direction = '0.5 0.5 0.5'
    point = '0.25 0.3 0.27'
  []
  [vx2]
    type = NekPointValue
    field = velocity_x_squared
    point = '0.25 0.3 0.27'
  []
  [vy2]
    type = NekPointValue
    field = velocity_y_squared
    point = '0.25 0.3 0.27'
  []
  [vz2]
    type = NekPointValue
    field = velocity_z_squared
    point = '0.25 0.3 0.27'
  []
  [vel]
    type = NekPointValue
    field = velocity
    point = '0.25 0.3 0.27'
  []
  [temp]
    type = NekPointValue
    field = temperature
    point = '0.25 0.3 0.27'
  []
  [p]
    type = NekPointValue
    field = pressure
    point = '0.25 0.3 0.27'
  []
  [scalar01]
    type = NekPointValue
    field = scalar01
    point = '0.25 0.3 0.27'
  []
  [scalar02]
    type = NekPointValue
    field = scalar02
    point = '0.25 0.3 0.27'
  []
  [scalar03]
    type = NekPointValue
    field = scalar03
    point = '0.25 0.3 0.27'
  []
  [unity]
    type = NekPointValue
    field = unity
    point = '0.25 0.3 0.27'
  []
  [usrwrk00]
    type = NekPointValue
    field = usrwrk00
    point = '0.25 0.3 0.27'
  []
  [usrwrk01]
    type = NekPointValue
    field = usrwrk01
    point = '0.25 0.3 0.27'
  []
  [usrwrk02]
    type = NekPointValue
    field = usrwrk02
    point = '0.25 0.3 0.27'
  []
[]

[Outputs]
  csv = true
[]

(test/tests/postprocessors/nek_point_value/usrwrk_units.i)

[Problem]
  type = NekRSProblem
  casename = 'brick'

  [FieldTransfers]
    [avg_flux]
      type = NekBoundaryFlux
      initial_flux_integral = 1
      usrwrk_slot = 1
      direction = to_nek
    []
    [heat_source]
      type = NekVolumetricSource
      initial_source_integral = 2
      usrwrk_slot = 2
      direction = to_nek
    []
  []

  [Dimensionalize]
    L = 5.0
    U = 0.2
    T = 10.0
    dT = 200.0
    rho = 1000
    Cp = 3000
  []
[]

[Mesh]
  type = NekRSMesh
  boundary = '1'
  volume = true
  scaling = 5.0
[]

[ICs]
  [flux]
    type = FunctionIC
    variable = avg_flux
    function = f
  []
  [source]
    type = FunctionIC
    variable = heat_source
    function = s
  []
[]

[Functions]
  [f]
    type = ParsedFunction
    expression = '(100+200*x+300*y+400*z)/4e4'
  []
  [s]
    type = ParsedFunction
    expression = '(100+200*x+300*y+500*z)/2.937500e+05/1.106383e+00*2'
  []
[]

[Executioner]
  type = Transient

  [TimeStepper]
    type = NekTimeStepper
  []
[]

[Postprocessors]
  [u01]
    type = NekPointValue
    field = usrwrk01
    point = '0.25 0.3 0.27'
  []
  [u02]
    type = NekPointValue
    field = usrwrk02
    point = '0.25 0.3 0.27'
  []
  [integral_flux]
    type = SideIntegralVariablePostprocessor
    variable = avg_flux
    boundary = '1'
  []
  [integral_source]
    type = ElementIntegralVariablePostprocessor
    variable = heat_source
  []
[]

[Outputs]
  csv = true
  hide = 'avg_flux_integral heat_source_integral'
[]

(test/tests/postprocessors/nek_point_value/points_nondimensional.i)

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

  [Dimensionalize]
    L = 5.0
    U = 0.2
    T = 10.0
    dT = 200.0
    rho = 1000
    Cp = 3000

    s01 = 15.0
    ds01 = 250.0
    s02 = 20.0
    ds02 = 300.0
    s03 = 5.0
    ds03 = 100.0
  []
[]

[Mesh]
  type = NekRSMesh
  volume = true
  scaling = 5.0
[]

[Executioner]
  type = Transient

  [TimeStepper]
    type = NekTimeStepper
  []
[]

[Postprocessors]
  [temp]
    type = NekPointValue
    field = temperature
    point = '0.25 0.3 0.27'
  []
  [vx]
    type = NekPointValue
    field = velocity_x
    point = '0.25 0.3 0.27'
  []
  [vy]
    type = NekPointValue
    field = velocity_y
    point = '0.25 0.3 0.27'
  []
  [vz]
    type = NekPointValue
    field = velocity_z
    point = '0.25 0.3 0.27'
  []
  [comp]
    type = NekPointValue
    field = velocity_component
    velocity_direction = '0.5 0.5 0.5'
    point = '0.25 0.3 0.27'
  []
  [vx2]
    type = NekPointValue
    field = velocity_x_squared
    point = '0.25 0.3 0.27'
  []
  [vy2]
    type = NekPointValue
    field = velocity_y_squared
    point = '0.25 0.3 0.27'
  []
  [vz2]
    type = NekPointValue
    field = velocity_z_squared
    point = '0.25 0.3 0.27'
  []
  [vel]
    type = NekPointValue
    field = velocity
    point = '0.25 0.3 0.27'
  []
  [p]
    type = NekPointValue
    field = pressure
    point = '0.25 0.3 0.27'
  []
  [scalar01]
    type = NekPointValue
    field = scalar01
    point = '0.25 0.3 0.27'
  []
  [scalar02]
    type = NekPointValue
    field = scalar02
    point = '0.25 0.3 0.27'
  []
  [scalar03]
    type = NekPointValue
    field = scalar03
    point = '0.25 0.3 0.27'
  []
  [unity]
    type = NekPointValue
    field = unity
    point = '0.25 0.3 0.27'
  []
  [usrwrk00]
    type = NekPointValue
    field = usrwrk00
    point = '0.25 0.3 0.27'
  []
  [usrwrk01]
    type = NekPointValue
    field = usrwrk01
    point = '0.25 0.3 0.27'
  []
  [usrwrk02]
    type = NekPointValue
    field = usrwrk02
    point = '0.25 0.3 0.27'
  []
[]

[Outputs]
  csv = true
[]

(test/tests/postprocessors/nek_point_value/points.i)

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

[Mesh]
  type = NekRSMesh
  volume = true
[]

[Executioner]
  type = Transient

  [TimeStepper]
    type = NekTimeStepper
  []
[]

[Postprocessors]
  [vx]
    type = NekPointValue
    field = velocity_x
    point = '0.25 0.3 0.27'
  []
  [vy]
    type = NekPointValue
    field = velocity_y
    point = '0.25 0.3 0.27'
  []
  [vz]
    type = NekPointValue
    field = velocity_z
    point = '0.25 0.3 0.27'
  []
  [comp]
    type = NekPointValue
    field = velocity_component
    velocity_direction = '0.5 0.5 0.5'
    point = '0.25 0.3 0.27'
  []
  [vx2]
    type = NekPointValue
    field = velocity_x_squared
    point = '0.25 0.3 0.27'
  []
  [vy2]
    type = NekPointValue
    field = velocity_y_squared
    point = '0.25 0.3 0.27'
  []
  [vz2]
    type = NekPointValue
    field = velocity_z_squared
    point = '0.25 0.3 0.27'
  []
  [vel]
    type = NekPointValue
    field = velocity
    point = '0.25 0.3 0.27'
  []
  [temp]
    type = NekPointValue
    field = temperature
    point = '0.25 0.3 0.27'
  []
  [p]
    type = NekPointValue
    field = pressure
    point = '0.25 0.3 0.27'
  []
  [scalar01]
    type = NekPointValue
    field = scalar01
    point = '0.25 0.3 0.27'
  []
  [scalar02]
    type = NekPointValue
    field = scalar02
    point = '0.25 0.3 0.27'
  []
  [scalar03]
    type = NekPointValue
    field = scalar03
    point = '0.25 0.3 0.27'
  []
  [unity]
    type = NekPointValue
    field = unity
    point = '0.25 0.3 0.27'
  []
  [usrwrk00]
    type = NekPointValue
    field = usrwrk00
    point = '0.25 0.3 0.27'
  []
  [usrwrk01]
    type = NekPointValue
    field = usrwrk01
    point = '0.25 0.3 0.27'
  []
  [usrwrk02]
    type = NekPointValue
    field = usrwrk02
    point = '0.25 0.3 0.27'
  []
[]

[Outputs]
  csv = true
[]

(test/tests/postprocessors/nek_point_value/unknown.i)

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

  [Dimensionalize]
    L = 5.0
    U = 0.2
    T = 10.0
    dT = 200.0
    rho = 1000
    Cp = 3000
  []
[]

[Mesh]
  type = NekRSMesh
  volume = true
  scaling = 5.0
[]

[Executioner]
  type = Transient

  [TimeStepper]
    type = NekTimeStepper
  []
[]

[Postprocessors]
  [usrwrk]
    type = NekPointValue
    field = usrwrk02
    point = '0.25 0.3 0.27'
  []
[]

Description
Example Input Syntax
Input Parameters
Input Files