Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “core modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

A New Global Core Plasma Model of the Plasmasphere

The Global Core Plasma Model (GCPM) is the first empirical model for thermal inner magnetospheric plasma designed to integrate previous models and observations into a continuous in value and gradient representation of typical total densities. New information about the plasmasphere, in particular, makes possible significant improvement. The IMAGE Mission Radio Plasma Imager (RPI) has obtained the first observations of total plasma densities along magnetic field lines in the plasmasphere and polar cap. Dynamics Explorer 1 Retarding Ion Mass Spectrometer (RIMS) has provided densities in temperatures in the plasmasphere for 5 ion species. These and other works enable a new more detailed empirical model of thermal in the inner magnetosphere that will be presented.

Gallagher, D. L.↗

An Iron-Rain Model for Core Formation on Asteroid 4 Vesta

Asteroid 4 Vesta is differentiated into a crust, mantle, and core, as demonstrated by studies of the eucrite and diogenite meteorites and by data from NASA's Dawn spacecraft. Most models for the differentiation and thermal evolution of Vesta assume that the metal phase completely melts within 20 degrees of the eutectic temperature, well before the onset of silicate melting. In such a model, core formation initially happens by Darcy flow, but this is an inefficient process for liquid metal and solid silicate. However, the likely chemical composition of Vesta, similar to H chondrites with perhaps some CM or CV chondrite, has 13-16 weight percent S. For such compositions, metal-sulfide melting will not be complete until a temperature of at least 1350 degrees Centigrade. The silicate solidus for Vesta's composition is between 1100 and 1150 degrees Centigrade, and thus metal and silicate melting must have substantially overlapped in time on Vesta. In this chemically and physically more likely view of Vesta's evolution, metal sulfide drops will sink by Stokes flow through the partially molten silicate magma ocean in a process that can be envisioned as "iron rain". Measurements of eucrites show that moderately siderophile elements such as Ni, Mo, and W reached chemical equilibrium between the metal and silicate phases, which is an important test for any Vesta differentiation model. The equilibration time is a function of the initial metal grain size, which we take to be 25-45 microns based on recent measurements of H6 chondrites. For these sizes and reasonable silicate magma viscosities, equilibration occurs after a fall distance of just a few meters through the magma ocean. Although metal drops may grow in size by merger with other drops, which increases their settling velocities and decreases the total core formation time, the short equilibration distance ensures that the moderately siderophile elements will reach chemical equilibrium between metal and silicate before metal drop merger becomes important. In this model, there must be at least 30 percent melting of the silicate phase when metal melting is complete, corresponding to a crust thickness of at least 30 kilometers on Vesta, consistent with Dawn gravity observations. Greater degrees of silicate melting and a correspondingly thicker crust are possible if Vesta accreted sufficiently rapidly.

Kiefer, Walter S.↗

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

Response functions and giant monopole resonances for light to medium-mass nuclei from the ab initio symmetry-adapted no-core–shell model

Using the ab initio symmetry-adapted no-core–shell model, we compute sum rules and response functions for light to medium-mass nuclei, starting from interactions that are derived in the chiral effective field theory. Specifically, we investigate electromagnetic transitions of monopole, dipole and quadrupole nature for 4 He, and explore dominant features of giant monopole resonances in symmetric nuclei such as the closed-shell 4 He and 16 O light nuclei, the intermediate-mass open-shell 20 Ne and the medium-mass closed-shell 40 Ca. Furthermore, for the NNLO opt chiral potential, we determine parameter-free monopole sum rules, which can provide information on the incompressibility of symmetric nuclear matter. Here, we report 213(10) MeV as an estimate for the compression modulus for infinite nuclear matter, which overlaps with the lower range of values often used in current astrophysical applications.

ab initio nuclear structure↗

Relevence of 3D Simulations and Sandwich Core Topology for the Modeling of Honeycomb Core Sandwich Panels Undergoing Mode I Delamination

A recently developed cohesive zone traction-separation law, which includes the effects of fiber bridging in a novel way, is extended from 2D to 3D. The proposed cohesive model is applied to low fidelity (i.e. homogenized core representation) and high fidelity (i.e. directly accounting for the core topology) finite element models of a composite panel comprised of carbon fiber reinforced plastic facesheets and a honeycomb sandwich core. This enables the investigation of 2D to 3D parameter transferability, width-dependent effects such as thumbnail-shaped crack growth, and the verification of plane strain / plane stress assumptions. A pronounced curvature of the initial interface-related crack front is observed, while the bridging-related crack front is straight. Furthermore, it is found that the cohesive parameters can easily be transferred from 2D to 3D under plane stress assumptions, but not under plane strain assumptions. The numerical predictions are compared to experimental load-displacement and R-curves.

Composites↗

Divertor-safe nonlinear burn control based on a SOLPS parameterized core-edge model for ITER

Abstract For ITER operations, the range of desirable burning-plasma regimes with high fusion power output will be restricted by various operational constraints. These constraints include the saturation of ITER’s various heating and fueling actuators such as the neutral beam injectors, the ion and electron cyclotron heating systems, the gas puffing system, and the deuterium–tritium pellet injectors. In addition to these actuator constraints, the H-mode power threshold, divertor detachment, and the heat load on the divertor targets may apply limitations to ITER’s operational space. In this work, Plasma Operation Contour (POPCON) plots that map the aforementioned constraints to the temperature-density space are used to investigate which constraints are most limiting towards accessing regimes with high fusion power output. The presented POPCON plots are based on a control-oriented core-edge model that couples the nonlinear density and energy response models for the core-plasma region with SOLPS4.3 parameterizations for conditions in the edge-plasma regions (scrape-off-layer and divertor). Using this control-oriented core-edge model, a nonlinear burn controller, which aims to regulate the plasma temperature and density in the core-plasma region, is constructed in this work. This controller is augmented with an online optimization scheme that governs the control references such that the plasma can be guided towards regimes with high fusion powers while protecting the divertor targets from dangerously high heat loads. A closed-loop simulation study illustrates the capability of this burn control scheme.

Physics↗

Integrating Remote Sensing Observations with NASA's GEOS-5 Modeling Framework in Support of Retrospective Analyses and Seasonal Prediction of Biosphere-Atmosphere CO2 Flux

The absorption of nearly half of all anthropogenic carbon dioxide (CO2) emissions by terrestrial and marine ecosystems has played a critical role in mitigating climate change. However, a persistent lack of understanding of these uptake processes impedes attribution, and thereby introduces large uncertainties into projections of the trajectory of carbon and climate in coming decades. A better understanding and quantification of the seasonal and interannual variability of carbon flux between the atmosphere and ocean and terrestrial biospheres is critically needed to improve predictions from Earth system models. We propose to bring together a diverse suite of remote sensing observations, mechanistic models, and NASA's Goddard Earth Observing System, version 5 (GEOS-5) model to address two objective 1) quantify the processes controlling the temporal variability of biosphere-atmosphere CO2 flux on seasonal and interannual timescales from 1989 to 2018 and 2) evaluating the predictability of the components of carbon flux on seasonal timescales. Because GEOS-5 has been developed to incorporate a variety of Earth system observations and includes a complete and physically consistent depiction of the atmosphere, land, and ocean carbon cycles, it represents a unique and ideal synthesis framework for this effort. By expanding the realism of processes that could be included in GEOS-5 seasonal forecasts, this effort supports NASA's goals toward an integrated Earth system modeling and prediction framework.We propose to integrate a number of observations into the GEOS-5 modeling system to refine retrospective estimates of both terrestrial and marine carbon flux and to test hypotheses regarding the effects of climate drivers on terrestrial carbon flux. This will include assimilating several satellite ocean color products in the NASA Ocean Biogeochemistry Model (NOBM) to improve the representation of ocean biology. We will also utilize remote sensing datasets to increase the temporal resolution of global land-use change estimates and fire emissions to better constrain the temporal variability of the terrestrial carbon cycle. This information will be incorporated into two NASA terrestrial biosphere models Catchment-CN, a dynamic vegetation model that is integrated within GEOS-5 and the Ecosystem Demography (ED) model, a core modeling resource supporting the Global Ecosystem Dynamics Investigation (GEDI) mission to map biomass. Fluxes will be evaluated through 1) comparison to top-down flux estimates inferred from surface CO2 observations and 2) comparison against available aircraft and satellite observations.By combining predictions of net ecosystem exchange and air-sea CO2 flux from land and ocean models with statistical predictions of land use change, fire, and fossil fuel emissions, we propose to create the first dynamical CO2 forecast on seasonal timescales. Because of the experimental nature of this endeavor, an important project objective will be to evaluate the limits to the predictability of different carbon cycle processes over lead times of 1-12 months using different satellite and surface data streams.

Ott, Lesley↗

Empirical Modeling of Core Plasma in the Inner Magnetosphere

This report presents first year progress on "Empirical Modeling of Core Plasma in the Inner Magnetosphere". In this subtask we are developing a three-dimensional model of plasma outflow from the high latitude ionosphere to the magnetosphere. Much of the main framework of the model has now been put in place.

Wilson, Gordon R.↗

Core-Formation Models and Extinct Nuclides

Zr and W isotope data are consistent with the Earth's core forming in a single event subsequent to about 113 Ma after the formation of the solar system. With continuous models of core formation the process can start early. Additional information is contained in the original extended abstract.

Jacobsen, S. B.↗

Modelling the core magnetic field of the earth

It is suggested that radial off-center dipoles located within the core of the earth be used instead of spherical harmonics of the magnetic potential in modeling the core magnetic field. The off-center dipoles, in addition to more realistically modeling the physical current systems within the core, are if located deep within the core more effective at removing long wavelength signals of either potential or field. Their disadvantage is that their positions and strengths are more difficult to compute, and such effects as upward and downward continuation are more difficult to manipulate. It is nevertheless agreed with Cox (1975) and Alldredge and Hurwitz (1964) that physical realism in models is more important than mathematical convenience. A radial dipole model is presented which agrees with observations of secular variation and excursions.

Harrison, C. G. A.↗

Plasmaspheric Density Troughs: Global IMAGE EUV Observations and Analysis via Global Core Plasma Modeling

To date, the IMAGE EUV camera has observed several plasmaspheric density trough features inside the plasmapause under a wide range of geomagnetic activity. From the perspective of EUV, a density trough feature appears as a channel of diminished pixel counts which spans a width of L-shell (DELTA L) and magnetic local time (MLT) inside the plasmapause. Plasmaspheric density troughs are found to be morphologically complex possessing considerable spatial and temporal variability. We present an analysis of the evolution of trough DELTA L and MLT extent as functions of associated D (sub ST) and K (sub p) history. Trough features range in size from 0.16 less than or equal to DELTA L less than or equal to 1.2 with azimuthal extent from 1500 less than or equal to MLT less than or equal to 1200. All cases of plasmaspheric density troughs studied to date appear to have evolved as a result of the inner edge of the afternoon/evening plasma drainage plume being wrapped around through the nightside plasmasphere. The structure of plasmaspheric density trough features is further probed by analyzing simulated EUV images produced by forward modeling artificially introduced regions of depleted density into both static and dynamic global core plasmaspheric models. Forward modeling suggests that (1) L-shell refilling of density troughs during storm recovery can be modeled as filling from the ionosphere toward the equator (i.e., bottom-up refilling), and (2) that an erosion process is operating within flux tubes beyond the outer L-shell wall of the observed density troughs.

Adrian, M. L.↗

A model for core formation in the early Earth

Two basic types exogenous models were proposed to account for siderophile and chalcophile element abundances in the Earth's upper mantle. The first model requires that the Earth be depleted in volatiles and that, after a core formation event which extracted the most siderophile elements into the core, additional noble siderophile elements (Pt, Ir, Au) were added as a late veneer and mixed into the mantle. The second model postulates a reduced Earth with approximately CI elemental abundances in which a primary core forming event depleted all siderophile elements in the mantle. The plausibility of models which require fine scale mixing of chondritic material into the upper mantle is analyzed. Mixing in liquids is more efficient, but large degrees of silicate partial melting will facilitate the separation of magma from residual solids. Any external events affecting the upper mantle of the Earth should also be evident in the Moon; but siderophile and chalcophile element abundance patterns inferred for the mantles of the Earth and Moon differ. There appear to be significant physical difficulties associated with chondritic veneer models.

Jones, J. H.↗

Study of the muon component in the core-corona model using CONEX 3D

The discrepancy between models and data regarding the muon content in air showers generated by ultra-high energy cosmic rays still needs to be solved. The CONEX simulation framework provides a flexible tool to assess the impact of different interaction properties and thus address the muon puzzle. In this work, we present the multidimensional extension of CONEX and show its performance compared to CORSIKA by discussing muon-related air-shower features for three experiments: KASCADE, IceTop, and the Pierre Auger Observatory. We also implement an effective version of the core-corona model to demonstrate the impact of the core effect, as observed at the LHC, on the muon content in air showers produced by ultra-high energy cosmic rays. At a primary energy of E$_{0}$ = 10$^{19}$ eV, we obtain an increase of 15% to 20% in the muon content.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

TRANSP-TGLF core predictive modeling of the JET DT baseline scenario

In recent years, an intense modeling activity has been focused on preparing and analyzing the second JET Deuterium–Tritium (D–T) experimental campaign DTE2. Among the numerous scientific outcomes of this campaign was the unique opportunity to test and validate the state of the art of modeling tools with fusion-relevant DT plasmas using the full metallic ITER-like wall in different scenarios. This work reports on the core predictive modeling of plasma density, electron and ion temperatures (n e , T e , T i ) performed using TRANSP (Pankin et al 2025 Comput. Phys. Commun. 312 109611) coupled with Trapped Gyro-Landau Fluid (TGLF)-SAT2 (Kinsey et al 2008 Phys. Plasmas 15 055908), (Staebler et al 2021 Nucl. Fusion 61 116007) for the JET D–T baseline scenario (I p = 3.5 MA, q 95 = 3, β N < 2, with pellet pacing) (Garzotti et al 2025 Plasma Phys. Control. Fusion 67 075011). The sensitivity to different input parameters, as the $\vec{E}$ x $\vec{B}$ shear parameterization and the values of the kinetic quantities at the boundary of the prediction domain (ρ = 0.85) has been assessed, identifying the confidence interval of the prediction results. In particular, the dependence of the electron density profile on the particle source parameters has been studied, identifying the ionization source as the main cause for the density gradient under-prediction obtained by TRANSP-TGLF, reported in (Hyun-Tae et al 2023 Nucl. Fusion 63 112004).

JET↗

A post-core-collapse model for the nucleus of M33

A model for the post-collapse core of M33 is proposed. By analogy with globular clusters, the formation of tight binary systems powers the reexpansion of the core, and it is expected that low-mass X-ray binaries form at a rate about equal to that of all galactic globular clusters combined. About a dozen such binaries should be present, and their combined emission may explain the large and enigmatic unresolved X-ray emission from the nucleus of M33. In addition, tidal formation of cataclysmic binaries may lead to a nova rate of more than about one per century. Numerous blue stragglers have probably formed through mergers caused by stellar collisions, but their density is probably too low to explain the blue color of the nucleus. Young stars are the likely cause for the color, and their presence may complicate the simple dynamical picture presented here.

Hernquist, Lars↗

Enhanced Core Noise Modeling for Turbofan Engines

This report describes work performed by MTC Technologies (MTCT) for NASA Glenn Research Center (GRC) under Contract NAS3-00178, Task Order No. 15. MTCT previously developed a first-generation empirical model that correlates the core/combustion noise of four GE engines, the CF6, CF34, CFM56, and GE90 for General Electric (GE) under Contract No. 200-1X-14W53048, in support of GRC Contract NAS3-01135. MTCT has demonstrated in earlier noise modeling efforts that the improvement of predictive modeling is greatly enhanced by an iterative approach, so in support of NASA's Quiet Aircraft Technology Project, GRC sponsored this effort to improve the model. Since the noise data available for correlation are total engine noise spectra, it is total engine noise that must be predicted. Since the scope of this effort was not sufficient to explore fan and turbine noise, the most meaningful comparisons must be restricted to frequencies below the blade passage frequency. Below the blade passage frequency and at relatively high power settings jet noise is expected to be the dominant source, and comparisons are shown that demonstrate the accuracy of the jet noise model recently developed by MTCT for NASA under Contract NAS3-00178, Task Order No. 10. At lower power settings the core noise became most apparent, and these data corrected for the contribution of jet noise were then used to establish the characteristics of core noise. There is clearly more than one spectral range where core noise is evident, so the spectral approach developed by von Glahn and Krejsa in 1982 wherein four spectral regions overlap, was used in the GE effort. Further analysis indicates that the two higher frequency components, which are often somewhat masked by turbomachinery noise, can be treated as one component, and it is on that basis that the current model is formulated. The frequency scaling relationships are improved and are now based on combustor and core nozzle geometries. In conjunction with the Task Order No. 10 jet noise model, this core noise model is shown to provide statistical accuracy comparable to the jet noise model for frequencies below blade passage. This model is incorporated in the NASA FOOTPR code and a user s guide is provided.

Stone, James R.↗

Ab initio leading order effective potential for elastic proton scattering based on the symmetry-adapted no-core shell model

Calculating microscopic optical potentials for elastic scattering at intermediate energies from light nuclei in an ab initio fashion within the Watson expansion has been established within the last few years. Based on the Watson expansion of the multiple scattering series, we employ a nonlocal translationally invariant nuclear density derived within the symmetry-adapted no-core shell model (SA-NCSM) framework from a chiral next-to-next-to-leading order (NNLO) nucleon-nucleon interaction and the very same interaction for a consistent full-folding calculation of the effective (optical) potential for nucleon-nucleus scattering for medium-heavy nuclei. The leading order effective (optical) folding potential is computed by integrating over a translationally invariant SA-NCSM one-body scalar density, spin-projected momentum distribution, and the Wolfenstein amplitudes 𝐴, 𝐶, and 𝑀. The resulting nonlocal potentials serve as input for a momentum space Lippmann-Schwinger equation. In the SA-NCSM, the model space is systematically up-selected using Sp⁡(3,ℝ) symmetry considerations. For the light nucleus of 6 He, we establish a systematic selection scheme in the SA-NCSM for scattering observables. Then, we apply this scheme to calculations of scattering observables, such as differential cross sections, analyzing powers, and spin rotation functions for elastic proton scattering from 20 Ne and 40 Ca in the energy regime between 65 and 200 MeV, and compare to available data. Furthermore, our calculations show that the leading order effective nucleon-nucleus potential in the Watson expansion of multiple scattering theory obtained from an up-selected SA-NCSM model space describes 40 Ca elastic scattering observables reasonably well to about 60 degrees in the center-of-mass frame, which coincides roughly with the validity of the NNLO chiral interaction used to calculate both the nucleon-nucleon amplitudes and the one-body scalar and spin nuclear densities.

Ab initio calculations↗