Search NASA⌕ Search

SEARCH · Search NASA

Results for “solar convection zone”

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.

16 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Radiation burnthrough measurements to infer opacity at conditions close to the solar radiative zone–convective zone boundary

Recent measurements at the Sandia National Laboratory of the x-ray transmission of iron plasma have inferred opacities much higher than predicted by theory, which casts doubt on modeling of iron x-ray radiative opacity at conditions close to the solar convective zone-radiative zone boundary. An increased radiative opacity of the solar mixture, in particular iron, is a possible explanation for the disagreement in the position of the solar convection zone-radiative zone boundary as measured by helioseismology and predicted by modeling using the most recent photosphere analysis of the elemental composition. In this report we present data from radiation burnthrough experiments, which do not support a large increase in the opacity of iron at conditions close to the base of the solar convection zone and provide a constraint on the possible values of both the mean opacity and the opacity in the x-ray range of the Sandia experiments. The data agree with opacity values from current state-of-the-art opacity modeling using the CASSANDRA opacity code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

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↗
DOE OSTI2021

Quantitative comparison of opacities calculated using the R -matrix and distorted-wave methods: Fe xvii

ABSTRACT We present here a detailed calculation of opacities for Fe xvii at the physical conditions corresponding to the base of the Solar convection zone. Many ingredients are involved in the calculation of opacities. We review the impact of each ingredient on the final monochromatic and mean opacities (Rosseland and Planck). The necessary atomic data were calculated with the R-matrix and the distorted-wave (DW) methods. We study the effect of broadening, of resolution, of the extent of configuration sets and of configuration interaction to understand the differences between several theoretical predictions as well as the existing large disagreement with measurements. New Dirac R-matrix calculations including all configurations up to the n = 4, 5, and 6 complexes have been performed as well as corresponding Breit–Pauli DW calculations. The DW calculations have been extended to include autoionizing initial levels. A quantitative contrast is made between comparable DW and R-matrix models. We have reached self-convergence with n = 6 R-matrix and DW calculations. Populations in autoionizing initial levels contribute significantly to the opacities and should not be neglected. The R-matrix and DW results are consistent under the similar treatment of resonance broadening. The comparison with the experiment shows a persistent difference in the continuum while the filling of the windows shows some improvement. This study defines our path to the next generation of opacities and opacity tables for stellar modelling.

Delahaye, F.↗
DOE OSTI2024

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↗
DOE OSTI2022

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↗
DOE OSTI2026

Evidence of 5 minute Oscillations from Parker Solar Probe

The Sun’s surface vibrates in characteristic 5 minute oscillations, known as p-modes, generated by sound waves trapped within the convection zone. Although these oscillations have long been hypothesized to reach into the solar wind, direct in situ evidence has remained elusive, even during previous close encounters by Parker Solar Probe (PSP). Here, we present strong in situ evidence of 5 minute oscillations in the upper solar corona, based on observations from PSP’s three closest perihelia. In two events at 9.9 solar radii (R ⊙ ), we identify statistically significant (∼6σ) 3.1–3.2 mHz peaks in the magnetic field power spectrum, each appearing as a large-amplitude, spherically polarized Alfvénic wave train lasting approximately 35 minutes. These results demonstrate that global solar oscillations can reach the solar wind.

Fast solar wind↗
DOE OSTI2021

Detection Limits of Low-mass, Long-period Exoplanets Using Gaussian Processes Applied to HARPS-N Solar Radial Velocities

Radial velocity (RV) searches for Earth-mass exoplanets in the habitable zone around Sun-like stars are limited by the effects of stellar variability on the host star. In particular, suppression of convective blueshift and brightness inhomogeneities due to photospheric faculae/plage and starspots are the dominant contribution to the variability of such stellar RVs. Gaussian process (GP) regression is a powerful tool for statistically modeling these quasi-periodic variations. We investigate the limits of this technique using 800 days of RVs from the solar telescope on the High Accuracy Radial velocity Planet Searcher for the Northern hemisphere (HARPS-N) spectrograph. These data provide a well-sampled time series of stellar RV variations. Into this data set, we inject Keplerian signals with periods between 100 and 500 days and amplitudes between 0.6 and 2.4 m s{sup −1}. We use GP regression to fit the resulting RVs and determine the statistical significance of recovered periods and amplitudes. We then generate synthetic RVs with the same covariance properties as the solar data to determine a lower bound on the observational baseline necessary to detect low-mass planets in Venus-like orbits around a Sun-like star. Our simulations show that discovering planets with a larger mass (∼0.5 m s{sup −1}) using current-generation spectrographs and GP regression will require more than 12 yr of densely sampled RV observations. Furthermore, even with a perfect model of stellar variability, discovering a true exo-Venus (∼0.1 m s{sup −1}) with current instruments would take over 15 yr. Therefore, next-generation spectrographs and better models of stellar variability are required for detection of such planets.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Development of improved higher-order correction for the NIF opacity spectrometer

X-ray opacity measurements on the National Ignition Facility (NIF) are in the process of reproducing earlier measurements from the Sandia Z Facility, in particular for oxygen and iron plasmas. These measurements have the potential to revise our understanding of the “solar problem” and of the hot degenerate Q class white dwarf structure by probing plasma conditions near the base of their convection zones. Accurate opacity measurements using soft x-ray Bragg crystal spectrometers require correction for higher-order diffraction effects. Extending prior work in this area [Dutra et al., Review of Scientific Instruments 93, 113527 (2022)], we have developed a new method to remove higher-order spectral components from NIF opacity spectrometer data. Furthermore, by modeling absorption and backlighting continuum spectra and subtracting the second- and third-order components from the measured data, we are able to perform this correction while avoiding imprinting first-order model line features onto the data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Magnetized Winds of M-type Stars and Star–Planet Magnetic Interactions: Uncertainties and Modeling Strategy

M-type stars are the most common stars in the Universe. They are ideal hosts for the search of exoplanets in the habitable zone (HZ), as their small size and low temperature make the HZ much closer-in than their solar twins. Harboring very deep convective layers, they also usually exhibit very intense magnetic fields. Understanding their environment, in particular their coronal and wind properties, is thus very important, as they might be very different from what is observed in the solar system. The mass-loss rate of M-type stars is poorly known observationally, and recent attempts to estimate it for some of them (e.g., TRAPPIST-1 and Proxima Centauri) can vary by an order of magnitude. In this work, we revisit the stellar wind properties of M dwarfs in the light of the latest estimates of $\dot{M}$ through Lyα absorption at the astropause and slingshot prominences. We outline a modeling strategy to estimate the mass-loss rate, radiative loss, and wind speed, with uncertainties, based on an Alfvén-wave-driven stellar wind model. We find that it is very likely that several TRAPPIST-1 planets lie within the Alfvén surface, which implies that these planets experience star–planet magnetic interactions (SPMIs). We also find that SPMIs between Proxima Cen b and its host star could be the reason for recently observed radio emissions.

M stars↗
DOE OSTI2020

Axion helioscopes as solar magnetometers

Axion helioscopes search for solar axions and axionlike particles via inverse Primakoff conversion in strong laboratory magnets pointed at the Sun. Anticipating the detection of solar axions, in this work, we determine the potential for the planned next-generation helioscope, the International Axion Observatory (IAXO), to measure or constrain the solar magnetic field. To do this we consider a previously neglected component of the solar axion flux at sub-keV energies arising from the conversion of longitudinal plasmons. This flux is sensitively dependent to the magnetic field profile of the Sun, with lower energies corresponding to axions converting into photons at larger solar radii. If the detector technology eventually installed in IAXO has an energy resolution better than 200 eV, then solar axions could become an even more powerful messenger than neutrinos of the magnetic field in the core of the Sun. For energy resolutions better than 10 eV, IAXO could access the inner 70% of the Sun and begin to constrain the field at the tachocline: the boundary between the radiative and convective zones. The longitudinal plasmon flux from a toroidal magnetic field also has an additional 2% geometric modulation effect which could be used to measure the angular dependence of the magnetic field.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

On the Importance of the Convective Urca Process in 3D Simulations of a Simmering White Dwarf

Type Ia supernovae are bright thermonuclear explosions that are important to numerous areas of astronomy. However, the origins of these events are poorly understood. One proposed setting is that of a near Chandrasekhar mass white dwarf that undergoes runaway carbon burning in the core. During the thousand years leading up to the explosion, the white dwarf undergoes a simmering phase where slow carbon burning heats the core and drives convection. A poorly understood aspect of this phase is the convective Urca process, which links convection with weak nuclear reactions. We use the low Mach number code MAESTROeX to perform full 3D simulations as is required to accurately capture the turbulent convection. We present simulations with and without the A=23 convective Urca process, which have relaxed to a steady state. We characterize the effects of the convective Urca process on the neutrino losses, the nuclear energy generation, and the convective boundary. We find that the size of the convection zone is substantially reduced by the convective Urca process, though convection still extends past the Urca shell. Our findings on the structure of the convective zone and the compositional changes can be used to inform 1D stellar models that track the longer-timescale evolution.

FOS: Physical sciences↗
DOE OSTI2023

A reduced-order modeling of a tubular solar reactor for long duration thermochemical energy storage

The storage of solar energy in a solid form, referred to as a “solar fuel”, can be achieved through a process known as endothermic solar thermochemistry. This process transforms the absorbed solar energy into a stable and retrievable form that can be stored for extended periods of time. This paper presents a low–order heat transfer model of a counter–current tubular falling bed reactor designed to produce thermally reduced magnesium manganese oxide pellets for long duration thermochemical energy storage. The energy required for the endothermic reduction was supplied by concentrated solar energy or renewable electricity via indirect heating of the gas and solid reactants flowing in a ceramic tube. The counter-current gas flow enhances the mixing of the solid particles with the heat recuperation zone, allowing the gas and particles to enter and exit the tubular reactor close to room temperature. Further, the reactor was vertically oriented and was heated circumferentially by an adjustable level heat flux along a finite segment of its length. The temperature distribution of the reactor in response to transient changes along the tube was modeled by considering conduction, convection, and radiation heat transfer. Governing equations for the heat transfer model were solved by discretizing the reactor tube into a finite number of control volumes and using an energy balance for the heat exchange between the reactor wall, gas, and particles within the control volume. The energy absorbed during this endothermic reaction was modeled numerically by fitting the data of the chemical conversion rate with the corresponding temperature of particles in the heating zone. The numerical model has been experimentally validated using a reactor prototype made of a 121.92 cm alumina tube heated by a 7kW electric tube–furnace. The alumina tube receives magnesium manganese oxide pellets of 3.66±0.516 mm in diameter from the top, and a counter–current gas flow from the bottom. The reactor wall temperature was monitored by six thermocouples installed along the reactor tube length. The experimental procedure was numerically simulated, and the temperature variation along the reactor tube was compared with a matrix of experimental runs for a range of particles mass flowrates (0.75–1.25g/s) and corresponding gas flowrates (36–65 SLPM). The reactor system was heated gradually from room temperature to a steady state temperature of 1673K, and then cooled down to room temperature. The heating and cooling processes were simulated, and the numerical and experimental results were compared throughout processes. The numerical model showed similar trends to the experimental results, with an error of 0.69 to 7.9% for the particle inlet and 0.7 to 7.9% for the gas inlet during steady-state operation. The proposed numerical model can be implemented as a simplified physical model to design a feedback control system to regulate reactor temperature.

14 SOLAR ENERGY↗
DOE OSTI2026

Solar Heat for Industrial Processes: Integration with Chemical Reactors

The integration of solar thermal systems with chemical reactors has been proposed as part of a larger effort to develop and deploy solar heat for industrial processes (SHIP) technologies. A strong motivation for SHIP processes and technologies is the potential for high thermal efficiency coupled with low-cost thermal energy storage (TES) which can enable commercial deployment of such systems. While there are different ways to categorize SHIP technologies, one important such distinction is between directly irradiated systems and indirect off-sun process driven by a SHIP system. While directly irradiated systems can provide high thermal efficiencies and high fluxes, they usually require complex engineering solutions due to the need for redesigning the established processes and unit operations. In most cases, it is also more challenging to couple such a process to a TES system, losing some of the benefits of SHIP. On the other hand, using a SHIP system to drive an industrial process off-sun can allow better integration with existing process chains, easier TES capabilities, and potential for more applications fitting a specific SHIP technology. However, the integration of SHIP systems with the industrial processes is not fully explored in detail, especially in the case of high-temperature processes such as reforming, cracking, cement manufacturing, and iron/steelmaking. Many of these systems require heating fluxes of >50 kW/m^2, supplied via combustion of hydrocarbons in a fire box and benefitting from radiative heat transfer between the flue gases and the reaction zones. As such, using SHIP systems for such processes is more complex than providing the same thermal input in the form of a heat transfer medium (HTM) entering the reactor, kiln, or furnace. Moreover, in case convective heat transfer using SHIP is envisioned, for example using supercritical CO2 as the HTM from a particle receiver, the thermal integration might be more challenging than initially envisioned: lower heat transfer coefficients and limited approach temperature might require large flow rates, causing a mismatch between the process thermal requirements and the thermal capacity of the SHIP system. In addition, even if the heat exchange between SHIP and reactor is effective, there is still a cold leg HTM at the reaction temperature or slightly below it. Chemical plants usually include a set of heat exchangers, heat recovery steam generators, and even power generation units - in a tightly integrated design - to recover the flue gases which are eventually vented. With SHIP systems mostly operating on a closed HTM loop, bottoming the cold leg is crucial. In this talk, we will present different modeling results for a variety of syngas production reactions, using catalytic and chemical looping processes, and discuss some of the challenges and design considerations for off-sun chemical reactors using SHIP systems.

14 SOLAR ENERGY↗
DOE OSTI2025

Multidimensional Nova Simulations with an Extended Buffer and Lower Initial Mixing Temperatures

A classical nova is a thermonuclear runaway initiated on a white dwarf accreting solar-like material from its stellar companion. Once the white dwarf accretes enough mass, the pressure at the base of the accreted layer reaches a critical point, leading to the ignition of the hydrogen fuel at their interface. This paper presents a set of two-dimensional CO classical nova simulations with an extended buffer zone of a fixed low density and temperature between the top of the accreted layer and the upper boundary, allowing us to capture the thermonuclear outburst in the domain. Our domain reduces the role of the upper-outflow boundary condition that has affected previous simulations and allows us to explore the nucleosynthesis evolution in detail. We also study the effects of the initial temperature perturbation and buffer size to explore their sensitivity in our simulations. Finally, we start our simulations with a lower temperature at the base of the accreted layer ( 7 × 10 7 K ) than previous work, allowing us to capture mixing earlier in the evolution, reducing the effects of the mixing-length-theory assumptions. This allows for a more realistic description of convective transport in our models.

Smith Clark, Alexander (ORCID:0000000159611680)↗
DOE OSTI2020

A Survey of Multizone Airflow Simulation Software

H-Canyon simulates the building and process ventilation to ensure positive contamination control for any modifications as well as evaluating off-normal conditions. Simulations can also be used to identify air quality of personnel occupied areas. The models can be classified into two groups, single and multi-zone models. As the name states, a single zone model assumes that the building can be described by a single, well-mixed zone. Whereas large buildings have structures that must be described as multizone. Multizone airflow models deal with the difficulty of flows in a building by distinguishing the effects of internal flow restrictions. Several models have been developed that take internal partitions into account. CONTAM (currently used) is a multizone airflow program created by NIST to determine airflows including infiltration, exfiltration, wind pressures acting outside the building, buoyancy effects caused by the indoor and outdoor air temperature difference. COMIS is a multizone air infiltration, ventilation, and contaminant transport model that was developed at Lawrence Berkeley National Laboratory. COMIS deals with applied pressure distribution, airflow components, HVAC systems, schedules, and contaminant transport. EnergyPlus is designed by DOE. It supports flexible component level configuration of HVAC models and simulates sub-hourly time steps to control fast dynamic systems. The program features include providing accurate environmental performance data including energy consumption and carbon emissions, and reporting solar gains on surfaces. Heavent is used to model industrial ventilation systems to control airborne contaminants. It's designed based on the American Conference of Governmental Industrial Hygienists (ACGIH) Ventilation manual. DesignBuilder is the most recognized and advanced user interface to EnergyPlus. It gives access to all of the most frequently required simulation capabilities covering building fabric, thermal mass, glazing, shading renewables, and HVAC analysis. The objective of this study is to identify all relevant software available including commercially licensed, freely licensed, and open source, that can be used to simulate large, complex process facility ventilation systems. All discovered software tools are evaluated for salient features and compared to the software tool currently used to simulate the H-Canyon ventilation system. The design criteria of a ventilation system has two basic requirements: Supplying fresh air for the occupants and changing the air in the room sufficiently; therefore, the smells, fumes, and contaminants will be removed. In this study, first, whether the software simulates airflow (infiltration, natural ventilation, and interzonal) or considers them to be inputs is discussed. Second, if indoor contaminant levels are considered; and finally suggestions on how airflow analysis can be more appropriately accounted for when performing airflow simulation of a building. All of the airflow software reviewed for this study can account for constant infiltration and natural ventilation rates that are not affected by changes in indoor and outdoor conditions. Modelling infiltration in EnergyPlus or DesignBuilder are applicable for small buildings and not appropriate for larger buildings like the H-Canyon. The effect of wind on external pressures, infiltration, and natural and mechanical ventilation, can also be calculated using all of these software. These programs have control options to switch between natural to mechanical ventilation modes or both (hybrid). Interzonal airflow can be applied which means one zone can be specified to receive a fixed amount of air from another zone. The following ventilation modeling software is recommended for further testing: CONTAM, COMIS, Heavent. The next step in evaluating the software will consist of creating test facility simulations using evaluation copies of the software.

42 ENGINEERING↗