Move Ref-Data & Add my first shot

data/referance-data/x2s4111/databox/x2s4111.sh

@@ -0,0 +1,5 @@
+# x2s4111.sh:
+# shell script for x2s4111
+# 12/10/2018
+
+shot_class.py --shot_name x2s4111 --condition_title "80%Ne condition" --plot --analyse --error --lowcut 100000.0 --filter --get_canny_triggers --do_not_run_pitot_theoretical_eq --plot_canny_plots

data/referance-data/x2s4111/databox/x2s4111.shotcfg

@@ -0,0 +1,37 @@
+# Config file for x2s3690
+
+#details of the condition
+
+driver_condition = 'He:0.8,Ar:0.2'
+test_gas = 'gasgiant_h285ne'
+secondary = False
+nozzle = True
+conehead = True
+area_ratio = 5.64
+piston = 'lwp-1.2mm'
+shock_over_model = True
+
+#fill pressures (in Pa)
+nominal_shock_tube_fill = 2000.0
+nominal_acc_tube_fill = 1.0
+
+actual_shock_tube_fill = nominal_shock_tube_fill
+actual_acc_tube_fill = nominal_acc_tube_fill
+
+#any notes about the condition
+
+notes = "None."
+
+canny_settings_dict = {}
+
+canny_settings_dict['at'] = {'sigma':3, 'derivative_threshold':0.001, 'post_suppression_threshold':0.04}
+# supposedly this is a way to copy dictionaries in Python...
+# copy the at dictionary and then add a start time which is past the mechanical noise
+canny_settings_dict['at7'] = dict(canny_settings_dict['at'])
+canny_settings_dict['at7']['start_time'] = 0.00180
+canny_settings_dict['at8'] = canny_settings_dict['at7']
+canny_settings_dict['sd'] = {'sigma':2, 'derivative_threshold':1.0, 'post_suppression_threshold':10.0}
+canny_settings_dict['sd3'] = {'sigma':2, 'derivative_threshold':0.1, 'post_suppression_threshold':10.0}
+canny_settings_dict['st'] = canny_settings_dict['at']
+canny_settings_dict['pt'] = {'sigma':3, 'derivative_threshold':0.001, 'post_suppression_threshold':0.07,'start_time':0.00180}
+canny_settings_dict['trigbox'] = {'sigma':1, 'derivative_threshold':10.0, 'post_suppression_threshold':100.0,'start_time':0.00150}

data/referance-data/x2s4111/databox/x2s4111.txt

@@ -0,0 +1,86 @@
+Project.............................. Yu Liu -- Gas giant Neon condition
+Run number........................... x2s4111
+Date................................. 13/10/2018
+Blame................................ Yu Liu, Chris James
+Reservoir...........................  4.85 MPa compressed air
+Driver............................... 80% He,20% Ar filled to 1103 mba. Flushed twice with He
+Primary diaphragm.................... 1.2 mm hot rolled steel scored diaphragm 
+Shock tube........................... 2000 Pa 80% Ne, 20% H2 by volume, flushed twice with He
+Secondary diaphragm.................. 12 micron Al
+Acceleration tube.................... 1 Pa lab air
+Trig................................. Triggering off pt2, trigger level set to 5 on box 
+
+EXPERIMENT
+
+Gas Giant Helium condition.
+
+Centre IR 1500 nm 
+
+Centre UV 525 nm
+
+High speed camera in use
+
+---------------------
+-- UV spectrometer --
+---------------------
+
+Centre wavelength = 525 nm
+
+Aperture = 20 mm
+
+Delay = 200 us
+
+Exposure = 20 us
+
+Slit width = 50 um
+
+Filter = 400 nm longpass filter
+
+Gain: 2800
+
+Grating: 150
+
+---------------------
+-- IR spectrometer --
+---------------------
+
+Centre wavelength = 1500 nm
+
+Aperture = 20 mm
+
+Delay = triggering off Andor
+
+Exposure = 20 us
+
+Slit width = 250 um
+
+Filter = 850 nm longpass filter FGL850S, I think
+
+-----------------------
+-- High speed camera --
+-----------------------
+
+Andover corp H-beta filter
+
+f-number: 8
+
+Framerate: 250 kHz (4 us between frames)
+
+Exposure : 1/2
+
+Pre-frames = 10
+
+Delay = 0 us
+
+Gain: 2
+
+-----------------------
+-- Thermocouple mount -
+-----------------------
+
+Channel: 1
+Gain: 100
+Resistance 3.6 - 3.7 Ohms (may change slightly from shot to sjot, but being checked...
+New thermocouple with ground wire
+
+

data/referance-data/x2s4111/databox/x2s4111.xconfig

@@ -0,0 +1,40 @@
+# x2s4111 config file
+# This files contains a text version of the National Instruments PXI unit configuration data.
+# Channels are represented as one line per configured signal with tab-separated values.
+#
+# Line content and types of values (first item is item 0):
+# i  item            type        comments
+# --------------------------------------------------------------------------------------------------------
+# 0  signal name     string    
+# 1  channel ID      string    legacy I/O channel designation
+# 2  external gain   float     1 if there is no external amplifier connected
+# 3  sensitivity     float     volts per measured unit (of the transducer)
+# 4  units           string    name of units for physical quantity (e.g. kPa)
+# 5  position        float     position of transducer (in mm)
+# 6  serial number   string    transducer identity 
+# 7  PXI channel ID  string    non-legacy channel ID and trigger source
+#---------------------------------------------------------------------------------------------------------
+
+sd1            	2 0 9	     1	7.321500E-4	kPa	0.000	113B24-LW30646 	PXI1Slot2/ai0	-1.000	1.000
+sd2            	2 7 9	     1	7.259100E-4	kPa	0.000	113B24-LW30648 	PXI1Slot2/ai7	-1.000	1.000
+sd3            	2 2 9	     1	6.820000E-4	kPa	0.000	112A24-2676    	PXI1Slot2/ai2	-1.000	1.000
+st1            	5 2 9	     1	1.239800E-2	kPa	0.000	112A24-9594    	PXI1Slot5/ai2	-1.000	1.000
+st2            	2 4 9	     1	1.187400E-2	kPa	0.000	112A24-19126   	PXI1Slot2/ai4	-1.000	1.000
+st3            	2 5 9	     1	1.482000E-2	kPa	0.000	112A24-35071   	PXI1Slot2/ai5	-1.000	1.000
+at1            	3 0 9	     1	1.443000E-2	kPa	0.000	112A22-35098   	PXI1Slot3/ai0	-1.000	1.000
+at2            	3 1 9	     1	1.450000E-2	kPa	0.000	112A22-9050    	PXI1Slot3/ai1	-1.000	1.000
+at3            	3 2 9	     1	1.447980E-2	kPa	0.000	112A22-19124   	PXI1Slot3/ai2	-1.000	1.000
+at4            	3 3 9	     1	1.435000E-2	kPa	0.000	112A22-34424   	PXI1Slot3/ai3	-1.000	1.000
+at5            	3 4 9	     1	1.393600E-2	kPa	0.000	112A22-35042   	PXI1Slot3/ai4	-1.000	1.000
+at6            	3 5 9	     1	1.442000E-2	kPa	0.000	112A22-34426   	PXI1Slot3/ai5	-1.000	1.000
+at7            	3 6 9	     1	1.451000E-2	kPa	0.000	112A22-35096   	PXI1Slot3/ai6	-1.000	1.000
+at8            	3 7 9	     1	1.490000E-2	kPa	0.000	112A22-35079   	PXI1Slot3/ai7	-1.000	1.000
+trigbox        	5 0 9	     1	1.000000E-3	mV	0.000	trigbox        	PXI1Slot5/ai0(Trigger)	-10.000	10.000
+photodiode     	5 1 9	     1	1.000000E-3	mV	0.000	photodiode     	PXI1Slot5/ai1	-10.000	10.000
+hpv1           	5 6 9	     1	1.100000E-3	mV	0.000	hpv1           	PXI1Slot5/ai6	-10.000	10.000
+pt1            	4 0 9	     1	1.464000E-2	kPa	0.000	112A22-36656   	PXI1Slot4/ai0	-10.000	10.000
+pt2            	4 1 9	     1	1.212200E-2	kPa	0.000	112A22-35073   	PXI1Slot4/ai1	-10.000	10.000
+pt3            	4 2 9	     1	1.448000E-2	kPa	0.000	112A22-36657   	PXI1Slot4/ai2	-10.000	10.000
+mwirout        	5 7 9	     1	1.000000E-3	mV	0.000	mwirout        	PXI1Slot5/ai7	-10.000	10.000
+tctest         	5 3 9	     1	1.000000E-3	mV	0.000	tctest         	PXI1Slot5/ai3	-5.000	5.000
+uvout          	5 4 9	     1	1.000000E-3	mV	0.000	uvout          	PXI1Slot5/ai4	-10.000	10.000

data/referance-data/x2s4111/databox/x2s4111_basic_data_summary.txt

@@ -0,0 +1,78 @@
+# Basic analysis performed using xlabs shot data analysis code version 15-Aug-2018
+------------------------------------------------------------
+Shot x2s4111 - 80%Ne condition
+------------------------------------------------------------
+# Curve fit shock speeds (with 1-sigma uncertainty)
+st1-at8 = 4883.79 +- 732.65 m/s (using sensors ['st1', 'st2', 'st3', 'at1', 'at2', 'at3', 'at4', 'at5', 'at6', 'at7', 'at8'])
+at1-at8 = 3261.69 +- 943.57 m/s (using sensors ['at1', 'at2', 'at3', 'at4', 'at5', 'at6', 'at7', 'at8'])
+at4-at8 = 1078.67 +- 730.94 m/s (using sensors ['at4', 'at5', 'at6', 'at7', 'at8'])
+------------------------------------------------------------
+# Trigger (shock arrival) times
+sd1 = 0.000637 +- 2.00e-07 s (0.0006368 - 0.0006372 s)
+sd2 = 0.0007052 +- 0.00e+00 s (0.0007052 - 0.0007052 s)
+sd3 = 0.0007814 +- 2.00e-07 s (0.0007812 - 0.0007816 s)
+st1 = 0.0010636 +- 0.00e+00 s (0.0010636 - 0.0010636 s)
+st2 = 0.0011468 +- 0.00e+00 s (0.0011468 - 0.0011468 s)
+st3 = 0.001313 +- 2.00e-07 s (0.0013128 - 0.0013132 s)
+at1 = 0.0014216 +- 0.00e+00 s (0.0014216 - 0.0014216 s)
+at2 = 0.0014506 +- 2.00e-07 s (0.0014504 - 0.0014508 s)
+at3 = 0.0014786 +- 2.00e-07 s (0.0014784 - 0.0014788 s)
+at4 = 0.0016018 +- 2.00e-07 s (0.0016016 - 0.001602 s)
+at5 = 0.0016406 +- 2.00e-07 s (0.0016404 - 0.0016408 s)
+at6 = 0.0016778 +- 2.00e-07 s (0.0016776 - 0.001678 s)
+at7 = 0.0019916 +- 0.00e+00 s (0.0019916 - 0.0019916 s)
+at8 = 0.0020326 +- 2.00e-07 s (0.0020324 - 0.0020328 s)
+pt1 = 0.0019722 +- 2.00e-07 s (0.001972 - 0.0019724 s)
+pt2 = 0.0019726 +- 2.00e-07 s (0.0019724 - 0.0019728 s)
+pt3 = 0.0019698 +- 2.00e-07 s (0.0019696 - 0.00197 s)
+------------------------------------------------------------
+# Experimental Time-of-flight Shock speeds
+sd1-sd2 = 3416.42 +- 46.06 m/s (1.35 %)
+sd1-sd3 = 3227.15 +- 22.44 m/s (0.70 %)
+sd2-sd3 = 3057.74 +- 42.00 m/s (1.37 %)
+st1-st3 = 4125.90 +- 13.94 m/s (0.34 %)
+st1-st2 = 6189.90 +- 45.18 m/s (0.73 %)
+st2-st3 = 3092.66 +- 19.31 m/s (0.62 %)
+at1-at3 = 6298.25 +- 73.43 m/s (1.17 %)
+at1-at2 = 6137.93 +- 142.35 m/s (2.32 %)
+at2-at3 = 6464.29 +- 151.64 m/s (2.35 %)
+at3-at4 = 6444.81 +- 34.41 m/s (0.53 %)
+at4-at6 = 6657.89 +- 56.81 m/s (0.85 %)
+at4-at5 = 6597.94 +- 110.70 m/s (1.68 %)
+at5-at7 = 886.04 +- 8.12 m/s (0.92 %)
+at5-at6 = 6720.43 +- 116.68 m/s (1.74 %)
+at6-at8 = 284.67 +- 7.98 m/s (2.80 %)
+at6-at7 = 194.39 +- 9.02 m/s (4.64 %)
+at7-at8 = 975.61 +- 69.96 m/s (7.17 %)
+------------------------------------------------------------
+Test section flow arrival time is 0.00197153333333 s (1971.53 microseconds).
+Found using transducers ['pt1', 'pt2', 'pt3'].
+------------------------------------------------------------
+trigbox arrival at 0.0019998 s (1999.80 microseconds).
+This is 28.27 microseconds after flow arrival.
+------------------------------------------------------------
+# Channel names
+sd1
+sd2
+sd3
+st1
+st2
+st3
+at1
+at2
+at3
+at4
+at5
+at6
+at7
+at8
+pt1
+pt2
+pt3
+hpv1
+mwirout
+photodiode
+tctest
+trigbox
+uvout
+------------------------------------------------------------

data/referance-data/x2s4111/databox/x2s4111_full_data_summary.txt

@@ -0,0 +1,116 @@
+# Full analysis performed using xlabs shot data analysis code version 15-Aug-2018
+------------------------------------------------------------
+Shot x2s4111 - 80%Ne condition
+Driver condition is He:0.8,Ar:0.2.
+Test gas is gasgiant_h285ne.
+The X2 nozzle is in use.
+------------------------------------------------------------
+Shock tube nominal fill pressure is 2000.0 Pa, actual fill pressure was 2000.0 Pa.
+Acceleration tube nominal fill pressure is 1.0 Pa, actual fill pressure was 1.0 Pa.
+------------------------------------------------------------
+Any notes about the shot are:
+None.
+------------------------------------------------------------
+# Curve fit shock speeds (with 1-sigma uncertainty)
+st1-at8 = 4883.79 +- 732.65 m/s (using sensors ['st1', 'st2', 'st3', 'at1', 'at2', 'at3', 'at4', 'at5', 'at6', 'at7', 'at8'])
+at1-at8 = 3261.69 +- 943.57 m/s (using sensors ['at1', 'at2', 'at3', 'at4', 'at5', 'at6', 'at7', 'at8'])
+at4-at8 = 1078.67 +- 730.94 m/s (using sensors ['at4', 'at5', 'at6', 'at7', 'at8'])
+------------------------------------------------------------
+# Trigger (shock arrival) times
+sd1 = 0.000637 +- 2.00e-07 s (0.0006368 - 0.0006372 s)
+sd2 = 0.0007052 +- 0.00e+00 s (0.0007052 - 0.0007052 s)
+sd3 = 0.0007814 +- 2.00e-07 s (0.0007812 - 0.0007816 s)
+st1 = 0.0010636 +- 0.00e+00 s (0.0010636 - 0.0010636 s)
+st2 = 0.0011468 +- 0.00e+00 s (0.0011468 - 0.0011468 s)
+st3 = 0.001313 +- 2.00e-07 s (0.0013128 - 0.0013132 s)
+at1 = 0.0014216 +- 0.00e+00 s (0.0014216 - 0.0014216 s)
+at2 = 0.0014506 +- 2.00e-07 s (0.0014504 - 0.0014508 s)
+at3 = 0.0014786 +- 2.00e-07 s (0.0014784 - 0.0014788 s)
+at4 = 0.0016018 +- 2.00e-07 s (0.0016016 - 0.001602 s)
+at5 = 0.0016406 +- 2.00e-07 s (0.0016404 - 0.0016408 s)
+at6 = 0.0016778 +- 2.00e-07 s (0.0016776 - 0.001678 s)
+at7 = 0.0019916 +- 0.00e+00 s (0.0019916 - 0.0019916 s)
+at8 = 0.0020326 +- 2.00e-07 s (0.0020324 - 0.0020328 s)
+pt1 = 0.0019722 +- 2.00e-07 s (0.001972 - 0.0019724 s)
+pt2 = 0.0019726 +- 2.00e-07 s (0.0019724 - 0.0019728 s)
+pt3 = 0.0019698 +- 2.00e-07 s (0.0019696 - 0.00197 s)
+------------------------------------------------------------
+# Experimental Time-of-flight Shock speeds
+sd1-sd2 = 3416.42 +- 46.06 m/s (1.35 %)
+sd1-sd3 = 3227.15 +- 22.44 m/s (0.70 %)
+sd2-sd3 = 3057.74 +- 42.00 m/s (1.37 %)
+st1-st3 = 4125.90 +- 13.94 m/s (0.34 %)
+st1-st2 = 6189.90 +- 45.18 m/s (0.73 %)
+st2-st3 = 3092.66 +- 19.31 m/s (0.62 %)
+at1-at3 = 6298.25 +- 73.43 m/s (1.17 %)
+at1-at2 = 6137.93 +- 142.35 m/s (2.32 %)
+at2-at3 = 6464.29 +- 151.64 m/s (2.35 %)
+at3-at4 = 6444.81 +- 34.41 m/s (0.53 %)
+at4-at6 = 6657.89 +- 56.81 m/s (0.85 %)
+at4-at5 = 6597.94 +- 110.70 m/s (1.68 %)
+at5-at7 = 886.04 +- 8.12 m/s (0.92 %)
+at5-at6 = 6720.43 +- 116.68 m/s (1.74 %)
+at6-at8 = 284.67 +- 7.98 m/s (2.80 %)
+at6-at7 = 194.39 +- 9.02 m/s (4.64 %)
+at7-at8 = 975.61 +- 69.96 m/s (7.17 %)
+------------------------------------------------------------
+Test section flow arrival time is 0.00197153333333 s (1971.53 microseconds).
+Found using transducers ['pt1', 'pt2', 'pt3'].
+------------------------------------------------------------
+trigbox arrival at 0.0019998 s (1999.80 microseconds).
+This is 28.27 microseconds after flow arrival.
+------------------------------------------------------------
+Now printing basic results of experimental equilibrium run of Pitot program.
+Vs1 = 3057.74 m/s.
+Vs2 = 6720.43 m/s.
+Nozzle entry pressure (p7) 125.928708702 Pa.
+Nozzle exit pressure (p8) 7.59 Pa (using area ratio of 5.64).
+Conehead pressure is 488.514417195 Pa.
+------------------------------------------------------------
+Now printing basic results of theoretical perfect gas run of Pitot program.
+Vs1 = 4521.09 m/s.
+Vs2 = 10323.69 m/s.
+Nozzle entry pressure (p7) 12.93 Pa.
+Nozzle exit pressure (p8) 0.80 Pa (using area ratio of 5.64).
+Conehead pressure is 184.50 Pa.
+------------------------------------------------------------
+# Experimental hold time data from the x-t diagram
+All hold time data is compared to the nominal solution on the experimental x-t diagram.
+------------------------------------------------------------
+Using ['st1', 'st2', 'st3'] the hold time is -5.93262347716e-05 s (-59.3262 microseconds).
+The slope and intercept used to find this are 0.000242345028717 and 0.000230156482825.
+From the linear regression, R**2 = 0.964187878845, the p value is 0.121205419603, the std err is 4.6705529942e-05.
+------------------------------------------------------------
+Using ['at1', 'at2', 'at3'] the hold time is 3.0503368026e-05 s (30.5034 microseconds).
+The slope and intercept used to find this are 0.000158762918594 and 0.000605669738024.
+From the linear regression, R**2 = 0.999776417807, the p value is 0.00951951700366, the std err is 2.37419417977e-06.
+------------------------------------------------------------
+Using ['at4', 'at5', 'at6'] the hold time is 6.32283189752e-05 s (63.2283 microseconds).
+The slope and intercept used to find this are 0.000150203086921 and 0.000667652193629.
+From the linear regression, R**2 = 0.999971823274, the p value is 0.00337930527904, the std err is 7.9731556264e-07.
+------------------------------------------------------------
+# Channel names
+sd1
+sd2
+sd3
+st1
+st2
+st3
+at1
+at2
+at3
+at4
+at5
+at6
+at7
+at8
+pt1
+pt2
+pt3
+hpv1
+mwirout
+photodiode
+tctest
+trigbox
+uvout
+------------------------------------------------------------

data/referance-data/x2s4111/databox/x2s4111eq_experimental-x-t-diagram-summary.txt

@@ -0,0 +1,19 @@
+# x-t diagram summary made using Pitot Version 15-Aug-2018.
+------------------------------------------------------------
+Shock will pass transducer sd1 at t = 0.000842778540772 s (842.778540772 microseconds).
+Shock will pass transducer sd2 at t = 0.000918978540772 s (918.978540772 microseconds).
+Shock will pass transducer sd3 at t = 0.000995178540772 s (995.178540772 microseconds).
+Shock will reach end of the shock tube at t = 0.00111781802575 s (1117.81802575 microseconds).
+------------------------------------------------------------
+Shock will pass transducer st1 at t = 0.00123879242575 s (1238.79242575 microseconds).
+Shock will pass transducer st2 at t = 0.00131542442575 s (1315.42442575 microseconds).
+Shock will pass transducer st3 at t = 0.00139190762575 s (1391.90762575 microseconds).
+Shock will pass transducer at1 at t = 0.00156704522575 s (1567.04522575 microseconds).
+Shock will pass transducer at2 at t = 0.00159353162575 s (1593.53162575 microseconds).
+Shock will pass transducer at3 at t = 0.00162046442575 s (1620.46442575 microseconds).
+Shock will pass transducer at4 at t = 0.00173861162575 s (1738.61162575 microseconds).
+Shock will pass transducer at5 at t = 0.00177670442575 s (1776.70442575 microseconds).
+Shock will pass transducer at6 at t = 0.00181390442575 s (1813.90442575 microseconds).
+Shock will pass transducer at7 at t = 0.00182298122575 s (1822.98122575 microseconds).
+Shock will pass transducer at8 at t = 0.00182893322575 s (1828.93322575 microseconds).
+Shock will reach end of the acceleration tube at t = 0.00188666762575 s (1886.66762575 microseconds).

data/referance-data/x2s4111/databox/x2s4111eq_experimental.csv

@@ -0,0 +1,28 @@
+Pitot Version,15-Aug-2018,expansion-tube mode
+Solver,equilibrium.
+Test,'experiment'
+Facility,x2.
+Test gas (state 1),gasgiant_h285ne,gamma,1.6078,R,476.310512747
+Accelerator gas (state 5),Air.
+Driver gas,{'Ar': 0.2, 'He': 0.8}.
+Vs1,3057.74 m/s,Ms1,6.40,Vs2,6720.43 m/s,Ms2,19.42
+state ,P          ,T        ,a     ,V        ,M     ,rho      ,pitot   ,stgn     
+      ,Pa         ,K        ,m/s   ,m/s      ,      ,kg/m^3   ,kPa      ,MPa      
+s4    ,1.550e+07  ,1993.0   ,1571  ,0.0      ,0.00  ,10.46800 ,15500.0 ,15.500   
+s3s   ,7.551e+06  ,1495.3   ,1361  ,1359.7   ,1.00  ,6.79690  ,15500.0 ,15.471   
+s1    ,2000.0     ,298.1    ,478   ,0.0      ,0.00  ,0.01408  ,2.0     ,0.002    
+s2    ,107510.0   ,3080.3   ,1377  ,2450.3   ,1.78  ,0.07089  ,457.6   ,0.568    
+s3    ,1.592e+06  ,802.6    ,997   ,2450.3   ,2.46  ,2.67010  ,14859.4 ,25.078   
+s5    ,1.0        ,298.1    ,346   ,0.0      ,0.00  ,1.17e-05 ,0.0     ,0.000    
+s6    ,499.0295   ,5301.3   ,1633  ,6342.5   ,3.88  ,2.08e-04 ,8.4     ,0.172    
+s7    ,125.9287   ,322.7    ,497   ,6720.4   ,13.52 ,8.19e-04 ,33.0    ,22.623   
+s8    ,7.585137   ,110.3    ,292   ,6759.8   ,23.11 ,1.44e-04 ,5.9     ,22.591   
+s10f  ,5018.354   ,17557.7  ,3691  ,5133.2   ,1.39  ,6.00e-04 ,13.7    ,0.010    
+s10e  ,5619.0     ,10132.5  ,2665  ,5748.8   ,2.16  ,9.65e-04 ,33.2    ,0.056    
+s10cf ,489.6711   ,1714.8   ,1153  ,6464.5   ,5.60  ,6.00e-04 ,22.4    ,0.237    
+s10c  ,488.5144   ,1650.2   ,1102  ,6466.8   ,5.87  ,6.21e-04 ,23.5    ,1.720    
+Ht,22.613, MJ/kg.
+h,-0.2343, MJ/kg.
+Tt, 13441, K.
+Ue,6725, m/s.
+Basic test time,61.40, microseconds

data/referance-data/x2s4111/databox/x2s4111eq_experimental.txt

@@ -0,0 +1,38 @@
+Pitot Version: 15-Aug-2018 doing an expansion tube calculation
+state 4 is the driver condition.
+state 1 is shock tube fill. state 5 is acceleration tube fill.
+state 2 is the shocked test gas.
+state 7 is expanded test gas entering the nozzle.
+state 8 is test gas exiting the nozzle (using area ratio of 5.64).
+state 10f is frozen shocked test gas flowing over the model.
+state 10e is equilibrium shocked test gas flowing over the model.
+state 10c is conditions over 15.0 degree conehead in the test section.
+Solver used is equilibrium.
+Test is 'experiment'
+Facility is x2.
+Driver gas is {'Ar': 0.2, 'He': 0.8} (by moles).
+Test gas (state 1) is gasgiant_h285ne (gamma = 1.6078, R = 476.310512747, {'H2': 0.15, 'Ne': 0.85} by moles).
+Accelerator gas (state 5) is Air.
+Vs1 = 3057.74 m/s, Ms1 = 6.40, Vs2 = 6720.43 m/s, Ms2 = 19.42
+state P          T        a     V        M     rho      pitot   stgn     
+      Pa         K        m/s   m/s            kg/m^3   kPa     MPa      
+s4    1.550e+07  1993.0   1571  0.0      0.00  10.46800 15500.0 15.500   
+s3s   7.551e+06  1495.3   1361  1359.7   1.00  6.79690  15500.0 15.471   
+s1    2000.0     298.1    478   0.0      0.00  0.01408  2.0     0.002    
+s2    107510.0   3080.3   1377  2450.3   1.78  0.07089  457.6   0.568    
+s3    1.592e+06  802.6    997   2450.3   2.46  2.67010  14859.4 25.078   
+s5    1.0        298.1    346   0.0      0.00  1.17e-05 0.0     0.000    
+s6    499.0295   5301.3   1633  6342.5   3.88  2.08e-04 8.4     0.172    
+s7    125.9287   322.7    497   6720.4   13.52 8.19e-04 33.0    22.623   
+s8    7.585137   110.3    292   6759.8   23.11 1.44e-04 5.9     22.591   
+s10f  5018.354   17557.7  3691  5133.2   1.39  6.00e-04 13.7    0.010    
+s10e  5619.0     10132.5  2665  5748.8   2.16  9.65e-04 33.2    0.056    
+s10cf 489.6711   1714.8   1153  6464.5   5.60  6.00e-04 22.4    0.237    
+s10c  488.5144   1650.2   1102  6466.8   5.87  6.21e-04 23.5    1.720    
+The total enthalpy (Ht) leaving the nozzle is 22.613 MJ/kg (H8 - h1).
+The total temperature (Tt) leaving the nozzle is 13441 K.
+The freestream enthalpy (h) leaving the nozzle is -0.2343 MJ/kg (h8 - h1).
+The flight equivalent velocity (Ue) is 6725 m/s.
+Basic test time = 61.40 microseconds
+Species in the shock layer at equilibrium (s10e) (by moles):
+{'e-': 0.046354, 'H': 0.20263, 'Ne+': 0.00020291, 'Ne': 0.70467, 'H+': 0.046151}

data/referance-data/x2s4111/databox/x2s4111pg_theoretical.csv

@@ -0,0 +1,27 @@
+Pitot Version,15-Aug-2018,expansion-tube mode
+Solver,perfect gas
+Test,'fulltheory-pressure'
+Facility,x2.
+Test gas (state 1),gasgiant_h285ne,gamma,1.6078,R,476.310512747
+Accelerator gas (state 5),Air.
+Driver gas,He:0.80,Ar:0.20.
+Vs1,4521.09 m/s,Ms1,9.46,Vs2,10323.69 m/s,Ms2,29.83
+state ,P          ,T        ,a     ,V        ,M     ,rho      ,pitot   ,stgn     
+      ,Pa         ,K        ,m/s   ,m/s      ,      ,kg/m^3   ,kPa      ,MPa      
+s4    ,1.550e+07  ,1993.0   ,1571  ,0.0      ,0.00  ,10.46941 ,15500.0 ,15.500   
+s3s   ,7.551e+06  ,1494.7   ,1360  ,1360.4   ,1.00  ,6.80010  ,15500.0 ,15.500   
+s1    ,2000.0     ,298.1    ,478   ,0.0      ,0.00  ,0.01408  ,2.0     ,0.002    
+s2    ,220307.2   ,7935.9   ,2465  ,3428.6   ,1.39  ,0.05828  ,721.4   ,0.748    
+s3    ,216358.6   ,360.9    ,668   ,3428.6   ,5.13  ,0.80696  ,8456.5  ,64.527   
+s5    ,1.0        ,298.1    ,346   ,0.0      ,0.00  ,1.17e-05 ,0.0     ,0.000    
+s6    ,1037.642   ,51852.0  ,4565  ,8593.4   ,1.88  ,6.97e-05 ,5.2     ,0.007    
+s7    ,12.92844   ,199.5    ,391   ,10227.2  ,26.16 ,1.36e-04 ,12.7    ,17.751   
+s8    ,0.7990541  ,69.7     ,231   ,10243.1  ,44.35 ,2.41e-05 ,2.2     ,17.750   
+s10f  ,1937.917   ,39436.9  ,5496  ,7851.8   ,1.43  ,1.03e-04 ,6.6     ,0.007    
+s10e  ,2287.0     ,12973.2  ,3227  ,9268.3   ,2.87  ,2.53e-04 ,22.1    ,0.093    
+s10c  ,184.495    ,3615.9   ,1664  ,9802.9   ,5.89  ,1.07e-04 ,9.2     ,0.119    
+Ht,52.172, MJ/kg.
+h,-0.2879, MJ/kg.
+Tt, 41705, K.
+Ue,10215, m/s.
+Basic test time,20.08, microseconds

data/referance-data/x2s4111/databox/x2s4111pg_theoretical.txt

@@ -0,0 +1,37 @@
+Pitot Version: 15-Aug-2018 doing an expansion tube calculation
+state 4 is the driver condition.
+state 1 is shock tube fill. state 5 is acceleration tube fill.
+state 2 is the shocked test gas.
+state 7 is expanded test gas entering the nozzle.
+state 8 is test gas exiting the nozzle (using area ratio of 5.64).
+state 10f is frozen shocked test gas flowing over the model.
+state 10e is equilibrium shocked test gas flowing over the model.
+state 10c is conditions over 15.0 degree conehead in the test section.
+Solver used is perfect gas.
+Test is 'fulltheory-pressure'
+Facility is x2.
+Driver gas is He:0.80,Ar:0.20.
+Test gas (state 1) is gasgiant_h285ne (gamma = 1.6078, R = 476.310512747).
+Accelerator gas (state 5) is Air.
+Vs1 = 4521.09 m/s, Ms1 = 9.46, Vs2 = 10323.69 m/s, Ms2 = 29.83
+state P          T        a     V        M     rho      pitot   stgn     
+      Pa         K        m/s   m/s            kg/m^3   kPa     MPa      
+s4    1.550e+07  1993.0   1571  0.0      0.00  10.46941 15500.0 15.500   
+s3s   7.551e+06  1494.7   1360  1360.4   1.00  6.80010  15500.0 15.500   
+s1    2000.0     298.1    478   0.0      0.00  0.01408  2.0     0.002    
+s2    220307.2   7935.9   2465  3428.6   1.39  0.05828  721.4   0.748    
+s3    216358.6   360.9    668   3428.6   5.13  0.80696  8456.5  64.527   
+s5    1.0        298.1    346   0.0      0.00  1.17e-05 0.0     0.000    
+s6    1037.642   51852.0  4565  8593.4   1.88  6.97e-05 5.2     0.007    
+s7    12.92844   199.5    391   10227.2  26.16 1.36e-04 12.7    17.751   
+s8    0.7990541  69.7     231   10243.1  44.35 2.41e-05 2.2     17.750   
+s10f  1937.917   39436.9  5496  7851.8   1.43  1.03e-04 6.6     0.007    
+s10e  2287.0     12973.2  3227  9268.3   2.87  2.53e-04 22.1    0.093    
+s10c  184.495    3615.9   1664  9802.9   5.89  1.07e-04 9.2     0.119    
+The total enthalpy (Ht) leaving the nozzle is 52.172 MJ/kg (H8 - h1).
+The total temperature (Tt) leaving the nozzle is 41705 K.
+The freestream enthalpy (h) leaving the nozzle is -0.2879 MJ/kg (h8 - h1).
+The flight equivalent velocity (Ue) is 10215 m/s.
+Basic test time = 20.08 microseconds
+Species in the shock layer at equilibrium (s10e) (by moles):
+{'e-': 0.21377, 'H': 0.026268, 'Ne+': 0.03494, 'Ne': 0.54618, 'H+': 0.17883}