diff --git a/Roccaraso/postprocessing/RPS/cleanData/ROC-01/engineFlightData.mat b/Roccaraso/postprocessing/RPS/cleanData/ROC-01/engineFlightData.mat deleted file mode 100644 index 3d993ba0656809dfec3908594186567bb60b3361..0000000000000000000000000000000000000000 Binary files a/Roccaraso/postprocessing/RPS/cleanData/ROC-01/engineFlightData.mat and /dev/null differ diff --git a/Roccaraso/postprocessing/RPS/cleanData/ROC-02/engineFlightData.mat b/Roccaraso/postprocessing/RPS/cleanData/ROC-02/engineFlightData.mat deleted file mode 100644 index ada981048b03e3b1971113f3d0ba9420636cc36a..0000000000000000000000000000000000000000 Binary files a/Roccaraso/postprocessing/RPS/cleanData/ROC-02/engineFlightData.mat and /dev/null differ diff --git a/Roccaraso/postprocessing/RPS/cleanData/ROC1/engineFlightData.mat b/Roccaraso/postprocessing/RPS/cleanData/ROC1/engineFlightData.mat new file mode 100644 index 0000000000000000000000000000000000000000..34ff92c1d33966780accb5848141b50f28b76a72 Binary files /dev/null and b/Roccaraso/postprocessing/RPS/cleanData/ROC1/engineFlightData.mat differ diff --git a/Roccaraso/postprocessing/RPS/cleanData/ROC2/engineFlightData.mat b/Roccaraso/postprocessing/RPS/cleanData/ROC2/engineFlightData.mat new file mode 100644 index 0000000000000000000000000000000000000000..cfd6deb231917c80950f74c2dbfb8ad51dc084f0 Binary files /dev/null and b/Roccaraso/postprocessing/RPS/cleanData/ROC2/engineFlightData.mat differ diff --git a/Roccaraso/postprocessing/RPS/configDataAnalyzer.m b/Roccaraso/postprocessing/RPS/configDataAnalyzer.m new file mode 100644 index 0000000000000000000000000000000000000000..5e45f0fb52fdfdf25262ec17e7b29681e7982524 --- /dev/null +++ b/Roccaraso/postprocessing/RPS/configDataAnalyzer.m @@ -0,0 +1,32 @@ +%% FLIGHT and ACQUISITION SETTINGS +settings.launchSite = 'Roccaraso'; % [-] Launch site name +settings.testName = 'ROC2'; % [-] Test name + +settings.frPT = 1000; % [Hz] PTs sampling rate +settings.frSrv = 10; % [Hz] SRVs sampling rate +settings.fileCSVpath = strcat('../../../', settings.launchSite, '/data/MOTORE/'); + +% Load mission parameters +run(strcat('engineData/missions/', settings.testName)); + +%% ROCKET SELECTION +% Choose the motor you want to fly +settings.engineName = 'Furia'; + +% Retrieve engine data +run(strcat('engineData/engines/', settings.engineName, '/', settings.engineName)); + +%% ODE SETTINGS +% !! Do not change unless you know what you are doing !! +settings.odeEngineSettings = odeset('RelTol', 1e-5, 'AbsTol', 1e-5, 'Events', @eventBurnOut); +settings.fsolveSettings = optimoptions('fsolve', 'Display', 'off'); +settings.odeRegressionRate = 'Marxman'; % choose between Marxman and Whitmore + +%% BISECTION SETTINGS +settings.bisection.maxIter = 20; +settings.bisection.xTol = 1e-3; +settings.bisection.absTol = 1e-3; + +%% PLOT SETTINGS +settings.plots = true; % [-] True if you want to see the plots + diff --git a/Roccaraso/postprocessing/RPS/dataAnalyzer.m b/Roccaraso/postprocessing/RPS/dataAnalyzer.m index bfce82bcd8c436db424a50ff6ccb8c251f07136e..b3f9c00e32f4fa8641288d81ef355cb38d87864e 100644 --- a/Roccaraso/postprocessing/RPS/dataAnalyzer.m +++ b/Roccaraso/postprocessing/RPS/dataAnalyzer.m @@ -1,6 +1,8 @@ -clear all; close all; clc +clear all; +close all; +clc -%% PATHS & PLOTS +%% PATHS & PLOTS PROPERTIES addpath(genpath('./')) % Settings the default plot properties @@ -15,49 +17,35 @@ set(0, 'defaultAxesColor', 'none'); set(groot,'DefaultAxesFontSize', 13); %% DATA & SETTINGS -settings.engine.cT = 1.3606; % [-] Mean thrust coefficient - corrected -settings.engine.lam = 0.9853; % [-] Nozzle 2D efficiency -settings.engine.rt = 14.9e-3; % [m] Throat radius -settings.engine.Isp = 173.3; % [s] Specific impulse - -settings.frPT = 1000; % [Hz] PTs sampling rate -settings.frSrv = 10; % [Hz] SRVs sampling rate -settings.g0 = 9.806; % [m/s] Acceleration of Gravity, constant -settings.launchSite = 'Roccaraso'; % [-] Launch site name -settings.testName = 'ROC-02'; % [-] Test name - -settings.fileCSVpath = strcat('../../../', settings.launchSite, '/data/MOTORE/'); +% Call settings function +configDataAnalyzer +% Load previously saved data if file exists (saved in dataLoader.m) if ~isfile(strcat('cleanData/', settings.testName, '/engineFlightData.mat')) flightData = dataLoader(settings); else load(strcat('cleanData/', settings.testName, '/engineFlightData.mat')); end - -%% CODE -f1 = figure; -yyaxis left -plot(flightData.timePTop, flightData.PTankTop) -hold on; grid minor -ylabel('Pressure [bar]') -plot(flightData.timePBot, flightData.PTankBot) -yyaxis right -stairs(flightData.timeSrvMain, flightData.srvMain) -xlabel('Time [s]'); ylabel('SRV Opening [-]'); title('Tank pressure values') -legend('Top tank', 'Bottom tank') - -f2 = figure; -plot(flightData.timePEng, flightData.PEng) -grid minor; xlabel('Time [s]'); ylabel('Pressure [bar]'); title('Engine pressure value') - -f3 = figure; -plot(flightData.timeThrust, flightData.Thrust) -grid minor; xlabel('Time [s]'); ylabel('Thrust [kgf]'); title('Expected thrust') - -f4 = figure; -plot(flightData.timeMp, flightData.mp) -grid minor; xlabel('Time [s]'); ylabel('Mass Flow [kg/s]'); title('Expected propellant mass flow rate') +%% DATA PROCESSING AND ANALYSIS +engine.acquisition = flightData; + +engine = tankEvacuation(engine, settings); + +%% PLOTTING +if settings.plots + f1 = figure; + yyaxis left + plot(flightData.timePTop, flightData.PTankTop) + hold on; grid minor + ylabel('Pressure [bar]') + plot(flightData.timePBot, flightData.PTankBot) + plot(flightData.timePEng, flightData.PEng) + yyaxis right + stairs(flightData.timeSrvMain, flightData.srvMain) + xlabel('Time [s]'); ylabel('SRV Opening [-]'); title('Tank pressure values') + legend('Top tank', 'Bottom tank', 'Combustion chamber') +end diff --git a/Roccaraso/postprocessing/RPS/engineData/engines/CEA_db.mat b/Roccaraso/postprocessing/RPS/engineData/engines/CEA_db.mat new file mode 100644 index 0000000000000000000000000000000000000000..2712efdac18d9d0c55b4166d865a2d70d69286e8 Binary files /dev/null and b/Roccaraso/postprocessing/RPS/engineData/engines/CEA_db.mat differ diff --git a/Roccaraso/postprocessing/RPS/engineData/engines/Furia/Furia.m b/Roccaraso/postprocessing/RPS/engineData/engines/Furia/Furia.m new file mode 100644 index 0000000000000000000000000000000000000000..09f094e2b292763bb2615892ff7fae2645ea4120 --- /dev/null +++ b/Roccaraso/postprocessing/RPS/engineData/engines/Furia/Furia.m @@ -0,0 +1,63 @@ +%{ + +ROCKET DATA FILE + +Project name: HRE FURIA +Year of launch: 2023 +Location of Launch: TBD +Date of launch: TBD + +Last update: 02/12/2022 + +Copyright © 2022, Skyward Experimental Rocketry, PRP department +All rights reserved + +SPDX-License-Identifier: GPL-3.0-or-later + +%} + +engine.name = 'Furia'; + +%% ENGINE PARAMETERS +engine.oxName = 'N2O'; % [-] Oxydizer name +engine.fuName = 'ABS'; % [-] Fuel name + +engine.grainGeometry = 'Cylinder'; % [-] Grain geometry + +%%% LOAD DATA +addpath(genpath('../../../oxydizers')); +addpath(genpath('../../../fuels')); +addpath(genpath('../../')); +engine.data.ox = eval(engine.oxName); +engine.data.fu = eval(engine.fuName); + +if exist('CEA_db.mat', 'file') + engine.data.CEA = load('CEA_db.mat'); +else + engine.data.CEA = []; +end + +%% ENGINE GEOMETRY + +%%% FEED LINE +engine.data.d_pipe = 1.25e-2; % [m] Diameter of the feedline pipes +engine.data.dInj = 1.5e-3; % [m] Diameter of injection holes +engine.data.Aox = engine.data.Ninj* ... + pi/4*engine.data.dInj.^2; % [m^2] Injection area + +%%% EFFICIENCIES +engine.data.etaCT = 0.96; +engine.data.etaCs = 0.85; + +%%% COMBUSTION CHAMBER +engine.data.rt = 0.0149; % [m] Throat radius +engine.data.epsilon = 3.4777; % [-] Expansion ratio +engine.data.lambda = 0.985; % [-] Loss due to nozzle divergence +engine.data.Rc_CM = 0.5; +engine.data.R_pre = 1.15; % [-] Ratio between pre-chamber length and port radius +engine.data.R_post = 2.35; % [-] ratio between post-chamber length and port radius + +%%% NOZZLE +engine.data.zopt = 0; % [m] Altitude of optimal expansion +engine.data.RAO_Lratio = 0.8; % [-] Rao length ratio + diff --git a/Roccaraso/postprocessing/RPS/engineData/fuels/ABS.m b/Roccaraso/postprocessing/RPS/engineData/fuels/ABS.m new file mode 100644 index 0000000000000000000000000000000000000000..592281101e5331f9e1e711675d75aee7b09f808b --- /dev/null +++ b/Roccaraso/postprocessing/RPS/engineData/fuels/ABS.m @@ -0,0 +1,49 @@ +function ABS = ABS() +% ------------------------------------------------------------------------- +% This function creates a dataset for curves of saturation for N2O using data +% from NIST database +% +% PROTOTYPE +% [N2O] = N2Osatcurve() +% +% INPUT +% +% OUTPUT +% ABS [struct] temperature [k], pressure [bar], densities [kg/m^3], +% specific entropies [J/K/Kg] and enthalpies [J/Kg] +% of N2O at saturation state +% +% CONTRIBUTORS +% Alberto Boffi +% +% VERSION +% 4-11-2022: v01.0 +%-------------------------------------------------------------------------% + +%% CEA PARAMETERS +ABS.CEA.name = 'ABS'; % [-] Fuel CEA database name +ABS.CEA.T = 293; % [K] Fuel temperature + +% Set the following to true if you want to input oxydizer chemical +% composition +ABS.CEA.elementsFlag = true; +% If the previous is true, set the following: +ABS.CEA.elements = {'C', 'H', 'N'}; +ABS.CEA.elConctr = [3.85 4.85 0.43]; +ABS.CEA.h = 62.63; + + +%% PERFORMANCE PARAMETERS +% Literature values +ABS.aLit = 1.1787e-5; % [m^(2n+1)/s^(n-1)/kg^n] Parameter a for Vieille's equation +ABS.nLit = 0.7990; % [-] Parameter n for Vieille's equation + +% Gox fitted Chimaera paramters calculated with SFT5 & SFT6 +ABS.a = 1.188 * 1e-4; % [m^(2n+1)/s^(n-1)/kg^n] Parameter a for Vieille's equation +ABS.n = 0.395; % [-] Parameter n for Vieille's equation + +ABS.rhofNom = 1007; % [kg/m^3] Nominal grain density +ABS.rhof = 987.6; % [kg/m^3] Measured grain density +ABS.hv = 3.07*1e6; % [J/kg] enthalpy of gassification +ABS.Tf = 713; % [K] gassification temperature +end \ No newline at end of file diff --git a/Roccaraso/postprocessing/RPS/engineData/missions/ROC1.m b/Roccaraso/postprocessing/RPS/engineData/missions/ROC1.m new file mode 100644 index 0000000000000000000000000000000000000000..9a39eb3c70e980ea0975bbd0e1478b38845ddf04 --- /dev/null +++ b/Roccaraso/postprocessing/RPS/engineData/missions/ROC1.m @@ -0,0 +1,18 @@ +% Environment +engine.data.T0 = 298.15; + +% Injector choice +engine.data.Ninj = 22; +engine.data.cd_inj = 0.7; + +% Residual web thickness +engine.data.resWeb = 0.3 + 1.4; % [mm] Residual fuel web thickness -- sum of SF1 and SF2 + +%%% FUEL +engine.data.grainGeometry = 'Cylinder'; % Choose between Cylinder and Star +engine.data.Lgrain = 0.510; % [m] Grain length +engine.data.ri = 0.0256; % [m] Grain core radius +engine.data.rf = 0.0369; % [m] Grain outer radius +engine.data.r0 = engine.data.ri; % [m] Initial regression coordinate - to be set equal to ri for cylinder +engine.data.webThickness = 0.0; % [m] Web thickness +engine.data.spikeNumber = 0; % [-] Number of spikes \ No newline at end of file diff --git a/Roccaraso/postprocessing/RPS/engineData/missions/ROC2.m b/Roccaraso/postprocessing/RPS/engineData/missions/ROC2.m new file mode 100644 index 0000000000000000000000000000000000000000..9a39eb3c70e980ea0975bbd0e1478b38845ddf04 --- /dev/null +++ b/Roccaraso/postprocessing/RPS/engineData/missions/ROC2.m @@ -0,0 +1,18 @@ +% Environment +engine.data.T0 = 298.15; + +% Injector choice +engine.data.Ninj = 22; +engine.data.cd_inj = 0.7; + +% Residual web thickness +engine.data.resWeb = 0.3 + 1.4; % [mm] Residual fuel web thickness -- sum of SF1 and SF2 + +%%% FUEL +engine.data.grainGeometry = 'Cylinder'; % Choose between Cylinder and Star +engine.data.Lgrain = 0.510; % [m] Grain length +engine.data.ri = 0.0256; % [m] Grain core radius +engine.data.rf = 0.0369; % [m] Grain outer radius +engine.data.r0 = engine.data.ri; % [m] Initial regression coordinate - to be set equal to ri for cylinder +engine.data.webThickness = 0.0; % [m] Web thickness +engine.data.spikeNumber = 0; % [-] Number of spikes \ No newline at end of file diff --git a/Roccaraso/postprocessing/RPS/engineData/oxydizers/N2O.m b/Roccaraso/postprocessing/RPS/engineData/oxydizers/N2O.m new file mode 100644 index 0000000000000000000000000000000000000000..5f1c6cca18c67631bf9c1b073d5a44e2dc70d0e1 --- /dev/null +++ b/Roccaraso/postprocessing/RPS/engineData/oxydizers/N2O.m @@ -0,0 +1,56 @@ +function [N2O] = N2O +% ------------------------------------------------------------------------- +% This function creates a dataset for curves of saturation for N2O using data +% from NIST database +% +% PROTOTYPE +% [N2O] = N2Osatcurve() +% +% INPUT +% +% OUTPUT +% N2O [struct] temperature [k], pressure [bar], densities [kg/m^3], +% specific entropies [J/K/Kg] and enthalpies [J/Kg] +% of N2O at saturation state +% +% CONTRIBUTORS +% Alberto Boffi +% +% VERSION +% 4-11-2022: v01.0 +%-------------------------------------------------------------------------% + +%% CEA PARAMETERS +N2O.CEA.name = 'N2O'; % [-] Oxydizer CEA database name +N2O.CEA.T = 293; % [K] Oxydizer temperature + +% Set the following to true if you want to input oxydizer chemical +% composition +N2O.CEA.elementsFlag = false; +% If the previous is true, set the following: +N2O.CEA.elements = {''}; +N2O.CEA.elConctr = []; +N2O.CEA.h = 0; + + +%% SATURATION CURVES +N2O.T = [182.33 184.43 186.33 188.08 189.7 191.22 192.64 193.98 195.25 196.46 197.61 198.71 199.77 200.78 201.76 202.7 203.61 204.49 205.35 206.17 206.98 207.76 208.52 209.27 209.99 210.7 211.39 212.06 212.72 213.37 214 214.62 215.23 215.82 216.41 216.98 217.55 218.1 218.65 219.18 219.71 220.23 220.74 221.25 221.74 222.23 222.71 223.19 223.66 224.12 224.58 225.03 225.47 225.91 226.35 226.78 227.2 227.62 228.03 228.44 228.85 229.25 229.64 230.03 230.42 230.81 231.19 231.56 231.93 232.3 232.67 233.03 233.39 233.74 234.09 234.44 234.79 235.13 235.47 235.81 236.14 236.47 236.8 237.12 237.45 237.77 238.08 238.4 238.71 239.02 239.33 239.63 239.94 240.24 240.54 240.83 241.13 241.42 241.71 242 242.28 242.57 242.85 243.13 243.41 243.69 243.96 244.24 244.51 244.78 245.05 245.31 245.58 245.84 246.1 246.36 246.62 246.88 247.13 247.39 247.64 247.89 248.14 248.39 248.64 248.88 249.13 249.37 249.61 249.85 250.09 250.33 250.56 250.8 251.03 251.26 251.5 251.73 251.96 252.18 252.41 252.64 252.86 253.08 253.31 253.53 253.75 253.97 254.19 254.4 254.62 254.84 255.05 255.26 255.48 255.69 255.9 256.11 256.32 256.52 256.73 256.94 257.14 257.35 257.55 257.75 257.95 258.15 258.35 258.55 258.75 258.95 259.14 259.34 259.54 259.73 259.92 260.12 260.31 260.5 260.69 260.88 261.07 261.26 261.44 261.63 261.82 262 262.19 262.37 262.55 262.74 262.92 263.1 263.28 263.46 263.64 263.82 264 264.17 264.35 264.53 264.7 264.88 265.05 265.23 265.4 265.57 265.74 265.92 266.09 266.26 266.43 266.6 266.76 266.93 267.1 267.27 267.43 267.6 267.76 267.93 268.09 268.26 268.42 268.58 268.75 268.91 269.07 269.23 269.39 269.55 269.71 269.87 270.03 270.18 270.34 270.5 270.65 270.81 270.97 271.12 271.27 271.43 271.58 271.74 271.89 272.04 272.19 272.34 272.5 272.65 272.8 272.95 273.1 273.24 273.39 273.54 273.69 273.84 273.98 274.13 274.28 274.42 274.57 274.71 274.86 275 275.14 275.29 275.43 275.57 275.71 275.86 276 276.14 276.28 276.42 276.56 276.7 276.84 276.98 277.12 277.25 277.39 277.53 277.67 277.8 277.94 278.08 278.21 278.35 278.48 278.62 278.75 278.89 279.02 279.15 279.29 279.42 279.55 279.69 279.82 279.95 280.08 280.21 280.34 280.47 280.6 280.73 280.86 280.99 281.12 281.25 281.38 281.5 281.63 281.76 281.89 282.01 282.14 282.27 282.39 282.52 282.64 282.77 282.89 283.02 283.14 283.27 283.39 283.51 283.64 283.76 283.88 284 284.13 284.25 284.37 284.49 284.61 284.73 284.85 284.97 285.09 285.21 285.33 285.45 285.57 285.69 285.81 285.93 286.04 286.16 286.28 286.4 286.51 286.63 286.75 286.86 286.98 287.1 287.21 287.33 287.44 287.56 287.67 287.79 287.9 288.01 288.13 288.24 288.35 288.47 288.58 288.69 288.81 288.92 289.03 289.14 289.25 289.36 289.48 289.59 289.7 289.81 289.92 290.03 290.14 290.25 290.36 290.47 290.57 290.68 290.79 290.9 291.01 291.12 291.22 291.33 291.44 291.55 291.65 291.76 291.87 291.97 292.08 292.18 292.29 292.4 292.5 292.61 292.71 292.82 292.92 293.03 293.13 293.23 293.34 293.44 293.54 293.65 293.75 293.85 293.96 294.06 294.16 294.26 294.36 294.47 294.57 294.67 294.77 294.87 294.97 295.07 295.17 295.27 295.37 295.47 295.57 295.67 295.77 295.87 295.97 296.07 296.17 296.27 296.37 296.46 296.56 296.66 296.76 296.86 296.95 297.05 297.15 297.24 297.34 297.44 297.53 297.63 297.73 297.82 297.92 298.01 298.11 298.2 298.3 298.39 298.49 298.58 298.68 298.77 298.87 298.96 299.05 299.15 299.24 299.33 299.43 299.52 299.61 299.71 299.8 299.89 299.98 300.07 300.17 300.26 300.35 300.44 300.53 300.62 300.71 300.81 300.9 300.99 301.08 301.17 301.26 301.35 301.44 301.53 301.62 301.71 301.79 301.88 301.97 302.06 302.15 302.24 302.33 302.41 302.5 302.59 302.68 302.77 302.85 302.94 303.03 303.11 303.2 303.29 303.37 303.46 303.55 303.63 303.72 303.8 303.89 303.98 304.06 304.15 304.23 304.32 304.4 304.49 304.57 304.66 304.74 304.82 304.91 304.99 305.07 305.16 305.24 305.33 305.41 305.49 305.57 305.66 305.74 305.82 305.9 305.99 306.07 306.15 306.23 306.31 306.39 306.48 306.56 306.64 306.72 306.8 306.88 306.96 307.04 307.12 307.2 307.28 307.36 307.44 307.52 307.6 307.68 307.76 307.83 307.91 307.99 308.07 308.15 308.22 308.3 308.38 308.46 308.53 308.61 308.69 308.77 308.84 308.92 308.99 309.07 309.15 309.22 309.3 309.37 309.45 309.52]; +N2O.P = [0.8784 0.99769 1.117 1.2363 1.3555 1.4748 1.5941 1.7134 1.8327 1.952 2.0713 2.1905 2.3098 2.4291 2.5484 2.6677 2.787 2.9063 3.0255 3.1448 3.2641 3.3834 3.5027 3.622 3.7413 3.8605 3.9798 4.0991 4.2184 4.3377 4.457 4.5763 4.6956 4.8148 4.9341 5.0534 5.1727 5.292 5.4113 5.5306 5.6498 5.7691 5.8884 6.0077 6.127 6.2463 6.3656 6.4848 6.6041 6.7234 6.8427 6.962 7.0813 7.2006 7.3198 7.4391 7.5584 7.6777 7.797 7.9163 8.0356 8.1548 8.2741 8.3934 8.5127 8.632 8.7513 8.8706 8.9898 9.1091 9.2284 9.3477 9.467 9.5863 9.7056 9.8248 9.9441 10.063 10.183 10.302 10.421 10.541 10.66 10.779 10.898 11.018 11.137 11.256 11.376 11.495 11.614 11.733 11.853 11.972 12.091 12.211 12.33 12.449 12.568 12.688 12.807 12.926 13.046 13.165 13.284 13.403 13.523 13.642 13.761 13.881 14 14.119 14.238 14.358 14.477 14.596 14.716 14.835 14.954 15.073 15.193 15.312 15.431 15.551 15.67 15.789 15.908 16.028 16.147 16.266 16.386 16.505 16.624 16.743 16.863 16.982 17.101 17.221 17.34 17.459 17.578 17.698 17.817 17.936 18.056 18.175 18.294 18.413 18.533 18.652 18.771 18.891 19.01 19.129 19.248 19.368 19.487 19.606 19.726 19.845 19.964 20.083 20.203 20.322 20.441 20.561 20.68 20.799 20.918 21.038 21.157 21.276 21.396 21.515 21.634 21.753 21.873 21.992 22.111 22.231 22.35 22.469 22.588 22.708 22.827 22.946 23.066 23.185 23.304 23.423 23.543 23.662 23.781 23.901 24.02 24.139 24.258 24.378 24.497 24.616 24.736 24.855 24.974 25.093 25.213 25.332 25.451 25.571 25.69 25.809 25.928 26.048 26.167 26.286 26.406 26.525 26.644 26.763 26.883 27.002 27.121 27.241 27.36 27.479 27.598 27.718 27.837 27.956 28.076 28.195 28.314 28.433 28.553 28.672 28.791 28.911 29.03 29.149 29.268 29.388 29.507 29.626 29.746 29.865 29.984 30.103 30.223 30.342 30.461 30.581 30.7 30.819 30.938 31.058 31.177 31.296 31.416 31.535 31.654 31.773 31.893 32.012 32.131 32.251 32.37 32.489 32.608 32.728 32.847 32.966 33.086 33.205 33.324 33.443 33.563 33.682 33.801 33.921 34.04 34.159 34.278 34.398 34.517 34.636 34.756 34.875 34.994 35.113 35.233 35.352 35.471 35.591 35.71 35.829 35.948 36.068 36.187 36.306 36.426 36.545 36.664 36.783 36.903 37.022 37.141 37.261 37.38 37.499 37.618 37.738 37.857 37.976 38.096 38.215 38.334 38.453 38.573 38.692 38.811 38.931 39.05 39.169 39.288 39.408 39.527 39.646 39.766 39.885 40.004 40.123 40.243 40.362 40.481 40.601 40.72 40.839 40.958 41.078 41.197 41.316 41.436 41.555 41.674 41.793 41.913 42.032 42.151 42.271 42.39 42.509 42.629 42.748 42.867 42.986 43.106 43.225 43.344 43.464 43.583 43.702 43.821 43.941 44.06 44.179 44.299 44.418 44.537 44.656 44.776 44.895 45.014 45.134 45.253 45.372 45.491 45.611 45.73 45.849 45.969 46.088 46.207 46.326 46.446 46.565 46.684 46.804 46.923 47.042 47.161 47.281 47.4 47.519 47.639 47.758 47.877 47.996 48.116 48.235 48.354 48.474 48.593 48.712 48.831 48.951 49.07 49.189 49.309 49.428 49.547 49.666 49.786 49.905 50.024 50.144 50.263 50.382 50.501 50.621 50.74 50.859 50.979 51.098 51.217 51.336 51.456 51.575 51.694 51.814 51.933 52.052 52.171 52.291 52.41 52.529 52.649 52.768 52.887 53.006 53.126 53.245 53.364 53.484 53.603 53.722 53.841 53.961 54.08 54.199 54.319 54.438 54.557 54.676 54.796 54.915 55.034 55.154 55.273 55.392 55.511 55.631 55.75 55.869 55.989 56.108 56.227 56.346 56.466 56.585 56.704 56.824 56.943 57.062 57.181 57.301 57.42 57.539 57.659 57.778 57.897 58.016 58.136 58.255 58.374 58.494 58.613 58.732 58.851 58.971 59.09 59.209 59.329 59.448 59.567 59.686 59.806 59.925 60.044 60.164 60.283 60.402 60.521 60.641 60.76 60.879 60.999 61.118 61.237 61.356 61.476 61.595 61.714 61.834 61.953 62.072 62.191 62.311 62.43 62.549 62.669 62.788 62.907 63.026 63.146 63.265 63.384 63.504 63.623 63.742 63.861 63.981 64.1 64.219 64.339 64.458 64.577 64.696 64.816 64.935 65.054 65.174 65.293 65.412 65.531 65.651 65.77 65.889 66.009 66.128 66.247 66.366 66.486 66.605 66.724 66.844 66.963 67.082 67.201 67.321 67.44 67.559 67.679 67.798 67.917 68.036 68.156 68.275 68.394 68.514 68.633 68.752 68.871 68.991 69.11 69.229 69.349 69.468 69.587 69.706 69.826 69.945 70.064 70.184 70.303 70.422 70.541 70.661 70.78 70.899 71.019 71.138 71.257 71.376 71.496 71.615 71.734 71.854 71.973 72.092 72.211 72.331 72.447]; + +N2O.rhoL = [1237.3 1231.2 1225.6 1220.4 1215.6 1211 1206.8 1202.7 1198.9 1195.2 1191.7 1188.3 1185.1 1182 1178.9 1176 1173.2 1170.4 1167.7 1165.1 1162.6 1160.1 1157.7 1155.3 1153 1150.8 1148.5 1146.4 1144.2 1142.1 1140.1 1138.1 1136.1 1134.1 1132.2 1130.3 1128.4 1126.6 1124.8 1123 1121.2 1119.5 1117.8 1116.1 1114.4 1112.8 1111.1 1109.5 1107.9 1106.3 1104.8 1103.2 1101.7 1100.2 1098.7 1097.2 1095.7 1094.3 1092.9 1091.4 1090 1088.6 1087.2 1085.8 1084.5 1083.1 1081.8 1080.4 1079.1 1077.8 1076.5 1075.2 1073.9 1072.6 1071.4 1070.1 1068.9 1067.6 1066.4 1065.2 1063.9 1062.7 1061.5 1060.3 1059.2 1058 1056.8 1055.6 1054.5 1053.3 1052.2 1051 1049.9 1048.8 1047.7 1046.5 1045.4 1044.3 1043.2 1042.1 1041 1040 1038.9 1037.8 1036.7 1035.7 1034.6 1033.6 1032.5 1031.5 1030.5 1029.4 1028.4 1027.4 1026.3 1025.3 1024.3 1023.3 1022.3 1021.3 1020.3 1019.3 1018.3 1017.3 1016.4 1015.4 1014.4 1013.4 1012.5 1011.5 1010.5 1009.6 1008.6 1007.7 1006.7 1005.8 1004.8 1003.9 1003 1002 1001.1 1000.2 999.26 998.34 997.42 996.51 995.59 994.68 993.78 992.87 991.97 991.07 990.17 989.28 988.38 987.49 986.61 985.72 984.84 983.96 983.08 982.2 981.33 980.45 979.58 978.72 977.85 976.99 976.12 975.26 974.41 973.55 972.7 971.84 970.99 970.15 969.3 968.45 967.61 966.77 965.93 965.09 964.26 963.42 962.59 961.76 960.93 960.1 959.27 958.45 957.62 956.8 955.98 955.16 954.34 953.53 952.71 951.9 951.09 950.28 949.47 948.66 947.85 947.05 946.24 945.44 944.64 943.84 943.04 942.24 941.45 940.65 939.86 939.06 938.27 937.48 936.69 935.9 935.11 934.33 933.54 932.75 931.97 931.19 930.41 929.63 928.85 928.07 927.29 926.51 925.74 924.96 924.19 923.41 922.64 921.87 921.1 920.33 919.56 918.79 918.02 917.25 916.49 915.72 914.96 914.19 913.43 912.67 911.9 911.14 910.38 909.62 908.86 908.1 907.34 906.59 905.83 905.07 904.32 903.56 902.81 902.05 901.3 900.55 899.79 899.04 898.29 897.54 896.79 896.04 895.29 894.54 893.79 893.04 892.29 891.54 890.8 890.05 889.3 888.56 887.81 887.06 886.32 885.57 884.83 884.08 883.34 882.6 881.85 881.11 880.37 879.62 878.88 878.14 877.4 876.65 875.91 875.17 874.43 873.69 872.95 872.21 871.47 870.72 869.98 869.24 868.5 867.76 867.02 866.28 865.54 864.8 864.06 863.32 862.58 861.84 861.1 860.36 859.62 858.88 858.14 857.4 856.66 855.92 855.18 854.44 853.7 852.96 852.22 851.47 850.73 849.99 849.25 848.51 847.77 847.03 846.28 845.54 844.8 844.06 843.31 842.57 841.83 841.08 840.34 839.59 838.85 838.1 837.36 836.61 835.87 835.12 834.37 833.63 832.88 832.13 831.38 830.63 829.88 829.13 828.38 827.63 826.88 826.13 825.38 824.63 823.87 823.12 822.37 821.61 820.86 820.1 819.34 818.59 817.83 817.07 816.31 815.55 814.79 814.03 813.27 812.51 811.74 810.98 810.21 809.45 808.68 807.92 807.15 806.38 805.61 804.84 804.07 803.3 802.52 801.75 800.98 800.2 799.42 798.65 797.87 797.09 796.31 795.53 794.74 793.96 793.17 792.39 791.6 790.81 790.02 789.23 788.44 787.65 786.85 786.06 785.26 784.46 783.67 782.87 782.06 781.26 780.46 779.65 778.84 778.03 777.22 776.41 775.6 774.78 773.97 773.15 772.33 771.51 770.69 769.86 769.04 768.21 767.38 766.55 765.71 764.88 764.04 763.2 762.36 761.52 760.67 759.83 758.98 758.13 757.28 756.42 755.57 754.71 753.84 752.98 752.12 751.25 750.38 749.5 748.63 747.75 746.87 745.99 745.1 744.21 743.32 742.43 741.54 740.64 739.73 738.83 737.92 737.01 736.1 735.18 734.26 733.34 732.41 731.49 730.55 729.62 728.68 727.74 726.79 725.84 724.89 723.93 722.97 722 721.03 720.06 719.08 718.1 717.12 716.13 715.13 714.13 713.13 712.12 711.11 710.09 709.07 708.04 707.01 705.98 704.93 703.88 702.83 701.77 700.71 699.64 698.56 697.48 696.39 695.29 694.19 693.08 691.96 690.84 689.71 688.58 687.43 686.28 685.12 683.95 682.78 681.59 680.4 679.2 677.99 676.77 675.54 674.3 673.05 671.79 670.52 669.24 667.94 666.64 665.32 663.99 662.65 661.3 659.93 658.55 657.16 655.74 654.32 652.88 651.42 649.95 648.45 646.94 645.41 643.86 642.29 640.7 639.09 637.45 635.79 634.11 632.4 630.66 628.89 627.09 625.26 623.4 621.5 619.56 617.58 615.57 613.5 611.39 609.23 607.02 604.75 602.42 600.02 597.55 595 592.36 589.63 586.8 583.86 580.79 577.58 574.21 570.66 566.9 562.9 558.62 553.99 548.94 543.37 537.1 529.9 521.31 510.48 495.1 452.01]; +N2O.rhoV = [2.6116 2.9395 3.2647 3.5875 3.9084 4.2275 4.5449 4.8609 5.1756 5.4891 5.8016 6.113 6.4235 6.7332 7.0422 7.3504 7.6579 7.9649 8.2713 8.5771 8.8825 9.1874 9.4919 9.7961 10.1 10.403 10.706 11.009 11.312 11.614 11.916 12.218 12.52 12.821 13.123 13.424 13.725 14.026 14.327 14.628 14.929 15.229 15.53 15.831 16.131 16.432 16.732 17.033 17.333 17.634 17.934 18.235 18.535 18.836 19.137 19.437 19.738 20.039 20.34 20.641 20.942 21.243 21.545 21.846 22.148 22.449 22.751 23.053 23.355 23.657 23.959 24.262 24.564 24.867 25.17 25.473 25.776 26.08 26.383 26.687 26.991 27.295 27.599 27.903 28.208 28.513 28.818 29.123 29.429 29.734 30.04 30.346 30.653 30.959 31.266 31.573 31.88 32.188 32.495 32.803 33.112 33.42 33.729 34.038 34.347 34.656 34.966 35.276 35.586 35.897 36.208 36.519 36.83 37.142 37.454 37.766 38.078 38.391 38.704 39.018 39.331 39.645 39.96 40.274 40.589 40.904 41.22 41.536 41.852 42.168 42.485 42.802 43.119 43.437 43.755 44.074 44.393 44.712 45.031 45.351 45.671 45.992 46.312 46.634 46.955 47.277 47.599 47.922 48.245 48.568 48.892 49.216 49.541 49.865 50.191 50.516 50.842 51.169 51.495 51.823 52.15 52.478 52.806 53.135 53.464 53.794 54.124 54.454 54.785 55.116 55.447 55.779 56.112 56.445 56.778 57.112 57.446 57.78 58.115 58.451 58.787 59.123 59.46 59.797 60.134 60.472 60.811 61.15 61.489 61.829 62.17 62.51 62.852 63.194 63.536 63.878 64.222 64.565 64.909 65.254 65.599 65.945 66.291 66.638 66.985 67.332 67.68 68.029 68.378 68.728 69.078 69.428 69.779 70.131 70.483 70.836 71.189 71.543 71.897 72.252 72.608 72.964 73.32 73.677 74.035 74.393 74.752 75.111 75.471 75.831 76.192 76.553 76.915 77.278 77.641 78.005 78.37 78.735 79.1 79.466 79.833 80.201 80.569 80.937 81.306 81.676 82.047 82.418 82.789 83.162 83.535 83.908 84.282 84.657 85.033 85.409 85.786 86.163 86.541 86.92 87.299 87.68 88.06 88.442 88.824 89.207 89.59 89.974 90.359 90.745 91.131 91.518 91.906 92.294 92.683 93.073 93.464 93.855 94.247 94.64 95.033 95.427 95.822 96.218 96.615 97.012 97.41 97.809 98.208 98.609 99.01 99.412 99.814 100.22 100.62 101.03 101.43 101.84 102.25 102.66 103.06 103.47 103.89 104.3 104.71 105.12 105.54 105.95 106.37 106.79 107.2 107.62 108.04 108.46 108.88 109.31 109.73 110.15 110.58 111 111.43 111.86 112.29 112.72 113.15 113.58 114.01 114.44 114.88 115.31 115.75 116.19 116.63 117.07 117.51 117.95 118.39 118.84 119.28 119.73 120.17 120.62 121.07 121.52 121.97 122.42 122.88 123.33 123.79 124.24 124.7 125.16 125.62 126.08 126.54 127.01 127.47 127.94 128.4 128.87 129.34 129.81 130.28 130.76 131.23 131.71 132.18 132.66 133.14 133.62 134.1 134.59 135.07 135.56 136.05 136.53 137.02 137.52 138.01 138.5 139 139.5 139.99 140.49 141 141.5 142 142.51 143.02 143.52 144.03 144.55 145.06 145.57 146.09 146.61 147.13 147.65 148.17 148.7 149.22 149.75 150.28 150.81 151.34 151.88 152.41 152.95 153.49 154.03 154.58 155.12 155.67 156.22 156.77 157.32 157.88 158.43 158.99 159.55 160.11 160.68 161.24 161.81 162.38 162.95 163.53 164.1 164.68 165.26 165.84 166.43 167.02 167.61 168.2 168.79 169.39 169.99 170.59 171.19 171.8 172.4 173.01 173.63 174.24 174.86 175.48 176.1 176.73 177.36 177.99 178.62 179.26 179.9 180.54 181.18 181.83 182.48 183.13 183.79 184.45 185.11 185.78 186.44 187.11 187.79 188.47 189.15 189.83 190.52 191.21 191.9 192.6 193.3 194.01 194.72 195.43 196.14 196.86 197.59 198.31 199.05 199.78 200.52 201.26 202.01 202.76 203.52 204.28 205.04 205.81 206.58 207.36 208.14 208.93 209.72 210.52 211.32 212.13 212.94 213.76 214.58 215.41 216.24 217.08 217.92 218.78 219.63 220.49 221.36 222.24 223.12 224.01 224.9 225.8 226.71 227.62 228.54 229.47 230.41 231.35 232.31 233.27 234.23 235.21 236.19 237.19 238.19 239.2 240.22 241.25 242.29 243.34 244.4 245.47 246.55 247.64 248.74 249.86 250.98 252.12 253.27 254.44 255.62 256.81 258.01 259.24 260.47 261.72 262.99 264.27 265.57 266.89 268.23 269.59 270.97 272.36 273.78 275.22 276.69 278.18 279.69 281.23 282.8 284.4 286.03 287.69 289.38 291.11 292.87 294.68 296.53 298.42 300.36 302.34 304.39 306.49 308.65 310.87 313.17 315.54 318 320.55 323.21 325.97 328.86 331.89 335.07 338.43 342 345.8 349.89 354.31 359.14 364.5 370.52 377.48 385.78 396.3 411.29 452.01]; + +N2O.sL = 1000*[-0.022012 -0.0024045 0.015246 0.031319 0.046093 0.059777 0.072533 0.084488 0.095745 0.10639 0.11648 0.12609 0.13527 0.14404 0.15245 0.16054 0.16832 0.17583 0.18308 0.19009 0.19687 0.20346 0.20985 0.21605 0.22209 0.22797 0.2337 0.23929 0.24474 0.25006 0.25527 0.26036 0.26534 0.27021 0.27499 0.27967 0.28426 0.28877 0.29319 0.29753 0.30179 0.30598 0.3101 0.31416 0.31814 0.32207 0.32593 0.32973 0.33348 0.33717 0.34081 0.3444 0.34794 0.35143 0.35487 0.35827 0.36162 0.36493 0.3682 0.37143 0.37462 0.37778 0.38089 0.38397 0.38702 0.39003 0.393 0.39595 0.39886 0.40174 0.4046 0.40742 0.41021 0.41298 0.41572 0.41843 0.42112 0.42378 0.42642 0.42903 0.43162 0.43418 0.43672 0.43924 0.44174 0.44422 0.44667 0.44911 0.45152 0.45391 0.45629 0.45864 0.46098 0.4633 0.4656 0.46789 0.47015 0.4724 0.47463 0.47685 0.47905 0.48123 0.4834 0.48555 0.48769 0.48981 0.49192 0.49402 0.4961 0.49816 0.50022 0.50226 0.50428 0.50629 0.5083 0.51028 0.51226 0.51422 0.51617 0.51811 0.52004 0.52196 0.52386 0.52576 0.52764 0.52951 0.53137 0.53322 0.53506 0.53689 0.53871 0.54052 0.54232 0.54411 0.5459 0.54767 0.54943 0.55118 0.55293 0.55466 0.55639 0.55811 0.55982 0.56152 0.56321 0.56489 0.56657 0.56824 0.5699 0.57155 0.57319 0.57483 0.57646 0.57808 0.5797 0.5813 0.5829 0.58449 0.58608 0.58766 0.58923 0.59079 0.59235 0.5939 0.59545 0.59699 0.59852 0.60004 0.60156 0.60308 0.60458 0.60608 0.60758 0.60907 0.61055 0.61203 0.6135 0.61496 0.61642 0.61788 0.61933 0.62077 0.62221 0.62364 0.62507 0.62649 0.62791 0.62932 0.63073 0.63213 0.63353 0.63492 0.6363 0.63769 0.63906 0.64044 0.6418 0.64317 0.64453 0.64588 0.64723 0.64857 0.64991 0.65125 0.65258 0.65391 0.65523 0.65655 0.65787 0.65918 0.66049 0.66179 0.66309 0.66438 0.66567 0.66696 0.66824 0.66952 0.6708 0.67207 0.67334 0.6746 0.67586 0.67712 0.67837 0.67962 0.68087 0.68211 0.68335 0.68458 0.68582 0.68704 0.68827 0.68949 0.69071 0.69193 0.69314 0.69435 0.69555 0.69676 0.69796 0.69915 0.70035 0.70154 0.70273 0.70391 0.70509 0.70627 0.70745 0.70862 0.70979 0.71096 0.71212 0.71328 0.71444 0.7156 0.71675 0.7179 0.71905 0.7202 0.72134 0.72248 0.72362 0.72475 0.72588 0.72701 0.72814 0.72927 0.73039 0.73151 0.73263 0.73374 0.73486 0.73597 0.73708 0.73818 0.73929 0.74039 0.74149 0.74258 0.74368 0.74477 0.74586 0.74695 0.74804 0.74912 0.7502 0.75128 0.75236 0.75344 0.75451 0.75558 0.75665 0.75772 0.75879 0.75985 0.76091 0.76197 0.76303 0.76409 0.76514 0.7662 0.76725 0.7683 0.76934 0.77039 0.77143 0.77248 0.77352 0.77455 0.77559 0.77663 0.77766 0.77869 0.77972 0.78075 0.78178 0.78281 0.78383 0.78485 0.78587 0.78689 0.78791 0.78893 0.78994 0.79096 0.79197 0.79298 0.79399 0.795 0.796 0.79701 0.79801 0.79902 0.80002 0.80102 0.80202 0.80301 0.80401 0.805 0.806 0.80699 0.80798 0.80897 0.80996 0.81095 0.81193 0.81292 0.8139 0.81489 0.81587 0.81685 0.81783 0.81881 0.81978 0.82076 0.82174 0.82271 0.82368 0.82466 0.82563 0.8266 0.82757 0.82854 0.8295 0.83047 0.83144 0.8324 0.83337 0.83433 0.83529 0.83626 0.83722 0.83818 0.83914 0.84009 0.84105 0.84201 0.84297 0.84392 0.84488 0.84583 0.84678 0.84774 0.84869 0.84964 0.85059 0.85154 0.85249 0.85344 0.85439 0.85534 0.85629 0.85723 0.85818 0.85913 0.86007 0.86102 0.86196 0.86291 0.86385 0.8648 0.86574 0.86668 0.86762 0.86857 0.86951 0.87045 0.87139 0.87233 0.87327 0.87421 0.87515 0.87609 0.87703 0.87797 0.87891 0.87985 0.88079 0.88173 0.88267 0.88361 0.88455 0.88549 0.88642 0.88736 0.8883 0.88924 0.89018 0.89112 0.89206 0.893 0.89394 0.89487 0.89581 0.89675 0.89769 0.89863 0.89957 0.90051 0.90146 0.9024 0.90334 0.90428 0.90522 0.90616 0.90711 0.90805 0.90899 0.90994 0.91088 0.91183 0.91277 0.91372 0.91467 0.91561 0.91656 0.91751 0.91846 0.91941 0.92036 0.92131 0.92227 0.92322 0.92417 0.92513 0.92608 0.92704 0.928 0.92896 0.92992 0.93088 0.93184 0.9328 0.93377 0.93473 0.9357 0.93667 0.93764 0.93861 0.93958 0.94055 0.94153 0.94251 0.94348 0.94446 0.94544 0.94643 0.94741 0.9484 0.94939 0.95038 0.95137 0.95236 0.95336 0.95436 0.95536 0.95636 0.95736 0.95837 0.95938 0.96039 0.96141 0.96242 0.96344 0.96446 0.96549 0.96652 0.96755 0.96858 0.96962 0.97066 0.9717 0.97275 0.97379 0.97485 0.9759 0.97696 0.97803 0.9791 0.98017 0.98124 0.98232 0.98341 0.9845 0.98559 0.98669 0.98779 0.9889 0.99001 0.99113 0.99226 0.99339 0.99452 0.99566 0.99681 0.99796 0.99912 1.0003 1.0015 1.0026 1.0038 1.005 1.0062 1.0074 1.0087 1.0099 1.0111 1.0124 1.0136 1.0149 1.0162 1.0175 1.0188 1.0201 1.0214 1.0228 1.0241 1.0255 1.0269 1.0283 1.0297 1.0311 1.0326 1.034 1.0355 1.037 1.0385 1.0401 1.0417 1.0433 1.0449 1.0466 1.0483 1.05 1.0518 1.0536 1.0554 1.0573 1.0592 1.0612 1.0632 1.0654 1.0675 1.0698 1.0721 1.0746 1.0771 1.0798 1.0826 1.0856 1.0888 1.0922 1.0959 1.1 1.1046 1.1099 1.1162 1.1241 1.1354 1.1682]; +N2O.sV = 1000*[2.0452 2.0286 2.0139 2.0008 1.9889 1.978 1.968 1.9588 1.9502 1.9421 1.9346 1.9274 1.9207 1.9143 1.9082 1.9024 1.8968 1.8915 1.8864 1.8815 1.8768 1.8723 1.8679 1.8637 1.8596 1.8556 1.8517 1.848 1.8444 1.8409 1.8374 1.8341 1.8308 1.8277 1.8246 1.8216 1.8186 1.8157 1.8129 1.8101 1.8074 1.8048 1.8022 1.7996 1.7971 1.7947 1.7923 1.7899 1.7876 1.7853 1.7831 1.7809 1.7787 1.7765 1.7744 1.7724 1.7703 1.7683 1.7664 1.7644 1.7625 1.7606 1.7587 1.7569 1.7551 1.7533 1.7515 1.7497 1.748 1.7463 1.7446 1.7429 1.7413 1.7397 1.738 1.7364 1.7349 1.7333 1.7318 1.7302 1.7287 1.7272 1.7258 1.7243 1.7228 1.7214 1.72 1.7186 1.7172 1.7158 1.7144 1.7131 1.7117 1.7104 1.7091 1.7078 1.7065 1.7052 1.7039 1.7026 1.7014 1.7001 1.6989 1.6976 1.6964 1.6952 1.694 1.6928 1.6916 1.6905 1.6893 1.6881 1.687 1.6858 1.6847 1.6836 1.6825 1.6813 1.6802 1.6791 1.678 1.677 1.6759 1.6748 1.6737 1.6727 1.6716 1.6706 1.6696 1.6685 1.6675 1.6665 1.6655 1.6644 1.6634 1.6624 1.6614 1.6605 1.6595 1.6585 1.6575 1.6566 1.6556 1.6546 1.6537 1.6527 1.6518 1.6508 1.6499 1.649 1.6481 1.6471 1.6462 1.6453 1.6444 1.6435 1.6426 1.6417 1.6408 1.6399 1.639 1.6381 1.6373 1.6364 1.6355 1.6346 1.6338 1.6329 1.6321 1.6312 1.6303 1.6295 1.6287 1.6278 1.627 1.6261 1.6253 1.6245 1.6236 1.6228 1.622 1.6212 1.6204 1.6195 1.6187 1.6179 1.6171 1.6163 1.6155 1.6147 1.6139 1.6131 1.6123 1.6116 1.6108 1.61 1.6092 1.6084 1.6077 1.6069 1.6061 1.6053 1.6046 1.6038 1.603 1.6023 1.6015 1.6008 1.6 1.5993 1.5985 1.5978 1.597 1.5963 1.5955 1.5948 1.594 1.5933 1.5926 1.5918 1.5911 1.5903 1.5896 1.5889 1.5882 1.5874 1.5867 1.586 1.5853 1.5845 1.5838 1.5831 1.5824 1.5817 1.581 1.5803 1.5795 1.5788 1.5781 1.5774 1.5767 1.576 1.5753 1.5746 1.5739 1.5732 1.5725 1.5718 1.5711 1.5704 1.5697 1.569 1.5683 1.5676 1.567 1.5663 1.5656 1.5649 1.5642 1.5635 1.5628 1.5622 1.5615 1.5608 1.5601 1.5594 1.5588 1.5581 1.5574 1.5567 1.556 1.5554 1.5547 1.554 1.5534 1.5527 1.552 1.5513 1.5507 1.55 1.5493 1.5487 1.548 1.5473 1.5467 1.546 1.5453 1.5447 1.544 1.5433 1.5427 1.542 1.5414 1.5407 1.54 1.5394 1.5387 1.538 1.5374 1.5367 1.5361 1.5354 1.5348 1.5341 1.5334 1.5328 1.5321 1.5315 1.5308 1.5302 1.5295 1.5288 1.5282 1.5275 1.5269 1.5262 1.5256 1.5249 1.5243 1.5236 1.5229 1.5223 1.5216 1.521 1.5203 1.5197 1.519 1.5184 1.5177 1.5171 1.5164 1.5158 1.5151 1.5145 1.5138 1.5131 1.5125 1.5118 1.5112 1.5105 1.5099 1.5092 1.5086 1.5079 1.5073 1.5066 1.506 1.5053 1.5046 1.504 1.5033 1.5027 1.502 1.5014 1.5007 1.5 1.4994 1.4987 1.4981 1.4974 1.4968 1.4961 1.4954 1.4948 1.4941 1.4935 1.4928 1.4921 1.4915 1.4908 1.4901 1.4895 1.4888 1.4882 1.4875 1.4868 1.4862 1.4855 1.4848 1.4842 1.4835 1.4828 1.4822 1.4815 1.4808 1.4801 1.4795 1.4788 1.4781 1.4774 1.4768 1.4761 1.4754 1.4747 1.4741 1.4734 1.4727 1.472 1.4713 1.4707 1.47 1.4693 1.4686 1.4679 1.4672 1.4665 1.4658 1.4652 1.4645 1.4638 1.4631 1.4624 1.4617 1.461 1.4603 1.4596 1.4589 1.4582 1.4575 1.4568 1.4561 1.4554 1.4547 1.454 1.4532 1.4525 1.4518 1.4511 1.4504 1.4497 1.4489 1.4482 1.4475 1.4468 1.446 1.4453 1.4446 1.4439 1.4431 1.4424 1.4417 1.4409 1.4402 1.4394 1.4387 1.438 1.4372 1.4365 1.4357 1.435 1.4342 1.4334 1.4327 1.4319 1.4312 1.4304 1.4296 1.4289 1.4281 1.4273 1.4265 1.4258 1.425 1.4242 1.4234 1.4226 1.4218 1.421 1.4202 1.4194 1.4186 1.4178 1.417 1.4162 1.4154 1.4146 1.4138 1.4129 1.4121 1.4113 1.4105 1.4096 1.4088 1.4079 1.4071 1.4063 1.4054 1.4045 1.4037 1.4028 1.402 1.4011 1.4002 1.3993 1.3985 1.3976 1.3967 1.3958 1.3949 1.394 1.3931 1.3922 1.3912 1.3903 1.3894 1.3885 1.3875 1.3866 1.3856 1.3847 1.3837 1.3828 1.3818 1.3808 1.3798 1.3788 1.3779 1.3769 1.3758 1.3748 1.3738 1.3728 1.3718 1.3707 1.3697 1.3686 1.3676 1.3665 1.3654 1.3643 1.3632 1.3621 1.361 1.3599 1.3588 1.3576 1.3565 1.3553 1.3541 1.3529 1.3517 1.3505 1.3493 1.3481 1.3469 1.3456 1.3443 1.3431 1.3418 1.3405 1.3391 1.3378 1.3364 1.3351 1.3337 1.3323 1.3308 1.3294 1.3279 1.3264 1.3249 1.3234 1.3218 1.3203 1.3186 1.317 1.3153 1.3137 1.3119 1.3102 1.3084 1.3065 1.3047 1.3028 1.3008 1.2988 1.2967 1.2946 1.2925 1.2902 1.2879 1.2855 1.2831 1.2805 1.2779 1.2751 1.2723 1.2693 1.2661 1.2627 1.2592 1.2554 1.2513 1.2469 1.242 1.2366 1.2304 1.223 1.2139 1.2012 1.1682]; + +N2O.hL = 1000*[-4.0502 -0.44503 2.8366 5.8554 8.6559 11.272 13.73 16.051 18.252 20.346 22.346 24.26 26.097 27.865 29.568 31.213 32.805 34.346 35.842 37.294 38.707 40.082 41.422 42.73 44.006 45.253 46.472 47.665 48.834 49.979 51.101 52.202 53.283 54.344 55.387 56.412 57.42 58.412 59.388 60.349 61.295 62.228 63.147 64.053 64.947 65.829 66.699 67.558 68.406 69.243 70.07 70.887 71.695 72.494 73.283 74.064 74.836 75.6 76.356 77.104 77.844 78.577 79.303 80.022 80.734 81.439 82.138 82.83 83.516 84.197 84.871 85.539 86.202 86.859 87.511 88.158 88.799 89.436 90.067 90.694 91.316 91.933 92.546 93.154 93.758 94.357 94.953 95.544 96.131 96.715 97.294 97.869 98.441 99.009 99.574 100.13 100.69 101.25 101.8 102.34 102.89 103.43 103.97 104.5 105.03 105.56 106.09 106.61 107.13 107.65 108.16 108.67 109.18 109.69 110.19 110.69 111.19 111.69 112.18 112.67 113.16 113.65 114.13 114.61 115.09 115.57 116.05 116.52 116.99 117.46 117.93 118.39 118.85 119.31 119.77 120.23 120.68 121.14 121.59 122.04 122.49 122.93 123.38 123.82 124.26 124.7 125.13 125.57 126 126.44 126.87 127.29 127.72 128.15 128.57 128.99 129.42 129.84 130.25 130.67 131.09 131.5 131.91 132.32 132.73 133.14 133.55 133.96 134.36 134.76 135.17 135.57 135.97 136.36 136.76 137.16 137.55 137.94 138.34 138.73 139.12 139.51 139.89 140.28 140.67 141.05 141.43 141.82 142.2 142.58 142.96 143.33 143.71 144.09 144.46 144.84 145.21 145.58 145.95 146.32 146.69 147.06 147.43 147.79 148.16 148.52 148.89 149.25 149.61 149.97 150.33 150.69 151.05 151.41 151.77 152.12 152.48 152.83 153.19 153.54 153.89 154.24 154.59 154.94 155.29 155.64 155.99 156.33 156.68 157.03 157.37 157.71 158.06 158.4 158.74 159.08 159.42 159.76 160.1 160.44 160.78 161.12 161.45 161.79 162.13 162.46 162.8 163.13 163.46 163.79 164.13 164.46 164.79 165.12 165.45 165.78 166.11 166.43 166.76 167.09 167.41 167.74 168.06 168.39 168.71 169.04 169.36 169.68 170 170.33 170.65 170.97 171.29 171.61 171.93 172.25 172.56 172.88 173.2 173.52 173.83 174.15 174.46 174.78 175.09 175.41 175.72 176.04 176.35 176.66 176.97 177.29 177.6 177.91 178.22 178.53 178.84 179.15 179.46 179.77 180.08 180.38 180.69 181 181.31 181.61 181.92 182.22 182.53 182.84 183.14 183.45 183.75 184.05 184.36 184.66 184.96 185.27 185.57 185.87 186.17 186.47 186.78 187.08 187.38 187.68 187.98 188.28 188.58 188.88 189.18 189.48 189.77 190.07 190.37 190.67 190.97 191.26 191.56 191.86 192.16 192.45 192.75 193.05 193.34 193.64 193.93 194.23 194.52 194.82 195.11 195.41 195.7 196 196.29 196.58 196.88 197.17 197.47 197.76 198.05 198.35 198.64 198.93 199.22 199.52 199.81 200.1 200.39 200.68 200.98 201.27 201.56 201.85 202.14 202.43 202.72 203.02 203.31 203.6 203.89 204.18 204.47 204.76 205.05 205.34 205.63 205.92 206.21 206.5 206.79 207.08 207.37 207.66 207.95 208.24 208.53 208.82 209.11 209.4 209.69 209.98 210.27 210.56 210.85 211.14 211.43 211.72 212.01 212.31 212.6 212.89 213.18 213.47 213.76 214.05 214.34 214.63 214.92 215.21 215.5 215.79 216.08 216.37 216.66 216.96 217.25 217.54 217.83 218.12 218.41 218.71 219 219.29 219.58 219.88 220.17 220.46 220.76 221.05 221.34 221.64 221.93 222.22 222.52 222.81 223.11 223.4 223.7 223.99 224.29 224.59 224.88 225.18 225.48 225.77 226.07 226.37 226.67 226.97 227.26 227.56 227.86 228.16 228.46 228.76 229.06 229.37 229.67 229.97 230.27 230.58 230.88 231.18 231.49 231.79 232.1 232.4 232.71 233.02 233.32 233.63 233.94 234.25 234.56 234.87 235.18 235.49 235.8 236.12 236.43 236.74 237.06 237.37 237.69 238.01 238.32 238.64 238.96 239.28 239.6 239.92 240.25 240.57 240.89 241.22 241.54 241.87 242.2 242.53 242.86 243.19 243.52 243.86 244.19 244.53 244.86 245.2 245.54 245.88 246.22 246.57 246.91 247.26 247.61 247.95 248.31 248.66 249.01 249.37 249.72 250.08 250.44 250.81 251.17 251.54 251.91 252.28 252.65 253.03 253.41 253.79 254.17 254.55 254.94 255.33 255.73 256.12 256.52 256.93 257.33 257.74 258.15 258.57 258.99 259.41 259.84 260.28 260.71 261.15 261.6 262.05 262.51 262.97 263.44 263.92 264.4 264.89 265.38 265.89 266.4 266.92 267.45 267.99 268.54 269.1 269.67 270.26 270.86 271.48 272.11 272.76 273.43 274.12 274.83 275.58 276.35 277.16 278 278.89 279.84 280.84 281.92 283.09 284.38 285.83 287.48 289.44 291.92 295.45 305.6]; +N2O.hV = 1000*[372.87 374.13 375.26 376.28 377.22 378.08 378.88 379.63 380.33 380.99 381.62 382.2 382.76 383.29 383.8 384.29 384.75 385.19 385.62 386.03 386.42 386.8 387.16 387.52 387.86 388.19 388.5 388.81 389.11 389.4 389.68 389.96 390.22 390.48 390.73 390.97 391.21 391.44 391.67 391.88 392.1 392.31 392.51 392.71 392.9 393.09 393.27 393.45 393.63 393.8 393.97 394.13 394.29 394.45 394.6 394.75 394.9 395.04 395.18 395.32 395.45 395.58 395.71 395.84 395.96 396.08 396.2 396.31 396.43 396.54 396.65 396.75 396.86 396.96 397.06 397.16 397.25 397.35 397.44 397.53 397.62 397.7 397.79 397.87 397.95 398.03 398.11 398.18 398.26 398.33 398.4 398.47 398.54 398.61 398.67 398.74 398.8 398.86 398.92 398.98 399.03 399.09 399.15 399.2 399.25 399.3 399.35 399.4 399.45 399.49 399.54 399.58 399.62 399.67 399.71 399.75 399.79 399.82 399.86 399.89 399.93 399.96 399.99 400.03 400.06 400.09 400.12 400.14 400.17 400.2 400.22 400.24 400.27 400.29 400.31 400.33 400.35 400.37 400.39 400.41 400.42 400.44 400.46 400.47 400.48 400.5 400.51 400.52 400.53 400.54 400.55 400.56 400.56 400.57 400.58 400.58 400.59 400.59 400.59 400.6 400.6 400.6 400.6 400.6 400.6 400.6 400.6 400.59 400.59 400.59 400.58 400.58 400.57 400.56 400.56 400.55 400.54 400.53 400.52 400.51 400.5 400.49 400.48 400.47 400.45 400.44 400.43 400.41 400.4 400.38 400.37 400.35 400.33 400.31 400.29 400.28 400.26 400.24 400.22 400.19 400.17 400.15 400.13 400.1 400.08 400.06 400.03 400.01 399.98 399.95 399.93 399.9 399.87 399.84 399.81 399.78 399.75 399.72 399.69 399.66 399.63 399.6 399.57 399.53 399.5 399.46 399.43 399.39 399.36 399.32 399.29 399.25 399.21 399.17 399.13 399.09 399.06 399.02 398.98 398.93 398.89 398.85 398.81 398.77 398.72 398.68 398.64 398.59 398.55 398.5 398.45 398.41 398.36 398.31 398.27 398.22 398.17 398.12 398.07 398.02 397.97 397.92 397.87 397.82 397.76 397.71 397.66 397.61 397.55 397.5 397.44 397.39 397.33 397.27 397.22 397.16 397.1 397.05 396.99 396.93 396.87 396.81 396.75 396.69 396.63 396.57 396.5 396.44 396.38 396.32 396.25 396.19 396.12 396.06 395.99 395.93 395.86 395.79 395.73 395.66 395.59 395.52 395.45 395.38 395.31 395.24 395.17 395.1 395.03 394.96 394.88 394.81 394.74 394.66 394.59 394.51 394.44 394.36 394.29 394.21 394.13 394.05 393.98 393.9 393.82 393.74 393.66 393.58 393.5 393.42 393.33 393.25 393.17 393.09 393 392.92 392.83 392.75 392.66 392.58 392.49 392.4 392.31 392.23 392.14 392.05 391.96 391.87 391.78 391.69 391.6 391.5 391.41 391.32 391.22 391.13 391.04 390.94 390.85 390.75 390.65 390.56 390.46 390.36 390.26 390.16 390.06 389.96 389.86 389.76 389.66 389.56 389.45 389.35 389.25 389.14 389.04 388.93 388.82 388.72 388.61 388.5 388.39 388.28 388.18 388.06 387.95 387.84 387.73 387.62 387.51 387.39 387.28 387.16 387.05 386.93 386.81 386.7 386.58 386.46 386.34 386.22 386.1 385.98 385.86 385.74 385.61 385.49 385.36 385.24 385.11 384.99 384.86 384.73 384.6 384.48 384.35 384.22 384.08 383.95 383.82 383.69 383.55 383.42 383.28 383.15 383.01 382.87 382.73 382.6 382.46 382.31 382.17 382.03 381.89 381.74 381.6 381.45 381.31 381.16 381.01 380.86 380.71 380.56 380.41 380.26 380.11 379.95 379.8 379.64 379.49 379.33 379.17 379.01 378.85 378.69 378.53 378.36 378.2 378.03 377.87 377.7 377.53 377.36 377.19 377.02 376.85 376.68 376.5 376.33 376.15 375.97 375.79 375.61 375.43 375.25 375.06 374.88 374.69 374.51 374.32 374.13 373.94 373.74 373.55 373.35 373.16 372.96 372.76 372.56 372.36 372.15 371.95 371.74 371.53 371.32 371.11 370.9 370.69 370.47 370.25 370.03 369.81 369.59 369.37 369.14 368.91 368.68 368.45 368.22 367.98 367.74 367.5 367.26 367.02 366.77 366.52 366.27 366.02 365.77 365.51 365.25 364.99 364.72 364.46 364.19 363.92 363.64 363.36 363.08 362.8 362.51 362.22 361.93 361.64 361.34 361.03 360.73 360.42 360.11 359.79 359.47 359.15 358.82 358.49 358.15 357.81 357.47 357.12 356.76 356.4 356.04 355.67 355.29 354.91 354.53 354.13 353.73 353.33 352.92 352.5 352.07 351.64 351.19 350.74 350.28 349.81 349.33 348.85 348.35 347.84 347.31 346.78 346.23 345.67 345.09 344.5 343.89 343.26 342.61 341.94 341.24 340.52 339.77 339 338.18 337.33 336.44 335.5 334.5 333.44 332.3 331.08 329.74 328.27 326.62 324.73 322.5 319.7 315.78 305.6]; + +%% OXYDIZER CRITICAL RATIO +% load fChockParamFun +% N2O.eta_cr = f_chock_param; +% clear fChockParamFun + +N2O.eta_cr = @(mox) 4.555e-6 * mox.^2 - 0.002462 * mox + 0.8302; + +end \ No newline at end of file diff --git a/Roccaraso/postprocessing/RPS/src/dataLoader.m b/Roccaraso/postprocessing/RPS/src/dataLoader.m index d47a651bdc192d0a4ed4c99abd1b6e794f108a54..e884a7dd4414bcc1dd3180d3f5ab78be24f728a8 100644 --- a/Roccaraso/postprocessing/RPS/src/dataLoader.m +++ b/Roccaraso/postprocessing/RPS/src/dataLoader.m @@ -26,7 +26,7 @@ rawData.timeSrvMain = rawData.timeSrv(rawData.srvOpening(:, 1) == 5); rawData.srvVent = rawData.srvOpening(rawData.srvOpening(:, 1) == 6, 2); rawData.timeSrvVent = rawData.timeSrv(rawData.srvOpening(:, 1) == 6); -if strcmp(settings.testName, 'ROC-01') +if strcmp(settings.testName, 'ROC1') % Defining the burning time tb = 6.2; @@ -51,7 +51,7 @@ if strcmp(settings.testName, 'ROC-01') flightData.srvMain = rawData.srvMain(is:js); flightData.srvVent = rawData.srvVent(is:js); -elseif strcmp(settings.testName, 'ROC-02') +elseif strcmp(settings.testName, 'ROC2') % Defining the burning time tb = 3.3; @@ -82,12 +82,4 @@ flightData.PTankTop = flightData.PTankTop + 1; flightData.PTankBot = flightData.PTankBot + 1; flightData.PEng = flightData.PEng + 1; -% Computing the thrust -flightData.Thrust = flightData.PEng * pi*settings.engine.rt^2 * settings.engine.cT * settings.engine.lam * 1e5; -flightData.timeThrust = flightData.timePEng; - -% Computing the propellant mass flow rate -flightData.mp = settings.engine.Isp * settings.g0 ./flightData.Thrust; -flightData.timeMp = flightData.timePEng; - save(strcat('cleanData/', settings.testName,'/engineFlightData.mat'), 'flightData') \ No newline at end of file diff --git a/Roccaraso/postprocessing/RPS/src/tankEvacuation.m b/Roccaraso/postprocessing/RPS/src/tankEvacuation.m new file mode 100644 index 0000000000000000000000000000000000000000..8e90e499b58e0fb049d3800c89be0a395b724464 --- /dev/null +++ b/Roccaraso/postprocessing/RPS/src/tankEvacuation.m @@ -0,0 +1,69 @@ +function engine = tankEvacuation(engine, settings) +%{ +tankEvacuation - this function outputs oxydizer mass flow rate usinìg NHNE +model + +CALLED SCRIPTS: + +CALLED FUNCTIONS: + +CALLED DATA FILES: + +REVISIONS: +- 0, 12/09/2022, Release, Alberto Boffi, Paolo Gnata, PRP +Department + +Copyright © 2022, Skyward Experimental Rocketry, PRP department +All rights reserved + +SPDX-License-Identifier: GPL-3.0-or-later + +%} +%% INPUT + +d_pipe = engine.data.d_pipe; +N2O = engine.data.ox; +pTank = engine.acquisition.PTankTop; +Aox = engine.data.Aox; +cd = engine.data.cd_inj; +alpha = 1; +time = engine.acquisition.timePTop; + +%% CALCULATIONS + +% estimate injection pressure from tank pressure +p1 = pTank; +% oxidizer mass flow rate initialization +mox = zeros(length(p1), 1); +% velocity initialization +v = 10; + +for i = 1 : length(p1) + rho1(i) = interp1(N2O.P, N2O.rhoL, p1(i)); + p1(i) = pTank(i) - (alpha/(2*Aox.^2*rho1(i)).*mox(i).^2)*1e-5; + s1(i) = interp1(N2O.P, N2O.sL, p1(i)); + h1(i) = interp1(N2O.P, N2O.hL, p1(i)); + p2(i) = p1(i) * N2O.eta_cr(p1(i)); + s2(i) = s1(i); + sV2(i) = interp1(N2O.P, N2O.sV, p2(i)); + sL2(i) = interp1(N2O.P, N2O.sL, p2(i)); + x2(i) = (s2(i)-sL2(i))/(sV2(i)-sL2(i)); + hV2(i) = interp1(N2O.P, N2O.hV, p2(i)); + hL2(i) = interp1(N2O.P, N2O.hL, p2(i)); + h2(i) = x2(i) * hV2(i) + (1-x2(i))*hL2(i); + rhoL2(i) = interp1(N2O.P, N2O.rhoL, p2(i)); + rhoV2(i) = interp1(N2O.P, N2O.rhoV, p2(i)); + rho2(i) = 1/(x2(i) * 1/rhoV2(i) + (1-x2(i))*1/rhoL2(i)); + k(i) = sqrt((p1(i)-p2(i))./(pTank(i)-p2(i))); + mox(i) = Aox * cd * (k(i)*sqrt(2*abs(p1(i)-p2(i))*1e5*rho1(i)) + ... + rho2(i)*sqrt(2*abs(h1(i)-h2(i))))/(1+k(i)); + v(i+1) = mox(i)./rho1(i)/(pi*d_pipe^2/4); +end + +%% OUTPUT + +% ozydizer mass flow rate +engine.performances.mox = mox; + + +