700 likes | 816 Views
Solar Surface Dynamics convection & waves. Bob Stein - MSU Dali Georgobiani - MSU Dave Bercik - MSU Regner Trampedach - MSU Aake Nordlund - Copenhagen Mats Carlsson - Oslo Viggo Hansteen - Oslo Andrew McMurry - Oslo Tom Bogdan - HAO O. Simulations. Computation. Solve
E N D
Solar Surface Dynamicsconvection & waves Bob Stein - MSU Dali Georgobiani - MSU Dave Bercik - MSU Regner Trampedach - MSU Aake Nordlund - Copenhagen Mats Carlsson - Oslo Viggo Hansteen - Oslo Andrew McMurry - Oslo Tom Bogdan -HAOO
Computation • Solve • Conservation equations • mass, momentum & internal energy • Induction equation • Radiative transfer equation • 3D, Compressible • EOS includes ionization • Open boundaries • Fix entropy of inflowing plasma at bottom
Method • Spatial derivatives - Finite difference • 6th order compact or 3rd order spline • Time advance - Explicit • 3rd order predictor-corrector or Runge-Kutta • Diffusion
Boundary Conditions • Periodic horizontally • Top boundary: Transmitting • Large zone, adjust <r> mass flux, ∂u/∂z=0, energy ≈ constant, drifts slowly with mean state • Bottom boundary: Open, but No net mass flux • (Node for radial modes so no boundary work) • Specify entropy of incoming fluid at bottom • (fixes energy flux) • Top boundary: B potential field • Bottom boundary: inflows advect 1G or 30G horizontal field, or B vertical
Wave Reflection Gravity wave Acoustic Wave
Radiation Transfer • LTE • Non-gray - multigroup • Formal Solution Calculate J - B by integrating Feautrier equations along one vertical and 4 slanted rays through each grid point on the surface.
Simplifications • Only 5 rays • 4 Multi-group opacity bins • Assume kLa kC
Advantage • Wavelengths with same t(z) are grouped together, so • integral over t and sum over l commute integral over t and sum over l commute
Energy Fluxes ionizationenergy 3X larger energy than thermal
Fluid Parcelsreaching the surface Radiate away their Energy and Entropy t Z r Q E S
Entropy Green & blue are low entropy downflows, red is high entropy upflows Low entropy plasma rains down from the surface
Stratified convective flow:diverging upflows, turbulent downflows Velocity arrows, temperature fluctuation image(red hot, blue cool)
Vorticity Downflows are turbulent, upflows are more laminar.
Velocity at Surface and Depth Horizontal scale of upflows increases with depth.
Velocity Distribution Up Down
Vorticity Distribution Down Up
Field Distribution observed simulation Both simulated and observed distributions are stretched exponentials.
Micropores David Bercik - Thesis
Strong Field Simulation • Initial Conditions • Snapshot of granular convection (6x6x3 Mm) • Impose 400G uniform vertical field • Boundary Conditions • Top boundary: B -> potential field • Bottom boundary: B -> vertical • Results • Micropores
Micropore Intensity image + B contours @ 0.5 kG intervals (black) + Vz=0 contours (red).
Solar velocity spectrum 3-D simulations (Stein & Nordlund) v ~ k-1/3 MDI correlation tracking (Shine) MDI doppler (Hathaway) TRACE correlation tracking (Shine) v ~ k
Line Profiles observed simulation Line profile without velocities. Line profile with velocities.
Convection produces line shifts, changes in line widths. No microturbulence, macroturbulence. Average profile is combination of lines of different shifts & widths. average profile
Spectrum of granulation Simulated intensity spectrum and distribution agree with observations after smoothing with telescope+seeing point spread function.
Stokes Image - Quiet SunSynthetic Observation - La Palma Telescope MTF + Moderate Seeing Stokes V Surface Intensity 6 Mm 6 Mm