Add data from x2s4111 and throw in some primitive graphing
This commit is contained in:
Cal Wing
parent
76012f919f
commit
395b36c361
BIN
data/x2s4111/databox/Thumbs.db
Normal file
BIN
data/x2s4111/databox/Thumbs.db
Normal file
Binary file not shown.
BIN
data/x2s4111/databox/experimental-x-t-diagram.pdf
Normal file
BIN
data/x2s4111/databox/experimental-x-t-diagram.pdf
Normal file
Binary file not shown.
BIN
data/x2s4111/databox/experimental-x-t-diagram.png
Normal file
BIN
data/x2s4111/databox/experimental-x-t-diagram.png
Normal file
Binary file not shown.
|
After Width: | Height: | Size: 114 KiB |
BIN
data/x2s4111/databox/shock-speed-plot.pdf
Normal file
BIN
data/x2s4111/databox/shock-speed-plot.pdf
Normal file
Binary file not shown.
BIN
data/x2s4111/databox/shock-speed-plot.png
Normal file
BIN
data/x2s4111/databox/shock-speed-plot.png
Normal file
Binary file not shown.
|
After Width: | Height: | Size: 140 KiB |
BIN
data/x2s4111/databox/x-t-diagram.pdf
Normal file
BIN
data/x2s4111/databox/x-t-diagram.pdf
Normal file
Binary file not shown.
BIN
data/x2s4111/databox/x-t-diagram.png
Normal file
BIN
data/x2s4111/databox/x-t-diagram.png
Normal file
Binary file not shown.
|
After Width: | Height: | Size: 121 KiB |
15024
data/x2s4111/databox/x2s4111.lvm
Normal file
15024
data/x2s4111/databox/x2s4111.lvm
Normal file
File diff suppressed because it is too large
Load Diff
5
data/x2s4111/databox/x2s4111.sh
Normal file
5
data/x2s4111/databox/x2s4111.sh
Normal file
@ -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
|
||||
37
data/x2s4111/databox/x2s4111.shotcfg
Normal file
37
data/x2s4111/databox/x2s4111.shotcfg
Normal file
@ -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}
|
||||
BIN
data/x2s4111/databox/x2s4111.tdms
Normal file
BIN
data/x2s4111/databox/x2s4111.tdms
Normal file
Binary file not shown.
BIN
data/x2s4111/databox/x2s4111.tdms_index
Normal file
BIN
data/x2s4111/databox/x2s4111.tdms_index
Normal file
Binary file not shown.
86
data/x2s4111/databox/x2s4111.txt
Normal file
86
data/x2s4111/databox/x2s4111.txt
Normal file
@ -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
|
||||
|
||||
|
||||
BIN
data/x2s4111/databox/x2s4111.xcfg
Normal file
BIN
data/x2s4111/databox/x2s4111.xcfg
Normal file
Binary file not shown.
40
data/x2s4111/databox/x2s4111.xconfig
Normal file
40
data/x2s4111/databox/x2s4111.xconfig
Normal file
@ -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
|
||||
78
data/x2s4111/databox/x2s4111_basic_data_summary.txt
Normal file
78
data/x2s4111/databox/x2s4111_basic_data_summary.txt
Normal file
@ -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
|
||||
------------------------------------------------------------
|
||||
8137
data/x2s4111/databox/x2s4111_full_data_summary.shot_obj
Normal file
8137
data/x2s4111/databox/x2s4111_full_data_summary.shot_obj
Normal file
File diff suppressed because one or more lines are too long
116
data/x2s4111/databox/x2s4111_full_data_summary.txt
Normal file
116
data/x2s4111/databox/x2s4111_full_data_summary.txt
Normal file
@ -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
|
||||
------------------------------------------------------------
|
||||
@ -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).
|
||||
28
data/x2s4111/databox/x2s4111eq_experimental.csv
Normal file
28
data/x2s4111/databox/x2s4111eq_experimental.csv
Normal file
@ -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
|
||||
|
38
data/x2s4111/databox/x2s4111eq_experimental.txt
Normal file
38
data/x2s4111/databox/x2s4111eq_experimental.txt
Normal file
@ -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}
|
||||
27
data/x2s4111/databox/x2s4111pg_theoretical.csv
Normal file
27
data/x2s4111/databox/x2s4111pg_theoretical.csv
Normal file
@ -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
|
||||
|
37
data/x2s4111/databox/x2s4111pg_theoretical.txt
Normal file
37
data/x2s4111/databox/x2s4111pg_theoretical.txt
Normal file
@ -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}
|
||||
35
main.py
Normal file
35
main.py
Normal file
@ -0,0 +1,35 @@
|
||||
# Cal Wing
|
||||
# Sep 2024
|
||||
|
||||
from nptdms import TdmsFile
|
||||
|
||||
from makeGraph import makeGraph, pltKeyClose, UQ_COLOURS as UQC
|
||||
|
||||
|
||||
DATA_FILE = './data/x2s4111/databox/x2s4111.tdms'
|
||||
|
||||
|
||||
def main():
|
||||
data = TdmsFile.read(DATA_FILE)
|
||||
|
||||
channels = data.groups()[0].channels()
|
||||
|
||||
graphData = {
|
||||
"plots": []
|
||||
}
|
||||
|
||||
for channel in channels[1:15]:
|
||||
graphData["plots"].append({
|
||||
"x": channels[0][:],
|
||||
"y": channel[:],
|
||||
"label": channel.name
|
||||
})
|
||||
|
||||
makeGraph(graphData, showPlot=True, doProgramBlock=False)
|
||||
|
||||
|
||||
|
||||
if __name__ == '__main__':
|
||||
main()
|
||||
|
||||
pltKeyClose()
|
||||
Loading…
x
Reference in New Issue
Block a user