Changes between Version 15 and Version 16 of FluxLimitedDiffusion


Ignore:
Timestamp:
03/19/13 14:16:39 (12 years ago)
Author:
Jonathan
Comment:

—

Legend:

Unmodified
Added
Removed
Modified

  • FluxLimitedDiffusion

    v15 v16  
    8282Krumholz et al. perform Implicit Radiative, Explicit Hydro, Explicit Radiative
    8383In AstroBEAR this would look like:
    84  * Prolongate, d, p, E, Edot
    85  * Overlap d, p, E
    86  * Do IR which updates e[0], and E[0] using d[1], e[1], E[1]
    87  * Update Edot[0] using pre IR and post IR E[0]
    88  * Ghost e[2*mbc], E[mbc+1], Edot[mbc+1]
    89  * Do first EH[mbc]
    90  * Do ER[mbc] --- Terms with grad E can be done without ghosting since EH did not change E.  The del dot vE term needs time centered face centered velocities which can be stored during the hydro update.
     84* Initialization
     85 * Prolongate, d, p, e, E, Edot
     86* Step 1
     87 * Overlap d, p, e, E and do physical BC's
     88 * Do IR which updates e,,0,,, and E,,0,, using d,,1,,, e,,1,,, E,,1,,
     89 * Update E,,1,, using Edot,,1,,
     90 * Update Edot,,0,, using pre IR and post IR E,,0,,
     91 * Ghost e,,2*mbc,,, E,,mbc+1,,, Edot,,mbc+1,,
     92 * Do first EH,,mbc,,
     93 * Do ER,,mbc,, --- Terms with grad E can be done without ghosting since EH did not change E.  The del dot vE term needs time centered face centered velocities which can be stored during the hydro update.
    9194 * Store Edot in child arrays to be prolongated
    92  * Overlap
    93  * Do IR which updates e[0], and E[0] using d[1], e[1], E[1]
    94  * Update Edot[0] using pre IR and post IR E[0]
    95  * Ghost e[mbc], E[1], Edot[1]
    96  * Do second EH[0]
    97  * Do ER[0] --- Terms with grad E can be done without ghosting since EH did not change E.  The del dot vE term needs time centered face centered velocities which can be stored during the hydro update.
     95* Step 2
     96 * Overlap d, p, e, E, and do physical BCs
     97 * Do IR which updates e,,0,,, and E,,0,, using d,,1,,, e,,1,,, E,,1,,
     98 * Update Edot,,0,, using pre IR and post IR E,,0,,
     99 * Update E,,1,, using Edot,,1,,
     100 * Ghost e,,mbc,,, E,,1,,, Edot,,1,,
     101 * Do second EH,,0,,
     102 * Do ER,,0,, --- Terms with grad E can be done without ghosting since EH did not change E.  The del dot vE term needs time centered face centered velocities which can be stored during the hydro update.
    98103
    99104