tionis said:
PAllen & Vanadium: do you guys agree with this plot?
First, a quick question that you've never answered (I'm afraid I'm fearing the worst). Can you view the reference link I posted, and did you read it?
You don't define your variables, but I think the idea is that BB is the spectral radiance, as per
http://en.wikipedia.org/wiki/Planck's_law
If you refer to the paper, which I *really hope* you're reading, you'll see that the paper uses the alternate formula from the wiki, where they use ##\lambda = \frac{c}{\nu}## instead of ##\nu##, ie in terms of wavelength rather than frequency.
This gives eq 12 of the above paper, which is equivalent to your own but in terms of wavelength rather than frequency.
eq 13 gives the actual received power / unit area, if the stars radius is a and it's distance is R.
Your transformation to the moving frame is incorrect, however.
The spectral radiance in the moving frame is D^3 times the spectral radiance in the stationary frame - see eq 11.
When you integrate this out, you find that the total energy recieved, integrating from lambda from 0 to infinity should scale as D^2. (D being the doppler factor). At least if you can get the integrals to work out - maple doesn't want to do them for me, and I don't even want to attempt to do them by hand.
The appendix to the paper gives a short derivation of why D^2 is correct for the total intensity. You can refer to the "photon arrival" thread for perhaps a clearer discussion, the gist is that the shift in energy per photon gives one factor of D, and that the photon arrival rate also increases by the doppler factor D, giving a total intensity increase of D^2. You can also try the Wiki article on "relativistic beaming"
http://en.wikipedia.org/wiki/Relativistic_beaming.
If you evaluate your integral that you give in the moving frame, you should see that the intensity does NOT increase as D^2, but much faster. (If you can get the integral to evalutate, that is - good luck with that!. )
But anyway that is where you appear to have "gone wrong".
[add]
WHile I can't do the integrals myself, I can point out that due to the Steffan-Boltman law
http://en.wikipedia.org/wiki/Stefan–Boltzmann_constant
one expects the total power emitted to be proportional to T^4. Therefore, when you multiply the temperature by D, you get a radiant power increase of D^4. To get the correct transformation law, you need to not only multiply the temperature by D, but divide the intensity by D^2, so that the radiant intensity scales by D^2.