(* ::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:: *)
celsiusToKelvin=Quantity[0.,"DegreesCelsius"]//UnitConvert//QuantityMagnitude; 
\[ScriptK]b=UnitConvert[Quantity["BoltzmannConstant"],"Joules"/"Kelvins"];
kb=\[ScriptK]b//UnitConvert//QuantityMagnitude;


(* ::Input::Initialization:: *)
\[ScriptM]Li=6.941 Quantity["AtomicMassUnit"];
mLi=\[ScriptM]Li//UnitConvert//QuantityMagnitude; (* in kg *)


(* ::Input::Initialization:: *)
\[ScriptL]Li =Quantity[147 10^3,"Joules"/"Moles"]/Quantity[6.9 10^-3,"Kilograms"/"Moles"];
lLi=\[ScriptL]Li//QuantityMagnitude;


(* ::Input::Initialization:: *)
pLi[TK_]:=Exp[26.89 - 18880/TK-0.4942 Log[TK]];


(* ::Input::Initialization:: *)
\[ScriptP]Li[TK_]:=Quantity[pLi[TK],"Pascals"];
\[ScriptP]Li[T\[DoubleStruckCapitalK]_Quantity]:=Quantity[pLi[QuantityMagnitude[T\[DoubleStruckCapitalK]]],"Pascals"];

\[ScriptN]eqLi[TC_]:=UnitConvert[\[ScriptP]Li[TC+celsiusToKelvin] /(\[ScriptK]b Quantity[(TC+celsiusToKelvin),"Kelvins"]),"Meters" ^-3];
\[ScriptN]eqLi[T\[DoubleStruckCapitalK]_Quantity]:=UnitConvert[\[ScriptP]Li[T\[DoubleStruckCapitalK]] /(\[ScriptK]b T\[DoubleStruckCapitalK])];

neqLi[TC_]:=\[ScriptN]eqLi[TC]//UnitConvert//QuantityMagnitude;


(* ::Input::Initialization:: *)
\[CapitalGamma]Li[T\[DoubleStruckCapitalK]_Quantity]:=\[ScriptN]eqLi[T\[DoubleStruckCapitalK]] Sqrt[\[ScriptK]b T\[DoubleStruckCapitalK]/(2\[Pi] \[ScriptM]Li)]//UnitConvert;
\[CapitalGamma]Li[TK_]:=\[CapitalGamma]Li[Quantity[TK,"Kelvins"]]


(* ::Input::Initialization:: *)
(* Latent heat flux through a surface at some density and temperature *)
\[CapitalGamma]\[ScriptL]Li[\[ScriptN]_,T\[DoubleStruckCapitalK]_]:=\[ScriptL]Li \[ScriptM]Li \[ScriptN] Sqrt[\[ScriptK]b T\[DoubleStruckCapitalK]/(2\[Pi] \[ScriptM]Li)]//UnitConvert



