PROPERTIES of (MEG) Mono Ethylene Glycol Mixture in Water

PROPERTIES of (MEG) Mono Ethylene Glycol Mixture in Water (Weight/Weight)

PROPERTIES of (MEG) Mono Ethylene Glycol Mixture in Water (Weight/Weight)

P = my_get_MEG_mix_props(tempC,pressBar,mixwwPerc)

tempC= temperature in degree C

pressBar= pressure in bar (unused, you can set it to 1)

mixwwPerc= mixture percent (Weight MEG / Weight Water) * 100

Look up tables obtained from:

INCOMPRESSIBLE FLUID

Bell, Ian H. and Wronski, Jorrit and Quoilin, Sylvain and Lemort, Vincent. Pure and Pseudo-pure Fluid Thermophysical Property Evaluation and the Open-Source Thermophysical Property Library CoolProp. 2014, Industrial & Engineering Chemistry Research. volume = {53}, number = {6}, pages = {2498--2508}, year = {2014}, doi = {10.1021/ie4033999},

AND

Melinder, Åke. (2007). Thermophysical properties of aqueous solutions used as secondary working fluids. Doctoral Thesis, Royal Institute of Technology, Stockholm.

Example usage:

clear all, close all, clc

itempC= linspace(-100,100, 21);

jmixwwPerc= linspace(0,100, 21);

for i= 1:numel(itempC)

for j= 1:numel(jmixwwPerc)

p= my_get_MEG_mix_props( itempC(i), 1, jmixwwPerc(j));

vv_perc(i,j)= p.vv_perc;

rho_kg_m3(i,j)= p.rho_kg_m3;

H_J_kg(i,j)= p.H_J_kg;

cp_J_kg_K(i,j)= p.cp_J_kg_K;

k_W_m_k(i,j)= p.k_W_m_k;

mu_Pa_s(i,j)= p.mu_Pa_s;

Prandtl(i,j)= p.Prandtl;

TFreeze_C(i,j)= p.TFreeze_C;

end

end

[X,Y]= meshgrid(jmixwwPerc,itempC);

subplot(2,3,1); plot3(X,Y,rho_kg_m3,'ko'), grid on

subplot(2,3,2); plot3(X,Y,H_J_kg,'ko'), grid on

subplot(2,3,3); plot3(X,Y,cp_J_kg_K,'ko'), grid on

subplot(2,3,4); plot3(X,Y,k_W_m_k,'ko'), grid on

subplot(2,3,5); plot3(X,Y,mu_Pa_s,'ko'), grid on

subplot(2,3,6); plot3(X,Y,Prandtl,'ko'), grid on