Version 10 (modified by 9 years ago) ( diff ) | ,
---|
A closer look at the 'Flux Limiter'
is the 'flux limiter'. It comes into the equations in the diffusion term:
This is because instead of solving the following conservation equation,
(which would require an additional equation for F), we make the 'flux-limited' approximation that,
This then turns the conservation equation into a diffusion equation:
The value of
controls whether the radiation is diffusing in the free-streaming limit ( ), i.e. at the speed of light, or is diffusing as it would in the optically thick limit ( ). The FLD approximation does well at these two limits, but not in between. Let's examine these two limits in more detail and then consider why it doesn't perform as well in between. Here is the functional form of ,
where,
Graphically, we have:
Optically thick limit
From the graph above, we see that as
, . How to interpret this? R is essentially the ratio of the optical depth , to the radiation's scale height, . Thus, as , we have:
That is, the radiation travels infinitesmally small distances before it is absorbed or scattered - and thus, we are in the optically thick regime.
The rad diffusion equation then becomes,
which is consistent with equation 6.59 in Drake's book describing optically thick, non-equilibrium radiation transfer.
Free streaming limit
In the other limit,
, we have:
That is, a photon travels an infinite distance before it interacts with another particle — i.e. we are in the free-streaming limit.
In this case, the rad diffusion equation becomes:
That is, the radiation diffuses instantly through the grid. Recall, how this radiative energy couples to the gas is given by the coupling term in the radiation equation, not shown here.
Estimating the diffusion time
Estimating
…Depending on the value of R, we will get different values of
. So, we start by approximating R.Note that,
so if we make the approximation,
(by dimensional arguments and assuming the gradient is negative), we have:
which integrates to:
or,
From this equation, it is clear that h is the scale-height. Thus, by writing R as:
we can set h to be the distance between the sink and box side (h=L), and imagine it as the scale height for the radiation.
Attachments (1)
- fld.png (28.2 KB ) - added by 9 years ago.
Download all attachments as: .zip