/*
Properties for R32.  
by Ian Bell


Thermo props from
Tillner-Roth, R. and Yokozeki, A.,
"An international standard equation of state for difluoromethane (R-32)
for temperatures from the triple point at 136.34 K to 435 K and pressures
up to 70 MPa,"
J. Phys. Chem. Ref. Data, 26(6):1273-1328, 1997.

*/

#if defined(_MSC_VER)
#define _CRTDBG_MAP_ALLOC
#define _CRT_SECURE_NO_WARNINGS
#include <stdlib.h>
#include <crtdbg.h>
#else
#include <stdlib.h>
#endif

#include <math.h>
#include "string.h"
#include "stdio.h"
#include <stdlib.h>
#include "CoolProp.h"

static const double a[]={
	 0.0,			//[0]
	 0.1046634e+1, 	//[1]
	-0.5451165,		//[2]
	-0.2448595e-2,	//[3]
	-0.4877002e-1,	//[4]
	 0.3520158e-1,	//[5]
	 0.1622750e-2,	//[6]
	 0.2377225e-4,	//[7]
	 0.29149e-1,	//[8]
	 0.3386203e-2,	//[9]
	-0.4202444e-2,	//[10]
	 0.4782025e-3,	//[11]
	-0.5504323e-2,	//[12]
	-0.2418396e-1,	//[13]
	 0.4209034,		//[14]
 	-0.4616537,		//[15]
	-0.1200513e+1,	//[16]
	-0.2591550e+1,	//[17]
	-0.1400145e+1,	//[18]
	 0.8263017		//[19]
};

static const int d[]={
	0,			//[0]
	1, 			//[1]
	2, 			//[2]
	5, 			//[3]
	1, 			//[4]
	1, 			//[5]
	3, 			//[6]
	8, 			//[7]
	4, 			//[8]
	4, 			//[9]
	4, 			//[10]
	8, 			//[11]
	3, 			//[12]
	5, 			//[13]
	1, 			//[14]
	1, 			//[15]
	3, 			//[16]
	1, 			//[17]
	2, 			//[18]
	3 			//[19]
};

static const double t[]={
	0.0,		//[0]
	1.0/4.0,	//[1]
	1.0,		//[2]
	-1.0/4.0,	//[3]
	-1.0,		//[4]
	2.0,		//[5]
	2.0,		//[6]
	3.0/4.0,	//[7]
	1.0/4.0,	//[8]
	18.0,		//[9]
	26.0,		//[10]
	-1.0,		//[11]
	25.0, 		//[12]
	7.0/4.0,	//[13]
	4.0,		//[14]
	5.0,		//[15]
	1.0,		//[16]
	3.0/2.0,	//[17]
	1.0,		//[18]
	1.0/2.0		//[19]
};

static const int c[]={
	0,			//[0]
	0,			//[1]
	0,			//[2]
	0,			//[3]
	0,			//[4]
	0,			//[5]
	0,			//[6]
	0,			//[7]
	0,			//[8]
	4,			//[9]
	3,			//[10]
	1,			//[11]
	4,			//[12]
	1,			//[13]
	2,			//[14]
	2,			//[15]
	1,			//[16]
	1,			//[17]
	1,			//[18]
	1			//[19]
};

static const double a0[]={
	-8.258096,	//[0]
	6.353098,	//[1]
	3.004486,	//[2]
	1.160761,	//[3]
	2.645151,	//[4]
	5.794987,	//[5]
	1.129475	//[6]
};
static const double n0[]={
	0.0,		//[0]
	0.0,		//[1]
	0.0,		//[2]
	2.2718538,	//[3]
	11.9144210,	//[4]
	5.1415638,	//[5]
	32.7682170	//[6]
};

R32Class::R32Class()
{
	std::vector<double> n_v(a,a+sizeof(a)/sizeof(double));
	std::vector<double> d_v(d,d+sizeof(d)/sizeof(int));
	std::vector<double> t_v(t,t+sizeof(t)/sizeof(double));
	std::vector<double> l_v(c,c+sizeof(c)/sizeof(int));
	std::vector<double> a0_v(a0,a0+sizeof(a0)/sizeof(double));
	std::vector<double> n0_v(n0,n0+sizeof(n0)/sizeof(double));

	phi_BC * phir_ = new phir_power(n_v,d_v,t_v,l_v,1,19);
	phirlist.push_back(phir_);

	// return log(delta)+a0[0]+a0[1]*tau+a0[2]*log(tau)+a0[3]*log(1-exp(-n0[3]*tau))+a0[4]*log(1-exp(-n0[4]*tau))+a0[5]*log(1-exp(-n0[5]*tau))+a0[6]*log(1-exp(-n0[6]*tau));
	phi_BC * phi0_lead_ = new phi0_lead(a0[0],a0[1]);
	phi_BC * phi0_logtau_ = new phi0_logtau(a0[2]);
	phi_BC * phi0_Planck_Einstein_ = new phi0_Planck_Einstein(a0_v,n0_v,3,6);

	phi0list.push_back(phi0_lead_);
	phi0list.push_back(phi0_logtau_);
	phi0list.push_back(phi0_Planck_Einstein_);

	// Critical parameters
	crit.rho = 424;
	crit.p = 5782;
	crit.T = 351.255;
	crit.v = 1.0/crit.rho;

	// Other fluid parameters
	params.molemass = 52.024;
	params.Ttriple = 136.34;
	params.accentricfactor = 0.2769;
	params.R_u = 8.314471;

	// Limits of EOS
	limits.Tmin = params.Ttriple;
	limits.Tmax = 435.0;
	limits.pmax = 70000.0;
	limits.rhomax = 27.4734*params.molemass;
	
	EOSReference.assign("Tillner-Roth, R. and Yokozeki, A.,"
						" \"An international standard equation of state for difluoromethane (R-32)"
						" for temperatures from the triple point at 136.34 K to 435 K and pressures"
						" up to 70 MPa,\""
						" J. Phys. Chem. Ref. Data, 26(6):1273-1328, 1997.");
	TransportReference.assign("Surface Tension: R. Heide, \"The surface tension of HFC refrigerants and mixtures\", Int J. Refrig. Vol. 20, No. 7, pp. 496-503, 1997");

	name.assign("R32");
}
double R32Class::psat(double T)
{
	double theta,phi,p0=5781.16;

	phi=T/crit.T;
	theta=1-phi;

	return p0*exp((-7.44892*theta+1.6886*pow(theta,3.0/2.0)-1.908*pow(theta,5.0/2.0)-2.810*theta*theta*theta*theta*theta)/phi);
}

double R32Class::rhosatL(double T)
{
	double theta, phi;
	phi=T/crit.T;
	theta=1-phi;

	return 424.0+434.55*pow(theta,1.0/4.0)+1296.53*pow(theta,2.0/3.0)-777.49*theta+366.84*pow(theta,5.0/3.0);
}
double R32Class::rhosatV(double T)
{
	double theta, phi;
	phi=T/crit.T;
	theta=1-phi;

	return crit.rho*exp(-1.969*pow(theta,1.0/3.0)-2.0222*pow(theta,2.0/3.0)-6.7409*pow(theta,4.0/3.0)-27.479*pow(theta,11.0/3.0));
}
double R32Class::ECS_psi_viscosity(double rhor)
{
	return 0.7954+5.426580e-2*rhor;
}
double R32Class::ECS_f_int(double T)
{
	return 0.000436654+0.00000178134*T;
}
double R32Class::ECS_chi_conductivity(double rhor)
{
	return 1.29424-0.0924549*rhor;
}
void R32Class::ECSParams(double *e_k, double *sigma)
{
	*e_k = 289.65;
	*sigma = 0.4098;
}
double R32Class::viscosity_Trho(double T, double rho)
{
	// Use propane as the reference
	Fluid * ReferenceFluid = new R290Class();
	ReferenceFluid->post_load();
	// Calculate the ECS
	double mu = viscosity_ECS_Trho(T, rho, ReferenceFluid);
	// Delete the reference fluid instance
	delete ReferenceFluid;
	return mu;
}
double R32Class::conductivity_Trho(double T, double rho)
{
	// Use propane as the reference
	Fluid * ReferenceFluid = new R290Class();
	ReferenceFluid->post_load();
	// Calculate the ECS
	double cond = conductivity_ECS_Trho(T, rho, ReferenceFluid);
	// Delete the reference fluid instance
	delete ReferenceFluid;
	return cond;
}
double R32Class::surface_tension_T(double T)
{
	return 0.0705581*pow(1-T/reduce.T,1.234);
}