Context Navigation

Meeting Update 1/21/2013

Worked with Christine in setting up her problem module.

Uses AmbientObjects, DiskObjects, RefinementObjects, ProjectionObjects, ProfileObjects, and implements the stream the ol' fashioned way - though we may want to make a StreamObjects or a JetObjects to do 'jets' from a boundary.

Ran into problems with the boundaries due to cooling time scales being short and DM Cooling has no equilibrium temps.
Solution was to have very high density contrast between disk (1e-7 g/cc) and outer edge of ambient (1e-20 g/cc).
Also wanted to have hydrostatic equilibrium without requiring temperature extremes - so set up isothermal solution… \rho=\rho_\infty e^{-\phi/\left( P_\infty / \rho_\infty \right )} T=P_\infty/\rho_\infty
First we added a constant pressure to p_inf to get the desired pressure at the disk edge - but since the density changes, this makes the ambient no longer isothermal.
IF we try to solve for P_\infty to give pressure matching at the disk edge we end up with a transcendental equation
P(R_{disk})=P_\infty e^{-\phi(R)/\left( P_\infty/\rho_\infty \right )}
If instead we fix T or P_\infty / \rho_\infty and solve for \rho_\infty we get an easier equation to solve.
P(R_{disk})=\rho_{\infty}(P_\infty/\rho_\infty) e^{-\phi(R)/\left( P_\infty/\rho_\infty \right )}

\rho_{\infty} = P(R_{disk})(\rho_\infty/P_\infty) e^{\phi(R)/\left( P_\infty/\rho_\infty \right )}