MFEMComplexIntegratedBC

Summary

Base class for objects applying complex-valued integrated boundary conditions to an MFEM FE problem.

Overview

MFEMComplexIntegratedBC applies boundary integrator(s) to the weak form equation that is labeled according to the test variable name returned from getTestVariableName(), similar to MFEMComplexKernel for domain integrators.

Similarly to how MFEMComplexKernel works, MFEMComplexIntegratedBC is a container for two MFEMIntegratedBC objects, one representing the real part of the integrated BC, and the other representing the imaginary part. These two can be set up using the RealComponent and ImagComponent sub-blocks on the script. The two integrators need not be the same, but they do need to be applied to the same variable.

Input Parameters

boundary-1 The list of boundaries (ids or names) from the mesh where this object applies. Defaults to all boundaries.
Default:-1
C++ Type:std::vector<BoundaryName>
Controllable:No
Description:The list of boundaries (ids or names) from the mesh where this object applies. Defaults to all boundaries.
variableVariable on which to apply the boundary condition
C++ Type:VariableName
Unit:(no unit assumed)
Controllable:No
Description:Variable on which to apply the boundary condition

Optional Parameters

control_tagsAdds user-defined labels for accessing object parameters via control logic.
C++ Type:std::vector<std::string>
Controllable:No
Description:Adds user-defined labels for accessing object parameters via control logic.
enableTrueSet the enabled status of the MooseObject.
Default:True
C++ Type:bool
Controllable:No
Description:Set the enabled status of the MooseObject.

Advanced Parameters

Child Objects

(framework/include/mfem/bcs/MFEMRWTE10IntegratedBC.h)

(framework/include/mfem/bcs/MFEMRWTE10IntegratedBC.h)

// 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

#ifdef MOOSE_MFEM_ENABLED

#pragma once

#include "MFEMComplexIntegratedBC.h"

class MFEMRWTE10IntegratedBC : public MFEMComplexIntegratedBC
{
public:
  static InputParameters validParams();

  MFEMRWTE10IntegratedBC(const InputParameters & parameters);

  void RWTE10(const mfem::Vector & x, std::vector<std::complex<mfem::real_t>> & E);
  void RWTE10Real(const mfem::Vector & x, mfem::Vector & v);
  void RWTE10Imag(const mfem::Vector & x, mfem::Vector & v);

  virtual mfem::LinearFormIntegrator * getRealLFIntegrator() override
  {
    return (getParam<bool>("input_port") ? new mfem::VectorFEBoundaryTangentLFIntegrator(*_u_real)
                                         : nullptr);
  }
  virtual mfem::LinearFormIntegrator * getImagLFIntegrator() override
  {
    return (getParam<bool>("input_port") ? new mfem::VectorFEBoundaryTangentLFIntegrator(*_u_imag)
                                         : nullptr);
  }
  virtual mfem::BilinearFormIntegrator * getRealBFIntegrator() override { return nullptr; }
  virtual mfem::BilinearFormIntegrator * getImagBFIntegrator() override
  {
    return new mfem::VectorFEMassIntegrator(_robin_coef_im.get());
  }

  /// Computes the unit vector in the direction of the cross product of two 3d vectors
  mfem::Vector normalizedCrossProduct(const mfem::Vector & va, const mfem::Vector & vb)
  {
    mfem::Vector vec;
    va.cross3D(vb, vec);
    return vec /= vec.Norml2();
  }

protected:
  mfem::real_t _mu;
  mfem::real_t _epsilon;
  mfem::real_t _omega;

  mfem::Vector _a1;
  mfem::Vector _a2;

  mfem::real_t _k0;
  mfem::real_t _kc;
  std::complex<mfem::real_t> _k;

  mfem::Vector _k_c;
  mfem::Vector _k_a;

  std::unique_ptr<mfem::ConstantCoefficient> _robin_coef_im{nullptr};
  std::unique_ptr<mfem::VectorFunctionCoefficient> _u_real{nullptr};
  std::unique_ptr<mfem::VectorFunctionCoefficient> _u_imag{nullptr};
};

#endif

Summary
Overview
Input Parameters
Child Objects