Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “equilibrium core”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Impact of U-10Mo HALEU fuel element tolerances on the Massachusetts Institute of Technology reactor safety and operational performance – Thermal hydraulics

The U.S. is coordinating efforts for the conversion of six U.S. High Performance Research Reactors (USHPRR) including one critical facility from highly enriched uranium (HEU) to low-enriched uranium (LEU). In order to continue the mission of these reactors, including the Massachusetts of Institute of Technology Reactor (MITR), and achieve similar performance, high-assay low-enriched uranium (HALEU) with a high-density metallic alloy of uranium with 10 wt% molybdenum (U-10Mo) is being evaluated. The impact of the fabrication specification and tolerances was assessed following the preliminary design of the MITR LEU fuel elements using the U-10Mo monolithic alloy. This research focuses on the analysis of fabrication specification impact on thermal hydraulics (TH) characteristic of the MITR LEU core as a function of the variation of the relevant fuel specification parameters (e.g., coolant channel gap thickness, fuel plate thickness, etc.). The analyses are performed based on an all-fresh LEU fuel conversion plan identified in a preliminary safety analysis report submitted to the Nuclear Regulatory Commission. The reactor power margin to the onset of nucleate boiling (ONB) is assessed under the limiting safety system settings (LSSS), where a scram occurs, to ensure there is sufficient margin to the reactor safety limit, which is defined by the onset of flow instability that occurs after the ONB. The best estimate plus uncertainty approach is employed to analyze this TH characteristic, which yields realistic results while maintaining adequate conservatism, utilizing a statistical uncertainty propagation method with the STAT7 code. The TH characteristic is analyzed as a function of the variability of the specification parameters resulting from the fabrication process. The main findings of this study show that the MITR core can meet the TH safety and operational requirements at the all-LEU initial core startup (cycle 1), selected transition cycles (most reactive cycle and most limiting cycle: cycle 3 and 5, respectively) and equilibrium (cycle 14) cores under all limiting fabrication parameter combinations considered. In addition, the analyses show that the dependency of the core power margin to ONB on those specification parameters that have the most direct impact on TH performance is non-linear but monotonically decreasing within the specification tolerances. The third order polynomial fit curves are reported in detail for selected limiting cases and can serve as a powerful tool for future MITR fuel management in cases such as when HALEU supply is established that may allow additional cycle length or other operational benefits.

Conversion↗
NASA NTRS1981

Models of the earth's core

Combined inferences from seismology, high-pressure experiment and theory, geomagnetism, fluid dynamics, and current views of terrestrial planetary evolution lead to models of the earth's core with five basic properties. These are that core formation was contemporaneous with earth accretion; the core is not in chemical equilibrium with the mantle; the outer core is a fluid iron alloy containing significant quantities of lighter elements and is probably almost adiabatic and compositionally uniform; the more iron-rich inner solid core is a consequence of partial freezing of the outer core, and the energy release from this process sustains the earth's magnetic field; and the thermodynamic properties of the core are well constrained by the application of liquid-state theory to seismic and labroatory data.

Stevenson, D. J.↗
NASA NTRS2015

The Combined Strength of Thermodynamics and Comparative Planetology: Application of Activity Models to Core Formation in Terrestrial Bodies

Recent models for accretion of terrestrial bodies involve metal-silicate equilibrium as the metallic core formed during growth. Most elements considered are either refractory or well studied elements for which effects of pressure, temperature, oxygen fugacity, and metallic liquid composition are well known. There are a large number of elements that are both siderophile and volatile, whose fate in such models is unknown, largely due to a lack of data at comparable conditions and com-positions (FeNi core with light elements such as S, C, Si, and O). We have focused on Ge, In, As, Sb and determined the effect of Si and C on metal-silicate partitioning, and developed a thermo-dynamic model that allows application of these new data to a wide range of planetary bodies. New experiments: We have previously carried out experiments with FeSi metallic liquid at C-saturated conditions at 1600 and 1800 C [4]. In a new series of experiments we investigate the effect of Si in carbon-free systems at 1600 C for comparison. Experiments were carried out at 1 GPa in MgO capsules using the same basaltic starting composition as in previous studies. The MgO capsule reacts with the silicate melt to form more MgO-rich liquids that have 22-26 wt% MgO. Experimental met-als and silicates were analyzed using a combination of electron microprobe analysis and laser ablation ICP-MS. Results: The new results can be interpreted by considering Ge as an example, in the simple exchange equilibrium Fe + GeO = FeO + Ge, where the equilibrium constant Kd can be examined as a function of Si content of the metal. The slope of lnKd vs. (1-XSi) for this new series allows derivation of the epsilon interaction parameter for each of these four elements and Si (both C-saturated and C-free).All four elements have positive epsilon values, indicating that Si causes a decrease in the partition coefficients; values are 6.6, 6.5, 27.8 and 25.2 for In, Ge, As, and Sb, respectively, at 1 GPa and 1600 C. As an example of how large the effect of Si can be, these epsilon values correspond to activity coefficients (gamma) for As of 0.01 when XSi = 0, and up to gamma = 23 when XSi = 0.2. Combining these new results with previous determinations [5,6] of epsilon parameters for S and C for these elements allows us calculate activity of Ge, In, As, and Sb in Fe-Ni-Si-S-C-O metallic liquids. We apply this new model to sever-al terrestrial bodies such as Earth (Si-rich core), Mars (S-rich core), Moon (S-, C-, and Si-poor core), and Vesta, and examine the resulting core and mantle concentrations of these elements. Mantle concentrations of these four elements are well explained for Earth and Mars in models that call for mid-mantle equilibration between Si-bearing and S-bearing FeNi cores, respectively. Modeling results for the Moon and Vesta will also be presented.

Righter, K.↗
NASA NTRSDate not supplied

Machine Learning Applications to Metal-Silicate Equilibria and their Insights into Core Formation

An extensive number of studies have experimentally investigated how elements distribute between metal and silicate phases, to better constrain core-mantle chemical equilibrium. Here, we present a new database compiling all (to our knowledge) experimental data on liquid metal-silicate partitioning from 118 peer-reviewed publications. We applied various machine learning techniques to gain further insights into these partitioning equilibria and their dependencies. We performed a network analysis to investigate the relationship between experiments and partition coefficients, which enables visualizing gaps in the experimental dataset and biases related to varying experimental conditions and analytical setup. In addition, semi-empirical thermodynamic models are commonly used to extrapolate these chemical reactions to the wide range of pressure, temperature and compositional conditions of planetary differentiation. These models are based on linear regressions that assume continuous relationship between partition coefficients and experimental variables. Here, we considered random forest regressions, which are algorithms based on ensembles of decision trees and does not consider continuous effects of each variable. The application of this regression significantly improves the prediction of metal-silicate partitioning for several elements including Ni, Si and Cr. We will show how this new approach improves our understanding of elemental exchange between metal and silicate and their implications for the Earth’s core formation.

siderophile element↗
NASA NTRS1974

The internal structure of the convective envelope of the sun

An axisymmetric hydrodynamic model of the sun consisting of a core at radiant equilibrium and a convective envelope is studied with consideration of gravitational forces, rotation, turbulence effects, and electromagnetic forces. The dimensionless parameters that control the equations of plasma motion within the sun and the boundary conditions are investigated; estimates are given for these parameters on the basis of experimental data and the existing theoretical models. Stability conditions are derived for convective equilibrium, developed cellular convection, and developed turbulent convection with consideration of rotation and turbulent Reynolds stresses. These conditions characterize the possible physical state of the matter in the convective envelope. Using estimates of the dimensionless parameters and models of some of them, a method of successive approximations is proposed for solution of the problem as a whole.

Yavorskaya, I. M.↗
NASA NTRS1974

N-Player Stochastic Differential Games

Conditions are described which guarantee that the control strategies adopted by N players constitute an efficient solution, an equilibrium, or a core solution. The system dynamics are described by an Ito equation, and all players have perfect information. It was found that when the set of instantaneous joint costs and velocity vectors is convex, the conditions are necessary.

Varaiya, P.↗
NASA NTRS1976

N-player stochastic differential games

The paper presents conditions which guarantee that the control strategies adopted by N players constitute an efficient solution, an equilibrium, or a core solution. The system dynamics are described by an Ito equation, and all players have perfect information. When the set of instantaneous joint costs and velocity vectors is convex, the conditions are necessary.

Varaiya, P.↗
NASA NTRS1985

Exercise thermoregulation in men after 6 hours of immersion

The present investigation is concerned with thermoregulation at rest and during exercise after water-immersion deconditioning, giving particular attention to the effects of fluid shifts and negative water balance on sweat rate and rectal temperature. Six healthy males 20-35 years old were used in the experiments. Rectal and mean skin temperature, skin heat conductance, heart rate, and total body sweat rate were measured during 70 min of supine leg exercise at 50 percent of peak O2 uptake. The data were taken after a 6-h control period in air and after immersion to the neck in water (34.5 C) for 6 h after overnight food and fluid restriction. Attention is given to end exercise heart rates and data during exercise. The obtained results suggest that, compared with control responses, the equilibrium level of core temperature during submaximal exercise is regulated at a higher level after immersion.

Greenleaf, J. E.↗
NASA NTRS1969

ELAS: A general-purpose computer program for the equilibrium problems of linear structures. Volume 2: Documentation of the program

A general purpose digital computer program for the in-core solution of linear equilibrium problems of structural mechanics is documented. The program requires minimum input for the description of the problem. The solution is obtained by means of the displacement method and the finite element technique. Almost any geometry and structure may be handled because of the availability of linear, triangular, quadrilateral, tetrahedral, hexahedral, conical, triangular torus, and quadrilateral torus elements. The assumption of piecewise linear deflection distribution insures monotonic convergence of the deflections from the stiffer side with decreasing mesh size. The stresses are provided by the best-fit strain tensors in the least squares at the mesh points where the deflections are given. The selection of local coordinate systems whenever necessary is automatic. The core memory is used by means of dynamic memory allocation, an optional mesh-point relabelling scheme and imposition of the boundary conditions during the assembly time.

Utku, S.↗
DOE OSTI2026

Reduced fast-ion transport calculations of m = n = 1 fishbone-like instabilities in MAST-U

Fast-ion transport associated with an m = n = 1 fishbone-like burst in MAST-U discharge 47128 is investigated using a reduced guiding-center-based transport model (ORBIT-Kick) constrained by multi-diagnostic measurements. The two-dimensional beam-emission spectroscopy system provides measurements of the core poloidal mode structure and fluctuation amplitude, while EFIT++ reconstructions constrained by the motional Stark effect diagnostic indicate a flat q-profile with q 0 > 1⁠, indicating the absence of a resonant q = 1 surface and supporting a pressure-driven infernal-mode interpretation. Analytic m = n = 1 displacement profiles consistent with the measured core mode structure and equilibrium constraints are used as the mode structure inputs to ORBIT-Kick. The calculations show that the dominant resonances occur between the mode and co-passing fast ions, producing redistribution localized near the magnetic axis. Synthetic neutron camera signals from TRANSP-Kick recover up to 90% of the experimentally observed neutron deficit at the time of peak mode amplitude, indicating that the measured m = n = 1 mode is a dominant contributor to core fast-ion transport. However, the synthetic neutron signals recover rapidly, whereas the measured neutron emission continues to decrease after the peak amplitude. In conclusion, the remaining discrepancy may arise from contributions not included in the present single-harmonic model, including higher-m and higher-n harmonics, multi-harmonic interactions, and additional transport mechanisms, motivating future diagnostic development and modeling efforts to resolve and incorporate these additional contributions.

Wong, Henry H. [University of California, Los Ange↗
NASA NTRS2004

On the Final Mass of Giant Planets

In the core accretion model of giant planet formation, when the core reaches critical mass, hydrostatic equilibrium is no longer possible and gas accretion ensues. If the envelope is radiative, the critical core mass is nearly independent of the boundary conditions and is roughly M(sub crit) ~ 10Mass of the Earth (with weak dependence on the rate of planetesimal accretion M(sub core) and the disk opacity k). Given that such a core may form at the present location of Jupiter in a time comparable to its Type I migration time (10(exp 5) - 10(exp 6) years) provided that the nebula was significantly enhanced in solids with respect to the MMSN and stall at this location in a weakly turbulent (alpha approximately less than 10(exp -4) disk, it may be appropriate to assume that such objects inevitably form and drive the evolution of late-phase T Tauri star disks. Here we investigate the final masses of giant planets in disks with one or more than one such cores. Although the presence of several planets would lead to Type II migration (due to the effective viscosity resulting from the planetary tidal torques), we ignore this complication for now and simply assume that each core has stalled at its location in the disk. Once a core has achieved critical mass, its gaseous accretion is governed by the given Kelvin-Helmholtz timescale.

Estrada, P. R.↗
DOE OSTI2025

High-Temperature Gas-Cooled Reactors Multiphysics Simulation Demonstration and Code Validation

This study presents a comprehensive benchmarking and verification effort of several thermal-hydraulic and multiphysics capabilities for high-temperature gas-cooled reactor applications. The first part of this effort focuses on the running-in verification of Griffin’s multiphysics capabilities, specifically for simulating the evolution of pebble-bed reactor cores from startup to equilibrium. Since Fiscal Year 2024, improvements and enhancements have been implemented in Griffin, including simplifying the process to specify streamlines and developing the online cross-section generation capability. In the absence of validation data, code-to-code comparisons are conducted with kugelpy, showing good agreement for integral quantities like k-eff predictions and predictions for maximum power density. However, accuracy issues are noted for more detailed quantities like the spatial distribution of fission rate densities which will require further work to address. The second part of this report presents an improved System Analysis Module (SAM) core channel model where the effects of cross flow are considered during the pressurized loss of forced cooling transient, resulting in an improved agreement of the predicted pebble temperature with respect to the predictions from the SAM 2D porous media model. Additionally, the wall channeling effect due to variable porosity at the near wall region of the core is also investigated. Furthermore, to demonstrate Griffin’s online cross-section generation capability, a Multiphysics simulation is performed by coupling Griffin to the SAM core channel model. In the third part of the report, as a part of the Organisation for Economic Co-operation and Development/Nuclear Energy Agency (OECD/NEA) thermal-hydraulic code validation benchmark activity for a high-temperature gas-cooled reactor, the High Temperature Test Facility (HTTF) is investigated first using the NekRS computational fluid dynamics (CFD) code to study the flow mixing phenomenon in the lower plenum of the facility. Then, code-to-code and code-to-data comparisons are performed for Test PG27, which is a pressurized conduction cooldown (PCC) test, using five different codes by six organizations from five countries. The different simulations show good agreements in terms of the general trend but there are differences in some results such as the peak temperatures of different regions and heat removal rate.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1992

A new class of g-modes in neutron stars

Because a neutron star is born hot, its internal composition is close to chemical equilibrium. In the fluid core, this implies that the ratio of the number densities of charged particles (protons and electrons) to neutrons is an increasing function of the mass density. This composition gradient stably stratifies the matter giving rise to a Brunt-Vaisala frequency N of about 500/s. Consequently, a neutron star core provides a cavity that supports gravity modes (g-modes). These g-modes are distinct from those previously identified with the thermal stratification of the surface layers and the chemical stratification of the crust. We compute the lowest-order, quadrupolar, g-modes for cold, Newtonian, neutron star models with M/solar M = 0.581 and M/solar M = 1.405, and show that the crustal and core g-modes have similar periods. We also discuss damping mechanisms and estimate damping rates for the core g-modes. Particular attention is paid to damping due to the emission of gravitational radiation.

Reisenegger, Andreas↗
NASA NTRS2012

Provenance and Concentration of Water in the Shergottite Mantle

The water content of the martian mantle is controversial. In particular, the role of water in the petrogenesis of the shergottites has been much debated. Although the shergottites, collectively, contain very little water [e.g., 1,2], some experiments have been interpreted to show that percent levels of water are required for the petrogenesis of shergottites such as Shergotty and Zagami [3]. In this latter interpretation, the general paucity of water in the shergottites and their constituent minerals is attributed to late-stage degassing. Y980459 (Y98) is a very primitive, perhaps even parental, martian basalt, with a one-bar liquidus temperature of approx.1400 C. Olivine is the liquidus phase, and olivine core compositions are in equilibrium with the bulk rock [e.g., 4]. Petrogenetically, therefore, Y98 has had a rather simple history and can potentially help constrain the role of water in martian igneous processes. In particular, once trapped, melt inclusions should not be affected by subsequent degassing.

Jones, J. H.↗
NASA NTRS1986

Fe-Mg-Mn relations of ureilite olivines and pyroxenes and the genesis of ureilites

Microprobe analyses of ureilite and pigeonite cores are studied. The Fe/Mn-Fe/Mg relationship in the olivine core is examined. It is observed that magnetic processs such as fractional crystallization and partial melting, and FeO reduction contribute to the olivine core composition. The study of the Mg/Mn and Fe/Mn distributions reveals that these two distributions are not in equilibrium in the olivine and pigeonite cores. The effect of a reducing agent, carbon, on the ureilite genesis is investigated. It is concluded that fractional crystallization and FeO reduction are the major processes of ureilite genesis.

Mittlefehldt, D. W.↗
NASA NTRS1987

Experiments pertaining to the formation and equilibration of planetary cores

The phase diagram of FeO was experimentally determined to pressures of 155 GPa and temperatures of 4000 K using shock wave and diamond-cell techniques. Researchers discovered a metallic phase of FeO at pressures greater than 70 GPa and temperatures exceeding 1000 K. The metallization of FeO at high pressures implies that oxygen can be present as the light alloying element of the Earth's outer core, in accord with the geochemical predictions of Ringwood. The high pressures necessry for this metallization suggest that the core has acquired its composition well after the initial stages of the Earth's accretion. The core forming alloy can react chemically with oxides such as those forming the mantle. The core and mantle may never have reached complete chemical equilibrium, however. If this is the case, the core-mantle boundary is likely to be a zone of active chemical reactions.

Jeanloz, Raymond↗
NASA NTRS1986

High-pressure metallization of FeO and implications for the earth's core

The phase diagram of FeO has been experimentally determined to pressures of 155 GPa and temperatures of 4000 K using shock-wave and diamond-cell techniques. A metallic phase of FeO is observed at pressures greater than 70 GPa and temperatures exceeding 1000 K. The metallization of FeO at high pressures implies that oxygen can be present as the light alloying element of the earth's outer core, in accord with the geochemical predictions of Ringwood (1977 and 1979). The high pressures necessary for this metallization suggest that the core has acquired its composition well after the initial stages of the earth's accretion. Direct experimental observations at elevated pressures and temperatures indicate that core-forming alloy can react chemically with oxides such as those forming the mantle. The core and mantle may never have reached complete chemical equilibrium, however. If this is the case, the core-mantle boundary is likely to be a zone of active chemical reactions.

Knittle, Elise↗
NASA NTRS1972

Evolution of stars with suppressed core convection

Stellar evolution on the upper main sequence was computed for models of stars with cores assumed to be in radiative equilibrium, up to the point of central helium ignition. The role of the Schonberg-Chandrasekhar limit for an isothermal core is found to be critical for the evolutionary tracks. Observational data are used to rule out the hypothesis of evolution with radiative cores (in upper main-sequence stars) and, by implication, of magnetic fields that are sufficiently strong to have suppressed the core convention.

Stothers, R. B.↗