Testvs19_NujicSourceTreatment.h

Go to the documentation of this file.

 
#pragma once
 
// BaseNumTypes include
#include "Sfx/Sfx_scalarField1D.h"
// DamBreak App inludes
#include "../dbpp_ListSectionsFlow.h"
#include "../dbpp_SrcNumericalTreatment.h"
 
// Hydraulic radius
//    HydRadius wR; 
 
//      for( int32 j = 0; j < NbSections-1; j++) {
//            Rh = wR( A[j]);
 
//      TermeSf = g * A[j] * n[j]*n[j] * Q[j]*fabs(Q[j])/(A[j]*A[j]*pow (Rh, 4./3.));
 
//        if( aIndex == 0)
//            wTermeS0 = (0.5*g*( mA[aIndex+1] + mA[aIndex])) * (mH[aIndex+1] - mH[aIndex])/mDx;
//        else
//          wTermeS0 = (0.5*g*( mA[aIndex+1] + mA[aIndex-1])) * (mH[aIndex+1] - mH[aIndex-1])/(2.*mDx);
// 
// //       S[j] = (TermeSf - TermeS0);
//          return wTermeS0;
//      }
 
namespace Sfx {
  class FieldLattice;
  class EMcNeilBndCnd;
  class PhysicalSystem;
}
 
namespace Testvs19 {
  template<typename T> class ListSectionsFlow;
}
 
namespace Testvs19
{
  class NujicSrcTermAlgorithm : public dbpp::SrcNumericalTreatment
  {
  public:
    enum class eSrcTermType
    {
      Hderivative = 0, 
      Zderivative = 1  
    };

    enum class eSrcDxType
    {
      central_2ndorder = 0, 
      central_4thorder = 1  
    };
  public:
    NujicSrcTermAlgorithm() : m_derivativeOrder{ 2 }, m_listSectFlow19{nullptr} {}
    NujicSrcTermAlgorithm( const dbpp::ListSectionsFlow& aSectFlow)
    : m_derivativeOrder{ 2 },  // default value
      m_listSectFlow{ aSectFlow }, //
      m_listSectFlow19{nullptr}
    {}

    NujicSrcTermAlgorithm( const Testvs19::ListSectionsFlow<dbpp::SectionFlow>* aSectFlow)
    : m_derivativeOrder{ 2 },  // default value
      m_listSectFlow19{aSectFlow}
    {}

    std::valarray<float64> TraitementTermeSourceS0( const Sfx::scalarField1D& aAfield,
      const Sfx::scalarField1D& aQfield) const;
    std::valarray<float64> TraitementTermeSourceS2( const Sfx::FieldLattice& aAfield, 
      const Sfx::FieldLattice& aHfield) const;
    std::valarray<float64> TraitementTermeSourceS2( const Sfx::PhysicalSystem* aPhysys, const Sfx::EMcNeilBndCnd& aBC);
    void setHderivativeOrder(int32 aDerivOrder) noexcept { m_derivativeOrder = aDerivOrder; }
    int32 getderivativeOrder() const noexcept { return m_derivativeOrder; }
    std::valarray<float64> getManningValues() const noexcept { return std::valarray<float64>(); }
    void setSectionsManningValues( const std::vector<float64>& aVecNvalues) noexcept {/*not implemented yet*/ }
    [[nodiscard]] bool isFictionLess() const noexcept { return true; }
    void setVariableTraitementS2(eSrcTermType aSrcTemrType) noexcept { m_srcType = aSrcTemrType; }  
    void setDxOrderSourceTerm(eSrcDxType aSrcDxType) noexcept { m_dxsrcType = aSrcDxType; }
    eSrcTermType getSrcTermType() const noexcept { return m_srcType; }
    eSrcDxType getSrcDxType() const noexcept { return m_dxsrcType; } 
    void setData( const dbpp::ListSectionsFlow& aData) noexcept { m_listSectFlow = aData; }
  private:
    int m_derivativeOrder;                 
    // Friction energy slope formulae
    std::valarray<float64> m_frictionCoeff; 
    eSrcTermType m_srcType{ eSrcTermType::Hderivative };    
    eSrcDxType m_dxsrcType{ eSrcDxType::central_2ndorder }; 
    dbpp::ListSectionsFlow m_listSectFlow;                  
    const Testvs19::ListSectionsFlow<dbpp::SectionFlow>* m_listSectFlow19; 
  };
} // End of namespace