Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Modeling of Metal AM”

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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Quantitative modeling of gate-tunable cavity exciton-polaritons in transition-metal dichalcogenides

We performed a comprehensive real-space simulation study of propagative exciton-polaritons (EPs) of tungsten disulfide (WS 2 ) inside a Fabry–Pérot cavity. Through our simulations, we were able to construct the dispersion relation and determine the propagation length (L p ) of EPs, both of which show sensitive dependence on the applied gate voltage. We found that L p increases at higher gate voltages when excited at the exciton energy but shows opposite gate dependence at the trion energy. These observations are direct evidence of the coupling between the cavity mode with excitons and trions of WS 2 . Furthermore, we demonstrate a proof-of-concept device based on the tunable cavity EPs of WS 2 , which is promising for application as photonic transistors, modulators, or electrical-to-optical converters. Furthermore, our work unveils the real-space transport properties of gate-tunable cavity EPs and paves the way for future applications of EPs in nanophotonic devices and circuits.

42 ENGINEERING↗
DOE OSTI2021

CFD/FEA Co-Simulation Framework for Analysis of the Thermal Barrier Coating Design and Its Impact on the HD Diesel Engine Performance

Thermal barrier coatings (TBCs) have been investigated both experimentally and through simulation for mixing controlled combustion (MCC) concepts as a method for reducing heat transfer losses and increasing cycle efficiency, but it is still a very active research area. Early studies were inconclusive, with different groups discovering obstacles to realizing the theoretical potential. Nuanced papers have shown that coating material properties, thickness, microstructure, and surface morphology/roughness all can impact the efficacy of the thermal barrier coating and must be accounted for. Adding to the complexities, a strong spatial and temporal heat flux inhomogeneity exists for mixing controlled combustion (diesel) imposed onto the surfaces from the impinging flame jets. In support of the United States Department of Energy SuperTruck II program goal to achieve 55% brake thermal efficiency on a heavy-duty diesel engines, this study sought to develop a deeper insight into the inhomogeneous heat flux from mixing controlled combustion on thermal barrier coatings and to infer concrete guidance for designing coatings. To that end, a co-simulation approach was developed that couples high-fidelity computational fluid dynamics (CFD) modeling of in-cylinder processes and combustion, and finite element analysis (FEA) modeling of the thermal barrier-coated and metal engine components to resolve spatial and temporal thermal boundary conditions. The models interface at the surface of the combustion chamber; FEA modeling predicts the spatially resolved surface temperature profile, while CFD develops insights into the effect of the thermal barrier coating on the combustion process and the boundary conditions on the gas side. The paper demonstrates the capability of the framework to estimate cycle impacts of the temperature swing at the surface, as well as identify critical locations on the piston/thermal barrier coating that exhibit the highest charge temperature and highest heat fluxes. In addition, the FEA results include predictions of thermal stresses, thus enabling insight into factors affecting coating durability. An example of the capability of the framework is provided to illustrate its use for investigating novel coatings and provide deeper insights to guide future coating design.

42 ENGINEERING↗
DOE OSTI2023

Weak-coupling theory of neutron scattering as a probe of altermagnetism

Inelastic neutron scattering provides a powerful probe of the magnetic excitations of quantum magnets. Altermagnets have recently emerged as a new class of magnets with vanishing net magnetization characteristic of antiferromagnets and with a spin-split electronic structure typical of ferromagnets. Here we introduce a minimal Hubbard model with two-sublattice orthorhombic anisotropy as a framework to study altermagnetism. Using unrestricted Hartree-Fock calculations, we find an altermagnetic state for this model that evolves from a metallic state to an insulating state with increasing Hubbard-U Coulomb repulsion. We then examine the inelastic neutron scattering response in these states using random-phase approximation calculations of the dynamic spin susceptibility χ"(q,ω). We find that the magnetic excitation spectrum depends on its chirality for q along certain directions in reciprocal space, an observation that may be used in inelastic neutron scattering experiments as a probe of altermagnetism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2016

Silicon Emitter Needle and Array Design for Indium Electrospray Arrays for Spacecraft Propulsion

A design for silicon microfabricated emitter arrays was developed for electrospray thrusters for indium propellant in compact architectures with scalable thrust, low mass and volume, high specific impulse and high efficiency operation. The design elements include tip height, height uniformity across the array, tip radii, axial groove depth, tip angle, emitter shank sidewall angle and number of emitters. They were derived from commercial liquid metal ion source designs and modeling, fabrication and test results. The most critical results of the emitter array design analysis suggest that microfabricated silicon emitter array emitters should have a height greater than 280 microns with a height uniformity of +/-10 microns, a tip half angle of 49° for uniform turn-on voltages and current across the array, low beam divergence and high mass utilization efficiency. Elements of the design were fabricated and demonstrated in single emitter and in 400 element emitter arrays to validate the design, fabrication and performance capability. The design height and uniformity of arrays was demonstrated for 85% of the emitters in a prototype array. Required tip and sidewall angles and groove depths have been microfabricated. Single microfabricated silicon emitters demonstrated better performance in current and voltage characteristics than commercially available liquid metal ion sources. Microfabricated emitters and arrays demonstrated the required current and stability performance to enable the MEP thruster development for indium propellant.

Marrese-Reading, Colleen M.↗
DOE OSTI2024

Modeling the Progenitor Stars of Observed Type IIP Supernovae

Abstract Type IIP supernovae (SNe IIP) are thought to originate from the explosion of massive stars >10M ⊙ . Their luminosity is primarily powered by the explosion energy and the radioactive decay energy of 56 Co, with the photosphere location regulated by hydrogen recombination. However, the physical connections between SNe IIP and their progenitor stars remain unclear. This paper presents a comprehensive study of SNe IIP and their progenitor stars by using the one-dimensional stellar evolution code,MESA. Our model grids consider the effects of stellar metallicity, mass, and rotation in the evolution of massive stars, as well as the explosion energy and 56 Ni production in modeling supernovae. To elucidate the observed SNe IIP and their origins, we compare their light curves (LCs) with our models. Furthermore, we investigate the impact of stellar parameters on LCs by considering stellar mass, metallicity, rotation, explosion energy, and 56 Ni production. We find that more massive stars exhibit longer plateaus due to increased photon diffusion time caused by massive ejecta. Higher metallicity leads to increased opacity and mass loss of progenitor stars. Rapid rotation affects internal stellar structures, enhancing convective mixing and mass loss, potentially affecting the plateau’s brightness and duration. Higher explosion energy results in brighter but shorter plateaus due to faster-moving ejecta. 56 Ni mass affects late-time luminosity and plateau duration, with larger masses leading to slower declines.

Astronomy & Astrophysics↗
DOE OSTI2021

The effect of aspect ratio on the mechanical behavior of Li metal in solid-state cells

The Li metal anode is necessary for solid-state batteries to dramatically exceed the energy density of state-of-the-art conventional Li-ion batteries. To enable the Li metal anode, the mechanical properties of Li metal under compression in a commercially relevant cell configuration must be determined. This work examines the compression creep behavior of Li metal coupled with a relevant model solid-state electrolyte (Li 6.5 La 3 Ta 0.5 Zr 1.5 O 12 ) and current collector (Ni) as a function of applied stress (1–10 MPa) and Li foil thickness (100–700 μm). At 1 MPa applied stress, the strain rate of all Li foil thicknesses under compression was almost four orders of magnitude lower than what is predicted from tensile creep tests. Steady-state creep behavior was observed for Li foil ≥350 μm at 5 MPa and ≥250 μm at 10 MPa. Conversely, sustained creep deformation was not observed for 100 and 150 μm Li foils at any applied stress. It is suggested that friction at the Li/LLZO and Li/Ni interfaces leads to the formation of hydrostatic stress that controls Li deformation at low aspect ratios. Finally, these results have implications for battery pack stack pressure requirements, particularly relating to void formation during stripping and Li metal filament penetration during plating.

25 ENERGY STORAGE↗
DOE OSTI2023

Improved melt model for power flow

Accelerators that drive z-pinch experiments transport current densities in excess of 1 MA/cm 2 in order to melt or ionize the target and implode it on axis. These high current densities stress the transmission lines upstream from the target, where rapid electrode heating causes plasma formation, melt, and possibly vaporization. These plasmas negatively impact accelerator efficiency by diverting some portion of the current away from the target, referred to as “current loss”. Simulations that are able to reproduce this behavior may be applied to improving the efficiency of existing accelerators and to designing systems operating at ever higher current densities. The relativistic particle-in-cell code CHICAGO ® is the primary code for modeling power flow on Sandia National Laboratories’ Z accelerator. We report here on new algorithms that incorporate vaporization and melt into the standard power-flow simulation framework. Taking a hybrid approach, the CHICAGO® kinetic/multi-fluid treatment has been expanded to include vaporization while the quasi-neutral equation-of-motion has been updated for melt at high current-densities. For vaporization, a new one-dimensional substrate model provides a more accurate calculation of electrode thermal, mass, and magnetic field diffusion as well as a means of emitting absorbed contaminants and vaporized metal ions. A quasi-fluid model has been implemented expressly to mimic the motion of imploding liners for accurate inductance histories. For melt, a multi-ion Hall-MHD option has been implemented and benchmarked against Alegra MHD. This new model is described with sufficient detail to reproduce these algorithms in any hybrid kinetic code. Physics results from the new code are also presented. A CHICAGO ® Hall-MHD simulation of a radial transmission line demonstrates that Hall physics, not included in Alegra, has no significant impact on the diffusion of electrode material. When surface contaminant desorption is mocked in as a hydrogen surface plasma, both the surface and bulk-material plasmas largely compress under the influence of the j × B force. Similar results are seen in Alegra, which also shows magnetic and material diffusion scaling with peak current. Test vaporization simulations using MagLIF and a power-flow experimental geometry show Fe + ions diffuse only a few hundred µm from the electrodes, so present models of Z power flow remain valid.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Development of neural network force fields for corrosion studies

To fully understand the chemistry and physics of corrosion, novel methods of simulation must be developed. One approach is designing machine learning (ML) algorithms integrated with density functional theory to develop adaptive force fields to gain insight into corrosion behavior namely at the surface of metal oxides. Current methods of modeling corrosion are slow due to the computational cost of resolving both reaction mechanics and mass transport processes. Machine learning methods can be implemented to obtain structure-activity relationships at both the molecular and bulk scale while still retaining the accuracy of density functional theory (DFT) and significantly decreasing the time needed for simulations of complex chemical processes in the various environments of corrosion. Multiscale models are needed for corrosion studies to fully understand its processes not only at the atomic length scale (chemical bonding, energies, and forces), but also at the nano and meso length scales (solid-state physics and material science processes). Current methods of study include DFT, molecular dynamics, and Monte Carlo. The limitation of DFT is that only a small number of atoms or molecules can be simulated at that level of theory. Density functional theory is used to study the electronic structure of atoms and molecules, and calculate the force component of each atom. However, these calculations are limited to about 1000 atoms. Custom periodic boundary conditions (PBC) can be used to describe the various environments and defects that affect the atomic forces to produce a large data set from which a training set can be derived. Machine learning can be utilized to overcome the barrier of modeling macroscopic and multi-scale processes from ab initio calculations through the development of adaptive force fields. Local environments determine the atomic forces of a given system, therefore adaptive force fields must be created to produce reliable quantum mechanical calculations. This can be achieved by developing a learning algorithm that uses the mapped atomic forces or fingerprint as an input to produce energies and magnetic moments as output. A systematic approach was used to begin to build a data set in order to accurately describe the atomic forces in various environments. In Figure 4 below, a simple PBC cell of Fe{sub 2}O{sub 3} was first optimized. A surface optimization was performed next, followed by a hydroxylated surface optimization. Once this calculation has converged, the adsorption of halide species to the hydroxylated surface will be investigated. TensorFlow is an open source platform for machine learning developed by Google. Using a high level application program interface (API) such as Keras allows for building and training ML models easily in a number of different environments and languages. For this project, a neural network was developed within Anaconda in Python. Future Work: Further development of reference data set; Refining neural network and learning algorithm; Fingerprinting atomic environment to enable mapping of atomic force components; Choosing appropriate training set from reference data; Learning from training set and enabling non-linear mapping of training set fingerprints and the atomic forces; Estimation of uncertainty to identify ranges of outside applicability; Testing and analysis of molecular dynamic simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Assessing the Role of Interfacial and Metal Sites in Pt/TiO2-Catalyzed Acetic Acid Hydrodeoxygenation

Thermochemical conversion of biomass to produce drop-in quality biofuels typically involves hydrodeoxygenation (HDO) steps following catalytic fast pyrolysis (CFP) to remove excess oxygen and create a more-stable bio-oil product. HDO reactions are performed by co-feeding the CFP vapor-phase product and H2 gas over a bi-functional catalyst. Noble-metal catalysts supported on reducible metal oxides (e.g., Pt/TiO2) are active and selective toward these HDO reactions. Griffin and co-workers showed that Pt/TiO2 catalysts promote the desired deoxygenation steps for m-cresol HDO while mitigating undesired C-C bond-breaking steps that reduce the overall value/energy density of the biofuel. Such model-compound studies to inform the design of improved catalysts for HDO chemistry are necessary to improve the overall process economics/efficiencies for biofuels production. One important class of bio-derived compounds that has not been studied extensively with respect to HDO chemistry, particularly at the atomic level, is carboxylic acids. Past experimental work indicates that Pt/C and Pt/TiO2 catalysts are selective toward C-C and C-O bond-dissociation products for acetic acid HDO, respectively. To understand the role of Pt/TiO2 active sites in this observed change in selectivity and guide catalyst development, the work in this presentation focuses on modeling the role of Pt-metal and Pt-TiO2-interfacial sites in promoting key C-C bond-breaking, C-O bond-breaking, and (de)hydrogenation steps in acetic acid HDO. Density functional theory (DFT) calculations were performed using the Vienna Ab initio Simulation Package (VASP). The exchange-correlation functional was approximated by the Perdew-Burke-Ernzerhof functional. Dispersion interactions were captured using the D3 method. Projector augmented-wave potentials described electron-ion interactions, and electron wavefunctions were expanded via a planewave basis with an energy cutoff of 500 eV. Activation barriers were calculated using the climbing image nudged elastic band method. Results and Discussion: To discern the role of Pt-metal and Pt-TiO2 interface sites in promoting acetic acid HDO chemistry, Pt(111) slab and anatase TiO2(101)-supported Pt-nanowire (PtNW/OH-TiO2) surface models were constructed, respectively. Because H2 is co-fed in HDO reactions, the anatase support was terminated with OH groups. Interfacial vacancies have been shown to facilitate Ru/TiO2-catalyzed phenol HDO; thus, an interfacial model with an OH vacancy was also considered (PtNW/OHv-TiO2). Pt-TiO2-interface sites stabilize adsorption of all studied acetic acid HDO surface intermediates relative to terrace Pt-metal sites, particularly when an interfacial-OH vacancy is present. Oxygenated species prefer to bind at the OH vacancy through the O atom, suggesting a preference for C-O over C-C bond cleavage at these sites. This hypothesis is supported by net-negative and net-positive average shifts in the reaction energy and activation energy barriers for C-O and C-C bond-breaking steps, respectively, at Pt-TiO2-interface sites relative to Pt-metal sites. Using the calculated energetics, the predicted minimum-energy pathway was determined for each surface. Pt(111) is predicted to follow decarboxylation, producing undesired methane and carbon dioxide. Conversely, PtNW/OH-TiO2 and PtNW/OHv-TiO2 are both predicted to produce desired acetaldehyde and ethane. The interfacial vacancy may also play a key role in facilitating the first C-O bond-breaking step in acetic acid HDO, lowering the barrier by 0.6 eV relative to the defect-free interface model. These results demonstrate the critical role of the Pt-TiO2 interface in the shift in acetic acid HDO selectivity experimentally observed on Pt/C and Pt/TiO2 catalysts. The results herein demonstrate the important role of the Pt-TiO2 interface and interfacial oxygen vacancies in improving the carbon efficiency for HDO reactions in CFP upgrading.

BIOMASS FUELS,INORGANIC, ORGANIC, PHYSICAL, AND AN↗
NASA NTRS1979

An exploratory investigation of the effect of a plastic coating on the profile drag of a practical-metal-construction sailplane airfoil

The Langley low-turbulence pressure tunnel to determine the effect of a plastic coating on the profile drag of a practical-metal-construction sailplane airfoil was investigated. The model was tested with three surface configurations: (1) filled, painted, and sanded smooth; (2) rough bare metal; and (3) plastic-coated. The results are compared with data for the design airfoil (Wortmann FX 67-K-170/17) from another low-turbulence wind tunnel. The investigation was conducted at Reynolds numbers based on airfoil chord of 1.1 x 10 to the 6th power, 2.2 x 10 to the 6th power, and 3.3 x 10 to the 6th power at a Mach number of 0.10.

Somers, D. M.↗
DOE OSTI2023

Utilization of Additive Manufacturing for the Rapid Prototyping of C-Band Radiofrequency Loads

Additive manufacturing is a versatile technique that shows promise in providing quick and dynamic manufacturing for complex engineering problems. Research has been ongoing into the use of additive manufacturing for potential applications in radiofrequency (RF) component technologies. Here, we present a method for developing an effective prototype load produced from 316L stainless steel on a direct metal laser sintering machine. The model was tested using simulation software to verify the validity of the design. The load structure was manufactured by an online digital manufacturing company, showing the viability of using easily accessible tools to manufacture RF structures. The produced load was able to produce an S11 value of −22.8 dB at a C-band frequency of 5.712 GHz while under a vacuum. In a high-power test, the load was able to terminate a peak power of 8.1 MW. The discussion includes future applications of the present method and how it will help to improve the implementation of future accelerator concepts.

36 MATERIALS SCIENCE↗
NASA NTRS2003

Data and Modeling of Dendrites Subject to A Step Change in Pressure (TDSE)

There is considerable interest in dendritic solidification because of the influence dendrites have in the determination of microstructure, and thereby in the physical properties of cast metals and alloys. Current theories and models of dendritic growth generally couple diffusion effects in the melt with the physics of the interface. Data and subsequent analysis prior of the tip growth speed and radii of thermal succinonitrile dendrites in the near-convection free, on-orbit, free-fall environment demonstrate that these theories yield predictions that are reasonably in agreement with the results of experiment. However, data and analysis for assessing the interfacial physics component of theory are not sufficiently detailed or definitive. To study fundamental aspects of dendritic interface stability, we are measuring and modeling the kinetics and morphology of dendrites as they evolve from one well-defined steady state at a pre-set supercooling, through a transient stage, to a different well-defined steady state. More specifically, we subject succinonitrile dendrites, growing under steady-state conditions, to a rapid change in pressure. This leads to a rapid change in thermal driving force from the corresponding change in both the equilibrium melting temperature due to the Clapeyron effect, and a change in the far-field temperature due to adiabatic temperature changes in the bulk liquid and solid. Subsequently, we observe transformations from a well-characterized initial state into a new steady-state. Initial data reveal that the dendrite tip velocity changes almost as fast as the pressure charges, while the tip radius changes occur more slowly, taking from 10 60 seconds depending on the size of the step change and the final supercooling. Computer modeling of this process shows both agreements and disagreements with the experimental data. In making these observations and measurements, we are gaining new understandings of interfacial dynamics and state-selection physics.

Koss, Matthew B.↗
NASA NTRS2005

Multiscale Modeling for the Analysis for Grain-Scale Fracture Within Aluminum Microstructures

Multiscale modeling methods for the analysis of metallic microstructures are discussed. Both molecular dynamics and the finite element method are used to analyze crack propagation and stress distribution in a nanoscale aluminum bicrystal model subjected to hydrostatic loading. Quantitative similarity is observed between the results from the two very different analysis methods. A bilinear traction-displacement relationship that may be embedded into cohesive zone finite elements is extracted from the nanoscale molecular dynamics results.

Glaessgen, Edward H.↗
DOE OSTI2021

AFRL Additive Manufacturing Modeling Series: Challenge 1, Characterization of Residual Strain Distribution in Additively-Manufactured Metal Parts Using Energy-Dispersive Diffraction

A machine component fabricated by additive manufacturing (AM) of metal powders often possesses steep residual stress gradients with significant magnitudes due to large temperature gradients that transpire within a localized area during fabrication. These processing-induced residual stresses can cause distortion of the part, and if they are of significant magnitude, they could induce cracking of the component, degrade printability, and/or diminish subsequent mechanical performance. The ability to predict these residual stresses imparted by AM is an important step in permeating AM technology for advanced manufacturing. Calibration and validation of AM process models used for prediction are, therefore, a critical step in understanding the origin and mitigating the challenges associated with residual stresses inherent to the AM process. In the present work, the residual strain distributions in components with simple geometries fabricated by a laser powder bed fusion (LBPF) process were characterized non-destructively using energy-dispersive X-ray diffraction in support of the US Air Force Research Laboratory Additive Manufacturing Modeling Challenge Series Groeber et al. (JOM 70:441-448, 2018). The measurement setup and approach are described in detail so that the data can be used as a benchmark to calibrate and validate models for the prediction of macro-scale residual stresses due to the LPBF process.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Text Mining for Process–Structure–Properties Relationships in Metals

With the advent of large language models (LLMs), the vast unstructured text within millions of academic papers is increasingly accessible for materials discovery—although significant challenges remain. While LLMs offer promising few- and zero-shot learning capabilities, particularly valuable in the materials domain where expert annotations are scarce, general-purpose LLMs often fail to address key materials-specific queries without further adaptation. To bridge this gap, fine-tuning LLMs on human-labeled data is essential for effective structured knowledge extraction (Liu in The Importance of Human-Labeled Data in the Era of LLMs, 2023). Here, in this study, we introduce a novel annotation schema designed to extract generic process–structure–properties relationships from scientific literature. We demonstrate the utility of this approach using a dataset of 128 abstracts, with annotations drawn from two distinct domains: high-temperature materials (Domain I) and uncertainty quantification in simulating materials microstructure (Domain II). Initially, we developed a conditional random field (CRF) model based on MatBERT—a domain-specific BERT variant—and evaluated its performance on Domain I. Subsequently, we compared this model with a fine-tuned LLM (GPT-4o from OpenAI) under identical conditions. Our results indicate that fine-tuning LLMs can significantly improve entity extraction performance over the BERT-CRF baseline on Domain I. However, when additional examples from Domain II were incorporated, the performance of the BERT-CRF model became comparable to that of the GPT-4o model. These findings underscore the potential of our schema for structured knowledge extraction and highlight the complementary strengths of both modeling approaches.

Materials science↗
NASA NTRS2007

Distortion and Residual Stress Control in Integrally Stiffened Structure Produced by Direct Metal Deposition

2-D thermo-mechanical model developed to characterize distortion and residual stresses in integral structure produced by DMD. Demonstrated as a tool to guide experimental development of DMD fabrication process for aero structures. Distortion and residual stresses are local to deposit. Most distortion develops during deposition of the first few layers; Little change in distortion or residual stresses after fifth deposit layer Most of distortion is localized just beneath the build. Thicker build plates and the use of build lands results in greatest decrease in levels of distortion. Pre-straining shown to reduce distortion. Difficult to implement, particularly for complex stiffener arrays. Clamp position has complex effect on distortion and stresses. Overall distortion reduced with decreasing clamp clearance. Larger clamp clearances induce bending. Use of pre-heat and active cooling show minor influence on panel distortion. Generate changes in thermal gradients in the build plate.

Lin, Shih-Yung↗