FluxAnodeCathodeSolver.cpp

Go to the documentation of this file.

/*~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
   This file is part of CEPS.
 
   CEPS is free software: you can redistribute it and/or modify
   it under the terms of the GNU General Public License as published by
   the Free Software Foundation, either version 3 of the License, or
   (at your option) any later version.
 
   CEPS is distributed in the hope that it will be useful,
   but WITHOUT ANY WARRANTY; without even the implied warranty of
   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
   GNU General Public License for more details.
 
   You should have received a copy of the GNU General Public License
   along with CEPS (see file LICENSE at root of project).
   If not, see <https://www.gnu.org/licenses/>.
 
 
   Copyright 2019-2024 Inria, Universite de Bordeaux
 
   Authors, in alphabetical order:
 
     Pierre-Elliott BECUE, Florian CARO, Yves COUDIERE(*), Andjela DAVIDOVIC, 
     Charlie DOUANLA-LONTSI, Marc FUENTES, Mehdi JUHOOR, Michael LEGUEBE(*), 
     Pauline MIGERDITICHAN, Valentin PANNETIER(*), Nejib ZEMZEMI.
     * : currently active authors
 
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*/
#include "pde/solver/static/FluxAnodeCathodeSolver.hpp"
#include "pde/problem/static/FluxAnodeCathodeProblem.hpp"
#include "geometry/selectors/AttributesSelector.hpp"
 
FluxAnodeCathodeSolver::FluxAnodeCathodeSolver(FluxAnodeCathodeProblem * problem)
  : LaplacianSolver(problem)
{
  FluxAnodeCathodeSolver::setupWithParameters(problem->getParameters());
}
 
FluxAnodeCathodeSolver::~FluxAnodeCathodeSolver()
{
}
 
void 
FluxAnodeCathodeSolver::setupWithParameters(InputParameters* params)
{
  FluxAnodeCathodeProblem*  pb        = getFluxAnodeCathodeProblem();
  BoundaryConditionManager* bcManager = pb->getBoundaryConditionManager();
 
  auto     bcFE        = m_fe->getBoundaryFiniteElements();
  CepsReal anodeMeas   = getDomainMeasure(bcFE,pb->getAnodeAttributes  ());
  CepsReal cathodeMeas = getDomainMeasure(bcFE,pb->getCathodeAttributes());
 
 
  bcManager->getNeumannBCs()->at("ANODE"  )->setScaleFactor(1./anodeMeas  );
  bcManager->getNeumannBCs()->at("CATHODE")->setScaleFactor(1./cathodeMeas);
}
 
void
FluxAnodeCathodeSolver::solve()
{
  LaplacianSolver::solve();
 
  FEMassIntegrator integrator(m_fe);
  CepsReal uMean = integrator.integrate(m_solution, m_fe->getMassMatrix()); // integrate 1.M.u
  CEPS_SAYS("  Solution vector with mean " << uMean);
}
 
FluxAnodeCathodeProblem*
FluxAnodeCathodeSolver::getFluxAnodeCathodeProblem() const
{
  return ceps::runtimeCast<FluxAnodeCathodeProblem*>(m_problem);
}
 
void
FluxAnodeCathodeSolver::initializeAssemblers()
{
  if (getFluxAnodeCathodeProblem()->requireNullMean())
  {
    using LhsAssembler = FENullMeanConstraintAssembler<FEDivKGradAssembler>;
 
    LhsAssembler* opAsb = ceps::getNew<LhsAssembler>(m_fe);
    
    // Set diffusion coefficient for the single unknown and for the lagrangian
    opAsb->setKForUnknown(m_problem->getUnknown(0),getFluxAnodeCathodeProblem()->getDiffusionTensor());
    opAsb->addLagrangian(m_problem->getUnknown(1));
 
    m_opAsb = opAsb;
  }
  else 
  {
    FEDivKGradAssembler* opAsb = ceps::getNew<FEDivKGradAssembler>(m_fe);
 
    // Set diffusion coefficient for the single unknown and for the lagrangian
    opAsb->setKForUnknown(m_problem->getUnknown(0),getFluxAnodeCathodeProblem()->getDiffusionTensor());
    m_opAsb = opAsb;
  }
 
  m_opAsb->setMatrix(m_lhs.get());
 
  m_bcAsb  = ceps::getNew<FEDivKGradBCAssembler>(m_fe);
  m_bcAsb->setMatrix(m_lhs.get());
  m_bcAsb->setVector(m_rhs.get());
}