Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Boundary identification”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Materials and Approaches for Mitigation of SOFC Cathode Degradation in SOFC Power Systems

Perovskites such as doped lanthanum manganite and cobaltite have been used successfully as a cathode material in SOFC power systems operating between 700-1000°C. The electrodes, however, remain prone to degradation in presence of Cr containing gas atmosphere under systems operating conditions. Cathode in SOFCs is considered as one of the largest contributors of electrochemical performance degradations, leading to significant increase in ohmic and non-ohmic losses as well as interfacial compound formation. These degradations are mostly irreversible due to permanent changes to chemical and morphological alterations in the electrodes, resulting in the blockage of triple phase boundary (TPB) sites and decrease in oxygen reduction reactions (ORR). Fundamental mechanistic understanding has been developed to elucidate the reaction pathways of electrode poisoning in electronically conducting electrodes (LSM) as well as mixed ionic and electronically conducting electrodes (LSCF). The overall goal of this program has been identification and development of cost effective solutions to mitigate electrode degradations under “real world” systems operating conditions. The overall technological innovations and impact from this work are: The research provides an understanding of the degradation mechanism and cost-effective approaches for implementation in SOFC systems; Develop mitigation processes utilizing low cost getters using conventional synthesis and fabrication route; Develop tailored high surface area powders and coatings for high Cr capture capacity under systems operating conditions; >Developed computational models to understand mechanistic behavior; Experimentally measure rate of Cr evaporation from chromia and alumina forming alloys under SOFC operating conditions; Developed conditions to minimize Cr evaporation from balance-of-plant alloys; This innovation will support the development of getter platform to capture both intrinsic and extrinsic impurities present in air; and, Applications under wide operating temperature range (500-900°C) as well for wide high-temperature electrochemical systems (SOFC, SOEC, and OTM). The work has been successfully developed, validated experimentally and implemented. Successful technology transfer of this knowledge with industrial partners has been accomplished and the getters have been independently tested and validated under their systems operating conditions. The innovation will also find application in related high temperature electrochemical systems such as OTM and SOEC for the prevention of Cr assisted performance degradation. The proposed approach for Cr capture can also be applied to oxy-combustion and other advanced combustion techniques for the reduction of Cr vapor in the exhaust gas stream. This fundamental knowledge gained also has led to the path forward for developing advanced getters for capturing more than one contaminant in air.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2025

Identification of a window for quiescent H-mode operation in MHD stability diagram of DIII-D plasmas

A window for quiescent H-mode (QH-mode) operation in DIII-D plasmas in the diagram of pedestal magnetohydrodynamic (MHD) stability was identified through linear MHD stability analysis, accounting for the effects of plasma rotation and ion diamagnetic drift. The operation window lies between the stability boundaries of the kink/peeling mode (K/PM) identified with and without the inclusion of plasma rotation effects. In this region, the mode remains unstable unless rotation effects are considered alongside the ion diamagnetic drift, which is consistently included in the analysis. The stabilization of the MHD mode, facilitated by the coupled effects of plasma rotation and ion diamagnetic drift, plays a crucial role in widening the window, enabling the attainment of the QH-mode state. Furthermore, the suppression of edge localized modes (ELM) can be achieved by controlling the pedestal structure to ensure the plasma state remains within the operation window. The location of the operation window in the stability diagram depends on the K/PM stability properties. Therefore, optimizing conditions for QH-mode requires adjustments based on changes in stability characteristics. A pressure pedestal and its associated bootstrap current density near the last closed flux surface are advantageous for situating the plasma state within the window. However, excessive current density can trigger ELMs. This trend was confirmed through comparisons of MHD stability diagrams between QH-mode and ELMy H-mode plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Automated extraction of interfacial dislocations and disconnections from atomistic data

Here we introduce Interfacial Line Defect Analysis (ILDA), a method for identifying and extracting interfacial dislocations and disconnections with little-to-no input from the user on the nature of the interface. By simply providing the orientations and coherency strains for the crystals in a coherent reference state of the interface, ILDA provides exact Burgers vectors. Alternatively, these orientations and strains can be estimated using local atomic deformation gradients, making ILDA fully automated and providing estimated Burgers vectors. ILDA also determines the step height associated with each defect line segment, in case the associated defect is a disconnection. The heart of the method is the identification of atoms residing at co-incidence sites between the two crystals and the construction of a surface mesh connecting these sites that is used to compose Burgers circuits and insert defect line segments. We demonstrate the performance of ILDA in two test cases: a twist grain boundary and a phase boundary.

36 MATERIALS SCIENCE↗
NASA NTRS1975

Delineation of the boundaries of a buried pre-glacial valley with LANDSAT-1 data

The continuity of a narrow meandering strip of Udoll (prairie) soils running east and west for approximately 40 miles across north central Indiana in an area predominantly of Udalfs (timber soils) was detected from LANDSAT-1 data taken on June 9, 1973. This data was processed through a clustering procedure and classified with resulting increased definition of the boundaries among soils grouped according to nine categories and vegetation to two categories of reflectance. This dark stretch of prairie soil is believed to have formed in the heavy textured, poorly drained glacial debris which filled a major pre-glacial tributary of the Teays River System. Ready identification and location of the valley has significance to soil survey and land classification people as a guide to soil classification and land use and to geologists as a guide to location of a potentially economically significant aquifer.

Peterson, J. B.↗
NASA NTRS1989

AVHRR for monitoring global tropical deforestation

Advanced Very High Resolution Radiometer (AVHRR) data have been used to assess the dynamics of forest trnsformations in three parts of the tropical belt. A large portion of the Amazon Basin has been systematically covered by Local Area Coverage (LAC) data in the 1985-1987 period. The analysis of the vegetation index and thermal data led to the identification and measurement of large areas of active deforestation. The Kalimantan/Borneo forest fires were monitored and their impact was evaluated using the Global Area Coverage (GAC) 4 km resolution data. Finally, High Resolution Picture Transmission (HRPT) data have provided preliminary information on current activities taking place at the boundary between the savanna and the forest in the Southern part of West Africa. The AVHRR approach is found to be a highly valuable means for carrying out deforestation assessments in regional and global perspectives.

Malingreau, J. P.↗
NASA NTRS2024

Critical Load Exceedances for North America and Europe using an Ensemble of Models and an Investigation of Causes for Environmental Impact Estimate Variability: An AQMEII4 Study

Exceedances of critical loads for deposition of sulphur (S) and nitrogen (N) to different ecosystems were estimated using European and North American ensembles of air quality models, under Phase 4 of the Air Quality Model Evaluation International Initiative (AQMEII4), to identify where risk of ecosystem harm is expected to occur based on model deposition estimates. The ensembles were driven by common emissions and lateral boundary condition inputs. Model output was regridded to common North American and Europe 0.125° resolution domains, which were then used to calculate critical load exceedances. New, targeted deposition diagnostics implemented in AQMEII4 allowed an unprecedented level of post-simulation analysis to be carried out and facilitated the identification of specific causes of model-to-model variability in critical load exceedance estimates. New datasets for North American critical loads for acidity for forest soil water and aquatic ecosystems were combined with the ensemble deposition predictions to show a substantial decrease in the area and number of locations in exceedance between 2010 and 2016 (forest soils: 13.2 % to 6.1 %; aquatic ecosystems: 21.2 % to 11.4 %). All models agreed in the direction of the ensemble exceedance change between 2010 and 2016. The North American ensemble also predicted a decrease in both severity and total area in exceedance between the years 2010 and 2016 for eutrophication-impacted ecosystems in the USA (sensitive epiphytic lichen: 81.5 % to 75.8 %). The exceedances for herbaceous community richness also decreased between 2010 and 2016, from 13.9 % to 3.9 %. The uncertainty associated with the North American eutrophication results is high; there were sharp differences between the models in both predictions of total N deposition and the change in N deposition, and hence in the predicted eutrophication exceedances between the two years. The European ensemble was used to predict relatively static exceedances of critical loads with respect to acidification (4.48 % to 4.32 % from 2009 to 2010) while eutrophication exceedance increased slightly (60.2 % to 62.2 %). While most models showed the same changes in critical load exceedances as the ensemble between the two years, the spatial extent and magnitude of exceedances varied significantly between the models. The reasons for this variation were examined in detail by first ranking the relative contribution of different sources of sulphur and nitrogen deposition in terms of deposited mass and model-to-model variability in that deposited mass, followed by their analysis using AQMEII4 diagnostics, along with evaluation of the most recent literature. All models in both the North American and European ensembles had net annual negative biases with respect to observed wet deposition of sulphate, nitrate and ammonium. Diagnostics and recent literature suggest that this bias may stem from insufficient cloud scavenging of aerosols and gases, and may be improved through the incorporation of multiphase hydrometeor scavenging within the modelling frameworks. The inability of North American models to predict the timing of the seasonal peak in wet ammonium ion deposition (observed maximum was in April, while all models predicted a June maximum) may also relate to the need for multiphase hydrometeor scavenging (absence of snow scavenging in all models employed here). High variability in the relative importance of particulate sulphate, nitrate and ammonium deposition fluxes between models was linked to the use of updated particle dry deposition parameterizations in some models. However, recent literature and further development of some of the models within the ensemble suggests these particulate biases may also be ameliorated via the incorporation of multiphase hydrometeor scavenging. Annual sulphur and nitrogen deposition prediction variability was linked to SO 2 and HNO 3 dry deposition parameterizations, and diagnostic analysis showed that the cuticle and soil deposition pathways dominate the deposition mass flux of these species. Further work improving parameterizations for these deposition pathways should reduce variability in model acidifying gas deposition estimates. The absence of base cation chemistry in some models was shown to be a major factor in positive biases in fine mode particulate ammonium and particle nitrate concentrations. Models employing ammonia bidirectional fluxes had both the largest and the smallest magnitude biases, depending on the model and bidirectional flux algorithm employed. A careful analysis of bidirectional flux models suggests that those with poor NH 3 performance may underestimate the extent of NH 3 emissions fluxes from forested areas. Based on these results, an increased process-research focus is therefore recommended for the following model processes and on observations which may assist in model evaluation and improvement: multiphase hydrometeor scavenging combined with updated particle dry deposition, cuticle and soil deposition pathway algorithms for acidifying gases, base cation chemistry and emissions, and NH 3 bidirectional fluxes. Comparisons with satellite observations suggest that oceanic NH 3 emissions sources should be included in regional chemical transport models. The choice of land use database employed within any given model was shown to significantly influence deposition totals in several instances, and employing a common land use database across chemical transport models and critical load calculations is recommended for future work.

critical loads↗
DOE OSTI2025

Water Mass Transformation Budgets in Finite‐Volume Generalized Vertical Coordinate Ocean Models

Water Mass Transformation (WMT) theory provides conceptual tools that in principle enable innovative analyses of numerical ocean models; in practice, however, these methods can be challenging to implement and interpret, and therefore remain under-utilized. Our aim is to demonstrate the feasibility of diagnosing all terms in the water mass budget and to exemplify their usefulness for scientific inquiry and model development by quantitatively relating water mass changes, overturning circulations, boundary fluxes, and interior mixing. We begin with a pedagogical derivation of key results of classical WMT theory. We then describe best practices for diagnosing each of the water mass budget terms from the output of Finite-Volume Generalized Vertical Coordinate (FV-GVC) ocean models, including the identification of a non-negligible remainder term as the spurious numerical mixing due to advection scheme discretization errors. We illustrate key aspects of the methodology through the analysis of a polygonal region of the Greater Baltic Sea in a regional demonstration simulation using the Modular Ocean Model v6 (MOM6). We verify the convergence of our WMT diagnostics by brute-force, comparing time-averaged (“offline”) diagnostics on various vertical grids to timestep-averaged (“online”) diagnostics on the native model grid. Finally, we briefly describe a stack of xarray-enabled Python packages for evaluating WMT budgets in FV-GVC models (culminating in the new xwmb package), which is intended to be model-agnostic and available for community use and development.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1978

Discrimination of rock classes and alteration products in southwestern Saudi Arabia with computer-enhanced LANDSAT data

Digital LANDSAT MSS data for an area in the southwestern Arabian Shield were computer-enhanced to improve discrimination of rock classes, and recognition of gossans associated with massive sulphide deposits. The test area is underlain by metamorphic rocks that are locally intruded by granites; these are partly overlain by sandstones. The test area further includes the Wadi Wassat and Wadi Qatan massive sulphide deposits, which are commonly capped by gossans of ferric oxides, silica, and carbonates. Color patterns and boundaries on contrast-stretched ratio color composite imagery, and on complementary images constructed using principal component and canonical analyses transformations, correspond exceptionally well to 1:100,000 scale field maps. A qualitative visual comparison of information content showed that the ratio enhancement provided the best overall image for identification of rock type and alteration products.

Blodget, H. W.↗
NASA NTRS1983

Low-level flow conditions hazardous to aircraft

Low level flow conditions known to be hazardous to aircraft during takeoff/climbout and approach/landing operations are turbulence, wind shear, and vertical motion. These conditions can and frequently do occur separately and in combinations. The identification and selection were completed of representative data cases to determine magnitude, frequency, duration, and simultaneity of occurrence of turbulence (gustiness and gust factor), wind shear (speed and direction), and vertical motion (updraft and downdraft), along with temperature inversions. New representations of temporal and spatial variations in the atmospheric boundary layer were developed. Efforts continued relative to low level flow conditions where published results imply strong vertical shear with virtually no horizontal shear and where order of magnitude analyses of the equations of motion for an aircraft illustrates that low values of horizontal shear (along the flight path) are much more hazardous than larger values of vertical wind shear (altitude).

Alexander, M. B.↗
DOE OSTI2021

Alpha Shape Analysis (ASA) Framework for Post- Clustering Property Determination in Atom Probe Tomographic Data

While application of clustering algorithms to atom probe tomography data have enabled quantification of solute clusters in terms of number density, size, and subcomposition there exist other properties (e.g., volume, surface area, and composition) that are better determined by defining an interface between the cluster and the surrounding matrix. The limitation in composition results from an ion selection step where the expected matrix ion types are omitted from the cluster search algorithm to enhance the contrast between the matrix and cluster and to reduce the complexity of the search. Previously, composition determination within solute clusters has utilized a secondary envelopment and erosion step on top of conventional methods such as maximum separation. In this work, we present a novel stochastic method that combines the particle identification fidelity of a conventional clustering algorithm with the analytical flexibility of mesh-based approaches through the generation of alpha shapes for each identified cluster. Furthermore, the corresponding mesh accounts for concave components of the clusters and determines the volume and surface area of the clusters; additionally, the mesh boundary is utilized to update the total composition according to the internal ions.

42 ENGINEERING↗
NASA NTRS1989

Analytical modeling of structure-soil systems for lunar bases

The study of the behavior of granular materials in a reduced gravity environment and under low effective stresses became a subject of great interest in the mid 1960's when NASA's Surveyor missions to the Moon began the first extraterrestrial investigation and it was found that Lunar soils exhibited properties quite unlike those on Earth. This subject gained interest during the years of the Apollo missions and more recently due to NASA's plans for future exploration and colonization of Moon and Mars. It has since been clear that a good understanding of the mechanical properties of granular materials under reduced gravity and at low effective stress levels is of paramount importance for the design and construction of surface and buried structures on these bodies. In order to achieve such an understanding it is desirable to develop a set of constitutive equations that describes the response of such materials as they are subjected to tractions and displacements. This presentation examines issues associated with conducting experiments on highly nonlinear granular materials under high and low effective stresses. The friction and dilatancy properties which affect the behavior of granular soils with low cohesion values are assessed. In order to simulate the highly nonlinear strength and stress-strain behavior of soils at low as well as high effective stresses, a versatile isotropic, pressure sensitive, third stress invariant dependent, cone-cap elasto-plastic constitutive model was proposed. The integration of the constitutive relations is performed via a fully implicit Backward Euler technique known as the Closest Point Projection Method. The model was implemented into a finite element code in order to study nonlinear boundary value problems associated with homogeneous as well as nonhomogeneous deformations at low as well as high effective stresses. The effect of gravity (self-weight) on the stress-strain-strength response of these materials is evaluated. The calibration of the model is performed via three techniques: (1) physical identification, (2) optimized calibration at the constitutive level, and (3) optimized calibration at the finite element level (Inverse Identification). Activities are summarized in graphic and outline form.

Macari-Pasqualino, Jose Emir↗
NASA NTRS1981

Study project of intrusive rocks: States of Espirito Santo and Rio de Janeiro, south and east of Minas Gerais and southeast of the state of Sao Paulo

The feasibility of mapping intrusive rocks in polycyclic and polymetamorphic areas using the logic method for photointerpretation of LANDSAT and radar imagery was investigated. The resolution, scale and spectral characteristics of the imagery were considered. Spectral characteristics of the intrusive rock units mapped using image 100 were investigated. It was determined that identification of acidic and basic intrusive bodies and determination of their relationships with principal structural directions using the logic method was feasible. Tectonic compartments were subdivided into units according to their predominant lithographic types, ignoring stratigraphy. The principal directions of various foliations, faults, megafolds, and fractural systems were defined. Delineation of the boundaries of intrusive bodies mapped using the spectral characteristics of Image 100 imagery ws determined to be more accurate than visual analysis. A 1:500,000 scale map of intrusions in the areas studied was generated.

Dejesusparada, N.↗
NASA NTRS1986

Utility of remotely sensed data for identification of soil conservation practices

Discussed are a variety of remotely sensed data sources that may have utility in the identification of conservation practices and related linear features. Test sites were evaluated in Alabama, Kansas, Mississippi, and Oklahoma using one or more of a variety of remotely sensed data sources, including color infrared photography (CIR), LANDSAT Thematic Mapper (TM) data, and aircraft-acquired Thermal Infrared Multispectral Scanner (TIMS) data. Both visual examination and computer-implemented enhancement procedures were used to identify conservation practices and other linear features. For the Kansas, Mississippi, and Oklahoma test sites, photo interpretations of CIR identified up to 24 of the 109 conservation practices from a matrix derived from the SCS National Handbook of Conservation Practices. The conservation practice matrix was modified to predict the possibility of identifying the 109 practices at various photographic scales based on the observed results as well as photo interpreter experience. Some practices were successfully identified in TM data through visual identification, but a number of existing practices were of such size and shape that the resolution of the TM could not detect them accurately. A series of computer-automated decorrelation and filtering procedures served to enhance the conservation practices in TM data with only fair success. However, features such as field boundaries, roads, water bodies, and the Urban/Ag interface were easily differentiated. Similar enhancement techniques applied to 5 and 10 meter TIMS data proved much more useful in delineating terraces, grass waterways, and drainage ditches as well as the features mentioned above, due partly to improved resolution and partly to thermally influenced moisture conditions. Spatially oriented data such as those derived from remotely sensed data offer some promise in the inventory and monitoring of conservation practices as well as in supplying parameter data for a variety of computer-implemented agricultural models.

Pelletier, R. E.↗
NASA NTRS2008

Services, Perspective and Directions of the Space Physics Data Facility

The multi-mission data and orbit services of NASA's Space Physics Data Facility (SPDF) project offer unique capabilities supporting science of the Heliophysics Great Observatory and that are highly complementary to other services now evolving in the international heliophysics data environment. The VSPO (Virtual Space Physics Observatory) service is an active portal to a wide rage of distributed data sources. CDAWeb (Coordinated Data Analysis Web) offers plots, listings and file downloads for current data from many missions across the boundaries of missions and instrument types. CDAWeb now includes extensive new data from STEREO and THEMIS, plus new ROCSAT IPEI data, the latest data from all four TIMED instruments and high-resolution data from all DE-2 experiments. SSCWeb, Helioweb and out 3D Animated Orbit Viewer (TIPSOD) provide position data and identification of spacecraft and ground conjunctions. OMNI Web, with its new extension to 1- and 5-minute resolution, provides interplanetary parameters at the Earth's bow shock. SPDF maintains NASA's CDF (Common Data Format) standard and a range of associated tools including format translation services. These capabilities are all now available through web services based APIs, one element in SPDF's ongoing work to enable heliophysics community development of Virtual discipline Observatories (e.g. VITMO). We will demonstrate out latest data and capabilities, review the lessons we continue to learn in what science users need and value in this class of services, and discuss out current thinking to the future role and appropriate focus of the SPDF effort in the evolving and increasingly distributed heliophysics data environment.

McGuire, Robert E.↗
NASA NTRS2012

Anomaly Detection in Test Equipment via Sliding Mode Observers

Nonlinear observers were originally developed based on the ideas of variable structure control, and for the purpose of detecting disturbances in complex systems. In this anomaly detection application, these observers were designed for estimating the distributed state of fluid flow in a pipe described by a class of advection equations. The observer algorithm uses collected data in a piping system to estimate the distributed system state (pressure and velocity along a pipe containing liquid gas propellant flow) using only boundary measurements. These estimates are then used to further estimate and localize possible anomalies such as leaks or foreign objects, and instrumentation metering problems such as incorrect flow meter orifice plate size. The observer algorithm has the following parts: a mathematical model of the fluid flow, observer control algorithm, and an anomaly identification algorithm. The main functional operation of the algorithm is in creating the sliding mode in the observer system implemented as software. Once the sliding mode starts in the system, the equivalent value of the discontinuous function in sliding mode can be obtained by filtering out the high-frequency chattering component. In control theory, "observers" are dynamic algorithms for the online estimation of the current state of a dynamic system by measurements of an output of the system. Classical linear observers can provide optimal estimates of a system state in case of uncertainty modeled by white noise. For nonlinear cases, the theory of nonlinear observers has been developed and its success is mainly due to the sliding mode approach. Using the mathematical theory of variable structure systems with sliding modes, the observer algorithm is designed in such a way that it steers the output of the model to the output of the system obtained via a variety of sensors, in spite of possible mismatches between the assumed model and actual system. The unique properties of sliding mode control allow not only control of the model internal states to the states of the real-life system, but also identification of the disturbance or anomaly that may occur.

Solano, Wanda M.↗
NASA NTRSDate not supplied

Exploring rock-regolith interfaces in Jezero crater with Mars 2020 SHERLOC

The Perseverance rover successfully landed in Jezero crater, Mars in February 2021 at the Octavia E. Butler landing site and began its mission to explore and sample an ancient crater lake basin. Principal goals of the Mars 2020 mission include characterizing the geology of Mars and seeking signs of ancient microbial life via the spacecraft cameras and spectroscopic instruments onboard. The Scanning Habitable Environments with Raman and Luminescence for Organics and Chemicals (SHERLOC) instrument is a deep UV Raman spectrometer that utilizes a 248.6nm pulsed laser. Part of SHERLOC is a color camera known as the Wide Angle Topographic Sensor for Operations and eNgineering (WATSON). The SHERLOC suite provides coordinated, spectroscopic and imaging capabilities at high spatial resolution, to detect minerals and organic molecules in microtextural context. By pairing high spatial resolution (~100 μm) resonance Raman and native fluorescence spectroscopy with microscopic imaging in a novel spacecraft capability, SHERLOC enables texture-specific molecular composition measurements of rock and regolith targets on Mars. Coordinated rock-regolith observations illuminate unique insights into weathering processes and thereby to primary properties of rocks in Jezero crater. This work describes the potential of rock-regolith interfaces to preserve unique records of geological processes in Jezero crater and can powerfully supplement observations of the more general rock record on Mars. Linking observations of local rock texture with associated regolith reveals important lithologic information based on the interrelationship between differential weathering behavior and mineralogy, grain size, and cement chemistry. Preliminary observations indicate that the polygonally fractured lithotype common near the Octavia E. Butler landing site may weather by granular disintegration and/or surface creep, a relation that can be uniquely observed at the rock-regolith interface. SHERLOC -specific observations of microtextural and elemental composition transitions presented here trace rock-regolith boundaries at multiple indurated surfaces adjacent to regolith. At these locales, grain-scale based examinations suggest chemical weathering could be related to a variably distributed coating or rind on dark rock targets that may be mafic in composition. Granule deposits overlying widely distributed fine-grained material are also observable. Mineral identifications of each phase are presented, with cross-scale comparisons to the remote insights gained by the SuperCam instrument.

Emily L. Cardarelli↗
NASA NTRSDate not supplied

Trajectory Engineering with Modular Patched Conics for Entry Systems and TPS (TEMPEST)

Brief Presenter Biography (35 word limit): Bohdan Wesely is an Aerospace Engineer in the Entry Systems and Technology Division at Ames. He has worked on a variety of projects for NASA including integrated TPS (thermal protection system) flight hardware deliveries and testing services for commercial partners. Introduction: TEMPEST is a new trajectory analysis framework that is designed to fill the gap between dedicated flight mechanics tools and aerothermal and TPS sizing tools. The project started as an SJSU master’s thesis and has since evolved into a general conceptual design tool capable of studying a wide variety of entry problems. Development is ongoing in the Entry Systems and Technology Division at NASA ARC. Why TEMPEST: Space missions involving entry into a planetary atmosphere involve a series of unique requirements across multiple disciplines. Whether it is traditional entry descent and landing (EDL), or aerocapture, the vehicle must navigate to its target landing location or orbit state, and the TPS must protect the payload during entry. The design process typically involves iterative handoffs between various flight mechanics, flow solver, and material response level tools. During the early conceptual phase, a wide variety of feasible trajectories are simulated in a Monte Carlo scenario which broadly satisfy the mission or landing requirements. Next, computational fluid dynamics (CFD), direct simulation Monte Carlo (DSMC), and other flow solver analyses are performed at various key trajectory points to generate an aero-database, heating and TPS design requirements also emerge at this stage. At this point, with updated aerodynamics from the various flow solvers, trajectories can be re-run, this in turn can change the required freestream conditions for the CFD tools, and as a project progresses, these analyses converge, and uncertainty is reduced. However, there is always a “hand-off” occurring between two inherently coupled phenomena. Analysis Description: One of the goals with TEMPEST is to use a variety of first principles estimation methods coupled with an atmosphere model to predict vehicle aerothermodynamics across the entire flight regime while propagating a 3 or 6 degree of freedom (DoF) trajectory. Aerodynamics methods include modified Newtonian, Maxwell and Cercignani- Lampis-Lord (CLL) for continuum, transitional, and free molecular flow regimes. Aerothermodynamics include boundary layer and reference enthalpy methods, and Mutation++ for non-equilibrium chemistry modeling. TEMPEST is also capable of stitching multiple trajectory segments together to study mission scenarios like multi-pass aerocapture and aero-gravity assists. Most of the program is implemented in MATLAB using modern system objects, it relies on several C++ shared libraries for supporting tools like Gmsh, the Global Reference Atmospheric Model (GRAM), and Mutation++. The various first principles aerothermal estimation methods are discretized across either a structured axisymmetric panel mesh or an unstructured tri-mesh generated from an open-source tool such as Gmsh, this allows solutions on the same mesh to be compared across tools such as CB-Aero. CFD Coupling. A physics-aware, gaussian process CFD anchoring scheme is proposed to adjust the various first principles methods as a CFD database is populated. One goal for this anchoring module is to inform the project where CFD should be run. Full knowledge of the entire trajectory, atmosphere, and aerothermodynamics allows for easier identification of high sensitivity areas and uncertainty quantification. While the first principles effects are well known and proven accurate in existing tools such as CB- Aero and Cart3D, a physics aware CFD anchoring scheme increases tool credibility across a project lifecycle. Material Response Modeling. Correct TPS sizing is critical for optimizing mass for science payloads and ensuring mission success. The process typically involves a thermal analysis along the trajectory with surface heating environments as a boundary condition. Several design constraints are maximum bondline temperature and maximum recession with various margining techniques. The material response tool FIAT, developed out of NASA Ames, is currently being integrated into the TEMPEST environment. TPS recession, shape change, mass loss, and mass property alteration are all factors that can perturb an entry trajectory. For missions like Mars 2020, recession was minimal and was safely handled separately as a post process. For missions such as Jupiter Galileo with a high TPS mass fraction or asteroid entries, recession plays a major role. The proposed fully coupled scheme is to use an epoch-based approach where the trajectory integration is halted after a recession threshold, the energy balance and FIAT are solved at each panel, the mesh, aerodynamics, and mass properties are updated, and the trajectory continues. Several computational tradeoffs have been made during the development of TEMPEST to limit the cost of a single trajectory and preserve its utility as a conceptual, rapid iteration tool. Conclusion: Development of TEMPEST is ongoing and the project is still in its infancy. This talk aims to showcase its unique capabilities to support future NASA entry systems missions.

Bohdan O Wesely↗
DOE OSTI2025

Fast meta-solvers for 3D complex-shape scatterers using neural operators trained on a non-scattering problem

Three-dimensional target identification using scattering techniques requires high accuracy solutions and very fast computations for real-time predictions in some critical applications. We first train a deep neural operator (DeepONet) to solve wave propagation problems described by the Helmholtz equation in a domain without scatterers but at different wavenumbers and with a complex absorbing boundary condition. We then design two classes of fast meta-solvers by combining DeepONet with either relaxation methods, such as Jacobi and Gauss-Seidel, or with Krylov methods, such as GMRES and BiCGStab, using the trunk basis of DeepONet as a coarse-scale preconditioner. We leverage the spectral bias of neural networks to account for the lower part of the spectrum in the error distribution while the upper part is handled inexpensively using relaxation methods or fine-scale preconditioners. The meta-solvers are then applied to solve scattering problems with different shape of scatterers, at no extra training cost. We first demonstrate that the resulting meta-solvers are shape-agnostic, fast, and robust, whereas the standard standalone solvers may even fail to converge without the DeepONet. We then apply both classes of meta-solvers to scattering from a submarine, a complex three-dimensional problem. We achieve very fast solutions, especially with the DeepONet-Krylov methods, which require orders of magnitude fewer iterations than any of the standalone solvers.

97 MATHEMATICS AND COMPUTING↗