MHDTests.hpp

Go to the documentation of this file.

// Class to provide various test cases for the Stokes problem
//
//
// Author: Liam Yemm (liam.yemm@monash.edu)
//
 
#include <mesh.hpp>
// #include <Math.hpp>
#include "basis.hpp"
 
using Math::PI;
 
#ifndef _MHD_TESTS_HPP
#define _MHD_TESTS_HPP
 
 
namespace HArDCore3D
{
 
    class MHDTests
    {
    public:
        MHDTests(
            size_t u_id, 
            size_t b_id,  
            size_t p_id,  
            double visc,  
            double diff  
        );
 
        FType<VectorRd> velocity();
 
        FType<MatrixRd> grad_velocity();
 
        FType<VectorRd> laplace_velocity();
 
        FType<VectorRd> magnetic();
 
        FType<MatrixRd> grad_magnetic();
 
        FType<VectorRd> laplace_magnetic();
 
        FType<VectorRd> u_dot_grad_u();
        FType<VectorRd> u_dot_grad_b();
        FType<VectorRd> b_dot_grad_u();
        FType<VectorRd> b_dot_grad_b();
 
        FType<double> pressure();
 
        FType<VectorRd> grad_pressure();
 
        FType<VectorRd> magnetic_source();
        FType<VectorRd> velocity_source();
 
    private:
        size_t m_u_id;
        size_t m_b_id;
        size_t m_p_id;
        double m_visc;
        double m_diff;
 
        void validate();
    };
}
 
 
#endif