changed rohacell31 to rohacell41S to agree with Brewer

example/cryo_input.toml

@@ -56,7 +56,7 @@ name = "TASOPT Model with cryo fuel and HX"
     tank_placement = "front" #Location of tank in fuselage. Choices are "front", "rear", or "both".
     insulation_segment_base_thickness = [0.01,0.01,0.01] #m
     insulation_thicknesses_design_indices = [1,2,3]
-    insulation_material = ["rohacell31","rohacell31","rohacell31"] #Options: "rohacell31", "polyurethane"
+    insulation_material = ["rohacell41S","rohacell41S","rohacell41S"] #Options: "rohacell41S", "polyurethane"
     skin_yield_strength = "172e6 Pa"
     tank_skin_density = 2825 #kg/m^3
     maximum_boiloff_rate = 0.15 #%/hour

src/structures/fuseW.jl

@@ -210,7 +210,7 @@ function fusew(Nland,Wfix,Wpay,Wpadd,Wseat,Wapu,Weng,
                   xcabin = 0.5 * (xshell1 + xshell2 - (ltank + 2.0*ft_to_m))
             end
       else
-            xcabin = 0.5 * (xshell1 + xshell2)
+            xcabin = 0.5 * (xshell1 + xshell2) #two or zero tanks
       end
       lcabin = xshell2 - xshell1 - nftanks * (ltank + 2.0*ft_to_m) #cabin length is smaller if there are fuel tanks
 

src/structures/tankWmech.jl

@@ -160,8 +160,8 @@ This function calculates the density of different insulation materials.
       - `ρ::Float64`: mass density (kg/m^3).
 """
 function insulation_density_calc(material)
-      if material == "rohacell31"
-            ρ = 32.0 #kg/m^3. From manufacturer sheet
+      if material == "rohacell41S"
+            ρ = 35.0 #kg/m^3. From Brewer (1991)
       elseif material == "polyurethane"
             ρ = 27.0 #kg/m^3. From Brewer (1991)
       end

src/structures/tankWthermal.jl

@@ -198,7 +198,8 @@ This function calculates the thermal conductivity of different insulation materi
       - `k::Float64`: thermal conductivity (W/(m K)).
 """
 function insulation_conductivity_calc(T, material)
-      if material == "rohacell31"
+      if material == "rohacell41S"
+            #Note: Brewer (1991) only had one point for rohacell thermal conductivity. They assumed the same curve as PVC
             k = 0.001579 + 1.283e-4 * T - 3.353e-7*T^2 + 8.487e-10 * T^3 # W/(m K), polynomial fit to Fig. 4.78 in Brewer (1991) between 20 and 320 K
       elseif material == "polyurethane"
             k = 2.114e-13 * T^5 - 1.639e-10 *T^4 + 4.438e-8 * T^3 - 5.222e-6*T^2 + 3.740e-4*T - 2.192e-3

test/unit_test_fueltank.jl

@@ -6,7 +6,7 @@ dRfuse = 0.3
 hconvgas = 0.0
 Tfuel = 20.0
 t_cond = [0.15,0.15,0.15]
-material_insul = ["rohacell31", "rohacell31", "polyurethane"]
+material_insul = ["rohacell41S", "rohacell41S", "polyurethane"]
 iinsuldes = [1,2,3]
 time_flight = 7*3600.0
 fstring = 0.1
@@ -34,7 +34,7 @@ m_boiloff = threshold_percent * Wfuel /(100 * gee) *time_flight/3600
                         wfb, nfweb, sigskin, material_insul, rhoskin, Wfuel, threshold_percent, clearance_fuse, AR, 
                         iinsuldes, ifuel, qfac)
 
-    outputs_size_check = (162741.4353849525, 0.49077274862420084, 12.790944804923033, 0.004247366632687734, 145.623998835008, 100000.0, 107.03363914373088, 0.004916887075394796, 0.0049105639181476314, 1.9092272513757993, 2141.366879896994, 21900.169398210935, 31321.95159534658, [10017.291267982546, 11061.33253815481, 10243.327789209226], 14.695255169223437, 62741.43538495249)
+    outputs_size_check = (164717.55636677786, 0.49077274862420084, 12.790944804923033, 0.004247366632687734, 145.623998835008, 100000.0, 107.03363914373088, 0.004916887075394796, 0.0049105639181476314, 1.9092272513757993, 2141.366879896994, 21900.169398210935, 33298.07257717196, [10956.41232435591, 12098.332463606823, 10243.327789209226], 14.695255169223437, 64717.55636677787)
     for i in 1:length(outputs_size)
         @test outputs_size[i] ≈ outputs_size_check[i]
     end
@@ -45,7 +45,7 @@ m_boiloff = threshold_percent * Wfuel /(100 * gee) *time_flight/3600
                   sigskin, material_insul, rhoskin,
                   Wfuel, m_boiloff, t_cond, clearance_fuse, AR)
 
-    outputs_mech_check = (162741.4353849525, 12.790944804923033, 0.004916887075394796, 1.9092272513757993, 145.623998835008, 62741.43538495249, 100000.0, 31321.95159534658, 0.0049105639181476314, 2141.366879896994, 21900.169398210935, [10017.291267982546, 11061.33253815481, 10243.327789209226], [15.832773779668218, 19.254692964997925, 23.016493195949316, 27.117221510580325], 14.695255169223437)
+    outputs_mech_check = (164717.55636677786, 12.790944804923033, 0.004916887075394796, 1.9092272513757993, 145.623998835008, 64717.55636677787, 100000.0, 33298.07257717196, 0.0049105639181476314, 2141.366879896994, 21900.169398210935, [10956.41232435591, 12098.332463606823, 10243.327789209226], [15.832773779668218, 19.254692964997925, 23.016493195949316, 27.117221510580325], 14.695255169223437)
     for i in 1:length(outputs_mech)
         @test outputs_mech[i] ≈ outputs_mech_check[i]
     end