Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “First principles”

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 1,243 records · Page 69
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

The Effects of Magnetic Nozzle Configurations on Plasma Thrusters

Over the course of eight years, the Ohio State University has performed research in support of electric propulsion development efforts at the NASA Lewis Research Center, Cleveland, OH. This research has been largely devoted to plasma propulsion systems including MagnetoPlasmaDynamic (MPD) thrusters with externally-applied, solenoidal magnetic fields, hollow cathodes, and Pulsed Plasma Microthrusters (PPT's). Both experimental and theoretical work has been performed, as documented in four master's theses, two doctoral dissertations, and numerous technical papers. The present document is the final report for the grant period 5 December 1987 to 31 December 1995, and summarizes all activities. Detailed discussions of each area of activity are provided in appendices: Appendix 1 - Experimental studies of magnetic nozzle effects on plasma thrusters; Appendix 2 - Numerical modeling of applied-field MPD thrusters; Appendix 3 - Theoretical and experimental studies of hollow cathodes; and Appendix 4 -Theoretical, numerical and experimental studies of pulsed plasma thrusters. Especially notable results include the efficacy of using a solenoidal magnetic field downstream of a plasma thruster to collimate the exhaust flow, the development of a new understanding of applied-field MPD thrusters (based on experimentally-validated results from state-of-the art, numerical simulation) leading to predictions of improved performance, an experimentally-validated, first-principles model for orificed, hollow-cathode behavior, and the first time-dependent, two-dimensional calculations of ablation-fed, pulsed plasma thrusters.

Turchi, P. J.↗
NASA NTRS2012

Prediction of Liquid Slosh Damping Using a High Resolution CFD Tool

Propellant slosh is a potential source of disturbance critical to the stability of space vehicles. The slosh dynamics are typically represented by a mechanical model of a spring mass damper. This mechanical model is then included in the equation of motion of the entire vehicle for Guidance, Navigation and Control analysis. Our previous effort has demonstrated the soundness of a CFD approach in modeling the detailed fluid dynamics of tank slosh and the excellent accuracy in extracting mechanical properties (slosh natural frequency, slosh mass, and slosh mass center coordinates). For a practical partially-filled smooth wall propellant tank with a diameter of 1 meter, the damping ratio is as low as 0.0005 (or 0.05%). To accurately predict this very low damping value is a challenge for any CFD tool, as one must resolve a thin boundary layer near the wall and must minimize numerical damping. This work extends our previous effort to extract this challenging parameter from first principles: slosh damping for smooth wall and for ring baffle. First the experimental data correlated into the industry standard for smooth wall were used as the baseline validation. It is demonstrated that with proper grid resolution, CFD can indeed accurately predict low damping values from smooth walls for different tank sizes. The damping due to ring baffles at different depths from the free surface and for different sizes of tank was then simulated, and fairly good agreement with experimental correlation was observed. The study demonstrates that CFD technology can be applied to the design of future propellant tanks with complex configurations and with smooth walls or multiple baffles, where previous experimental data is not available.

Yang, H. Q.↗
NASA NTRS2019

Risk Analysis Associated with Loss of Toxic Gases During Orion Landing and Recovery Operations

Mission, landing and recovery operations for the Orion crew module involve reentry into the Earth's atmosphere and the deployment of three Nomex parachutes to slow the descent before landing along the west coast of the United States. Orion may have residual fuel (hydrazine, N2H4) or coolant (ammonia, NH3) on board which are both highly toxic to crew in the event of exposure. These risks were evaluated using a first principles analysis approach through fluid dynamics modeling. Plume calculations were first performed with the ANSYS Fluent computational fluid dynamics code. Data were then extracted at locations relevant to crew safety such as the snorkel fan inlet and the egress hatch. Mixing calculations were performed to quantify exposure concentrations within the crew bay before and during egress and departure. Finally, results included herein were used to inform the Orion post-landing Concept of Operations (ConOps) so that strategies could be formulated to maintain crew safety in the event of the loss of fuel or coolant.

Orion↗
DOE OSTI2025

Particle-in-cell simulations of expanding high energy density plasmas with laser ray tracing

The design and analysis of high energy density (HED) laser experiments typically rely on radiation hydrodynamics simulations. However, some laser–plasma interaction regimes are not collisional and cannot be adequately modeled with hydrodynamics. For example, strongly driven magnetic reconnection and magnetized collisionless shock experiments possess extended hydrodynamic or even kinetic properties, necessitating first-principles kinetic simulations. In this paper, we present the benchmarking and first results obtained with a laser-ray-tracing and inverse bremsstrahlung absorption module implemented in the particle-in-cell code PSC. The simulation results are compared to radiation hydrodynamic simulations using the FLASH code as well as analytical estimates. We successfully benchmark the energy deposition model and overall hydrodynamic evolution of the systems. We also consider possible kinetic effects that may be expected from laser-target ablation in the HED regime, including non-local transport and two-temperature effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2013

Using Outcrop Exposures on the Road to Yellowknife Bay to Build a Stratigraphic Column, Gale Crater, Mars

Since landing in Gale Crater on August 5, 2012, the Curiosity rover has driven 450 m east, descending approximately 15 m in elevation from the Bradbury landing site to Yellowknife Bay. Outcrop exposure along this drive has been discontinuous, but isolated outcrops may represent windows into underlying inplace stratigraphy. This study presents an inventory of outcrops targeted by Curiosity (Figs. 1-2), grouped by lithological properties observed in Mastcam and Navcam imagery. Outcrop locations are placed in a stratigraphic context using orbital imagery and first principles of stratigraphy. The stratigraphic models presented here represent an essential first step in understanding the relative age relationships of lithological units encountered at the Curiosity landing site. Such observations will provide crucial context for assessing habitability potential of ancient Gale crater environments and organic matter preservation.

stratigraphic column↗
DOE OSTI2025

Two-Dimensional Ferroelectric Altermagnets: From Model to Material Realization

Multiferroic altermagnets offer new opportunities for magnetoelectric coupling and electrically tunable spintronics. However, due to intrinsic symmetry conflicts between altermagnetism and ferroelectricity, achieving their coexistence, known as ferroelectric altermagnets (FEAM), remains an outstanding challenge, especially in two-dimensional (2D) systems. Here, we propose a universal, symmetry-based design principle for 2D FEAM, supported by tight-binding models and first-principles calculations. We show that lattice distortions can break spin equivalence and introduce the necessary rotation-related symmetry, enabling altermagnetism with electrically reversible spin splitting. Guided by this framework, we identify a family of 2D vanadium oxyhalides and sulfide halides as promising FEAM candidates. In these compounds, pseudo Jahn-Teller distortions and Peierls-like dimerization cooperatively establish the required symmetry conditions. Here, we further propose the magneto-optical Kerr effect as an experimental probe to confirm FEAM and its electric spin reversal. Furthermore, our findings provide a practical framework for 2D FEAM and advancing electrically controlled spintronic devices.

2D materials↗
DOE OSTI2024

Self-consistent interaction of linear gravitational and electromagnetic waves in non-magnetized plasma

Here, this paper explores the hybridization of linear metric perturbations with linear electromagnetic (EM) perturbations in non-magnetized plasma for a general background metric. The local wave properties are derived from first principles for inhomogeneous plasma, without assuming any symmetries of the background metric. First, we derive the effective (“oscillation-center”) Hamiltonian that governs the average dynamics of plasma particles in a prescribed quasimonochromatic wave that involves metric perturbations and EM fields simultaneously. Then, using this Hamiltonian, we derive the backreaction of plasma particles on the wave itself and obtain gauge-invariant equations that describe the resulting self-consistent gravito-electromagnetic (GEM) waves in a plasma. The transverse tensor modes of gravitational waves are found to have no interaction with the plasma and the EM modes in the geometrical-optics limit. However, for $longitudinal$ GEM modes with large values of the refraction index, the interplay between gravitational and EM interactions in plasma can have a strong effect. In particular, the dispersion relation of the Jeans mode is significantly affected by electrostatic interactions. As a spin-off, our calculation also provides an alternative resolution of the so-called Jeans swindle.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Fully Kinetic Shearing-box Simulations of Magnetorotational Turbulence in 2D and 3D. I. Pair Plasmas

The magnetorotational instability (MRI) is a fundamental mechanism determining the macroscopic dynamics of astrophysical accretion disks. In collisionless accretion flows around supermassive black holes, MRI-driven plasma turbulence cascading to microscopic (i.e., kinetic) scales can result in enhanced angular-momentum transport and redistribution, nonthermal particle acceleration, and a two-temperature state where electrons and ions are heated unequally. However, this microscopic physics cannot be captured with standard magnetohydrodynamic (MHD) approaches typically employed to study the MRI. In this work, we explore the nonlinear development of MRI turbulence in a pair plasma, employing fully kinetic particle-in-cell (PIC) simulations in two and three dimensions. First, we thoroughly study the axisymmetric MRI with 2D simulations, explaining how and why the 2D geometry produces results that differ substantially from 3D MHD expectations. We then perform the largest (to date) 3D simulations, for which we employ a novel shearing-box approach, demonstrating that 3D PIC models can reproduce the mesoscale (i.e., MHD) MRI dynamics in sufficiently large runs. With our fully kinetic simulations, we are able to describe the nonthermal particle acceleration and angular-momentum transport driven by the collisionless MRI. Since these microscopic processes ultimately lead to the emission of potentially measurable radiation in accreting plasmas, our work is of prime importance to understand current and future observations from first principles, beyond the limitations imposed by fluid (MHD) models. While in this first study we focus on pair plasmas for simplicity, our results represent an essential step toward designing more realistic electron–ion simulations, on which we will focus in future work.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Exploring the interfacial structure and crystallinity for direct growth of Mn 3 Sn(0001) on sapphire (0001) by molecular beam epitaxy

The Kagome antiferromagnet Mn 3 Sn has garnered a lot of attention due to the presence of interesting properties such as anomalous Hall effects and Nernst effects. Until now, few papers have been reported to grow using molecular beam epitaxy with a buffer layer or on other substrates. In this paper, we discuss the synthesis of crystalline Mn 3 Sn layers, prepared on Al 2 O 3 (0001) without a buffer layer using molecular beam epitaxy. The growth is monitored in-situ using reflection high energy electron diffraction and measured ex-situ using X-ray diffraction, Rutherford back-scattering, and cross-sectional scanning transmission electron microscopy. The samples were deposited at 524 ± 5 °C, with an Mn: Sn atomic flux ratio of 3.2:1 for 90 minutes. Orientation relationships between the Mn 3 Sn films and the sapphire substrates are determined from in-plane and out-of-plane measurements. Our analysis indicates that the resulting film is predominantly c−plane oriented. Lastly, the samples prepared in this way were found to be discontiguous, showing a 3-dimensional morphology. According to first-principles calculations, the Mn 3 Sn exhibits a displaced Kagome structure in the very first stages of growth, for 2ML and 4ML growth on Al 2 O 3 (0001). This result is corroborated by calculating the surface formation energies and explains the observed RHEED patterns.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Symmetry progression and band vs Mott character of CdPS 3 under pressure

Complex chalcogenides are renowned for their tunable electronic, magnetic, and optical properties under external stimuli. The MPX 3 family (M = Mn, Ni, Co, V; X = S, Se) is a platform for many exciting discoveries—especially under compression—although CdPS3 is thought to be different because the Cd center possesses a filled 4d shell, which precludes Mottness. Here, we combine synchrotron-based infrared absorbance and Raman scattering spectroscopies with diamond anvil cell techniques, complementary lattice dynamics calculations, and an analysis of the energy landscape to reveal a series of structural phase transitions in CdPS 3 . We find four distinct pressure-driven transitions, with low frequency modes detectable over the full 35 GPa range of our investigation. A group–subgroup analysis along with our first-principles calculations allows us to partially unravel the space group sequence. For instance, the first critical pressure is a monoclinic C2/m to trigonal $\overline{R}$3 transition at 10 GPa. Despite the softness and overall sensitivity to pressure, we do not locate an insulator-to-metal transition in this pressure range, indicating that the energy scale for gap closure is significantly higher than expected. We discuss these findings in terms of force-induced color change and Mott vs band character in this system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

High-pressure behavior of 3.65 Å phase: Insights from Raman spectroscopy

Abstract The 3.65 Å phase [MgSi(OH)6] is a hydrous phase that is predicted to be stable in a simplified MgO-SiO2-H2O (MSH) ternary system at pressures exceeding 9 GPa. Along cold subduction zones, it is likely to transport water, bound in its crystalline lattice, into the Earth’s interior. The 3.65 Å phase consists of Mg and Si octahedral sites attached to the hydroxyl group that forms a hydrogen bond and is predicted to undergo pressure-induced symmetrization of the hydrogen bond. Therefore, in this study, we investigate the high-pressure behavior of the 3.65 Å phase using Raman spectroscopy. We have conducted five distinct compressions up to ~60 GPa using two different pressure-transmitting media—alcohol mixture and neon. At ambient conditions, we identified vibrational modes using complementary first-principles simulations based on density functional perturbation theory. Upon compression, we note that the first derivative of the vibrational modes in the lattice region stiffens, i.e., b1lattice > 0. In contrast, the hydroxyl region softens, i.e., b1OH > 0. This is indicative of the strengthening of hydrogen bonding upon compression. We noticed a significant broadening of vibrational modes related to hydroxyl groups that are indicative of proton disorder. However, within the maximum pressures explored in this study, we did not find evidence for pressure-induced symmetrization of the hydrogen bonds. We used the pressure derivative of the vibrational modes to determine the ratio of the bulk moduli and their pressure derivative. We note that the smaller bulk moduli of hydrous phases compared to the major mantle phases are compensated by significantly larger pressure derivatives of the bulk moduli for the hydrous phases. This leads to a significant reduction in the elasticity contrast between hydrous and major mantle phases. Consequently, the detection of the degree of mantle hydration is likely to be challenging at greater depths.

Geochemistry & Geophysics↗
NASA NTRS2012

Performance Characterization and Simulation of Amine-Based Vacuum Swing Sorption Units for Spacesuit Carbon Dioxide and Humidity Control

Controlling carbon dioxide (CO2) and water (H2O) vapor concentrations in a space suit is critical to ensuring an astronauts safety, comfort, and capability to perform extra-vehicular activity (EVA) tasks. Historically, this has been accomplished using lithium hydroxide (LiOH) and metal oxide (MetOx) canisters. Lithium hydroxide is a consumable material that requires priming with water before it becomes effective at removing carbon dioxide. MetOx is regenerable through a power-intensive thermal cycle but is significantly heavier on a volume basis than LiOH. As an alternative, amine-based vacuum swing beds are under aggressive development for EVA applications. The vacuum swing units control atmospheric concentrations of both CO2 and H2O through fully-regenerative process. The current concept, referred to as the rapid cycle amine (RCA), has resulted in numerous laboratory prototypes. Performance of these prototypes have been assessed experimentally and documented in previous reports. To support developmental e orts, a first principles model has also been established for the vacuum swing sorption technology. For the first time in several decades, a major re-design of Portable Life Support System (PLSS) for the extra-vehicular mobility unit (EMU) is underway. NASA at Johnson Space Center built and tested an integrated PLSS test bed of all subsystems under a variety of simulated EVA conditions of which the RCA prototype played a significant role. The efforts documented herein summarize RCA test performance and simulation results for single and variable metabolic rate experiments in an integrated context. In addition, a variety of off-nominal tests were performed to assess the capability of the RCA to function under challenging circumstances. Tests included high water production experiments, degraded vacuum regeneration, and deliberate valve/power failure and recovery.

Swickrath, Michael J.↗
DOE OSTI2026

Locating the QCD critical point through contours of constant entropy density

We propose a new method to investigate the existence and location of the conjectured high-temperature critical point of strongly interacting matter via contours of constant entropy density. By approximating these lines as a power series in the baryon chemical potential 𝜇 𝐵 , one can extrapolate them from first-principle results at zero net-baryon density, and use them to locate the quantum chromodynamics (QCD) critical point, including the associated first-order and spinodal lines. As a proof of principle, we employ currently available continuum-extrapolated lattice data from the Wuppertal-Budapest collaboration to find a critical point at a temperature and a baryon chemical potential of 𝑇 𝑐 = 114.3 ± 6.9 MeV and 𝜇 𝐵,𝑐 = 602.1 ± 62.1 MeV, respectively, at expansion order 𝒪⁡(𝜇$^{2}_{𝐵}$). We advocate for a more precise determination of the required expansion coefficients via lattice QCD simulations as a means of pinpointing the location of the critical endpoint in the phase diagram of strongly interacting matter.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Importance of van der Waals interactions for ab initio studies of topological insulators

We investigate the lattice and electronic structures of the bulk and surface of the prototypical layered topological insulators Bi 2 Se 3 and Bi 2 Te 3 using ab initio density functional methods, and systematically compare the results of different methods of including van der Waals (vdW) interactions. We show that the methods utilizing semi-empirical energy corrections yield accurate descriptions of these materials, with the most precise results obtained by properly accounting for the long-range tail of the vdW interactions. The bulk lattice constants, distances between quintuple layers and the Dirac velocity of the topological surface states (TSS) are all in excellent agreement with experiment. In Bi 2 Te 3 , hexagonal warping of the energy dispersion leads to complex spin textures of the TSS at moderate energies, while in Bi 2 Se 3 these states remain almost perfectly helical away from the Dirac point, showing appreciable signs of hexagonal warping at much higher energies, above the minimum of the bulk conduction band. Furthermore, our results establish a framework for unified and systematic self-consistent first principles calculations of topological insulators in bulk, slab and interface geometries, and provides the necessary first step toward ab initio modeling of topological heterostructures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Magnetic contributions to corundum-eskolaite and corundum-hematite phase equilibria: A DFT cluster expansion study

Magnetic contributions have the potential to significantly influence predicted phase stability within alloy and mineral mixing phase diagrams, yet have been historically challenging to incorporate due to a significant increase to phase space sampling. Here, in this work, we employ a computational protocol that includes spin orientation as an additional configurational component within multi-component cluster expansions between magnetic and non-magnetic metal oxide alloys [calculated using density functional theory (DFT) and the generalized gradient approximation]. This approach was used to determine the effect of magnetic contributions to corundum-eskolaite and corundum-hematite phase equilibria from first principles. Two-component cluster expansions of the magnetic components of eskolaite and hematite were first performed showing the ability of this method to properly calculate their respective magnetic properties. Two-component cluster expansions were then performed for non-magnetic Al(III) and ferromagnetic Cr(III) and Fe(III), and phase diagrams were calculated for later comparison. Finally, a non-magnetic Al(III) and “up” and “down” magnetic configurations for anti-ferromagnetic Cr(III) and Fe(III) were performed. Magnetic contributions to the calculated phase diagram for the corundumeskolaite system were shown to be inconsequential but are vital for accurate determination of the corundum-hematite solvus.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1989

Automatic determination of fault effects on aircraft functionality

The problem of determining the behavior of physical systems subsequent to the occurrence of malfunctions is discussed. It is established that while it was reasonable to assume that the most important fault behavior modes of primitive components and simple subsystems could be known and predicted, interactions within composite systems reached levels of complexity that precluded the use of traditional rule-based expert system techniques. Reasoning from first principles, i.e., on the basis of causal models of the physical system, was required. The first question that arises is, of course, how the causal information required for such reasoning should be represented. The bond graphs presented here occupy a position intermediate between qualitative and quantitative models, allowing the automatic derivation of Kuipers-like qualitative constraint models as well as state equations. Their most salient feature, however, is that entities corresponding to components and interactions in the physical system are explicitly represented in the bond graph model, thus permitting systematic model updates to reflect malfunctions. Researchers show how this is done, as well as presenting a number of techniques for obtaining qualitative information from the state equations derivable from bond graph models. One insight is the fact that one of the most important advantages of the bond graph ontology is the highly systematic approach to model construction it imposes on the modeler, who is forced to classify the relevant physical entities into a small number of categories, and to look for two highly specific types of interactions among them. The systematic nature of bond graph model construction facilitates the process to the point where the guidelines are sufficiently specific to be followed by modelers who are not domain experts. As a result, models of a given system constructed by different modelers will have extensive similarities. Researchers conclude by pointing out that the ease of updating bond graph models to reflect malfunctions is a manifestation of the systematic nature of bond graph construction, and the regularity of the relationship between bond graph models and physical reality.

Feyock, Stefan↗
DOE OSTI2024

Theoretical Investigation of the Adsorbate and Potential–Induced Stability of Cu Facets During Electrochemical CO 2 and CO Reduction

The activity and product selectivity of electrocatalysts for reactions like the carbon dioxide reduction reaction (CO 2 RR) are intimately dependent on the catalyst's structure and composition. While engineering catalytic surfaces can improve performance, discovering the key sets of rational design principles remains challenging due to limitations in modeling catalyst stability under operating conditions. Herein, we perform first-principles density functional calculations adopting implicit solvation methods with potential control to study the influence of adsorbates and applied potential on the stability of different facets of model Cu electrocatalysts. Using coverage dependencies extracted from microkinetic models, we describe an approach for calculating potential and adsorbate-dependent contributions to surface energies under reaction conditions, where Wulff constructions are used to understand the morphological evolution of Cu electrocatalysts under CO 2 RR conditions. Here we identify that CO*, a key reaction intermediate, exhibits higher kinetically and thermodynamically accessible coverages on (100) relative to (111) facets, which can translate into an increased relative stabilization of the (100) facet during CO 2 RR. Our results support the known tendency for increased (111) faceting of Cu nanoparticles under more reducing conditions and that the relative increase in (100) faceting observed under CO 2 RR conditions is likely attributed to differences in CO* coverage between these facets.

30 DIRECT ENERGY CONVERSION↗
NASA NTRS1985

Physical concepts in the development of constitutive equations

Proposed viscoplastic material models include in their formulation observed material response but do not generally incorporate principles from thermodynamics, statistical mechanics, and quantum mechanics. Numerous hypotheses were made for material response based on first principles. Many of these hypotheses were tested experimentally. The proposed viscoplastic theories and the experimental basis of these hypotheses must be checked against the hypotheses. The physics of thermodynamics, statistical mechanics and quantum mechanics, and the effects of defects, are reviewed for their application to the development of constitutive laws.

Cassenti, B. N.↗