390 likes | 560 Views
Convection Simulations. Robert Stein Ake Nordlund Dali Georgobiani David Benson Werner Schafenberger. Objectives. Understand supergranulation & the magnetic network Find signatures of emerging magnetic flux Validate and refine local helioseismic methods. Competed Simulations.
E N D
Convection Simulations Robert Stein Ake Nordlund Dali Georgobiani David Benson Werner Schafenberger
Objectives • Understand supergranulation & the magnetic network • Find signatures of emerging magnetic flux • Validate and refine local helioseismic methods
Competed Simulations • Hydrodynamic solar surface convection 96 Mm wide x 20 Mm deepDuration: 20.9 hours (bootstrapped from 48 Mm wide)Resolution: 100 km horizontal, 12-75 km vertical • 48 Mm wide x 20 Mm deepDuration: 60 hours (boostrapped from 24 Mm wide)Resolution: 100 km horizontal, 12-75 km vertical
Numerical Method • Spatial differencing • 6th-order f.d. • staggered • Time advancement • 3rd order Runga-Kutta • Equation of state • tabular • including ionization • H, He + abundant elements • Radiative transfer • 3D, LTE • 4 bin opacity distrib. fxn • Quenching
20 Mm Computational Domain 96 Mm 96 Mm Computational Domain for the CFD Simulations of Solar Convection
Mean Atmosphere Temperature, Density and Pressure (K) (105 dynes/cm2) (10-7 gm/cm2)
Mean Atmosphere Ionization of He, He I and He II
rms Velocity, O-burning convection Courtesy Dave Arnett Vertical Horizontal
Fractional area in Upflows, O burning Courtesy Dave Arnett
Vertical velocity:scan from temperature minimum to 20 Mm depthsize of cellular structures increases with depth
Streamlines: red down, blue up;15 hours; 48 x 20 Mm, vertical scale is depthsome downflows are halted, others merge into larger structures
Tracer particles: 48 x 20 Mm, 15 hours,vertical scale is grid not depth(white is fastest, dark is slowest)
Horizontal Velocity Spectrum -200 km 0 Mm 2 Mm 4 Mm 8 Mm 16 Mm
Upflows at surface come from small area at bottom (left)Downflows at surface converge to supergranule boundaries (right)
Upflows at surface come from small area at bottom (left)Downflows at surface converge to supergranule boundaries (right)
Mass Conservation Rising fluidmust turn over and descend within about ascale height to conserve mass. The actualmixing(entrainment) length in the simulation is1.8 HP.
Vorticity:11.75 hours,Finite Time Lyapunov Exponent Field, subdomain 21 Mm wide x 19 Mm high x 0.5 Mm thick, (from 48x 20 Mm simulation)
Energy balancebuoyancy work ~ dissipation(~ div FKE @ surface)
Wave generation & propagation Courtesy Junwei Zhao
k-w Diagram simulation MDI
Horizontal velocity (2-3 Mm depth):simulation (left), travel time inversion (right)
Vertical velocity (image) (km/s) & Velocity vectors in slice
Vertical velocity (image) (km/s) & Magnetic field linesin slice