Numerical Simulation of Magnetoconvection in a Stellar Envelope

Total Page:16

File Type:pdf, Size:1020Kb

Numerical Simulation of Magnetoconvection in a Stellar Envelope Center for Turbulence Research 281 Annual Research Briefs 2002 Numerical simulation of magnetoconvection in a stellar envelope By S. D. Ustyugov AND A. N. Andrianov y 1. Introduction One of the main problems in the physics of the Sun is the interaction between con- vection and magnetic ¯elds. Convective motions in stellar envelopes span a substantial radial distance from the convectively unstable region into the adjacent stable zones and can distort the magnetic ¯eld into concentrated flux sheets and tubes in which the mag- netic pressure is comparable to the gas pressure (see Hurlburt & Toomre (1988)). On the other hand, magnetic ¯elds can be su±ciently strong to suppress convection on granu- lar and supergranular scales. The dynamics of such couplings is strongly nonlinear, and the flow spans many scale heights in the vertical direction and weakens the stable ther- mal strati¯cation. Studies of convection in the presence of magnetic ¯elds have already shown e®ects arising from compressibility, revealing distinctive asymmetries between up- ward and downward flows. Convection penetrates into the underlying stable layers as downward-directed plumes (see Hurlburt et.al. (1986)). Here we consider compressible convection in three spatial dimensions in the presence of an externally imposed magnetic ¯eld. We study penetrative convection within relatively simple con¯gurations consisting of two layers with well-posed boundary conditions. We include all di®usive processes and do not model the unresolved scales. 2. Initial model and equations We shall consider penetrative convection in a compressible stellar envelope in the pres- ence of an imposed magnetic ¯eld. We assume that this envelope experiences uniform gravitational acceleration directed downward, and possesses constant thermal conductiv- ity, a constant magnetic di®usivity, and a constant shear viscosity. The strati¯cation consists of an unstable layer bounded below by a stable polytrope. This con¯guration mimics the lower part of the solar convection zone and the upper radiative zone. The initial distribution of temperature, density, and pressure are expressed as m® m®+1 z ¡ z® ½ T P T T = T® + ; = ; = (2.1) K® ½® µT® ¶ P® µT® ¶ for the unstable (® = 1) and stable (® = 2) layers. Here the polytropic indices m® are given by m® + 1 = gK®=RFT , where R is the gas constant, K® is the coe±cient of thermal conductivity, and FT is the total energy flux. Convective instability occurs when m® < mc, where mc, the critical polytropic index, depends on the ratio of speci¯c heats y Keldysh Institute of Applied Mathematics, Moscow, Russia 282 S. D. Ustyugov & A. N. Andrianov γ. For perfect a monoatomic gas (γ = 5=3), mc = 3=2. We take m1 = 1 in the unstable layer and m2 = 3 in the stable layer. We have solved the equations of three-dimensional, fully compressible nonlinear mag- netoconvection. We use cartesian coordinates x = x1; y = x2; z = x3, where the positive z-axis points downwards. The equations of motion are as follows: conservation of mass, @½ @ + (½ui) = 0; (2.2) @t @xi conservation of momentum, @½ui @ 1 + ±ij P + BkBk + ½uiuj ¡ ¿ij ¡ BiBj ¡ ½g±i3 = 0; (2.3) @t @xj · µ 2 ¶ ¸ conservation of energy, @½E @ 1 P @T + ½uj uiui + ² + ¡ gx3 ¡ ui¿ji ¡ K @t @xj · µ2 ½ ¶ @xj ¸ (2.4) @ @Bi @Bj ¡ ´Bi ¡ + Bi (Bjui ¡ Biuj) = 0; @xj · µ @xj @xi ¶ ¸ and the induction equation, @Bi @ @Bi @Bj + ujBi ¡ uiBj ¡ ´ ¡ = 0: (2.5) @t @xj · µ @xj @xi ¶¸ These are augmented by the equation of state for a perfect gas, P = R½T; (2.6) and expressions for the total energy per unit mass, 1 1 BiBi E = ² + u u ¡ gx3 + ; (2.7) 2 i i 2 ½ and the viscous stress tensor, @ui @uj 2 @ul ¿ij = ¹ + ¡ ±ij : (2.8) µ@xj @xi 3 @xl ¶ Here ½, T , P , and ² are density, temperature, pressure, and speci¯c internal energy, and Bi denotes components of the magnetic ¯eld. We write the equations in dimensionless form and take our unit of length to be the depth of the unstable layer d. The time unit d= (Rd¯1) is related to the sound travel time across the unstable layer, where ¯1 is the initial temperature gradient within the unstable layer. The magnetic ¯eld is scaled by B^, the value of the initially uniform vertical magnetic ¯eld. The evolution of the system depends on some physical and geometrical parameters. We de¯ne the computational box to extend in a horizontal plane from x = y = 0 to x = y = A, where A is the aspect ratio of the unstable layer, and vertically from z = z1 Magnetoconvection in a stellar envelope 283 to z = z3. The depth of the unstable layer is unity in our units. The density contrast Â2 is the ratio of the density at the bottom of the unstable layer (z = z2) to that at the top: ½ (z2) Â2 = (2.9) ½ (z1) The Prandtl number, which gives the ratio of viscous to thermal di®usivity, is de¯ned as ¹c σ = p (2.10) K The degree of instability in the unstable layer is given by the Rayleigh number ¡1 ¡1 2 2m1¡1 Ra = (m1 + 1) (mc ¡ m1) (mc + 1) σ1 ¸ z ; (2.11) where ¸ is the ratio of the sound speed to the thermal di®usion time and σ1 is the Prandtl number in the unstable layer. Since Ra depends on vertical position in the unstable layer, we evaluate it at the center, taking 1 R^ = R z1 + : (2.12) a a µ 2¶ as a nominal Rayleigh number. The strength of the imposed magnetic ¯eld and its e®ect on convective stability can be measured by the Chandrasekhar number d2 Q = B2 (2.13) ¹´ The magnetic Prandtl number (the ratio of magnetic to thermal di®usivity) is de¯ned as ³ = ´½cp=K (2.14) 3. Numerical method For the numerical simulation we have used an explicit TVD method, second order ac- curate in space with the time advance using a three-step Runge-Kutta scheme (see Yee et.al. (1990), Shu (1989)). Our algorithm uses equations in conservative form for advec- tion, wave propagation, and di®usion processes. The 3-D compressible time-dependent MHD equation with source term can be written as @U @F @G @H + + + = S; (3.1) @t @x @y @z The dependent variable U is vector of conserved variables, and F; G; H are vectors of flux in the three space directions. Let Ut = L(U)i;j;k; (3.2) 284 S. D. Ustyugov & A. N. Andrianov be the semi-discrete form of (3.1) at point (i; j; k), where L is the spatial discretization operator 1 1 L = ¡ F +1 2 ¡ F ¡1 2 ¡ G +1 2 ¡ G ¡1 2 i;j;k ¢x i = i = ¢y j = j = ¡ ¢ ¡ ¢ (3.3) 1 ¡ H +1 2 ¡ H ¡1 2 + S ; ¢z k = k = i;j;k ¡ ¢ The numerical flux, for example in the x direction, is de¯ned as 1 F +1 2 = F +1 + F + P +1 2© +1 2 : (3.4) i = 2 i i i = i = ¡ ¢ The last term Pi+1=2©i+1=2 is a nonlinear dissipation, and the quantity Pi+1=2 is the right eigenvector matrix @F=@U using, for example, Roe's approximate average state. l The components Ái+1=2 of vector ©i+1=2 can be written as 1 Ál = Ã(al )(gl + gl) ¡ Ã(al + γl )®l ; (3.5) i+1=2 2 i+1=2 i+1 i i+1=2 i+1=2 i+1=2 1 γl = Ã(al )(gl ¡ gl)=®l if ®l 6= 0; (3.6) i+1=2 2 i+1=2 i+1 i i+1=2 i+1=2 l l γi+1=2 = 0 if ®i+1=2 = 0: (3.7) l The ai+1=2; l = 1; :::; 7, are the characteristic speeds of @F=@U evaluated at some l symmetric average of Ui and Ui+1. The function à is an entropy correction to j ai+1=2 j. l ¡1 The ®i+1=2 are elements of Ri+1=2(Ui+1 ¡ Ui). l The limiter function gi can be expressed as l l l 1 l l g = minmod(2a ¡ ; 2a ; (a ¡ + a )): (3.8) i i 1=2 i+1=2 2 i 1=2 i+1=2 For the thermal conductivity we have applied an explicit multistage numerical scheme. This scheme is absolutely stable and is described in detail by Zhukov et.al. (1993). The numerical scheme has been tested widely and used successfully in earlier work (see Ustyu- gov & Andrianov (2002)). We have carried out computations in a cube with 128 mesh points in each direction on 256 processors of the parallel machine with cluster architec- ture MBC 1000M. The computer program was parallelized using the NORMA system (see Andrianov et.al. (2001)). The Norma program was then compiled using Fortran and the MPI library. A typical simulation requires on this cluster 50 hours to advance the solution well past its start-up transient. 4. Boundary and Initial conditions We impose periodic boundary conditions with period A for all variables in x and y. On the upper and lower surfaces we use stress-free boundary conditions: Bx = By = @Bz=@z = uz = @ux=@z = @uy=@z = 0; (4.1) Magnetoconvection in a stellar envelope 285 i.e., vertical velocities, together with the horizontal components of both the viscous and magnetic stresses, vanish at the upper and lower surfaces. The temperature is held ¯xed on the upper surface, and its vertical derivative is imposed on the lower. These conditions require T = T1 at z = z1; (4.2) @T=@z = 1=K2 at z = z3; (4.3) We introduce a small-amplitude velocity perturbation over several wavenumbers within the unstable layer.
Recommended publications
  • Structure and Energy Transport of the Solar Convection Zone A
    Structure and Energy Transport of the Solar Convection Zone A DISSERTATION SUBMITTED TO THE GRADUATE DIVISION OF THE UNIVERSITY OF HAWAI'I IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF DOCTOR OF PHILOSOPHY IN ASTRONOMY December 2004 By James D. Armstrong Dissertation Committee: Jeffery R. Kuhn, Chairperson Joshua E. Barnes Rolf-Peter Kudritzki Jing Li Haosheng Lin Michelle Teng © Copyright December 2004 by James Armstrong All Rights Reserved iii Acknowledgements The Ph.D. process is not a path that is taken alone. I greatly appreciate the support of my committee. In particular, Jeff Kuhn has been a friend as well as a mentor during this time. The author would also like to thank Frank Moss of the University of Missouri St. Louis. His advice has been quite helpful in making difficult decisions. Mark Rast, Haosheng Lin, and others at the HAO have assisted in obtaining data for this work. Jesper Schou provided the helioseismic rotation data. Jorgen Christiensen-Salsgaard provided the solar model. This work has been supported by NASA and the SOHOjMDI project (grant number NAG5-3077). Finally, the author would like to thank Makani for many interesting discussions. iv Abstract The solar irradiance cycle has been observed for over 30 years. This cycle has been shown to correlate with the solar magnetic cycle. Understanding the solar irradiance cycle can have broad impact on our society. The measured change in solar irradiance over the solar cycle, on order of0.1%is small, but a decrease of this size, ifmaintained over several solar cycles, would be sufficient to cause a global ice age on the earth.
    [Show full text]
  • The Sun and the Solar Corona
    SPACE PHYSICS ADVANCED STUDY OPTION HANDOUT The sun and the solar corona Introduction The Sun of our solar system is a typical star of intermediate size and luminosity. Its radius is about 696000 km, and it rotates with a period that increases with latitude from 25 days at the equator to 36 days at poles. For practical reasons, the period is often taken to be 27 days. Its mass is about 2 x 1030 kg, consisting mainly of hydrogen (90%) and helium (10%). The Sun emits radio waves, X-rays, and energetic particles in addition to visible light. The total energy output, solar constant, is about 3.8 x 1033 ergs/sec. For further details (and more accurate figures), see the table below. THE SOLAR INTERIOR VISIBLE SURFACE OF SUN: PHOTOSPHERE CORE: THERMONUCLEAR ENGINE RADIATIVE ZONE CONVECTIVE ZONE SCHEMATIC CONVECTION CELLS Figure 1: Schematic representation of the regions in the interior of the Sun. Physical characteristics Photospheric composition Property Value Element % mass % number Diameter 1,392,530 km Hydrogen 73.46 92.1 Radius 696,265 km Helium 24.85 7.8 Volume 1.41 x 1018 m3 Oxygen 0.77 Mass 1.9891 x 1030 kg Carbon 0.29 Solar radiation (entire Sun) 3.83 x 1023 kW Iron 0.16 Solar radiation per unit area 6.29 x 104 kW m-2 Neon 0.12 0.1 on the photosphere Solar radiation at the top of 1,368 W m-2 Nitrogen 0.09 the Earth's atmosphere Mean distance from Earth 149.60 x 106 km Silicon 0.07 Mean distance from Earth (in 214.86 Magnesium 0.05 units of solar radii) In the interior of the Sun, at the centre, nuclear reactions provide the Sun's energy.
    [Show full text]
  • 198Lapj. . .243. .945G the Astrophysical Journal, 243:945-953
    .945G .243. The Astrophysical Journal, 243:945-953, 1981 February 1 . © 1981. The American Astronomical Society. All rights reserved. Printed in U.S.A. 198lApJ. CONVECTION AND MAGNETIC FIELDS IN STARS D. J. Galloway High Altitude Observatory, National Center for Atmospheric Research1 AND N. O. Weiss Sacramento Peak Observatory2 Received 1980 March 10; accepted 1980 August 18 ABSTRACT Recent observations have demonstrated the unity of the study of stellar and solar magnetic fields. Results from numerical experiments on magnetoconvection are presented and used to discuss the concentration of magnetic flux into isolated ropes in the turbulent convective zones of the Sun or other late-type stars. Arguments are given for siting the solar dynamo at the base of the convective zone. Magnetic buoyancy leads to the emergence of magnetic flux in active regions, but weaker flux ropes are shredded and dispersed throughout the convective zone. The observed maximum field strengths in late-type stars should be comparable with the field (87rp)1/2 that balances the photospheric pressure. Subject headings: convection — hydrodynamics — stars: magnetic — Sun : magnetic fields I. INTRODUCTION In § II we summarize the results obtained for simplified The Sun is unique in having a magnetic field that can be models of hydromagnetic convection. The structure of mapped in detail, but magnetic activity seems to be a turbulent magnetic fields is then described in § III. standard feature of stars with deep convective envelopes. Theory and observation are brought together in § IV, Fields of 2000 gauss or more have been found in two where we discuss the formation of isolated flux tubes in late-type main sequence stars (Robinson, Worden, and the interior of the Sun.
    [Show full text]
  • Single and Binary Evolution of Population III Stars and Their Supernova Explosions T
    Mon. Not. R. Astron. Soc. 384, 1533–1543 (2008) doi:10.1111/j.1365-2966.2007.12810.x Single and binary evolution of Population III stars and their supernova explosions T. M. Lawlor,1 T. R. Young,2 T. A. Johnson3 and J. MacDonald4 1Department of Physics, Penn State University, Brandywine, Media, PA 19063, USA 2Department of Physics, University of North Dakota, Grand Forks, ND 58202, USA 3Department of Math, Science & Technology, University of Minnesota, Crookston, Crookston, MN 56716 4Department of Physics and Astronomy, University of Delaware, Newark, DE 19716, USA Accepted 2007 November 30. Received 2007 November 30; in original form 2007 October 12 ABSTRACT We present stellar evolution calculations for Population III stars for both single- and binary- star evolutions. Our models include 10- and 16.5-M single stars and a 10-M model star that undergoes an episode of accretion resulting in a final mass of 16.1 M. For comparison, we present the evolution of a solar heavy element abundance model. We use the structure from late-stage evolution models to calculate simulated supernova light curves. Light curve comparisons are made between accretion and non-accretion progenitor models, and models for single-star evolution of comparable masses. Where possible, we make comparisons to previous works. Similar investigations have been carried out, but primarily for solar or near-solar heavy metal abundance stars and not including both the evolution and the supernova explosions in one work. Key words: binaries: general – stars: evolution – supernovae: general. low-mass, very low Z object He 0107–5240 ([Fe/H] =−5.3) as an 1 INTRODUCTION example of a possible Pop III star that makes it difficult to rule out It is commonly believed that the first stars, Population III (Pop III) low-mass fragments.
    [Show full text]
  • Radiation Zone
    AST 100 General Astronomy: 5. What’s inside the Sun? Stars & Galaxies From the Center Outwards • Core: Hydrogen fusing into helium, releasing energy in ANNOUNCEMENTS the form of gamma rays, neutrinos, and positrons Midterm I on Tue, Sept. 29 it will cover class material up to today (included) • Core temp = 15 million K, hot & dense from the gravitational Quiz #2 today, end of class weight of all that mass Meanderings of outbound photons Our gamma-ray photons random walk outwards (getting Radiation Zone redirected with every step), gradually cooling • Gamma ray photons leave the core and move into an area known as the Takes hundreds of Radiation Zone thousands to a million – Neutrinos? They leave right years to get out!! away, no interaction – Positrons? Quickly find electrons in the core to annihilate with. • Photons only travel about 1 mm before being redirected in another direction T=10 million K. Photosphere Convection Zone • At the top of the convection • Eventually, gas is cool zone, the densities are now enough (2 million K at the low enough that our boundary) and becomes photons can zoom away. turbulent – Now downgraded all the – No longer just redirects way to visible energies photons, now absorbs them • Photosphere is the “visible surface” of the Sun • Convection: hotter • T = only 5800 K regions rise, cooler • Photons free - seen at Earth regions sink 8 min later • Blackbody spectrum (T= • Energy continues to work 5800 K) + absorption from its way out - nearly 1 cooler gasses just on top million years to get out Granulation: turbulent convection Granulation Movie Appearance of the photosphere • Typical granulations last only 8 - 15 minutes • Movie covers 35 min Taken by G.
    [Show full text]
  • Solar Photosphere and Chromosphere
    SOLAR PHOTOSPHERE AND CHROMOSPHERE Franz Kneer Universit¨ats-Sternwarte G¨ottingen Contents 1 Introduction 2 2 A coarse view – concepts 2 2.1 The data . 2 2.2 Interpretation – first approach . 4 2.3 non-Local Thermodynamic Equilibrium – non-LTE . 6 2.4 Polarized light . 9 2.5 Atmospheric model . 9 3 A closer view – the dynamic atmosphere 11 3.1 Convection – granulation . 12 3.2 Waves ....................................... 13 3.3 Magnetic fields . 15 3.4 Chromosphere . 16 4 Conclusions 18 1 1 Introduction Importance of solar/stellar photosphere and chromosphere: • photosphere emits 99.99 % of energy generated in the solar interior by nuclear fusion, most of it in the visible spectral range • photosphere/chromosphere visible “skin” of solar “body” • structures in high atmosphere are rooted in photosphere/subphotosphere • dynamics/events in high atmosphere are caused by processes in deep (sub-)photospheric layers • chromosphere: onset of transport of mass, momentum, and energy to corona, solar wind, heliosphere, solar environment chromosphere = burning chamber for pre-heating non-static, non-equilibrium Extent of photosphere/chromosphere: • barometric formula (hydrostatic equilibrium): µg dp = −ρgdz and dp = −p dz (1) RT ⇒ p = p0 exp[−(z − z0)/Hp] (2) RT with “pressure scale height” Hp = µg ≈ 125 km (solar radius R ≈ 700 000 km) • extent: some 2 000 – 6 000 km (rugged) • “skin” of Sun In following: concepts, atmospheric model, dyanmic atmosphere 2 A coarse view – concepts 2.1 The data radiation (=b energy) solar output: measure radiation at Earth’s position, distance known 4 2 10 −2 −1 ⇒ F = σTeff, = L /(4πR ) = 6.3 × 10 erg cm s (3) ⇒ Teff, = 5780 K .
    [Show full text]
  • Sub-Surface Convection Zones in Hot Massive Stars and Their Observable Consequences
    A&A 499, 279–290 (2009) Astronomy DOI: 10.1051/0004-6361/200911643 & c ESO 2009 Astrophysics Sub-surface convection zones in hot massive stars and their observable consequences M. Cantiello1, N. Langer1,2, I. Brott1,A.deKoter1,3,S.N.Shore4,J.S.Vink5, A. Voegler1, D. J. Lennon6, and S.-C. Yoon7 1 Astronomical Institute, Utrecht University, Princetonplein 5, 3584 CC, Utrecht, The Netherlands e-mail: [email protected] 2 Argelander-Institut für Astronomie der Universität Bonn, Auf dem Hügel 71, 53121 Bonn, Germany 3 Astronomical Institute Anton Pannekoek, University of Amsterdam, Kruislaan 403, 1098 SJ, Amsterdam, The Netherlands 4 Dipartmento di Fisica “Enrico Fermi”, Università di Pisa, via Buonarroti 2, Pisa 56127 and INFN – Sezione di Pisa, Italy 5 Armagh Observatory, College Hill, Armagh, BT61 9DG, Northern Ireland, UK 6 Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA 7 Department of Astronomy & Astrophysics, University of California, Santa Cruz, High Street, Santa Cruz, CA 95064, USA Received 12 January 2009 / Accepted 3 March 2009 ABSTRACT Context. We study the convection zones in the outer envelope of hot massive stars which are caused by opacity peaks associated with iron and helium ionization. Aims. We determine the occurrence and properties of these convection zones as function of the stellar parameters. We then confront our results with observations of OB stars. Methods. A stellar evolution code is used to compute a grid of massive star models at different metallicities. In these models, the mixing length theory is used to characterize the envelope convection zones. Results. We find the iron convection zone (FeCZ) to be more prominent for lower surface gravity, higher luminosity and higher initial 3.2 3.9 4.2 metallicity.
    [Show full text]
  • Photons in the Radiative Zone: a Search for the Beginning Which Way Is Out? an A-Maz-Ing Model TEACHER GUIDE
    The Sun and Solar Wind: Photons in the Radiative Zone: A Search for the Beginning Which Way is Out? An A-Maz-ing Model TEACHER GUIDE BACKGROUND INFORMATION Students try to work their way out of a circular maze, thereby modeling the movement of a photon as it travels through the radiative zone of the sun. Classroom discussion after they complete the activity is focused on the Standard Solar Model and its importance in further scientific studies of the sun. STANDARDS ADDRESSED Grades 5-8 Science As Inquiry Understandings about scientific inquiry Science and Technology Understandings about science and technology Physical Science Properties and changes of properties in matter Transfer of energy History and Nature of Science Science as a human endeavor Nature of science and scientific knowledge History of science and historical perspectives Grades 9-12 Science As Inquiry Understandings about scientific inquiry Science and Technology Understandings about science and technology Earth and Space Science The origin and evolution of the Earth system Physical Science Properties and changes of properties in matter Transfer of energy History and Nature of Science Science as a human endeavor Nature of science and scientific knowledge History of science and historical perspectives GENESIS 1 MATERIALS For each student (or pair of students) Copy of Student Activity “Photons in the Radiative Zone: Which Way is Out?” A protractor A straight edge Copy of “Standard Model of the Sun” Copy of Student Text “Models in Science” PROCEDURE Alternative Strategy Tips 1. Before class, make copies of the Student Activity, “Photons in the Radiative Zone: Which Way is Out?” If you have not already done so, Prepare “black boxes” for each group of students by placing a small, makes copies of the Handout “Standard Model of the Sun” and Student unbreakable object in a box and Text, “Models in Science”.
    [Show full text]
  • The Solar Atmosphere
    The Solar Atmosphere Viggo H. Hansteen and Mats Carlsson Institute of Theoretical Astrophysics, University of Oslo PB 1029 Blindern, 0315 Oslo, Norway [email protected] 1 Introduction Looking at the solar photosphere, we see the top of the convection zone in the form of granulation: Hot gas rising from the solar interior as part of the energy transport process reaches a position where the opacity is no longer sufficient to prevent the escape of radiation. The gas expands, radiates, and cools and in so doing loses its buoyancy and descends. These motions, ultimately driven by the requirement that the energy generated by nuclear fusion in the Sun’s core be transported in the most efficient manner, represent a vast reservoir of “mechanical” energy flux. Looking closer, we see that granulation is not the only phenomenon visible at the solar surface. The quiet and semi-quiet photosphere is also threaded by magnetic fields that appear as bright points, as well as darker micro-pores and pores. These small scale magnetic structures are, while able to modify photospheric emission, subject to granular flows and seem to be passively carried by the convective motions. Convective flows are also known to generate the perturbations that drive solar oscillations. Oscillations, sound waves, with frequencies mainly in the band centered roughly at 3 mHz or 5 minutes are omnipresent in the solar photosphere and are collectively known as p-modes (‘p’ for pressure). These p-modes are a subject in their own right and studies of their properties have given solar physicists a unique tool in gathering information on solar struc- ture — the variation of the speed of sound cs, the rotation rate, and other important quantities — at depths far below those accessible through direct observations.
    [Show full text]
  • Dynamic Models of the Sun from the Convection Zone to the Chromosphere
    Convection in Astrophysics Proceedings IAU Symposium No. 239, 2006 c 2007 International Astronomical Union F. Kupka, I. W. Roxburgh & K. L. Chan, eds. doi:10.1017/S1743921307000105 Dynamic models of the sun from the convection zone to the chromosphere Sven Wedemeyer-B¨ohm1 1Kiepenheuer-Institut f¨ur Sonnenphysik, Sch¨oneckstr. 6, 79104 Freiburg, Germany email: [email protected] Abstract. The chromosphere in internetwork regions of the quiet Sun was regarded as a static and homogeneous layer for a long time. Thanks to advances in observations and numerical mod- elling, the wave nature of these atmospheric regions received increasing attention during the last decade. Recent three-dimensional radiation magnetohydrodynamic simulations with CO5BOLD feature the chromosphere of internetwork regions as a dynamic and intermittent phenomenon. It is a direct product of interacting waves that form a mesh-like pattern of hot shock fronts and cool post-shock regions. The waves are excited self-consistently at the top of the convec- tion zone. In the middle chromosphere above an average height of 1000 km, plasma beta gets larger than one and magnetic fields become more important. The model chromosphere exhibits a magnetic field that is much more homogeneous than in the layers below and evolves much faster. That includes fast propagating (MHD) waves. Further improvements of the simulations like time-dependent hydrogen ionisation are currently in progress. This class of models is capable of explaining apparently contradicting diagnostics such as carbon monoxide and UV emission at the same time. Keywords. Convection, hydrodynamics, MHD, radiative transfer, shock waves, Sun: chromo- sphere, granulation, photosphere, radio radiation, UV radiation.
    [Show full text]
  • Quiz 4 1)From the Center Outward, Which of the Following Lists The
    Quiz 4 1)From the center outward, which of the following lists the "layers" of the Sun in the correct order? A)core, convection zone, radiation zone, photosphere, chromosphere, corona B)core, radiation zone, convection zone, corona, chromosphere, photosphere * C)core, radiation zone, convection zone, photosphere, chromosphere, corona D)core, corona, radiation zone, convection zone, photosphere, chromosphere E)core, convection zone, radiation zone, corona, chromosphere, photosphere 2)Sunspots are cooler than the surrounding solar surface because A)they are regions where convection carries cooler material downward. B)magnetic fields trap ionized gases that absorb light. C)magnetic fields lift material from the surface of the Sun, cooling off the material faster. * D)strong magnetic fields slow convection and prevent hot plasma from entering the region. E)there is less fusion occurring there. 3)How do human-built nuclear power plants on Earth generate energy? A)chemical reactions B)converting kinetic energy into electricity * C)nuclear fission D)nuclear fusion E)converting gravitational potential energy into electricity 4)Studies of sunquakes, or helioseismology, have revealed that A)"sunquakes" are caused by similar processes that create earthquakes on the earth. B)neutrinos from the solar core reach the solar surface easily. C)the Sun vibrates only on the surface. D)the Sun generates energy by nuclear fusion. * E)our mathematical models of the solar interior are fairly accurate. 5)Why are neutrinos so difficult to detect? A)because they move at, or close to, the speed of light B)because there are so rare * C)because they rarely interact with matter D)because they have no mass E)We don't know: this is the essence of the solar neutrino problem.
    [Show full text]
  • Physics of the Solar Corona
    other process is needed to close the chain Fig. 2 — A possible sequence of processes related to the origin of the solar cycle and allow for self-excitation. It is tempting Region Physical Description Processes Related to Magnetic Field and to assume a certain flow pattern that twists Activity Cycle the toroidal field in such a way that the net effect is a meridional component of the Interface Zone Transition between ω(core) and Generation of toroidal field by differential rota­ ω(envelope)· Shear zone with pos- tion. Field is kept down by topological pump­ field. This leads to the formulation of the sible turbulence driven by instabi­ ing until it has grown strong enough to escape kinematical dynamo problem: Can we find lities due to differential rotation. owing to buoyancy and magnetic instabilities ; a velocity field v(x, t) in such a way, that formation of flux tubes. Helical flows produced the electrical field vxB/c drives a current j in that phase regenerate the poloidal field. which amplifies the magnetic field B? Convection Zone Nearly adiabatically stratified re­ Flux tubes are transported, shredded and refor­ In 1955, E.N. Parker presented a concep­ gion. Energy is transported almost med by convection and buoyancy. Small frag­ tual picture of a possible dynamo process entirely by convection. The veloci­ ments can be carried around by convection and based on the influence of rotation on con­ ty field may consist of large-scale appear later as small active regions randomly vective motions due to the Coriolis force. eddies as well as small-scale tur­ distributed over the surface.
    [Show full text]