Planetary Atmospheres

Updated boundary conditions to go from rho = 1e-9 to tau = 1e4

Inner boundary: (tau=1e4)

  • r=7.815240e+09
  • rho= 4.972392e-03
  • T= 2.306330e+03
  • M_enclosed=2.173174e+30
  • P=4.145923e+08

Point particle

  • M_point=2.1632e30 = M_enclosed - V(r)*rho

Outer Boundary (rho < 1e-9)

  • r=8.181350e+09
  • rho=9.488437e-10
  • T=1.115419e+03
  • P=3.826197e+01

Also, calculated the effective aspect ratio for simulation region: 'w' = 75.672

which gives Ncells = Nr3w2 = 5726 Nr3 = (18 Nr)3

So if we want 100 zones in radius, we would be effectively simulating an 18003 sim.

Here are the previous values I used

Inner boundary: (tau=10)

  • r=7.9568e9
  • rho=3.8901e-5
  • T=1.4682e3
  • M_enclosed=2.1733e30
  • P=2.0648e6 (2.0648 atm)

Point particle

  • M_point=2.1732e30 = M_enclosed - V(r)*rho

Outer Boundary (rho < 1e-9)

  • r=8.1814e9 (3% of inner boundary)
  • rho=9.4884e-10
  • T=1.1154e+03
  • P=38.2620

PLoit showing mesh

Plot showing density (click for movie)

Plot showing temperature (click for movie)

Attachments (6)

Comments

No comments.