SEARCH · Search NASA
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.
This core model proposal updates modeling structures, key technology assumptions, and data sources for a detailed representation of hydrogen in GCAM’s energy system, including production, transmission and distribution, and consumption in hydrogen end-use sectors (buildings, transportation, industry).
The goal of this core model proposal (CMP) is to break out Ukraine from the Europe_Eastern region. This work aims to establish Ukraine as an independent region in the GCAM core (region 14) while moving Belarus and Moldova to region 15 (Europe_Non_EU). We have: 1) Updated several mappings to recode region 14 (formerly Europe_Eastern) as Ukraine and moved Belarus and Moldova to region 15 (Europe_Non_EU); 2) Updated several assumptions in the raw data files which provide information by region to reframe Ukraine as the 14th region, including coefficients, base year values, share weight interpolation values and rules, pipeline networks for gas trade, elasticities, shares, etc. 3) Changed documentation and in-code comments at several places referring to fixed 32 regions in GCAM to indicate that GCAM can have any number of regions; 4) Updated code base in gcamdata to dynamically process data for Ukraine given special cases.
This Core Model Proposal (CMP) comprehensively restructures and updates the macroeconomic and socioeconomic modules in gcamdata. It includes visualizations of key data inputs, accounting identities, and data flows in the context of GCAM-Macro-KLEM. Major improvements include: (1) updating the Penn World Table (PWT) to version 10 and incorporating a new source, the Global Macro Database (GMD); (2) updating the SSP socioeconomics database from version 3.0.1 to 3.2; (3) introducing SSP-specific differentiation of employment and labor force data; (4) improving data integration between national accounts (from PWT, GMD, and GTAP) and GDP/population data from external sources; and (5) general data cleaning and structural refinements. We document the data sources and key assumptions used throughout the processing. These updates establish the foundation for the forthcoming KLEAM version of GCAM-macro.
This Core Model Proposal (CMP) updates the Shared Socioeconomic Pathway (SSP) database to a recent version (v3.0.1; released in 2024) within GCAM. Currently, GCAM relies on socioeconomic drivers, including population, GDP, and labor force projections, from the original SSP database version released in 2013. These projections, provided by independent socioeconomic dynamic models (e.g., multi-dimensional demographic models and macroeconomic models of convergence growth), may need regular updates when (1) near-term observations become available and (2) there are updates and advancements in the socioeconomic modeling. Timely updates of socioeconomic drivers in global economic equilibrium and multisector dynamic modeling will ensure (1) alignment of historical years and near-term projections with observations, enhancing base year calibrations, including calibration parameters and labor productivity, and (2) improvement of long-term projections with updated socioeconomic drivers, which set the scale of the economy. This CMP updates the SSP data (from v2013 to v2024) and also fixes/reconciles historical GDP data sources in GCAM. We investigate the impact of these updates on GCAM projections.
GCAM's current base year – the last historical year used for calibration – 2015, is considerably lagged compared to the present year, indicating out-of-date assumptions about near-past state of the resources, commodities, and technologies represented in GCAM. This proposal makes data and code changes to enable updating the GCAM base year easier by making the processing code more flexible and robust.
Hector V3.2.0 is the version documented in Dorheim et. al (accepted in GMD), the changes between the previous version coupled with GCAM were in response to the reviewer feedback. We corrected aerosol forcing coefficients based on Zelinka et al. (2023), enabled the permafrost module to be on by default, and recalibrated the model. These changes mean that we had to update the hector-gcam.ini file, it also causes some changes in Hector output behavior (described below) which may have implications on GCAM runs. Ultimately Hector is cooler by about 0.15 degrees, although this is scenario dependent.
The lunar model proposed helps to account for the offset of the center of gravity from the center of the optical figure, the moments of inertia of the Moon, the 'mascons,' the localization of the maria basins on the near side of the Moon, the igneous nature of rocks, and the remanent magnetism. In the proposed model the Moon has a core whose center is offset from the center of the outside spheroid towards the earth. Such a core will be formed if the Moon were entirely molten at some time in its past, and on solidification was synchronous with the earth.
A model is proposed for the formation and evolution of the Earth's core which provides an adequate energy source for maintaining the geodynamo. A modified inhomogeneous accretion model is proposed which leads to initial iron and refractory enrichment at the center of the planet. The probable heat source for melting of the core is the decay of Al. The refractory material is emplaced irregularly in the lowermost mantle with uranium and thorium serving as a long lived heat source. Fluid motions in the core are driven by the differential heating from above and the resulting cyclonic motions may be the source of the geodynamo.
Equations providing numerical values of the geomagnetic field spherical harmonic spatial power spectrum as defined by Lowes (1966, 1974) are obtained and this power spectrum is related to various other power spectra. Equations relating the spherical harmonic spatial power spectrum to average great circle power spectra for components of the vector magnetic field in the radial direction, along the great circle track and perpendicular to the first two directions are derived under the assumption that the sources of the field are internal. A statistical model for the crustal and core geomagnetic fields is proposed and used to derive equations for the expected main and crustal spherical harmonic power spectra. The model equations are then compared with observations to determine a scale factor which is then used to obtain an estimate for the core radius and a great circle power spectrum for the field component perpendicular to the great circle and radial directions which are in good agreement with observations. The predicted spherical harmonic power spectrum for the crustal field is found to be consistent with POGO satellite and aircraft data. Other possible models for the crustal and core geomagnetic fields are also briefly considered.
A model for the formation and evolution of the earth's core, which provides an adequate energy source for maintaining the geodynamo, is proposed. A modified inhomogeneous accretion model is proposed which leads to initial iron and refractory enrichment at the center of the planet. The probable heat source for melting of the core is the decay of Al-26. The refractory material is emplaced irregularly in the lowermost mantle with uranium and thorium serving as a long-lived heat source. Fluid motions in the core are driven by the differential heating from above and the resulting cyclonic motions may be the source of the geodynamo.
This report summarizes the work completed in FY-2024 to analyze the power reactivity decrement (PRD) concepts of the ARC-100 core. The PRD has been traditionally defined by the reactivity change from a hot zero power (HZP) to a particular power. Consequently, the PRD accounts for the core reactivity changes due to increase coolant temperature gradient axially and radially across the core, and increased fuel temperature. The coolant temperature gradient leads to sodium and structure density changes, to radial core expansion from assembly flowering and bowing (due to axial and radial temperature gradients within the assemblies), and to control rod driveline thermal expansion. The fuel temperature increase associated with coolant temperature and power increases leads to Doppler effect and axial thermal expansion. In this work, the normal operating Hot Full Power (HFP) state is the only power level of interest, so analyses focus on the PRD calculated from HZP to HFP. The PRD has been used to assess the reactor safety features asymptotically in unprotected accident scenarios, including the loss of heat sink (LOHS), loss of flow (LOF), and transient overpower (TOP) without scram. The PRD concept relies on the “global” reactivity coefficients A, B, and C that are estimated based on “individual” reactivity effects (Doppler, sodium density, etc.). The objectives of this work are: 1) to improve and verify the methodology used to calculate the ABC coefficients used in the PRD, 2) to assess if the PRD can be used to reliably identify abnormal events. This report fulfills the FY-2024 scope of WBS#1.15.8.3 activity, “ANL0120 – Predictive Analyses of PRD as function of deformation”. The PRD concept is described in Sections 2. Additional effort in refining the methodology for core bowing modeling is performed in Sections 3. Then, two verification exercises are proposed in Section 4 to benchmark these coefficients based on direct neutronic-only calculations and on dynamic core transient simulations. Finally, the PRD approach is assessed for the detection of several unexpected events, such as primary flow perturbation or improper fuel loading, in Section 5.
This document is the compilation of the milestone portion to a larger end of project NEUP report. The executive summary of the modeling portion is provided below: In the event of cladding rupture during a postulated LOCA in a pressurized water reactor, fuel particles, along with fission gases, can be expelled into the reactor core from the fractured fuel rod, a phenomenon referred to as fuel dispersal. The initial stage of fuel dispersal is strongly influenced by the high-pressure ejection of fuel fragments, the size and geometry of the ruptured cladding, and the depressurization history of the fuel rod during the postulated LOCA transient. Depending on the location of the burst orifice relative to the quench front, the dispersal event represents an intricate three-phase flow and heat transfer phenomenon, where high-temperature fuel particles carried by the fission gases interact with the coolant within the narrow subchannels of the fuel assemblies, inducing localized phase change. Given the unique multiphysics nature of this phenomena, the current study develops a dedicated computational framework to predict the mass distribution and cooling of dispersing fuel particles, facilitating post-accident assessment and management of the fuel assemblies. Considering the scale of nuclear reactor applications, a continuum three-fluid model is proposed for simulating the transport of solids within the reactor core. With high-temperature fuel fragments within the liquid media, nucleation sites inducing phase changes are dispersed within the flow domain. Coupled with the fact that the transient dispersal event occurs on different time scales than other three-phase flow applications, this study derives a time-averaged three-fluid flow model without losing generality. The assumptions regarding the continuum treatment of the solid phase and the modeling of fuel dispersal behavior are incorporated to simplify the governing equations and derive applicable closure relations. The computational validation of the model was conducted using adiabatic experimental results obtained from ongoing research at Oregon State University, focusing on characterizing fuel dispersal behavior during simulated LOCA conditions. Settlement characteristics of the solids, quantified by the probability distribution of equivalent particles, closely matched the probability density functions reported in experimental studies. The transport of fuel particles within a scaled 5 × 5 lattice of a pressurized-water reactor rod bundle geometry was modeled through a two-fluid Eulerian framework. The required boundary conditions were evaluated from the fuel performance code BISON in a postulated large-break LOCA scenario. The modeling framework considered solid fuel particles as granular matter, interacting with the gaseous dry steam phase and fission gases through the governing interfacial momentum exchange between the participating fluids. The simulation results provided the volume fraction of the solids obtained at the bottom surface of the enclosing tank geometry. Postulated LOCA leading to fuel dispersal phenomena involves the strong coupling between fuel thermomechanics, cladding deformation, thermal-hydraulics, and fuel particle transport. Incorporation of such a strong coupling in numerical simulation is performed by coupling the multiphysics solvers. In the case of fuel dispersal, a strong coupled simulation can be performed by coupling the BISON code for fuel performance, the TRACE code for system-level thermal hydraulics, and fuel particle transport in Multiphysics Object-Oriented Simulation Environment (MOOSE). For such intricate infrastructure, the MOOSE Framework eases the data transfer between codes. The recent version of MOOSE has incorporated the Navier-Stokes module for the fluid flow. An exploratory exercise was done to gain familiarity with finite volume capabilities in the MOOSE framework to incorporate the Spalart-Allmaras (SA) turbulence model. New finite-volume and auxiliary kernels were introduced to assemble the SA transport equation, compute turbulent viscosity, and evaluate wall distance and diagnostic turbulence terms, fully integrated with existing Navier-Stokes modules. A turbulent lid-driven cavity at a Reynolds number of approximately 10,000 is used for verification. MOOSE shows the robust solver convergence and produces the turbulent features. But it underpredicts the velocity profile and turbulent quantities, emphasizing the need to develop improved SA near-wall treatments (e.g., low-Re corrections or wall functions) as a key direction for future work.
It is proposed that a gravitationally bound supermassive object is located in the 'core' of each of the Cyg A radio clouds and that these objects are responsible for the prolonged ejection of magnetic plasma which produces the synchrotron radio emission. The objects are considered to be analogous to the massive rotators, or spinars, which have been used to model the central radio and optical activity of quasars. The formation of Cyg A is described in terms of the breakup of such an object in the central optical galaxy, and the physical properties of both the compact and the giant clouds are discussed, including the effects of intergalactic drag. It is shown that the spinars are capable of effectively replenishing the extended radio clouds with relativistic electrons, even if they are embedded within the compact clouds. The spinars are modeled quantitatively, and their possible detection is considered. The requirements of the proposed model are compared with those of several other models for bright double radio galaxies.
Accurate reconstruction of the neutron flux distribution within a reactor core is essential for safe and efficient reactor operation. Traditional power shape synthesis in Light Water Reactors relies on hundreds of in-core detectors. However, this approach becomes impractical for Advanced Reactors and Microreactors due to limited space and harsh environments. To address this challenge, we propose a data-driven methodology that combines high-fidelity modeling with real-time ex-core sensor measurements, enabling the reconstruction of core power distribution while minimizing the reliance on intrusive in-core instrumentation. This project began in FY24 and achieved two initial milestones: (1) the definition of a three-year development plan for a Digital Twin framework and (2) the development of high-fidelity neutronics models of the Purdue University Reactor One (PUR-1) using both MCNP6 and OpenMC. The PUR-1 reactor, a zero-power facility, was selected due to its suitability for neutronics-focused modeling and the availability of experimental data for validation. Both models were benchmarked using neutron flux measurements obtained from irradiated gold foils, which were strategically placed within the core during a dedicated campaign in July 2024. This report marks the continuation and completion of those foundational tasks. The OpenMC model has been refined (improved geometric accuracy, expanded cross-section libraries, and refined sampling) and validated using additional experimental data. An updated sensor design—based on quadrupole configuration—was designed to measure both ex-core flux and its spatial gradient. These measurements will serve as inputs to a neural network-based reconstruction algorithm. Finally, the methodology was demonstrated on a two-dimensional test case representative of the heterogeneous material composition of the PUR-1 reactor core. A neural network implementation of the Kirchhoff-Helmholtz integral equation was employed to solve the boundary value problem using peripheral sensor measurements. The preliminary results confirm the strong potential of the proposed approach for accurate and minimally invasive neutron flux reconstruction.
The physical structure of Jupiter is discussed on the basis of data obtained by Pioneers 10 and 11. It is argued that the elemental composition of Jupiter is similar to that of the sun, and it is shown that this argument is supported by measurements of the planet's density and H/He ratio. Jupiter's shape and gravitational field are described, and a model of the planet is proposed in which there are a small iron-silicate core at the center, a very thick liquid-hydrogen stratum divided into metallic (inner) and molecular (outer) layers, and a gaseous atmosphere. According to this model, the excess heat radiated by Jupiter is simply a remnant of the heat generated when the planet coalesced from the solar nebula. The appearance of the planetary disk is described together with the Jovian magnetic field, and the Great Red Spot is shown to be a cyclonic disturbance similar to a hurricane. Effects of the Galilean satellites on the magnetic field are considered.
Models for the upper photospheres and lower chromospheres of Alpha Centauri A and B are derived from high-dispersion spectrograms of the Ca II K-line emission cores and damping wings. Effective temperatures, surface gravities, and ages consistent with the measured broadband colors, metallicities, and absolute magnitudes are estimated for the two stars. The spectrograms are calibrated by fitting the far-wing K profiles with synthetic fluxes based on radiative-equilibrium models, and the model atmospheres are obtained through a partial-redistribution analysis of the line cores and damping wings. These models are compared with those previously proposed for Procyon, Arcturus, and the sun. Many features in all the models are found to be quite similar, and some evidence suggests that Alpha Cen A may be significantly older and more evolved than the sun, even though both stars exhibit similar chromospheric properties.
This paper presents a multiscale modeling framework (MMF) to model moist atmospheric limited-area weather. The MMF resolves large-scale convection using a coarse grid while simultaneously resolving local features through numerous fine local grids and coupling them seamlessly. Both large- and small-scale processes are modeled using the compressible Navier-Stokes equations within the Nonhydrostatic Unified Model of the Atmosphere (NUMA), and are discretized using a continuous element-based Galerkin method (spectral elements) with high-order basis functions. Consequently, the large-scale and small-scale models share the same dynamical core but have the flexibility to be adjusted individually. The proposed MMF method is tested in 2D and 3D idealized limited-area weather problems involving storm clouds produced by squall line and supercell simulations. Numerical results from the MMF showed enhanced representation of cloud processes compared to the coarse model.
Models of the internal structure of Mars consistent with the mass, radius and moment of inertia of the planet are constructed. The models assume that the radius of the core is between 0.36 and 0.60 of the radius of the planet, that the zero-pressure density of the mantle is between 3.54 and 3.49 g/cu cm, and that the planet contains 25 to 28% iron. Meteorite models of Mars containing 25 wt % iron and 12 wt % core are also proposed. It is maintained that Mars in contrast to the earth is an incompletely differentiated planet with a core substantially richer in sulfur than the core of the earth. The absence of a magnetic field on Mars is possibly linked with lack of lunar precessional torque and the small size and high resistivity of the Martian core.