wiki:u/adebrech/Matlab/RadiationPressure

Version 1 (modified by adebrech, 7 years ago) ( diff )

Lyman- Line Transfer

This is basically identical to the line transfer for ionizing radiation, except:

and the energy per photon is instead the momentum per photon,

.

To test, I've put neutral hydrogen of uniform density and temperature 104 K in a uniform gravitational field opposed by radiation pressure. To balance, equate volumetric force:

with . Clearly the value of nH doesn't affect the result, and I've chosen a small enough length scale that absorption doesn't matter and we can take across the grid.

With a flux , we get a result of (tweaking some decimals - flux in physics.data for this is actually ):

http://www.pas.rochester.edu/~adebrech/PlanetIonization/radPress/visit0002.png

after ~0.1 seconds, when acceleration under only gravity would give . Radiation pressure and gravity are clearly not perfectly matched, but it also seems clear that radiation pressure is working as it should.

Note: See TracWiki for help on using the wiki.