| 15 | | Dobbs Dixon '08 |
| 16 | | Domain |
| 17 | | * Radius goes from 7.95e9 cm to 8.65e9 cm (1.06 -1.2 R_j) |
| 18 | | * Single Temperature (assumes Planck opacity is >> ) |
| 19 | | * They use [[latex($\lambda = \frac{2+R}{6+3R+R^2}$)]] instead of [[latex($\frac{1}{R}\left ( \coth{R}-\frac{1}{R} \right )$)]] |
| 20 | | [[Image(Limiter.png, width=400)]] |
| 21 | | * [[latex($\gamma=7/5$)]] |
| 22 | | * [[latex($\chi=2.3$)]] |
| 23 | | * Roseland opacities found using Pollack et al 1985 for low T with Alexander & Ferguson for high T |
| 24 | | * Initial Density and Temperature profiles taken from 1D model (B1 of Bodenheimer 2001) (.63 M,,J,, for 4.5 Gyr with a 1200 degree outer boundary) |
| 25 | | * Outer boundary set to T= max(1200 K * [cos(phi)*cos(theta)]^(1/4), 100 K) |
| | 24 | == Module supports == |
| 30 | | |
| 31 | | |
| 32 | | Dobbs Dixon'10 |
| 33 | | * Added two temperature model |
| 34 | | * Added viscosity |
| 35 | | * Added stellar heating using optical depth and high energy planck opacity |
| 36 | | |
| 37 | | |
| 38 | | Current differences between Ian's code and AstroBEAR |
| 39 | | * AstroBEAR folds in gas/radiation energy coupling into implicit solve... Hypre has both kappa,,R,, and kappa,,P,, terms - with sub-cycling for changes in source function. |
| 40 | | * AstroBEAR also assumes 'kappa' is constant over implicit update... assumes [[latex($\frac{\partial \kappa}{\partial T} = 0$)]]... |
| 41 | | * No Kappa_A terms (stellar heating) in AstroBEAR |
| 42 | | * Need to develop communication patterns for calculating optical depth from star to get heating term correct. Could use initial spherical profile to estimate optical depth as a function of position? Not sure how this is done in Ian's code |
| 43 | | * Viscosity may need to be incorporated later... Shule is developing this. |
| 44 | | |
| 45 | | * Second AMR implicit solve may need to be added later (ie Howell and Greenough 2002) |
| 46 | | |
| 47 | | == Initial Params == |
| 48 | | |
| 49 | | * Density Profile, Temperature Profile |
| 50 | | * Opacity Tables |
| 51 | | * Planet held together by self gravity |
| 52 | | * Rotating frame about COM |
| 53 | | * |
| | 32 | == Still working on == |
| | 33 | * Basic testing of hydrostatic equilibrium for planet |
| | 34 | * Line transfer for stellar heating |
| | 35 | * Second AMR implicit solve may need to be added later (ie Howell and Greenough 2002) |