imhd_to_twofluid.h

Go to the documentation of this file.

#pragma once
 
// wxm includes
#include "apps/wmapplication.h"
 
namespace wxm
{
namespace apps
{
namespace five_moment
{
namespace bc
{
class imhd_to_twofluid : public WmApplication
{
public:
    imhd_to_twofluid();
 
    ~imhd_to_twofluid() override;
 
    void setup(const WxCryptSet& wxc) override;
 
    const std::vector<int>& getInputVariableIndexes(int flag) const override
    {
        return _input_variables;
    }
    const std::vector<int>& getOutputVariableIndexes(int flag) const override
    {
        return _out_variables;
    }
 
    // 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;
 
    void evaluate_function(const real* q_in,
                           const real* aux,
                           const solverVariables_t* pFV,
                           real* q_out) const override;
 
protected:
    std::vector<int> _input_variables;
    std::vector<int> _out_variables;
 
private:
    // mass ratios of two species
    real _Aa, _Ab;
    // charges of two species
    real _Za, _Zb;
    // Za/Zb*Tb/Ta
    real _ztheta_ba;
    real _skin_depth_norm;
 
    // assumptions
    // can be "minimal" or "consistent"
    // "minimal" is based on section 4.1 on Frontiers paper:
    // https://www.frontiersin.org/articles/10.3389/fphy.2018.00105/full
    // "consistent" is based on section 4.2 and is further described in Iman's General
    // Exam
    // This app was originally written only for "minimal" so that is default
    std::string _assumptions;
 
    imhd_to_twofluid& operator=(const imhd_to_twofluid& var);
    imhd_to_twofluid(const imhd_to_twofluid& var);
};
} // namespace bc
} // namespace five_moment
} // namespace apps
} // namespace wxm