/*===============================================================
                        Introduction
  This is the computation of J.A.Beattie - O.C.Bridgeman's model,
which is revised from ver de Waals equation. The form of this 
equation is like this:
		
	p=R*T[1-c/(V*T^3)]/V^2*[V+B0*(1-b/V)]-A0*(1-a/V)/V^2

This equation is much acute than the RK equation. But it is full
of complicated items and requires more original data. This 
equation can be still be very accurate when the special volume is 
more than two times of that on critical point.
---------------------------------------------------------------*/

#pragma once

#ifndef BBGas
#define BBGas
#endif

#include "GasBase.h"

namespace XZGas
{  
    class BBGasModel:public GasModel
    {
    public:
        BBGasModel(const GasData &gas);
        ~BBGasModel();

        inline double Cp(double TT,double pp);
		inline double Cv(double TT,double vv);
	
		inline double Gama(double TT,double pp);
		inline double h(double TT,double pp);
		inline double u(double TT,double vv);
		inline double s(double TT,double pp);
		
        inline double p(double TT,double vv);
        inline double v(double TT,double pp);
        inline double T(double pp,double vv);	

		inline double Alphav(double pp,double vv);
		inline double KT(double TT,double vv);
		inline double Beta(double pp,double vv);
    
		inline void ChangeGas(const GasData &newgas);

    protected:
        inline double dpdvT(double TT,double vv);
		inline double dpdTv(double TT,double vv);
		
	private:
		double A0;
        double a;
        double B0;
        double b;
        double c;                   
    };
}