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) |