Champ_front_contact_VEF.h

/****************************************************************************
* Copyright (c) 2024, 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.
*
*****************************************************************************/
 
#ifndef Champ_front_contact_VEF_included
#define Champ_front_contact_VEF_included
 
#include <Ch_front_var_instationnaire_dep.h>
#include <Domaine_VEF.h>
#include <TRUST_Ref.h>
#include <Front_VF.h>
#include <Motcle.h>
 
class Equation_base;
class Milieu_base;
class Domaine_dis_base;
class Domaine_Cl_dis_base;
class Champ_front_contact_VEF;
class Champ_Inc_base;
class Front_dis_base;
class Equation_base;
class Milieu_base;
 
/*! @brief class Champ_front_contact_VEF Enables scalar coupling (temperature or concentration) between problems by computing
 *
 *      the wall value that equalizes the fluxes.
 *      calculer_coeffs_echange stores grad_num_local and grad_fro_local such that (up to sign?)
 *      flux = grad_num_local + Scal_paroi*grad_fro_local
 *      mettre_a_jour retrieves the grad* from the opposite side and computes Scal_paroi
 *      so as to equalize the fluxes.
 *
 * @sa Champ_front_var_instationnaire Champ_Inc
 */
class Champ_front_contact_VEF : public Ch_front_var_instationnaire_dep
{
 
  Declare_instanciable(Champ_front_contact_VEF);
 
public:
  int initialiser(double temps, const Champ_Inc_base& inco) override;
  Champ_front_base& affecter_(const Champ_front_base& ch) override ;
  void mettre_a_jour(double temps) override;
  void calculer_coeffs_echange(double temps) override;
  void creer(const Nom&, const Nom&,const Nom&, const Nom&, const Motcle&);
  const Champ_Inc_base& inconnue() const;
  const Equation_base& equation() const;
  const Domaine_dis_base& domaine_dis() const override;
  const Domaine_Cl_dis_base& domaine_Cl_dis() const ;
  const Frontiere_dis_base& front_dis() const ;
  void associer_ch_inc_base(const Champ_Inc_base&);
  const Nom& nom_bord_oppose() const;
  const Milieu_base& milieu() const;
 
  virtual void calcul_grads_transf(double temps);
  virtual void calcul_grads_locaux(double temps);
  void calcul_coeff_amort();
  void recupere_gradients_autre_pb(const Front_VF&, DoubleVect&, DoubleVect&);
  DoubleVect& trace_face_raccord(const Front_VF& fr_vf,const DoubleVect& y,DoubleVect& x);
  void remplir_connect_bords();
  void verifier_scalaire_bord(double temps);
  void associer_front_vf_et_ch_front_autre_pb();
 
  void connectivity_failed(const Domaine_VEF&,  int&,  int&,
                           const Domaine_VEF&,  int&,  int&, Nom&, Nom&) ;
  void remplir_faces_coin();
  void remplir_elems_voisin_bord();
  void test_faces_coin();
 
  inline DoubleVect& set_gradient_num_local();
  inline const DoubleVect& get_gradient_num_local() const;
  inline DoubleVect& set_gradient_fro_local();
  inline  const DoubleVect& get_gradient_fro_local() const;
 
  inline  DoubleVect& set_gradient_num_transf();
  inline  const DoubleVect& get_gradient_num_transf() const;
  inline  DoubleVect& set_gradient_fro_transf();
  inline  const DoubleVect& get_gradient_fro_transf() const;
 
  inline  DoubleVect& set_coeff_amort_num();
  inline  const DoubleVect& get_coeff_amort_num() const;
  inline DoubleVect& set_coeff_amort_denum();
  inline  const DoubleVect& get_coeff_amort_denum() const;
 
  inline const IntVect& get_elems_voisin_bord() const;
  inline const Nom& get_nom_bord1() const;
  inline const IntVect& get_connect_bords() const;
  inline const IntVect& get_faces_coin() const;
 
protected :
  int is_conduction=0; // Are we in a conduction problem?
  OBS_PTR(Champ_Inc_base) l_inconnue1;
  OBS_PTR(Champ_Inc_base) l_inconnue2;
  OBS_PTR(Champ_Inc_base) l_inconnue;
  OBS_PTR(Front_VF) fr_vf_autre_pb;
  OBS_PTR(Champ_front_contact_VEF) ch_fr_autre_pb;
 
  Nom nom_bord1;
  Nom nom_bord2;
  Nom nom_bord ;
  Nom nom_pb1;
  Nom nom_pb2;
 
  IntVect connect_bords; // Connectivity array of the boundary faces of the two conjugate boundaries (sequential)
  int connect_est_remplit = -1;
 
  DoubleVect gradient_num_local;
  DoubleVect gradient_fro_local;
  DoubleVect gradient_num_transf;
  DoubleVect gradient_fro_transf;
  DoubleVect Scal_moy;
 
  DoubleVect coeff_amort_num;
  DoubleVect coeff_amort_denum;
 
  IntVect elems_voisin_bord_; //Contains the index of the neighboring element for each boundary face
  IntVect faces_coin;         //Takes the value 1 if the boundary face is a corner face (0 otherwise)
  int verification_faces_coin=0;
};
 
 
 
 
 
inline DoubleVect& Champ_front_contact_VEF::set_gradient_num_local()
{
  return gradient_num_local;
}
 
inline const DoubleVect& Champ_front_contact_VEF::get_gradient_num_local() const
{
  return gradient_num_local;
}
 
inline DoubleVect& Champ_front_contact_VEF::set_gradient_fro_local()
{
  return gradient_fro_local;
}
 
inline const DoubleVect& Champ_front_contact_VEF::get_gradient_fro_local() const
{
  return gradient_fro_local;
}
 
 
inline  DoubleVect& Champ_front_contact_VEF::set_gradient_num_transf()
{
  return gradient_num_transf;
}
 
inline  const DoubleVect& Champ_front_contact_VEF::get_gradient_num_transf() const
{
  return gradient_num_transf;
}
 
inline  DoubleVect& Champ_front_contact_VEF::set_gradient_fro_transf()
{
  return gradient_fro_transf;
}
 
inline  const DoubleVect& Champ_front_contact_VEF::get_gradient_fro_transf() const
{
  return gradient_fro_transf;
}
 
inline  DoubleVect& Champ_front_contact_VEF::set_coeff_amort_num()
{
  return coeff_amort_num;
}
 
inline  const DoubleVect& Champ_front_contact_VEF::get_coeff_amort_num() const
{
  return coeff_amort_num;
}
 
inline DoubleVect& Champ_front_contact_VEF::set_coeff_amort_denum()
{
  return coeff_amort_denum;
}
 
inline const DoubleVect& Champ_front_contact_VEF::get_coeff_amort_denum() const
{
  return coeff_amort_denum;
}
 
inline const IntVect& Champ_front_contact_VEF::get_elems_voisin_bord() const
{
  return elems_voisin_bord_;
}
 
inline const Nom& Champ_front_contact_VEF::get_nom_bord1() const
{
  return nom_bord1;
}
 
inline const IntVect& Champ_front_contact_VEF::get_connect_bords() const
{
  return connect_bords;
}
 
inline const IntVect& Champ_front_contact_VEF::get_faces_coin() const
{
  return faces_coin;
}
 
#endif