SEARCH · Engineering Papers

Results for “Calculational Margin”

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.

SourceResearch connectionFrom yearThrough year

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

At least 487 records · Page 27

All sources NASA NTRS DOE OSTI

NASA NTRS2022

Quantifying Feedback from Narrow Line Region Outflows in Nearby Active Galaxies. IV. The Effects of Different Density Estimates on the Ionized Gas Masses and Outflow Rates

Active galactic nuclei (AGN) can launch outflows of ionized gas that may influence galaxy evolution, and quantifying their full impact requires spatially resolved measurements of the gas masses, velocities, and radial extents. We previously reported these quantities for the ionized narrow-line region outflows in six low-redshift AGN, where the gas velocities and extents were determined from Hubble Space Telescope long-slit spectroscopy. However, calculating the gas masses required multicomponent photoionization models to account for radial variations in the gas densities, which span ∼6 orders of magnitude. To simplify this method for larger samples with less spectral coverage, we compare these gas masses with those calculated from techniques in the literature. First, we use a recombination equation with three different estimates for the radial density profiles. These include constant densities, those derived from [S ii], and power-law profiles based on constant values of the ionization parameter (U). Second, we use single-component photoionization models with power-law density profiles based on constant U, and allow U to vary with radius based on the [O iii]/Hβ ratios. We find that assuming a constant density of nH = 10(^2) cm(^−3) overestimates the gas masses for all six outflows, particularly at small radii where the outflow rates peak. The use of [S ii] marginally matches the total gas masses, but also overestimates at small radii. Overall, single-component photoionization models where U varies with radius are able to best match the gas mass and outflow rate profiles when there are insufficient emission lines to construct detailed models.

Mitchell Revalski↗

DOE OSTI2023

Can upscaling ground nadir SIF to eddy covariance footprint improve the relationship between SIF and GPP in croplands?

Ground solar-induced chlorophyll fluorescence (SIF) is important for the mechanistic understanding of the dynamics of vegetation gross primary production (GPP) at fine spatiotemporal scales. However, eddy covariance (EC) observations generally cover larger footprint areas than ground SIF observations (a bare fiber with nadir), and this footprint mismatch between nadir SIF and GPP could complicate the canopy SIF-GPP relationships. Here, we upscaled nadir SIF observations to EC footprint and investigated the change in SIF-GPP relationships after the upscaling in cropland. We included 13 site-years data in our study, with seven site-years corn, four site-years soybeans, and two site-years miscanthus, all located in the US Corn Belt. All sites’ crop nadir SIF observations collected from the automated FluoSpec2 system (a hemispheric-nadir system) were upscaled to the GPP footprint-based SIF using vegetation indices (VIs) calculated from high spatiotemporal satellite reflectance data. We found that SIF-GPP relationships were not substantially changed after upscaling nadir SIF to GPP footprint at our crop sites planted with corn, soybean, and miscanthus, with R 2 change after the upscaling ranging from -0.007 to 0.051 and root mean square error (RMSE) difference from -0.658 to 0.095 umol m –2 s –1 relative to original nadir SIF-GPP relationships across all the site-years. The variation of the SIF-GPP relationship within each species across different site-years was similar between the original nadir SIF and upscaled SIF. Different VIs, EC footprint models, and satellite data led to marginal differences in the SIF-GPP relationships when upscaling nadir SIF to EC footprint. Furthermore, our study provided a methodological framework to correct this spatial mismatch between ground nadir SIF and GPP observations for croplands and potentially for other ecosystems. Our results also demonstrated that the spatial mismatch between ground nadir SIF and GPP might not significantly affect the SIF-GPP relationship in cropland that are largely homogeneous.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Confronting interacting dark radiation scenarios with cosmological data

Dark radiation (DR) is generally predicted in new physics scenarios that address fundamental puzzles of the Standard Model or tensions in the cosmological data. Cosmological data have the sensitivity to constrain not only the energy density of DR but also whether it is interacting. In this paper, we present a systematic study of five types of interacting DR (free-streaming, fluid, decoupling, instantaneous decoupling, and recoupling DR) and their impact on cosmological observables. We modify the Boltzmann hierarchy to describe all these types of interacting DR under the relaxation time approximation. We, for the first time, robustly calculate the collision terms for recoupling scalar DR and provide a better estimation of the recoupling transition redshift. We demonstrate the distinct features of each type of DR on the cosmic microwave background and matter power spectra. We perform Markov-chain Monte Carlo scans using the Planck 2018 data and baryon acoustic oscillation data. Assuming no new physics in the standard model neutrino sector, we find no statistically significant constraints on the couplings of DR, although there is a slight preference for the fluidlike limit of all the cases. In the case of instantaneous decoupling DR, this limit corresponds to a late transition redshift around recombination. Furthermore, the ΔN eff constraint varies marginally depending on the type of DR.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Impact of uranium oxide (UO 2 ) fuel with molybdenum (Mo) inserts on pressurized water reactor performance and safety

This work investigates nuclear reactor performance and safety characteristics of UO 2 with high thermal conductivity Mo insert structures by using multiphysics modeling techniques. Additionally, the purpose of this study is to use scoping analyses to quantify the impact of using Mo inserts from neutronic and heat transfer standpoints. Attention is given to reactor performance parameters, such as cycle length, maximum fuel temperature, temperature gradients in the fuel, and stored energy in the fuel. The finite-element code BISON and the Monte Carlo particle transport code Serpent were used to perform sensitivity analyses on the Mo insert geometry to optimize the insert design and inform larger scale modeling that required the homogenization of the UO 2 and Mo. Although BISON is often used as a fuel performance analysis tool, it is used in this context for heat transfer analysis only. Fuel performance optimization is outside the scope of the current study, but would be important for future work focused on this concept. The results showed that the insert had little impact on neutronic performance and that homogenizing the UO 2 and Mo was acceptable for reactor physics calculations. Reactivity temperature coefficients calculated using homogeneous UO 2 -Mo were shown to be relatively similar to UO 2 , but higher Mo content and 235 U enrichment can reduce the worth of soluble boron and control rods. The effect of insert geometry on heat transfer was much greater, and an approximately 15–20% difference in maximum fuel temperature was predicted between the best and worst performing heat transfer geometries. Furthermore, thermal conductivity calibration based on the finite element analysis results was performed to improve the accuracy of temperature predictions in reactor analysis models that homogenized the UO 2 -Mo fuel. Compared with UO 2 in a pressurized water reactor (PWR), the optimized UO 2 -Mo design increased the margin to fuel melt by 13–32% across the fuel cycle, but it requires the 235U enrichment to exceed 5% to match the cycle length of conventional UO 2 .

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS2000

Depositional facies and aqueous-solid geochemistry of travertine-depositing hot springs (Angel Terrace, Mammoth Hot Springs, Yellowstone National Park, U.S.A.)

Petrographic and geochemical analyses of travertine-depositing hot springs at Angel Terrace, Mammoth Hot Springs, Yellowstone National Park, have been used to define five depositional facies along the spring drainage system. Spring waters are expelled in the vent facies at 71 to 73 degrees C and precipitate mounded travertine composed of aragonite needle botryoids. The apron and channel facies (43-72 degrees C) is floored by hollow tubes composed of aragonite needle botryoids that encrust sulfide-oxidizing Aquificales bacteria. The travertine of the pond facies (30-62 degrees C) varies in composition from aragonite needle shrubs formed at higher temperatures to ridged networks of calcite and aragonite at lower temperatures. Calcite "ice sheets", calcified bubbles, and aggregates of aragonite needles ("fuzzy dumbbells") precipitate at the air-water interface and settle to pond floors. The proximal-slope facies (28-54 degrees C), which forms the margins of terracette pools, is composed of arcuate aragonite needle shrubs that create small microterracettes on the steep slope face. Finally, the distal-slope facies (28-30 degrees C) is composed of calcite spherules and calcite "feather" crystals. Despite the presence of abundant microbial mat communities and their observed role in providing substrates for mineralization, the compositions of spring-water and travertine predominantly reflect abiotic physical and chemical processes. Vigorous CO2 degassing causes a +2 unit increase in spring water pH, as well as Rayleigh-type covariations between the concentration of dissolved inorganic carbon and corresponding delta 13C. Travertine delta 13C and delta 18O are nearly equivalent to aragonite and calcite equilibrium values calculated from spring water in the higher-temperature (approximately 50-73 degrees C) depositional facies. Conversely, travertine precipitating in the lower-temperature (< approximately 50 degrees C) depositional facies exhibits delta 13C and delta 18O values that are as much as 4% less than predicted equilibrium values. This isotopic shift may record microbial respiration as well as downstream transport of travertine crystals. Despite the production of H2S and the abundance of sulfide oxidizing microbes, preliminary delta 34S data do not uniquely define the microbial metabolic pathways present in the spring system. This suggests that the high extent of CO2 degassing and large open-system solute reservoir in these thermal systems overwhelm biological controls on travertine crystal chemistry.

Non-NASA Center↗

NASA NTRS2007

Can the Equivalent Sphere Model Approximate Organ Doses in Space?

For space radiation protection it is often useful to calculate dose or dose,equivalent in blood forming organs (BFO). It has been customary to use a 5cm equivalent sphere to. simulate the BFO dose. However, many previous studies have concluded that a 5cm sphere gives very different dose values from the exact BFO values. One study [1] . concludes that a 9 cm sphere is a reasonable approximation for BFO'doses in solar particle event environments. In this study we use a deterministic radiation transport [2] to investigate the reason behind these observations and to extend earlier studies. We take different space radiation environments, including seven galactic cosmic ray environments and six large solar particle events, and calculate the dose and dose equivalent in the skin, eyes and BFO using their thickness distribution functions from the CAM (Computerized Anatomical Man) model [3] The organ doses have been evaluated with a water or aluminum shielding of an areal density from 0 to 20 g/sq cm. We then compare with results from the equivalent sphere model and determine in which cases and at what radius parameters the equivalent sphere model is a reasonable approximation. Furthermore, we address why the equivalent sphere model is not a good approximation in some cases. For solar particle events, we find that the radius parameters for the organ dose equivalent increase significantly with the shielding thickness, and the model works marginally for BFO but is unacceptable for the eye or the skin. For galactic cosmic rays environments, the equivalent sphere model with an organ-specific constant radius parameter works well for the BFO dose equivalent, marginally well for the BFO dose and the dose equivalent of the eye or the skin, but is unacceptable for the dose of the eye or the skin. The ranges of the radius parameters are also being investigated, and the BFO radius parameters are found to be significantly, larger than 5 cm in all cases, consistent with the conclusion of an earlier study [I]. The radius parameters for the dose equivalent in GCR environments are approximately between 10 and I I cm for the BFO, 3.7 to 4.8 cm for the eye, and 3.5 to 5.6 cm for the skin; while the radius parameters are between 10 and 13 cm for the BFO dose.

Lin, Zi-Wei↗

DOE OSTI2024

Theoretical modeling of a bottom-raised oscillating surge wave energy converter structural loadings and power performances

Here, this study presents theoretical formulations to evaluate the fundamental parameters and performance characteristics of a bottom-raised oscillating surge wave energy converter (OSWEC) device. Employing a flat plate assumption and potential flow formulation in elliptical coordinates, closed-form equations for the added mass, radiation damping, and excitation forces/torques in the relevant pitch-pitch and surge-pitch directions of motion are developed and used to calculate the system's response amplitude operator and the forces and moments acting on the foundation. The model is benchmarked against numerical simulations using WAMIT and WEC-Sim, showcasing excellent agreement. The sensitivity of plate thickness on the analytical hydrodynamic solutions is investigated over several thickness-to-width ratios ranging from 1:80 to 1:10. The results show that as the thickness of the benchmark OSWEC increases, the deviation of the analytical hydrodynamic coefficients from the numerical solutions grows from 3% to 25%. Differences in the excitation forces and torques, however, are contained within 12%. While the flat plate assumption is a limitation of the proposed analytical model, the error is within a reasonable margin for use in the design space exploration phase before a higher-fidelity (and thus more computationally expensive) model is employed. A parametric study demonstrates the ability of the analytical model to quickly sweep over a domain of OSWEC dimensions, illustrating the analytical model's utility in the early phases of design.

13 HYDRO ENERGY↗

DOE OSTI2022

Uncertainty propagation in pore water chemical composition calculation using surrogate models

Performance assessment in deep geological nuclear waste repository systems necessitates an extended knowledge of the pore water chemical conditions prevailing in host-rock formations. In the last two decades, important progress has been made in the experimental characterization and thermodynamic modeling of pore water speciation, but the influence of experimental artifacts and uncertainties of thermodynamic input parameters are seldom evaluated. In this respect, we conducted an uncertainty propagation study in a reference geochemical model describing the pore water chemistry of the Callovian-Oxfordian clay formation. Nineteen model input parameters were perturbed, including those associated to experimental characterization (leached anions, exchanged cations, cation exchange selectivity coefficients) and those associated to generic thermodynamic databases (solubilities). A set of 13 quantities of interest were studied by the use of polynomial chaos expansions built non-intrusively with a least-squares forward stepwise regression approach. Training and validation sets of simulations were carried out using the geochemical speciation code PHREEQC. The statistical results explored the marginal distribution of each quantity of interest, their bivariate correlations as well as their global sensitivity indices. The influence of the assumed distributions for input parameters uncertainties was evaluated by considering two parametric domain sizes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Scrambling versus relaxation in Fermi and non-Fermi liquids

In this work, we compute the Lyapunov exponent characterizing quantum scrambling in a family of generalized Sachdev- Ye-Kitaev models, which can be tuned between different low temperature states from Fermi liquids, through non-Fermi liquids, to fast scramblers. The analytic calculation, controlled by a small coupling constant and large $\textit{N}$, allows us to clarify the relations between the quasiparticle relaxation rate $\textit{1/τ}$ and the Lyapunov exponent $λ_L$ characterizing scrambling. In the Fermi liquid states we find that the quasiparticle relaxation rate dictates the Lyapunov exponent. In non-Fermi liquids, where $1/τ \gg T$, we find that $λ_L$ is always $\textit{T}$ linear with a prefactor that is independent of the coupling constant in the limit of weak coupling. Instead it is determined by a scaling exponent that characterizes the relaxation rate. $λ_L$ approaches the general upper bound 2$\textit{πT}$ at the transition to a fast scrambling state. Finally, in a marginal Fermi liquid state the exponent is linear in temperature with a prefactor that vanishes as a nonanalytic function ~$g\text{ln}(1/g)$ of the coupling constant $\textit{g}$.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1987

Dendritic growth in a supercooled alloy melt

A simple model which describes the growth of an 'array' of dendrites into a supercooled, binary, alloy melt is presented. Solute diffusion is calculated by superposing the solutions given by Flemings and Zener, and also, by superposing the solutions given by Ivantsov and Flemings. A general expression for the transport solution is suggested from which all other dendrite growth models presented earlier may be obtained as special cases. It is shown that both 'free' and 'constrained' growth may be described by a single transport solution, which indicates that (1) both thermal and solutal effects will be important during 'free' growth in dilute alloys, (2) only solutal effects are predominant during 'free' growth in concentrated alloys and during 'constrained' growth. An examination of the relevant dimensionless parameters also suggests that all dendrite growth models, regardless of the assumptions used to determine the tip radius (marginal stability, minimum undercooling, maximum velocity, minimum entropy production) should predict the experimentally observed extrema in tip radius and growth velocity in dilute alloys, during 'free' dendritic growth. Experimental data in binary H2O-NaCl and succinonitrile-acetone solutions are shown to be in good agreement with the model.

Laxmanan, V.↗

NASA NTRS2013

Component Testing of the J-2X Augmented Spark Igniter (ASI)

In support of the development of the J-2X engine, 201 low pressure, liquid oxygen / liquid hydrogen (LOX/LH2) J-2X Augmented Spark Igniter (ASI) subsystem ignition tests were conducted at Marshall Space Flight Center (MSFC). The main objective of these tests was to start the ASI within the anticipated J-2X engine start box, as well as outside of it, to check for ignition margin. The setup for the J-2X ASI component testing simulated, as much as possible, the tank-head start-up configuration of the ASI within the J-2X Engine. The ignition tests were divided into 124 vacuum start tests to simulate altitude start on a flight engine, and 77 sea-level start tests to simulate the first set of ground tests for the J-2X Engine at Stennis Space Center (SSC). Other ignition parameters that were varied included propellant tank pressures, oxidizer temperature entering the ASI oxidizer feedline, oxidizer valve timing, spark igniter condition (new versus damaged), and oxidizer and fuel feedline orifice sizes. Propellant blowdowns using venturis sized to simulate the ASI resistance allowed calculation of transient propellant mass flow rates as well as global mixture ratio for all ignition tests. Global mixture ratio within the ASI at the time of ignition varied from 0.2 to 1.2. Detailed electronics data obtained from an instrumented ignition lead allowed characterization of the breakdown voltage, sustaining voltage and energy contained in each spark as the ASI propellants ignited. Results indicated that ignition always occurred within the first five sparks when both propellants were present in the ASI chamber.

Osborne, Robin J.↗

DOE OSTI2023

Evaluation of Thermal Neutron Scattering Cross Section of Uranium Silicide with Ab Initio Lattice Dynamics

Uranium silicide (U 3 Si 2 ) is a candidate material for the high-density nuclear fuel in commercial light water reactors [1], [2]. Its higher uranium density, 11.3 g-U/cm3, compared to that of uranium dioxide (UO 2 ), 9.7 g-U/cm3, can improve the performance of a nuclear reactor while using low enriched uranium (LEU) and diversify the choice of cladding materials [1]–[3]. It also has a higher thermal conductivity than UO 2 , which can reduce the thermal stress on the material caused by a temperature gradient across the fuel pellet and provide a larger margin for some postulated accidents [1], [2], [4]. Furthermore, compared to U3Si, another high-density fuel candidate, it has better resistance to in-pile swelling due to less irradiation-induced rapid amorphization [1], [3]. Corresponding to its importance in nuclear engineering, many previous studies have reported the properties of U3Si2. Experiments showed that U 3 Si 2 is a paramagnetic (PM) metal, where a slight linear increase in magnetic susceptibility was measured with increasing temperature [5], [6]. In addition, thermodynamic quantities such as thermal expansion coefficient, heat capacity, and thermal conductivity were experimentally determined over a wide temperature range [1], [7], [8]. In several computational studies, ab initio atomistic simulations based on density functional theory (DFT) were performed to calculate various properties including elastic constants, electronic density of states (DOS), and phonon dispersion curves [9]–[12]. Nevertheless, thermal neutron scattering cross sections, which are critical to the prediction of the parameters in reactor physics that are ultimately related to reactor criticality, have not yet been evaluated for U3Si2. The scattering cross section can be calculated from the phonon DOS, or the energy spectrum of lattice vibrations, of the crystalline system [13], [14]. However, there is also no experimental data available for the phonon DOS of U 3 Si 2 . While some computational studies reported the phonon DOS and/or dispersion curves from ab initio simulations [9]–[12], the accuracy cannot be guaranteed because it is unclear whether the spin-polarization behavior of PM U 3 Si 2 was properly described. In the present study, the thermal neutron scattering cross section for U 3 Si 2 is evaluated for the first time by calculating the phonon DOS for U3Si2 from ab initio lattice dynamics (AILD) simulations based on DFT. First, U 3 Si 2 is modeled based on the experimental structure, and AILD simulations are performed on the modeled U3Si2 to optimize the structure. Next, AILD simulations are performed for supercells with atomic displacement to calculate Hellmann-Feynman forces. Based on the calculated forces, partial phonon DOSs for U and Si are obtained, and the thermal neutron scattering law (TSL) for U 3 Si 2 is finally evaluated. To verify the accuracy of the calculations in the present study, the calculation results are compared with experimental data on the structure and heat capacity of U3Si2 [1], [7], [8], [15].

Geometry Optimization↗

NASA NTRSDate not supplied

Magnetohydrodynamics (MHD) Aerocapture System for Enabling Faster-Larger Planetary Science & Human Exploration Missions

Since our completing the NIAC Phase I NIAC on this Advanced Aerocapture System, NASA Langley Research Center has funded or supported a number of studies and code enhancements through its Center Innovation Fund (CIF) and NASA’s NSTGRO and Internship Programs to mature the analysis capabilities and quantify the merits of the MHD Aerocapture System technology. These efforts have resulted in a plug and play analysis capability for assessing MHD aerocapture system performance for arrival at many planetary bodies of interest. Our efforts have especially focused on the potential mass savings for improving the capacity for science observations at Neptune and Triton. A re-cent Forbes article published “‘Orbital mechanics is probably going to decide for us whether we go to Uranus or Neptune because we need to flyby Jupiter,’ said Kunio Sayanagi at Hampton University, Virginia, who also worked on the Neptune Odyssey proposal…. Exactly when a mission can be sent to Uranus, or Neptune, depends on the relative position of Jupiter, which can help give a spacecraft a gravitational slingshot. That drastically shortens the cruise phase.” [1] Since shortening the cruise phase is important for these science missions, any mass savings enabled by the MHD Aerocapture System could be reallocated to increasing Thermal Protection System mass to allow faster arrival speeds and/or for onboarding additional payloads such as science instruments, batteries, or propellant for conducting more science for longer durations in the desired orbits. The analysis steps and codes for conducting trades and sizing vehicles for aerocapture are as follows: Step 1 is to conduct aeroheating analysis using LAURA of the selected entry vehicle shape to identify locations on the forebody where ionization and flow velocity are sufficient for producing Lo-rentz forces. LAURA is a multiblock structured grid finite-volume CFD solver developed at the NASA Langley Research Center. [2] LAURA has been used for aerothermal analysis support of the entry, de-scent and landing (EDL) phase of interplanetary missions over the last three decades [3-7]. Step 2 is to port the LAURA results into CFDWARP to calcu-late electrical and thermal conductivities of ionized flow for sizing MHD patch system and calculating Lorentz forces needed for controls analysis. CFDWARP is a CFD code that uses advanced nu-merical methods that enable the simulation of the full coupling between the aerodynamics, the magne-tohydrodynamics, and the non-neutral plasma sheaths. CFDWARP has the unique capability to simulate efficiently the non-neutral sheaths (near the electrodes) in coupled form with the quasi-neutral bulk MHD flow [8-11]. Step 3 is to link re-sults from LAURA and CFDWARP into POST2 for calculating entry trajectories and comparing MHD control results with other aerodynamic control strategies. The Program to Optimize Simulated Tra-jectories II (POST2) is a generalized point mass, discrete parameter targeting and optimization pro-gram. POST2 provides the capability to target and optimize point mass trajectories for multiple pow-ered or un-powered vehicles near an arbitrary rotat-ing, oblate planet [12]. Step 4: TPS sizing was per-formed using the Fully Implicit Ablation and Ther-mal-response code (FIAT) tool which computes the transient one-dimensional thermal response and surface thermochemistry of a multilayer stackup of thermal protection, bonding, and structural materi-als subject to aeroheating on one surface [13]. The sizing and margining methodology used was based on the approach documented by Mahzari and Milos [14] for the dual-layered heatshield for extreme entry environment technology (DL-HEEET) TPS concept. TPS analysis utilizes trajectory information from POST2. Using this step-wise plug and play MHD Aerocapture performance assessment process, our analysis targets a Neptune aerocapture trajectory that will place the spacecraft in an observation orbit for Triton. [15]. Magnetohydrodynamic (MHD) control of a 4.5-meter diameter MSL-style capsule resulted in TPS mass savings of nearly 2000 kg when using an MHD system mass of under 200 kg. The flight path for a vehicle using the MHD control strategy has a much lower heat rate and heat load compared to the conventional aerodynamic aerocapture strategies known as bank angle con-trolled (BAC) and direct force controlled (DFC). Both BAC and DFC have heat rates significantly greater than 1500 W/cm2 typically used as an upper limit for PICA. Thus, DL-HEEET TPS concept was required for the BAC and DFC control strategies. However, considering the more benign environ-ments for the MHD case, additional TPS concepts with improved mass efficiency were also assessed. PICA was considered for the MHD controlled strat-egy since the maximum heat rate was well within the limits (<1500 W/cm2) of PICA. TPS sizing re-sulted in a significant mass reduction. The PICA layer for this sizing case was about 7.8 cm. As a point of reference, the Mars 2020 mission, which used this same PICA concept, had a PICA thickness of 3.18 cm [16]. The trajectories used for the TPS sizing originat-ed from the POST2 simulations. The current, I, to an electromagnet configuration can be manipulated to allow for active control of the vehicle. Manipula-tion of the current, I, changes the magnetic field, B, which affects the Lorentz force and therefore the MHD drag force on the vehicle. Our analysis in-cluded both open-loop and close-loop control. Closed-loop control will enable improved overall performance when taking into account mission level uncertainties, such as interplanetary delivery errors and atmospheric modeling uncertainties. The open-loop and closed-loop MHD control cases do not dip as deep into the atmosphere as the aerodynamic cases. Three types of aerodynamic-only approach-es are investigated: bank angle modulation (BAM), director force control (DFC), and Drag Modulated. BAM and DFC make use of vehicle aerodynamic angles to steer the vehicle. Thus, changing the aer-odynamic forces acting on the vehicle for control, aerodynamic drag modulated case requires a vary-ing drag area to modulate the drag force. The MHD drag modulated case modulates MHD generated drag force that adds to the aerodynamic drag. This higher atmospheric activation of drag forces by the MHD patch results in significantly less heat flux on the vehicle. The MHD technology will enable shorter cruise times and deceleration of larger payloads for increasing the capacity for science at the Ice Giants or for returning astronauts to Earth from cislunar space or from Mars. The purpose of this presentation is to provide more details about this work and to highlight plans for further research and development including a flight demonstration.

R. W. Moses↗

DOE OSTI2023

Ab Initio Simulation of Structure and Properties in Ni-Based Superalloys: Haynes282 and Inconel740

The electronic structure, interatomic bonding, and mechanical properties of two supercell models of Ni-based superalloys are calculated using ab initio density functional theory methods. The alloys, Haynes282 and Inconel740, are face-centered cubic lattices with 864 atoms and eleven elements. These multi-component alloys have very complex electronic structure, bonding and partial-charge distributions depending on the composition and strength of the local bonding environment. We employ the novel concept of total bond order density (TBOD) and its partial components (PBOD) to ascertain the internal cohesion that controls the intricate balance between the propensity of metallic bonding between Ni, Cr and Co, and the strong bonds with C and Al. We find Inconel740 has slightly stronger mechanical properties than Haynes282. Both Inconel740 and Haynes282 show ductile natures based on Poisson’s ratio. Poisson’s ratio shows marginal correlation with the TBOD. Comparison with more conventional high entropy alloys with equal components are discussed.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Electron-phonon coupling and spin fluctuations in the Ising superconductor NbSe2

Abstract Ising superconductivity, observed in NbSe 2 and similar materials, has generated tremendous interest. Recently, attention was called to the possible role that spin fluctuations (SF) play in this phenomenon, in addition to the dominant electron–phonon coupling (EPC); the possibility of a predominantly triplet state was discussed and led to a conjecture of viable singlet–triplet Leggett oscillations. However, these hypotheses have not been put to a quantitative test. In this paper, we report first principle calculations of the EPC and also estimate coupling with SF, including full momentum dependence. We find that: (1) EPC is strongly anisotropic, largely coming from the $$K-{K}^{{\prime} }$$ K − K ′ scattering, and therefore excludes triplet symmetry even as an excited state; (2) superconductivity is substantially weakened by SF, but anisotropy remains as above; and, (3) we do find the possibility of a Leggett mode, not in a singlet–triplet but in an s ++ – s ± channel.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Attention-Augmented Parametric Kernel Graph Neural Network (APKGNN) for Node Classification

We present a new graph neural network, the Attention-based Parametric-Kernel augmented Graph Neural Network (APKGNN), developed for node classification tasks. Despite extensive work on modeling multi-faceted relationships between connected nodes of a graph, the effect of attention on edge features mapped to relationships has not yet been analyzed through learning representation. This study derives such an attention vector by first calculating node features corresponding to endpoints of an edge and then aggregating these with extracted local intrinsic patches of a given graph to generate augmented local patch vectors. This process uses a parametric kernel based on Gaussian mixture models (GMMs) to embed local neighborhoods of the graph in local patches. The patch vectors then convolve with the above node features to produce an updated node representation. We show that this new learning representation (APKGNN) achieves higher node classification accuracy on tasks - both standard benchmarks (Cora, PubMed, Citeseer) and new experimental short text corpora where nodes correspond to text documents and words. This implementation of the GNN convolution layer outperforms state-of-the-art (SOTA) algorithms, achieving higher training, validation, and test accuracy by a significant margin on three standard benchmark data sets under both SOTA experimental settings and those for new testbeds.

Bose, Avishek↗

NASA NTRS1982

A comparative study of computational methods in cosmic gas dynamics

Many theoretical investigations of fluid flows in astrophysics require extensive numerical calculations. The selection of an appropriate computational method is, therefore, important for the astronomer who has to solve an astrophysical flow problem. The present investigation has the objective to provide an informational basis for such a selection by comparing a variety of numerical methods with the aid of a test problem. The test problem involves a simple, one-dimensional model of the gas flow in a spiral galaxy. The numerical methods considered include the beam scheme, Godunov's method (G), the second-order flux-splitting method (FS2), MacCormack's method, and the flux corrected transport methods of Boris and Book (1973). It is found that the best second-order method (FS2) outperforms the best first-order method (G) by a huge margin.

Van Albada, G. D.↗

NASA NTRS1985

Fault type predictions from stress distributions on planetary surfaces - Importance of fault initiation depth

The prediction of fault type on planetary surfaces from model stresses calculated at depth is discussed. These fault-type predictions yield different faults than those predicted using the surface criteria commonly employed in geophysical models. For elastic-plate flexure models of mascon loading on the moon, stresses calculated at the surface predict the occurrence of strike-slip faulting at the radial distance where grabens are found. Normal faults bounding lunar grabens and thrust faults responsible for wrinkle ridges are analyzed. It is found that the former initiate at the mechanical discontinuity that separates the breccia of the megaregolith from in situ fractured rock and that the latter initiate at the mechanical discontinuity between basalt layers and the underlying basin floor. The difference between elastic constants for the outer few kilometers of brecciated megaregolith and the underlying lunar lithosphere are evaluated. Superposing nonisotropic stresses resulting from the weight of overburden to the depth of the relevant mechanical discontinuity yield stresses that predict wrinkle ridges in the basin centers and grabens outside the basin margin, and eliminate the predicted zone of strike-slip faults.

Golombek, M. P.↗