Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Material Point Method”

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 793 records · Page 44

Anharmonic Eigenvectors and Acoustic Phonon Disappearance in Quantum Paraelectric SrTiO 3

Pronounced anomalies in the SrTiO 3 dynamical structure factor, S(Q,E), including the disappearance of acoustic phonon branches at low temperatures, were uncovered with inelastic neutron scattering (INS) and simulations. The striking effect reflects anharmonic couplings between acoustic and optic phonons and the incipient ferroelectric instability near the quantum critical point. It is rationalized using a first-principles renormalized anharmonic phonon approach, pointing to nonlinear Ti-O hybridization causing unusual changes in real-space phonon eigenvectors, frequencies, group velocities, and scattering phase space. Lastly, our method is general and establishes how T dependences beyond the harmonic regime, assessed by INS mapping of large reciprocal-space volumes, provide real-space insights into anharmonic atomic dynamics near phase transitions.

36 MATERIALS SCIENCE↗

Quantification of LEU Holdup using gamma ray imaging and inverse transport solver

Holdup is the residual amount of special nuclear material (SNM) remaining in a processing facility after the bulk materials have been cleaned out. In commercial uranium processing facilities, quantification of holdup is a major challenge because of the highly variable shapes and sizes of the deposits. Any method that attempts to generalize and calibrate deposit shapes in order to quantify holdup will be prone to high uncertainties. Uncertainties on the order of ±50% are typical in holdup results. In international safeguards applications, a ±50% uncertainty can result in a large amount of material unaccounted for (MUF) thereby increasing the difficulty of detecting material diversion and facility misuse. An imaging-based methodology has been developed with the objective of significantly reducing this uncertainty by using the true deposit shape, instead of relying on oversimplified geometric assumptions. The project is a collaboration between ORNL, Y-12, and the University of Tennessee, Knoxville, TN. Uranium sources of known masses were measured using the Germanium Gamma-ray Imager (GeGI), a high-resolution imaging spectrometer, creating a pixelated map for each spectral bin. Two different gamma imaging methods are employed in this work: coded aperture imaging and Compton imaging. A validated MonteCarlo model of the detector has been developed using the GEANT4 code for determining the intrinsic response of the detector, its enclosure, and the coded aperture mask. An inverse transport solver based on the Markov Chain Monte-Carlo approach known as Differential Evolution Adaptive Metropolis (DREAM) is employed to use the measurement data from the image pixels (coded aperture or Compton) to solve for the mass of 235 U in the deposit. A reliable method based on the DREAM solver has been developed to flag the infinite thickness condition of a uranium deposit. The project team is working towards improving the image reconstruction for Compton imaging so that a better localization of the source can be achieved. Besides treating the coded aperture and Compton imaging methods independently, the project is also evaluating a combined method that uses the Compton scatter data from a coded aperture measurement. GEANT4 simulations are being performed to evaluate the combined approach. The impact on the DREAM optimization as the source thickness progressively approaches infinite thickness is being evaluated. A number of uranium sources available at ORNL have been measured, and the DREAM results have been tested and validated for the coded aperture imaging. A similar effort will be carried out to validate the Compton based method once the development of algorithms for better localization are complete. The imaging based quantification is very amenable to unattended monitoring of holdup accumulation at key measurement points. A proof of concept measurement has been completed to demonstrate this capability The current work used the high energy resolution imager GeGI. However, the approach and methodologies are applicable to other imagers such as the cadmium zin telluride (CZT) based imager manufactured by H3D, Inc.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Unveiling the Finite Temperature Physics of Hydrogen Chains via Auxiliary Field Quantum Monte Carlo

The ability to accurately predict the finite temperature properties and phase diagrams of realistic quantum solids is central to uncovering new phases and engineering materials with novel properties ripe for device applications. Nonetheless, there remain comparatively few many-body techniques capable of elucidating the finite temperature physics of solids from first principles. Here, we take a significant step toward developing such a technique by generalizing our previous, exact fully ab initio finite temperature Auxiliary Field Quantum Monte Carlo (FT-AFQMC) method to model periodic solids and employing it to uncover the finite temperature physics of periodic hydrogen chains. Our chains’ unit cells consist of 10 hydrogen atoms modeled in a minimal basis, and we sample 5 k-points from the first Brillouin zone to arrive at a supercell consisting of 50 orbitals and 50 electrons. Based upon our calculations of these chains’ many-body energies, free energies, entropies, heat capacities, double and natural occupancies, and charge and spin correlation functions, we outline their metal–insulator and magnetic ordering as a function of both H–H bond distance and temperature. At low temperatures approaching the ground state, we observe both metal–insulator and ferromagnetic–antiferromagnetic crossovers at bond lengths between 0.5 and 0.75 Å. We then demonstrate how this low-temperature ordering evolves into a metallic phase with decreasing magnetic order at higher temperatures. In order to contextualize our results, we compare the features we observe to those previously seen in one-dimensional, half-filled Hubbard models at finite temperature and in ground state hydrogen chains. Interestingly, we identify signatures of the Pomeranchuk effect in hydrogen chains for the first time and show that spin and charge excitations that typically arise at distinct temperatures in the Hubbard model are indistinguishably coupled in these systems. Beyond qualitatively revealing the many-body phase behavior of hydrogen chains in a numerically exact manner without invoking the phaseless approximation, our efforts shed light on the further theoretical developments that will be required to construct the phase diagrams of the more complex transition metal, lanthanide, and actinide solids of longstanding interest to physicists.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermal-Mechanical Analysis of an Additive Manufacturing Ceramic Heat Exchanger for High-Temperature Recuperator in a sCO2 Power System

Supercritical CO2 (sCO2) Brayton power cycle can be configured in a closed-loop power system and has a potentially high cycle efficiency. Compactness and high efficiency of a sCO2 power block make the sCO2 Brayton cycle a versatile power cycle in broad applications. While many heat sources are of sufficient intensity to produce high temperature working fluids to achieve high cycle efficiency, the thermal-mechanical stability of traditional materials (e.g., steels and nickel-based superalloys) used in construction of heat exchangers and turbine components limits the operating conditions and thus thermodynamic efficiency of the system. This effort seeks to establish the viability of ceramic heat exchanger technologies for the most extreme operating conditions envisioned for power generation and other high temperature processes. Heat exchangers constructed from ultra-high temperature ceramics, a class of extreme environment materials featuring melting points (Tmp.) above 3000 degrees C, is particularly appealing for sCO2 Brayton cycles given their ultra-low creep rates and very high retained strength at low homologous temperatures (i.e., T < 0.5 Tmp., or at least 1500 degrees C). To translate these materials properties to ultra-high temperature heat exchangers, innovations are required in ceramic manufacturing techniques to realize the complex architectures featured in compact heat exchangers with high power density. With appropriate processing, ZrB2-SiC based compositions can be sintered to near full density and shaped into complex topologies via ceramic additive manufacturing methods. This paper analyzes heat exchanger designs and explores thermal-mechanical implications of the operating environment. Thermal flow, heat transfer, and conjugate mechanical analyses provide insights into benefits and risks associated with the design approach.

additive manufacturing↗

Uncertainty Analysis of Slug Calorimeters in the HyMETS Arc-Jet Facility

The objective of this work is to perform an uncertainty analysis of the deduced stagnation heat flux environment on a slug calorimeter for conditions that span the performance envelope of the Hypersonic Materials Environmental Test System arc-jet facility located at NASA Langley Research Center. Analytical solutions are developed for boundary-value problems on the slug element accounting for non-ideal effects, including spatial variation in the slug heat flux, multi-dimensional thermal conduction, and back-face losses, which departs from the state-of-the-art method derived from the American Society of Testing and Materials. Boundary-value problem definitions are informed by preliminary finite element thermal analysis of the slug calorimeter assembly (including both slug and housing) and just the slug element. The analytical solutions are presented in a general sense and in a truncated form from error analysis. Results are shown in optimizing and validating the analytical models against available slug back-face thermal data. The optimization results indicate that the appropriate epistemic uncertainty of the deduced stagnation heat flux on the slug calorimeter is at most±2.5% for both a high-and low-enthalpy test condition. In addition, a numerical approach is used to determine the aleatory (probabilistic) uncertainty component in the slug stagnation heat flux by applying a marching least-squares slope routine through the steady-state portion of the slug back-face thermal response. Results indicate a compromise between the number of samples and the filter frequency of slug back-face thermal data points when evaluating the standard deviation of the deduced stagnation heat flux statistics. When combining the mixed uncertainty, both aleatory and epistemic, the interval of uncertainty in the deduced stagnation heat flux is determined to be up to ±4%, which is at least a 60% reduction from the standard uncertainty used in the state-of-the-art method.

uncertainty↗

Application of the generalized reduced gradient method to conceptual aircraft design

The complete aircraft design process can be broken into three phases of increasing depth: conceptual design, preliminary design, and detail design. Conceptual design consists primarily of developing general arrangements and selecting the configuration that optimally satisfies all mission requirements. The result of the conceptual phase is a conceptual baseline configuration that serves as the starting point for the preliminary design phase. The conceptual design of an aircraft involves a complex trade-off of many independent variables that must be investigated before deciding upon the basic configuration. Some of these variables are discrete (number of engines), some represent different configurations (canard vs conventional tail) and some may represent incorporation of new technologies (aluminum vs composite materials). At Lockheed-Georgia, the sizing program is known as GASP (Generalized Aircraft Sizing Program). GASP is a large program containing analysis modules covering the many different disciplines involved fin defining the aricraft, such as aerodynamics, structures, stability and control, mission performance, and cost. These analysis modules provide first-level estimates the aircraft properties that are derived from handbook, experimental, and historical sources.

Gabriele, G. A.↗

Mitigate Stress Corrosion Cracking (SCC) in High-strength Al Castings (CRADA 508)

Development of high strength, corrosion resistant, aluminum alloys is important to a number of industries, including automotive, aerospace, power, etc. Understanding how processing methods such as casting, chemistry, additives and recycle content affect an alloy’s microstructure and performance in real world environments is crucial to their implementation. One major concern is how the repeated exposure to saline solutions can be detrimental to the structural integrity of aluminum alloy components. This study was undertaken as a collaboration between Eck Industries and PNNL to investigate the stress corrosion cracking (SCC) and corrosion behavior of aluminum alloys. This work examines the SCC response of three different materials set. First is cast A206 Al alloy produced by three different casting techniques (i.e., permanent mold, and sand casting with and without chill) and with or without nano forming additives as potential grain refiners. The second materials set comprises Al alloy tubes extruded by PNNL’s ShAPE process and produced using raw material with three different ratios of twitch and pre-consumer Al 6061 scrap. Finally, the third materials set comprises two Al-Si-Mg-Fe alloys with controlled Si% and Fe% to mimic recycle grade Al alloys. The corrosion behavior of these materials was characterized using SCC testing and scanning electrochemical cell microscopy (SECCM) technique; mechanical behavior was characterized using tensile testing, and microstructure was characterized using optical and electron microscopy techniques. Of the A206 Al alloys, the samples with nano forming additives performed poorly compared to those without. The best SCC performance (i.e. least degradation in mechanical properties) was demonstrated by A206 alloy sand cast with chill while the permanent mold sample showed the worst SCC performance. The SECCM was performed on A206 alloy sand cast with chill to understand the effect of microstructure and determine the location-dependent corrosion properties on grain boundaries (GB) and inside grains (IG) on the sample. The corrosion potential and current measured on GB and IG at various points on the sample established that the GB locations are more cathodic and corrosive than the IG locations. Further work needs to be conducted to investigate the mechanisms behind the observed dependence of SCC on casting technique and nano forming additives. Among the ShAPE extruded tubes, SECCM data suggests the decreasing order of tendency to corrode as follows: 100% twitch > 75% twitch > 50% twitch. Finally, both variants of the Al-Si-Mg-Fe alloys showed a large volume fraction of Fe- and Si-containing intermetallics and additional work is needed to discern differences in their respective corrosion responses. In summary, correlating local electrochemical behavior (e.g. via SECCM) with micro/nanostructure, in conjunction with bulk mechanical behavior, can help identify the right fabrication method and alloy chemistry to minimize SCC issues in Al alloys.

36 MATERIALS SCIENCE↗

Foundational insights into the mechanical and molecular evolution of porcine skin gelatin during gelation

Gelatin is a widely used material in biomedical fields, particularly in regenerative medicine and tissue engineering, due to its biocompatibility and versatile properties. While prior research has explored methods to enhance gelatin's mechanical strength and stability, fundamental studies on gelatin, specifically its curing process, mechanical stiffness, and chemical evolution during gelation, remain limited. This study uses ultrasonic testing and Fourier Transform Infrared Spectroscopy (FTIR) to examine gelatin's stiffness and molecular changes during gelation. Samples of 175 and 300 Porcine Skin Bloom Strength Gelatin at concentrations of 2% and 6% (w/v) were analyzed. Through transmission ultrasonic testing helped identify key transition points in gelation, with higher concentrations exhibiting delayed transitions. FTIR revealed that C-N bond formation peaks early while N-H bond deformation persists. A correlation emerged between sound speed and peak absorbance, suggesting that changes in molecular mobility may contribute to the observed sound speed behavior during periods of active bond formation. However, as gelation continues, fewer bonding components may be available, potentially decreasing molecular movement and contributing to the observed increase in sound speed. These findings provide insights into gelatin's mechanical and chemical evolution, offering a framework for improved control over its gelation kinetics. Swept-Frequency Acoustic Interferometry (SFAI) was performed at the end of the curing process to measure the sound speed, enabling the calculation of the bulk moduli of the gelatin samples. The combined use of ultrasonic and FTIR testing provides a non-destructive method for characterizing gelatin and other biomaterials. This approach advances understanding of gelatin curing behavior and supports the development of safer biomaterials with tailored mechanical properties for various applications such as tissue engineering and regenerative medicine.

Biomaterials↗

Surface Change Detection Using Large Footprint Laser Altimetry

Laser altimeters provide a precise and accurate method for mapping topography at fine horizontal and vertical scales. A laser altimeter provides range by measuring the roundtrip flight time of a short pulse of laser light from the laser altimeter instrument to the target surface. The range is then combined with laser beam pointing knowledge and absolute position knowledge to provide an absolute measurement of the surface topography. Newer generations of laser altimeters measure the range by recording the shape and time of the outgoing and received laser pulses. The shape of the return pulse can also provide unique information about the vertical structure of material such as vegetation within each laser footprint. Distortion of the return pulse is caused by the time-distributed reflections adding together and representing the vertical distribution of surfaces within the footprint. Larger footprints (10 - 100m in diameter) can support numerous target surfaces and thus provide the potential for producing complex return pulses. Interpreting the return pulse from laser altimeters has evolved from simple timing between thresholds, range-walk corrections, constant-fraction discriminators, and multi-stop time interval units to actual recording of the time varying return pulse intensity - the return waveform. Interpreting the waveform can be as simple as digitally thresholding the return pulse, calculating a centroid, to fitting one or more gaussian pulse-shapes to the signal. What we present here is a new technique for using the raw recorded return pulse as a raw observation to detect centimeter-level vertical topographic change using large footprint airborne and spaceborne laser altimetry. We use the correlation of waveforms from coincident footprints as an indication of the similarity in structure of the waveforms from epoch to epoch, and assume that low correlation is an indicator of vertical structure or elevation change. Thus, using vertically and horizontally geolocated waveforms as raw observables (i.e., waveforms tied to a common reference ellipsoid), we assess whether epoch-to-epoch vertical ground motion results in a decrease in the correlation of coincident waveforms over time, and whether this can be used to quantify the magnitude of the deformation. Results of computer models and an example over an area of eroded beachfront will be presented.

Blair, J. Bryan↗

Modeling Progressive Damage Using Local Displacement Discontinuities Within the FEAMAC Multiscale Modeling Framework

A method for performing progressive damage modeling in composite materials and structures based on continuum level interfacial displacement discontinuities is presented. The proposed method enables the exponential evolution of the interfacial compliance, resulting in unloading of the tractions at the interface after delamination or failure occurs. In this paper, the proposed continuum displacement discontinuity model has been used to simulate failure within both isotropic and orthotropic materials efficiently and to explore the possibility of predicting the crack path, therein. Simulation results obtained from Mode-I and Mode-II fracture compare the proposed approach with the cohesive element approach and Virtual Crack Closure Techniques (VCCT) available within the ABAQUS (ABAQUS, Inc.) finite element software. Furthermore, an eccentrically loaded 3-point bend test has been simulated with the displacement discontinuity model, and the resulting crack path prediction has been compared with a prediction based on the extended finite element model (XFEM) approach.

Ranatunga, Vipul↗

Intentional carbon doping reveals CH as an abundant charged impurity in nominally undoped synthetic WS 2 and WSe 2

Understanding the physical properties and controlling the generation of intrinsic and extrinsic defects is central to the technological adoption of 2D materials in devices. Here we identify a charged carbon-hydrogen complex at a chalcogen site (CH X ) as a common, charged impurity in synthetically grown transition metal dichalcogenides (TMDs). This conclusion is drawn by comparing high resolution scanning probe microscopy measurements of nominally undoped and intentionally carbon doped TMD samples. While CH impurity densities in undoped CVD-grown WS 2 and MOCVD-grown WSe 2 can range anywhere from parts per million to parts per thousand, CH densities in the percentage levels were selectively generated by a post-synthetic methane plasma treatment. Our study indicates that methane plasma treatment is a selective and clean method for the controlled introduction of a charged carbon-hydrogen complex at a surface chalcogen site, adefect that is commonly present in synthetic TMDs.

36 MATERIALS SCIENCE↗

Nuclear–electronic orbital approach to quantization of protons in periodic electronic structure calculations

The nuclear–electronic orbital (NEO) method is a well-established approach for treating nuclei quantum mechanically in molecular systems beyond the usual Born–Oppenheimer approximation. In this work, we present a strategy to implement the NEO method for periodic electronic structure calculations, particularly focused on multicomponent density functional theory (DFT). The NEO-DFT method is implemented in an all-electron electronic structure code, FHI-aims, using a combination of analytical and numerical integration techniques as well as a resolution of the identity scheme to enhance computational efficiency. After validating this implementation, proof-of-concept applications are presented to illustrate the effects of quantized protons on the physical properties of extended systems, such as two-dimensional materials and liquid–semiconductor interfaces. Specifically, periodic NEO-DFT calculations are performed for a trans-polyacetylene chain, a hydrogen boride sheet, and a titanium oxide–water interface. The zero-point energy effects of the protons as well as electron–proton correlation are shown to noticeably impact the density of states and band structures for these systems. These developments provide a foundation for the application of multicomponent DFT to a wide range of other extended condensed matter systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

SCALE 6.2.4 Validation: Radiation Shielding

For safe and reliable use of computer codes by the community, accuracy must be clearly evaluated. In particular, the nuclear reactor engineering and licensing field needs accurate tools for radiation shielding modeling. Monaco with Automated Variance Reduction using Importance Calculations (MAVRIC) is one such tool, with built-in variance reduction methods distributed within the SCALE code, and its validity is demonstrated in this report for the released version 6.2.4. Representative benchmarks corresponding to shielding analysis are selected for the validation study. Typical experimental results analyzed from those benchmarks include neutron fluxes, detector count rates, detector energy response functions, neutron and gamma doses, foil neutron activation rates and activities, neutron leakage fluxes, and skyshine dose rates. Thousands of points of comparison between experiment and calculation are presented in this work. Other than rare outliers typically explained by either a lack of information or large uncertainties in the experiment conditions, material, or dimensions, MAVRIC agrees well with the experiment results. MAVRIC is also compared to Monte Carlo N-Particle (MCNP) calculations when available, and both codes generally produce good agreements within estimated uncertainties. The selected benchmarks are obtained from reliable sources such as the International Criticality Safety Benchmark Evaluation Project Handbook (ICSBEP Handbook), the Shielding Integral Benchmark Archive & Database (SINBAD), and other shielding validation work found in the literature. Additional datapoints and benchmarks will be added to future versions of this report to incrementally expand the shielding validation suite incrementally.

61 RADIATION PROTECTION AND DOSIMETRY↗

Radiographic Moire

A method for the x-ray inspection of materials making use of the Moire effect is described. The Moire effect results when two patterns are superimposed; a third pattern is produced. Any change in either of the first two patterns creates a change in the third. Moire inspection is common with visible light. This invention allows the technique to be extended to locations inaccessible to visual inspection. A first pattern of high radio contrast material is attached to or included in the sample. X-rays are projected through the sample. A second pattern is imposed at the observation point, either before or after the formation of the x-ray image. The two patterns interact to create a third, Moire, pattern. As the material is stressed the Moire pattern changes, the degree of change indicating the degree of stress.

Madaras, Eric Irvine↗

Radiographic moire

A method for the x-ray inspection of materials making use of the Moire effect is described. The Moire effect results when two patterns are superimposed, a third pattern is produced. Any change in either of the first two patterns creates a change in the third. Moire inspection is common with visible light, this invention allows the technique to be extended to locations inaccessible to visual inspection. A first pattern of high radio contrast material is attached to or included in the sample. X-rays are projected through the sample. A second pattern is imposed at the observation point, either before or after the formation of the x-ray image. The two patterns interact to create a third, Moire, pattern. As the material is stressed the Moire pattern changes, the degree of change indicating the degree of stress.

Madaras, Eric I.↗

Devices and methods for creating plasma channels for laser plasma acceleration

This disclosure provides systems, methods, and apparatus related to devices and methods for creating hollow, near-hollow, and parabolic plasma channels. In one aspect, a device includes a block of material and a cooling system. The block of material defines a channel having a cylindrical shape and having a first open end and a second open end. An axis of the channel lies along a straight line. The block of material further defines a first gas port and a second gas port. The first gas port and the second gas port are in fluid communication with channel. The cooling system is operable to cool the channel to below the freezing point of a gas.

Leemans, Wim↗

A method for rapid liquid-solid phase solubility measurements using the protein lysozyme

Using hen's egg white lysozyme crystals as the test material, a simple system was developed for rapidly and unambiguously determining solubilities in (aqueous) solutions. The system is based upon a maximization of the solid surface area available for solute transfer to or from the solution, and a minimization of both the solution volume which must be equilibrated and the distance over which diffusive solute exchange occurs. This technique is further enhanced by using duplicate test systems which differ only in that one approaches equilibrium from an oversaturated solution, while the other from an undersaturated solution. Thus, the resulting data pair brackets the solubility value. In practical terms, the data points are found to usually be within 3 percent of each other, and individual solubility data points may usually be made at this resolution within 8-24 h depending upon the temperature change made since the previous determination.

Pusey, Marc L.↗

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↗