Simulation of Propellant Loading using GFSSP - National Council of ...

LiquidVolumeInput Data for Integrated Ground System, LH 2 Tank and Flare Stack Modelof Propellant LoadingLH2 Storage Tank PressureAmbient TemperatureLH2 Propellant Load46.3 psia85 ° F48593 lbPre-Chill Valve C v 16Slow Fill & Topping Valve C v 12Fast Fill Valve C v 140Replenish Valve C v 5.64Vent Valve Area 20.94 in 2Vent Valve C d 0.552Ground System Pipe Length and Volume 1910 ft / 879 ft 3Flare Stack Pipe Length and Volume 1305 ft /1605 ft 3Tank Volume 11,620 ft 3Ground System Pipe Mass29314 lbTank Mass8742 lbFoam MassMetal (Al-Li) thicknessFoam (BX-265) thickness (Tank Barrel)Foam (BX-265) thickness (Dome)673 lb0.1934 in1 in0.5 inCommon Bulkhead Conductance 0.045 Btu/hr-ft 2 -F6

Summary Result for LH2 LoadingDesign ParametersGround Supply VaporGround Supply Vapor QualityQuality = 50%= 2%Pre-chill Time (after start) 129 Minutes 129 Minutes5% Tank Fill Time (after pre-chill) 23 Minutes 23 Minutes95% Tank Fill Time (after pre-chill) 73 Minutes 73 Minutes100% Tank Fill Time (after pre-chill) 87 Minutes 87 MinutesTank Chill-down Time (after start) 194 Minutes 194 MinutesMaximum Tank Pressure (pre-chill) 15.94 psia 15.94 psiaMaximum Ullage Pressure (Replenish) 15.5 psia 14.85 psiaMaximum Vent Flowrate 0.95 lb/sec 0.67 lb/secAmount of GH2 Vented 4069 lb 3681 lbMinimum Foam Surface Temperature 6.5 F 6.2 F7

Time Step Study and Effect of CPU TIMEDerived Parametersfrom SimulationBase Run (Time step= 0.1 s)Run1 (Time step =0.05 s)Run 2(Time step = 0.2s)Pre-chilled Time (min) 128 138 117Time taken for 5% Fill(min)Time taken for 98%Fill (min)23.2 22.8 23.472.7 72.5 72.9Time taken for 100%86.88 86.66 87.0Fill (min)Computational Time11 hours 26 mins 23 hours 16 mins 6 hours 17 minsneeded for Simulation

Tank Inflow rate and Vent flow rate16F2627 LBM/SEC Vent Valve F910 LBM/SEC Pipe 910WinPlot v4.60 rc1Inflow128Fastfill4SlowfillPrechillTopping0Vent FlowReplenish-40 2000 4000 6000 8000 10000 12000 14000TIME SECONDS9:49:43AM 03/02/20099

Pressure in LH2 Tank18P10 PSIA Node 10P26 PSIA Node 26P22 PSIAP28 PSIA Node 28WinPlot v4.60 rc117Common Bulkhead16End of PrechillMid Barrel15UllageAtmosphere140 4000 8000 12000 16000TIME SECONDS3:26:58PM 03/02/200910

Propellant Temperature and Quality in LH2 TankHydrogen TemperatureQuality (Vapor Fraction)200T10 DEG_F Node T22 DEG_F T26 DEG_F Node 2XV10 - Node 10 XV22 - XV26 - Node 26WinPlot v4.60 rc1 2.0WinPlot v4.60 rc11001.50-100-200-300-400Common BulkheadMid BarrelUllage10 1.00.50.0-0.5Condensation BeginsCommon Bulkhead Mid BarrelCondensation EndsUllage-5000 2000 4000 6000 8000 10000 12000 14000TIME SECONDS 2:17:40PM 03/03/2009-1.00 2000 4000 6000 8000 10000 12000 14000TIME SECONDS 2:25:10PM 03/03/200911

Pressure and Facility Temperature in the Ground System60P1 PSIA Node 1 P8 PSIA Node 8 P9 PSIA Node 9WinPlot v4.60 rc1100DEG_F S Node 11 DEG_F S Node 15 DEG_F S Node 17WinPlot v4.60 rc150040Storage Tank-100 Cross Country LinePdSl Pad Slope30Valve Upstream20Valve Downstream100 2000 4000 6000 8000 10000 12000 14000TIME SECONDS2:56:57PM 03/03/2009-200-300Storage Tank-400-5000 2000 4000 6000 8000 10000 12000 14000TIME SECONDS3:02:10PM 03/03/200912

Ullage Pressure & Flowrate in Vent Line2.0LBM/SEC Vent ValveLBM/SEC Vent ValveWinPlot v4.60 rc11.5PSIDLBM/SECVent Flow (50% Quality)Gauge Pressure (50% quality)10 1.00.50.0-0.5Gauge Pressure (2% Quality)Vent tFl Flow (2%Quality)-1.00 2000 4000 6000 8000 10000 12000 14000TIME SECONDS9:17:05AM 03/04/200913

CONCLUSION• The Simulation Model will be carried on at Alabama A & M University by thestudents in the senior design students for further optimization and parametricstudy.• The model simulates all phases of loading: Pre-chill, Slow fill, Fast fill, Toppingand Replenish• A new solver technique developed at MSFC will be implemented and tested in thesenior design course.• The computation time is cut by 50% by optimizing the time step. However, furtherstudy is needed to optimize further.• STUDENTS WORK• Improve computational efficiency of numerical simulation by introducing fastsolver and using fast compiler• PARAMETRIC STUDY• TEST WITH BROYDEN SOLVER (replacing the existing Newton Raphson).14

REFERENCES and ACKNOWLEDGEMENTREFERENCES:1. Majumdar, Aloke and LeClair Andre, “LO2 and LH2 Tank Loading and Boil-Off AnalysisReport”, report no. MPS-PR-08, NASA MSFC, March 21, 2008Acknowledgments:Thanks to My mentor Dr. Alok Majumdar, Ms. Melissa Van Dyke, Thermal AnalysisBranch Chief, Dr. Andre Le’Clair, and Mr. Rick Moore, Marshall Space Center,Huntsville, Alabama.Special Thanks to ESMD Space Grant Faculty Fellowship Program for all the supporttowards this work