CD_ComputationalGeometry.H

Go to the documentation of this file.

/*
 * SPDX-FileCopyrightText: 2021-2026 SINTEF Energy Research
 *
 * SPDX-License-Identifier: GPL-3.0-or-later
 */
 
#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_beginDomain);
 
  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