-
Leptocline as a Shallow Substructure of Near-Surface Shear Layer in 3D Radiative Hydrodynamic Simulations
Authors:
Irina N. Kitiashvili,
Alexander G. Kosovichev,
Alan A. Wray,
Viacheslav M. Sadykov,
Gustavo Guerrero
Abstract:
Understanding effects driven by rotation in the solar convection zone is essential for many problems related to solar activity, such as the formation of differential rotation, meridional circulation, and others. We analyze realistic 3D radiative hydrodynamics simulations of solar subsurface dynamics in the presence of rotation in a local domain 80 Mm wide and 25 Mm deep, located at 30 degrees lati…
▽ More
Understanding effects driven by rotation in the solar convection zone is essential for many problems related to solar activity, such as the formation of differential rotation, meridional circulation, and others. We analyze realistic 3D radiative hydrodynamics simulations of solar subsurface dynamics in the presence of rotation in a local domain 80 Mm wide and 25 Mm deep, located at 30 degrees latitude. The simulation results reveal the development of a shallow 10-Mm deep substructure of the Near-Surface Shear Layer (NSSL), characterized by a strong radial rotational gradient and self-organized meridional flows. This shallow layer ("leptocline") is located in the hydrogen ionization zone associated with enhanced anisotropic overshooting-type flows into a less unstable layer between the H and HeII ionization zones. We discuss current observational evidence of the presence of the leptocline and show that the radial variations of the differential rotation and meridional flow profiles obtained from the simulations in this layer qualitatively agree with helioseismic observations.
△ Less
Submitted 27 August, 2022; v1 submitted 2 March, 2022;
originally announced March 2022.
-
Modeling Stellar Jitter for the Detection of Earth-Mass Exoplanets via Precision Radial Velocity Measurements
Authors:
Samuel Granovsky,
Irina N. Kitiashvili,
Alan Wray
Abstract:
The detection of Earth-size exoplanets is a technological and data analysis challenge. Future progress in Earth-mass exoplanet detection is expected from the development of extreme precision radial velocity measurements. Increasing radial velocity precision requires developing a new physics-based data analysis methodology to discriminate planetary signals from host-star-related effects, taking ste…
▽ More
The detection of Earth-size exoplanets is a technological and data analysis challenge. Future progress in Earth-mass exoplanet detection is expected from the development of extreme precision radial velocity measurements. Increasing radial velocity precision requires developing a new physics-based data analysis methodology to discriminate planetary signals from host-star-related effects, taking stellar variability and instrumental uncertainties into account. In this work, we investigate and quantify stellar disturbances of the planet-hosting solar-type star HD121504 from 3D radiative modeling obtained with the StellarBox code. The model has been used for determining statistical properties of the turbulent plasma and obtaining synthetic spectroscopic observations for several Fe I lines at different locations on the stellar disk to mimic high-resolution spectroscopic observations.
△ Less
Submitted 14 February, 2022;
originally announced February 2022.
-
3D Modeling of the Structure and Dynamics of a Main-Sequence F-type Star
Authors:
Irina N. Kitiashvili,
Alan A. Wray
Abstract:
Current state-of-the-art computational modeling makes it possible to build realistic models of stellar convection zones and atmospheres that take into account chemical composition, radiative effects, ionization, and turbulence. The standard 1D mixing-length-based evolutionary models are not able to capture many physical processes of the stellar interior dynamics. Mixing-length models provide an in…
▽ More
Current state-of-the-art computational modeling makes it possible to build realistic models of stellar convection zones and atmospheres that take into account chemical composition, radiative effects, ionization, and turbulence. The standard 1D mixing-length-based evolutionary models are not able to capture many physical processes of the stellar interior dynamics. Mixing-length models provide an initial approximation of stellar structure that can be used to initialize 3D radiative hydrodynamics simulations which include realistic modeling of turbulence, radiation, and other phenomena.
In this paper, we present 3D radiative hydrodynamic simulations of an F-type main-sequence star with 1.47 solar mass. The computational domain includes the upper layers of the radiation zone, the entire convection zone, and the photosphere. The effects of stellar rotation is modeled in the f-plane approximation. These simulations provide new insight into the properties of the convective overshoot region, the dynamics of the near-surface, highly turbulent layer, and the structure and dynamics of granulation. They reveal anti-solar-type differential rotation and latitudinal dependence of the tachocline location.
△ Less
Submitted 26 July, 2021;
originally announced July 2021.
-
Challenges and Advances in Modeling of the Solar Atmosphere: A White Paper of Findings and Recommendations
Authors:
Serena Criscuoli,
Maria Kazachenko,
Irina Kitashvili,
Alexander Kosovichev,
Valentín Martínez Pillet,
Gelu Nita,
Viacheslav Sadykov,
Alan Wray
Abstract:
The next decade will be an exciting period for solar astrophysics, as new ground- and space-based instrumentation will provide unprecedented observations of the solar atmosphere and heliosphere. The synergy between modeling effort and comprehensive analysis of observations is crucial for the understanding of the physical processes behind the observed phenomena. However, the unprecedented wealth of…
▽ More
The next decade will be an exciting period for solar astrophysics, as new ground- and space-based instrumentation will provide unprecedented observations of the solar atmosphere and heliosphere. The synergy between modeling effort and comprehensive analysis of observations is crucial for the understanding of the physical processes behind the observed phenomena. However, the unprecedented wealth of data on one hand, and the complexity of the physical phenomena on the other, require the development of new approaches in both data analysis and numerical modeling. In this white paper, we summarize recent numerical achievements to reproduce structure, dynamics, and observed phenomena from the photosphere to the low corona and outline challenges we expect to face for the interpretation of future observations.
△ Less
Submitted 7 January, 2021; v1 submitted 30 December, 2020;
originally announced January 2021.
-
Connecting Atmospheric Properties and Synthetic Emission of Shock Waves Using 3D RMHD Simulations of Quiet Sun
Authors:
Viacheslav M. Sadykov,
Irina N. Kitiashvili,
Alexander G. Kosovichev,
Alan A. Wray
Abstract:
We analyze the evolution of shock waves in high-resolution 3D radiative MHD simulations of the quiet Sun and their synthetic emission characteristics. The simulations model the dynamics of a 12.8x12.8x15.2 Mm quiet-Sun region (including a 5.2 Mm layer of the upper convection zone and a 10 Mm atmosphere from the photosphere to corona) with an initially uniform vertical magnetic field of 10 G, natur…
▽ More
We analyze the evolution of shock waves in high-resolution 3D radiative MHD simulations of the quiet Sun and their synthetic emission characteristics. The simulations model the dynamics of a 12.8x12.8x15.2 Mm quiet-Sun region (including a 5.2 Mm layer of the upper convection zone and a 10 Mm atmosphere from the photosphere to corona) with an initially uniform vertical magnetic field of 10 G, naturally driven by convective flows. We synthesize the Mg II and C II spectral lines observed by the IRIS satellite and EUV emission observed by the SDO/AIA telescope. Synthetic observations are obtained using the RH1.5D radiative transfer code and temperature response functions at both the numerical and instrumental resolutions. We found that the Doppler velocity jumps of the C II 1334.5 A IRIS line and a relative enhancement of the emission in the 335 A SDO/AIA channel are the best proxies for the enthalpy deposited by shock waves into the corona (with Kendall's $τ$ correlation coefficients of 0.59 and 0.38, respectively). The synthetic emission of the lines and extreme ultraviolet passbands are correlated with each other during the shock wave propagation. All studied shocks are mostly hydrodynamic (i.e., the magnetic energy carried by horizontal fields is $\leq{}$2.6% of the enthalpy for all events) and have Mach numbers > 1.0-1.2 in the low corona. The study reveals the possibility of diagnosing energy transport by shock waves into the solar corona, as well as their other properties, by using IRIS and SDO/AIA sensing observations.
△ Less
Submitted 12 March, 2021; v1 submitted 13 August, 2020;
originally announced August 2020.
-
Origin of Deep Acoustic Sources Associated with Solar Magnetic Structures
Authors:
Irina N. Kitiashvili,
Alexander G. Kosovichev,
Nagi N. Mansour,
Alan A. Wray,
Timothy A. Sandstrom
Abstract:
It is generally accepted that solar acoustic (p) modes are excited by near-surface turbulent motions, in particular, by downdrafts and interacting vortices in intergranular lanes. Recent analysis of Solar Dynamics Observatory data by (Zhao et al., 2015) revealed fast-moving waves around sunspots, which are consistent with magnetoacoustic waves excited approximately 5 Mm beneath the sunspot. We ana…
▽ More
It is generally accepted that solar acoustic (p) modes are excited by near-surface turbulent motions, in particular, by downdrafts and interacting vortices in intergranular lanes. Recent analysis of Solar Dynamics Observatory data by (Zhao et al., 2015) revealed fast-moving waves around sunspots, which are consistent with magnetoacoustic waves excited approximately 5 Mm beneath the sunspot. We analyzed 3D radiative MHD simulations of solar magnetoconvection with a self-organized pore-like magnetic structure, and identified more than 600 individual acoustic events both inside and outside of this structure. By performing a case-by-case study, we found that surrounding the magnetic structure, acoustic sources are associated with downdrafts. Their depth correlates with downdraft speed and magnetic fields. The sources often can be transported into deeper layers by downdrafts. The wave front shape, in the case of a strong or inclined downdraft, can be stretched along the downdraft. Inside the magnetic structure, excitation of acoustic waves is driven by converging flows. Frequently, strong converging plasma streams hit the structure boundaries, causing compressions in its interior that excite acoustic waves. Analysis of the depth distribution of acoustic events shows the strongest concentration at 0.2 - 1 Mm beneath the surface for the outside sources and 2.5 - 3 Mm for the excitation event inside the structure.
△ Less
Submitted 14 October, 2018;
originally announced October 2018.
-
Dynamics of Turbulent Convection and Convective Overshoot in a Moderate Mass Star
Authors:
Irina N. Kitiashvili,
Alexander G. Kosovichev,
Nagi N. Mansour,
Alan A. Wray
Abstract:
Continued progress in observational stellar astrophysics requires a deep understanding of the underlying convection dynamics. We present results of realistic 3D radiative hydrodynamic simulations of the outer layers of a moderate mass star (1.47 Msun), including the full convection zone, the overshoot region, and the top layers of the radiative zone. The simulation results show that the surface gr…
▽ More
Continued progress in observational stellar astrophysics requires a deep understanding of the underlying convection dynamics. We present results of realistic 3D radiative hydrodynamic simulations of the outer layers of a moderate mass star (1.47 Msun), including the full convection zone, the overshoot region, and the top layers of the radiative zone. The simulation results show that the surface granulation has a broad range of scales, from 2 to 12 Mm, and that large granules are organized in well-defined clusters, consisting of several granules. Comparison of the mean structure profiles from 3D simulations with the corresponding 1D standard stellar model shows an increase of the stellar radius by ~800 km, as well as significant changes in the thermodynamic structure and turbulent properties of the ionization zones. Convective downdrafts in the intergranular lanes between granulation clusters reach speeds of more than 20 km/s, penetrate through the whole convection zone, hit the radiative zone, and form a 8 Mm thick overshoot layer. Contrary to semi-empirical overshooting models, our results show that the 3D dynamic overshoot region consists of two layers: a nearly adiabatic extension of the convection zone and a deeper layer of enhanced subadiabatic stratification. This layer is formed because of heating caused by the braking of the overshooting convective plumes. This effect has to be taken into account in stellar modeling and the interpretation of asteroseismology data. In particular, we demonstrate that the 3D model can qualitatively explain deviations from the standard solar model found by helioseismology.
△ Less
Submitted 22 December, 2015;
originally announced December 2015.
-
Simulations of Stellar Magnetoconvection using the Radiative MHD Code `StellarBox'
Authors:
Alan A. Wray,
Khalil Bensassi,
Irina N. Kitiashvili,
Nagi N. Mansour,
Alexander G. Kosovichev
Abstract:
Realistic numerical simulations, i.e., those that make minimal use of ad hoc modeling, are essential for understanding the complex turbulent dynamics of the interiors and atmospheres of the Sun and other stars and the basic mechanisms of their magnetic activity and variability. The goal of this paper is to present a detailed description and test results of a compressible radiative MHD code, `Stell…
▽ More
Realistic numerical simulations, i.e., those that make minimal use of ad hoc modeling, are essential for understanding the complex turbulent dynamics of the interiors and atmospheres of the Sun and other stars and the basic mechanisms of their magnetic activity and variability. The goal of this paper is to present a detailed description and test results of a compressible radiative MHD code, `StellarBox', specifically developed for simulating the convection zones, surface, and atmospheres of the Sun and moderate-mass stars. The code solves the three-dimensional, fully coupled compressible MHD equations using a fourth-order Padé spatial differentiation scheme and a fourth-order Runge-Kutta scheme for time integration. The radiative transfer equation is solved using the Feautrier method for bi-directional ray tracing and an opacity-binning technique. A specific feature of the code is the implementation of subgrid-scale MHD turbulence models. The data structures are automatically configured, depending on the computational grid and the number of available processors, to achieve good load balancing. We present test results and illustrate the code's capabilities for simulating the granular convection on the Sun and a set of main-sequence stars. The results reveal substantial changes in the near-surface turbulent convection in these stars, which in turn affect properties of the surface magnetic fields. For example, in the solar case initially uniform vertical magnetic fields tend to self-organize into compact (pore-like) magnetic structures, while in more massive stars such structures are not formed and the magnetic field is distributed more-or-less uniformly in the intergranular lanes.
△ Less
Submitted 28 July, 2015;
originally announced July 2015.
-
Realistic modeling of local dynamo processes on the Sun
Authors:
I. N. Kitiashvili,
A. G. Kosovichev,
N. N. Mansour,
A. A. Wray
Abstract:
Magnetic fields are usually observed in the quiet Sun as small-scale elements that cover the entire solar surface (the `salt and pepper' patterns in line-of-sight magnetograms). By using 3D radiative MHD numerical simulations we find that these fields result from a local dynamo action in the top layers of the convection zone, where extremely weak 'seed' magnetic fields (e.g., from a $10^{-6}$ G) c…
▽ More
Magnetic fields are usually observed in the quiet Sun as small-scale elements that cover the entire solar surface (the `salt and pepper' patterns in line-of-sight magnetograms). By using 3D radiative MHD numerical simulations we find that these fields result from a local dynamo action in the top layers of the convection zone, where extremely weak 'seed' magnetic fields (e.g., from a $10^{-6}$ G) can locally grow above the mean equipartition field, to a stronger than 2000~G field localized in magnetic structures. Our results reveal that the magnetic flux is predominantly generated in regions of small-scale helical downflows. We find that the local dynamo action takes place mostly in a shallow, about 500~km deep, subsurface layer, from which the generated field is transported into the deeper layers by convective downdrafts. We demonstrate that the observed dominance of vertical magnetic fields at the photosphere and horizontal fields above the photosphere can be explained by small-scale magnetic loops produced by the dynamo. Such small-scale loops play an important role in the structure and dynamics of the solar atmosphere and that their detection in observations is critical for understanding the local dynamo action on the Sun.
△ Less
Submitted 29 June, 2015;
originally announced June 2015.
-
Mechanism of local dynamo action on the Sun
Authors:
I. N. Kitiashvili,
A. G. Kosovichev,
N. N. Mansour,
A. A. Wray
Abstract:
In the quiet Sun, magnetic fields are usually observed as small-scale magnetic elements, `salt and pepper', covering the entire solar surface. By using 3D radiative MHD numerical simulations we demonstrate that these fields are a result of local dynamo action in the top layers of the convection zone, where extremely weak `seed' magnetic fields can locally grow above the mean equipartition field (e…
▽ More
In the quiet Sun, magnetic fields are usually observed as small-scale magnetic elements, `salt and pepper', covering the entire solar surface. By using 3D radiative MHD numerical simulations we demonstrate that these fields are a result of local dynamo action in the top layers of the convection zone, where extremely weak `seed' magnetic fields can locally grow above the mean equipartition field (e.g., from a $10^{-6}$ G `seed' field to more than 1000 G magnetic structures). We find that the local dynamo action takes place only in a shallow, about 500 km deep, subsurface layer, from which the generated field is transported into deeper layers by convection downdrafts. We demonstrate that the observed dominance of vertical magnetic fields at the photosphere and the horizontal fields above the photosphere can be explained by multi-scale magnetic loops produced by the dynamo.
△ Less
Submitted 3 December, 2013;
originally announced December 2013.
-
Ubiquitous Solar Eruptions Driven by Magnetized Vortex Tubes
Authors:
I. N. Kitiashvili,
A. G. Kosovichev,
S. K. Lele,
N. N. Mansour,
A. A. Wray
Abstract:
The solar surface is covered by high-speed jets transporting mass and energy into the solar corona and feeding the solar wind. The most prominent of these jets have been known as spicules. However, the mechanism initiating these eruptions events is still unknown. Using realistic numerical simulations we find that small-scale eruptions are produced by ubiquitous magnetized vortex tubes generated by…
▽ More
The solar surface is covered by high-speed jets transporting mass and energy into the solar corona and feeding the solar wind. The most prominent of these jets have been known as spicules. However, the mechanism initiating these eruptions events is still unknown. Using realistic numerical simulations we find that small-scale eruptions are produced by ubiquitous magnetized vortex tubes generated by the Sun's turbulent convection in subsurface layers. The swirling vortex tubes (resembling tornadoes) penetrate into the solar atmosphere, capture and stretch background magnetic field, and push surrounding material up, generating quasiperiodic shocks. Our simulations reveal a complicated high-speed flow patterns, and thermodynamic and magnetic structure in the erupting vortex tubes. We found that the eruptions are initiated in the subsurface layers and are driven by the high-pressure gradients in the subphotosphere and photosphere, and by the Lorentz force in the higher atmosphere layers.
△ Less
Submitted 31 December, 2012;
originally announced January 2013.
-
Turbulent Kinetic Energy Spectra of Solar Convection from NST Observations and Realistic MHD Simulations
Authors:
I. N. Kitiashvili,
V. I. Abramenko,
P. R. Goode,
A. G. Kosovichev,
S. K. Lele,
N. N. Mansour,
A. A. Wray,
V. B. Yurchyshyn
Abstract:
Turbulent properties of the quiet Sun represent the basic state of surface conditions, and a background for various processes of solar activity. Therefore understanding of properties and dynamics of this `basic' state is important for investigation of more complex phenomena, formation and development of observed phenomena in the photosphere and atmosphere. For characterization of the turbulent pro…
▽ More
Turbulent properties of the quiet Sun represent the basic state of surface conditions, and a background for various processes of solar activity. Therefore understanding of properties and dynamics of this `basic' state is important for investigation of more complex phenomena, formation and development of observed phenomena in the photosphere and atmosphere. For characterization of the turbulent properties we compare kinetic energy spectra on granular and sub-granular scales obtained from infrared TiO observations with the New Solar Telescope (Big Bear Solar Observatory) and from 3D radiative MHD numerical simulations ('SolarBox' code). We find that the numerical simulations require a high spatial resolution with 10 - 25 km grid-step in order to reproduce the inertial (Kolmogorov) turbulence range. The observational data require an averaging procedure to remove noise and potential instrumental artifacts. The resulting kinetic energy spectra show a good agreement between the simulations and observations, opening new perspectives for detailed joint analysis of more complex turbulent phenomena on the Sun, and possibly on other stars. In addition, using the simulations and observations we investigate effects of background magnetic field, which is concentrated in self-organized complicated structures in intergranular lanes, and find an increase of the small-scale turbulence energy and its decrease at larger scales due to magnetic field effects.
△ Less
Submitted 22 June, 2012;
originally announced June 2012.
-
Dynamics of Magnetized Vortex Tubes in the Solar Chromosphere
Authors:
I. N. Kitiashvili,
A. G. Kosovichev,
N. N. Mansour,
A. A. Wray
Abstract:
We use 3D radiative MHD simulations to investigate the formation and dynamics of small-scale (less than 0.5 Mm in diameter) vortex tubes spontaneously generated by turbulent convection in quiet-Sun regions with initially weak mean magnetic fields. The results show that the vortex tubes penetrate into the chromosphere and substantially affect the structure and dynamics of the solar atmosphere. The…
▽ More
We use 3D radiative MHD simulations to investigate the formation and dynamics of small-scale (less than 0.5 Mm in diameter) vortex tubes spontaneously generated by turbulent convection in quiet-Sun regions with initially weak mean magnetic fields. The results show that the vortex tubes penetrate into the chromosphere and substantially affect the structure and dynamics of the solar atmosphere. The vortex tubes are mostly concentrated in intergranular lanes and are characterized by strong (near sonic) downflows and swirling motions that capture and twist magnetic field lines, forming magnetic flux tubes that expand with height and which attain magnetic field strengths ranging from 200 G in the chromosphere to more than 1 kG in the photosphere. We investigate in detail the physical properties of these vortex tubes, including thermodynamic properties, flow dynamics, and kinetic and current helicities, and conclude that magnetized vortex tubes provide an important path for energy and momentum transfer from the convection zone into the chromosphere.
△ Less
Submitted 26 January, 2012;
originally announced January 2012.
-
Vortex Tubes of Turbulent Solar Convection
Authors:
I. N. Kitiashvili,
A. G. Kosovichev,
N. N. Mansour,
S. K. Lele,
A. A. Wray
Abstract:
Investigation of the turbulent properties of solar convection is extremely important for understanding the multi-scale dynamics observed on the solar surface. In particular, recent high-resolution observations have revealed ubiquitous vortical structures, and numerical simulations have demonstrated links between vortex tube dynamics and magnetic field organization and have shown the importance of…
▽ More
Investigation of the turbulent properties of solar convection is extremely important for understanding the multi-scale dynamics observed on the solar surface. In particular, recent high-resolution observations have revealed ubiquitous vortical structures, and numerical simulations have demonstrated links between vortex tube dynamics and magnetic field organization and have shown the importance of vortex tube interactions in the mechanisms of acoustic wave excitation on the Sun. In this paper we investigate the mechanisms of the formation of vortex tubes in highly-turbulent convective flows near the solar surface by using realistic radiative hydrodynamic LES simulations. Analysis of data from the simulations indicates two basic processes of vortex tube formation: 1) development of small-scale convective instability inside convective granules, and 2) a Kelvin-Helmholtz type instability of shearing flows in intergranular lanes. Our analysis shows that vortex stretching during these processes is a primary source of generation of small-scale vorticity on the Sun.
△ Less
Submitted 11 January, 2012; v1 submitted 27 December, 2011;
originally announced December 2011.
-
LES of turbulent convection in solar-type stars and formation of large-scale magnetic structures
Authors:
I. Rogachevskii,
N. Kleeorin,
I. N. Kitiashvili,
A. G. Kosovichev,
A. A. Wray,
N. N. Mansour
Abstract:
In this study we investigate the effects of turbulent convection on formation of large-scale inhomogeneous magnetic structures by means of Large-Eddy Simulation (LES) for convection in solar-type stars. The main idea of this study is the implementation of a new subgrid-scale model for the effective Lorentz force in a three-dimensional nonlinear radiative magnetohydrodynamics (MHD) code developed f…
▽ More
In this study we investigate the effects of turbulent convection on formation of large-scale inhomogeneous magnetic structures by means of Large-Eddy Simulation (LES) for convection in solar-type stars. The main idea of this study is the implementation of a new subgrid-scale model for the effective Lorentz force in a three-dimensional nonlinear radiative magnetohydrodynamics (MHD) code developed for simulating the upper solar convection zone and lower atmosphere. To this end we derived the energy budget equations, which include the effects of the subgrid-scale turbulence on the Lorentz-force, and implemented the new subgrid-scale turbulence model (TELF-Model) in a three-dimensional nonlinear MHD LES code. Using imposed initial vertical and horizontal uniform magnetic fields in LES with the TELF-Model, we have shown that the magnetic flux tubes formation is started when the initial mean magnetic field is larger than a threshold value (about 100 G). This is in agreement with the theoretical studies by Rogachevskii and Kleeorin (2007). We have determined the vertical profiles of the velocity and magnetic fluctuations, total MHD energy and anisotropy of turbulent magneto-convection, kinetic and current and cross helicities.
△ Less
Submitted 6 February, 2011;
originally announced February 2011.
-
Numerical MHD Simulations of Solar Magnetoconvection and Oscillations in Inclined Magnetic Field Regions
Authors:
I. N. Kitiashvili,
A. G. Kosovichev,
N. N. Mansour,
A. A. Wray
Abstract:
The sunspot penumbra is a transition zone between the strong vertical magnetic field area (sunspot umbra) and the quiet Sun. The penumbra has a fine filamentary structure that is characterized by magnetic field lines inclined toward the surface. Numerical simulations of solar convection in inclined magnetic field regions have provided an explanation of the filamentary structure and the Evershed ou…
▽ More
The sunspot penumbra is a transition zone between the strong vertical magnetic field area (sunspot umbra) and the quiet Sun. The penumbra has a fine filamentary structure that is characterized by magnetic field lines inclined toward the surface. Numerical simulations of solar convection in inclined magnetic field regions have provided an explanation of the filamentary structure and the Evershed outflow in the penumbra. In this paper, we use radiative MHD simulations to investigate the influence of the magnetic field inclination on the power spectrum of vertical velocity oscillations. The results reveal a strong shift of the resonance mode peaks to higher frequencies in the case of a highly inclined magnetic field. The frequency shift for the inclined field is significantly greater than that in vertical field regions of similar strength. This is consistent with the behavior of fast MHD waves.
△ Less
Submitted 30 November, 2010; v1 submitted 24 November, 2010;
originally announced November 2010.
-
Excitation of acoustic waves by vortices in the quiet Sun
Authors:
I. N. Kitiashvili,
A. G. Kosovichev,
N. N. Mansour,
A. A. Wray
Abstract:
Five-minutes oscillations is one of the basic properties of solar convection. Observations show mixture of a large number of acoustic wave fronts propagating from their sources. We investigate the process of acoustic waves excitation from the point of view of individual events, by using realistic 3D radiative hydrodynamic simulation of the quiet Sun. The results show that the excitation events are…
▽ More
Five-minutes oscillations is one of the basic properties of solar convection. Observations show mixture of a large number of acoustic wave fronts propagating from their sources. We investigate the process of acoustic waves excitation from the point of view of individual events, by using realistic 3D radiative hydrodynamic simulation of the quiet Sun. The results show that the excitation events are related to dynamics vortex tubes (or swirls) in the intergranular lanes. These whirlpool-like flows are characterized by very strong horizontal velocities (7 - 11 km/s) and downflows (~ 7 km/s), and are accompanied by strong decreases of the temperature, density and pressure at the surface and in a ~ 0.5-1 Mm deep layer below the surface. High-speed whirlpool flows can attract and capture other vortices. According to our simulation results, the processes of the vortex interaction, such as vortex annihilation, can cause the excitation of acoustic waves.
△ Less
Submitted 16 November, 2010;
originally announced November 2010.
-
Modeling of the subgrid-scale term of the filtered magnetic field transport equation
Authors:
G. Balarac,
A. G. Kosovichev,
O. Brugière,
A. A. Wray,
N. N. Mansour
Abstract:
Accurate subgrid-scale turbulence models are needed to perform realistic numerical magnetohydrodynamic (MHD) simulations of the subsurface flows of the Sun. To perform large-eddy simulations (LES) of turbulent MHD flows, three unknown terms have to be modeled. As a first step, this work proposes to use a priori tests to measure the accuracy of various models proposed to predict the SGS term appear…
▽ More
Accurate subgrid-scale turbulence models are needed to perform realistic numerical magnetohydrodynamic (MHD) simulations of the subsurface flows of the Sun. To perform large-eddy simulations (LES) of turbulent MHD flows, three unknown terms have to be modeled. As a first step, this work proposes to use a priori tests to measure the accuracy of various models proposed to predict the SGS term appearing in the transport equation of the filtered magnetic field. It is proposed to evaluate the SGS model accuracy in term of "structural" and "functional" performance, i.e. the model capacity to locally approximate the unknown term and to reproduce its energetic action, respectively. From our tests, it appears that a mixed model based on the scale-similarity model has better performance.
△ Less
Submitted 27 October, 2010;
originally announced October 2010.
-
Mechanism of spontaneous formation of stable magnetic structures on the Sun
Authors:
I. N. Kitiashvili,
A. G. Kosovichev,
A. A. Wray,
N. N. Mansour
Abstract:
One of the puzzling features of solar magnetism is formation of long-living compact magnetic structures; such as sunspots and pores, in the highly turbulent upper layer of the solar convective zone. We use realistic radiative 3D MHD simulations to investigate the interaction between magnetic field and turbulent convection. In the simulations, a weak vertical uniform magnetic field is imposed in a…
▽ More
One of the puzzling features of solar magnetism is formation of long-living compact magnetic structures; such as sunspots and pores, in the highly turbulent upper layer of the solar convective zone. We use realistic radiative 3D MHD simulations to investigate the interaction between magnetic field and turbulent convection. In the simulations, a weak vertical uniform magnetic field is imposed in a region of fully developed granular convection; and the total magnetic flux through the top and bottom boundaries is kept constant. The simulation results reveal a process of spontaneous formation of stable magnetic structures, which may be a key to understanding of the magnetic self-organization on the Sun and formation of pores and sunspots. This process consists of two basic steps: 1) formation of small-scale filamentary magnetic structures associated with concentrations of vorticity and whirlpool-type motions, and 2) merging of these structures due to the vortex attraction, caused by converging downdrafts around magnetic concentration below the surface. In the resulting large-scale structure maintained by the converging plasma motions, the magnetic field strength reaches ~1.5 kG at the surface and ~6 kG in the interior; and the surface structure resembles solar pores. The magnetic structure remains stable for the whole simulation run of several hours with no sign of decay.
△ Less
Submitted 13 April, 2010;
originally announced April 2010.
-
Explanation of the sea-serpent magnetic structure of sunspot penumbrae
Authors:
I. N. Kitiashvili,
L. R. Bellot Rubio,
A. G. Kosovichev,
N. N. Mansour,
A. Sainz Dalda,
A. A. Wray
Abstract:
Recent spectro-polarimetric observations of a sunspot showed the formation of bipolar magnetic patches in the mid penumbra and their propagation toward the outer penumbral boundary. The observations were interpreted as being caused by sea-serpent magnetic fields near the solar surface (Sainz Dalda & Bellot Rubio 2008). In this Letter, we develop a 3D radiative MHD numerical model to explain the…
▽ More
Recent spectro-polarimetric observations of a sunspot showed the formation of bipolar magnetic patches in the mid penumbra and their propagation toward the outer penumbral boundary. The observations were interpreted as being caused by sea-serpent magnetic fields near the solar surface (Sainz Dalda & Bellot Rubio 2008). In this Letter, we develop a 3D radiative MHD numerical model to explain the sea-serpent structure and the wave-like behavior of the penumbral magnetic field lines. The simulations reproduce the observed behavior, suggesting that the sea-serpent phenomenon is a consequence of magnetoconvection in a strongly inclined magnetic field. It involves several physical processes: filamentary structurization, high-speed overturning convective motions in strong, almost horizontal magnetic fields with partially frozen field lines, and traveling convective waves. The results demonstrate a correlation of the bipolar magnetic patches with high-speed Evershed downflows in the penumbra. This is the first time that a 3D numerical model of the penumbra results in downward directed magnetic fields, an essential ingredient of sunspot penumbrae that has eluded explanation until now.
△ Less
Submitted 26 February, 2010;
originally announced March 2010.
-
Traveling Waves of Magnetoconvection and the Origin of the Evershed Effect in Sunspots
Authors:
Irina N. Kitiashvili,
Alexander G. Kosovichev,
Alan A. Wray,
Nagi N. Mansour
Abstract:
Discovered in 1909 the Evershed effect represents strong mass outflows in sunspot penumbra, where the magnetic field of sunspots is filamentary and almost horizontal. These flows play important role in sunspots and have been studied in detail using large ground-based and space telescopes, but the basic understanding of its mechanism is still missing. We present results of realistic numerical sim…
▽ More
Discovered in 1909 the Evershed effect represents strong mass outflows in sunspot penumbra, where the magnetic field of sunspots is filamentary and almost horizontal. These flows play important role in sunspots and have been studied in detail using large ground-based and space telescopes, but the basic understanding of its mechanism is still missing. We present results of realistic numerical simulations of the Sun's subsurface dynamics, and argue that the key mechanism of this effect is in non-linear magnetoconvection that has properties of traveling waves in the presence of strong, highly inclined magnetic field. The simulations reproduce many observed features of the Evershed effect, including the high-speed "Evershed clouds", the filamentary structure of the flows, and the non-stationary quasi-periodic behavior. The results provide a synergy of previous theoretical models and lead to an interesting prediction of a large-scale organization of the outflows.
△ Less
Submitted 23 April, 2009;
originally announced April 2009.
-
Realistic Numerical Modeling of Solar Magnetoconvection and Oscillations
Authors:
Irina Kitiashvili,
Laetitia Jacoutot,
Alexander Kosovichev,
Alan Wray,
Nagi Mansour
Abstract:
We have developed 3D, compressible, non-linear radiative MHD simulations to study the influence of the magnetic field of various strength and geometry on the turbulent convective cells and on the excitation mechanisms of the acoustic oscillations. The results reveal substantial changes of the granulation structure with increased magnetic field, and a frequency-dependent reduction in the oscillat…
▽ More
We have developed 3D, compressible, non-linear radiative MHD simulations to study the influence of the magnetic field of various strength and geometry on the turbulent convective cells and on the excitation mechanisms of the acoustic oscillations. The results reveal substantial changes of the granulation structure with increased magnetic field, and a frequency-dependent reduction in the oscillation power. These simulation results reproduce the enhanced high-frequency acoustic emission observed at the boundaries of active region ("acoustic halo" phenomenon). In the presence of inclined magnetic field the solar convection develops filamentary structure with flows concentrated along magnetic filaments, and also exhibits behavior of running magnetoconvective waves, resembling recent observations of the sunspot penumbra dynamics from Hinode/SOT.
△ Less
Submitted 28 January, 2009; v1 submitted 27 January, 2009;
originally announced January 2009.
-
Solar convection and oscillations in magnetic regions
Authors:
L. Jacoutot,
A. G. Kosovichev,
A. Wray,
N. N. Mansour
Abstract:
The goal of this research is to investigate how magnetic field affects the dynamics of granular convection and excitation of solar oscillations by means of realistic numerical simulations. We have used a 3D, compressible, non-linear radiative magnetohydrodynamics code developed at the NASA Ames Research Center. This code takes into account several physical phenomena: compressible fluid flow in a…
▽ More
The goal of this research is to investigate how magnetic field affects the dynamics of granular convection and excitation of solar oscillations by means of realistic numerical simulations. We have used a 3D, compressible, non-linear radiative magnetohydrodynamics code developed at the NASA Ames Research Center. This code takes into account several physical phenomena: compressible fluid flow in a highly stratified medium, sub-grid scale turbulence models, radiative energy transfer between the fluid elements, and a real-gas equation of state. We have studied the influence of the magnetic field of various strength on the convective cells and on the excitation mechanisms of the acoustic oscillations by calculating spectral properties of the convective motions and oscillations. The results reveal substantial changes of the granulation structure with increased magnetic field, and a frequency-dependent reduction in the oscillation power in a good agreement with solar observations. These simulations suggest that the enhanced high-frequency acoustic emission at the boundaries of active region ("acoustic halo" phenomenon) is caused by the changes of the spatial-temporal spectrum of the turbulent convection in magnetic field, resulting in turbulent motions of smaller scales and higher frequencies than in quiet Sun regions.
△ Less
Submitted 24 May, 2008;
originally announced May 2008.
-
Numerical simulation of excitation of solar oscillation modes for different turbulent models
Authors:
L. Jacoutot,
A. G. Kosovichev,
A. Wray,
N. N. Mansour
Abstract:
The goal of this research is to investigate how well various turbulence models can describe physical properties of the upper convective boundary layer of the Sun. An accurate modeling of the turbulence motions is necessary for understanding the excitation mechanisms of solar oscillation modes. We have carried out realistic numerical simulations using several different physical Large Eddy Simulat…
▽ More
The goal of this research is to investigate how well various turbulence models can describe physical properties of the upper convective boundary layer of the Sun. An accurate modeling of the turbulence motions is necessary for understanding the excitation mechanisms of solar oscillation modes. We have carried out realistic numerical simulations using several different physical Large Eddy Simulation (LES) models (Hyperviscosity approach, Smagorinsky, and dynamic models) to investigate how the differences in turbulence modeling affect the damping and excitation of the oscillations and their spectral properties and compare with observations. We have first calculated the oscillation power spectra of radial and non-radial modes supported by the computational box with the different turbulence models. Then we have calculated the work integral input to the modes to estimate the influence of the turbulence model on the depth and strength of the oscillation sources. We have compared these results with previous studies and with the observed properties of solar oscillations. We find that the dynamic turbulence model provides the best agreement with the helioseismic observations.
△ Less
Submitted 11 October, 2007;
originally announced October 2007.