MultiAppGeometricInterpolationTransfer.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

// MOOSE includes
#include "MultiAppConservativeTransfer.h"
#include "MooseVariableFieldBase.h"

#include "libmesh/mesh_base.h"

class MultiAppCoordTransform;
namespace libMesh
{
template <unsigned int>
class InverseDistanceInterpolation;
}

class MultiAppGeometricInterpolationTransfer : public MultiAppConservativeTransfer
{
public:
  static InputParameters validParams();

  MultiAppGeometricInterpolationTransfer(const InputParameters & parameters);

  virtual void execute() override;

protected:
  void fillSourceInterpolationPoints(
      FEProblemBase & from_problem,
      const MooseVariableFieldBase & from_var,
      const MultiAppCoordTransform & from_app_transform,
      std::unique_ptr<libMesh::InverseDistanceInterpolation<Moose::dim>> & idi);

  void interpolateTargetPoints(
      FEProblemBase & to_problem,
      MooseVariableFieldBase & to_var,
      NumericVector<Real> & to_solution,
      const MultiAppCoordTransform & to_app_transform,
      const std::unique_ptr<libMesh::InverseDistanceInterpolation<Moose::dim>> & idi);

  void
  subdomainIDsNode(MooseMesh & mesh, const Node & node, std::set<subdomain_id_type> & subdomainids);

  void computeTransformation(const MooseMesh & mesh,
                             std::unordered_map<dof_id_type, Point> & transformation);

  unsigned int _num_points;
  Real _power;
  MooseEnum _interp_type;
  Real _radius;
  // How much we want to shrink gap
  Real _shrink_gap_width;
  // Which mesh we want to shrink
  MooseEnum _shrink_mesh;
  // Which gap blocks want to exclude during solution transfers
  std::vector<SubdomainName> _exclude_gap_blocks;
  // How small we can consider two points are identical
  Real _distance_tol;

private:
  bool usesMooseAppCoordTransform() const override { return true; }
};