DC add permanent discharge calculations

TS-DC-schematic/TS-DC-schematic.tex

@@ -1,7 +1,7 @@
 \documentclass[paper=A4]{article}
 \usepackage[utf8]{inputenc}
 
-\usepackage[a4paper, left=2cm, right=2cm, top=2cm, bottom=2cm]{geometry}
+\usepackage[a4paper, left=2cm, rightr2cm, top=2cm, bottom=2cm]{geometry}
 
 \usepackage{siunitx}
 \sisetup{
@@ -18,6 +18,34 @@
 \usepackage{multicol}
 \usepackage{amsmath, amssymb}
 \usepackage{hyperref}
+\usepackage{listings}
+\usepackage{xcolor}
+
+\definecolor{codegreen}{rgb}{0,0.6,0.4}
+\definecolor{codegray}{rgb}{0.5,0.5,0.5}
+\definecolor{codepurple}{rgb}{0.6,0,0.82}
+\definecolor{backcolour}{rgb}{0.92,0.92,0.92}
+
+\lstdefinestyle{mystyle}{
+    backgroundcolor=\color{backcolour},
+    commentstyle=\color{codegreen},
+    keywordstyle=\color{orange},
+    numberstyle=\tiny\color{codegray},
+    stringstyle=\color{codepurple},
+    basicstyle=\ttfamily\footnotesize,
+    breakatwhitespace=false,         
+    breaklines=true,                 
+    captionpos=b,                    
+    keepspaces=true,                 
+    numbers=left,                    
+    numbersep=5pt,                  
+    showspaces=false,                
+    showstringspaces=false,
+    showtabs=false,                  
+    tabsize=2
+}
+
+\lstset{style=mystyle}
 
 \hyphenpenalty=10000	%to stop cutting words in a paragraph
 
@@ -53,7 +81,7 @@ As seen in the Schematic, for our discharge circuitry a PTC (PTCEL13R251NxE) is
 	\label{fig:PTC_T_R}
 \end{wrapfigure}
 
-To Calculate how much discharge attempts we can have before the discharge reaches 5s, we first see at what temperature the PTC the resistance of \SI{13123}{\ohm} have. We can get this from the datasheet of the PTC (fig. \ref{fig:PTC_T_R}).
+To Calculate how much discharge attempts we can have before the discharge reaches 5s, we first see at what temperature the PTC the resistance of \SI{13123}{\ohm} have. We can get this from the datasheet \cite{ptc_datasheet} of the PTC (fig. \ref{fig:PTC_T_R}).
 
 As seen in the datasheet, the temperature is about \SI{165}{\celsius}. Assuming that the temperature rises instantly after the discharge, and the heat dissipation is negligible (since the thermal time constant $\tau_{th}$ is \SI{130}{\second}). As we are calculating the lest amount of precharge allowed, we assume the temperature of the environment to be \SI{45}{\celsius}. Since the thermal capacity $C_{th}$ is \SI{1.45}{\joule/\kelvin}, we can see that the total energy that we can generate from discharge is $E = \SI{120}{\kelvin} \cdot \SI{1.45}{\joule/\kelvin} = \SI{174}{\joule}$.
 
@@ -67,21 +95,34 @@ We then calculate how much energy in produced in one discharge:
 Therefore, the total amount of discharge we can carry out before the time exceeds \SI{5}{\second} is $\SI{174}{\joule} / \SI{16.26}{\joule} = 10.7$ --- ten times.
 
 \section*{Permanent TS Voltage}
-We can find the equilibrium temperature by finding the temperature at which the heat loss is equal to the power emitted. Using the model for LTSpice from Vishay, we can see what the heat loss is at equilibrium. (fig. ) We can see that it is about \SI{2.56}{\watt} Since the dissipation factor is \SI{19.5}{\milli\watt/\kelvin}, we can calculate what the temperature difference $\Delta T$ is.
+We can find the equilibrium temperature by finding the temperature at which the heat loss is equal to the power emitted. To find that, we first convert the graph provided in the datasheet (fig. \ref{fig:PTC_T_R}) to a Look Up Table (LUT), a python script (listed below) is then created with the two function listed below to find the equilibrium point. (Here is the $DF$ the dissipation factor, and with this PTC it is \SI{19.5}{\milli\watt/\kelvin}).
 
 \begin{align}
-	(T_{eq} - \SI{45}{\celsius}) & \cdot \SI{19.5}{\milli\watt/\kelvin} = \SI{2.56}{\watt} \\
-	T_{eq}                       & \approx \SI{176.3}{\celsius}
+	(T_{eq} - T_{amb})  \cdot DF = & P_{dissipated} \\
+	V_{TS} ^ 2 / R_{PTC} =         & P_{created}
 \end{align}
 
+After the execution of the script, we can see that the power dissipation at equilibrium is about \SI{1.84}{\watt}. The equilibrium temperature and the corresponding resistance calculated is then \SI{139}{\celsius} and \SI{88.5}{\kilo\ohm} accordingly. We can see that this is smaller then the maximum temperature rated at \SI{165}{\celsius}.
+
+To find whether the MOSFET STB10LN80K5 can survive the permanent TS voltage, we first have to calculate the current going through it.
+
+\begin{align}
+	I & = V/R = \SI{403.2}{\volt}/\SI{88.5}{\kilo\ohm} \\
+	  & = \SI{4.56}{\milli\ampere}
+\end{align}
+
+Since the MOSFET drain current $I_D$ is rated for 8A, it will work under permanent TS voltage. \cite{mosfet_datasheet}
+
+\section*{Python script}
+\lstinputlisting[language=python]{./Documents/ptc.py}
+
 \bibliographystyle{plain}
 
 \renewcommand\refname{Reference}
 \begin{thebibliography}{00}
 
-	\bibitem{lem_datasheet} \textit{Current Transducer HO-S series}. \href{https://www.lem.com/sites/default/files/products_datasheets/ho-50__250-s_v7.pdf}{https://www.lem.com/sites/default/files/products\_datasheets/ho-50\_\_250-s\_v7.pdf}, 03.2022
-	\bibitem{adz_datasheet} \textit{ADZ SME 200bar}. \href{https://www.adz.de/fileadmin/user_upload/downloads/produkte/SME/ADZ_SME_OperatingManual_DE-EN_A.pdf}{https://www.adz.de/fileadmin/user\_upload/downloads/produkte/SME/ADZ\_SME\_
-		OperatingManual\_DE-EN\_A.pdf}, 12.2010
+	\bibitem{ptc_datasheet} \textit{Vishay PTCEL13R251NxE}. \href{https://www.vishay.com/docs/29165/ptcel_series.pdf}{https://www.vishay.com/docs/29165/ptcel\_series.pdf}, 09.2024
+	\bibitem{mosfet_datasheet} \textit{ST STB10LN80K5}. \href{https://www.st.com/resource/en/datasheet/stb10ln80k5.pdf}{https://www.st.com/resource/en/datasheet/stb10ln80k5.pdf}, 02.2016
 
 \end{thebibliography}