Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “co-simulation”

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.

810 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

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↗
NASA NTRSDate not supplied

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↗
DOE OSTI2024

Simulating hindered grain boundary diffusion using the smoothed boundary method

Abstract Grain boundaries can greatly affect the transport properties of polycrystalline materials, particularly when the grain size approaches the nanoscale. While grain boundaries often enhance diffusion by providing a fast pathway for chemical transport, some material systems, such as those of solid oxide fuel cells and battery cathode particles, exhibit the opposite behavior, where grain boundaries act to hinder diffusion. To facilitate the study of systems with hindered grain boundary diffusion, we propose a model that utilizes the smoothed boundary method to simulate the dynamic concentration evolution in polycrystalline systems. The model employs domain parameters with diffuse interfaces to describe the grains, thereby enabling solutions with explicit consideration of their complex geometries. The intrinsic error arising from the diffuse interface approach employed in our proposed model is explored by comparing the results against a sharp interface model for a variety of parameter sets. Finally, two case studies are considered to demonstrate potential applications of the model. First, a nanocrystalline yttria-stabilized zirconia solid oxide fuel cell system is investigated, and the effective diffusivities are extracted from the simulation results and are compared to the values obtained through mean-field approximations. Second, the concentration evolution during lithiation of a polycrystalline battery cathode particle is simulated to demonstrate the method’s capability.

Materials Science↗
NASA NTRSDate not supplied

The Use of Immersion Rigs for High Temperature Hydrogen Exposure Testing Within the Nuclear Thermal Rocket Element Environmental Simulator (NTREES): NTREES-5 Hydrogen Heat Exchanger (N5HHX-01)

The Nuclear Thermal Rocket Element Environmental Simulator (NTREES) facility at NASA’s Marshall Space Flight Center (MSFC) subjects materials and components typically found on nuclear thermal propulsion (NTP) systems to non-nuclear simulations of prototypic NTP reactor thermohydraulic environments. With 1.2 MW of inductive power available, the currently employed methods to simulate nuclear heating in NTREES consist of 1) inductively heating a test specimen, which typically cannot be performed on ceramic materials, or 2) heating gas that envelops the test specimen. The latter of these has provided a basis for developing so-called “immersion rigs” to be used in NTREES. Furthermore, immersion rigs are especially useful in the case of smaller, technically challenging test specimens. One such device, the “NTREES-5 Hydrogen Heat Exchanger” (N5HHX-01), designed and deployed to test small material coupons to ~0.5 g/s of 2800 K, 500 psi Hydrogen, is discussed herein. To achieve this desired delivery temperature level, the most suitable material for the N5HHX’s construction was pure tungsten, which imposed budgetary and manufacturing challenges. Also discussed in this paper is the design of the N5HHX rig’s sample container, which became the first iteration of improving NTREES small material coupon test operations by permitting rapid changeout of different material samples.

Reactor Testing↗
DOE OSTI2026

Passive Confirmation of the Presence Of High-Explosive Material Via Neutron Transmission Spectroscopy

The goal of this project was to explore a novel approach for identifying presence and potential types of high explosives (HE) in treaty-controlled items with passive neutron sources by using passive neutron transmission spectroscopy. We predominantly look to measure relative elemental abundance of Carbon, Hydrogen, Oxygen, and Nitrogen (CHON) since most relevant materials are certain mix of these elements, as shown in Table 1. Previously, most studies exploring potential identification of CHON elemental content of targets in the vicinity of passive neutron source were focused solely on gamma spectroscopy using high resolution gamma detectors. Using neutron transmission spectroscopy with pulse-shape discrimination (PSD) capable organic scintillators is not only a novel idea in of itself, but arguably necessary for accurate identification of the CHON elemental content of an interrogated target. While high resolution gamma spectroscopy in principle can determine a presence of hydrogen and nitrogen, it is effectively blind to a quantitatively measuring areal densities of C, O. The initially proposed approach, described in detail in the project proposal, was to leverage previous Monte Carlo study on neutron transmission spectroscopy using slab shaped target interrogated with well-collimated (“pencil”) neutron beam. In this study, a simulated test CO 2 target was interrogated by a pencil neutron beam and transmitted neutron spectra were measured using PSD capable organic scintillator using MLEM based unfolding technique. To establish relative elemental content of carbon and oxygen, the unfolded neutron spectrum was fit with a parametrized combination of their respective elemental spectral templates. The individual spectral template was calculated by convoluting ENDF neutron crosssection with the known resolution of the PSD detector used in the study. This approach in the studied configuration successfully quantitatively established relative elemental composition of carbon and oxygen.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRSDate not supplied

Effect of Barely Visible Impact Damage on Thermoplastic Welded Structure

Spirit Aerosystems recently conducted screening tests on thermoplastic skin/stringer panels manufactured with Automated Fiber Placement (AFP) and stamp forming which were joined using Spirit’s Co-Fusion process. Skin and stringer gauge were both set to 0.168 cm. The resulting panel design is 48 cm tall, 49 cm wide, and has 3 stringers welded at a 21.5 cm pitch. Barely Visible Impact Damage (BVID) and Residual Compressive Strength were chosen to conduct a damage tolerance study. Zee stringers, made with IM8/PEKK-FC material, were Co-Fusion welded to a T1100/TC1225 skin. Out of these 5 panels, 1 was chosen to be a BVID survey panel, 1 for pristine compressive strength determination, and 3 to be impacted to create BVID followed by testing for residual compressive strength. On the survey panel, impacts ranged in energy from 40J – 65J and dent depth measurements were conducted at 24 hours, 48 hours, and two weeks after impacting. The residual compressive strength panel failed near the predicted load of 214 MPa. For the residual compressive strength panels, BVID was created at two different mid-bay locations: 1) striking the OML coincident with the stringer and 2) striking the skin IML 2.5 cm away from the stringer flange. Pulse Echo Ultrasonic Testing was used to inspect the panels before and after creating BVID.

HiCAM↗
DOE OSTI2026

Cave-based chronology of valley incision in the glaciated Northern Carpathians

While glaciers leave a mark on mountain landscapes, the importance of glacial erosion in mountains at a global scale remains insufficiently constrained despite over a century of research aimed at quantifying the rates and magnitudes of such erosion. Measuring rates of landscape change across the Quaternary in repeatedly glaciated mountain landscapes requires the ability to accurately partition the relative role of glacial versus fluvial processes in forming relief; however, the modification of landforms by surface erosion hinders quantification of process rates over the last few million years. This study presents 32 paired 26 Al and 10 Be cosmogenic nuclide burial ages from allochthonous, quartz-bearing sediment deposited within a paleophreatic karst system in the northern Carpathians. These data quantify Pliocene-to-present timing of valley incision and karst system reactivation, as well as erosion rates over time. Post-burial, muon production-correction yields 24 fully resolved ages, six minimum ages, and two unresolvable ages. Burial ages are interpreted to document the tempo of valley incision via phreatic passage formation and reveal that ∼80% of current valley relief (∼700 m) had developed by approximately 1.6 Ma. Paleoerosion rates derived from this dataset suggest little change in the processes generating sediment washed into the karst networks since the Pliocene. These results suggest that river incision dominated relief creation in Carpathian valleys that bear the imprint of glaciation. Our findings suggest the importance of investigating other mid-latitude mountainous regions worldwide to determine if they have similar erosion and topographic histories.

Carpathians↗
DOE OSTI2023

CFD Simulation of the Dosing Behavior within the Atomic Layer Deposition Feeding System

The effective operation of atomic layer deposition (ALD) feeding system is the premise of realizing specific ALD processes. In the present work, a detailed computational fluid dynamics (CFD) model of the feeding system has been developed and validated, which accounts for the roles of ALD valves and manifolds. A numerical simulation of the compressible fluid flow and heat/mass transfer within the feeding system was conducted. The dosing amounts and the spatiotemporal distributions of the precursors can be accurately predicted using the CFD model, as validated by experimental results. Different precursors, operating conditions, and structures of the feeding system were simulated and analyzed to examine the operating flexibility of the feeding system. The simulation results can be adopted as the upstream boundary conditions for simulations of the ALD process in the reaction chamber. The substrate-scale simulation indicates that the effect of the feeding system on the film deposition is highly related to the surface kinetics of ALD. The present work can serve as a guide for the development and optimization of different ALD-based processes via proper operation and even the design of the feeding system.

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

Ka-Band Monopulse Antenna-Pointing Systems Analysis and Simulation

NASA's Deep Space Network (DSN) has been using both 70-m and 34-m reflector antennas to communicate with spacecraft at S-band (2.3 GHz) and X-band (8.45 GHz). To improve the quality of telecommunication and to meet future mission requirements, JPL has been developing 34-m Ka-band (32-GHz) beamwave guide antennas. Presently, antenna pointing operates in either the open-loop mode with blind pointing using navigation predicts or the closed-loop mode with conical scan (conscan). Pointing accuracy under normal conscan operating conditions is in the neighborhood of 5 mdeg. This is acceptable at S- and X-bands, but not enough at Ka-band. Due to the narrow beamwidth at Ka-band, it is important to improve pointing accuracy significantly (approx. 2 mdeg). Monopulse antenna tracking is one scheme being developed to meet the stringent pointing-accuracy requirement at Ka-band. Other advantages of monopulse tracking include low sensitivity to signal amplitude fluctuations as well as single-pulse processing for acquisition and tracking. This article presents system modeling, signal processing, simulation, and implementation of Ka-band monopulse tracking feed for antennas in NASA/DSN ground stations.

V Y Lo↗
DOE OSTI2026

Development of Ternary Transition Metal Oxide Catalysts for Oxygen Evolution Reaction

Electrochemical water splitting, a promising method for green hydrogen production, is currently hindered by the high cost of precious-group metal (PGM)-based catalysts for the oxygen evolution reaction (OER). This study addresses this challenge by advancing the development of catalysts for OER, focusing on the development of high-performance PGM-free catalysts using nickel-iron-cobalt (NiFeCo)-based aerogels. The catalysts were synthesized via a sol-gel method and critical point drying to achieve a highly porous structure with an exceptionally high surface area. The designed catalyst structure provides an ideal platform for maximizing catalytic active sites and enhancing mass transport kinetics. Co has been systematically incorporated into the current PGM-free state-of-the-art, NiFeOx catalyst, and the metal ratios have been optimized. In addition to the experimental studies, density functional theory calculations were performed to study the material’s properties of this ternary catalyst and the effect of Co addition on enhancing OER catalysis.

Mizrahi, Michal [Bar-Ilan University, Ramat Gan, I↗
DOE OSTI2024

Factors Influencing Gas Evolution from High‐Nickel Layered Oxide Cathodes in Lithium‐Based Batteries

Abstract Gas evolution from high‐nickel layered oxide cathodes (>90% Ni) remains a major issue for their practical application. Gaseous species, such as CO 2 , O 2 , and CO, that are evolved at high states of charge (SOC) worsen the overall safety of batteries, as pressure build‐up within the cell may lead to cell rupture. Since these gasses are produced during cathode degradation, tracking the formation of gasses is also important in diagnosing cathode failure. Online electrochemical mass spectrometry (OEMS) is a powerful in situ technique to study gas evolution from the cathode during high‐voltage charge. However, the differences in the OEMS experimental setups between different groups make it challenging to compare results between groups. In this perspective, the various factors that influence gas evolution based on the OEMS results collected in this group are presented. The focus is on the conditions that lead to gas release, with a particular emphasis on reactive oxygen formation and subsequent chemical reactions with the electrolyte. Promising strategies, such as electrolytes, compositional tuning, and surface coatings that are effective at suppressing gas evolution from the cathode are highlighted. Critical insights into mitigating cathode degradation and gas evolution are provided to guide the development of safer, high‐energy batteries.

Chemistry↗
DOE OSTI2026

An integrated in-situ coordination strategy enabling high-performance layered cathodes for sodium-ion batteries

O3-type layered transition metal oxide cathodes hold tremendous potential in sodium-ion batteries (SIBs) due to their low cost and high energy density. However, the structure instability associated with detrimental phase transitions and severe interface parasitic reactions exacerbate the material's electrochemical performance degradation. Herein, we develop an integrated in-situ coordination strategy via heteroatomic modulation inducing coherent epitaxial layer to collaboratively enhance the overall framework robustness from surface to bulk. The theoretical calculation and multiple in/ex-situ characterizations demonstrate the charge density around oxygen is redistributed, which promotes the electron localization, thus widening the NaO 2 lattice space and accelerating the Na + transport dynamics. Furthermore, the formed strengthened oxygen bond energy effectively distributes the long-range coordination of Mn 3+ O 6 octahedron, thereby alleviating Jahn-Teller distortion and local stress. Importantly, the in-situ formed conformal buffer layer dramatically relieves the adverse interface side reactions, facilitating the construction of robust cathode-electrolyte interface, which ameliorate the whole structure stability of designed materials. Consequently, the optimized NFMZ@NZO-1.0 exhibits the excellent cycling stability with 80.2% capacity retention after 300 cycles at 1C, and delivers a high discharge capacity of 107.1 mAh g −1 at 10C. In conclusion, this distinctive coupling strategy provides valuable insights for developing high-performance layered cathode materials in SIBs.

Coherent epitaxial layer↗
NASA NTRS2026

Isotopic Evidence for A Cold and Distant Origin of 3I/Atlas

Interstellar objects provide the only directly observable samples of icy planetesimals formed around other stars, and can therefore provide insight into the diversity of physical and chemical conditions occurring during exoplanet formation1−3. Here we report isotopic measurements of the interstellar comet 3I/ATLAS, which reveal an elemental composition unlike any Solar System body. The water in 3I/ATLAS is enriched in deuterium, at a level of D/H = (0.98 ± 0.06)%, which is more than an order of magnitude higher than in known comets, while its range of 12 C/ 13 C ratios (141–191 for CO 2 and 123–172 for CO) exceeds typical values found in the Solar System, as well as nearby interstellar clouds and protoplanetary disks. Such extreme isotopic signatures indicate formation at temperatures ≲ 30 K in a relatively metal-poor environment. When interpreted with respect to models for Galactic chemical evolution, the carbon isotopic composition implies that 3I/ATLAS may have accreted as long ago as 12 billion years, following a period of intense, early star formation. 3I/ATLAS thus represents a preserved fragment of an ancient planetary system.

Martin A Cordiner↗
DOE OSTI2025

Oxidant-assisted methane pyrolysis

Adding small amounts of CO 2 or H 2 O to methane pyrolysis boosts hydrogen and carbon yieldsviacyclic catalyst regeneration, enabling efficient low-carbon hydrogen and crystalline carbon production.

Chemistry↗
NASA NTRS2019

Physics-Based Modeling and Simulation of Emerging Battery Technologies for Aerospace

Recently there is a growing interest in the aviation sector to reduce air and noise pollution. Electrochemical energy storage devices such as batteries coupled with a distributed electric propulsion system can reduce noise concerns as well as emissions and allow the concept of Urban Air Mobility to come to fruition. The battery performance needs to improve considerably compared from current state-of-art Li-ion battery (about 200Wh/Kg) to realize all-electric passenger jet for short flights (up to 690 miles). NASA is exploring lithium-oxygen battery chemistry to power hybrid (battery powered electrical system) and all-electric aircraft for short distance and long-distance flights. Li-O2 is one of the advanced Li-ion technologies that promise to provide specific energy of more than 750Wh/Kg. For this presentation, we present our work on improving power density of Li-O2 batteries through the use of multiphysics simulations. Next, a path is outlined to port these model to simulate performance for a new battery chemistry for space application, Li-CO2. Li-CO2 uses carbon dioxide as the active material instead of oxygen. Although this technology is in its early development, the offers two benefits: it can be used as a CO2 scrubber, as oxygen is one of the by-products on charging, and as a backup or a standalone battery for various Mars or Venus missions, where the carbon dioxide content in the atmosphere is high and need battery to operate at higher temperatures.

Mehta, Mohit↗
DOE OSTI2025

Numerical simulation of frost formation and heat transfer on fin-and-tube heat exchangers in turbulent cross-flow

Frost formation in fin-and-tube heat exchangers in turbulent cross-flow presents significant challenges in industrial refrigeration applications, affecting heat transfer efficiency and operational reliability. The purpose of this work is to investigate frost deposition and growth on a staggered bank of a fin-and-tube freezer coil under turbulent forced convection conditions. The focus here is on investigating conditions that closely replicate real-world scenarios in large walk-in industrial freezers. Using a direct numerical simulation approach, we examine the flow dynamics and thermal behaviour in the presence of frost, considering turbulent regimes characterized by a Reynolds number in the range 1050 ≤ R e D , avg ≤ 4800 , with the characteristic length being the outer diameter of the tube and the velocity being the bulk fluid velocity between the plates (fins). Computational fluid dynamics simulations are employed to resolve the interactions between turbulent airflow and the frost layer. Our approach incorporates a modified immersed boundary method and a slow-time acceleration technique to address the complex dynamic interface between the continuously evolving frost layer and the flowing air stream. Our findings indicate that frost forms more on the sides of the finned surfaces (plates) and less on the tubes themselves. This article is part of the theme issue ‘Heat and mass transfer in frost and ice’.

Science & Technology - Other Topics↗
DOE OSTI2026

Entropy-defect synergy for dual luminescence mechanism in spinel: Time-resolved anti-counterfeiting and fingerprint visualization

Multimodal luminescent materials, while promising for anti-counterfeiting, often lack dynamic time-dependent responses and controllable spatial distribution, limiting their encryption capabilities in the spatiotemporal dimension. Here, this work presents a coordinated control strategy based on entropy and defect engineering, and uses a backpropagation (BP) neural network for material screening to successfully prepare spinel Mg 0.8 (Fe 0.04 Co 0.04 Ni 0.04 Cu 0.04 Zn 0.04 )Cr 2 O 4 (MgA 5 CO) phosphors with time-dependent dynamic luminescence behavior. This phosphor simultaneously activated the d-d transition luminescence (∼618 nm) derived from Co 2+ /Cr 3+ and the defect luminescence (∼398 nm) related to zinc vacancies (V Zn ) in a single-phase solid solution. The phosphor exhibits a time-dependent color evolution from pink to purple under fixed-wavelength excitation, due to the different excited-state dynamics and decay lifetimes associated with the d-d transition and defect luminescence. Structural characterization and spectral analysis confirmed the existence of V Zn and its significant role in defect luminescence process. The fluorescent and dynamic luminescent properties of entropy-based spinel oxide enable its use in advanced anti-counterfeiting applications like fingerprint recognition and color-changing dedicated anti-counterfeiting mark, showing promise in high-end and time-dynamic anti-counterfeiting fields. This research not only developed a new type of fluorescent dynamic anti-counterfeiting material, but also provided a new idea for constructing advanced optical functional materials with multiple luminescence mechanisms.

Defect project↗