NodeValueAtXFEMInterface.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* 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

#include "GeneralUserObject.h"
#include "ElementPairLocator.h"
#include "MooseVariableFE.h"

class XFEM;
class InterfaceMeshCutUserObjectBase;

class NodeValueAtXFEMInterface : public GeneralUserObject
{
public:
  static InputParameters validParams();

  NodeValueAtXFEMInterface(const InputParameters & parameters);

  virtual ~NodeValueAtXFEMInterface() {}

  virtual void initialize() override;
  virtual void execute() override;
  virtual void finalize() override;

  std::map<unsigned int, Real> getValueAtPositiveLevelSet() const
  {
    return _values_positive_level_set_side;
  };

  std::map<unsigned int, Real> getValueAtNegativeLevelSet() const
  {
    return _values_negative_level_set_side;
  };

  std::map<unsigned int, RealVectorValue> getGradientAtPositiveLevelSet() const
  {
    return _grad_values_positive_level_set_side;
  };

  std::map<unsigned int, RealVectorValue> getGradientAtNegativeLevelSet() const
  {
    return _grad_values_negative_level_set_side;
  };

  unsigned int numberNodes() const { return _nodes.size(); };

protected:
  const Elem * getElemContainingPoint(const Node & p, bool positive_level_set);

  MooseMesh & _mesh;

  std::vector<Point> _nodes;

  std::unique_ptr<libMesh::PointLocatorBase> _pl;

  std::shared_ptr<XFEM> _xfem;

  const ElementPairLocator::ElementPairList * _elem_pairs;

  const InterfaceMeshCutUserObjectBase * _mesh_cut;

  MooseVariableFEBase * _var;

  const unsigned int _level_set_var_number;

  const libMesh::System & _system;

  const NumericVector<Number> & _solution;

  std::map<unsigned int, Real> _values_positive_level_set_side;

  std::map<unsigned int, Real> _values_negative_level_set_side;

  std::map<unsigned int, RealVectorValue> _grad_values_positive_level_set_side;

  std::map<unsigned int, RealVectorValue> _grad_values_negative_level_set_side;
};