% AUTHOR: Courtney Hiatt
% DATE: 3/26/24
% This code is the driver for the material properties function and outputs
% mass, volume, and cost requirements for different materials. 


clc
clear all

LowestTemp = -1; %deg C
MinEnergyRequirement = 84153; %kJ 
MaxEnergyRequirement = 119593; %kJ

%Properties of water
Water = 'Water';
WaterHeatofFusion = 334; %kJ/kg
WaterSpecificHeat = 4.187; %kJ/kgC
TempofFusion = 0; %deg C
WaterDensityMatrix = [-50 -40 -35 -30 -25 -20 -15 -10 -5 0 1 4 10 15 20 25 30 35 40 45 50 55 60 65 70; 921.6 920.8 920.4 920 919.6 919.4 919.4 918.9 917.5 916.2 999.90 999.97 999.70 999.10 998.21 997.05 995.65 994.03 992.22 990.21 998.04 985.69 983.21 980.55 977.76]; %kg/m3
WaterCost = 0.0002189;
MaterialProperties(Water, WaterHeatofFusion, TempofFusion, WaterSpecificHeat, WaterDensityMatrix, MinEnergyRequirement, MaxEnergyRequirement, LowestTemp, WaterCost)


%Properties of concrete
Concrete = 'Concrete';
ConcreteHeatofFusion = 0; %kJ/kg
ConcreteSpecificHeat = 1; %kJ/kgC
ConcreteDensityMatrix = [-100,0, 80, 95, 180; 2300, 2300, 2300, 2300, 2254];
ConcreteTempofFusion = 1200; %deg C
ConcreteCost = 0.10;
MaterialProperties(Concrete, ConcreteHeatofFusion, ConcreteTempofFusion, ConcreteSpecificHeat, ConcreteDensityMatrix, MinEnergyRequirement, MaxEnergyRequirement, LowestTemp, ConcreteCost)