R Notebook

bind_rows(data.frame(Component="Electrolyte Cost $/Ah",
          High=CostAh,
          Medium=CostAh*2/3,
          Low=CostAh/3),
          componentcosts)

Costs<-findCost(0.738,150) #SOC Range, Flow rate

pal2 <- brewer.pal(5,"Blues")%>%rev

powerComponentColours<-data.frame(Component=c("Bipolar Plate",
                                               "PVC Frame",
                                               "Membrane",
                                               "Felt",
                                               "Al End Plate"),colour=pal2)
othercolours<-data.frame(Component=c("PCS",
                                     "Pumps",
                                     "Labor",
                                     "Chemicals"),
                         colour=c("Gold2",
                                  "Green3",
                                  "Red4",
                                  "Red3"))

Costs%>%
  melt(id.var="Component")%>%
  mutate(Scenario2=ifelse(variable=="High","Low","Medium"),
         Scenario2=ifelse(variable=="Low","High",Scenario2),
         Scenario2=factor(Scenario2,levels=c("Low","Medium","High")))%>%
  mutate(Component=factor(Component,levels=c("Bipolar Plate",
                                             "PVC Frame",
                                             "Nf 211",
                                             "Carbon Felt",
                                             "Al Endplate",
                                             "PCS",
                                             "Pumps",
                                             "Labor",
                                             "Chemicals")))%>%
  ggplot(aes(x=Scenario2,y=value,fill=Component))+
  geom_bar(stat="identity", colour="black",size=1)+theme_bw()+geom_hline(yintercept = 250,size=1.2)+ylab("Total Cost $/kWh")+xlab("Volume")+
  theme(strip.background = element_blank(),
        axis.title = element_text(face="bold",size=16),
        axis.text = element_text(face="bold",size=10,colour="black"),
        strip.text =  element_text(face="bold",size=16,colour="black"),
        legend.text =   element_text(face="bold",size=14,colour="black"),
        legend.title =  element_text(face="bold",size=14,colour="black"),
        legend.position="bottom")+
    guides(fill = guide_legend(nrow = 3))+
  scale_y_continuous(expand = expand_scale(mult = c(0,0.05), add = 0), limits = c(0, NA))+
  scale_fill_manual(values = bind_rows(powerComponentColours,othercolours)$colour,name="")

million pie charts

DHPS earlier gen

"C:\\Users\\craw038\\OneDrive - PNNL\\Raw Data - DHPS - 200 cm2.xlsx"%>%
            read_excel(skip = 6,sheet=3)%>%
            filter(Step_Index==4,
                   Cycle_Index==3)%>%
            transmute(Voltage=`Voltage(V)`,
                      DisCap=`Discharge_Capacity(Ah)`,
                      Current=round(-`Current(A)`/200*1000),
                      Flow=400)%>%
  mutate(SOC=DisCap/Cap,
         SOC=(1-max(SOC))/2+SOC,
         SOC=1-SOC)->data3

Costs<-Costs%>%mutate(Type="Mixed")%>%bind_rows(findCost(0.73,150)%>%mutate(Type="Unmixed")) #SOC Range, Flow rate


Costs%>%
  melt(id.var=c("Component","Type"))%>%
  mutate(Scenario2=ifelse(variable=="High","Low","Medium"),
         Scenario2=ifelse(variable=="Low","High",Scenario2),
         Scenario2=factor(Scenario2,levels=c("Low","Medium","High")))%>%
  mutate(Component=factor(Component,levels=c("Bipolar Plate",
                                             "PVC Frame",
                                             "Nf 211",
                                             "Carbon Felt",
                                             "Al Endplate",
                                             "PCS",
                                             "Pumps",
                                             "Labor",
                                             "Chemicals")))%>%
  ggplot(aes(x=Scenario2,y=value,fill=Component))+
  geom_bar(stat="identity", colour="black",size=1)+theme_bw()+geom_hline(yintercept = 250,size=1.2)+ylab("Total Cost $/kWh")+xlab("Volume")+
  facet_wrap(.~Type)+
  theme(strip.background = element_blank(),
        axis.title = element_text(face="bold",size=16),
        axis.text = element_text(face="bold",size=10,colour="black"),
        strip.text =  element_text(face="bold",size=16,colour="black"),
        legend.text =   element_text(face="bold",size=14,colour="black"),
        legend.title =  element_text(face="bold",size=14,colour="black"),
        legend.position="bottom")+
    guides(fill = guide_legend(nrow = 3))+
  scale_y_continuous(expand = expand_scale(mult = c(0,0.05), add = 0), limits = c(0, NA))+
  scale_fill_manual(values = bind_rows(powerComponentColours,othercolours)$colour,name="")