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;
};