CD_MFHelmholtzElectrostaticEBBCImplem.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_MFHELMHOLTZELECTROSTATICEBBCIMPLEM_H
#define CD_MFHELMHOLTZELECTROSTATICEBBCIMPLEM_H
 
// Std includes
#include <limits>
 
// Our includes
#include <CD_MFHelmholtzElectrostaticEBBC.H>
#include <CD_NamespaceHeader.H>
 
Real
MFHelmholtzElectrostaticEBBC::getElectrodePotential(
  const RealVect& a_pos) const // NOLINT(readability-convert-member-functions-to-static)
{
  CH_TIME("MFHelmholtzElectrostaticEBBC::getElectrodePotential(RealVect)");
 
  // Find closest electrode
  int  closestElectrode = 0;
  Real minDist          = std::numeric_limits<Real>::infinity();
 
  const std::vector<std::pair<Electrode, ElectrostaticEbBc::BcFunction>>& electrodeBCs = m_electrostaticBCs.getBcs();
 
  for (int i = 0; i < electrodeBCs.size(); i++) {
    const RefCountedPtr<BaseIF>& impFunc = electrodeBCs[i].first.getImplicitFunction();
 
    const Real curDist = std::abs(impFunc->value(a_pos));
 
    if (curDist < minDist) {
      closestElectrode = i;
      minDist          = curDist;
    }
  }
 
  // Return potential of closest electrode. Again, following the (perhaps odd) convention that the FieldSolver is "time dependent" but the operator
  // factory is not, we've passed the time in by reference to the functions in m_electrostaticBCs. Calling those functions simply ignore the time-argument
  // and uses FieldSolver::m_dt instead. So, we can just use a dummy dt here.
  constexpr Real dummyDt = 0.0;
 
  return electrodeBCs[closestElectrode].second(a_pos, dummyDt);
}
 
#include <CD_NamespaceFooter.H>
 
#endif