Initialisation

I will list out all relevant general parameters that I have initialised for my galaxy(in normalised form):

1. Resolution of $z$ is $dz = 0.04$
2. Resolution of $t$ is $dt = 0.00079$
3. $\alpha = sin(\pi z)$
4. I solve my equations for $r = 0.4kpc$
5. When working on Task 3, I evaluate $B_{tot}$ at $z=0$
6. I work with Dirichlet Boundary conditions for Task 3

For completeness, I will list all scaling parameters again:

$$ \begin{array}{|c | c | c |}\hline \text{Scaling Constant} & \text{ Value } & \text{Unit} \\ \hline R_0 & 10 & kpc \\ \hline h_0 & 0.5 & kpc \\ \hline B_0 & 1e-6 & G \\ \hline \alpha_0 & 2000 & cm s^{-1}\\ \hline \Omega_0 & 100 & km s^{-1} kpc^{-1} \\ \hline S_0 & -50 & km s^{-1} kpc^{-1} \\ \hline \eta_T & 1e26 & cm^2 s^{-1} \\ \hline t_0 = h_0^2/\eta_T & 0.745 & Gyr \\ \hline r_\Omega & 2 & kpc \\ \hline \end{array} $$

Task 2 Results

First, I will review some results from Task 2, where I solved the $\alpha-\omega$ dynamo equations. I was able to match my results with those provided by the Teaching assistant. I evolved my system for around $15t_0$ or $11.9 Gyr$.

The seed field I used is given in figure (1). I performed many simulations with different dynamo numbers to understand the evolution of $B_{tot}$.

Figure 1

Figure 2: Red - Just a constant line, to help distinguish growing and decaying solutions