Search NASA⌕ Search

SEARCH · Search NASA

Results for “SOLAR MODELING”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

Solar Models in Light of New High Metallicity Measurements from Solar Wind Data

We study the impact of new metallicity measurements, from solar wind data, on the solar model. The “solar modeling problem” refers to the persisting discrepancy between helioseismological observations and predictions of solar models computed implementing state-of-the-art photospheric abundances. We critically reassess the problem, in particular considering the new set of abundances of von Steiger and Zurbuchen, determined through the in situ collection of solar wind samples from polar coronal holes. This new set of abundances indicates a solar metallicity Z{sub ⊙}⩾0.0196±0.0014, significantly higher than the currently established value. The new values hint at an abundance of volatile elements (i.e., C, N, O, Ne) close to previous results of Grevesse and Sauval, whereas the abundance of refractory elements (i.e., Mg, Si, S, Fe) is considerably increased. Using the Linear Solar Model formalism, we determine the variation of helioseismological observables in response to the changes in elemental abundances, in order to explore the consistency of these new measurements with constraints from helioseismology. We find that for observables that are particularly sensitive to the abundance of volatile elements, in particular the radius of the convective zone boundary (CZB) and the sound speed around the radius of CZB, improved agreement over previous models is obtained. Conversely, the high abundance of refractories correlates with a higher core temperature, resulting in an overproduction of neutrinos and a huge increase in the surface helium abundance. We conclude that the “solar modeling problem” remains unsolved.

79 ASTRONOMY AND ASTROPHYSICS↗

Solar Modeling and Opacities

In 2004, improved analyses of solar spectra resulted in a downward revision of the abundances of elements heavier than hydrogen and helium, in particular, the elements C, N, and O. Solar models evolved using these lower abundances showed discrepancies with inferences from solar oscillation observations (helioseismology). This problem has not been solved satisfactorily in the intervening 16 years. How serious/important is this problem? This talk will give an overview of how standard solar models are calculated and of the constraints from helioseismology. The talk will also present results of some attempts to change input physics or assumptions in the standard solar model, including opacity modifications, to try to resolve the discrepancies. I would like to have a discussion with the participants after the talk about which directions appear most promising to resolve the problem.

79 ASTRONOMY AND ASTROPHYSICS↗

FIRST NEW SOLAR MODELS WITH OPAS OPACITY TABLES

Stellar seismology appears more and more as a powerful tool for a better determination of the fundamental properties of solar-type stars. However, the particular case of the Sun is still challenging. For about a decade now, the helioseismic sound-speed determination has continued to disagree with the standard solar model (SSM) prediction, questioning the reliability of this model. One of the sources of uncertainty could be in the treatment of the transport of radiation from the solar core to the surface. In this Letter, we use the new OPAS opacity tables, recently available for solar modeling, to address this issue. We discuss first the peculiarities of these tables, then we quantify their impact on the solar sound-speed and density profiles using the reduced OPAS tables taken on the grids of the OPAL ones. We use the two evolution codes, Modules for Experiments in Stellar Astrophysics and Code Liégeois d’Evolution Stellaire, that led to similar conclusions in the solar radiative zone. In comparison to commonly used OPAL opacity tables, the new solar models are computed for the most recent photospheric composition with OPAS tables and present improvements to the location of the base of the convective zone and to the description of the solar radiative zone in comparison to the helioseismic observations, even if the differences in the Rosseland mean opacity do not exceed 6%. We finally carry out a comparison to a solar model computed with the OP opacity tables.

79 ASTRONOMY AND ASTROPHYSICS↗

An Expanded Set of Los Alamos OPLIB Tables in MESA: Type-1 Rosseland-mean Opacities and Solar Models

Abstract We present a set of 1194 Type-1 Rosseland-mean opacity tables for four different metallicity mixtures. These new Los Alamos OPLIB atomic radiative opacity tables are an order of magnitude larger in number than any previous opacity table release, and span regimes where previous opacity tables have not existed. For example, the new set of opacity tables expands the metallicity range to Z = 10 −6 to Z = 0.2, which allows improved accuracy of opacities at low and high metallicity, increases the table density in the metallicity range Z = 10 −4 to Z = 0.1 to enhance the accuracy of opacities drawn from interpolations across neighboring metallicities, and adds entries for hydrogen mass fractions between X = 0 and X = 0.1 including X = 10 −2 , 10 −3 , 10 −4 , 10 −5 , 10 −6 that can improve stellar models of hydrogen deficient stars. We implement these new OPLIB radiative opacity tables in MESA and find that calibrated solar models agree broadly with previously published helioseismic and solar neutrino results. We find differences between using the new 1194 OPLIB opacity tables and the 126 OPAL opacity tables range from ≈20% to 80% across individual chemical mixtures, up to ≈8% and ≈15% at the bottom and top of the solar convection zone respectively, and ≈7% in the solar core. We also find differences between standard solar models using different opacity table sources that are on par with altering the initial abundance mixture. We conclude that this new, open-access set of OPLIB opacity tables does not solve the solar modeling problem, and suggest the investigation of physical mechanisms other than the atomic radiative opacity.

79 ASTRONOMY AND ASTROPHYSICS↗

A New Generation of Standard Solar Models

We compute a new generation of standard solar models (SSMs) that includes recent updates on some important nuclear reaction rates and a more consistent treatment of the equation of state. Models also include a novel and flexible treatment of opacity uncertainties based on opacity kernels, required in light of recent theoretical and experimental works on radiative opacity. Two large sets of SSMs, each based on a different canonical set of solar abundances with high and low metallicity (Z), are computed to determine model uncertainties and correlations among different observables. We present detailed comparisons of high- and low-Z models against different ensembles of solar observables, including solar neutrinos, surface helium abundance, depth of the convective envelope, and sound speed profile. A global comparison, including all observables, yields a p-value of 2.7σ for the high-Z model and 4.7σ for the low-Z one. When the sound speed differences in the narrow region of 0.65<r/R{sub ⊙}<0.70 are excluded from the analysis, results are 0.9σ and 3.0σ for high- and low-Z models respectively. These results show that high-Z models agree well with solar data but have a systematic problem right below the bottom of the convective envelope linked to steepness of molecular weight and temperature gradients, and that low-Z models lead to a much more general disagreement with solar data. We also show that, while simple parametrizations of opacity uncertainties can strongly alleviate the solar abundance problem, they are insufficient to substantially improve the agreement of SSMs with helioseismic data beyond that obtained for high-Z models due to the intrinsic correlations of theoretical predictions.

79 ASTRONOMY AND ASTROPHYSICS↗

Application of New Los Alamos OPLIB Opacities in Solar Modeling Using the Mesa Code

The Sun provides the only local laboratory to test our understanding of stellar physics. The production, then, of valid dynamical models of the Sun is critical to how we think about other stars. Despite improvements in modeling methodologies and abundance measurements, solar models continue to reproduce errors in convection zone depth and surface helium abundance compared to helioseismic observations. Using the 1D stellar evolution code MESAstar, this work first evaluates the effect of changes to timestep resolution and mixing length theory on model outputs. A standard test case is then constructed with both GS98 and AGSS09 abundances using OP and OPAL opacities before implementing the new Los Alamos OPLIB opacities. It is shown that, when using GS98 abundances, OPLIB opacities improve agreement with helioseismic inference of the convection zone base radius, R czb , by 1.68-σ and 1.19-σ compared to OP and OPAL, respectively, while worsening agreement with the inferred envelope helium mass fraction, Y surf , by 0.63-σ and 0.98-σ for OP and OPAL.

14 SOLAR ENERGY↗

In-depth analysis of solar models with high-metallicity abundances and updated opacity tables

Context. As a result of the high-quality constraints available for the Sun, we are able to carry out detailed combined analyses using neutrino, spectroscopic, and helioseismic observations. These studies lay the ground for future improvements of the key physical components of solar and stellar models because ingredients such as the equation of state, the radiative opacities, or the prescriptions for macroscopic transport processes of chemicals are then used to study other stars in the Universe. Aims. We study the existing degeneracies in solar models using the recent high-metallicity spectroscopic abundances by comparing them to helioseismic and neutrino data and discuss the effect on their properties of changes in the micro and macro physical ingredients. Methods. We carried out a detailed study of solar models computed with a high-metallicity composition from the literature based on averaged 3D models that were claimed to resolve the solar modelling problem. We compared these models to helioseismic and neutrino constraints. Results. The properties of the solar models are significantly affected by the use of the recent OPLIB opacity tables and the inclusion of macroscopic transport. The properties of the standard solar models computed using the OPAL opacities are similar to those for which the OP opacities were used. We show that a modification of the temperature gradient just below the base of the convective zone is required to remove the discrepancies in solar models, particularly in the presence of macroscopic mixing. This can be simulated by a localised increase in the opacity of a few percent. Conclusions. We conclude that the existing degeneracies and issues in solar modelling are not removed by using an increase in the solar metallicity, in contradiction to what has been suggested in the recent literature. Therefore, standard solar models cannot be used as an argument for a high-metallicity composition. While further work is required to improve solar models, we note that direct helioseismic inversions indicate a low metallicity in the convective envelope, in agreement with spectroscopic analyses based on full 3D models.

74 ATOMIC AND MOLECULAR PHYSICS↗

SOLAR MODELS WITH NEW LOW METAL ABUNDANCES

In the past decade, the photospheric abundances of the Sun had been revised several times by many observers. The standard solar models constructed with the new low-metal abundances disagree with helioseismic results and detected neutrino fluxes. The solar model problem has puzzled some stellar physicists for more than 10 years. Rotation, enhanced diffusion, convection overshoot, and magnetic fields are used to reconcile the new abundances with helioseismology. The too low helium subsurface abundance in enhanced diffusion models can be improved by the mixing caused by rotation and magnetic fields. The problem of the depth of the convective zone in rotating models can be resolved by convection overshoot. Consequently, the Asplund–Grevesse–Sauval rotation model including overshooting (AGSR) reproduces the seismically inferred sound-speed and density profiles and the convection zone depth as well as the Grevesse and Sauval model computed before. But this model fails to reproduce the surface helium abundance, which is 0.2393 (2.6 σ away from the seismic value), and neutrino fluxes. The magnetic model called AGSM keeps the agreement of the AGSR and improves the prediction of the surface helium abundance. The observed separation ratios r {sub 02} and r {sub 13} are reasonably reproduced by AGSM. Moreover, neutrino fluxes calculated by this model are not far from the detected neutrino fluxes and the predictions of previous works.

79 ASTRONOMY AND ASTROPHYSICS↗

Assessing the WRF-Solar Model Performance Using Satellite-Derived Irradiance from the National Solar Radiation Database

Abstract WRF-Solar is a numerical weather prediction model specifically designed to meet the increasing demand for accurate solar irradiance forecasting. The model provides flexibility in the representation of the aerosol–cloud–radiation processes. This flexibility can be argued to make it more difficult to improve the model’s performance because of the necessity of inspecting different configurations. To alleviate this situation, WRF-Solar has a reference configuration to use as a benchmark in sensitivity experiments. However, the scarcity of high-quality ground observations is a handicap to accurately quantify the model performance. An alternative to ground observations are satellite irradiance retrievals. Herein we analyze the adequacy of the National Solar Radiation Database (NSRDB) to validate the WRF-Solar performance using high-quality global horizontal irradiance (GHI) observations across the contiguous United States (CONUS). Based on the sufficient performance of NSRDB, we further analyze the WRF-Solar forecast errors across the CONUS, the growth of the forecasting errors as a function of the lead time, and sensitivities to the grid spacing and the representation of the radiative effects of unresolved clouds. Our results based on WRF-Solar forecasts spanning 2018 reveal a 7% median degradation of the mean absolute error (MAE) from the first to the second daytime period. Reducing the grid spacing from 9 to 3 km leads to a 4% improvement in the MAE, whereas activating the radiative effects of unresolved clouds is desirable over most of the CONUS even at 3 km of grid spacing. A systematic overestimation of the GHI is found. These results illustrate the potential of GHI retrievals to contribute to increasing the WRF-Solar performance.

14 SOLAR ENERGY↗

{sup 26}Al production: The Allende meteorite (Chihuahua) stellar nucleosynthesis and solar models

In 1969 a meteorite fell near the small town of Allende, state of Chihuahua in the north of Mexico. Its study yielded information that changed the current understanding of the solar model. In particular traces of {sup 26}Al were found. Abundances of that isotope had been seen in the universe and were related to regions of active heavy nucleosynthesis. Its presence on the solar system was unexpected. It is now understood that cosmic rays induce nuclear reactions on materials to produce {sup 26}Al, on Earth this is well known and it is the basis of many environmental studies, so it is not only the product of some high metalicity star collapse. Taking advantage of the recently reinforced laboratory infrastructure of the Instituto de Física, at UNAM in Mexico City, we proposed to measure the cross section for {sup 26}Al production via some of the most likely reactions, from the nuclear physics point of view (highest Q-values). In this paper the study of the {sup 28}Si(d,α){sup 26} Al nuclear reaction is shown. A target is prepared by a mixture of silicon and aluminum powders. It is irradiated with a deuteron beam (≈1 µA current) at the MV CN-Van de Graaff accelerator laboratory. The number of projectiles is deduced by Rutherford Backscattering Spectrometry (RBS). The produced {sup 26}Al nuclei are then counted at the Accelerator Mass Spectrometry Laboratory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Influences of Cloud Microphysics on the Components of Solar Irradiance in the WRF-Solar Model

An accurate forecast of Global Horizontal solar Irradiance (GHI) and Direct Normal Irradiance (DNI) in cloudy conditions remains a major challenge in the solar energy industry. This study focuses on the impact of cloud microphysics on GHI and its partition into DNI and Diffuse Horizontal Irradiance (DHI) using the Weather Research and Forecasting model specifically designed for solar radiation applications (WRF-Solar) and seven microphysical schemes. Three stratocumulus (Sc) and five shallow cumulus (Cu) cases are simulated and evaluated against measurements at the US Department of Energy’s Atmospheric Radiation Measurement (ARM) user facility, Southern Great Plains (SGP) site. Results show that different microphysical schemes lead to spreads in simulated solar irradiance components up to 75% and 350% from their ensemble means in the Cu and Sc cases, respectively. The Cu cases have smaller microphysical sensitivity due to a limited cloud fraction and smaller domain-averaged cloud water mixing ratio compared to Sc cases. Cloud properties also influence the partition of GHI into DNI and DHI, and the model simulates better GHI than DNI and DHI due to a non-physical error compensation between DNI and DHI. The microphysical schemes that produce more accurate liquid water paths and effective radii of cloud droplets have a better overall performance.

54 ENVIRONMENTAL SCIENCES↗

CORONAL JETS SIMULATED WITH THE GLOBAL ALFVÉN WAVE SOLAR MODEL

This paper describes a numerical modeling study of coronal jets to understand their effects on the global corona and their contribution to the solar wind. We implement jets into a well-established three-dimensional, two-temperature magnetohydrodynamic (MHD) solar corona model employing Alfvén-wave dissipation to produce a realistic solar-wind background. The jets are produced by positioning a compact magnetic dipole under the solar surface and rotating the boundary plasma around the dipole's magnetic axis. The moving plasma drags the magnetic field lines along with it, ultimately leading to a reconnection-driven jet similar to that described by Pariat et al. We compare line-of-sight synthetic images to multiple jet observations at EUV and X-ray bands, and find very close matches in terms of physical structure, dynamics, and emission. Key contributors to this agreement are the greatly enhanced plasma density and temperature in our jets compared to previous models. These enhancements arise from the comprehensive thermodynamic model that we use and, also, our inclusion of a dense chromosphere at the base of our jet-generating regions. We further find that the large-scale corona is affected significantly by the outwardly propagating torsional Alfvén waves generated by our polar jet, across 40° in latitude and out to 24 R {sub ⊙}. We estimate that polar jets contribute only a few percent to the steady-state solar-wind energy outflow.

79 ASTRONOMY AND ASTROPHYSICS↗

Jackson County Community Solar Modeling Results [Slides]

Jackson County, via the U.S. Department of Energy's Communities LEAP (Local Energy Action Program) requested support on the topic of Community Solar in their area. NREL completed a technical and financial modeling exercise to determine the potential for developing community solar in the region. The analysis performed determined that the current policy and economic environment, in large part due to Illinois's Solar For All program, is advantageous to the deployment and growth of community solar in Illinois. Across all 6 scenarios modeled, including current minimum and maximum financial incentives, the economics yielded positive net present values for both the project owner and its subscribers. The analysis shown in this presentation is only one of the first steps in deploying community solar and suggests Jackson County should pursue next steps. Such actions could include refining project economics, partnering with local utilities, policy makers, and developers to identify potential projects and outreach opportunities, and educating and engaging community members.

14 SOLAR ENERGY↗

pvlib python: a python package for modeling solar energy systems

pvlib python is a community-supported open source tool that provides a set of functions and classes for simulating the performance of photovoltaic energy systems. pvlib python aims to provide reference implementations of models relevant to solar energy, including for example algorithms for solar position, clear sky irradiance, irradiance transposition, DC power, and DC-to-AC power conversion. pvlib python is an important component of a growing ecosystem of open source tools for solar energy (William F. Holmgren, Hansen, Stein, & Mikofski, 2018).

14 SOLAR ENERGY↗

Is There A Solar Model Solution to The Faint Young Sun Paradox? [Slides]

We have calculated possible solar rotation histories and determined reasonable estimates for the total mass shed off by the star over its lifetime. We predict that all solar-type stars spin down very early in their evolution, are therefore inactive, and as a result have low mass loss. The Sun was probably 0.1−0. 2% more massive than it is today. Better solutions to the Faint Young Sun Paradox, from more advanced models of early Earth and Mars, are needed.

79 ASTRONOMY AND ASTROPHYSICS↗

Modeling solar wind with boundary conditions from interplanetary scintillations

Interplanetary scintillations make it possible to create three-dimensional, time- dependent distributions of the solar wind velocity. Combined with the magnetic field observations in the solar photosphere, they help perform solar wind simulations in a genuinely time-dependent way. Interplanetary scintillation measurements from the Ooty Radio Astronomical Observatory in India provide directions to multiple stars and may assure better resolution of transient processes in the solar wind. In this paper, we present velocity distributions derived from Ooty observations and compare them with those obtained with the Wang-Sheeley-Arge (WSA) model. We also present our simulations of the solar wind flow from 0.1 AU to 1 AU with the boundary conditions based on both Ooty and WSA data.

79 ASTRONOMY AND ASTROPHYSICS↗

Solar Models with Dynamic Screening and Early Mass Loss Tested by Helioseismic, Astrophysical, and Planetary Constraints

The faint young Sun paradox remains an open question. Presented here is one possible solution to this apparent inconsistency, a more massive early Sun. Based on the conditions on early Earth and Mars, a luminosity constraint is set as a function of mass. For this, we examine helioseismic constraints of these alternative mass-losing models. Additionally, we explore a dynamic electron screening correction in an effort to improve helioseismic agreement in the core of the early mass-losing model.

79 ASTRONOMY AND ASTROPHYSICS↗

MODELING SOLAR CORONAL BRIGHT-POINT OSCILLATIONS WITH MULTIPLE NANOFLARE HEATED LOOPS

Intensity oscillations of coronal bright points (BPs) have been studied for the past several years. It has been known for a while that these BPs are closed magnetic loop-like structures. However, the initiation of such intensity oscillations is still an enigma. There have been many suggestions to explain these oscillations, but so far modeling such BPs has not been explored. Using a multithreaded nanoflare heated loop model we study the behavior of such BPs in this work. We compute typical loop lengths of BPs using potential field-line extrapolation of available data, and set this as the length of our simulated loops. We produce intensity-like observables through forward modeling and analyze the intensity time series using wavelet analysis, as was done by previous observers. The result reveals similar intensity oscillation periods reported in past observations. It is suggested these oscillations are actually shock wave propagations along the loop. We also show that if one considers different background subtractions, one can extract adiabatic standing modes from the intensity time-series data as well, both from the observed and simulated data.

79 ASTRONOMY AND ASTROPHYSICS↗