Interlayer

library(tidyverse)
library(gganimate)
library(readxl)

interlayer_ <- read_excel(file.path("interlayer.xlsx"))

interlayer <- interlayer_ %>% 
  pivot_longer(cols = c(colnames(interlayer_)[2],colnames(interlayer_)[9:15]),names_to = c("composition"),values_to=c("CTE"))

new <- interlayer %>% filter(composition %in% c("Zr","Cr")) %>%
  rename("comp"="composition")

p <- interlayer %>% 
  filter(!composition %in% c("Zr","Cr")) %>% 
  ggplot(aes(Temp,CTE,colour = composition))+
  geom_text(aes(label = ifelse(Temp == 1600,composition,"")),hjust = 0, nudge_y = 0.2,nudge_x = 20)+
  geom_text(data = new,aes(colour = comp,label = ifelse(Temp == 1600,comp,"")),hjust = 0, nudge_y = 0.2,nudge_x = 20)+
  theme_minimal()+
  labs(title="HEA CTE",
       x="Temperature (K)",
       y=expression(paste(alpha, " ppm ", "(",K^-1,")")))+
  geom_line(size = 2)+
  geom_line(size = 2,data = new,aes(Temp, CTE, colour = comp))+
  theme(legend.position = "none")+
  lims(x=c(0,2100),y=c(0,18))+
  annotate("text",x=300+273.15,y=1,angle=90,label="PWR Operation Range",size=3,hjust=0.15)+
  geom_vline(aes(xintercept=273.15+290),lty=2)+
   geom_vline(aes(xintercept=273.15+325),lty=2)

p

anim <- p + transition_states(composition,
                    transition_length = 10,
                    state_length = 1)+
  ggtitle('HEA Linear Coefficients of Thermal Expansion',
          subtitle = 'Showing {closest_state}')+
   ease_aes('cubic-in-out')
animate(anim)