Engineering PapersSearch

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 91 records · Page 5

The Effect of Metal Composition on Fe-Ni Partition Behavior between Olivine and FeNi-Metal, FeNi-Carbide, FeNi-Sulfide at Elevated Pressure

Metal-olivine Fe-Ni exchange distribution coefficients were determined at 1500 C over the pressure range of 1 to 9 GPa for solid and liquid alloy compositions. The metal alloy composition was varied with respect to the Fe/Ni ratio and the amount of dissolved carbon and sulfur. The Fe/Ni ratio of the metal phase exercises an important control on the abundance of Ni in the olivine. The Ni abundance in the olivine decreases as the Fe/Ni ratio of the coexisting metal increases. The presence of carbon (up to approx. 3.5 wt.%) and sulfur (up to approx. 7.5 wt.%) in solution in the liquid Fe-Ni-metal phase has a minor effect on the partitioning of Fe and Ni between metal and olivine phases. No pressure dependence of the Fe-Ni-metal-olivine exchange behavior in carbon- and sulfur-free and carbon- and sulfur-containing systems was found within the investigated pressure range. To match the Ni abundance in terrestrial mantle olivine, assuming an equilibrium metal-olivine distribution, a sub-chondritic Fe/Ni-metal ratio that is a factor of 17 to 27 lower than the Fe/Ni ratios in estimated Earth core compositions would be required, implying higher Fe concentrations in the core forming metal phase. A simple metal-olivine equilibrium distribution does not seem to be feasible to explain the Ni abundances in the Earth's mantle. An equilibrium between metal and olivine does not exercise a control on the problem of Ni overabundance in the Earth's mantle. The experimental results do not contradict the presence of a magma ocean at the time of terrestrial core formation, if olivine was present in only minor amounts at the time of metal segregation.

Holzheid, Astrid

Examining transport and integrated modeling predictive capabilities for negative-triangularity scenarios

This paper investigates the predictive capabilities of TGYRO and TGLF models in assessing the performance of negative triangularity (NT) plasmas compared to positive triangularity (PT) plasmas in fusion devices. TGYRO predicts kinetic profiles, while TGLF analyzes turbulent transport. The study reveals that TGYRO reasonably predicts NT profiles similar to PT, although it overpredicts the high-power scenarios where there is increased experimental MHD activity. TGLF analysis finds reduced linear growth rates in NT and altered flux spectra relative to PT. Additionally, the TGLF SAT0 saturation model is observed to predict high-k transport and a reduction of particle transport with the electron temperature gradient. These findings are further corroborated by core-pedestal modeling using the Stability Transport Equilibrium Pedestal workflow, showing stronger confinement improvements in NT, particularly at higher power densities for the SAT0 saturation model. Furthermore, the study underscores the importance of accurately capturing turbulence saturation mechanisms for NT in order to project its performance accurately in fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Conceptual study on using Doppler backscattering to measure magnetic pitch angle in tokamak plasmas

We introduce a new approach to measure the magnetic pitch angle profile in tokamak plasmas with Doppler backscattering (DBS), a technique traditionally used for measuring flows and density fluctuations. The DBS signal is maximised when its probe beam's wavevector is perpendicular to the magnetic field at the cutoff location, independent of the density fluctuations [Hillesheim \emph{et al} 2022 \emph{Nucl. Fusion} \textbf{55} 073024]. Hence, if one could isolate this effect, DBS would then yield information about the magnetic pitch angle. By varying the toroidal launch angle, the DBS beam reaches cutoff with different angles with respect to the magnetic field, but with other properties remaining similar. Hence, the toroidal launch angle which gives maximum backscattered power is thus that which is matched to the pitch angle at the cutoff location, enabling inference of the magnetic pitch angle. We performed systematic scans of the DBS toroidal launch angle for repeated DIII-D tokamak discharges. Experimental DBS data from this scan were analysed and combined with Gaussian beam-tracing simulations using the Scotty code [Hall-Chen \emph{et al} 2022 \emph{Plasma Phys. Control. Fusion} \textbf{64} 095002]. The pitch-angle inferred from DBS is consistent with that from magnetics-only and motional-Stark-effect-constrained (MSE) equilibrium reconstruction in the edge. In the core, the pitch angles from DBS and magnetics-only reconstructions differ by one to two degrees, while simultaneous MSE measurements were not available. The uncertainty in these measurements was under a degree; we show that this uncertainty is primarily due to the error in toroidal steering, the number of toroidally separated measurements, and shot-to-shot repeatability. We find that the error of pitch-angle measurements can be reduced by optimising the poloidal launch angle and initial beam properties. Since DBS has high spatial and temporal resolutions, is non-perturbative, does not require neutral beams, and is likely robust to neutron damage of and debris on the first mirrors, using DBS to measure the pitch angle in future fusion energy systems is especially appealing.

beam tracing

Strong Coupling of Hydrodynamics and Reactions in Nuclear Statistical Equilibrium for Modeling Convection in Massive Stars

We build on the simplified spectral deferred corrections (SDC) coupling of hydrodynamics and reactions to handle the case of nuclear statistical equilibrium (NSE) and electron/positron captures/decays in the cores of massive stars. Our approach blends a traditional reaction network on the grid with a tabulated NSE state from a very large, ${\mathcal O }(100)$ nuclei network. We demonstrate how to achieve second-order accuracy in the simplified-SDC framework when coupling NSE to hydrodynamics, with the ability to evolve the star on the hydrodynamics time step. We discuss the application of this method to convection in massive stars leading up to core collapse. We also show how to initialize the initial convective state from a 1D model in a self-consistent fashion. All of these developments are done in the publicly available Castro simulation code and the entire simulation methodology is fully GPU-accelerated.

Explosive nucleosynthesis

Framework for phase transitions between the Maxwell and Gibbs constructions at finite temperature

The characteristics of the hadron-to-quark first-order phase transition differ depending on whether charge neutrality is locally or globally fulfilled. In 𝛽-equilibrated matter, these two possibilities correspond to the Maxwell and Gibbs constructions. Recently, we presented a new framework in which a continuously varying parameter allows one to describe a first-order phase transition in intermediate scenarios to the two extremes of fully local and fully global charge neutrality. In this work, we extend the previous framework to finite temperatures and out-of-𝛽 equilibrium conditions, making it available for simulations of core-collapse supernovae and binary neutron star mergers. We investigate its impact on key thermodynamic quantities across a range of baryon densities, temperatures, and electron fractions. We find that when matter is not in 𝛽 equilibrium, the pressure in the mixed phase is not constant even for the case of fully local charge neutrality. Moreover, we compute the thermal index using three different approaches, demonstrating that the finite-temperature extension of an equation of state using a constant thermal index can be ill defined when applied to the mixed phase.

QCD phase transitions

Kinetic and Thermodynamic Characterization of the Cobalt and Manganese Oxyhydroxide Cores Formed in Horse Spleen Ferritin

Horse spleen ferritin (HoSF) containing 800-1500 cobalt or 250-1200 manganese atoms as Co(O)OH and Mn(O)OH mineral cores within the HoSF interior (Co-HoSF and Mn-HoSF) was synthesized, and the chemical reactivity, kinetics of reduction, and the reduction potentials were measured. Microcoulometric and chemical reduction of HoSF containing the M(O)OH mineral core (M = Co or Mn) was rapid and quantitative with a reduction stoichiometry of 1.05+/-0.10 e/M forming a stable M(OH)2 mineral core. At pH 9.0, ascorbic acid (AH2), a two-electron reductant, effectively reduced the mineral cores; however, the reaction was incomplete and rapidly reached equilibrium. The addition of excess AH2 shifted the reaction to completion with a M(3+)/AH2 stoichiometry of 1.9-2.1, consistent with a single electron per metal atom reduction. The rate of reaction between M(0)OH and excess AH2 was measured by monitoring the decrease in mineral core absorbance with time. The reaction was first order in each reactant with second-order rate constants of 0.53 and 4.74/M/min, respectively, for Co- and Mn-HoSF at pH 9.0. From the variation of absorbance with increasing AH2 concentration, equilibrium constants at pH 9.0 of 5.0+/-1.9 for Co-HoSF and 2.9+/-0.9 for Mn-HoSF were calculated for 2M(O)OH + AH2 = 2M(OH)2 f D, where AH2 and D are ascorbic acid and dehydroascorbic acid, respectively. Consistent with these equilibrium constants, the standard potential for the reduction of Co(III)-HoSF is 42 mV more positive than that of the ascorbic acid reaction, while the standard potential of Mn(III)-HoSF is 27 mV positive relative to AH2. Fe(2+) in solution with Co- and Mn-HoSF under anaerobic conditions was oxidized to form Fe(O)OH within the HoSF interior, resulting in partial displacement of the Co or Mn by iron.

Zhang, Bo

Partnership Center for High-fidelity Boundary Plasma Simulation

The University of Colorado Boulder contributed to the DOE Partnership Center for High‑fidelity Boundary Plasma Simulation (award DE‑SC0018271) by advancing gyrokinetic modeling of core–edge plasma interactions in tokamaks. Using the GEM and XGC codes, the project investigated turbulence spreading from the plasma edge to the core, finding that while edge‑driven trapped electron mode turbulence enhances core ion heat flux in DIII‑D L‑mode plasmas, it does not fully explain experimentally observed transport shortfalls, indicating a need for flux‑driven models. The team also demonstrated that edge micro‑instabilities are highly sensitive to small variations in magnetic equilibrium, particularly near the pedestal and edge, highlighting the importance of accurately characterizing equilibrium uncertainties. Significant progress was made in spatially coupling the core GEM and edge XGC simulations, with successful validation against standalone XGC results in realistic tokamak geometries, laying the groundwork for future extensions to kinetic electrons and electromagnetic effects.

42 ENGINEERING

Comparison of dynamic crystallization techniques on Apollo 15 quartz normative basalts

The paper compares dynamic crystallization techniques used on Apollo 15 quartz normative basalts which show differences related to their method. The variation of pyroxene core major element chemistry with cooling rate is confirmed, and it is shown that the Fe loss from liquid is shifted toward the En-Di joint; it is suggested that the nucleation temperature and breakdown of equilibrium partitioning as a function of cooling rate dictate the core composition. It is concluded that the different crystallization behavior that results from the Fe-capsule and Fe-loop container methods is related to the nucleation kinetics.

Lofgren, G. E.

Dopant, coating, and grating effects in silica optical fibers under extreme neutron irradiation

Radiation effects on fiber Bragg gratings (FBGs) have been studied but data gaps related to the effects of displacement damage resulting from high fast neutron fluence remain. In this work, Type-I and Type-II FBGs inscribed in optical fibers with various core dopants (Ge and F) and fiber coatings (acrylate, polyimide) were monitored in situ during 75 days of neutron irradiation to a peak fast (>0.1 MeV) neutron fluence of 3 X 10 21 n fast /cm 2 . Further, the reflected intensity of the Type-I FBGs inscribed in a Ge-doped core fiber decreased by > 30 dB within 7 h of irradiation (10 19 n fast /cm 2 ), whereas Type-II FBGs inscribed in a pure silica core fiber eventually approached >40 dB attenuation after accumulating a fast neutron fluence on the order of 10 20 n fast /cm 2 . Type-II FBGs inscribed in F-doped core fiber improved stability: the attenuation approached an equilibrium value in the range of 10 to 20 dB.

47 OTHER INSTRUMENTATION

FLARE: field line analysis and reconstruction for 3D boundary plasma modeling

The FLARE code is a magnetic mesh generator that is integrated within a suite of tools for the analysis of the magnetic geometry in toroidal fusion devices. A magnetic mesh is constructed from field line segments and permits fast reconstruction of field lines in 3D boundary plasma codes such as EMC3-EIRENE. Both intrinsically non-axisymmetric configurations (stellarators) and those with symmetry breaking perturbations of an axisymmetric equilibrium (tokamaks) are supported. The code itself is written in Modern Fortran with MPI support for parallel computing, and it incorporates object-oriented programming for the definition of the magnetic field and the material surface geometry. Extended derived types for a number of different magnetohydrodynamic equilibrium and plasma response models are implemented. The core element of FLARE is a field line tracer with adaptive step-size control, and this is integrated into tools for the construction of Poincaré maps and invariant manifolds of X-points. A collection of high-level procedures that generate output files for visualization is build on top of that. The analysis modules are build with Python frontends that facilitate customization of tasks and/or scripting of parameter scans.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Plant Reload Optimization (prlo)

The PRLO framework is built on a modular and extensible architecture that tightly couples advanced evolutionary optimization algorithms with nuclear fuel depletion solvers (i.e., nuclear physics neutronics code). It supports exploring complex, high-dimensional design spaces constrained by user-specified operational, safety, and economic constraints. Objectives such as minimizing fresh fuel enrichment, flattening radial and axial power distributions, and maximizing discharge burnup are evaluated. PRLO’s equilibrium cycle optimization capability enables the identification of core configurations that maintain fuel cycle sustainability over extended planning horizons. Its integration with the RAVEN platform facilitates optimization of loading patterns or fuel shuffling schemes across multiple cycles. The interface with SIMULATE, a licensed industry-standard nodal code developed by Studsvik, ensures accurate neutronic and thermal-hydraulic feedback for reactor core design. PRLO’s automated workflow engine supports iterative design refinement, enabling utilities to streamline core design processes and meet evolving performance and regulatory targets.

Kim, Junyung [Idaho National Laboratory] (00090005

Broadband photoexcitation of lasers

Existing techniques are reviewed and novel approaches for broadband photoexcitation of lasers are discussed. This subject is of interest because of the emergence of future intense near equilibrium and blackbody sources, such as the gas core reactor, and electric discharge plasma devices, which are the key to the development of intense broadband radiation sources. The possible use of radiation from solar concentrators in space also fits into this general area. The varied use of photoexcitation of lasers are discussed ranging over direct optical pumping, photodissociation, and photoreionization. In addition, the new uses of laser transitions and chemical reactions is shown. This is accomplished through multiline operation of high pressure tunable lasers which can use the broadband radiation for multiphoton excitation.

Hess, R. V.

Thermoregulatory responses of rats to varying environmental temperatures

The peripheral (tail) circulatory responses of six male albino rats were measured at ambient temperatures between 5 and 40 C, using impedance plethysmography. Each animal was anesthetized, instrumented, and placed in a thermal environmental chamber to reach equilibrium. Tail blood flow, respiration rate, heart rate, core temperature, and tail skin temperature were then monitored at each ambient temperature. The mean tail blood flow was significantly higher (p less than 0.05) at 5 C than at 10 C. The mean blood flow increased significantly (p less than 0.01) at each of the temperatures above 10 C. Tail skin temperature and internal body (core) temperature increased significantly with increasing ambient temperature.

Berry, J. J.

Geologic History of Asteroid 4 Vesta

Some types of meteorites - most irons, stony irons, some achondrites - hail from asteroids that were heated to the point where magmatism occurred within a very few million years of the formation of the earliest solids in the solar system. The largest clan of achondrites, the howardite, eucrite and diogenite (HED) meteorites, represent the crust of their parent asteroid]. Diogenites are cumulate harzburgites and orthopyroxenites from the lower crust whilst eucrites are basalts, diabases and cumulate gabbros from the upper crust. Howardites are impact-engendered breccias mostly of diogenites and eucrites. There remains only one large asteroid with a basaltic crust, 4 Vesta, which is thought to be the source of the HED clan. Differentiation models for Vesta are based on HED compositions. Proto-Vesta consisted of chondritic materials containing Al-26, a potent, short-lived heat source. Inferences from compositional data are that Vesta was melted to high degree (≥50%) allowing homogenization of the silicate phase and separation of a metallic core. Convection of the silicate magma ocean allowed equilibrium crystallization, forming a harzburgitic mantle. After convective lockup occurred, melt collected between the mantle and the cool thermal boundary layer and underwent fractional crystallization forming an orthopyroxene-rich (diogenite) lower crust. The initial thermal boundary layer of chondritic material was replaced by a mafic upper crust through impact disruption and foundering. The mafic crust thickened over time as additional residual magma intrudes and penetrates the mafic crust forming plutons, dikes, sills and flows of cumulate and basaltic eucrite composition. This magmatic history may have taken only 2-3 Myr. This magma ocean scenario is at odds with a model of heat and magma transport that indicates that small degrees of melt would be rapidly expelled from source regions, precluding development of a magma ocean. Constraints from radiogenic Mg-26 distibutions suggest that the parent asteroid of HEDs was much smaller than Vesta. Thus, first-order questions regarding asteroid differentiation remain.

Mittlefehldt, David W.

Core saturation in a moving medium

A numerical technique for solving the line transfer equation of a two-level atom in static equilibrium is presented. Complete redistribution of emitted photons is assumed, as is saturation at the line core. Emission intensity is calculated either by a generalized Eddington-Barber relation, a first-order differential equation for the specific intensity, or by a formal transfer integral. Sample calculations are performed of the line transfer equation in a semi-infinite atmosphere with a constant Planck function of the collision parameter and for the Mg II resonance line in a model solar atmosphere experiencing shocks. Attention is focused on the line wings in the latter problem. The first order differential equation approach yields the best intensity values and temperature structure.

Kalkofen, W.

Effect of the fluid core on changes in the length of day due to long period tides

The long period luni-solar tidal potential is known to cause periodic changes in the earth's rotation rate. It is found that the effect of a dissipationless fluid outer core is to reduce the amplitudes of these tidal perturbations by about 11 percent. When the fluid core effect is added to Agnew and Farrell's (1978) estimate of the effect of an equilibrium ocean, the result is in accord with observation. The effects of dissipative processes within the fluid core are also examined. Out-of-phase perturbations are found which could be as large as about 10 ms at 18.6 yr. It is concluded, however, that the poorly understood decade fluctuations in the earth's rotation rate will prohibit observation of this effect.

Wahr, J. M.

Protracted core formation and impact disruptions shaped the earliest outer Solar System planetesimals

The distinct compositions of metallic cores from noncarbonaceous (NC) and carbonaceous (CC) iron meteorite parent bodies (IMPBs) reflect differences in accretion and differentiation histories of earliest inner and outer Solar System planetesimals. Compared to NC IMPBs, CC IMPBs have smaller, sulfur (S)–poor, highly siderophile element (HSE)–enriched cores and younger core formation ages. However, the origins of these differences remain debated. Using equilibrium partitioning models between the S-poor solid and the S-rich liquid metal, we argue that HSE enrichment in IID, IIF, IIIF, and IVB cores resulted from a multistage evolutionary sequence: (i) segregation of S-rich, HSE-depleted protocores during initial planetesimal heating; (ii) collisional disruption before S-poor metal segregation; (iii) reaccretion of mantle fragments into daughter planetesimals; and (iv) further 26 Al-driven heating producing HSE-enriched, S-poor cores. We suggest that iron meteorites from these CC IMPBs originate from such second-generation cores. Accounting for “missing” S-rich protocores helps reconcile several NC-CC IMPB differences and highlights the role of early collisional processing in shaping planetesimal chemical evolution.

Geosciences

Submillimeter observations of W3

A 35 arcsec resolution, 400 micron continuum map of the core of the W3 molecular cloud complex is presented. The map shows prominent peaks centered approximately at the positions of the infrared sources IRS 4 and IRS 5 with stronger-emission at the western (IRS 4) peak. The visual extinction (inferred from the submillimeter measurements) through the center of both peaks is large (A sub v = 200). A total mass for the core of 900 solar masses is estimated, consistent with dynamical estimates if the observed velocity gradient is a result of rotation. The peak H2 density is approximately 2.8 x 10 to the 5th/cu cm. The observations exclude the possibility that the shell H II regions, W3A-C, could be in pressure equilibrium with the surrounding cloud. The observations establish that the core of W3 contains several massive condensations surrounding the youngest stars in the region.

Jaffe, D. T.