DamBreakProb.h

Go to the documentation of this file.

 
#pragma once
 
// C++ include
#include <vector>
 
#include "Sfx_SimulationMgr.h"
// App include
#include "../Nov_VS2019/dbpp_RhsPtr2FuncFlux.h" 
#include "../Nov_VS2019/dbpp_SweRhsAlgorithm.h" 
// SfxBase19 library include
#include "include/Sfx_DefineTypes.h"
 
namespace dbpp
{
  // alias 
  using Time = float64;

  class DamBreakProblem
  {
  public:
    enum class eProbType
    {
      ProblemA = 0,     
      ProblemB = 1      
    };
 
  public:

    DamBreakProblem( float64 aPhi0, float64 aPhi1);
    float64 computeVelocity( float64 aCoord, float64 aTime);
    std::vector<float64> computeVelocity( const std::vector<float64>& aListofCoord, float64 aTime);
    float64 computeWaterDepth(float64 aCoord, float64 aTime);
    std::vector<float64> computeWaterDepth( const std::vector<float64>& aListofCoord, float64 aTime);
    void setIC(eProbType aTestProblem = eProbType::ProblemA);
    void timeLoop() 
    {
      /*not implemented yet*/ 
      //while (SfxTimePrm) {/*advance simulation by step*/ }
    }

    void init() 
    {
      // create solver(this) (ctor from DamBreakProblem)
      // create rhs algorithm (factory) and set from config file
      // set solver rhs discretization
      // set simulation start and stop time SfxTimePrm
    }

    virtual SweRhsAlgorithm* getRHSdiscretization() 
    { 
      auto w_calculFF = Sfx::SimulationMgr::getSingleton().getPtr2FAlgorithm(std::string{ "Whatever" });
      return new RhsPtr2FuncFlux{ std::string{"DamBreakProblem"}, w_calculFF, new SrcNumericalTreatment{}, false };
    }
    // Plot profile for velocity and water depth
//      void plotCurve();
 
  private:
    float64 mS;                
    float64 mU2;               
    float64 mC2;               
    float64 mSqrtG;             
    float64 mPhi1;             
    float64 mPhi0;              
    eProbType mTestProblem;    // Type of the initial profile */
    std::vector<float64> mU0;  // Initial water velocity */
    std::vector<float64> mH0;  // Initial water level */
 
    // ODESolver* m_odeSlvr1D;
     
    // Math equations that model the wave propagation 
    //  StVenant1DSys* mCLSystem;
 
    // Initial condition to solve the DAM-BREAK problem
    // InitialCnd mInitialCnd;
 
    // Boundary condition
    //  BoundaryCnd mBndCnd;
 };
} // End of namespace