Engineering PapersSearch

SEARCH · Engineering Papers

Results for “approximation methods”

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.

619 records · Page 8

Simulations of Yarn Micro-Mechanics of Woven Heat Shield Materials

Carbon and phenolic fibers are commonly used in ablative thermal protection materials, such as 3-dimensional Mid-Density Carbon Phenolic (3MDCP), a 3D-woven composite comprised of mixed-fiber yarn bundles. Predicting the micro-mechanical response and fracture of twisted yarns composed of brittle and ductile fibers requires a modeling approach that captures per-fiber yielding, fiber fracture, and inter-fiber friction and contact. This work presents an extended bonded particle model (BPM) for discrete element method (DEM) simulation of fiber and yarn mechanics, implemented in LAMMPS. The model builds on the incremental bond formulation of Guo et al. and introduces a piecewise elasto-plastic constitutive law for axial extension, enabling representation of fibers that yield before failure. 3MDCP yarns were constructed using measured fiber radius distributions and helical twist geometry. Tensile simulations of single-ply 3MDCP yarns show good agreement with vender stress–strain results. Fiber breakage models also show details on yarn breakage propagration, centered radially in the yarn. Yarn breakage of multi-ply 3MDCP also matched experimental observations in per-ply breakage; however, predicted yarn breakage strength were found higher than experimental observations.

Discrete Element Method

Simulations of Yarn Micro-Mechanics of Woven Heat Shield Materials (3MDCP)

Carbon and phenolic fibers are commonly used in ablative thermal protection materials, such as 3-dimensional Mid-Density Carbon Phenolic (3MDCP), a 3D-woven composite comprised of mixed-fiber yarn bundles. Predicting the micro-mechanical response and fracture of twisted yarns composed of brittle and ductile fibers requires a modeling approach that captures per-fiber yielding, fiber fracture, and inter-fiber friction and contact. This work presents an extended bonded particle model (BPM) for discrete element method (DEM) simulation of fiber and yarn mechanics, implemented in LAMMPS. The model builds on the incremental bond formulation of Guo et al. and introduces a piecewise elasto-plastic constitutive law for axial extension, enabling representation of fibers that yield before failure. 3MDCP yarns were constructed using measured fiber radius distributions and helical twist geometry. Tensile simulations of single-ply 3MDCP yarns show good agreement with vender stress–strain results. Fiber breakage models also show details on yarn breakage propagration, centered radially in the yarn. Yarn breakage of multi-ply 3MDCP also matched experimental observations in per-ply breakage; however, predicted yarn breakage strength were found higher than experimental observations.

Woven

230Th/234U radiochronometry for uranium materials by alpha spectrometry for nuclear forensics analysis

We have developed an alpha spectrometry method for 230Th/234U radiochronometry to determine the model separation age for uranium materials. This method offers an alternative to the more commonly used, but significantly more resource intensive multi-collector inductively coupled plasma mass spectrometry (MC-ICP-MS) while providing similar accuracy and precision. The development included radiochemical separation of uranium and thorium, calibration of the 232U tracer and the alpha spectrometry measurement. The method was validated by analyzing a certified reference material, CRM 125-A, which is a uranium oxide pellet assay, isotopic, and radiochronometric standard. The results were compared to those obtained by our routine MC-ICP-MS 230Th/234U radiochronometry analysis method to assess the performance parameters for alpha spectrometry compared to the state-of-the-art technique.

Macsik, Zsuzsanna [Los Alamos National Laboratory

NASA’s SLS (Space Launch System) Rocket Ready for Artemis II Lunar Mission

In early 2026, NASA will launch the Artemis II mission, an approximately 10-day long lunar mission that will fly three NASA astronauts and one CSA (Canadian Space Agency) astronaut on a free-return trajectory around the Moon, following a one-day checkout of their Orion spacecraft in Earth orbit. The mission will be the first to launch astronauts aboard NASA’s Orion spacecraft and on top of the agency’s SLS (Space Launch System) rocket. The mission will also deploy four 12U CubeSats as secondary payloads from the Orion stage adapter, following Orion separation and departure. The payloads, developed by four of NASA’s international partners, will perform a variety of science and technology investigations. The SLS and Orion for the mission are currently stacked in the Vehicle Assembly Building (VAB) at NASA’s Kennedy Space Center in Florida and are undergoing final preparations to rollout to Launch Pad 39B for a tanking test before launch. In addition to preparations for the Artemis II mission, significant progress is being made on the SLS for the Artemis III mission, which is targeted to return astronauts to the lunar surface no later than 2029. Major components of the core stage and solid rocket boosters are already at NASA Kennedy undergoing build-up for the mission. Data from the Artemis II launch and mission, as well as progress to subsequent missions, as available, will be presented.

Bruce Askins

NASA’s SLS (Space Launch System) Rocket Ready for Artemis II Lunar Mission

In early 2026, NASA will launch the Artemis II mission, an approximately 10-day long lunar mission that will fly three NASA astronauts and one CSA (Canadian Space Agency) astronaut on a free-return trajectory around the Moon, following a one-day checkout of their Orion spacecraft in Earth orbit. The mission will be the first to launch astronauts aboard NASA’s Orion spacecraft and on top of the agency’s SLS (Space Launch System) rocket. The mission will also deploy four 12U CubeSats as secondary payloads from the Orion stage adapter, following Orion separation and departure. The payloads, developed by four of NASA’s international partners, will perform a variety of science and technology investigations. The SLS and Orion for the mission are currently stacked in the Vehicle Assembly Building (VAB) at NASA’s Kennedy Space Center in Florida and are undergoing final preparations to rollout to Launch Pad 39B for a tanking test before launch. In addition to preparations for the Artemis II mission, significant progress is being made on the SLS for the Artemis III mission, which is targeted to return astronauts to the lunar surface no later than 2029. Major components of the core stage and solid rocket boosters are already at NASA Kennedy undergoing build-up for the mission. Data from the Artemis II launch and mission, as well as progress to subsequent missions, as available, will be presented.

Bruce R Askins

Electrochemical Deconstruction of Ortho-Phthalate Plasticizers and Recovery of Plasticizing Moieties

Electrochemical methods for the direct reduction of alkyl esters have been understudied but provide a potential advantage for addressing specific waste streams given that such strategies often require only a minimal chemical profile. One such waste stream that could benefit from electrochemical remediation is poly(vinyl chloride) (PVC) plastics. PVC recycling often faces challenges due to the complexity of such waste. Solvent-based recycling methods pose several advantages for addressing post-use PVC; however, such methods are complicated by the high amounts of toxic legacy plasticizers (ortho-phthalates) present in PVC. A potential solution to addressing such phthalates is to address them electrochemically, coupled with recovery of the resulting valuable products. Presented herein is an optimized electrochemical method for ester activation that is leveraged to separate and recover the aromatic and alkyl components of ortho-phthalate plasticizers. The resulting aliphatic alcohols may be reused to prepare other non-toxic plasticizers. This electro-degradation method is demonstrated on six phthalates that have been identified to pose health concerns in addition to plasticizers recovered from commercial samples of PVC.

di-(2-ethylhexyl) phthalate

Up-and-Coming IMCs

While the good oxidation and environmental resistance, high melting points, and comparatively low densities of such ordered intermetallics as Ti 3 Al, NiAl, FeAl, and NbAl 3 render them good candidates for advanced aerospace structures, their poor toughness at low temperatures and low strength at elevated temperatures have prompted the development of fiber-reinforced intermetallic-matrix composites (IMCs) with more balanced characteristics. Fabrication methods for continuous-fiber IMCs under development include the P/M 'powder cloth' method, the foil/fiber method, and thermal spraying. The ultimate success of IMCs depends on fibers truly compatible with the matrix materials.

Randy Bowman

Computing the Critical Temperature of the Affine-Transformed $D=3$ Ising Model Using Masked Autoregressive Flow

The simple Ising model provides a rich environment to build and study lattice field theories. As part of an ongoing project to construct a conformal field theory (CFT) on an arbitrarily curved manifold, in this work we develop methods to measure the critical temperature $β_c$ of the affine-transformed Ising model on the face-centered cubic (FCC) lattice. The main challenge in this endeavor is finding a computationally efficient and accurate method of interpolating and extrapolating Monte Carlo observables with respect to coupling coefficients and temperature. Herein, we compare two such methods. A traditional statistical approach uses the multiple histogram (MH) method, while a newer machine learning approach uses a masked autoregressive flow (MAF) to estimate the underlying probability density function of a set of observables. While the MH method is specifically designed to interpolate and extrapolate Monte Carlo observables, we find that MAF is a viable alternative for measuring $β_c$ with a computational cost that scales more favorably. Furthermore, we comment on additional advantages of MAF relevant to our work, such as extrapolating in system volume.

Svenson, Kai [Texas U.]

Quantum graph models for transport in filamentary switching

The formation of metallic nanofilaments bridging two electrodes across an insulator is a mechanism for resistive switching. Examples of such phenomena include atomic synapses, which constitute a distinct class of memristive devices the behavior of which is closely tied to the properties of the filament. Until recently, experimental investigation of the low-temperature regime and quantum transport effects has been limited. However, with growing interest in understanding the true impacts of the filament on device conductance, comprehending quantum effects has become crucial for quantum neuromorphic hardware. Here, we discuss quantum transport resulting from filamentary switching in a narrow region where the continuous approximation of the contact is not valid, and only a few atoms are involved. In this scenario, the filament can be represented by a graph depicting the adjacency of atoms and the overlap between atomic orbitals. Using the theory of quantum graphs with locally diffusive node scattering, we calculate the scattering amplitude of charge carriers on this graph and explore the interplay between filamentary formation and quantum transport effects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC

Photophoretic Propulsion Enabling Mesosphere Exploration NIAC Phase I Final Report

This Phase I report presents a comprehensive study on photophoretic flyers—innovative, ultralight, solar-powered vehicles that harness photophoretic forces generated via Knudsen pumping to achieve sustained flight in the mesosphere (50–80 km altitude). By integrating advanced materials such as nanocardboard— characterized by its extremely low areal density (~1 g/m²) and high bending stiffness—with ultrathin light-absorbing coatings, the project converts incident solar radiation directly into a directed thrust. Extensive experimental investigations, coupled with high-fidelity computational fluid dynamics (CFD) simulations using ANSYS Fluent, validate the concept across various three-dimensional geometries, including spherical, conical, and rocket-shaped configurations. These simulations bridge the gap between free-molecular and continuum flow regimes, demonstrating that optimized designs can generate lift forces sufficient to support kilogram-scale payloads even in low-pressure environments. At the heart of this innovation is the use of Knudsen pumping, where temperature gradients across porous surfaces induce directional gas flow, creating a modest overpressure that provides lift. The report introduces an analytical framework that interpolates between the well-known low-Reynolds number drag regime and the high-Reynolds number momentum theory. This model accurately predicts lift based on design parameters such as microchannel dimensions, porous wall geometry, areal density, and nozzle exit area. For instance, simulations indicate that 10-meter-scale structures with carefully engineered porous walls can achieve the necessary pressure differential to support scientifically significant payloads (~1 kg). The study also explores a hybrid propulsion approach that combines solar buoyancy with photophoretic lift. Initially, solar heating creates a buoyant force that elevates the flyer to mesospheric altitudes. Once in the optimal pressure range, the photophoretic mechanism—powered by Knudsen pumping—takes over as the primary source of lift, ensuring stable, long-duration flight. This dual-mode operation not only facilitates the deployment of photophoretic flyers but also broadens the potential applications for mesospheric exploration. In addition to propulsion, the report investigates the integration of photophoretic thrusters for trajectory control of existing research balloons in the upper stratosphere. By dynamically adjusting the nozzle orientation and controlling the flow-through velocity, these thrusters provide precise maneuverability, enabling the flyers to counteract atmospheric disturbances and adjust their flight paths in real time. For example, a photophoretic thruster approximately 7.5 by 7.5 meters in size could be unfolded below a payload gondola of a 60 million-cubic-foot zero-pressure balloon. Such a thruster can provide horizontal speed control of approximately 1 m/s using only sunlight and no moving parts (except those needed to track the Sun and control the jet direction). Importantly, photophoretic thrusters operate more efficiently at higher altitudes, which is complementary to known trajectory control techniques, such as propellers and tethered wings, which are more effective at lower altitudes. Finally, the report identifies three scientific research thrusts where mesospheric aircraft technology can have a profound impact: atmospheric tides, characterization of gravity waves, and investigation of mesospheric instabilities. Overall, the findings of this Phase I project represent a significant advancement in photophoretic propulsion technology. By demonstrating that large-scale, ultralight structures can be powered solely by solar radiation—via carefully engineered Knudsen pumping—this work lays a robust foundation for scalable, near-space flight architectures. Future refinements in material fabrication, structural optimization, and integrated trajectory control are expected to further enhance performance, paving the way for operational demonstrations that could revolutionize atmospheric science, remote sensing, and communication networks.

Knudsen Pump

High Capacity Battery Cell Flight Qualified

The High Capacity Battery Cell project is an effort equally funded by the NASA Lewis Research Center and Hughes Space and Communications Company (a unit of Hughes Aircraft Company) to develop and flight qualify a higher capacity nickel hydrogen battery for continuing use on commercial spacecraft. The larger diameter, individual pressure vessel cell will provide approximately twice the power, while occupying the same volume, as the current state-of-the-art nickel hydrogen cell. These cells are also anticipated to reduce battery cost by 20 percent. The battery is currently booked for use on 26 spacecraft, with the first flight scheduled in 1997. A strong requirement for batteries with higher power levels (6 to 12 kW), long life, and reduced cost was identified in studies of the needs of commercial communications spacecraft. With the design developed in this effort, the higher power level was accommodated without having to modify the rest of the existing spacecraft bus. This design scaled-up the existing state-of-the-art nickel hydrogen battery cell from a 3.5-in., 50-Ahr cell to a 5.5-in., 350-Ahr cell. An improvement in cycle life was also achieved by the use of the 26-percent KOH electrolyte design developed by NASA Lewis. The cell design was completed, and flight batteries were built and flight qualified by Hughes Space and Communications Company with input from NASA Lewis. Two batteries were shipped in September 1996 to undergo life cycle testing under the purview of NASA Lewis.

McKissock, Barbara I.

Reduced Order Models for Liquid Hydrogen Pooling and Vaporization Supported by Experiments

In the event of a leak of liquid hydrogen, a pool can form that vaporizes, disperses, and eventually dilutes to a non-flammable mixture. In this work, we describe fast-running models for the pooling and vaporization of liquid hydrogen in a steady cross-wind. Several pooling models from the literature are compared to solve for the flow and extent of the pool. The size of the pool can serve as the source for a separate dispersion model, which builds upon the existing one-dimensional Gaussian plume model in HyRAM+. Additional terms for the effects of a cross-wind on momentum and entrainment were added so that the model could handle the effects of a cross-wind on a low-speed flow. The models are compared to experimental data on pooling extent and downwind dispersion for steady flow rates of liquid hydrogen in a steady cross-wind. In the two compared experiments, liquid flow rates of 15 and 45 g/s were spilled onto concrete in cross-winds of approximately 1.8 m/s. The rate of growth of the pool and the downwind concentration boundaries are compared to the models, showing good agreement, although additional tuning is needed. These models can contribute to the advancement of codes and standards for liquid hydrogen systems.

dispersion

A First-Principles Study of the Structural and Thermo-Mechanical Properties of Tungsten-Based Plasma-Facing Materials

Tungsten (W) and tungsten alloys are being considered as leading candidates for structural and functional materials in future fusion energy devices. The most attractive properties of tungsten for the design of magnetic and inertial fusion energy reactors are its high melting point, high thermal conductivity, low sputtering yield, and low long-term disposal radioactive footprint. Despite these relevant features, there is a lack of understanding of how the structural and mechanical properties of W-based alloys are affected by the temperature in fusion power plants. In this work, we present a study on the thermo-mechanical properties of five W-based plasma-facing materials. First-principles density functional theory (DFT) calculations are combined with the quasi-harmonic approximation (QHA) theory to investigate the electronic, structural, mechanical, and thermal properties of these W-based alloys as a function of temperature. The coefficient of thermal expansion, temperature-dependent elastic constants, and several elastic parameters, including bulk and Young’s modulus, are calculated. Our work advances the understanding of the structural and thermo-mechanical behavior of W-based materials, thus providing insights into the design and selection of candidate plasma-facing materials in fusion energy devices.

42 ENGINEERING

Rapid Assessment of Sulfate Resistance in Mortar and Concrete

Extensive research has been conducted on the sulfate attack of concrete structures; however, the need to adopt the use of more sustainable materials is driving a need for a quicker test method to assess sulfate resistance. This work presents accelerated methods that can reduce the time required for assessing the sulfate resistance of mixtures by 70%. Class F fly ash has historically been used in concrete mixtures to improve sulfate resistance. However, environmental considerations and the evolving energy industry have decreased its availability, requiring the identification of economically viable and environmentally friendly alternatives to fly ash. Another challenge in addressing sulfate attack durability issues in concrete is that the standard sulfate attack test (ASTM C1012) is time-consuming and designed for only standard mortars (not concrete mixtures). To expedite the testing process, accelerated testing methods for both mortar and concrete mixtures were adopted from previous work to further the development of the accelerated tests and to assess the feasibility of testing the sulfate resistance of mortar and concrete mixtures rapidly. This study also established criteria for interpreting sulfate resistance for each of the test methods used in this work. A total of 14 mortar mixtures and four concrete mixtures using two types of Portland cement (Type I and Type I/II) and various supplementary cementitious materials (SCMs) were evaluated in this study. The accelerated testing methods significantly reduced the evaluation time from 12 months to 21 days for mortar mixtures and from 6 months to 56 days for concrete mixtures. The proposed interpretation method for mortar accelerated test results showed acceptable consistency with the ACI 318-19 interpretations for ASTM C1012 results. The interpretation methods proposed for the two concrete sulfate attack tests demonstrated excellent consistency with the ASTM C1012 results from mortar mixtures with the same cementitious materials combinations. Metakaolin was shown to improve sulfate resistance for both mortar and concrete mixtures, while silica fume and natural pozzolan had a limited impact. Using 15% metakaolin in mortar or concrete mixtures with Type I/II cement provided the best sulfate resistance.

Chemistry

Comparison Testing and Analysis of Braided Cords Constructed from Varying Kevlar Fiber Types

Material advancements move quickly, resulting in material obsolescence issues that require parachute engineers to understand how the material properties change with a material upgrade and if there are expected impacts to overall performance of the parachute. The objective of this test was to measure the elongation of different Kevlar types at both the cord-level and at the yarn-level to determine if a change in Kevlar type impacted parachute performance. An Instron machine with a pin-to-pin connection was used to measure load and displacement time histories for various Kevlar cord and yarn samples in a cyclic load profile. Twenty-four cord samples and twenty-four yarn samples were tested; each sample was cycled five times to the relevant maximum load (either to breaking strength or to the 30% peak load expected in flight). In addition to the Instron measurements, a photogrammetric technique was also implemented to measure elongation of the samples using button targets. Modulus test data shows approximately a 10% difference in elongation between Kevlar 29AP and Kevlar 129, and another 10% difference in elongation between Kevlar 29 and Kevlar 29AP. Although this is a measurable change in modulus, the effect on parachute performance is likely small, particularly for non-primary structure elements of relatively small length. This test campaign provides an example of the type of testing and analysis that parachute engineers could do to understand how material upgrades affect performance of their systems.

Modulus Testing

Elastic-plastic models for multi-site damage

This paper presents recent developments in advanced analysis methods for the computation of stress site damage. The method of solution is based on the p-version of the finite element method. Its implementation was designed to permit extraction of linear stress intensity factors using a superconvergent extraction method (known as the contour integral method) and evaluation of the J-integral following an elastic-plastic analysis. Coarse meshes are adequate for obtaining accurate results supported by p-convergence data. The elastic-plastic analysis is based on the deformation theory of plasticity and the von Mises yield criterion. The model problem consists of an aluminum plate with six equally spaced holes and a crack emanating from each hole. The cracks are of different sizes. The panel is subjected to a remote tensile load. Experimental results are available for the panel. The plasticity analysis provided the same limit load as the experimentally determined load. The results of elastic-plastic analysis were compared with the results of linear elastic analysis in an effort to evaluate how plastic zone sizes influence the crack growth rates. The onset of net-section yielding was determined also. The results show that crack growth rate is accelerated by the presence of adjacent damage, and the critical crack size is shorter when the effects of plasticity are taken into consideration. This work also addresses the effects of alternative stress-strain laws: The elastic-ideally-plastic material model is compared against the Ramberg-Osgood model.

Actis, Ricardo L.

Polynomial-time preparation of low-temperature Gibbs states for two-dimensional toric code

In this work, we propose a polynomial-time algorithm for preparing the Gibbs state of the two-dimensional toric code Hamiltonian at any temperature, starting from any initial state, significantly improving upon prior estimates that suggested exponential scaling with inverse temperature. We prove that fast mixing at low temperature for the two-dimensional toric code can be achieved by augmenting local jump operators with simple global jump operators, which enable efficient transitions between logical sectors. To establish tight lower bounds on the spectral gap, we introduce a new reduction method that eventually maps the problem to estimating the spectral gap of a perturbed graph Laplacian on a stair graph. Our proof also shows that the Lindblad dynamics with a digitally implemented low-temperature local Davies generator is able to efficiently drive the quantum state toward the ground state manifold.

97 MATHEMATICS AND COMPUTING

Reduced-order CFD modeling of cryogenic hydrogen isotope extrusion for pellet fueling

This study presents a reduced-order model (ROM) for computational fluid dynamics (CFD) simulations of cryogenic hydrogen isotope extrusions, focusing on protium (H₂) and deuterium (D₂) piston extruders. Using a 2D axisymmetric ROM in ANSYS-Polyflow, significant computational savings were achieved (runtime reduced from 9∼24 h to 3∼5 min), with extrusion force discrepancies between the 2D ROM and 3D models being on the order of 1%. Parametric studies identified optimal cutoff shear rates in the viscosity model (0.01/s for H₂ and 0.001/s for D₂), providing recommendations for future simulations. Finally, a comprehensive comparison of ROM results with experimental data was performed across varying geometries, cryogenic materials, temperatures, extrusion lengths, and piston velocities. Predictions at low extrusion temperatures met the objective of providing quick and efficient solutions with an acceptable extrusion force error of approximately 10% or less, validating the effectiveness of the 2D ROM approach. However, at high temperatures closer to the triple point, extrusion force error grows, which necessitates developing an improved model that accounts for temperature effects, e.g. melting. Nevertheless, the findings still represent a significant improvement in efficiency of CFD modeling of cryogenic hydrogenic extrusion. The ROM framework can also be extended to tritium (T2) and screw extruders, which will ultimately provide a fast and effective tool for optimizing pellet injector design for ITER and future reactor systems.

Fan, Joy [ORNL] (ORCID:0000000229751735)