(* ::Package:: *)

(************************************************************************)
(* This file was generated automatically by the Mathematica front end.  *)
(* It contains Initialization cells from a Notebook file, which         *)
(* typically will have the same name as this file except ending in      *)
(* ".nb" instead of ".m".                                               *)
(*                                                                      *)
(* This file is intended to be loaded into the Mathematica kernel using *)
(* the package loading commands Get or Needs.  Doing so is equivalent   *)
(* to using the Evaluate Initialization Cells menu command in the front *)
(* end.                                                                 *)
(*                                                                      *)
(* DO NOT EDIT THIS FILE.  This entire file is regenerated              *)
(* automatically each time the parent Notebook file is saved in the     *)
(* Mathematica front end.  Any changes you make to this file will be    *)
(* overwritten.                                                         *)
(************************************************************************)



(* ::Input::Initialization:: *)
(* some abbreviations for units *)
gpcc=Quantity["Grams"/("Centimeters")^3];
jpkgk=Quantity[1.0,"Joules"/("Kilograms""Kelvins")];
wpmk=Quantity["Watts"/("Meters""Kelvins")];


(* ::Input::Initialization:: *)
\[Rho]SS=8. gpcc;


(* ::Input::Initialization:: *)
\[Rho]Mullite=2.8 gpcc ;


(* ::Input::Initialization:: *)
\[ScriptC]SS=490 jpkgk;


(* ::Input::Initialization:: *)
(* Temperature is in Kelvin *)
\[ScriptC]MulliteT[T_]:=Quantity[a+10^-5 b T + 10^4 c T^-2+ 10 d T^-0.5+10^6 e T^-3,"Joules"/("Grams""Kelvins")]/.{a->1.58816,b->-1.2254,c->-2.2240,d->-1.1142,e->2.487}//Evaluate;
\[ScriptC]Mullite=\[ScriptC]MulliteT[950];(*worst case scenario for highest specific heat *);
\[ScriptC]MulliteTFactor[TK_]:=\[ScriptC]MulliteT[TK]/\[ScriptC]Mullite


(* ::Input::Initialization:: *)
\[Epsilon]SS=0.3;
\[Epsilon]CeramicHeaterHighEmissivity=0.85;
\[Epsilon]CeramicHeaterUncoatedSurfaces=0.6; (* at 800C *)


(* ::Input::Initialization:: *)
\[Rho]ElTemperatureFactorKanthal[Tc_]:=Interpolation[{{20,1},{250,1.01},{500,1.03},{800,1.05},{1000,1.06},{1200,1.06}},InterpolationOrder->1][Tc];


(* ::Input::Initialization:: *)
materialLib[\[Rho]_,\[ScriptC]_,\[ScriptK]_,\[Epsilon]_]:=<|"\[Rho]"->\[Rho],"\[ScriptC]"->\[ScriptC],"\[ScriptK]"->\[ScriptK],"\[Epsilon]"->\[Epsilon],"\[ScriptS]"->\[Rho] \[ScriptC]|>;
materialLib["SS"]:=materialLib[\[Rho]SS,\[ScriptC]SS,I,\[Epsilon]SS];

materialLib["CeramicHeater"]:=materialLib[\[Rho]Mullite,\[ScriptC]Mullite,I,\[Epsilon]CeramicHeaterHighEmissivity];
materialLib["CeramicHeaterUncoated"]:=materialLib[\[Rho]Mullite,\[ScriptC]Mullite,I,\[Epsilon]CeramicHeaterUncoatedSurfaces];
materialLib["Arbitrary\[Epsilon]0.1"]:=materialLib[I,I,I,0.1];
materialLib["Arbitrary\[Epsilon]0.2"]:=materialLib[I,I,I,0.2];


(* ::Input::Initialization:: *)
sheetGage[{"SS",g_}]:=sheetGageSS[g];

sheetGageSS["16th"]:=Quantity[0.0625,"Inches"];
