CD_ComputationalGeometry.H

Go to the documentation of this file.

/* chombo-discharge
 * Copyright © 2021 SINTEF Energy Research.
 * Please refer to Copyright.txt and LICENSE in the chombo-discharge root directory.
 */
 
#ifndef CD_ComputationalGeometry_H
#define CD_ComputationalGeometry_H
 
// Chombo includes
#include <BaseIF.H>
#include <MFIndexSpace.H>
#include <Box.H>
#include <RealVect.H>
#include <ProblemDomain.H>
 
// Our includes
#include <CD_Dielectric.H>
#include <CD_Electrode.H>
#include <CD_MultiFluidIndexSpace.H>
#include <CD_RealBox.H>
#include <CD_NamespaceHeader.H>
 
class ComputationalGeometry
{
public:
  ComputationalGeometry();
 
  virtual ~ComputationalGeometry();
 
  const Vector<Dielectric>&
  getDielectrics() const;
 
  const Vector<Electrode>&
  getElectrodes() const;
 
  Real
  getGasPermittivity() const;
 
  void
  useScanShop(const ProblemDomain a_startDomain);
 
  void
  useChomboShop();
 
  void
  setDielectrics(const Vector<Dielectric>& a_dielectrics);
 
  void
  setElectrodes(const Vector<Electrode>& a_electrodes);
 
  void
  setGasPermittivity(const Real a_eps0);
 
  const RefCountedPtr<MultiFluidIndexSpace>&
  getMfIndexSpace() const;
 
  const RefCountedPtr<BaseIF>&
  getGasImplicitFunction() const;
 
  const RefCountedPtr<BaseIF>&
  getSolidImplicitFunction() const;
 
  const RefCountedPtr<BaseIF>&
  getImplicitFunction(const phase::which_phase a_phase) const;
 
  virtual void
  buildGeometries(const ProblemDomain a_finestDomain,
                  const RealVect      a_probLo,
                  const Real          a_finestDx,
                  const int           a_nCellMax,
                  const int           a_maxGhostEB,
                  const int           a_maxCoarsen = -1);
 
protected:
  constexpr static Real s_thresh = 1.E-15;
 
  RefCountedPtr<MultiFluidIndexSpace> m_multifluidIndexSpace;
 
  Real m_eps0;
 
  bool m_useScanShop;
 
  ProblemDomain m_scanDomain;
 
  int m_maxGhostEB;
 
  Vector<Dielectric> m_dielectrics;
 
  Vector<Electrode> m_electrodes;
 
  RefCountedPtr<BaseIF> m_implicitFunctionGas;
 
  RefCountedPtr<BaseIF> m_implicitFunctionSolid;
 
  void
  buildGasGeometry(GeometryService*&   a_geoserver,
                   const ProblemDomain a_finestDomain,
                   const RealVect      a_probLo,
                   const Real          a_finestDx);
 
  void
  buildSolidGeometry(GeometryService*&   a_geoserver,
                     const ProblemDomain a_finestDomain,
                     const RealVect      a_probLo,
                     const Real          a_finestDx);
};
 
#include <CD_NamespaceFooter.H>
 
#endif