Pb_Euler.cpp

/****************************************************************************
* Copyright (c) 2026, CEA
* All rights reserved.
*
* Redistribution and use in source and binary forms, with or without modification, are permitted provided that the following conditions are met:
* 1. Redistributions of source code must retain the above copyright notice, this list of conditions and the following disclaimer.
* 2. Redistributions in binary form must reproduce the above copyright notice, this list of conditions and the following disclaimer in the documentation and/or other materials provided with the distribution.
* 3. Neither the name of the copyright holder nor the names of its contributors may be used to endorse or promote products derived from this software without specific prior written permission.
*
* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED.
* IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS;
* OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
*
*****************************************************************************/
 
#include <Milieu_composite_Euler.h>
#include <Pb_Euler.h>
 
Implemente_instanciable(Pb_Euler, "Pb_Euler", Pb_Fluide_base);
// XD Pb_Euler Pb_base Pb_Euler INHERITS_BRACE A problem that allows the resolution of N-phases Euler equations
// XD attr milieu_composite_Euler bloc_lecture milieu_composite_Euler OPT The composite medium associated with the
// XD_CONT problem.
// XD attr momentum_euler momentum_euler qdm_euler REQ Momentum conservation equation for a multi-phase Euler problem
// XD_CONT where the unknown is the velocity
// XD attr masse_euler masse_euler density_euler REQ Mass consevation equation for a multi-phase Euler problem where the
// XD_CONT unknown is the density
// XD attr energy_euler energy_euler energie_euler REQ Internal energy conservation equation for a multi-phase Euler
// XD_CONT problem where the unknown is the temperature
// XD attr fraction_euler fraction_euler fraction_euler REQ Void fraction conservation equation for a multi-phase Euler
// XD_CONT problem where the unknown is the alpha
 
Sortie& Pb_Euler::printOn(Sortie& os) const { return Pb_Fluide_base::printOn(os); }
 
Entree& Pb_Euler::readOn(Entree& is)
{
  if (!discretisation().is_coloc())
    Process::exit("Error: Pb_Euler is only available for Coloc discretization !!! Update your data file ...\n");
 
  if (Objet_U::dimension != 2)
    Process::exit("Error: Pb_Euler is currently only available in 2D ...\n");
 
  return Pb_Fluide_base::readOn(is);
}
 
void Pb_Euler::typer_lire_milieu(Entree& is)
{
  le_milieu_.resize(1);
  is >> le_milieu_[0];
 
  if (!sub_type(Milieu_composite_Euler, le_milieu_[0].valeur()))
    {
      Cerr << "Error: Fluid of type " << le_milieu_[0]->le_type() << " is not compatible with " << que_suis_je() << " problem which accepts only Milieu_composite_Euler medium !!!" << finl;
      Cerr << "Check your datafile!" << finl;
      Process::exit();
    }
 
  noms_phases_ = ref_cast(Milieu_composite_Euler, le_milieu_[0].valeur()).noms_phases();
  associer_milieu_base(le_milieu_[0].valeur());
  Probleme_base::discretiser_equations();
  for (int i = 0; i < nombre_d_equations(); i++)
    equation(i).associer_milieu_equation();
 
  equation(0).milieu().discretiser((*this), la_discretisation_.valeur());
}
 
Entree& Pb_Euler::lire_equations(Entree& is, Motcle& mot)
{
  Cerr << "Pb_Euler : Reading of the equations ..." << finl;
  for (int i = 0; i < nombre_d_equations(); i++)
    {
      is >> mot;
      is >> getset_equation_by_name(mot);
    }
 
  is >> mot; // XXX read last word just after equations !
 
  Cerr << "Pb_Euler : Reading of the equations => OK" << finl;
 
  //correction of data read from .data
  Cerr << "Pb_Euler : Initializing mass's equation unknown " << eq_masse_.inconnue().le_nom() << finl;
  eq_masse_.init_alpha_rho(); // init unknown because the data file initializes rho !!
 
  Cerr << "Pb_Euler : Initializing momentum's equation unknown " << eq_qdm_.inconnue().le_nom() << finl;
  eq_qdm_.init_alpha_rho_u(); // init unknown + phase velocity because the data file initializes velocity !!
 
  Cerr << "Pb_Euler : Initializing energy's equation unknown " << eq_energie_.inconnue().le_nom() << finl;
  eq_energie_.init_energie_tot(); // init unknown because the data file initializes nothing (we don't know the value) ...
 
  return is;
}
 
const Equation_base& Pb_Euler::equation(int i) const
{
  if (i == 0)
    return eq_qdm_;
  else if (i == 1)
    return eq_masse_;
  else if (i == 2)
    return eq_energie_;
  else if (i == 3)
    return eq_fraction_;
  else
    {
      Cerr << "Pb_Euler::equation() : Wrong equation number" << i << "!" << finl;
      Process::exit();
    }
  return eq_qdm_;
}
 
Equation_base& Pb_Euler::equation(int i)
{
  if (i == 0)
    return eq_qdm_;
  else if (i == 1)
    return eq_masse_;
  else if (i == 2)
    return eq_energie_;
  else if (i == 3)
    return eq_fraction_;
  else
    {
      Cerr << "Pb_Euler::equation() : Wrong equation number" << i << "!" << finl;
      Process::exit();
    }
  return eq_qdm_; //to return something
}
 
void Pb_Euler::associer_milieu_base(const Milieu_base& mil)
{
  for (int i = 0; i < nombre_d_equations(); i++)
    equation(i).associer_milieu_base(mil);
}
 
void Pb_Euler::preparer_calcul()
{
  Pb_Fluide_base::preparer_calcul();
  mettre_a_jour(schema_temps().temps_courant());
}
 
void Pb_Euler::mettre_a_jour(double temps)
{
  Probleme_base::mettre_a_jour(temps);
  equation_qdm().mettre_a_jour_p_c();
}