\section{MatLab}
\label{V:Matlab}

\subsection{Kode}
I det følgende er de ulike matlabskriptene vist.

\subsubsection{Fordamper.m}
\label{V:Mford}
Beregner entalpidifferanse over fordamper ved ulike fordampertemperaturer. Gitt: T og p før ventil, T etter ventil (fordampertemperatur) og grad av overheting. Bruker skriptene \begin{verbatim}pH_co2_pr.m\end{verbatim} og \begin{verbatim}pH_co2_sw.m\end{verbatim} i \ref{V:pHpr} og \ref{V:pHsw} for beregning av volum på metningskurven.
\lstinputlisting{Vedlegg/Matlab/Fordamper.m}

\newpage
\subsubsection{Kompressor.m}
\label{V:Mkomp}
Beregner entalpidifferanse over kompressor. Gitt: T og p før kompressor, p etter kompressor samt at entropi er konstant.
\lstinputlisting{Vedlegg/Matlab/Kompressor.m}

\newpage
\subsubsection{co2 pr.m}
\label{V:pr}
Fra databasen DIPPR96\cite{dippr} hentes data for \co2 for å generere en standardtilstand, en ideell tilstandslikning og et reelt residualledd til tilstandslikningen (i dette tilfellet Peng-Robinson). Standardtilstanden genereres ved hjelp av kjemisk potensial, dannelsesentalpi og standard entropi for \co2. I den reelle tilstandslikningen benyttes verdier som kritisk trykk og temperatur.



\begin{verbatim}
### Peng-Robinson CO2 equation of state

n   = ['carbon dioxide']

co2 = Surface.new('co2_pr') * (
         Helmholtz.new() * (
           StandardState.new() * (
             MuT_cp.new(:ig,:dippr,:dippr96) * (
               MuT_hs.new(:ig,:h0,:dippr96) +
               MuT_hs.new(:ig,:s0,:dippr96).mixture(n)
             )
           ) +
           EquationOfState.new() * (
             ModTVN_ideal.new(:fluid,:idealgas).mixture(n)
           ) +
           EquationOfState.new() * (
             ModTVN.new(:fluid,:pr,:dippr96,:home).mixture(n)
           )
         )
       )

co2.LaTeX_report!

:sw12polar,:sw03
\end{verbatim}

\lstinputlisting{Vedlegg/Matlab/co2_pr.m}


\newpage
\subsubsection{co2 sw.m}
\label{V:sw}
\lstinputlisting{Vedlegg/Matlab/co2_sw.m}

\newpage
\subsubsection{pH co2 pr.m}
\label{V:pHpr}
\lstinputlisting{Vedlegg/Matlab/pH_co2_pr.m}

\newpage
\subsubsection{pH co2 sw.m}
\label{V:pHsw}
\lstinputlisting{Vedlegg/Matlab/pH_co2_sw.m}

\newpage
\subsection{Entalpidifferanser}
\subsubsection{Fordamper}
\label{V:Hford}
I figur \ref{fig:VMat1}-\ref{fig:VMat5} er entalpidifferanse ved ulik grad av overheting plottet mot ulike fordampertemperaturer. 

\begin{figure}[!ht]
\centering
\includegraphics[width=15cm,height=!]{Figurer/PNG/Matlab/ResFordOH0.png}\caption{Entalpidifferanse over fordamper, ingen overheting.}
\label{fig:VMat1}
\end{figure}

\begin{figure}[!ht]
\centering
\includegraphics[width=15cm,height=!]{Figurer/PNG/Matlab/ResFordOH5.png}\caption{Entalpidifferanse over fordamper, overheting: 5K.}
\label{fig:VMat2}
\end{figure}

\begin{figure}[!ht]
\centering
\includegraphics[width=15cm,height=!]{Figurer/PNG/Matlab/ResFordOH10.png}\caption{Entalpidifferanse over fordamper, overheting: 10K.}
\label{fig:VMat3}
\end{figure}

\begin{figure}[!ht]
\centering
\includegraphics[width=15cm,height=!]{Figurer/PNG/Matlab/ResFordOH15.png}\caption{Entalpidifferanse over fordamper, overheting: 15K.}
\label{fig:VMat4}
\end{figure}

\begin{figure}[!ht]
\centering
\includegraphics[width=15cm,height=!]{Figurer/PNG/Matlab/ResFordOH20.png}\caption{Entalpidifferanse over fordamper, overheting: 20K.}
\label{fig:VMat5}
\end{figure}
