Sfx_ImposeBnd.h

Go to the documentation of this file.

#pragma once
 
// SfxBase19 Library includes
#include "Sfx/Sfx_NodalValues.h"
#include "Sfx/Sfx_FieldLattice.h"
// App include
#include "../Nov_VS2019/dbpp_SectionFlow.h"
 
// forward declaration
namespace Sfx { 
    class PhysicalSystem; 
    class StateVector;
}
 
namespace Sfx
{
    class EMcNeilBndCnd
    {
    public:
        using pair_node = std::pair<Nodal_Value,Nodal_Value>;     // deprecated
 
    public:
    // DESIGN NOTE
    /*
         * @brief Its important to note that B.C. (characteristic) is applied
     *  to the physical system. This is a prototype version impl. 
        */
    void applyBC( const Sfx::PhysicalSystem* aPhySys, double aDt);

        pair_node getBC() const { return std::make_pair(m_upstream,m_downstream);}
        void setBC( Sfx::NumArray<float64>& aU1, Sfx::NumArray<float64>& aU2);
    private:
        //std::set<dbpp::SectionFlow> m_setSectFlow;  totally ridiculus!!!
        Nodal_Value m_upstream;      
        Nodal_Value m_downstream;    
    private:
        struct UpstreamCelerity
        {
            float64 CN{};
            float64 CB{};
        };

        struct DownstreamCelerity
        {
            float64 CA{};
            float64 CP{};
        };

      DownstreamCelerity computeAval( const Sfx::StateVector& aVal,
                          const dbpp::SectionFlow& aSectAval,
                          float64 aDt);
        UpstreamCelerity computeAmont( const Sfx::StateVector& aVal,
                           const dbpp::SectionFlow& aSectAval,
                           float64 aDt);
    };
} // End of namespace