kh_instability_kinetic.h

Go to the documentation of this file.

#ifndef wmapplication_kh_instability_kinetic_fluid_vars_h
#define wmapplication_kh_instability_kinetic_fluid_vars_h
 
// WARPXM includes
#include "apps/wmapplication.h"
 
namespace wxm
{
namespace apps
{
namespace functions
{
namespace five_moment
{
class Kinetics2d2vKelvinHelmholtzInstabilityFluidVariables : public WmApplication
{
public:
    Kinetics2d2vKelvinHelmholtzInstabilityFluidVariables();
 
    ~Kinetics2d2vKelvinHelmholtzInstabilityFluidVariables();
 
    void setup(const WxCryptSet& wxc) override;
 
    const std::vector<int>& getOutputVariableIndexes(int flag) const override
    {
        return _output_fluid;
    }
 
    void evaluate_function(const real* q,
                           const real* aux,
                           const solverVariables_t* pSV,
                           real* result) const override;
 
    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:
    std::vector<int> _output_fluid;
 
    real _gamma;
    real _mass;
    real _d_width; // shear_layer_width
 
    real _perturbation_wavelength;
    real _perturbation_amplitude;
    std::vector<real> _perturbation_center;
 
    std::vector<real> _n_fit_xvec;
    std::vector<real> _n_fit_acoeffs;
    std::vector<real> _n_fit_bcoeffs;
 
    std::vector<real> _vy_fit_xvec;
    std::vector<real> _vy_fit_acoeffs;
    std::vector<real> _vy_fit_bcoeffs;
 
    std::vector<real> _Ty_fit_xvec;
    std::vector<real> _Ty_fit_acoeffs;
    std::vector<real> _Ty_fit_bcoeffs;
 
    mutable std::map<real, real> _n_fit_cache;
    mutable std::map<real, real> _vy_fit_cache;
    mutable std::map<real, real> _Ty_fit_cache;
 
    real _Tx;
 
private:
    Kinetics2d2vKelvinHelmholtzInstabilityFluidVariables&
    operator=(const Kinetics2d2vKelvinHelmholtzInstabilityFluidVariables& var);
    Kinetics2d2vKelvinHelmholtzInstabilityFluidVariables(
        const Kinetics2d2vKelvinHelmholtzInstabilityFluidVariables& var);
};
 
} // namespace five_moment
} // namespace functions
} // namespace apps
} // namespace wxm
 
#endif