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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Electrochemical Activation of Ni–Fe Oxides for the Oxygen Evolution Reaction in Alkaline Media

The oxygen evolution reaction (OER) is essential to many key electrochemical devices, including H 2 O electrolyzers, CO 2 electrolyzers, and metal−air batteries. NiFe oxides have been historically identified as active for the OER, though they have been less studied in their more commercially relevant bulk oxide forms, such as NiFe 2 O 4 . Past works have demonstrated that the initial starting phase of Ni(Fe) precatalysts can influence their activation to the Ni(Fe)OOH active phase, including the rate and degree of conversion, pointing to the necessity of understanding activation protocols and in situ characteristics of catalyst materials at the device level. In this work, we investigate the characteristics of commercially relevant NiFe bulk oxides (NiFe 2 O 4 and a physical mixture of NiO and γ-Fe 2 O 3 ) during multiple activation procedures. Our results demonstrate that significant performance enhancement is observed for these bulk oxides regardless of the Fe incorporation in the initial form (i.e., atomically or macroscopically integrated), leading to significant performance enhancement (up to 30×) over time on stream. We hypothesize that this activation is due to the formation of NiFeOOH active sites on the surface, supported by in situ cyclic voltammetry and Raman spectroscopy results. We further show that not only the starting material but also the method of activation influences the number of Ni(Fe)OOH active sites formed and suggest that these sites can be quantified from the Ni 2+ to Ni 3+ redox transition using cyclic voltammetry. Broadly, this work demonstrates the necessity of in situ characterization of catalyst materials for cell-level design and testing.

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

Thermal and chemical convection in planetary mantles

Melting of the upper mantle and extraction of melt result in the formation of a less dense depleted mantle. This paper describes series of two-dimensional models that investigate the effects of chemical buoyancy induced by these density variations. A tracer particles method has been set up to follow as closely as possible the chemical state of the mantle and to model the chemical buoyant force at each grid point. Each series of models provides the evolution with time of magma production, crustal thickness, surface heat flux, and thermal and chemical state of the mantle. First, models that do not take into account the displacement of plates at the surface of Earth demonstrate that chemical buoyancy has an important effect on the geometry of convection. Then models include horizontal motion of plates 5000 km wide. Recycling of crust is taken into account. For a sufficiently high plate velocity which depends on the thermal Rayleigh number, the cell's size is strongly coupled with the plate's size. Plate motion forces chemically buoyant material to sink into the mantle. Then the positive chemical buoyancy yields upwelling as depleted mantle reaches the interface between the upper and the lower mantle. This process is very efficient in mixing the depleted and undepleted mantle at the scale of the grid spacing since these zones of upwelling disrupt the large convective flow. At low spreading rates, zones of upwelling develop quickly, melting occurs, and the model predicts intraplate volcanism by melting of subducted crust. At fast spreading rates, depleted mantle also favors the formation of these zones of upwelling, but they are not strong enough to yield partial melting. Their rapid displacement toward the ridge contributes to faster large-scale homogenization.

Dupeyrat, L.↗
NASA NTRS2002

A Two-Dimensional Linear Bicharacteristic FDTD Method

The linear bicharacteristic scheme (LBS) was originally developed to improve unsteady solutions in computational acoustics and aeroacoustics. The LBS has previously been extended to treat lossy materials for one-dimensional problems. It is a classical leapfrog algorithm, but is combined with upwind bias in the spatial derivatives. This approach preserves the time-reversibility of the leapfrog algorithm, which results in no dissipation, and it permits more flexibility by the ability to adopt a characteristic based method. The use of characteristic variables allows the LBS to include the Perfectly Matched Layer boundary condition with no added storage or complexity. The LBS offers a central storage approach with lower dispersion than the Yee algorithm, plus it generalizes much easier to nonuniform grids. It has previously been applied to two and three-dimensional free-space electromagnetic propagation and scattering problems. This paper extends the LBS to the two-dimensional case. Results are presented for point source radiation problems, and the FDTD algorithm is chosen as a convenient reference for comparison.

Beggs, John H.↗
NASA NTRS2007

Reentry Thermal Analysis of a Generic Crew Exploration Vehicle Structure

Comparative studies were performed on the heat-shielding characteristics of honeycomb-core sandwich panels fabricated with different materials for possible use as wall panels for the proposed crew exploration vehicle. Graphite/epoxy sandwich panel was found to outperform aluminum sandwich panel under the same geometry due to superior heat-shielding qualities and lower material density. Also, representative reentry heat-transfer analysis was performed on the windward wall structures of a generic crew exploration vehicle. The Apollo low Earth orbit reentry trajectory was used to calculate the reentry heating rates. The generic crew exploration vehicle has a graphite/epoxy composite honeycomb sandwich exterior wall and an aluminum honeycomb sandwich interior wall, and is protected with the Apollo thermal protection system ablative material. In the thermal analysis computer program used, the TPS ablation effect was not yet included; however, the results from the nonablation heat-transfer analyses were used to develop a "virtual ablation" method to estimate the ablation heat loads and the thermal protection system recession thicknesses. Depending on the severity of the heating-rate time history, the virtual ablation period was found to last for 87 to 107 seconds and the ablation heat load was estimated to be in the range of 86 to 88 percent of the total heat load for the ablation time period. The thermal protection system recession thickness was estimated to be in the range of 0.08 to 0.11 inches. For the crew exploration vehicle zero-tilt and 18-degree-tilt stagnation points, thermal protection system thicknesses of h = {0.717, 0.733} inches were found to be adequate to keep the substructural composite sandwich temperature below the limit of 300 F.

Ko, William L.↗
NASA NTRSDate not supplied

MEDLI2: MISP Inferred Aerothermal Environment and Flow Transition Assessment

The Mars Entry, Descent, and Landing Instrumentation 2 (MEDLI2) sensor suite on the Mars2020 mission contained multiple sensors on the aeroshell to measure the aerothermal environment during entry into the Martian atmosphere. These sensors performed superbly and successfully returned forebody and aftbody heating measurements. Analysis of MEDLI2 data indicated flow transitioning from a laminar to turbulent state on the heatshield. No evidence of flow transition was observed on the backshell. Two methods were used to estimate flow transition times on the heatshield: (1) temperature gradient of near-surface thermocouple data and (2) heat flux gradient from an inverse reconstruction approach using thermocouple data and material response modeling. Both methods produced similar transition times with an estimated accuracy of ±1 s. To assess various transition criteria, transition parameters were evaluated at each sensor location using flow field solutions from computational fluid dynamics (CFD) simulations. The idea was to use conservative values inferred from MEDLI2 data as transition criteria for other Mars missions. To test this hypothesis, MEDLI data from the Mars Science Laboratory (MSL) mission was used to compare predicted vs. actual flow transition times. The comparisons suggest smooth wall transition criteria are not well-suited in modeling the rapid progression of a turbulent transition front. Transition criteria containing a roughness element parameter agreed better with the flight data. In summary, critical transition values derived from MEDLI2 data may be used as a starting point in constructing a flow transition model for future Mars missions.

Chun Y Tang↗
DOE OSTI2024

Identifying Topological Defects in Lamellar Phases through Contour Analysis of Complex Wave Fields

Lamellar phases frequently contain structural imperfections that significantly affect their behaviors and properties. Our previous research successfully reconstructed real-space configurations of defective lamellar phases from diffuse scattering patterns, indicating the presence of phase vortices as a potential method for identifying topological defects disrupting the smectic ordering. Here, this report presents a mathematical framework using regularized wave fields to represent defective lamellar structures in real space. Phase singularities, resulting from the interference of random waves and indicating lamellar order disruption, are identified through a contour integral. These wave fields, derived from coherent scattering in reciprocal space, were validated via computational benchmarks analyzing small-angle neutron scattering data from AOT surfactant solutions, facilitating further statistical analysis of the defects. Our study highlights the potential to extract meaningful information about topological defects in lyotropic phases by inversely analyzing experimentally measured two-point static correlations. Our method allows for detailed structural analysis of various lyotropic phases, both particulate and nonparticulate, in their quiescent states and facilitates quantitative investigation of defects’ role in phase transitions. By integrating small-angle scattering, deep learning, and vortex tangle analysis, our comprehensive approach shows promise in addressing complex challenges in the structural analysis of soft matter systems.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Electronic structure and optical properties of halide double perovskites from a Wannier-localized optimally-tuned screened range-separated hybrid functional

Halide double perovskites are a chemically diverse and growing class of compound semiconductors that are promising for optoelectronic applications. However, the prediction of their fundamental gaps and optical properties with density functional theory (DFT) and ab initio many-body perturbation theory has been a significant challenge. Recently, a nonempirical Wannier-localized optimally tuned screened range-separated hybrid (WOT-SRSH) functional has been shown to accurately produce the fundamental band gaps of a wide set of semiconductors and insulators, including lead halide perovskites. Here, in this study, we apply the WOT-SRSH functional to five halide double perovskites and compare the results with those obtained from other known functionals and previous GW calculations. We also use the approach as a starting point for GW calculations and we compute the band structures and optical absorption spectrum for Cs 2 AgBiBr 6 , using both time-dependent DFT and the GW-Bethe-Salpeter equation approach. We show that the WOT-SRSH functional leads to accurate fundamental and optical band gaps, as well as optical absorption spectra, consistent with spectroscopic measurements, thereby establishing WOT-SRSH as a viable method for the accurate prediction of optoelectronic properties of halide double perovskites.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Growth of Carbon Nanotubes Using Mixed-Metal Catalysts That Include Heavy Refractory Metals as Catalyst Stabilizers

Here, we report new findings from studies of a novel technique to augment the growth of CNTs. In this method, the catalytic metals used for CNT growth are combined with high-melting-point heavy refractory metals, which serve to stabilize the metal catalyst nanoparticles, from which the CNTs nucleate and grow. Our study investigates the combination of an iron catalyst with tungsten and osmium refractory stabilizers. We find that both tungsten and osmium enhance the lifetime of the iron catalyst particles by a factor of approximately two. Tungsten stabilizer resulted in production of CNTs longer by a factor of 1.25 compared to those produced using a pure iron catalyst, while the osmium stabilizer did not significantly affect the ultimate length of the CNTs due to the decreased CNT growth rate with a combined Fe/Os catalyst. These results suggest a possible new route toward producing CNTs with lengths sufficient for materials applications.

Park, Wonjung [Univ. of Tampa, FL (United States)]↗
DOE OSTI2026

From molecular to macroscopic: predicting liquid–liquid phase equilibria and small-angle scattering of mixtures of organic liquids from atomistic simulation using Kirkwood–Buff theory

Macroscopic phase equilibria between solutions define the functionality of many biological and industrial processes, yet they are challenging to predict due to the inherent complexity of liquids containing large molecules. This work introduces an approach for the purely predictive calculation of such phase equilibria in temperature-composition space from molecular dynamics (MD) simulations at one temperature in the single-phase region. We use an approach developed previously to obtain the entropic and enthalpic contributions to the free energy of mixing from the atomic-scale information given by MD simulations via Kirkwood–Buff theory. This allows us to accurately estimate the free energy of mixing as a function of temperature, and thus obtain liquid–liquid phase equilibria, including liquid–liquid critical points, associated binodal and spinodal lines, and composition fluctuations across a region of temperature and composition. Results for binary malonamide–alkane systems are validated by comparison to a direct experimental probe of the fluctuations: the small angle X-ray scattering intensity near zero wavenumber. The MDKB → Phase method demonstrated here provides a significant improvement in predicting liquid–liquid equilibria and free energy as a function of temperature for our systems of interest compared to conventional thermodynamic models. The accurate performance of this purely predictive approach lies in its preservation of atomistic details when determining thermodynamic properties. Furthermore, its inherent extensibility to multi-component systems will likely make the MDKB → Phase approach a valuable general tool for connecting molecular interactions to macroscopic phase equilibria and for the computational screening of materials for targeted thermodynamic behavior.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Development of Bismuth and Platinum Bi-Metallic Nanoparticles to Enhance Melt Wire Temperature Resolution

Advanced manufacturing (AM) based on direct-write (DW) technologies has emerged as the predominant enabler for the fabrication of active and passive sensors to be deployed in harsh operating environments seen in a nuclear reactor. Recently, Idaho National Laboratory and Boise State University have recently established capabilities to incorporate advanced manufacturing (AM) methods to accelerate, modernize and enhance functionality of nuclear sensors and instrumentation to achieve the goal of enhancing the safety and efficiency of nuclear reactors. A significant thrust of this work includes the development of nuclear relevant feedstock materials compatible with a variety of direct-write processes for the development, fabrication and testing of AM sensors for peak temperature detection and neutron flux monitoring. For this report, bismuth and bismuth/platinum (BiPt) bi-metallic nanoparticles were synthesized using wet chemical approaches to develop bi-metallic feedstock materials to enhance the sensitivity of melt wires for peak temperature detection. Nanoparticle characterization was performed with Transmission Electron Microscopy (TEM) for nanoparticle composition and size, Transmission X-Ray Diffraction Fluorescence (TXRF) for nanoparticle composition, and Differential Scanning Calorimetry to elucidate the melting point of BiPt bi-metallic nanoparticles. Preliminary results indicate bi-metallic BiPt nanoparticles as a viable pathway for fabricating high resolution AM melt wires.

36 MATERIALS SCIENCE↗
DOE OSTI2025

An unwanted guest in the electrochemical oxidation of high-voltage Li-ion battery electrolytes: the life of highly reactive protons

Lithium-ion batteries (LIBs) are central to the urgent societal need to decarbonize both transportation and energy storage on the grid. Unfortunately, despite their attractive energy/power density, as well as high coulombic and energy efficiencies, further improvement of this technology – especially their durability – is desperately needed. To support these efforts, our study focuses on fundamental understanding of the decomposition pathways for LIB electrolytes at the cathode–electrolyte interface (CEI), as the nature of these reactions directly controls the extent to which cell capacity and voltage decays in these systems. In this study, we employ electrochemical methods, coupled with product analysis using NMR spectroscopy and mass spectrometry, to determine the decomposition mechanisms in both model and technologically relevant electrolytes. Remarkably, we discovered the electrochemical formation of protons with high chemical activity, comparable to known superacids, at potentials relevant to practical Li-ion batteries. Their reactivity toward every individual component of the CEI provides a unified thermochemical origin for a myriad of side reactions that are commonly associated with the electrochemical reaction. In particular, electrochemically generated protons react with intact EC molecules to form CO 2 and other short and long chain ethers. They also undergo an acid–base reaction with LiPF 6 , to form the weaker acid HF, and with the cathode active material, leaching transition metals into the electrolyte. Collectively, the results of this study all point to the urgent need to either mitigate this proton formation or introduce benign harvesting additives via new electrolyte design strategies.

Ilic, Stefan [Argonne National Laboratory (ANL), A↗
DOE OSTI2021

Bayesian automated weighting of aggregated DFT, MD, and experimental data for candidate thermodynamic models of aluminum with uncertainty quantification

Atomic-scale modeling methods such as density functional theory (DFT) and molecular dynamics (MD) can predict the thermodynamic properties of materials at a lower cost than experimental measurements. However, their regular usage in thermodynamic model construction is hampered by the lack of quantitative agreement with experimental measurements and the lack of uncertainty estimates on the data. To make regular usage of this atomistic simulation data, it is important to assess whether the atomistic simulation datasets, by themselves or in combination with experimental measurements, result in the same physics-informed models best supported by experimental measurements alone. Here, models of aluminum thermodynamic properties are discussed using three data sources: atomistic calculations (DFT and MD), experiments, and a combination of atomistic calculations and experiments. The study shows that, after ensuring self-consistency in predicting key invariant points, both experimental measurements and atomistic calculations can significantly contribute to an optimal model.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Automated tabletop exfoliation and identification of monolayer graphene flakes

Over the past two decades, graphene has been intensively studied because of its remarkable mechanical, optical, and electronic properties. Initial studies were enabled by manual “Scotch Tape” exfoliation; nearly two decades later, this method is still widely used to obtain chemically pristine flakes of graphene and other 2D van der Waals materials. Unfortunately, the yield of large, pristine flakes with uniform thickness is inconsistent. Thus, significant time and effort are required to exfoliate and locate flakes suitable for fabricating multilayer van der Waals heterostructures. Here, we describe a relatively affordable tabletop device (the “eXfoliator”) that can reproducibly control key parameters and largely automate the exfoliation process. In a typical exfoliation run, the eXfoliator produces 3 or more large (≥ 400 μm 2 ) high-quality graphene monolayer flakes, allowing new users to produce such flakes at a rate comparable to manual exfoliation by an experienced user. Furthermore, we use an automated mapping system and a computer vision algorithm to locate candidate flakes. Our results provide a starting point for future research efforts to identify more precisely which parameters matter for the success of exfoliation and to optimize them.

47 OTHER INSTRUMENTATION↗
NASA NTRS1992

Lunar oxygen - The reduction of glass by hydrogen

The direct reduction of volcanic glass by hydrogen has been proposed as a method of extracting oxygen from the lunar soil. Experiments using lunar simulant glasses reacted with flowing hydrogen gas have demonstrated reduction at temperatures from 1000 to 1200 C. For melted samples ferrous iron was reduced to the metal, which formed large crystals at the expense of the glass. Samples held below the melting point rapidly devitrified, and iron was formed from submicrometer crystals of ilmenite and pyroxene. Weight losses of 3.6 - 4.5 percent, depending on glass composition, were achieved in 3 hours at 1100 C. A lunar oxygen plant operating at this efficiency and utilizing Apollo 17 orange glass as a feedstock could produce 50 kg of oxygen per ton of soil. The processes of reduction and sintering of lunar soil are synergistic, and could be combined to produce both oxygen and construction material at a moon base.

Allen, Carlton C.↗
DOE OSTI2020

Composite Structured Surfaces for Durable Dropwise Condensation

Dropwise condensation of steam on hydrophobic substrates has a 10X higher heat transfer coefficient compared to filmwise condensation. To promote dropwise condensation, low surface energy hydrophobic coatings (polymers) are typically utilized. The low intrinsic thermal conductivity (k < 1 W/(m · K)) of polymers, coupled with high heat transfer coefficient of dropwise condensation (100 kW/(m 2 · K)), necessitates that the coating be thin (< 1µm) in order to avoid reducing the overall heat exchanger conductance. However, thin polymeric films easily degrade. The two opposing requirements result in the need for optimization between the durability (thick coating) and the heat transfer (thin coating). To enable high thermal conductivity in thicker coatings, we develop metal-polymer structured surfaces. By using porous structures as inter-connected heat-conducting backbones that are filled with hydrophobic materials, we enable tuning of the coating effective thermal conductivity and surface energy. Three metal structures were studied; micro/nanowires, inverse opals, and sintered spheres. Heat transfer performance was calculated using three-dimensional finite element method simulations with two distinct boundary conditions; convection at the walls and isothermal walls. Interestingly, the overall conductance shows up to 40% difference depending on the boundary condition used in calculating the composite coating effective thermal conductivity. We use our model to predict the heat transfer performance as a function of metal fraction by volume and by surface area for condensation. By coupling our thermal simulations with a previously verified analytical model for predicting wetting behavior on heterogeneous surfaces, we propose a regime map to predict dropwise-to-filmwise transition. Furthermore, our work not only forms a starting point for the development of durable dropwise condensing surfaces, it identifies important considerations needed for computing effective thermal conductivity of composites.

42 ENGINEERING↗
DOE OSTI2021

Identification of Known and Novel Monomers for Poly(hydroxyurethanes) from Biobased Materials

Sourcing polymers from biobased materials is desirable for a transition to a more sustainable world. However, finding new monomers sourced from biobased molecules requires extensive experimental effort, due to the fact that most biomolecules lack the functional groups or properties needed to serve as a monomer in a particular polymeric system without chemical modification. In this work, we demonstrate a computational screening method to obtain monomers of a non-isocyanate polyurethane (NIPU), polyhydroxyurethane (PHU), from biomolecules. Our method generates candidate monomers from a list of 15 selected biobased starting molecules by performing chemical transformations that are standard to chemical catalysis and organic synthesis on these starting molecules and then screening the results for molecules with moieties that can be used as monomers for a PHU. Product molecules from two generations of reaction, a total number of 20,057, were analyzed automatically to identify molecules that have two or more epoxy or cyclic carbonate functional groups as candidates for further derivatization. These 21 candidate molecules were then subjected to another two generations of reaction, producing a set of 17,913 additional molecules. Further scrutiny of this set for desired functional groups characterizing PHUs and querying against the PubChem database was carried out, and 15 known molecules that are commercially available were identified, with the desirable features that they are predicted to be derived from biobased starting molecules in four or fewer steps. Moreover, molecule space generated from this process was primarily comprised of molecules that were not known to PubChem, indicating the framework also has the ability to generate novel candidate monomers that originate from biobased starting points.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Molecular Rotations, Multiscale Order, Hyperuniformity, and Signatures of Metastability during the Compression/Decompression Cycles of Amorphous Ices

We model, via large-scale molecular dynamics simulations, the isothermal compression of low-density amorphous ice (LDA) to generate high-density amorphous ice (HDA) and the corresponding decompression extending to negative pressures to recover the low-density amorphous phase (LDA HDA ). Both LDA and HDA are nearly hyperuniform and are characterized by a dynamical HBN, showing that amorphous ices are nonstatic materials and implying that nearly hyperuniformity can be accommodated in dynamical networks. In correspondence with both the LDA-to-HDA and the HDA-to-LDA HDA phase transitions, the (partial) activation of rotational degrees of freedom activates a cascade effect that induces a drastic change in the connectivity and a pervasive reorganization of the HBN topology which, ultimately, break the samples’ hyperuniform character. Key to this effect is the rapid rate at which changes occur, and not their magnitude. The inspection of structural properties from the short- to the long-range shows that signatures of metastability are present at all length-scales, hence providing further solid evidence in support of the liquid–liquid critical point scenario. LDA and LDA HDA differ in terms of HBN and structural properties, implying that they are distinct low-density glasses. Our work unveils the role of molecular rotations in the phase transitions between amorphous ices and shows how the unfreezing of rotational degrees of freedom generates a cascade effect that propagates over multiple length-scales. Our findings greatly improve our basic understanding of water and amorphous ices and can potentially impact the field of molecular network-forming materials at large.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

A probabilistic-phase field model for the fracture of brittle materials

Abstract We develop a computational method to determine the failure probability of brittle materials under general mechanical loading conditions. The method is a combination of two parts: (1) numerical simulations of materials with multiple cracks using phase field theory, where the complete fracture process is viewed as ‘damage percolation’ along critical paths or clusters of cracks, rather than the traditional weak-link failure mechanism of Weibull, and (2) an extension of the Batdorf statistical theory of fracture to finite domains, where it is implemented within the finite element framework. The results of phase-field simulations at the ‘percolation threshold’ are used as failure data in the Batdorf theory to determine the overall probability of failure. The input to this approach is the size distribution of cracks in a pristine material. An example is shown, where alumina samples that were previously tested by Abe and coworkers (Abeet al2003J. Am. Ceram. Soc.861019–21) in four-point loading are compared to the results of our numerical simulations. The approach developed here has the advantage of being extendable to more complex thermomechanical loading.

Materials Science↗