Champ_Face_VDF.h

 
/****************************************************************************
* Copyright (c) 2023, 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_Face_VDF_included
#define Champ_Face_VDF_included
 
#include <Champ_Face_VDF_implementation.h>
#include <Champ_Face_base.h>
#include <Domaine_VDF.h>
 
class Domaine_Cl_VDF;
 
/*! @brief class Champ_Face_VDF
 *
 * @brief This class represents a vector field for which only the normal components
 *   to the faces are computed in VDF. There is therefore only one degree of
 *   freedom per face and the nb_comp_ attribute of a Champ_Face object
 *   equals 1. All components of the field can nevertheless be imposed on
 *   the boundary. If n is the total number of faces in the domain and nb_faces_bord
 *   is the number of boundary faces, the value array associated with the field
 *   is built as follows:
 *       - n values to represent the normal components at faces
 *       - nb_faces_bord*dimension to store the imposed values
 *         on the boundary faces
 *  Note: this class is specific to the VDF module.
 *
 * @sa Champ_Inc_base
 */
class Champ_Face_VDF : public Champ_Face_base, public Champ_Face_VDF_implementation
{
  Declare_instanciable(Champ_Face_VDF);
public:
  int fixer_nb_valeurs_nodales(int) override;
 
  inline const Domaine_VDF& domaine_vdf() const override { return ref_cast(Domaine_VDF, le_dom_VF.valeur()); }
  Champ_base& affecter_(const Champ_base&) override;
  virtual const Champ_Proto& affecter(const double x1, const double x2);
  virtual const Champ_Proto& affecter(const double x1, const double x2, const double x3);
  virtual const Champ_Proto& affecter(const DoubleTab&);
  void verifie_valeurs_cl() override;
  int compo_normale_sortante(int) const;
 
  inline const DoubleTab& tau_diag() const { return tau_diag_; }
  inline const DoubleTab& tau_croises() const { return tau_croises_; }
  inline DoubleTab& tau_diag() { return tau_diag_; }
  inline DoubleTab& tau_croises() { return tau_croises_; }
  void dimensionner_tenseur_Grad();
  void calculer_rotationnel_ordre2_centre_element(DoubleTab&) const;
  int imprime(Sortie&, int) const override;
  DoubleTab& trace(const Frontiere_dis_base&, DoubleTab&, double, int distant) const override;
  void mettre_a_jour(double temps) override;
 
  void calculer_dscald_centre_element(DoubleTab&) const;
 
  void calcul_critere_Q(DoubleTab&, const Domaine_Cl_VDF&);
  void calcul_grad_u(const DoubleTab&, DoubleTab&, const Domaine_Cl_VDF&);
  void calcul_y_plus(DoubleTab&, const Domaine_Cl_VDF&);
  void calcul_y_plus_diphasique(DoubleTab& , const Domaine_Cl_VDF& );
 
  DoubleTab& calcul_duidxj(const DoubleTab&, DoubleTab&) const;
  DoubleTab& calcul_duidxj(const DoubleTab&, DoubleTab&, const Domaine_Cl_VDF&) const;
  DoubleVect& calcul_S_barre(const DoubleTab&, DoubleVect&, const Domaine_Cl_VDF&) const;
  DoubleTab& calcul_S_barre_Multiphase(const DoubleTab&,  DoubleTab&, const Domaine_Cl_VDF&) const;
  DoubleVect& calcul_S_barre_sans_contrib_paroi(const DoubleTab&, DoubleVect&, const Domaine_Cl_VDF&) const;
  void calculer_dercov_axi(const Domaine_Cl_VDF&);
 
  // inline methods
  inline DoubleTab& valeur_aux_faces_post(DoubleTab& result) const override
  {
    return Champ_Face_VDF_implementation::valeur_aux_faces_post_impl(domaine_vdf(), result);
  }
 
  inline DoubleVect& valeur_a_elem(const DoubleVect& position, DoubleVect& val, int le_poly) const override
  {
    return Champ_Face_VDF_implementation::valeur_a_elem(position, val, le_poly);
  }
 
  inline double valeur_a_elem_compo(const DoubleVect& position, int le_poly, int ncomp) const override
  {
    return Champ_Face_VDF_implementation::valeur_a_elem_compo(position, le_poly, ncomp);
  }
 
  inline DoubleTab& valeur_aux_elems(const DoubleTab& positions, const IntVect& les_polys, DoubleTab& tab_valeurs) const override
  {
    return Champ_Face_VDF_implementation::valeur_aux_elems(positions, les_polys, tab_valeurs);
  }
 
  inline DoubleTab& valeur_aux_elems_passe(const DoubleTab& positions, const IntVect& les_polys, DoubleTab& tab_valeurs) const override
  {
    return Champ_Face_VDF_implementation::valeur_aux_elems_passe(positions, les_polys, tab_valeurs);
  }
 
  inline DoubleVect& valeur_aux_elems_compo(const DoubleTab& positions, const IntVect& les_polys, DoubleVect& tab_valeurs, int ncomp) const override
  {
    return Champ_Face_VDF_implementation::valeur_aux_elems_compo(positions, les_polys, tab_valeurs, ncomp);
  }
 
  inline DoubleTab& valeur_aux_sommets(const Domaine& dom, DoubleTab& val) const override
  {
    return Champ_Face_VDF_implementation::valeur_aux_sommets(dom, val);
  }
 
  inline DoubleVect& valeur_aux_sommets_compo(const Domaine& dom, DoubleVect& val, int comp) const override
  {
    return Champ_Face_VDF_implementation::valeur_aux_sommets_compo(dom, val, comp);
  }
 
  inline DoubleTab& remplir_coord_noeuds(DoubleTab& positions) const override
  {
    return Champ_Face_VDF_implementation::remplir_coord_noeuds(positions);
  }
 
  inline int remplir_coord_noeuds_et_polys(DoubleTab& positions, IntVect& polys) const override
  {
    return Champ_Face_VDF_implementation::remplir_coord_noeuds_et_polys(positions, polys);
  }
 
  /* utility for computing source terms: computes the vector v_e + n_f (v_f - v_e. n_f)
     returns: the vector, its norm and the derivatives of the latter with respect to v_e and v_f
  */
  inline double v_norm(const DoubleTab& val, const DoubleTab& val_f, int e, int f, int k, int l, double *v_ext, double *dnv) const
  {
    const Domaine_VDF& domaine = domaine_vdf();
    int d, D = dimension, N = val_f.line_size();
    const DoubleTab& nf = domaine.face_normales();
    const DoubleVect& fs = domaine.face_surfaces();
 
    double scal = 0, vf = f >= 0 ? val_f(f, k) - (l >= 0 ? val_f(f, l) : 0) : 0, v_temp[3], *v = v_ext ? v_ext : v_temp;
    for (d = 0; d < D; d++)
      v[d] = val(e, N*d+k ) - (l >= 0 ? val(e, N*d+l ) : 0);
 
    if (f >= 0)
      for (d = 0, scal = domaine.dot(v, &nf(f, 0)) / fs(f); d < D; d++) v[d] += (vf - scal) * nf(f, d) / fs(f);
 
    double nv = sqrt(domaine.dot(v, v));
 
    if (dnv)
      for (d = 0; d < D; d++) dnv[d] = nv ? (v[d] - (f >= 0 ? vf * nf(f, d) / fs(f) : 0)) / nv : 0;
 
#ifdef _COMPILE_AVEC_PGCC
    if (dnv) dnv[3] = f >= 0 && nv ? vf / (nv+DMINFLOAT) : 0; // nvc++ sucks
#else
    if (dnv) dnv[3] = f >= 0 && nv ? vf / nv : 0;
#endif
    return nv;
  }
 
private:
  double val_imp_face_bord_private(int face, int comp) const;
  double val_imp_face_bord_private(int face, int comp1, int comp2) const;
  inline const Champ_base& le_champ() const override { return *this; }
  inline Champ_base& le_champ() override { return *this; }
 
  DoubleTab tau_diag_;       // diagonal terms of the Grad tensor
  DoubleTab tau_croises_;    // off-diagonal terms of the Grad tensor
};
 
double Champ_Face_coeff_frottement_face_bord(const int, const int , const Domaine_Cl_VDF& zclo);
double Champ_Face_coeff_frottement_grad_face_bord(const int, const int , const Domaine_Cl_VDF& zclo);
double Champ_Face_get_val_imp_face_bord_sym(const DoubleTab& tab_valeurs, const double temp,int face,int comp, const Domaine_Cl_VDF& zclo);
double Champ_Face_get_val_imp_face_bord( const double temp,int face,int comp, const Domaine_Cl_VDF& zclo) ;
double Champ_Face_get_val_imp_face_bord( const double temp,int face,int comp, int comp2, const Domaine_Cl_VDF& zclo) ;
 
#endif /* Champ_Face_VDF_included */