water.m

classdef water
    % WATER Physical properties of water dependent on temperature,
    %    from 1 to 99 degrees Celsius.
    %
    %    Water properties:
    %        .density [kg/m^3]
    %        .viscosity [m^2/s]
    %        .bulk_modulus [Pa]
    %
    %    Examples:
    %        % Compute water properties for a single temperature value
    %        water.density(27)
    %        water.viscosity(27)
    %        water.bulk_modulus(27)
    %        % Plot kinematic viscosity for multiple temperature values
    %        t = 1:50;
    %        plot(t,water.viscosity(t),'*')
    %
    %    Programmed by Ildeberto de los Santos Ruiz
    
    properties (Constant)
        nu = [1.7264;1.6706;1.6176;1.5671;1.5189;1.473;1.4292;1.3874;1.3475;1.3093;1.2728;1.2379;1.2044;1.1724;1.1417;1.1122;1.0839;1.0568;1.0307;1.0056;0.98155;0.95838;0.93607;0.9146;0.89391;0.87398;0.85476;0.83622;0.81834;0.80108;0.78441;0.76831;0.75275;0.73771;0.72317;0.70909;0.69548;0.68229;0.66953;0.65716;0.64518;0.63356;0.6223;0.61137;0.60077;0.59048;0.58049;0.57079;0.56137;0.55222;0.54332;0.53467;0.52626;0.51807;0.51011;0.50237;0.49482;0.48748;0.48033;0.47336;0.46658;0.45996;0.45351;0.44722;0.44109;0.43511;0.42928;0.42358;0.41802;0.4126;0.4073;0.40212;0.39707;0.39213;0.3873;0.38258;0.37797;0.37346;0.36905;0.36474;0.36052;0.35639;0.35235;0.3484;0.34453;0.34074;0.33703;0.33339;0.32983;0.32634;0.32292;0.31957;0.31629;0.31307;0.30992;0.30682;0.30379;0.30081;0.29789]*1e-6;
        rho = [1;0.99995;0.9999;0.99983;0.99976;0.99968;0.9996;0.9995;0.9994;0.99929;0.99918;0.99906;0.99893;0.99879;0.99864;0.99849;0.99833;0.99816;0.99798;0.9978;0.99761;0.99741;0.99721;0.99699;0.99677;0.99654;0.99631;0.99606;0.99581;0.99556;0.99529;0.99502;0.99474;0.99445;0.99415;0.99385;0.99354;0.99322;0.9929;0.99256;0.99222;0.99187;0.99152;0.99115;0.99078;0.99041;0.99002;0.98963;0.98923;0.98882;0.9884;0.98798;0.98755;0.98711;0.98667;0.98621;0.98575;0.98529;0.98481;0.98433;0.98384;0.98334;0.98283;0.98232;0.9818;0.98127;0.98074;0.98019;0.97964;0.97909;0.97852;0.97795;0.97737;0.97678;0.97619;0.97558;0.97497;0.97435;0.97373;0.9731;0.97246;0.97181;0.97115;0.97049;0.96982;0.96914;0.96846;0.96776;0.96706;0.96636;0.96564;0.96492;0.96419;0.96345;0.96271;0.96195;0.96119;0.96043;0.95965]*1e3;
        K = [2.0424;2.051;2.0595;2.0678;2.0759;2.0839;2.0917;2.0993;2.1068;2.1141;2.1212;2.1281;2.1349;2.1416;2.148;2.1544;2.1605;2.1665;2.1724;2.178;2.1836;2.189;2.1942;2.1993;2.2042;2.209;2.2136;2.2181;2.2224;2.2266;2.2306;2.2345;2.2383;2.2419;2.2454;2.2487;2.2519;2.255;2.2579;2.2607;2.2633;2.2659;2.2683;2.2705;2.2726;2.2747;2.2765;2.2783;2.2799;2.2814;2.2828;2.284;2.2851;2.2862;2.287;2.2878;2.2885;2.289;2.2894;2.2897;2.2899;2.29;2.29;2.2898;2.2896;2.2892;2.2888;2.2882;2.2875;2.2867;2.2858;2.2849;2.2838;2.2826;2.2813;2.2799;2.2784;2.2769;2.2752;2.2734;2.2716;2.2696;2.2676;2.2654;2.2632;2.2609;2.2585;2.256;2.2535;2.2508;2.2481;2.2453;2.2424;2.2394;2.2364;2.2332;2.23;2.2268;2.2234]*1e9;
    end
    methods (Static)
        function nu = viscosity(t)
            % Kinematic Viscosity, nu (m^2/s)
            nu = interp1(1:99,water.nu,t,'spline');
        end
        function rho = density(t)
            % Density, rho (kg/m^3)
            rho = interp1(1:99,water.rho,t,'spline');
        end
        function K = bulk_modulus(t)
            % Bulk Modulus, K (Pa)
            K = interp1(1:99,water.K,t,'spline');
        end
    end
end