Changes between Version 3 and Version 4 of u/adebrech/Matlab/IonizationFront
- Timestamp:
- 09/08/17 12:39:43 (7 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
u/adebrech/Matlab/IonizationFront
v3 v4 51 51 where dx is the cell size in cm, and we are solving for X(c), the ionization fraction in the current cell. A modified logistic curve was then fit approximately (i.e. by hand) to the ionization fraction solution and used as the initial conditions for a simulation, which was allowed to relax to true equilibrium. A better fit for initial conditions could be obtained with a minimization routine, but I think it is more instructive to see the simulation have some significant relaxation. 52 52 53 [[Image(http://www.pas.rochester.edu/~adebrech/PlanetIonization/ionization_eq.png,link=http://www.pas.rochester.edu/~adebrech/PlanetIonization/ionization_eq.gif)]] 54 53 55 One thing to note is that the mean free path of a photon in the ionization front is not particularly comparable to the actual width of the ionization front in this case, since the ionization front cannot be considered thin on our scale.