flowthrough_wall_gradients.h

Go to the documentation of this file.

#pragma once
 
// Wm includes
#include "apps/wmapplication.h"
 
// STL includes
#include <memory>
 
namespace wxm
{
namespace apps
{
namespace rmhd
{
namespace bc
{
/*
 *
 * Sets BC for Flowthrough endwall gradients
 * From r=0 to r=0.03 and from r = 0.05 to r = 0.1, use conducting wall
 * From r = 0.03 to r = 0.05 use an outflow
 * Transition between the two using a sine wave function to mitigate possible numerical
 * instability
 * Assume the wall normal is z (x in 2d axisymmetric) direction
 *
 */
class FlowthroughWallGradients : public WmApplication
{
public:
    FlowthroughWallGradients();
    ~FlowthroughWallGradients();
    void setup(const WxCryptSet& wxc) override;
 
    const std::vector<int>& getInputVariableIndexes(int flag) const override
    {
        return _grad_q_idcs;
    }
 
    const std::vector<int>& getOutputVariableIndexes(int flag) const override
    {
        return _grad_q_idcs;
    }
 
    const std::vector<int>& getAuxiliaryVariableIndexes(int flag) const override
    {
        return _q_idcs;
    }
 
    // redefinition of bc
    void bc_q(const real* q_in,
              const real* aux_in,
              const real* aux_out,
              const solverVariables_t* pFV,
              real* q_out) const override;
 
protected:
    real _gas_gamma;
 
    real _r1;
    real _r2;
    real _w;
    real _w_over_2;
 
    std::vector<int> _grad_q_idcs;
    std::vector<int> _q_idcs;
 
    // pointer to subapplication
    std::unique_ptr<WmApplication> _cw_app;
    std::unique_ptr<WmApplication> _ft_app;
 
private:
    FlowthroughWallGradients& operator=(const FlowthroughWallGradients& var);
    FlowthroughWallGradients(const FlowthroughWallGradients& var);
};
 
} // namespace bc
} // namespace rmhd
} // namespace apps
} // namespace wxm