GaussianEnergyFluxBC.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 "ADIntegratedBC.h"
 
class GaussianEnergyFluxBC : public ADIntegratedBC
{
public:
  static InputParameters validParams();
  static InputParameters beamParams();
 
  GaussianEnergyFluxBC(const InputParameters & params);
 
  template <typename T, typename PointType>
  static ADReal beamFlux(const T & flux_obj, const PointType & flux_obj_location);
 
protected:
  virtual ADReal computeQpResidual() override;
 
  const Real _P0;
  const Real _R;
  const Function & _x_beam_coord;
  const Function & _y_beam_coord;
  const Function & _z_beam_coord;
};
 
template <typename T, typename PointType>
ADReal
GaussianEnergyFluxBC::beamFlux(const T & flux_obj, const PointType & flux_obj_location)
{
  const Point origin(0, 0, 0);
  const RealVectorValue beam_coords{flux_obj._x_beam_coord.value(flux_obj._t, origin),
                                    flux_obj._y_beam_coord.value(flux_obj._t, origin),
                                    flux_obj._z_beam_coord.value(flux_obj._t, origin)};
  const auto r = (flux_obj_location - beam_coords).norm();
  const auto R2 = flux_obj._R * flux_obj._R;
 
  using std::exp;
  return -2 * flux_obj._P0 / (libMesh::pi * R2) * exp(-2 * r * r / R2);
}