SEARCH · Engineering Papers

Results for “Matrix completion”

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 811 records · Page 45

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Rapid High-Resolution Analysis of Polysaccharide-Lignin Interactions in Secondary Plant Cell Walls Using Proton-Detected Solid-State NMR

The plant secondary cell wall, a complex matrix composed of cellulose, hemicellulose, and lignin, is crucial for the mechanical strength and water-proofing properties of plant tissues, and serves as a primary source of biomass for biorenewable energy and biomaterials. Structural analysis of these polymers and their interactions within the secondary cell wall has been heavily relying on 13 C-based solid-state NMR techniques. In this study, we explore the application of 1 H-detected solid-state NMR techniques for rapid, high-resolution structural characterization of polysaccharides and lignin, demonstrated on the stems of hardwood eucalyptus. We explored the use of synthesized 2D spectra to resolve central 1 H resonances and the combined application of 3D hCCH and hCHH experiments for complete resonance assignment and unambiguous identification of lignin-carbohydrate interactions. Our findings emphasize the central role of acetylated three-fold xylan conformers, rather than two-fold, in stabilizing the carbohydrate-lignin interface, with glucuronic acid sidechains in eucalyptus glucuronoxylan colocalizing with lignin, revised cellulose-lignin interactions involving uncoated microfibril surfaces, and pectin-lignin interactions indicative of early-stage lignification. These results present a novel approach for rapid structural analysis of lignocellulosic biomaterials without the need for solubilization or extraction.

09 BIOMASS FUELS↗

DOE OSTI2023

Processing Temperature Effect on Failure Behavior of Unidirectional Thermoplastic Polymer-fiber-reinforced Polymers (PFRPs)

This study aimed to examine the influence of processing temperature on the mechanical and morphological characteristics of unidirectional (UD) thermoplastic polymer-fiber-reinforced polymers (PFRPs), utilizing UHMWPE-polymer-fiber-reinforced HDPE composites as an example. To achieve this objective, UD thermoplastic PFRPs were produced through filament winding and hot pressing techniques, employing five distinct processing temperatures within the range spanning from the melting onset temperature of the polymer matrix to the melting peak temperature of the reinforcing polymer fiber. By conducting the uni-axial tensile tests on the UD UHMWPE-polymer-fiber-reinforced HDPE composites fabricated at different processing temperatures, it was shown that the optimal longitudinal tensile properties of the composites, comparable with glass-fiber-reinforced thermoset or thermoplastic polymers (GFRPs or GFRTPs) can be achieved when the processing temperature does not significantly exceeds the melting onset temperature of the reinforcing polymer fiber. However, when the processing temperature exceeds this threshold, the mechanical properties of the composites are significantly reduced, as evidenced by the transition in the failure morphology from the presence of significant splitting cracks to the occurrence of plastic necking due to the complete melting of the reinforcing polymer fibers in a composite. Particularly, the ductile behavior of the optimal thermoplastic PFRPs investigated in this study is superior than that of carbon-fiber-reinforced thermoset or thermoplastic polymers (CFRPs or CFRTPs) and even GFRPs and GFRTPs. This study not only provides valuable insights into the proper fabrication of high-performance thermoplastic PFRPs but also offers useful experimental data that can aid in the validation and development of computational models, particularly those related to processing modeling.

Qiao, Yao↗

DOE OSTI2024

Transverse momentum-dependent heavy-quark fragmentation at next-to-leading order

The transverse momentum-dependent fragmentation functions (TMD FFs) of heavy (bottom and charm) quarks, which we recently introduced, are universal building blocks that enter predictions for a large number of observables involving final-state heavy quarks or hadrons. They enable the extension of fixed-order subtraction schemes to quasi-collinear limits, and are of particular interest in their own right as probes of the nonperturbative dynamics of hadronization. In this paper we calculate all TMD FFs involving heavy quarks and the associated TMD matrix element in heavy-quark effective theory (HQET) to next-to-leading order in the strong interaction. Our results confirm the renormalization properties, large-mass, and small-mass consistency relations predicted in our earlier work. We also derive and confirm a prediction for the large-z behavior of the heavy-quark TMD FF by extending, for the first time, the formalism of joint resummation to capture quark mass effects in heavy-quark fragmentation. Our final results in position space agree with those of a recent calculation by another group that used a highly orthogonal organization of singularities in the intermediate momentum-space steps, providing a strong independent cross check. As an immediate application, we present the complete quark mass dependence of the energy-energy correlator (EEC) in the back-to-back limit at $\mathcal{O}\left({\alpha}_s\right)$.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Power and isotope effects in the ITER baseline scenario with tungsten and tungsten-equivalent radiators in DIII-D

Abstract Experiments in DIII-D document the ITER Baseline Scenario (IBS) at q 95 ∼ 3 and P IN / P LH ∼ 1–2, in both deuterium and hydrogen utilizing Kr and Xe as Tungsten-equivalent radiators. The power threshold for H-mode operation ( P LH ) was determined experimentally without added impurities and found to be about a factor of two higher than the scaling law. In recent IBS experiments in deuterium, intrinsic levels of metals such as Tungsten (W) or molybdenum and inconel are present that reduce the pedestal pressure by 20%–25%. A complete radiative collapse of deuterium IBS plasmas occurs at W core concentrations C W = 10 −5 . Simulations show that for core temperatures expected for ITER, the plasmas would not have a radiative collapse at C W = 1 × 10 −5 , moreover Q = 8–10 would still be achieved for C W up to 3 × 10 −5 . In contrast to deuterium, the IBS in hydrogen is not affected by intrinsic high-Z impurities, indicating that hydrogen H-modes in ITER may not inform the D-T phase with respect to W accumulation and discharge survival. Compared to deuterium, the pedestal pressure in hydrogen is ∼25% lower, with much higher ELM frequency of 150 Hz, decreasing with input power. Krypton was injected in a matrix scan of input power and impurity flow in IBS hydrogen discharges. Krypton impurity density profiles in hydrogen are similar to deuterium plasmas, but at Kr flows that are 2–3 times higher for the same input power. Krypton is transported into the core and affects the whole radius; at the highest injection rates a radiative collapse occurs at core radiation fractions of 0.3–0.35, consistent with the expected maximum W radiation fraction for ITER core plasmas. Comparing the results with previous International Tokamak Physics Activity database studies of the IBS confirms that at higher radiation fraction due to high-Z impurities, a drop in H 98 of >10% is observed. On the other hand, the results using Kr as a W-equivalent radiator indicate that metal (W) devices at lower core temperatures than ITER may provide overly pessimistic performance extrapolations to ITER for deuterium-tritium operation. The new DIII-D results support a more attractive option for the ITER Research Plan with a short hydrogen phase for system commissioning, transitioning to deuterium operations as soon as possible to provide relevant conditions for deuterium-tritium operations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Assessment of the Phase-to-Ground Fault Apparent Admittance Method with Phase/Ground Boundaries to Detect Types of Electrical Faults for Protective Relays Using Signature Library and Simulated Events

Protective relays in electric power grids recognize the types of electrical faults in a few seconds. The most common detection method to detect the types of electrical faults is based on measuring the angle between the zero and negative sequence currents. However, it is not completely accurate because the phase-to-phase-ground and phase-to-ground electrical faults could have the same detection conditions. Therefore, engineers need to plot the events after an electrical fault to observe the nature of the incidents in detail. In this study, the phase-to-ground fault apparent (PGFA) admittance method with phase/ground boundaries identified the types of electrical faults located in distribution power lines and feeders. This method was based on measuring the PGFA admittance magnitudes for the faulted and nonfaulted phases, resulting in greater than zero and near zero, respectively. The PGFA admittance algorithm was built with MATLAB/Simulink software and tested with signature library and grid simulation events. The PGFA method with phase/ground boundaries was evaluated with the confusion matrix. The measured and predicted values matched in more than 90% of the tests, and the PGFA admittance method with phase/ground boundaries presented an accuracy of 94.3% and a precision of 100%.

42 ENGINEERING↗

NASA NTRS1995

Robust integration schemes for generalized viscoplasticity with internal-state variables. Part 1: Theoretical developments and applications

This two-part report is concerned with the development of a general framework for the implicit time-stepping integrators for the flow and evolution equations in generalized viscoplastic models. The primary goal is to present a complete theoretical formulation, and to address in detail the algorithmic and numerical analysis aspects involved in its finite element implementation, as well as to critically assess the numerical performance of the developed schemes in a comprehensive set of test cases. On the theoretical side, the general framework is developed on the basis of the unconditionally-stable, backward-Euler difference scheme as a starting point. Its mathematical structure is of sufficient generality to allow a unified treatment of different classes of viscoplastic models with internal variables. In particular, two specific models of this type, which are representative of the present start-of-art in metal viscoplasticity, are considered in applications reported here; i.e., fully associative (GVIPS) and non-associative (NAV) models. The matrix forms developed for both these models are directly applicable for both initially isotropic and anisotropic materials, in general (three-dimensional) situations as well as subspace applications (i.e., plane stress/strain, axisymmetric, generalized plane stress in shells). On the computational side, issues related to efficiency and robustness are emphasized in developing the (local) interative algorithm. In particular, closed-form expressions for residual vectors and (consistent) material tangent stiffness arrays are given explicitly for both GVIPS and NAV models, with their maximum sizes 'optimized' to depend only on the number of independent stress components (but independent of the number of viscoplastic internal state parameters). Significant robustness of the local iterative solution is provided by complementing the basic Newton-Raphson scheme with a line-search strategy for convergence. In the present first part of the report, we focus on the theoretical developments, and discussions of the results of numerical-performance studies using the integration schemes for GVIPS and NAV models.

Saleeb, Atef F.↗

NASA NTRS1995

Robust integration schemes for generalized viscoplasticity with internal-state variables. Part 2: Algorithmic developments and implementation

This two-part report is concerned with the development of a general framework for the implicit time-stepping integrators for the flow and evolution equations in generalized viscoplastic models. The primary goal is to present a complete theoretical formulation, and to address in detail the algorithmic and numerical analysis aspects involved in its finite element implementation, as well as to critically assess the numerical performance of the developed schemes in a comprehensive set of test cases. On the theoretical side, the general framework is developed on the basis of the unconditionally-stable, backward-Euler difference scheme as a starting point. Its mathematical structure is of sufficient generality to allow a unified treatment of different classes of viscoplastic models with internal variables. In particular, two specific models of this type, which are representative of the present start-of-art in metal viscoplasticity, are considered in applications reported here; i.e., fully associative (GVIPS) and non-associative (NAV) models. The matrix forms developed for both these models are directly applicable for both initially isotropic and anisotropic materials, in general (three-dimensional) situations as well as subspace applications (i.e., plane stress/strain, axisymmetric, generalized plane stress in shells). On the computational side, issues related to efficiency and robustness are emphasized in developing the (local) interative algorithm. In particular, closed-form expressions for residual vectors and (consistent) material tangent stiffness arrays are given explicitly for both GVIPS and NAV models, with their maximum sizes 'optimized' to depend only on the number of independent stress components (but independent of the number of viscoplastic internal state parameters). Significant robustness of the local iterative solution is provided by complementing the basic Newton-Raphson scheme with a line-search strategy for convergence. In the present second part of the report, we focus on the specific details of the numerical schemes, and associated computer algorithms, for the finite-element implementation of GVIPS and NAV models.

Li, Wei↗

DOE OSTI2023

Assessment of the Phase-to-Ground Fault Apparent Admittance Method with Phase/Ground Boundaries to Detect Types of Electrical Faults for Protective Relays Using Signature Library and Simulated Events

Protective relays in electric power grids recognize the types of electrical faults in a few seconds. The most common detection method to detect the types of electrical faults is based on measuring the angle between the zero and negative sequence currents. However, it is not completely accurate because the phase to phase ground and phase to ground electrical faults could have the same detection conditions. Therefore, engineers need to plot the events after an electrical fault to observe the nature of the incidents in detail. In this study, the phase to ground fault apparent (PGFA) admittance method with phase/ground boundaries identified the types of electrical faults located in distribution power lines and feeders. This method was based on measuring the PGFA admittance magnitudes for the faulted and non-faulted phases, resulting in greater than zero and near zero, respectively. The PGFA admittance algorithm was built with MATLAB/ Simulink software and tested with signature library and grid simulation events. The PGFA method with phase/ground boundaries was evaluated with the confusion matrix. The measured and predicted values matched in more than 90% of the tests, and the PGFA admittance method with phase/ground boundaries presented an accuracy of 94.3% and a precision of 100%.

Piesciorovsky, Emilio↗

DOE OSTI2025

Nonequilibrium steady-state thermoelectrics of Kondo-correlated quantum dots

The transport across a Kondo-correlated quantum dot coupled to two leads with independent temperatures and chemical potentials is studied using a controlled nonperturbative, and in this sense numerically exact, treatment based on a hybrid numerical renormalization group combined with time-dependent density matrix renormalization group (NRG-tDMRG). In the Kondo regime, for sufficiently large fixed voltage bias V ≳ T K , with T K the Kondo temperature, we find a peak in the conductance vs the temperature gradient Δ⁢T = T R - T L across left and right lead. Focusing then on zero voltage bias but finite ΔT far beyond linear response, we reveal the dependence of the characteristic zero-bias conductance on the individual lead temperatures. Here, we find that the finite-Δ⁢T data behaves quantitatively similar to linear response with an effective equilibrium temperature derived from the different lead temperatures. The regime of sign changes in the Seebeck coefficient, signaling the presence of Kondo correlations, and its dependence on the individual lead temperatures provide a complete picture of the Kondo regime in the presence of finite-temperature gradients. The results from the zero-bias conductance and Seebeck coefficient studies unveil an approximate “Kondo circle” in the T L /T R plane as the regime within which the Kondo correlations dominate. We also study the heat current and the corresponding heat conductance vs finite Δ⁢T. We provide a polynomial fit for our numerical results for the thermocurrent as a function of the individual lead temperatures, which may be used to fit experimental data in the Kondo regime.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Exploring the interference between the atmospheric and solar neutrino oscillation subamplitudes

The interference between the atmospheric and solar neutrino oscillation subamplitudes is said to be responsible for C P violation () in neutrino appearance channels. More precisely, is generated by the interference between the parts of the neutrino oscillation amplitude that are C P even and C P odd: even or odd when the neutrino mixing matrix is replaced with its complex conjugate. This is the interference term, as it gives a contribution to the oscillation probability, the square of the amplitude, which is opposite in sign for neutrinos and antineutrinos and is unique. For this interference to be nonzero, at least two subamplitudes are required. There are, however, other interference terms, which are even under the above exchange, and these are the C P conserving () interference terms. In this paper, we explore in detail these interference terms and show that they cannot be uniquely defined, as one can move pieces of the amplitude from the atmospheric subamplitude to the solar subamplitude and vice versa. This freedom allows one to move the interference terms around, but does not let you eliminate them completely. We also show that there is a reasonable definition of the atmospheric and solar subamplitudes for the appearance channels such that in neutrino disappearance probability there is no atmospheric-solar interference term. However, with this choice, there is a interference term within the atmospheric sector. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS2011

Supersonic Retropropulsion Experimental Results from the NASA Langley Unitary Plan Wind Tunnel

A new supersonic retropropulsion experimental effort, intended to provide code validation data, was recently completed in the Langley Research Center Unitary Plan Wind Tunnel Test Section 2 over the Mach number range from 2.4 to 4.6. The experimental model was designed using insights gained from pre-test computations, which were instrumental for sizing and refining the model to minimize tunnel wall interference and internal flow separation concerns. A 5-in diameter 70-deg sphere-cone forebody with a roughly 10-in long cylindrical aftbody was the baseline configuration selected for this study. The forebody was designed to accommodate up to four 4:1 area ratio supersonic nozzles. Primary measurements for this model were a large number of surface pressures on the forebody and aftbody. Supplemental data included high-speed Schlieren video and internal pressures and temperatures. The run matrix was developed to allow for the quantification of various sources of experimental uncertainty, such as random errors due to run-to-run variations and bias errors due to flow field or model misalignments. Preliminary results and observations from the test are presented, while detailed data and uncertainty analyses are ongoing.

Berry, Scott A.↗

NASA NTRS1980

Completion of evaluation of manufacturing processes for B/Al composites containing 0.2mm diameter boron fibers

Four fabricators produced a total of 54 B/1100 Al, B/6061 Al, and B/2024 Al panels for evaluation. The 8 ply unidirectional, 45 to 50 volume percent, panels were made using 0.20 mm diameter boron fibers which were obtained from a single supplier. Hot press consolidation was carried out in vacuum except for one set of dry woven tape panels which were hot pressed in air. A single testing contractor conducted nondestructive inspection, metallography, fractography and mechanical property tests. The mechanical property tests included 21 and 260 C tensile tests and 21 C shear tests. Panel quality, as measured by nondestructive evaluation, was generally good as were the 21 C tensile properties. The panels hot pressed in air delaminated in the shear tests. Shear strength values were lower in these panels. But tensile strengths were not affected by the delaminations because of the relation between the tensile loading direction and the delaminations. Composite tensile strength was found to be proportional to the volume percent boron and the aluminum matrix rather than to the tape used or fabrication technique. Suitability of these composites for 260 C service was confirmed by tensile tests.

Moore, T. J.↗

NASA NTRS2003

High-Time-Resolution Study of Magnetic Holes in the Solar Wind

The objectives of this investigation are to determine the internal plasma structure of kinetic-scale and larger scale magnetic holes, and to determine their stability, their source mechanism(s), and their spatial extent. It is also of importance to determine the relationship between kinetic-scale holes and long-duration holes. As smaller and smaller magnetic depressions are investigated in order to make this a complete study, a robust criterion is necessary for distinguishing magnetic holes from random or unresolvable fluctuations in the interplanetary magnetic field. In order to resolve this ambiguity, we obtained from the MFI experiments magnetic field measurements from the WIND spacecraft at a time resolution of 46 to 184 ms over certain periods. We have also devised a measure of certainty for magnetic hole detections. The certainty factor, q, is defined as the difference between the mean magnetic field in the hole and the local magnetic field, in units of the local standard deviation of the field strength. For fullest generality, it is necessary to calculate this q over the range of available scales of interest, from 60 ms up to 300 s. This technique establishes a two dimensional matrix of relative probabilities that a hole of some duration (d) might exist in the data set at a given time (t). In identifying q-peaks in time and duration, we also come upon a natural method for distinguishing holes with internal structure from multiple holes in close proximity or holes nested inside of others. If two q-peaks are more than a half-width apart, they are simply said to be separate events.

Lazarus, Alan↗

NASA NTRS2018

The Second European Service Module (ESM-2) Evolutions, Production and Challenges

This paper presents an overview of the Second European Service Module (ESM-2), the second in a series of European Service Modules produced as part of the Barter agreement between NASA and ESA for the Orion Program. The European Industrial consortium is led by the ESA prime contractor Airbus Defence and Space in Bremen. ESA and Airbus signed the ESM-2 contract on 16 February 2017, for this key element of the Orion Exploration Mission 2 (EM-2). EM-2 is the first crewed mission for Orion and will take astronauts farther into the solar system than humanity has ever travelled. EM-2 will also be a historic mission for Europe, as the ESM-2 will be the first European spacecraft to be part of a human transportation system carrying humans beyond low Earth orbit. ESM-2 is mainly a recurring production following ESM-1. Nevertheless, there are a number of important changes being implemented, for example, to incorporate upgrades to further enhance safety and reliability. The challenging delivery schedule for ESM-2 has driven the need to commence manufacturing prior to completion of the qualification on ESM-1. In addition, some requirement deviations and non-compliances approved for ESM-1 have resulted in modifications for ESM-2. In order to manage the competing constraints effectively, the ESM-2 Team has put in place a number of novel approaches to manage schedule, risk, and technical changes. Airbus has set up multi-functional teams according to an approach known as "Major Spacecraft Deliveries" consisting of quality assurance, engineering and procurement. The risk of starting manufacturing prior to qualification is managed through a special risk share agreement. This agreement necessitates rigorous risk reviews across the board for all manufacturing, assembly, integration and test milestones. The ESM-2 changes are managed by Configuration Management, but Airbus has also introduced the Technical Baseline Matrix to provide a transparent top-level overview of the changes from ESM-1 to ESM-2. The tool provides the basis for ESM-2 design and development needs, decisions, as well as the input for the Orion EM-2 Critical Design Review (CDR). The main technical evolutions, status of the production and the novel management approaches for ESM-2 are presented and discussed in the paper.

Orion European Service Module Programme↗

NASA NTRS2007

Comparison of Fluka-2006 Monte Carlo Simulation and Flight Data for the ATIC Detector

We have performed a detailed Monte Carlo (MC) simulation for the Advanced Thin Ionization Calorimeter (ATIC) detector using the MC code FLUKA-2006 which is capable of simulating particles up to 10 PeV. The ATIC detector has completed two successful balloon flights from McMurdo, Antarctica lasting a total of more than 35 days. ATIC is designed as a multiple, long duration balloon flight, investigation of the cosmic ray spectra from below 50 GeV to near 100 TeV total energy; using a fully active Bismuth Germanate(BGO) calorimeter. It is equipped with a large mosaic of.silicon detector pixels capable of charge identification, and, for particle tracking, three projective layers of x-y scintillator hodoscopes, located above, in the middle and below a 0.75 nuclear interaction length graphite target. Our simulations are part of an analysis package of both nuclear (A) and energy dependences for different nuclei interacting in the ATIC detector. The MC simulates the response of different components of the detector such as the Si-matrix, the scintillator hodoscopes and the BGO calorimeter to various nuclei. We present comparisons of the FLUKA-2006 MC calculations with GEANT calculations and with the ATIC CERN data and ATIC flight data.

Gunasingha, R.M.↗

DOE OSTI2024

Extracting Higher Central Charge from a Single Wave Function

A (2+1)D topologically ordered phase may or may not have a gappable edge, even if its chiral central charge c_ is vanishing. Recently, it was discovered that a quantity regarded as a "higher" version of chiral central charge gives a further obstruction beyond c_ to gapping out the edge. Here, in this Letter, we show that the higher central charges can be characterized by the expectation value of the partial rotation operator acting on the wave function of the topologically ordered state. This allows us to extract the higher central charge from a single wave function, which can be evaluated on a quantum computer. Our characterization of the higher central charge is analytically derived from the modular properties of edge conformal field theory, as well as the numerical results with the ν = 1/2 bosonic Laughlin state and the non-Abelian gapped phase of the Kitaev honeycomb model, which corresponds to U(1) 2 and Ising topological order, respectively. The Letter establishes a numerical method to obtain a set of obstructions to the gappable edge of (2+1)D bosonic topological order beyond c_, which enables us to completely determine if a (2+1)D bosonic Abelian topological order has a gappable edge or not. We also point out that the expectation values of the partial rotation on a single wave function put a constraint on the low-energy spectrum of the bulk-boundary system of (2+1)D bosonic topological order, reminiscent of the Lieb-Schultz-Mattis-type theorems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Comparing models for a unitary black hole S matrix

This paper compares features, challenges, and puzzles of different models for a unitary black hole S -matrix, focusing on both recent nonisometric models, as well as “nonviolent unitarization,” which is based on new quantum interactions of a black hole. As a foundation for comparison, the description of real-time Hawking evolution is first overviewed, including leading effects of gravitational dressing and backreaction. Connection is then made to qubit models for evolution, and some technology is outlined to facilitate their description. Important features of both nonisometric models and nonviolent unitarization are investigated in qubit models, which illustrate essential differences between the respective approaches. The nonisometric models present puzzles for understanding evolution of internal outgoing excitations, which can be excited by interactions such as particle decay. Qubit models for nonviolent unitarization are further developed, and nicely illustrate aspects of that approach. Some remaining questions in generalizing to more complete models for evolution are discussed. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

NASA NTRS2023

Extended Kalman Filter Performance on the Artemis-1 Mission

The Artemis Program is NASA’s campaign to explore the Moon and beyond. Artemis-1, the uncrewed exoLEO test flight of the Orion spacecraft, was completed in 2022. There were four navigation Extended Kalman Filters (EKFs) that are part of the Orion navigation system. The Atmospheric Extended Kalman Filter (ATMEKF) estimates the vehicle position, velocity, and attitude (referred to as the vehicle state) during the ascent and entry phases of flight. Once Orion is outside of Earths atmosphere, the Earth Orbit Extended Kalman Filter (EOEKF) and CisLunar Extended Kalman Filter (CLEKF) estimate the translational states, depending on the phase of flight, while the Attitude Extended Kalman Filter (ATTEKF) estimates the rotational state of the vehicle. The Kalman filters propagate the vehicle state forward in time using a combination of dynamics models and the output data from the Inertial Measurement Unit (IMU). The filters update the vehicle states and associated uncertainties, in the form of the covariance matrix, using pseudorange measurements from GPS (in ATMEKF/EOEKF), optical navigation measurements of the Earth or Moon (in CLEKF), and star tracker measurements (in ATTEKF). Simultaneously, the Kalman filters estimate error sources in the sensors, which are included in the state vectors as Exponentially Correlated Random Variables (ECRVs). This paper will summarize the performance of these filters during the Artemis-1 mission.

Artemis-1↗