intraspecies_cyl_source.h

Go to the documentation of this file.

#pragma once
 
// Wm includes
#include "apps/wmapplication.h"
 
#include "apps/mhd_two_temperature/app_setup_helper.h"
 
// STL includes
#include <memory>
 
namespace wxm
{
namespace apps
{
namespace mhd_two_temperature
{
namespace viscousmhd
{
class IntraspeciesCyl : public WmApplication
{
public:
    IntraspeciesCyl();
 
    ~IntraspeciesCyl() override;
 
    void setup(const WxCryptSet& wxc) override;
 
    const std::vector<int>& getInputVariableIndexes(int flag) const override
    {
        return q_indexes_;
    }
    const std::vector<int>& getAuxiliaryVariableIndexes(int flag) const override
    {
        return grad_q_indexes_;
    }
    const std::vector<int>& getOutputVariableIndexes(int flag) const override
    {
        return q_indexes_;
    }
 
    real source(const real* q,
                const real* aux,
                const elementGeometry_t* pEG,
                real* source) const override;
 
protected:
    std::unique_ptr<wxm::apps::mhd_two_temperature::AppSetupHelper> setup_helper_;
 
    real gas_gamma_;
 
    real Zi_;
 
    real Ai_;
 
    real const_ion_mu_;
 
    real const_ion_k_;
    real const_electron_k_;
 
    real min_density_floor_;
    real min_pressure_floor_;
 
    int radial_index_;
    int ir_, ith_, iz_;
 
    std::vector<int> q_indexes_;
    std::vector<int> grad_q_indexes_;
 
private:
    IntraspeciesCyl& operator=(const IntraspeciesCyl& var) = delete;
    IntraspeciesCyl(const IntraspeciesCyl& var) = delete;
};
 
} // namespace viscousmhd
} // namespace mhd_two_temperature
} // namespace apps
} // namespace wxm