-
The Stagger Code for Accurate and Efficient, Radiation-Coupled MHD Simulations
Authors:
Robert F. Stein,
Åke Nordlund,
Remo Collet,
Regner Trampedach
Abstract:
We describe the Stagger Code for simulations of magneto-hydrodynamic (MHD) systems. This is a modular code with a variety of physics modules that will let the user run simulations of deep stellar atmospheres, sunspot formation, stellar chromospheres and coronae, proto-stellar disks, star formation from giant molecular clouds and even galaxy formation. The Stagger Code is efficiently and highly par…
▽ More
We describe the Stagger Code for simulations of magneto-hydrodynamic (MHD) systems. This is a modular code with a variety of physics modules that will let the user run simulations of deep stellar atmospheres, sunspot formation, stellar chromospheres and coronae, proto-stellar disks, star formation from giant molecular clouds and even galaxy formation. The Stagger Code is efficiently and highly parallelizable, enabling such simulations with large ranges of both spatial and temporal scales. We, describe the methodology of the code, and present the most important of the physics modules, as well as its input and output variables. We show results of a number of standard MHD tests to enable comparison with other, similar codes. In addition, we provide an overview of tests that have been carried out against solar observations, ranging from spectral line shapes, spectral flux distribution, limb darkening, intensity and velocity distributions of granulation, to seismic power-spectra and the excitation of p modes. The Stagger Code has proven to be a high fidelity code with a large range of uses.
△ Less
Submitted 3 May, 2024;
originally announced May 2024.
-
Link between trees of fragmenting granules and deep downflows in MHD simulation
Authors:
T. Roudier,
J. M. Malherbe,
R. F. Stein,
Z. Frank
Abstract:
Trees of fragmenting granules (TFG) and associated flows are suspected to playa major role in the formation of the network in the quiet Sun. We investigate the counterparts, in terms of dynamics, of surface structures detectable by high resolution observations in deeper layers up to 15 Mm, which are only available from numerical simulations. The first aim is to demonstrate that TFG can be evidence…
▽ More
Trees of fragmenting granules (TFG) and associated flows are suspected to playa major role in the formation of the network in the quiet Sun. We investigate the counterparts, in terms of dynamics, of surface structures detectable by high resolution observations in deeper layers up to 15 Mm, which are only available from numerical simulations. The first aim is to demonstrate that TFG can be evidenced either from surface intensitites, vertical (Vz), or Doppler (Vdop) velocities. The second is to show that horizontal flows, which are derived from intensities or Vz/Vdop flows, are in good agreement, and that this is the case for observations and numerical simulations. The third objective is to apply this new Vz-based method to a 3D simulation to probe relationships between horizontal surface flows, TFG, and deep vertical motions. The TFG were detected after oscillation filtering of intensities or Vz/Vdop flows, using a segmentation and labelling technique. Surface horizontal flows were derived from local correlation tracking (LCT) and from intensities or Vz/Vdop flows. These methods were applied to Hinode observations, 2D surface results of a first simulation, and 3D Vz data of a second simulation. We find that TFG and horizontal surface flows (provided by the LCT) can be detected either from intensities or Vz/Vdop component, for high resolution observations and numerical simulations. We apply this method to a 3D run providing the Vz component in depth. This reveals a close relationship between surface TFG (5 Mm mesoscale) and vertical downflows 5 Mm below the surface. We suggest that the dynamics of TFG form larger scales (the 15-20 Mm supergranulation) associated with 15 Mm downflowing cells below the surface. The TFG and associated surface flows seem to be essential to understanding the formation and evolution of the network at the meso and supergranular scale.
△ Less
Submitted 10 January, 2019;
originally announced January 2019.
-
The asteroseismic surface effect from a grid of 3D convection simulations. I. Frequency shifts from convective expansion of stellar atmospheres
Authors:
Regner Trampedach,
Magnus J. Aarslev,
Günter Houdek,
Remo Collet,
Jørgen Christensen-Dalsgaard,
Robert F. Stein,
Martin Asplund
Abstract:
We analyse the effect on adiabatic stellar oscillation frequencies of replacing the near-surface layers in 1D stellar structure models with averaged 3D stellar surface convection simulations. The main difference is an expansion of the atmosphere by 3D convection, expected to explain a major part of the asteroseismic surface effect; a systematic overestimation of p-mode frequencies due to inadequat…
▽ More
We analyse the effect on adiabatic stellar oscillation frequencies of replacing the near-surface layers in 1D stellar structure models with averaged 3D stellar surface convection simulations. The main difference is an expansion of the atmosphere by 3D convection, expected to explain a major part of the asteroseismic surface effect; a systematic overestimation of p-mode frequencies due to inadequate surface physics.
We employ pairs of 1D stellar envelope models and 3D simulations from a previous calibration of the mixing-length parameter, alpha. That calibration constitutes the hitherto most consistent matching of 1D models to 3D simulations, ensuring that their differences are not spurious, but entirely due to the 3D nature of convection. The resulting frequency shift is identified as the structural part of the surface effect. The important, typically non-adiabatic, modal components of the surface effect are not included in the present analysis, but relegated to future papers.
Evaluating the structural surface effect at the frequency of maximum mode amplitude, $ν_{\rm max}$, we find shifts from $δν$=-0.8 microHz for giants at $\log g$=2.2 to -35 microHz for a ($T_{\rm eff}=6901$ K, $\log g$=4.29) dwarf. The fractional effect $δν(ν_{\rm max})/ν_{\rm max}$, ranges from -0.1% for a cool dwarf (4185 K, 4.74) to -6% for a warm giant (4962 K, 2.20).
△ Less
Submitted 8 November, 2016;
originally announced November 2016.
-
Improvements to stellar structure models, based on a grid of 3D convection simulations. II. Calibrating the mixing-length formulation
Authors:
Regner Trampedach,
Robert F. Stein,
Jørgen Christensen-Dalsgaard,
Åke Nordlund,
Martin Asplund
Abstract:
We perform a calibration of the mixing length of convection in stellar structure models against realistic 3D radiation-coupled hydrodynamics (RHD) simulations of convection in stellar surface layers, determining the adiabat deep in convective stellar envelopes.
The mixing-length parameter $α$ is calibrated by matching averages of the 3D simulations to 1D stellar envelope models, ensuring identic…
▽ More
We perform a calibration of the mixing length of convection in stellar structure models against realistic 3D radiation-coupled hydrodynamics (RHD) simulations of convection in stellar surface layers, determining the adiabat deep in convective stellar envelopes.
The mixing-length parameter $α$ is calibrated by matching averages of the 3D simulations to 1D stellar envelope models, ensuring identical atomic physics in the two cases. This is done for a previously published grid of solar-metallicity convection simulations, covering from 4200 K to 6900 K on the main sequence, and 4300-5000 K for giants with logg=2.2.
Our calibration results in an $α$ varying from 1.6 for the warmest dwarf, which is just cool enough to admit a convective envelope, and up to 2.05 for the coolest dwarfs in our grid. In between these is a triangular plateau of $α$ ~ 1.76. The Sun is located on this plateau and has seen little change during its evolution so far. When stars ascend the giant branch, they largely do so along tracks of constant $α$, with $α$ decreasing with increasing mass.
△ Less
Submitted 6 October, 2014;
originally announced October 2014.
-
Improvements to Stellar Structure Models, Based on a Grid of 3D Convection Simulations. I. $T(τ)$-Relations
Authors:
Regner Trampedach,
Robert F. Stein,
Jørgen Christensen-Dalsgaard,
Åke Nordlund,
Martin Asplund
Abstract:
Relations between temperature, T, and optical depth, tau, are often used for describing the photospheric transition from optically thick to optically thin in stellar structure models. We show that this is well justified, but also that currently used T(tau) relations are often inconsistent with their implementation. As an outer boundary condition on the system of stellar structure equations, T(tau)…
▽ More
Relations between temperature, T, and optical depth, tau, are often used for describing the photospheric transition from optically thick to optically thin in stellar structure models. We show that this is well justified, but also that currently used T(tau) relations are often inconsistent with their implementation. As an outer boundary condition on the system of stellar structure equations, T(tau) relations have an undue effect on the overall structure of stars. In this age of precision asteroseismology, we need to re-assess both the method for computing and for implementing T(tau) relations, and the assumptions they rest on. We develop a formulation for proper and consistent evaluation of T(tau) relations from arbitrary 1D or 3D stellar atmospheres, and for their implementation in stellar structure and evolution models. We extract radiative T(tau) relations, as described by our new formulation, from 3D simulations of convection in deep stellar atmospheres of late-type stars from dwarfs to giants. These simulations employ realistic opacities and equation of state, and account for line-blanketing. For comparison, we also extract T(tau) relations from 1D MARCS model atmospheres using the same formulation. T(tau)-relations from our grid of 3D convection simulations display a larger range of behaviours with surface gravity, compared with those of conventional theoretical 1D hydrostatic atmosphere models. Based on this, we recommend no longer to use scaled solar T(tau) relations. Files with T(tau) relations for our grid of simulations are made available to the community, together with routines for interpolating in this irregular grid. We also provide matching tables of atmospheric opacity, for consistent implementation in stellar structure models.
△ Less
Submitted 1 May, 2014; v1 submitted 1 May, 2014;
originally announced May 2014.
-
Models of solar surface dynamics: impact on eigenfrequencies and radius
Authors:
Laurent Piau,
Remo Collet,
Robert F. Stein,
Regner Trampedach,
Pierre Morel,
Sylvaine Turck-Chieze
Abstract:
We study the effects of different descriptions of the solar surface convection on the eigenfrequencies of p-modes. 1-D evolution calculations of the whole Sun and 3-D hydrodynamic and magnetohydrodynamic simulations of the current surface are performed. These calculations rely on realistic physics. Averaged stratifications of the 3-D simulations are introduced in the 1-D solar evolution or in the…
▽ More
We study the effects of different descriptions of the solar surface convection on the eigenfrequencies of p-modes. 1-D evolution calculations of the whole Sun and 3-D hydrodynamic and magnetohydrodynamic simulations of the current surface are performed. These calculations rely on realistic physics. Averaged stratifications of the 3-D simulations are introduced in the 1-D solar evolution or in the structure models. The eigenfrequencies obtained are compared to those of 1-D models relying on the usual phenomenologies of convection and to observations of the MDI instrument aboard SoHO. We also investigate how the magnetic activity could change the eigenfrequencies and the solar radius, assuming that, 3 Mm below the surface, the upgoing plasma advects a 1.2 kG horizontal field.
All models and observed eigenfrequencies are fairly close below 3 mHz. Above 3 mHz the eigenfrequencies of the phenomenological convection models are above the observed eigenfrequencies. The frequencies of the models based on the 3-D simulations are slightly below the observed frequencies. Their maximum deviation is ~ 3 μHz at 3 mHz but drops below 1 μHz at 4 mHz. Replacing the hydrodynamic by the magnetohydrodynamic simulation increases the eigenfrequencies. The shift is negligible below 2.2 mHz and then increases linearly with frequency to reach ~ 1.7 μHz at 4 mHz. The impact of the simulated activity is a 14 milliarcsecond shrinking of the solar layers near the optical depth unity.
△ Less
Submitted 27 September, 2013;
originally announced September 2013.
-
A Grid of 3D Stellar Atmosphere Models of Solar Metallicity: I. General Properties, Granulation and Atmospheric Expansion
Authors:
Regner Trampedach,
Martin Asplund,
Remo Collet,
Åke Nordlund,
Robert F. Stein
Abstract:
Present grids of stellar atmosphere models are the workhorses in interpreting stellar observations, and determining their fundamental parameters. These models rely on greatly simplified models of convection, however, lending less predictive power to such models of late type stars.
We present a grid of improved and more reliable stellar atmosphere models of late type stars, based on deep, 3D, con…
▽ More
Present grids of stellar atmosphere models are the workhorses in interpreting stellar observations, and determining their fundamental parameters. These models rely on greatly simplified models of convection, however, lending less predictive power to such models of late type stars.
We present a grid of improved and more reliable stellar atmosphere models of late type stars, based on deep, 3D, convective, stellar atmosphere simulations. This grid is to be used in general for interpreting observations, and improve stellar and asteroseismic modeling.
We solve the Navier Stokes equations in 3D and concurrent with the radiative transfer equation, for a range of atmospheric parameters, covering most of stellar evolution with convection at the surface. We emphasize use of the best available atomic physics for quantitative predictions and comparisons with observations.
We present granulation size, convective expansion of the acoustic cavity, asymptotic adiabat, as function of atmospheric parameters. These and other results are also available in electronic form.
△ Less
Submitted 7 March, 2013;
originally announced March 2013.
-
The Stagger-grid: A Grid of 3D Stellar Atmosphere Models - I. Methods and General Properties
Authors:
Z. Magic,
R. Collet,
M. Asplund,
R. Trampedach,
W. Hayek,
A. Chiavassa,
R. F. Stein,
Å. Nordlund
Abstract:
We present the Stagger-grid, a comprehensive grid of time-dependent, 3D hydrodynamic model atmospheres for late-type stars with realistic treatment of radiative transfer, covering a wide range in stellar parameters. This grid of 3D models is intended for various applications like stellar spectroscopy, asteroseismology and the study of stellar convection. In this introductory paper, we describe the…
▽ More
We present the Stagger-grid, a comprehensive grid of time-dependent, 3D hydrodynamic model atmospheres for late-type stars with realistic treatment of radiative transfer, covering a wide range in stellar parameters. This grid of 3D models is intended for various applications like stellar spectroscopy, asteroseismology and the study of stellar convection. In this introductory paper, we describe the methods used for the computation of the grid and discuss the general properties of the 3D models as well as their temporal and spatial averages (<3D>). All our models were generated with the Stagger-code, using realistic input physics for the equation of state (EOS) and for continuous and line opacities. Our ~220 grid models range in Teff from 4000 to 7000K in steps of 500K, in log g from 1.5 to 5.0 in steps of 0.5 dex, and [Fe/H] from -4.0 to +0.5 in steps of 0.5 and 1.0 dex. We find a tight scaling relation between the vertical velocity and the surface entropy jump, which itself correlates with the constant entropy value of the adiabatic convection zone. The range in intensity contrast is enhanced at lower metallicity. The granule size correlates closely with the pressure scale height sampled at the depth of maximum velocity. We compare the <3D> models with widely applied 1D models, as well as with theoretical 1D hydrostatic models generated with the same EOS and opacity tables as the 3D models, in order to isolate the effects of using self-consistent and hydrodynamic modeling of convection, rather than the classical mixing length theory approach. For the first time, we are able to quantify systematically over a broad range of stellar parameters the uncertainties of 1D models arising from the simplified treatment of physics, in particular convective energy transport. In agreement with previous findings, we find that the differences can be significant, especially for metal-poor stars.
△ Less
Submitted 12 July, 2013; v1 submitted 11 February, 2013;
originally announced February 2013.
-
On the Formation of Active Regions
Authors:
Robert F. Stein,
Åke Nordlund
Abstract:
Magneto-convection can produce an active region without an initial coherent flux tube. A simulation was performed where uniform, untwisted, horizontal magnetic field of 1 kG strenght was advected into the bottom of a computational domain 48 Mm wide by 20 Mm deep. The up and down convective motions produce a hierarchy of magnetic loops with a wide range of scales, with smaller loops riding "piggy b…
▽ More
Magneto-convection can produce an active region without an initial coherent flux tube. A simulation was performed where uniform, untwisted, horizontal magnetic field of 1 kG strenght was advected into the bottom of a computational domain 48 Mm wide by 20 Mm deep. The up and down convective motions produce a hierarchy of magnetic loops with a wide range of scales, with smaller loops riding "piggy back" in a serpentine fashion on larger loops. When a large loop approaches the surface it produces an small active region with a compact leading spot and more diffuse following spots.
△ Less
Submitted 17 July, 2012;
originally announced July 2012.
-
Turbulent Diffusion in the Photosphere as Derived from Photospheric Bright Point Motion
Authors:
V. I. Abramenko,
V. Carbone,
V. Yurchyshyn,
P. R. Goode,
R. F. Stein,
F. Lepreti,
V. Capparelli,
A. Vecchio
Abstract:
On the basis of observations of solar granulation obtained with the New Solar Telescope (NST) of Big Bear Solar Observatory, we explored proper motion of bright points (BPs) in a quiet sun area, a coronal hole, and an active region plage. We automatically detected and traced bright points (BPs) and derived their mean-squared displacements as a function of time (starting from the appearance of each…
▽ More
On the basis of observations of solar granulation obtained with the New Solar Telescope (NST) of Big Bear Solar Observatory, we explored proper motion of bright points (BPs) in a quiet sun area, a coronal hole, and an active region plage. We automatically detected and traced bright points (BPs) and derived their mean-squared displacements as a function of time (starting from the appearance of each BP) for all available time intervals. In all three magnetic environments, we found the presence of a super-diffusion regime, which is the most pronounced inside the time interval of 10-300 seconds. Super-diffusion, measured via the spectral index, $γ$, which is the slope of the mean-squared displacement spectrum, increases from the plage area ($γ=1.48$) to the quiet sun area ($γ=1.53$) to the coronal hole ($γ=1.67$). We also found that the coefficient of turbulent diffusion changes in direct proportion to both temporal and spatial scales. For the minimum spatial scale (22 km) and minimum time scale (10 sec), it is 22 and 19 km$^{2}$ s$^{-1}$ for the coronal hole and the quiet sun area, respectively, whereas for the plage area it is about 12 km$^{2}$ s$^{-1}$ for the minimum time scale of 15 seconds. We applied our BP tracking code to 3D MHD model data of solar convection (Stein et al. 2007) and found the super-diffusion with $γ=1.45$. An expression for the turbulent diffusion coefficient as a function of scales and $γ$ is obtained.
△ Less
Submitted 18 November, 2011;
originally announced November 2011.
-
The Mass Mixing Length in Convective Stellar Envelopes
Authors:
Regner Trampedach,
Robert F. Stein
Abstract:
The scale length over which convection mixes mass in a star can be calculated as the inverse of the vertical derivative of the unidirectional (up or down) mass flux. This is related to the mixing length in the mixing length theory of stellar convection. We give the ratio of mass mixing length to pressure scale height for a grid of 3D surface convection simulations, covering from 4300\,K to 6900\,K…
▽ More
The scale length over which convection mixes mass in a star can be calculated as the inverse of the vertical derivative of the unidirectional (up or down) mass flux. This is related to the mixing length in the mixing length theory of stellar convection. We give the ratio of mass mixing length to pressure scale height for a grid of 3D surface convection simulations, covering from 4300\,K to 6900\,K on the main-sequence, and up to giants at $\log g = 2.2$, all for solar composition. These simulations also confirm what is already known from solar simulations, that convection doesn't proceed by discrete convective elements, but rather as a continuous, slow, smooth, warm upflow and turbulent, entropy deficient, fast down drafts. This convective topology also results in mixing on a scale as that of the classic mixing length formulation, and is simply a consequence of mass conservation on flows in a stratified atmosphere.
△ Less
Submitted 5 February, 2011;
originally announced February 2011.
-
Emerging Flux Simulations and Proto-Active Regions
Authors:
Robert F. Stein,
Anders Lagerfjärd,
Åke Nordlund,
Dali Georgobiani
Abstract:
The emergence of minimally structured (uniform and horizontal) magnetic field from a depth of 20 Mm has been simulated. The field emerges first in a mixed polarity pepper and salt pattern, but then collects into separate, unipolar concentrations and produces pores. The field strength was then artificially increased to produce spot-like structures. The field strength at continuum optical depth unit…
▽ More
The emergence of minimally structured (uniform and horizontal) magnetic field from a depth of 20 Mm has been simulated. The field emerges first in a mixed polarity pepper and salt pattern, but then collects into separate, unipolar concentrations and produces pores. The field strength was then artificially increased to produce spot-like structures. The field strength at continuum optical depth unity peaks at 1 kG, with a maximum of 4 kG. Where the vertical field is strong, the spots persist (at present an hour of solar time has been simulated). Where the field is weak, the spot gets filled in and disappears. Stokes profiles have been calculated and processed with the Hinode annular mtf, the slit diffraction and frequency smoothing. These data are available at steinr.pa.msu.edu/~bob/stokes.
△ Less
Submitted 4 February, 2011;
originally announced February 2011.
-
Solar Flux Emergence Simulations
Authors:
R. F. Stein,
A. Lagerfjärd,
Å. Nordlund,
D. Georgobiani
Abstract:
We simulate the rise through the upper convection zone and emergence through the solar surface of initially uniform, untwisted, horizontal magnetic flux with the same entropy as the non-magnetic plasma that is advected into a domain 48 Mm wide from from 20 Mm deep. The magnetic field is advected upward by the diverging upflows and pulled down in the downdrafts, which produces a hierarchy of loop…
▽ More
We simulate the rise through the upper convection zone and emergence through the solar surface of initially uniform, untwisted, horizontal magnetic flux with the same entropy as the non-magnetic plasma that is advected into a domain 48 Mm wide from from 20 Mm deep. The magnetic field is advected upward by the diverging upflows and pulled down in the downdrafts, which produces a hierarchy of loop like structures of increasingly smaller scale as the surface is approached. There are significant differences between the behavior of fields of 10 kG and 20 or 40 kG strength at 20 Mm depth. The 10 kG fields have little effect on the convective flows and show little magnetic buoyancy effects, reaching the surface in the typical fluid rise time from 20 Mm depth of 32 hours. 20 and 40 kG fields significantly modify the convective flows, leading to long thin cells of ascending fluid aligned with the magnetic field and their magnetic buoyancy makes them rise to the surface faster than the fluid rise time. The 20 kG field produces a large scale magnetic loop that as it emerges through the surface leads to the formation of a bipolar pore-like structure.
△ Less
Submitted 24 December, 2009;
originally announced December 2009.
-
Impact of the physical processes in the modeling of HD49933
Authors:
L. Piau,
S. Turck-Chieze,
V. Duez,
R. F. Stein
Abstract:
Context : On its asteroseismic side, the initial run of CoRoT was partly devoted to the solar like star HD49933.The eigenmodes of this F dwarf have been observed with unprecedented accuracy.
Aims : We investigate quantitatively the impact of changes in the modeling parameters like mass and composition. More importantly we investigate how a sophisticated physics affects the seismological pictur…
▽ More
Context : On its asteroseismic side, the initial run of CoRoT was partly devoted to the solar like star HD49933.The eigenmodes of this F dwarf have been observed with unprecedented accuracy.
Aims : We investigate quantitatively the impact of changes in the modeling parameters like mass and composition. More importantly we investigate how a sophisticated physics affects the seismological picture of HD49933. We consider the effects of diffusion, rotation and the changes in convection efficiency.
Methods : We use the CESAM stellar evolution code coupled to the ADIPLS adiabatic pulsation package to build secular models and their associated oscillation frequencies. We also exploited the hydrodynamical code STAGGER to perform surface convection calculations. The seismic variables used in this work are : the large frequency separation, the derivative of the surface phase shift,and the eigenfrequencies $\rm ν_{\ell=0,n=14}$ and $\rm ν_{\ell=0,n=27}$.
Results : Mass and uncertainties on the composition have much larger impacts on the seismic variables we consider than the rotation. The derivative of the surface phase shift is a promising variable for the determination of the helium content. The seismological variables of HD49933 are sensitive to the assumed solar composition and also to the presence of diffusion in the models.
△ Less
Submitted 24 July, 2009;
originally announced July 2009.
-
Supergranulation Scale Connection Simulations
Authors:
R. F. Stein,
A. Nordlund,
D. Georgobiani,
D. Benson,
W. Schaffenberger
Abstract:
Results of realistic simulations of solar surface convection on the scale of supergranules (96 Mm wide by 20 Mm deep) are presented. The simulations cover only 10% of the geometric depth of the solar convection zone, but half its pressure scale heights. They include the hydrogen, first and most of the second helium ionization zones. The horizontal velocity spectrum is a power law and the horizon…
▽ More
Results of realistic simulations of solar surface convection on the scale of supergranules (96 Mm wide by 20 Mm deep) are presented. The simulations cover only 10% of the geometric depth of the solar convection zone, but half its pressure scale heights. They include the hydrogen, first and most of the second helium ionization zones. The horizontal velocity spectrum is a power law and the horizontal size of the dominant convective cells increases with increasing depth. Convection is driven by buoyancy work which is largest close to the surface, but significant over the entire domain. Close to the surface buoyancy driving is balanced by the divergence of the kinetic energy flux, but deeper down it is balanced by dissipation. The damping length of the turbulent kinetic energy is 4 pressure scale heights. The mass mixing length is 1.8 scale heights. Two thirds of the area is upflowing fluid except very close to the surface. The internal (ionization) energy flux is the largest contributor to the convective flux for temperatures less than 40,000 K and the thermal energy flux is the largest contributor at higher temperatures. This data set is useful for validating local helioseismic inversion methods. Sixteen hours of data are available as four hour averages, with two hour cadence, at steinr.msu.edu/~bob/96averages, as idl save files. The variables stored are the density, temperature, sound speed, and three velocity components. In addition, the three velocity components at 200 km above mean continuum optical depth unity are available at 30 sec. cadence.
△ Less
Submitted 4 November, 2008;
originally announced November 2008.
-
Helioseismic Holography of Simulated Solar Convection and Prospects for the Detection of Small-Scale Subsurface Flows
Authors:
D. C. Braun,
A. C. Birch,
D. Benson,
R. F. Stein,
Å. Nordlund
Abstract:
We perform helioseismic holography on realistic solar convection simulations and compare the observed travel-time perturbations with the expected travel times from the horizontal flows in the simulations computed from forward models under the assumption of the Born approximation. We demonstrate reasonable agreement between the observed and model travel times which reinforces the validity of heli…
▽ More
We perform helioseismic holography on realistic solar convection simulations and compare the observed travel-time perturbations with the expected travel times from the horizontal flows in the simulations computed from forward models under the assumption of the Born approximation. We demonstrate reasonable agreement between the observed and model travel times which reinforces the validity of helioseismic holography in the detection of subsurface horizontal flows. From the variation of the signal-to-noise ratio with depth, we conclude that the helioseismic detection of individual flow structures with spatial scales of supergranulation or smaller is not possible for depths below about 5 Mm below the surface over time scales less than a day. Approximately half of the observed signal originates within the first 2 Mm below the surface. A consequence of this is a rapid decrease (and reversal in some cases) of the travel-time perturbations with depth due to the contribution to the measurements of oppositely directed surface flows in neighboring convective cells. This confirms an earlier interpretation of similar effects reported from observations.
△ Less
Submitted 1 August, 2007;
originally announced August 2007.
-
Validation of Time-Distance Helioseismology by Use of Realistic Simulations of Solar Convection
Authors:
Junwei Zhao,
Dali Georgobiani,
Alexander G. Kosovichev,
David Benson,
Robert F. Stein,
Åke Nordlund
Abstract:
Recent progress in realistic simulations of solar convection have given us an unprecedented opportunity to evaluate the robustness of solar interior structures and dynamics obtained by methods of local helioseismology. We present results of testing the time-distance method using realistic simulations. By computing acoustic wave propagation time and distance relations for different depths of the…
▽ More
Recent progress in realistic simulations of solar convection have given us an unprecedented opportunity to evaluate the robustness of solar interior structures and dynamics obtained by methods of local helioseismology. We present results of testing the time-distance method using realistic simulations. By computing acoustic wave propagation time and distance relations for different depths of the simulated data, we confirm that acoustic waves propagate into the interior and then turn back to the photosphere. This demonstrates that in the numerical simulations properties of acoustic waves (p-modes) are similar to the solar conditions, and that these properties can be analyzed by the time-distance technique. For the surface gravity waves (f-mode), we calculate perturbations of their travel times, caused by localized downdrafts, and demonstrate that the spatial pattern of these perturbations (representing so-called sensitivity kernels) is similar to the patterns obtained from the real Sun, displaying characteristic hyperbolic structures. We then test the time-distance measurements and inversions by calculating acoustic travel times from a sequence of vertical velocities at the photosphere of the simulated data, and inferring a mean 3D flow fields by performing inversion based on the ray approximation. The inverted horizontal flow fields agree very well with the simulated data in subsurface areas up to 3 Mm deep, but differ in deeper areas. Due to the cross-talk effects between the horizontal divergence and downward flows, the inverted vertical velocities are significantly different from the mean convection velocities of the simulation dataset.
△ Less
Submitted 19 December, 2006;
originally announced December 2006.
-
Excitation of solar-like oscillations across the HR diagram
Authors:
R. Samadi,
D. Georgobiani,
R. Trampedach R.,
M. -J. Goupil,
R. F. Stein,
A. Nordlund
Abstract:
We extend semi-analytical computations of excitation rates for solar oscillation modes to those of other solar-like oscillating stars to compare them with recent observations. Numerical 3D simulations of surface convective zones of several solar-type oscillating stars are used to characterize the turbulent spectra as well as to constrain the convective velocities and turbulent entropy fluctuatio…
▽ More
We extend semi-analytical computations of excitation rates for solar oscillation modes to those of other solar-like oscillating stars to compare them with recent observations. Numerical 3D simulations of surface convective zones of several solar-type oscillating stars are used to characterize the turbulent spectra as well as to constrain the convective velocities and turbulent entropy fluctuations in the uppermost part of the convective zone of such stars. These constraints, coupled with a theoretical model for stochastic excitation, provide the rate 'P' at which energy is injected into the p-modes by turbulent convection. These energy rates are compared with those derived directly from the 3D simulations. The excitation rates obtained from the 3D simulations are systematically lower than those computed from the semi-analytical excitation model. We find that Pmax, the excitation rate maximum, scales as (L/M)^s where s is the slope of the power law and L and M are the mass and luminosity of the 1D stellar model built consistently with the associated 3D simulation. The slope is found to depend significantly on the adopted form of the eddy time-correlation ; using a Lorentzian form results in s=2.6, whereas a Gaussian one gives s=3.1. Finally, values of Vmax, the maximum in the mode velocity, are estimated from the computed power laws for Pmax and we find that Vmax increases as (L/M)^sv. Comparisons with the currently available ground-based observations show that the computations assuming a Lorentzian eddy time-correlation yield a slope, sv, closer to the observed one than the slope obtained when assuming a Gaussian. We show that the spatial resolution of the 3D simulations must be high enough to obtain accurate computed energy rates.
△ Less
Submitted 24 November, 2006;
originally announced November 2006.
-
Velocities measured in small scale solar magnetic elements
Authors:
Øystein Langangen,
Mats Carlsson,
Luc Rouppe van der Voort,
R. F. Stein
Abstract:
We have obtained high resolution spectrograms of small scale magnetic structures with the Swedish 1-m Solar Telescope. We present Doppler measurements at $0\farcs{2}$ spatial resolution of bright points, ribbons and flowers and their immediate surroundings, in the C {\small{I}} 5380.3 Å line (formed in the deep photosphere) and the two Fe {\small{I}} lines at 5379.6 Å and 5386.3 Å. The velocity…
▽ More
We have obtained high resolution spectrograms of small scale magnetic structures with the Swedish 1-m Solar Telescope. We present Doppler measurements at $0\farcs{2}$ spatial resolution of bright points, ribbons and flowers and their immediate surroundings, in the C {\small{I}} 5380.3 Å line (formed in the deep photosphere) and the two Fe {\small{I}} lines at 5379.6 Å and 5386.3 Å. The velocity inside the flowers and ribbons are measured to be almost zero, while we observe downflows at the edges. These downflows are increasing with decreasing height. We also analyze realistic magneto-convective simulations to obtain a better understanding of the interpretation of the observed signal. We calculate how the Doppler signal depends on the velocity field in various structures. Both the smearing effect of the non-negligible width of this velocity response function along the line of sight and of the smearing from the telescope and atmospheric point spread function are discussed. These studies lead us to the conclusion that the velocity inside the magnetic elements are really upflow of the order 1--2 km s${}^{-1}$ while the downflows at the edges really are much stronger than observed, of the order 1.5--3.3 km s${}^{-1}$.
△ Less
Submitted 23 November, 2006;
originally announced November 2006.
-
Local helioseismology and correlation tracking analysis of surface structures in realistic simulations of solar convection
Authors:
Dali Georgobiani,
Junwei Zhao,
Alexander Kosovichev,
David Benson,
Robert F. Stein,
Åke Nordlund
Abstract:
We apply time-distance helioseismology, local correlation tracking and Fourier spatial-temporal filtering methods to realistic supergranule scale simulations of solar convection and compare the results with high-resolution observations from the SOHO Michelson Doppler Imager (MDI). Our objective is to investigate the surface and sub-surface convective structures and test helioseismic measurements…
▽ More
We apply time-distance helioseismology, local correlation tracking and Fourier spatial-temporal filtering methods to realistic supergranule scale simulations of solar convection and compare the results with high-resolution observations from the SOHO Michelson Doppler Imager (MDI). Our objective is to investigate the surface and sub-surface convective structures and test helioseismic measurements. The size and grid of the computational domain are sufficient to resolve various convective scales from granulation to supergranulation. The spatial velocity spectrum is approximately a power law for scales larger than granules, with a continuous decrease in velocity amplitude with increasing size. Aside from granulation no special scales exist, although a small enhancement in power at supergranulation scales can be seen. We calculate the time-distance diagram for f- and p-modes and show that it is consistent with the SOHO/MDI observations. From the simulation data we calculate travel time maps for surface gravity waves (f-mode). We also apply correlation tracking to the simulated vertical velocity in the photosphere to calculate the corresponding horizontal flows. We compare both of these to the actual large-scale (filtered) simulation velocities. All three methods reveal similar large scale convective patterns and provide an initial test of time-distance methods.
△ Less
Submitted 31 January, 2007; v1 submitted 9 August, 2006;
originally announced August 2006.
-
Simulation of quiet Sun waves in the Ca II infrared triplet
Authors:
A. Pietarila,
H. Socas-Navarro,
T. Bogdan,
M. Carlsson,
R. F. Stein
Abstract:
The Ca II infrared triplet around 8540 A is a good candidate for observing chromospheric magnetism. We show results from combining a radiation hydrodynamic simulation with a Stokes synthesis code. The simulation shows interesting time-varying behavior of the Stokes V profiles as waves propagate through the lines' formation regions. Disappearing and reappearing lobes in the Stokes V profiles as w…
▽ More
The Ca II infrared triplet around 8540 A is a good candidate for observing chromospheric magnetism. We show results from combining a radiation hydrodynamic simulation with a Stokes synthesis code. The simulation shows interesting time-varying behavior of the Stokes V profiles as waves propagate through the lines' formation regions. Disappearing and reappearing lobes in the Stokes V profiles as well as profile asymmetries are closely related to the atmospheric velocity gradients.
△ Less
Submitted 26 October, 2005;
originally announced October 2005.
-
Effect of the radiative background flux in convection
Authors:
A. Brandenburg,
K. L. Chan,
A. Nordlund,
R. F. Stein
Abstract:
Numerical simulations of turbulent stratified convection are used to study models with approximately the same convective flux, but different radiative fluxes. As the radiative flux is decreased, for constant convective flux: the entropy jump at the top of the convection zone becomes steeper, the temperature fluctuations increase and the velocity fluctuations decrease in magnitude, and the distan…
▽ More
Numerical simulations of turbulent stratified convection are used to study models with approximately the same convective flux, but different radiative fluxes. As the radiative flux is decreased, for constant convective flux: the entropy jump at the top of the convection zone becomes steeper, the temperature fluctuations increase and the velocity fluctuations decrease in magnitude, and the distance that low entropy fluid from the surface can penetrate increases. Velocity and temperature fluctuations follow mixing length scaling laws.
△ Less
Submitted 18 August, 2005;
originally announced August 2005.
-
Spectrum and amplitudes of internal gravity waves excited by penetrative convection in solar-type stars
Authors:
B. Dintrans,
A. Brandenburg,
A. Nordlund,
R. F. Stein
Abstract:
The excitation of internal gravity waves by penetrative convective plumes is investigated using 2-D direct simulations of compressible convection. The wave generation is quantitatively studied from the linear response of the radiative zone to the plumes penetration, using projections onto the g-modes solutions of the associated linear eigenvalue problem for the perturbations. This allows an accu…
▽ More
The excitation of internal gravity waves by penetrative convective plumes is investigated using 2-D direct simulations of compressible convection. The wave generation is quantitatively studied from the linear response of the radiative zone to the plumes penetration, using projections onto the g-modes solutions of the associated linear eigenvalue problem for the perturbations. This allows an accurate determination of both the spectrum and amplitudes of the stochastically excited modes. Using time-frequency diagrams of the mode amplitudes, we then show that the lifetime of a mode is around twice its period and that during times of significant excitation up to 40% of the total kinetic energy may be contained into g-modes.
△ Less
Submitted 11 April, 2005; v1 submitted 7 February, 2005;
originally announced February 2005.
-
Excitation rates of p modes: mass luminosity relation across the HR diagram
Authors:
R. Samadi,
D. Georgobiani,
R. Trampedach,
M. J. Goupil,
R. F. Stein,
A. Nordlund
Abstract:
We compute the rates P at which energy is injected into the p modes for a set of 3D simulations of outer layers of stars.
We found that Pmax - the maximum in P - scales as (L/M)^s where s is the slope of the power law, L and M are the luminosity and the mass of the 1D stellar models associated with the simulations. The slope is found to depend significantly on the adopted representation for the…
▽ More
We compute the rates P at which energy is injected into the p modes for a set of 3D simulations of outer layers of stars.
We found that Pmax - the maximum in P - scales as (L/M)^s where s is the slope of the power law, L and M are the luminosity and the mass of the 1D stellar models associated with the simulations. The slope is found to depend significantly on the adopted representation for the turbulent eddy-time correlation function, chi_k. According to the expected performances of COROT, it will likely be possible to measure Pmax as a function of L/M and to constrain the properties of stellar turbulence as the turbulent eddy time-correlation.
△ Less
Submitted 2 October, 2004;
originally announced October 2004.
-
Observational manifestations of solar magneto-convection -- center-to-limb variation
Authors:
Mats Carlsson,
Robert F. Stein,
Ake Nordlund,
Goran Scharmer
Abstract:
We present the first center-to-limb G-band images synthesized from high resolution simulations of solar magneto-convection. Towards the limb the simulations show "hilly" granulation with dark bands on the far side, bright granulation walls and striated faculae, similar to observations. At disk center G-band bright points are flanked by dark lanes. The increased brightness in magnetic elements is…
▽ More
We present the first center-to-limb G-band images synthesized from high resolution simulations of solar magneto-convection. Towards the limb the simulations show "hilly" granulation with dark bands on the far side, bright granulation walls and striated faculae, similar to observations. At disk center G-band bright points are flanked by dark lanes. The increased brightness in magnetic elements is due to their lower density compared with the surrounding intergranular medium. One thus sees deeper layers where the temperature is higher. At a given geometric height, the magnetic elements are cooler than the surrounding medium. In the G-band, the contrast is further increased by the destruction of CH in the low density magnetic elements. The optical depth unity surface is very corrugated. Bright granules have their continuum optical depth unity 80 km above the mean surface, the magnetic elements 200-300 km below. The horizontal temperature gradient is especially large next to flux concentrations. When viewed at an angle, the deep magnetic elements optical surface is hidden by the granules and the bright points are no longer visible, except where the "magnetic valleys" are aligned with the line of sight. Towards the limb, the low density in the strong magnetic elements causes unit line-of-sight optical depth to occur deeper in the granule walls behind than for rays not going through magnetic elements and variations in the field strength produce a striated appearance in the bright granule walls.
△ Less
Submitted 7 June, 2004;
originally announced June 2004.
-
Stochastic excitation of gravity waves by overshooting convection in solar-type stars
Authors:
Boris Dintrans,
Axel Brandenburg,
Ake Nordlund,
R. F. Stein
Abstract:
The excitation of gravity waves by penetrative convective plumes is investigated using 2D direct simulations of compressible convection. The oscillation field is measured by a new technique based on the projection of our simulation data onto the theoretical g-modes solutions of the associated linear eigenvalue problem. This allows us to determine both the excited modes and their corresponding am…
▽ More
The excitation of gravity waves by penetrative convective plumes is investigated using 2D direct simulations of compressible convection. The oscillation field is measured by a new technique based on the projection of our simulation data onto the theoretical g-modes solutions of the associated linear eigenvalue problem. This allows us to determine both the excited modes and their corresponding amplitudes accurately.
△ Less
Submitted 3 March, 2004;
originally announced March 2004.
-
Numerical 3D constraints on convective eddy time-correlations : consequences for stochastic excitation of solar p modes
Authors:
R. Samadi,
A. Nordlund,
R. F. Stein,
M. J. Goupil,
I. Roxburgh
Abstract:
A 3D simulation of the upper part of the solar convective zone is used to obtain information on the frequency component, chi_k, of the correlation product of the turbulent velocity field. This component plays an important role in the stochastic excitation of acoustic oscillations. A time analysis of the solar simulation shows that a gaussian function does not correctly reproduce the nu-dependenc…
▽ More
A 3D simulation of the upper part of the solar convective zone is used to obtain information on the frequency component, chi_k, of the correlation product of the turbulent velocity field. This component plays an important role in the stochastic excitation of acoustic oscillations. A time analysis of the solar simulation shows that a gaussian function does not correctly reproduce the nu-dependency of chi_k inferred from the 3D simuation in the frequency range where the acoustic energy injected into the solar p modes is important (nu ~ 2 - 4 mHz). The nu-dependency of chi_k is fitted with different analytical functions which can then conveniently be used to compute the acoustic energy supply rate P injected into the solar radial oscillations. With constraints from a 3D simulation, adjustment of free parameters to solar data is no longer necessary and is not performed here. The result is compared with solar seismic data. Computed values of P obtained with the analytical function which fits best chi_k are found ~ 2.7 times larger than those obtained with the gaussian model and reproduce better the solar seismic observations. This non-gaussian description also leads to a Reynolds stress contribution of the same order as the one arising from the advection of the turbulent fluctuations of entropy by the turbulent motions. Some discrepancy between observed and computed P values still exist at high frequency and possible causes for this discrepancy are discussed.
△ Less
Submitted 28 April, 2003; v1 submitted 25 April, 2003;
originally announced April 2003.
-
Numerical constraints on the model of stochastic excitation of solar-type oscillations
Authors:
R. Samadi,
A. Nordlund,
R. F. Stein,
M. J. Goupil,
I. Roxburgh
Abstract:
Analyses of a 3D simulation of the upper layers of a solar convective envelope provide constraints on the physical quantities which enter the theoretical formulation of a stochastic excitation model of solar p modes, for instance the convective velocities and the turbulent kinetic energy spectrum. These constraints are then used to compute the acoustic excitation rate for solar p modes, P. The r…
▽ More
Analyses of a 3D simulation of the upper layers of a solar convective envelope provide constraints on the physical quantities which enter the theoretical formulation of a stochastic excitation model of solar p modes, for instance the convective velocities and the turbulent kinetic energy spectrum. These constraints are then used to compute the acoustic excitation rate for solar p modes, P. The resulting values are found ~5 times larger than the values resulting from a computation in which convective velocities and entropy fluctuations are obtained with a 1D solar envelope model built with the time-dependent, nonlocal Gough (1977) extension of the mixing length formulation for convection (GMLT). This difference is mainly due to the assumed mean anisotropy properties of the velocity field in the excitation region. The 3D simulation suggests much larger horizontal velocities compared to vertical ones than in the 1D GMLT solar model. The values of P obtained with the 3D simulation constraints however are still too small compared with the values inferred from solar observations. Improvements in the description of the turbulent kinetic energy spectrum and its depth dependence yield further increased theoretical values of P which bring them closer to the observations. It is also found that the source of excitation arising from the advection of the turbulent fluctuations of entropy by the turbulent movements contributes ~ 65-75 % to the excitation and therefore remains dominant over the Reynolds stress contribution. The derived theoretical values of P obtained with the 3D simulation constraints remain smaller by a factor ~3 compared with the solar observations. This shows that the stochastic excitation model still needs to be improved.
△ Less
Submitted 10 March, 2003;
originally announced March 2003.
-
Consequences of the non gaussian character of the stochastic excitation for solar-type oscillations
Authors:
R. Samadi,
A. Nordlund,
R. F. Stein,
M. J. Goupil,
I. Roxburgh
Abstract:
Our recent study, based on a 3D numerical simulation, suggested that the stochastic excitation of solar p-modes has a non gaussian nature and a non gaussian model was proposed which results in a good agreement between the computed rate $P$ at which p-modes are excited and the solar seismic observational constraints.
In the present work, we study some consequences of this non-gaussian model on…
▽ More
Our recent study, based on a 3D numerical simulation, suggested that the stochastic excitation of solar p-modes has a non gaussian nature and a non gaussian model was proposed which results in a good agreement between the computed rate $P$ at which p-modes are excited and the solar seismic observational constraints.
In the present work, we study some consequences of this non-gaussian model on the rate P at which p-modes are excited in intermediate mass stars (1 Mo < M < 2 Mo). The non gaussian model changes substantially the spectrum of P at high frequency compared with the commonly assumed gaussian model. The largest effects on P are found for stars in the mass range ~ 1.5-1.6 Mo. They are found large enough that observatins in solar-like oscillating stars observed with space based seismology missions (e.g. COROT) will be able to confirm or reject the gaussian nature of stochastic excitation.
△ Less
Submitted 1 October, 2002;
originally announced October 2002.
-
Solar Surface Magneto-Convection
Authors:
Robert F. Stein,
Aake Nordlund
Abstract:
Magneto-convection simulations on meso-granule and granule scales near the solar surface are used to study small scale dynamo activity, the emergence and disappearance of magnetic flux tubes, and the formation and evolution of micropores.
From weak seed fields, convective motions produce highly intermittent magnetic fields in the intergranular lanes which collect over the boundaries of the und…
▽ More
Magneto-convection simulations on meso-granule and granule scales near the solar surface are used to study small scale dynamo activity, the emergence and disappearance of magnetic flux tubes, and the formation and evolution of micropores.
From weak seed fields, convective motions produce highly intermittent magnetic fields in the intergranular lanes which collect over the boundaries of the underlying meso-granular scale cells. Instances of both emerging magnetic flux loops and magnetic flux disappearing from the surface occur in the simulations. We show an example of a flux tube collapsing to kG field strength and discuss how the nature of flux disappearance can be investigated. Observed stokes profiles of small magnetic structures are severely distorted by telescope diffraction and seeing.
Because of the strong stratification, there is little recycling of plasma and field in the surface layers. Recycling instead occurs by exchange with the deep layers of the convection zone. Plasma and field from the surface descend through the convection zone and rise again toward the surface. Because only a tiny fraction of plasma rising up from deep in the convection zone reaches the surface due to mass conservation, little of the magnetic energy resides in the near surface layers. Thus the dynamo acting on weak incoherent fields is global, rather than a local surface dynamo.
△ Less
Submitted 23 September, 2002;
originally announced September 2002.
-
Characterizing the dynamic properties of the solar turbulence with 3-D simulations: Consequences in terms of p-mode excitation
Authors:
R. Samadi,
A. Nordlund,
R. F. Stein,
M. J. Goupil,
I. Roxburgh
Abstract:
A 3D simulation of the upper part of the solar convective zone is used to derive constraints about the averaged and dynamic properties of solar turbulent convection. Theses constraints are then used to compute the acoustic energy supply rate P(nu) injected into the solar radial oscillations according to the theoretical expression in Samadi & Goupil (2001). The result is compared with solar seism…
▽ More
A 3D simulation of the upper part of the solar convective zone is used to derive constraints about the averaged and dynamic properties of solar turbulent convection. Theses constraints are then used to compute the acoustic energy supply rate P(nu) injected into the solar radial oscillations according to the theoretical expression in Samadi & Goupil (2001). The result is compared with solar seismic data. Assuming, as it is usually done, a gaussian model for the frequency (nu) component chi_k(nu) of the model of turbulence, it is found that the computed P(nu) is underestimated compared with the solar seismic data by a factor ~ 2.5. A frequency analysis of the solar simulation shows that the gaussian model indeed does not correctly model chi_k(nu) in the frequency range where the acoustic energy injected into the solar p-modes is important (nu ~ 2 - 4 mHz). One must consider an additional non-gaussian component for chi_k(nu) to reproduce its behavior. Computed values of P obtained with this non-gaussian component reproduce better the solar seismic observations. This non-gaussian component leads to a Reynolds stress contribution of the same order than the one arising from the advection of the turbulent fluctuations of entropy by the turbulent motions.
△ Less
Submitted 30 August, 2002;
originally announced August 2002.
-
What Causes P-mode Asymmetry Reversal?
Authors:
Dali Georgobiani,
Robert F. Stein,
Aake Nordlund
Abstract:
The solar acoustic p-mode line profiles are asymmetric. Velocity spectra have more power on the low-frequency sides, whereas intensity profiles show the opposite sense of asymmetry. Numerical simulations of the upper convection zone have resonant p-modes with the same asymmetries and asymmetry reversal as the observed modes. The temperature and velocity power spectra at optical depth…
▽ More
The solar acoustic p-mode line profiles are asymmetric. Velocity spectra have more power on the low-frequency sides, whereas intensity profiles show the opposite sense of asymmetry. Numerical simulations of the upper convection zone have resonant p-modes with the same asymmetries and asymmetry reversal as the observed modes. The temperature and velocity power spectra at optical depth $τ_{\rm cont} = 1$ have the opposite asymmetry as is observed for the intensity and velocity spectra. At a fixed geometrical depth, corresponding to $<τ_{\rm cont}>=1$, however, the temperature and velocity spectra have the same asymmetry. This indicates that the asymmetry reversal is produced by radiative transfer effects and not by correlated noise.
△ Less
Submitted 9 May, 2002;
originally announced May 2002.
-
Dynamic Hydrogen Ionization
Authors:
Mats Carlsson,
Robert F. Stein
Abstract:
We investigate the ionization of hydrogen in a dynamic Solar atmosphere. We show that the time scale for ionization/recombination can be estimated from the eigenvalues of a modified rate matrix where the optically thick Lyman transitions that are in detailed balance have been excluded. We find that the time scale for ionization/recombination is dominated by the slow collisional leakage from the…
▽ More
We investigate the ionization of hydrogen in a dynamic Solar atmosphere. We show that the time scale for ionization/recombination can be estimated from the eigenvalues of a modified rate matrix where the optically thick Lyman transitions that are in detailed balance have been excluded. We find that the time scale for ionization/recombination is dominated by the slow collisional leakage from the ground state to the first excited state. Throughout the chromosphere the time scale is long ($10^3$-$10^5$ s), except in shocks where the increased temperature and density shorten the time scale for ionization/recombination, especially in the upper chromosphere. Because the relaxation time scale is much longer than dynamic time scales, hydrogen ionization does not have time to reach its equilibrium value and its fluctuations are much smaller than the variation of its statistical equilibrium value appropriate for the instantaneous conditions. The ionization state tends to represent the higher temperature of the shocks, and the mean electron density is up to a factor of six higher than the electron density calculated in statistical equilibrium from the mean atmosphere. The simulations show that a static picture and a dynamic picture of the chromosphere are fundamentally different and that time variations are crucial for our understanding of the chromosphere itself and the spectral features formed there.
△ Less
Submitted 16 February, 2002;
originally announced February 2002.
-
Solar Convection and Magneto-Convection Simulations
Authors:
Robert F. Stein,
David Bercik,
Aake Nordlund
Abstract:
Magneto-convection simulations with two scenarios have been performed: In one, horizontal magnetic field is advected into the computational domain by fluid entering at the bottom. In the other, an initially uniform vertical magnetic field is imposed on a snapshot of non-magnetic convection and allowed to evolve. In both cases, the field is swept into the intergranular lanes and the boundaries of…
▽ More
Magneto-convection simulations with two scenarios have been performed: In one, horizontal magnetic field is advected into the computational domain by fluid entering at the bottom. In the other, an initially uniform vertical magnetic field is imposed on a snapshot of non-magnetic convection and allowed to evolve. In both cases, the field is swept into the intergranular lanes and the boundaries of the underlying mesogranules. The largest field concentrations at the surface reach pressure balance with the surrounding gas. They suppress both horizontal and vertical flows, which reduces the heat transport. They cool, become evacuated and their optical depth unity surface is depressed by several hundred kilometers. Micropores form, typically where a small granule disappears and surrounding flux tubes squeeze into its previous location.
△ Less
Submitted 5 December, 2001;
originally announced December 2001.
-
Magnetohydrodynamic turbulence in warped accretion discs
Authors:
U. Torkelsson,
G. I. Ogilvie,
A. Brandenburg,
J. E. Pringle,
AA. Nordlund,
R. F. Stein
Abstract:
Warped, precessing accretion discs appear in a range of astrophysical systems, for instance the X-ray binary Her X-1 and in the active nucleus of NGC4258. In a warped accretion disc there are horizontal pressure gradients that drive an epicyclic motion. We have studied the interaction of this epicyclic motion with the magnetohydrodynamic turbulence in numerical simulations. We find that the turb…
▽ More
Warped, precessing accretion discs appear in a range of astrophysical systems, for instance the X-ray binary Her X-1 and in the active nucleus of NGC4258. In a warped accretion disc there are horizontal pressure gradients that drive an epicyclic motion. We have studied the interaction of this epicyclic motion with the magnetohydrodynamic turbulence in numerical simulations. We find that the turbulent stress acting on the epicyclic motion is comparable in size to the stress that drives the accretion, however an important ingredient in the damping of the epicyclic motion is its parametric decay into inertial waves.
△ Less
Submitted 3 March, 2001;
originally announced March 2001.
-
Solar Oscillations and Convection: II. Excitation of Radial Oscillations
Authors:
R. F. Stein,
Å. Nordlund
Abstract:
Solar p-mode oscillations are excited by the work of stochastic, non-adiabatic, pressure fluctuations on the compressive modes. We evaluate the expression for the radial mode excitation rate derived by Nordlund and Stein (Paper I) using numerical simulations of near surface solar convection. We first apply this expression to the three radial modes of the simulation and obtain good agreement betw…
▽ More
Solar p-mode oscillations are excited by the work of stochastic, non-adiabatic, pressure fluctuations on the compressive modes. We evaluate the expression for the radial mode excitation rate derived by Nordlund and Stein (Paper I) using numerical simulations of near surface solar convection. We first apply this expression to the three radial modes of the simulation and obtain good agreement between the predicted excitation rate and the actual mode damping rates as determined from their energies and the widths of their resolved spectral profiles. We then apply this expression for the mode excitation rate to the solar modes and obtain excellent agreement with the low l damping rates determined from GOLF data. Excitation occurs close to the surface, mainly in the intergranular lanes and near the boundaries of granules (where turbulence and radiative cooling are large). The non-adiabatic pressure fluctuations near the surface are produced by small instantaneous local imbalances between the divergence of the radiative and convective fluxes near the solar surface. Below the surface, the non-adiabatic pressure fluctuations are produced primarily by turbulent pressure fluctuations (Reynolds stresses). The frequency dependence of the mode excitation is due to effects of the mode structure and the pressure fluctuation spectrum. Excitation is small at low frequencies due to mode properties -- the mode compression decreases and the mode mass increases at low frequency. Excitation is small at high frequencies due to the pressure fluctuation spectrum -- pressure fluctuations become small at high frequencies because they are due to convection which is a long time scale phenomena compared to the dominant p-mode periods.
△ Less
Submitted 2 August, 2000;
originally announced August 2000.
-
Solar Oscillations and Convection: I. Formalism for Radial Oscillations
Authors:
Å. Nordlund,
R. F. Stein
Abstract:
We present a formalism for investigating the interaction between pmode oscillations and convection by analyzing realistic, three-dimensional simulations of the near-surface layers of the solar convection zone. By choosing suitable definitions for fluctuations and averages, we obtain a separation that retains exact equations. The equations for the horizontal averages contain one part that corresp…
▽ More
We present a formalism for investigating the interaction between pmode oscillations and convection by analyzing realistic, three-dimensional simulations of the near-surface layers of the solar convection zone. By choosing suitable definitions for fluctuations and averages, we obtain a separation that retains exact equations. The equations for the horizontal averages contain one part that corresponds directly to the wave equations for a 1-D medium, plus additional terms that arise from the averaging and correspond to the turbulent pressure gradient in the momentum equation and the divergence of the convective and kinetic energy fluxes in the internal energy equation. These terms cannot be evaluated in closed form, but they may be measured in numerical simulations. The additional terms may cause the mode frequencies to shift, relative to what would be obtained if only the terms corresponding to a 1-D medium were retained---most straightforwardly by changing the mean stratification, and more subtly by changing the effective compressibility of the medium. In the presence of time dependent convection, the additional terms also have a stochastic time dependence, that acts as a source of random excitation of the coherent modes. In the present paper, we derive an expression for the excitation power and test it by applying it to a numerical experiment of sufficient duration for the excited modes to be spectrally resolved.
△ Less
Submitted 23 June, 2000;
originally announced June 2000.
-
Line formation in solar granulation: II. The photospheric Fe abundance
Authors:
M. Asplund,
AA. Nordlund,
R. Trampedach,
R. F. Stein
Abstract:
The solar photospheric Fe abundance has been determined using realistic ab initio 3D, time-dependent, hydrodynamical model atmospheres. The study is based on the excellent agreement between the predicted and observed line profiles directly rather than equivalent width, since the intrinsic Doppler broadening from the convective motions and oscillations provide the necessary non-thermal broadening…
▽ More
The solar photospheric Fe abundance has been determined using realistic ab initio 3D, time-dependent, hydrodynamical model atmospheres. The study is based on the excellent agreement between the predicted and observed line profiles directly rather than equivalent width, since the intrinsic Doppler broadening from the convective motions and oscillations provide the necessary non-thermal broadening. Thus, three of the four hotly debated parameters (equivalent widths, microturbulence and damping enhancement factors) in the center of the recent solar Fe abundance dispute regarding FeI lines no longer enter the analysis, leaving the transition probabilities as the main uncertainty. Both FeI (using the samples of lines of both the Oxford and Kiel studies) and FeII lines have been investigated, which give consistent results: log FeI = 7.44 +- 0.05 and log FeII = 7.45 +- 0.10. Also the wings of strong FeI lines return consistent abundances, log FeII = 7.42 +- 0.03, but due to the uncertainties inherent in analyses of strong lines we give this determination lower weight than the results from weak and intermediate strong lines. In view of the recent slight downward revision of the meteoritic Fe abundance log Fe = 7.46 +- 0.01, the agreement between the meteoritic and photospheric values is very good, thus appearingly settling the debate over the photospheric Fe abundance from FeI lines.
△ Less
Submitted 15 May, 2000;
originally announced May 2000.
-
Line formation in solar granulation: I. Fe line shapes, shifts and asymmetries
Authors:
M. Asplund,
AA. Nordlund,
R. Trampedach,
C. Allende Prieto,
R. F. Stein
Abstract:
Realistic ab-initio 3D, radiative-hydrodynamical convection simulations of the solar granulation have been applied to FeI and FeII line formation. In contrast to classical analyses based on hydrostatic 1D model atmospheres the procedure contains no adjustable free parameters but the treatment of the numerical viscosity in the construction of the 3D, time-dependent, inhomogeneous model atmosphere…
▽ More
Realistic ab-initio 3D, radiative-hydrodynamical convection simulations of the solar granulation have been applied to FeI and FeII line formation. In contrast to classical analyses based on hydrostatic 1D model atmospheres the procedure contains no adjustable free parameters but the treatment of the numerical viscosity in the construction of the 3D, time-dependent, inhomogeneous model atmosphere and the elemental abundance in the 3D spectral synthesis. However, the numerical viscosity is introduced purely for numerical stability purposes and is determined from standard hydrodynamical test cases with no adjustments allowed to improve the agreement with the observational constraints from the solar granulation. The non-thermal line broadening is mainly provided by the Doppler shifts arising from the convective flows in the solar photosphere and the solar oscillations. The almost perfect agreement between the predicted temporally and spatially averaged line profiles for weak Fe lines with the observed profiles and the absence of trends in derived abundances with line strengths, seem to imply that the micro- and macroturbulence concepts are obsolete in these 3D analyses. Furthermore, the theoretical line asymmetries and shifts show a very satisfactory agreement with observations with an accuracy of typically 50-100 m/s on an absolute velocity scale. The remaining minor discrepancies point to how the convection simulations can be refined further.
△ Less
Submitted 15 May, 2000;
originally announced May 2000.
-
The effects of numerical resolution on hydrodynamical surface convection simulations and spectral line formation
Authors:
M. Asplund,
H. -G. Ludwig,
AA. Nordlund,
R. F. Stein
Abstract:
The computationally demanding nature of radiative-hydrodynamical simulations of stellar surface convection warrants an investigation of the sensitivity of the convective structure and spectral synthesis to the numerical resolution and dimension of the simulations, which is presented here. With too coarse a resolution the predicted spectral lines tend to be too narrow, reflecting insufficient Dop…
▽ More
The computationally demanding nature of radiative-hydrodynamical simulations of stellar surface convection warrants an investigation of the sensitivity of the convective structure and spectral synthesis to the numerical resolution and dimension of the simulations, which is presented here. With too coarse a resolution the predicted spectral lines tend to be too narrow, reflecting insufficient Doppler broadening from the convective motions, while at the currently highest affordable resolution the line shapes have converged essentially perfectly to the observed profiles. Similar conclusions are drawn from the line asymmetries and shifts. In terms of abundances, weak FeI and FeII lines show a very small dependence (~0.02 dex) while for intermediate strong lines with significant non-thermal broadening the sensitivity increases (~0.10 dex). Problems arise when using 2D convection simulations to describe an inherent 3D phenomenon, which translates to inaccurate atmospheric velocity fields and temperature and pressure structures. In 2D the theoretical line profiles tend to be too shallow and broad compared with the 3D calculations and observations, in particular for intermediate strong lines. In terms of abundances, the 2D results are systematically about 0.1 dex lower than for the 3D case for FeI lines. Furthermore, the predicted line asymmetries and shifts are much inferior in 2D. Given these shortcomings and computing time considerations it is better to use 3D simulations of even modest resolution than high-resolution 2D simulations.
△ Less
Submitted 15 May, 2000;
originally announced May 2000.
-
The response of a turbulent accretion disc to an imposed epicyclic shearing motion
Authors:
Ulf Torkelsson,
Gordon I. Ogilvie,
Axel Brandenburg,
James E. Pringle,
Ake Nordlund,
Robert F. Stein
Abstract:
We excite an epicyclic motion, whose amplitude depends on the vertical position, $z$, in a simulation of a turbulent accretion disc. An epicyclic motion of this kind may be caused by a warping of the disc. By studying how the epicyclic motion decays we can obtain information about the interaction between the warp and the disc turbulence. A high amplitude epicyclic motion decays first by exciting…
▽ More
We excite an epicyclic motion, whose amplitude depends on the vertical position, $z$, in a simulation of a turbulent accretion disc. An epicyclic motion of this kind may be caused by a warping of the disc. By studying how the epicyclic motion decays we can obtain information about the interaction between the warp and the disc turbulence. A high amplitude epicyclic motion decays first by exciting inertial waves through a parametric instability, but its subsequent exponential damping may be reproduced by a turbulent viscosity. We estimate the effective viscosity parameter, $α_{\rm v}$, pertaining to such a vertical shear. We also gain new information on the properties of the disc turbulence in general, and measure the usual viscosity parameter, $α_{\rm h}$, pertaining to a horizontal (Keplerian) shear. We find that, as is often assumed in theoretical studies, $α_{\rm v}$ is approximately equal to $α_{\rm h}$ and both are much less than unity, for the field strengths achieved in our local box calculations of turbulence. In view of the smallness ($\sim 0.01$) of $α_{\rm v}$ and $α_{\rm h}$ we conclude that for $β= p_{\rm gas}/p_{\rm mag} \sim 10$ the timescale for diffusion or damping of a warp is much shorter than the usual viscous timescale. Finally, we review the astrophysical implications.
△ Less
Submitted 9 May, 2000;
originally announced May 2000.
-
Simulations of Oscillation Modes of the Solar Convection Zone
Authors:
D. Georgobiani,
A. G. Kosovichev,
R. Nigam,
A. Nordlund,
R. F. Stein
Abstract:
We use the three-dimensional hydrodynamic code of Stein and Nordlund to realistically simulate the upper layers of the solar convection zone in order to study physical characteristics of solar oscillations. Our first result is that the properties of oscillation modes in the simulation closely match the observed properties. Recent observations from SOHO/MDI and GONG have confirmed the asymmetry o…
▽ More
We use the three-dimensional hydrodynamic code of Stein and Nordlund to realistically simulate the upper layers of the solar convection zone in order to study physical characteristics of solar oscillations. Our first result is that the properties of oscillation modes in the simulation closely match the observed properties. Recent observations from SOHO/MDI and GONG have confirmed the asymmetry of solar oscillation line profiles, initially discovered by Duvall et al. In this paper we compare the line profiles in the power spectra of the Doppler velocity and continuum intensity oscillations from the SOHO/MDI observations with the simulation. We also compare the phase differences between the velocity and intensity data. We have found that the simulated line profiles are asymmetric and have the same asymmetry reversal between velocity and intensity as observed. The phase difference between the velocity and intensity signals is negative at low frequencies and jumps in the vicinity of modes as is also observed. Thus, our numerical model reproduces the basic observed properties of solar oscillations, and allows us to study the physical properties which are not observed.
△ Less
Submitted 23 December, 1999;
originally announced December 1999.
-
3D hydrodynamical model atmospheres of metal-poor stars: Evidence for a low primordial Li abundance
Authors:
M. Asplund,
AA. Nordlund,
R. Trampedach,
R. F. Stein
Abstract:
Realistic 3-dimensional (3D), radiative hydrodynamical surface convection simulations of the metal-poor halo stars HD 140283 and HD 84937 have been performed. Due to the dominance of adiabatic cooling over radiative heating very low atmospheric temperatures are encountered. The lack of spectral lines in these metal-poor stars thus causes much steeper temperature gradients than in classical 1D hy…
▽ More
Realistic 3-dimensional (3D), radiative hydrodynamical surface convection simulations of the metal-poor halo stars HD 140283 and HD 84937 have been performed. Due to the dominance of adiabatic cooling over radiative heating very low atmospheric temperatures are encountered. The lack of spectral lines in these metal-poor stars thus causes much steeper temperature gradients than in classical 1D hydrostatic model atmospheres where the temperature of the optically thin layers is determined by radiative equilibrium. The modified atmospheric structures cause changes in the emergent stellar spectra. In particular, the primordial Li abundances may have been overestimated by 0.2-0.35 dex with 1D model atmospheres. However, we caution that our result assumes local thermodynamic equilibrium (LTE), while the steep temperature gradients may be prone to e.g. over-ionization.
△ Less
Submitted 6 May, 1999;
originally announced May 1999.
-
Test of Convective Frequency Effects with SOI/MDI High-Degree Data
Authors:
C. S. Rosenthal,
J. Christensen-Dalsgaard,
A. G. Kosovichev,
Å. Nordlund,
J. Reiter,
E. J. Rhodes Jr,
J. Schou,
R. F. Stein,
R. Trampedach
Abstract:
Advances in hydrodynamical simulations have provided new insights into the effects of convection on the frequencies of solar oscillations. As more accurate observations become available, this may lead to an improved understanding of the dynamics of convection and the interaction between convection and pulsation (\cite{Rosenthal+99}). Recent high-resolution observations from the SOI/MDI instrumen…
▽ More
Advances in hydrodynamical simulations have provided new insights into the effects of convection on the frequencies of solar oscillations. As more accurate observations become available, this may lead to an improved understanding of the dynamics of convection and the interaction between convection and pulsation (\cite{Rosenthal+99}). Recent high-resolution observations from the SOI/MDI instrument on the SOHO spacecraft have provided the so-far most-detailed observations of high-degree modes of solar oscillations, which are particularly sensitive to the near-surface properties of the Sun. Here we present preliminary results of a comparison between these observations and frequencies computed for models based on realistic simulations of near-surface convection. Such comparisons may be expected to help in identifying the causes of the remaining differences between the observed frequencies and those of solar models.
△ Less
Submitted 7 July, 1998;
originally announced July 1998.
-
Convective contributions to the frequencies of solar oscillations
Authors:
C. S. Rosenthal,
J. Christensen-Dalsgaard,
A. Nordlund,
R. F. Stein,
R. Trampedach
Abstract:
Differences between observed and theoretical eigenfrequencies of the Sun have characteristics which identify them as arising predominantly from properties of the oscillations in the vicinity of the solar surface: in the super-adiabatic, convective boundary layer and above. These frequency differences may therefore provide useful information about the structure of these regions, precisely where t…
▽ More
Differences between observed and theoretical eigenfrequencies of the Sun have characteristics which identify them as arising predominantly from properties of the oscillations in the vicinity of the solar surface: in the super-adiabatic, convective boundary layer and above. These frequency differences may therefore provide useful information about the structure of these regions, precisely where the theory of solar structure is most uncertain.
In the present work we use numerical simulations of the outer part of the Sun to quantify the influence of turbulent convection on solar oscillation frequencies. Separating the influence into effects on the mean model and effects on the physics of the modes, we find that the main model effects are due to the turbulent pressure that provides additional support against gravity, and thermal differences between average 3-D models and 1-D models. Surfaces of constant pressure in the visible photosphere are elevated by about 150 km, relative to a standard envelope model.
As a result, the turning points of high-frequency modes are raised, while those of the low-frequency modes remain essentially unaffected. The corresponding gradual lowering of the mode frequencies accounts for most of the frequency difference between observations and standard solar models. Additional effects are expected to come primarily from changes in the physics of the modes, in particular from the modulation of the turbulent pressure by the oscillations.
△ Less
Submitted 16 March, 1999; v1 submitted 17 March, 1998;
originally announced March 1998.
-
Dynamics of Magnetic Flux Elements in the Solar Photosphere
Authors:
A. A. van Ballegooijen,
P. Nisenson,
R. W. Noyes,
M. G. Lofdahl,
R. F. Stein,
A. Nordlund,
V. Krishnakumar
Abstract:
The interaction of magnetic fields and convection is investigated in the context of the coronal heating problem. We study the motions of photospheric magnetic elements using filtergrams obtained at the Swedish Vacuum Solar Telescope at La Palma. We use potential-field modeling to extrapolate the magnetic and velocity fields to larger height. We find that the velocity in the chromosphere can be l…
▽ More
The interaction of magnetic fields and convection is investigated in the context of the coronal heating problem. We study the motions of photospheric magnetic elements using filtergrams obtained at the Swedish Vacuum Solar Telescope at La Palma. We use potential-field modeling to extrapolate the magnetic and velocity fields to larger height. We find that the velocity in the chromosphere can be locally enhanced at the separatrix surfaces between neighboring flux tubes. The predicted velocities are several km/s, significantly larger than those of the photospheric flux tubes, which may have important implications for coronal heating. s
△ Less
Submitted 27 February, 1998;
originally announced February 1998.
-
Does a Non-Magnetic Solar Chromosphere Exist?
Authors:
Mats Carlsson,
Robert F. Stein
Abstract:
Enhanced chromospheric emission which corresponds to an outwardly increasing semiempirical temperature structure can be produced by wave motion without any increase in the mean gas temperture. Hence, the sun may not have a classical chromosphere in magnetic field free internetwork regions. Other significant differences between the properties of dynamic and static atmospheres should be considered…
▽ More
Enhanced chromospheric emission which corresponds to an outwardly increasing semiempirical temperature structure can be produced by wave motion without any increase in the mean gas temperture. Hence, the sun may not have a classical chromosphere in magnetic field free internetwork regions. Other significant differences between the properties of dynamic and static atmospheres should be considered when analyzing chromospheric observations.
△ Less
Submitted 9 November, 1994;
originally announced November 1994.