CD_McPhotoImplem.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_MCPHOTOIMPLEM_H
#define CD_MCPHOTOIMPLEM_H
 
// Our includes
#include <CD_McPhoto.H>
#include <CD_NamespaceHeader.H>
 
template <class P, class Ret, Ret (P::*MemberFunc)() const>
void
McPhoto::depositPhotons(EBAMRCellData&        a_phi,
                        ParticleContainer<P>& a_photons,
                        const DepositionType& a_deposition) const noexcept
{
  CH_TIME("McPhoto::depositPhotons(ParticleContainer)");
  if (m_verbosity > 5) {
    pout() << m_name + "::depositPhotons(ParticleContainer)" << endl;
  }
 
  // a_phi contains only weights, i.e. not divided by kappa
  this->depositKappaConservative<P, Ret, MemberFunc>(a_phi, a_photons, a_deposition, m_coarseFineDeposition);
 
  // Compute m_depositionNC = sum(kappa*Wc)/sum(kappa)
  this->depositNonConservative(m_depositionNC, a_phi);
 
  // Compute hybrid deposition, including mass difference
  this->depositHybrid(a_phi, m_massDiff, m_depositionNC);
 
  // Redistribute
  if (m_blendConservation) {
    Vector<RefCountedPtr<EBFluxRedistribution>>& redistOps = m_amr->getRedistributionOp(m_realm, m_phase);
    for (int lvl = 0; lvl <= m_amr->getFinestLevel(); lvl++) {
      const Real     scale     = 1.0;
      const Interval variables = Interval(0, 0);
      const bool     hasCoar   = lvl > 0;
      const bool     hasFine   = lvl < m_amr->getFinestLevel();
 
      if (hasCoar) {
        redistOps[lvl]->redistributeCoar(*a_phi[lvl - 1], *m_massDiff[lvl], scale, variables);
      }
 
      redistOps[lvl]->redistributeLevel(*a_phi[lvl], *m_massDiff[lvl], scale, variables);
 
      if (hasFine) {
        redistOps[lvl]->redistributeFine(*a_phi[lvl + 1], *m_massDiff[lvl], scale, variables);
      }
    }
  }
 
  // Average down and interpolate
  m_amr->conservativeAverage(a_phi, m_realm, m_phase);
  m_amr->interpGhost(a_phi, m_realm, m_phase);
}
 
template <class P, class Ret, Ret (P::*MemberFunc)() const>
void
McPhoto::depositKappaConservative(EBAMRCellData&             a_phi,
                                  ParticleContainer<P>&      a_particles,
                                  const DepositionType       a_deposition,
                                  const CoarseFineDeposition a_coarseFineDeposition) const noexcept
{
  CH_TIME("McPhoto::depositKappaConservative");
  if (m_verbosity > 5) {
    pout() << m_name + "::depositKappaConservative" << endl;
  }
 
  CH_assert(a_phi[0]->nComp() == 1);
 
  m_amr->depositParticles<P, Ret, MemberFunc>(a_phi,
                                              m_realm,
                                              m_phase,
                                              a_particles,
                                              a_deposition,
                                              a_coarseFineDeposition,
                                              false);
}
 
#include <CD_NamespaceFooter.H>
 
#endif