Changes between Version 1 and Version 2 of u/adebrech/Matlab/RadiationPressure


Ignore:
Timestamp:
10/30/17 09:34:01 (7 years ago)
Author:
adebrech
Comment:

—

Legend:

Unmodified
Added
Removed
Modified

  • u/adebrech/Matlab/RadiationPressure

    v1 v2  
    1 == Lyman-$\alpha$ Line Transfer ==
     1= Lyman-$\alpha$ Line Transfer =
    22
    33This is basically identical to the line transfer for ionizing radiation, except:
    … …  
    2424
    2525$v \approx 0.1 \text{ cm/s}$ after ~0.1 seconds, when acceleration under only gravity would give $v = g t = 980(0.1) = 98 \text{ cm/s}$. Radiation pressure and gravity are clearly not perfectly matched, but it also seems clear that radiation pressure is working as it should.
     26
     27== Dynamic tests ==
     28
     29|| Lyman-$\alpha$ flux || 5.1d13 phot/cm^2^/s ||
     30|| lScale || 1 ||
     31
     32__Slab__
     33
     34|| Slab density || 1 CD ||
     35|| Ambient density || 1d-4 CD ||
     36|| Slab extent || (0,0.3) ||
     37|| Domain || {(0,1),(0,1)} ||
     38
     39__Clump__
     40
     41|| Clump density || 1 CD ||
     42|| Ambient density || 1d-4 CD ||
     43|| Clump extent || Circle around (0,0) with R = 0.5 ||
     44|| Domain || {(-1,1),(-1,1)} ||
     45
     46__Slab with ionization__
     47
     48|| Slab density || $1d4\cdot10^{-20(0.5-y)^2}$ CD ||
     49|| Ionizing flux || 2d8 || <-- Need to test value on smaller length scale