Commit 9bbc8924 authored by Jan Schäfer's avatar Jan Schäfer
Browse files

AUSMPlus special BC

parent f9cf85e1

src/Examples/flows/BoundaryConditions/Boiler/AUSMPlusBoundaryConditionsBoiler.h

0 → 100644
+153 −0
Original line number Diff line number Diff line 
#include <TNL/Functions/FunctionAdapter.h>



#include "DensityBoundaryConditionBoiler.h"

#include "MomentumXBoundaryConditionBoiler.h"

#include "MomentumYBoundaryConditionBoiler.h"

#include "MomentumZBoundaryConditionBoiler.h"

#include "AUSMPlusEnergyBoundaryConditionBoiler.h"



namespace TNL {



template< typename Mesh,

          typename Function,

          typename Real = typename Mesh::RealType,

          typename Index = typename Mesh::IndexType >

class BoundaryConditionsBoiler

{

   public:

      typedef Mesh MeshType;

      typedef Real RealType;

      typedef Index IndexType;

      typedef Function FunctionType;

      typedef Functions::MeshFunction< Mesh > MeshFunctionType;

      typedef typename Mesh::DeviceType DeviceType;



      typedef TNL::Operators::DensityBoundaryConditionsBoiler< MeshType, FunctionType, RealType, IndexType > DensityBoundaryConditionsType;

      typedef TNL::Operators::MomentumXBoundaryConditionsBoiler< MeshType, FunctionType, RealType, IndexType > MomentumXBoundaryConditionsType;

      typedef TNL::Operators::MomentumYBoundaryConditionsBoiler< MeshType, FunctionType, RealType, IndexType > MomentumYBoundaryConditionsType;

      typedef TNL::Operators::MomentumZBoundaryConditionsBoiler< MeshType, FunctionType, RealType, IndexType > MomentumZBoundaryConditionsType;

      typedef TNL::Operators::AUSMPlusEnergyBoundaryConditionsBoiler< MeshType, FunctionType, RealType, IndexType > EnergyBoundaryConditionsType;

      typedef CompressibleConservativeVariables< MeshType > CompressibleConservativeVariablesType;



      typedef Pointers::SharedPointer< DensityBoundaryConditionsType > DensityBoundaryConditionsTypePointer;

      typedef Pointers::SharedPointer< MomentumXBoundaryConditionsType > MomentumXBoundaryConditionsTypePointer;

      typedef Pointers::SharedPointer< MomentumYBoundaryConditionsType > MomentumYBoundaryConditionsTypePointer;

      typedef Pointers::SharedPointer< MomentumZBoundaryConditionsType > MomentumZBoundaryConditionsTypePointer;

      typedef Pointers::SharedPointer< EnergyBoundaryConditionsType > EnergyBoundaryConditionsTypePointer;

      typedef Pointers::SharedPointer< CompressibleConservativeVariablesType > CompressibleConservativeVariablesPointer;

      typedef Pointers::SharedPointer< MeshType > MeshPointer;

      typedef Pointers::SharedPointer< MeshFunctionType, DeviceType > MeshFunctionPointer;



      static void configSetup( Config::ConfigDescription& config,

                               const String& prefix = "" )

      {

      }



      bool setup( const MeshPointer& meshPointer,

                  const Config::ParameterContainer& parameters,

                  const String& prefix = "" )

      {

         this->densityBoundaryConditionsPointer->setup( meshPointer, parameters, prefix);

         this->momentumXBoundaryConditionsPointer->setup( meshPointer, parameters, prefix);

         this->momentumYBoundaryConditionsPointer->setup( meshPointer, parameters, prefix);

         this->momentumZBoundaryConditionsPointer->setup( meshPointer, parameters, prefix);

         this->energyBoundaryConditionsPointer->setup( meshPointer, parameters, prefix);

         return true;

      }



      void setCompressibleConservativeVariables(const CompressibleConservativeVariablesPointer& compressibleConservativeVariables)

      {

         this->densityBoundaryConditionsPointer->setCompressibleConservativeVariables(compressibleConservativeVariables);

         this->momentumXBoundaryConditionsPointer->setCompressibleConservativeVariables(compressibleConservativeVariables);

         this->momentumYBoundaryConditionsPointer->setCompressibleConservativeVariables(compressibleConservativeVariables);

         this->momentumZBoundaryConditionsPointer->setCompressibleConservativeVariables(compressibleConservativeVariables);

         this->energyBoundaryConditionsPointer->setCompressibleConservativeVariables(compressibleConservativeVariables);

      }



      void setTimestep(const RealType timestep)

      {

         this->densityBoundaryConditionsPointer->setTimestep(timestep);

         this->momentumXBoundaryConditionsPointer->setTimestep(timestep);

         this->momentumYBoundaryConditionsPointer->setTimestep(timestep);

         this->momentumZBoundaryConditionsPointer->setTimestep(timestep);

         this->energyBoundaryConditionsPointer->setTimestep(timestep);   

      }



      void setGamma(const RealType gamma)

      {

         this->densityBoundaryConditionsPointer->setGamma(gamma);

         this->momentumXBoundaryConditionsPointer->setGamma(gamma);

         this->momentumYBoundaryConditionsPointer->setGamma(gamma);

         this->momentumZBoundaryConditionsPointer->setGamma(gamma);

         this->energyBoundaryConditionsPointer->setGamma(gamma);

      }



      void setPressure(const MeshFunctionPointer& pressure)

      {

         this->densityBoundaryConditionsPointer->setPressure(pressure);

         this->momentumXBoundaryConditionsPointer->setPressure(pressure);

         this->momentumYBoundaryConditionsPointer->setPressure(pressure);

         this->momentumZBoundaryConditionsPointer->setPressure(pressure);

         this->energyBoundaryConditionsPointer->setPressure(pressure);

      }



      void setSpeed(const RealType cavitySpeed)

      {

         this->momentumXBoundaryConditionsPointer->setCavitySpeed(cavitySpeed);

         this->momentumYBoundaryConditionsPointer->setCavitySpeed(cavitySpeed);

         this->momentumZBoundaryConditionsPointer->setCavitySpeed(cavitySpeed);

         this->energyBoundaryConditionsPointer->setCavitySpeed(cavitySpeed);

      }



      void setHorizontalThrottleSpeed(const RealType startSpeed,

                                      const RealType finalSpeed,

                                      const RealType time,

                                      const RealType finalTime )

      {



      }



      void setVerticalThrottleSpeed(const RealType startSpeed,

                                    const RealType finalSpeed,

                                    const RealType time,

                                    const RealType finalTime )

      {



      }



      DensityBoundaryConditionsTypePointer& getDensityBoundaryCondition()

      {

         return this->densityBoundaryConditionsPointer;

      }



      MomentumXBoundaryConditionsTypePointer& getMomentumXBoundaryCondition()

      {

         return this->momentumXBoundaryConditionsPointer;

      }



      MomentumYBoundaryConditionsTypePointer& getMomentumYBoundaryCondition()

      {

         return this->momentumYBoundaryConditionsPointer;

      }



      MomentumZBoundaryConditionsTypePointer& getMomentumZBoundaryCondition()

      {

         return this->momentumZBoundaryConditionsPointer;

      }



      EnergyBoundaryConditionsTypePointer& getEnergyBoundaryCondition()

      {

         return this->energyBoundaryConditionsPointer;

      }





   protected:

      DensityBoundaryConditionsTypePointer densityBoundaryConditionsPointer;

      MomentumXBoundaryConditionsTypePointer momentumXBoundaryConditionsPointer;

      MomentumYBoundaryConditionsTypePointer momentumYBoundaryConditionsPointer;

      MomentumZBoundaryConditionsTypePointer momentumZBoundaryConditionsPointer;

      EnergyBoundaryConditionsTypePointer energyBoundaryConditionsPointer;



};



} //namespace TNL

src/Examples/flows/BoundaryConditions/Boiler/AUSMPlusEnergyBoundaryConditionBoiler.h

0 → 100644
+853 −0

File added.

Preview size limit exceeded, changes collapsed.

src/Examples/flows/BoundaryConditions/BoilerModel/AUSMPlusBoundaryConditionsBoilerModel.h

0 → 100644
+194 −0
Original line number Diff line number Diff line 
#include <TNL/Functions/FunctionAdapter.h>



#include "DensityBoundaryConditionBoilerModel.h"

#include "MomentumXBoundaryConditionBoilerModel.h"

#include "MomentumYBoundaryConditionBoilerModel.h"

#include "MomentumZBoundaryConditionBoilerModel.h"

#include "AUSMPlusEnergyBoundaryConditionBoilerModel.h"



namespace TNL {



template< typename Mesh,

          typename Function,

          typename Real = typename Mesh::RealType,

          typename Index = typename Mesh::IndexType >

class BoundaryConditionsBoilerModel

{

   public:

      typedef Mesh MeshType;

      typedef Real RealType;

      typedef Index IndexType;

      typedef Function FunctionType;

      typedef Functions::MeshFunction< Mesh > MeshFunctionType;

      typedef typename Mesh::DeviceType DeviceType;



      typedef TNL::Operators::DensityBoundaryConditionsBoilerModel< MeshType, FunctionType, RealType, IndexType > DensityBoundaryConditionsType;

      typedef TNL::Operators::MomentumXBoundaryConditionsBoilerModel< MeshType, FunctionType, RealType, IndexType > MomentumXBoundaryConditionsType;

      typedef TNL::Operators::MomentumYBoundaryConditionsBoilerModel< MeshType, FunctionType, RealType, IndexType > MomentumYBoundaryConditionsType;

      typedef TNL::Operators::MomentumZBoundaryConditionsBoilerModel< MeshType, FunctionType, RealType, IndexType > MomentumZBoundaryConditionsType;

      typedef TNL::Operators::AUSMPlusEnergyBoundaryConditionsBoilerModel< MeshType, FunctionType, RealType, IndexType > EnergyBoundaryConditionsType;

      typedef CompressibleConservativeVariables< MeshType > CompressibleConservativeVariablesType;



      typedef Pointers::SharedPointer< DensityBoundaryConditionsType > DensityBoundaryConditionsTypePointer;

      typedef Pointers::SharedPointer< MomentumXBoundaryConditionsType > MomentumXBoundaryConditionsTypePointer;

      typedef Pointers::SharedPointer< MomentumYBoundaryConditionsType > MomentumYBoundaryConditionsTypePointer;

      typedef Pointers::SharedPointer< MomentumZBoundaryConditionsType > MomentumZBoundaryConditionsTypePointer;

      typedef Pointers::SharedPointer< EnergyBoundaryConditionsType > EnergyBoundaryConditionsTypePointer;

      typedef Pointers::SharedPointer< CompressibleConservativeVariablesType > CompressibleConservativeVariablesPointer;

      typedef Pointers::SharedPointer< MeshType > MeshPointer;

      typedef Pointers::SharedPointer< MeshFunctionType, DeviceType > MeshFunctionPointer;



      static void configSetup( Config::ConfigDescription& config,

                               const String& prefix = "" )

      {

         config.addEntry< double >( prefix + "vertical-angle", "Vertical angle of throttle in degrees", 0 );

         config.addEntry< double >( prefix + "horizontal-angle", "Horizontal angle of throttle in degrees", 45 );

      }



      bool setup( const MeshPointer& meshPointer,

                  const Config::ParameterContainer& parameters,

                  const String& prefix = "" )

      {

         this->densityBoundaryConditionsPointer->setup( meshPointer, parameters, prefix);

         this->momentumXBoundaryConditionsPointer->setup( meshPointer, parameters, prefix);

         this->momentumYBoundaryConditionsPointer->setup( meshPointer, parameters, prefix);

         this->momentumZBoundaryConditionsPointer->setup( meshPointer, parameters, prefix);

         this->energyBoundaryConditionsPointer->setup( meshPointer, parameters, prefix);

         setZAngle(parameters.getParameter< double >( prefix + "vertical-angle" ) * M_PI / 180.0 );

         setXYAngle(parameters.getParameter< double >( prefix + "horizontal-angle" ) * M_PI / 180.0 );

         return true;

      }



      void setCompressibleConservativeVariables(const CompressibleConservativeVariablesPointer& compressibleConservativeVariables)

      {

         this->densityBoundaryConditionsPointer->setCompressibleConservativeVariables(compressibleConservativeVariables);

         this->momentumXBoundaryConditionsPointer->setCompressibleConservativeVariables(compressibleConservativeVariables);

         this->momentumYBoundaryConditionsPointer->setCompressibleConservativeVariables(compressibleConservativeVariables);

         this->momentumZBoundaryConditionsPointer->setCompressibleConservativeVariables(compressibleConservativeVariables);

         this->energyBoundaryConditionsPointer->setCompressibleConservativeVariables(compressibleConservativeVariables);

      }



      void setTimestep(const RealType timestep)

      {

         this->densityBoundaryConditionsPointer->setTimestep(timestep);

         this->momentumXBoundaryConditionsPointer->setTimestep(timestep);

         this->momentumYBoundaryConditionsPointer->setTimestep(timestep);

         this->momentumZBoundaryConditionsPointer->setTimestep(timestep);

         this->energyBoundaryConditionsPointer->setTimestep(timestep);   

      }



      void setGamma(const RealType gamma)

      {

         this->densityBoundaryConditionsPointer->setGamma(gamma);

         this->momentumXBoundaryConditionsPointer->setGamma(gamma);

         this->momentumYBoundaryConditionsPointer->setGamma(gamma);

         this->momentumZBoundaryConditionsPointer->setGamma(gamma);

         this->energyBoundaryConditionsPointer->setGamma(gamma);

      }



      void setZAngle(const RealType zAngle)

      {

         this->densityBoundaryConditionsPointer->setZAngle(zAngle);

         this->momentumXBoundaryConditionsPointer->setZAngle(zAngle);

         this->momentumYBoundaryConditionsPointer->setZAngle(zAngle);

         this->momentumZBoundaryConditionsPointer->setZAngle(zAngle);

         this->energyBoundaryConditionsPointer->setZAngle(zAngle);

      }



      void setXYAngle(const RealType xYAngle)

      {

         this->densityBoundaryConditionsPointer->setXYAngle(xYAngle);

         this->momentumXBoundaryConditionsPointer->setXYAngle(xYAngle);

         this->momentumYBoundaryConditionsPointer->setXYAngle(xYAngle);

         this->momentumZBoundaryConditionsPointer->setXYAngle(xYAngle);

         this->energyBoundaryConditionsPointer->setXYAngle(xYAngle);

      }



      void setPressure(const MeshFunctionPointer& pressure)

      {

         this->densityBoundaryConditionsPointer->setPressure(pressure);

         this->momentumXBoundaryConditionsPointer->setPressure(pressure);

         this->momentumYBoundaryConditionsPointer->setPressure(pressure);

         this->momentumZBoundaryConditionsPointer->setPressure(pressure);

         this->energyBoundaryConditionsPointer->setPressure(pressure);

      }



      void setSpeed(const RealType cavitySpeed)

      {



      }



      void setHorizontalThrottleSpeed(const RealType startSpeed,

                                      const RealType finalSpeed,

                                      const RealType time,

                                      const RealType finalTime )

      {

         RealType horizontalThrottleSpeed = 0;

         if(time <= finalTime)

            if( time != 0 && finalTime != 0 )

               horizontalThrottleSpeed = startSpeed + ( finalSpeed - startSpeed ) * ( time / finalTime );

            else

               horizontalThrottleSpeed = 0;

         else

            horizontalThrottleSpeed = finalSpeed;

         this->momentumXBoundaryConditionsPointer->setHorizontalThrottleSpeed(horizontalThrottleSpeed);

         this->momentumYBoundaryConditionsPointer->setHorizontalThrottleSpeed(horizontalThrottleSpeed);

         this->momentumZBoundaryConditionsPointer->setHorizontalThrottleSpeed(horizontalThrottleSpeed);

         this->energyBoundaryConditionsPointer->setHorizontalThrottleSpeed(horizontalThrottleSpeed);

      }



      void setVerticalThrottleSpeed(const RealType startSpeed,

                                    const RealType finalSpeed,

                                    const RealType time,

                                    const RealType finalTime )

      {

         RealType verticalThrottleSpeed = 0;

         if(time <= finalTime)

            if( time != 0 && finalTime != 0 )

               verticalThrottleSpeed = startSpeed + ( finalSpeed - startSpeed ) * ( time / finalTime );

            else

               verticalThrottleSpeed = 0;

         else

            verticalThrottleSpeed = finalSpeed;

         this->momentumXBoundaryConditionsPointer->setVerticalThrottleSpeed(verticalThrottleSpeed);

         this->momentumYBoundaryConditionsPointer->setVerticalThrottleSpeed(verticalThrottleSpeed);

         this->momentumZBoundaryConditionsPointer->setVerticalThrottleSpeed(verticalThrottleSpeed);

         this->energyBoundaryConditionsPointer->setVerticalThrottleSpeed(verticalThrottleSpeed);

      }



      DensityBoundaryConditionsTypePointer& getDensityBoundaryCondition()

      {

         return this->densityBoundaryConditionsPointer;

      }



      MomentumXBoundaryConditionsTypePointer& getMomentumXBoundaryCondition()

      {

         return this->momentumXBoundaryConditionsPointer;

      }



      MomentumYBoundaryConditionsTypePointer& getMomentumYBoundaryCondition()

      {

         return this->momentumYBoundaryConditionsPointer;

      }



      MomentumZBoundaryConditionsTypePointer& getMomentumZBoundaryCondition()

      {

         return this->momentumZBoundaryConditionsPointer;

      }



      EnergyBoundaryConditionsTypePointer& getEnergyBoundaryCondition()

      {

         return this->energyBoundaryConditionsPointer;

      }





   protected:

      DensityBoundaryConditionsTypePointer densityBoundaryConditionsPointer;

      MomentumXBoundaryConditionsTypePointer momentumXBoundaryConditionsPointer;

      MomentumYBoundaryConditionsTypePointer momentumYBoundaryConditionsPointer;

      MomentumZBoundaryConditionsTypePointer momentumZBoundaryConditionsPointer;

      EnergyBoundaryConditionsTypePointer energyBoundaryConditionsPointer;



};



} //namespace TNL

src/Examples/flows/BoundaryConditions/BoilerModel/AUSMPlusEnergyBoundaryConditionBoilerModel.h

0 → 100644
+1072 −0

File added.

Preview size limit exceeded, changes collapsed.

src/Examples/flows/BoundaryConditions/BoilerModel/DensityBoundaryConditionBoilerModel.h

+10 −19
Original line number Diff line number Diff line
@@ -265,20 +265,11 @@ class DensityBoundaryConditionsBoilerModel< Meshes::Grid< 2, MeshReal, Device, M 
         const auto& neighborEntities = entity.getNeighborEntities();

         const IndexType& index = entity.getIndex();

         if( entity.getCoordinates().x() == 0 )

         {

            return u[ neighborEntities.template getEntityIndex< 0, 0 >() ];

         }

            return u[ neighborEntities.template getEntityIndex< 1, 0 >() ];

         if( entity.getCoordinates().x() == entity.getMesh().getDimensions().x() - 1 )

         {

            if (entity.getCoordinates().y() < 0.835 * ( entity.getMesh().getDimensions().y() - 1 ))

               return u[ neighborEntities.template getEntityIndex< 0, 0 >() ];

            else

            return u[ neighborEntities.template getEntityIndex< -1, 0 >() ];

         }

         if( entity.getCoordinates().y() == 0 )

         {

            return u[ neighborEntities.template getEntityIndex< 0, 0 >() ];

         }

            return u[ neighborEntities.template getEntityIndex< 0, 1 >() ];

         // The following line is commented to avoid compiler warning

         //if( entity.getCoordinates().y() == entity.getMesh().getDimensions().y() - 1 )

         {
@@ -427,7 +418,7 @@ class DensityBoundaryConditionsBoilerModel< Meshes::Grid< 3, MeshReal, Device, M 
         const IndexType& index = entity.getIndex();

         if( entity.getCoordinates().x() == 0 )

         {

            return u[ neighborEntities.template getEntityIndex< 0, 0, 0 >() ];

            return u[ neighborEntities.template getEntityIndex< 1, 0, 0 >() ];

         }

         if( entity.getCoordinates().x() == entity.getMesh().getDimensions().x() - 1 )

         {
@@ -435,24 +426,24 @@ class DensityBoundaryConditionsBoilerModel< Meshes::Grid< 3, MeshReal, Device, M 
            if( entity.getCoordinates().z() > 0.835 * ( entity.getMesh().getDimensions().z() - 1 ) )

                 return u[ neighborEntities.template getEntityIndex< -1, 0, 0 >() ];

            else

               return u[ neighborEntities.template getEntityIndex< 0, 0, 0 >() ];

               return u[ neighborEntities.template getEntityIndex< -1, 0, 0 >() ];

         }

         if( entity.getCoordinates().y() == 0 )

         {

            return u[ neighborEntities.template getEntityIndex< 0, 0, 0 >() ];

            return u[ neighborEntities.template getEntityIndex< 0, 1, 0 >() ];

         }

         if( entity.getCoordinates().y() == entity.getMesh().getDimensions().y() - 1 )

         {

            return u[ neighborEntities.template getEntityIndex< 0, 0, 0 >() ];

            return u[ neighborEntities.template getEntityIndex< 0, -1, 0 >() ];

         }

         if( entity.getCoordinates().z() == 0 )

         {

            return u[ neighborEntities.template getEntityIndex< 0, 0, 0 >() ];

            return u[ neighborEntities.template getEntityIndex< 0, 0, 1 >() ];

         }

         // The following line is commented to avoid compiler warning

         //if( entity.getCoordinates().z() == entity.getMesh().getDimensions().z() - 1 )

         {

            return u[ neighborEntities.template getEntityIndex< 0, 0, 0 >() ];

            return u[ neighborEntities.template getEntityIndex< 0, 0, -1 >() ];

         }   

      }