\documentstyle[11pt]{article}
\setlength{\textheight}{222mm}
\setlength{\textwidth}{160mm}
\hoffset=-1.7cm
\voffset=-1cm
\setlength{\parskip}{1ex}


\voffset=-2.0truecm
\textheight=25truecm

\begin{document}

\section{Notes on blob 020608}

\begin{itemize}

\item Check the mass of the blob to be sensible; change ampllnrho 
accordingly

\item Check about the surrounding temperature; pressure of 55 would
  give about 2700 K.

\item c1: constant flux of temperature/entropy, bcs also for magnetic
  field set!; cT: constant T (clearest)

\item Add magnetic fields

\item Can we get a rid of the ambiguous thermal conduction for now, and
set the boundary values using cs2bot/cs2top?. It seems so, at least the code
does not complain, yet.

\end{itemize}

\section{Notes on blob 030608}

The blob has a density distribution of the form 
\begin{equation}
\rho(r) = \rho_{ampl} e^{-\frac{r}{R}}
\end{equation}
Where $\rho_{ampl}$ is the amplitude of density and R is a kind of a radius of 
blob.

We can integrate the total mass of the blob in a given radius
\begin{equation}
M(r) = \int_{V}\rho_{ampl} e^{-\frac{r'}{R}}dV 
= 4\pi\rho_{ampl} \int_{0}^{r}dr' (r')^{2} e^{-\frac{r'}{R}}
\end{equation}
And we get:
\begin{equation}
M(r) = 4\pi\rho_{ampl}\Big[2R^{3} - (r^{2}+ 2Rr + 2R^{2})Re^{-\frac{r}{R}}\Big]
\end{equation}

\section{Notes on blob 050608}

(Miikka writes:)

Today I tested the blob using sound speed, blob radius, density and pressure, 
that may be roughly about what what it should be (although there are still 
some uncertanties in calculating the sound speed and pressure that
 I would like 
to talk about). It is very high (roughly 
one tenth of the speed of light).

The problem still seems to be that allthough the boundary seems very hot, the 
energy does not seem to transfer into the system from the boundary -- at least 
efficiently. (The good thing is that the code does not crash)

I'm leaving the simulation running overnight to see if the situation 
changes over time.


I'm also going to write a python-program to calculate all of the parameters.
It takes too much time to calculate them every time by hand.

\section{Notes on blob 060608}

(Miikka writes:)

The simulation I mentioned yesterday crashed. $u_{rms}$ rose to infinity and
the time step went too short.

\section{Notes on blob 170608}

(Miikka writes:)

I tried today different runs with diffent sound speeds and densities. There did
not seem to be much differences when using different sound speeds and 
densities to the overall behavior of the system.

There was though one difference. When very high sound speeds were used, the
time step was very small(too small).

\end{document}
