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 37 records · Page 2

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↗

Improved solar models constructed with a formulation of convection for stellar structure and evolution calculations without the mixing-length theory approximations

We have updated a previous attempt to incorporate within a solar model a treatment of convection based upon numerical simulations of convection rather than mixing-length theory (MLT). We have modified our formulation of convection for a better treatment of the kinetic energy flux. Our solar model has been updated to include a complete range of OPAL opacities, the Debye-Hueckel correction to the equation of state, helium diffusion due to gravitational settling, and atmospheres by Kurucz. We construct a series of models using both MLT and our revised formulation of convection and the compared results to measurements of the solar radius, the solar luminosity, and the depth of the solar convection zone as inferred from helioseismology. We find X(solar) = 0.702 +/- 0.005, Y(solar) = 0.278 +/- 0.005, and Z(solar) = 0.0193 +/- 0.0005.

Lydon, Thomas J.↗

Ultrafast High Accuracy PCRTM_SOLAR Model for Cloudy Atmosphere

An ultrafast high accuracy PCRTM_SOLAR model is developed based on PCA compression and principal component-based radiative transfer model (PCRTM). A fast algorithm for simulation of multi-scattering properties of cloud and/or aerosols is integrated into the fast infrared PCRTM. We completed radiance simulation and training for instruments, such as IASI, AIRS, CrIS, NASTI and SHIS, under diverse conditions. The new model is 5 orders faster than 52-stream DISORT with very high accuracy for cloudy sky radiative transfer simulation. It is suitable for hyperspectral remote data assimilation and cloudy sky retrievals.

Yang, Qiguang↗

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↗

Ayres' bifurcated solar model

The assumptions of Ayres' model of the upper solar atmosphere are examined. It is found that the bistable character of his model is postulated - through the assumptions concerning the opacity sources and the effect of mechanical waves, which are allowed to destroy the CO molecules but not to heat the gas. The neglect of cooling by metal lines is based on their reduced local cooling rate, but it ignores the increased depth over which this cooling occurs. Thus, the bifurcated model of the upper solar atmosphere consists of two models, one cold at the temperature minimum, with a kinetic temperature of 2900 K, and the other hot, with a temperature of 4900 K.

Kalkofen, W.↗

Solar models, neutrino experiments, and helioseismology

The event rates and their recognized uncertainties are calculated for 11 solar neutrino experiments using accurate solar models. These models are also used to evaluate the frequency spectrum of the p and g oscillations modes of the sun. It is shown that the discrepancy between the predicted and observed event rates in the Cl-37 and Kamiokande II experiments cannot be explained by a 'likely' fluctuation in input parameters with the best estimates and uncertainties given in the present study. It is suggested that, whatever the correct solution to the solar neutrino problem, it is unlikely to be a 'trival' error.

Bahcall, John N.↗

The response of the adiabatic exponent Gamma(1) to modifications of solar models

Guenther et al. (1989) investigated the effet of varying the abundance of elemental Ne on solar models and found that the changes in some aspects of these models in the ionization zones of helium and hydrogen were nonmonotonic functions of the assumed Ne abundance. In this paper, the results of Guenther et al. are examined, and it is found that the apparent nonmonotonic behavior of the models was caused by slight errors in the radii of their models, combined with the fact that the adiabatic exponent Gamma(1) was compared at fixed distance r from the center of the model. When the models are compared at fixed fractional radius r/R, the differences in the values of Gamma(1) were found to behave in a much more regular way.

Christensen-Dalsgaard, J.↗

Comparison of measured and modeled solar EUV flux and its effect on the E-F1 region ionosphere

Two different photochemical schemes are used to investigate the response of the E-F1 region ionosphere to different solar EUV flux models, and the results are compared with incoherent scatter radar electron density measurements taken at Millstone Hill. The latest EUV flux model (Tobiska, 1991), which incorporates more recent measurements, has generally more flux at short wavelengths compared to the Hinteregger et al. (1981) flux model based on AE-E satellite data. This results in better agreement with the measurements in the E-F1 region and above. The Tobiska flux model gives a smaller E region peak density, due to the influence of low Ly-beta flux in the November 10, 1988 rocket measurements of Woods and Rottman (1990). The photoionization and photoabsorption cross sections of Conway (1988) give results in somewhat better agreement with observations than the cross sections of Torr et al. (1979). For the zenith angles considered (daytime conditions), the Chapman function method for calculating photoabsorption yields results in satisfactory agreement with a more rigorous calculation method using a formula from Rees (1989).

Buonsanto, M. J.↗

{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↗

The effect of the Mihalas, Hummer, and Daeppen equation of state and the molecular opacity on the standard solar model

Improvements to the Yale Rotating Stellar Evolution Code (YREC) by incorporating the Mihalas-Hummer-Daeppen equation of state, an improved opacity interpolation routine, and the effects of molecular opacities, calculated at Los Alamos, have been made. the effect of each of the improvements on the standard solar model has been tested independently by computing the corresponding solar nonradial oscillation frequencies. According to these tests, the Mihalas-Hummer-Daeppen equation of state has very little effect on the model's low l p-mode oscillation spectrum compared to the model using the existing analytical equation of state implemented in YREC. On the other hand, the molecular opacity does improve the model's oscillation spectrum. The effect of molecular opacity on the computed solar oscillation frequencies is much larger than that of the Mihalas-Hummer-Daeppen equation of state. together, the two improvements to the physics reduce the discrepancy with observations by 10 microHz for the low l modes.

Kim, Y.-C.↗

Can binary stars test solar models?

The position in the H-R diagram of the approximately solar-mass component of the Hyades eclipsing binary, HD 27130, is compared with the predictions of stellar structure theory. The stellar models are calibrated by matching a model with the solar heavy element composition and age to the solar radius and luminosity. The comparison to the Hyades binary then is a test of the prediction that the initial solar luminosity was only about 0.7 times the present solar luminosity. The agreement is satisfactory, lending a measure of confidence to the solar model employed, provided that the initial helium abundance of the Hyades stars is not greater than that of the sun and is not less by more than about 0.03 in Y. Unless the model is grossly incorrect, the inference of Stromgren, Olsen, and Gustafsson (1982) from the 'Hyades anomaly' in intermediate-band photometry that Y(Hyades) is less than Y(solar) by 0.1 or 0.15 is rejected by the observed properties of HD 27130.

Popper, D. M.↗

Modeling solar coronal streamers

Coronal streamer models must now make the transition from research projects to applications tools. SOHO (Solar and Heliospheric Observatory) will be making measurements on plasma parameters in and around streamers and a quantitative model will be required for interpreting the data. The reason for this is that a streamer is inherently a magnetohydrodynamic phenomenon; its properties are determined both by the magnetic field and the dynamics of the plasma. The purpose of the model will be to analyze the energetics of the streamer and place the results into the context of the magnetic field. Current streamer modeling, and why it has taken so long for these models to develop into useful tools, is described and indications on where further development is needed, are given.

Suess, Steven T.↗

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↗

A Solar Model with g-Modes

Good evidence is assembled showing that the Suit's core arid surface vary on time scales from a month to a decade arid that a number of scales are similar. The most plausible source for numerous long time scales and periodicities is long-lived global oscillations. This suggests g-modes (oscillations restored mainly by buoyancy) because they particularly affect the core and base of the convective envelope, which then indirectly modulates the surface. Also, standing g-modes have rotational properties that match many observed periodicities. But the standard solar model (SSM) has a static core and excites few if any g-modes. making new interior structures worth exploring. The model outlined here assumes two well mixed shells near 0.18 and 0.68 R, (13 = solar radius) where sound speed data shows sharp deviations from the SSM. Mixing is sustained by flows driven by the oscillations. The shells form a cavity that excludes g-modes from their main damping region below 0.1 R, assisting their net excitation and increasing their oscillation periods by at least a factor of two and probably much more. In terms of the solar luminosity L, the modes transport up through the cavity a power approx. 0.004 L as a lower limit and 0.11 L as all upper limit. The modes dissipate energy in the outer shell and cool the inner shell, asymmetrically in each case, and this stimulates occasional convective events whose response time is typically 0.8 years longer near the inner shell. Such events cool the core and reduce neutrino flux while heating the envelope and increasing solar activity. This gives a physical basis for a well mixed Sun with low neutrino flux and basis for the observed anticorrelation and lag of neutrino behind surface activity.

Wolff, Charles L.↗

Solar models of low neutrino-counting rate - The depleted Maxwellian tail

Evolutionary sequences for the sun are presented which confirm that the Cl-37 neutrino counting rate will be greatly reduced if the high-energy tail of the Maxwellian distribution of relative energies is progressively depleted. Thermonuclear reaction rates and pressure are reevaluated for a distribution function modified by the correction factor suggested by Clayton (1974), and the effect of the results on solar models calculated with a simple Henyey code is discussed. It is shown that if the depletion is characterized by a certain exponential dependence on the distribution function, the counting rate will fall below 1 SNU for a distribution function of not less than 0.01. Suggestions are made for measuring the distribution function in the sun by means of neutrino spectroscopy and photography.

Clayton, D. D.↗