Demonstrate contains point

framework/include/auxkernels/ContainsPointAux.h

@@ -0,0 +1,26 @@
+//* This file is part of the MOOSE framework
+//* https://www.mooseframework.org
+//*
+//* All rights reserved, see COPYRIGHT for full restrictions
+//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
+//*
+//* Licensed under LGPL 2.1, please see LICENSE for details
+//* https://www.gnu.org/licenses/lgpl-2.1.html
+
+#pragma once
+
+// MOOSE includes
+#include "AuxKernel.h"
+
+class ContainsPointAux : public AuxKernel
+{
+public:
+  static InputParameters validParams();
+
+  ContainsPointAux(const InputParameters & parameters);
+
+protected:
+  virtual Real computeValue() override;
+
+  const Point _point;
+};

framework/src/auxkernels/ContainsPointAux.C

@@ -0,0 +1,35 @@
+//* This file is part of the MOOSE framework
+//* https://www.mooseframework.org
+//*
+//* All rights reserved, see COPYRIGHT for full restrictions
+//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
+//*
+//* Licensed under LGPL 2.1, please see LICENSE for details
+//* https://www.gnu.org/licenses/lgpl-2.1.html
+
+// MOOSE includes
+#include "ContainsPointAux.h"
+#include "MooseMesh.h"
+
+registerMooseObject("MooseApp", ContainsPointAux);
+
+InputParameters
+ContainsPointAux::validParams()
+{
+  InputParameters params = AuxKernel::validParams();
+  params.addRequiredParam<Point>("point", "The point we're checking containment of");
+  return params;
+}
+
+ContainsPointAux::ContainsPointAux(const InputParameters & parameters)
+  : AuxKernel(parameters), _point(getParam<Point>("point"))
+{
+  if (isNodal())
+    paramError("variable", "This AuxKernel only supports Elemental fields");
+}
+
+Real
+ContainsPointAux::computeValue()
+{
+  return _current_elem->contains_point(_point);
+}

modules/tensor_mechanics/child.i

@@ -41,6 +41,16 @@ gamma = 0.5
     []
     initial_condition = 1.0
   []
+  [contains_point]
+    order = CONSTANT
+    family = MONOMIAL
+    [AuxKernel]
+      type = ContainsPointAux
+      point = '.104 .356 0'
+      execute_on = 'initial timestep_end'
+      use_displaced_mesh = true
+    []
+  []
 []
 
 [Modules/TensorMechanics/DynamicMaster]