Search NASA⌕ Search

SEARCH · Search NASA

Results for “3D Integration”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

DECOVALEX-2023: Task C Final Report

The Full-scale Emplacement (FE) heater experiment at the Mont Terri Underground Rock Laboratory (URL) was designed and conducted by Nagra to replicate an emplacement tunnel of Nagra’s reference repository design at 1:1 scale. Alongside testing the technical feasibility of constructing disposal tunnels, emplacing waste containers in the tunnels and then backfilling them, the main goals of the FE experiment are (1) to obtain a better understanding of the coupled effects of induced thermo-hydro-mechanical (THM) processes that may occur and (2) to validate existing coupled THM models (Müller et al., 2017). A key aspect of ensuring safety for repositories located in low-permeability rock involves minimizing any damage to the rock itself, thereby preserving its integrity and promoting a stable environment Amongst a number of processes that could damage the rock is the increase in pore pressure due to thermal loading caused by heat emitted from the waste. To reduce the potential damage of the rock, it is important to analyse the evolution of heat over time due to the heat load of the containers and assess possible consequences by coupled THM models. The aim of Task C of DECOVALEX-2023 was to build 3D numerical models of the FE experiment, focussing in particular on the heating induced pore pressure change in the Opalinus Clay. Data from a large number of sensors were available from the FE experiment for model comparison. These sensors measured temperature and relative humidity in the bentonite around the heaters, and temperature, pressure and displacement/strain in the surrounding Opalinus clay. Data were available from the start of excavation (April 2012) up to August 2020 for most sensors (more than 5 years from the start of heating in December 2014). To fulfil the overall aim of the task, the work was broken down into a number of steps, starting with simpler models to build confidence in each team’s approach and then moving to more complex models that better represent the FE experiment. Step 0 consisted of 2D benchmark models, gradually increasing the number of processes that are represented from thermal (T) only models in Step 0a, to coupled thermal hydraulic (TH) models in Step 0b with a representation of changing porosity, to coupled thermo-hydro-mechanical (THM) models in Step 0c, where porosity changes are calculated by the mechanical model. A detailed specification of processes, parameters, initial and boundary conditions was provided for this step, with the ambition that all teams would work towards close agreement in their model results, thus building confidence in the model implementations. vi It was not straightforward to achieve agreement between the teams, so additional steps (Step 0b2, 0b3, 0c2, 0c3) were added along with derivation of some analytical solutions against which the models could be compared. The reasons for the differences between teams were investigated and found to be caused primarily by different conceptual model assumptions (including temperature dependence of the thermal expansion of water), different model formulations (including porosity evolution) and differences in modelled domain sizes, boundary conditions and grid discretisation. This demonstrates that comparisons between multiple modelling teams and/or comparison with analytical results and experimental data are highly beneficial in providing an indication of uncertainty in model predictions. At the conclusion of Step 0, almost all teams had achieved a close agreement in model results and those that had not achieved an agreement knew the reason for this. Step 1 moved from 2D models to 3D models of the FE experiment without adding technical features like shotcrete or EDZ, and only considering the heating phase. Initially the 3D model was tightly specified to continue to build confidence in the model implementations (Step 1a). The results of Step 1a were compared to the data from the FE-experiment without the teams seeing the data. The teams were then provided with a sub-set of the data from the FE-experiment and invited to consider how best to use the large dataset for model comparison (Step 1b). Teams were then asked to use the data provided to calibrate their models, only changing material property values rather than adding features or processes to their models (Step 1c). In Step 1, teams were asked to only model the heating phase of the experiment, so pressure in the Opalinus Clay was reported as change in pressure since the initial conditions were specified rather than modelled. The change from 2D to 3D models was accompanied by an increase in the dispersion of results between the teams. Some of this was resolved during the task, but some remained and is potentially due to model discretisation. Calibration of parameters was useful in improving the fit of the models to the data but the remaining differences indicated that the models were missing features or processes. In Step 2, the teams were asked to update their models with additional features and processes as well as calibrating parameters to try and improve the fit of the models to the data. Teams were encouraged to represent ventilation of the open FE tunnel prior to backfilling with heaters and bentonite and in Step 2, the absolute pressure in the Opalinus Clay was compared between the teams. Teams took different approaches, but there was consideration of adding shotcrete and an EDZ into the model, representing stress change during excavation and different approaches to modelling ventilation of the FE tunnel. Overall, the documented results showed a very good agreement for temperature. The results for porewater pressure evolution showed a significant improvement for most teams compared to Step 1c with a good agreement to the measurements for several teams whereas some teams overpredicted the pressure increase and others overpredicted the drainage effect especially for the sensors close to the heater. Step 3 was an opportunity for teams to use the models developed in Step 1 and Step 2 to make predictions about the temperature and pressure changes that will be expected at the FE experiment over the next few years in light of the planned changes in thermal output of the heaters.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

DeepLynx Ecosystem 2025

Poor data integration and governance continue to plague complex engineering projects, resulting in missed cost, schedule, and performance targets. Departments operate in isolated systems with manual data exchange, creating fragmented information that compounds errors and leads to significant delays and cost overruns. The DeepLynx ecosystem addresses these challenges through an open-source, modular data management platform that transforms fragmented project data into an integrated digital thread. Built on a federated microservice architecture, the ecosystem comprises seven specialized tools centered around DeepLynx Nexus, a unified data catalog with hierarchical organization and graph-based navigation capabilities. The ecosystem includes: DeepLynx Stream for real-time timeseries data ingestion from industrial sources; DeepLynx Ingest for governed data uploads with formal review workflows; DeepLynx Lattice for ontology-based entity and relationship extraction; DeepLynx Run for workflow orchestration and secure AI/ML compute; DeepLynx Visualize for 3D digital twin visualization; and DeepLynx Insight for AI-assisted document analysis with traceable, grounded responses. Deployable in cloud, on-premise, or hybrid environments using containerized Docker applications and Helm charts, the DeepLynx ecosystem provides flexible infrastructure that adapts to organizational requirements. By consolidating project data into a unified data lake with role-based access controls and OAuth2 authentication, DeepLynx enables digital thread and digital twin capabilities that improve decision-making, reduce risk, and support complex engineering workflows throughout the project lifecycle.

42 - ENGINEERING↗
DOE OSTI2026

Physicochemical properties of digital light processing 3D-Printed alumina and mullite ceramics

Alumina ceramics fabricated using conventional techniques such as uniaxial pressing or injection molding are popular due to their low density, excellent insulation, and mechanical properties. However, these methods often limit the fabrication of complex geometries with high dimensional accuracy due to tooling constraints and limited design freedom. To overcome these limitations, this study employed Digital Light Processing (DLP) additive manufacturing (AM), which enables the production of precise structures. In addition, the optimization of sintering parameters to enhance the densification and performance of alumina and mullite ceramics was investigated, with a specific focus on how varying sintering temperatures and hold times affect part shrinkage, geometric accuracy, and material integrity. Flexural strength of both alumina and mullite specimens was clearly influenced by the way the layers was stacked. When layers were arranged across the direction of the applied load (Z = 4; XZ), the strength was higher than when they were stacked along the same direction as the load (Z = 3; XY). Weibull analysis based on these results showed high modulus values across all samples, indicating good reliability in the flexural strength measurements. This reliability, combined with the clear influence of sintering parameters, highlights how processing conditions effect the material properties and mechanical performance of parts produced through DLP additive manufacturing. In conclusion, these findings open new avenues for ceramic AM across various applications, enhancing the potential for innovation in fields such as aerospace, biomedical engineering, and energy.

36 MATERIALS SCIENCE↗
DOE OSTI2023

3P Program: Phenotyping X Prediction = Productivity (Final Scientific/Technical Report)

The goal of the 3P Program was to establish integrated, real-time phenotyping and to analyze above- and below-ground plant architecture and total carbon partitioning and allocation to predict heterosis and develop superior crop hybrids by fully leveraging the Sorghum gene pool. There were two overarching themes: 1) the development of a new crop improvement approach utilizing advances in high-throughput phenotyping (HTP), computing, and genomics for public dissemination and 2) leveraging this platform for sorghum crop improvement and commercialization. The Clemson team worked on creating genomic resources and using both statistical learning and high-throughput phenotyping in genomics-assisted breeding. Research was broadly interested in the genetics of carbon partitioning, with the aim of improving crop performance and achieving sustainability. The technology and resources created can be readily found in the public domain and serve to advance scientific understanding of crop genomics and breeding. Genomic prediction was able to identify top crosses to be made, and a hybrid prediction pipeline is in place to drive year-over-year genetic gain. Roots have long been ignored by plant breeders and agronomists, not because they are unimportant but because they are hard to measure. This is an untapped white space of potential insight and innovation. To address this, Hi Fidelity Genetics developed the RootTracker to measure roots in the field on a continuous basis. A database system called RootTracker Tracker was developed to handle data coming from the RootTrackers. In using this device, valuable data was observed for plant breeding, hydrochemical development, and other agricultural biology applications. Carnegie Mellon’s goal was developing new techniques to generate high-resolution 3D models of plants from data collected in the field. The idea was that more useful and more informative phenotypes could be extracted by resolving small features, such as seeds and flowers, and that by modeling in 3D, the spatial structure of plants could be examined. To achieve this, multiple images collected by a new small format structured light stereo imager were fused together. A sorghum panicle modeling pipeline was developed to allow the collection and processing of data. Carolina Seed Systems is an agricultural technology company focused on decarbonizing the agricultural system. Their technology pipeline serves to drive fundamental progress towards creation and distribution of carbon negative crops. The genomic and the engineering technology developed through the 3P Program was leveraged to deliver both value and sustainability from the grower to the consumer. Promising sorghum hybrids were scaled up and commercialized. The overall goal of our research was to integrate, create, and deploy genetic and engineering concepts and technologies to enhance crop productivity in a sustainable fashion. The combination of public and private partners allowed the basic research and hypothesis testing to be quickly accelerated for commercial application by the companies yet maintained that the core framework and academic insights remain in the public domain for continued market disruption, competition, and innovation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Toward Lithium Recovery Using Modular (and Membraneless) Phase Separation and Extraction (MPSE) Technology with Ionic Liquid (IL) Solvents: Effect of Coatings

In this work, a single-channel slope-plate Modular and Membrane-less Phase Separation and Extraction (MPSE) device was fabricated using 3D printing. The slope-plate was designed with an insertable glass slide to enable surface modification with coatings. Self-assembled monolayers (SAMs) and perfluoropolyether (Zdol) coatings were applied to tailor surface wettability and enhance phase separation. Using a model biphasic system of water and hexadecane, both coatings significantly improved the separation efficiency compared to the uncoated device by promoting selective wetting. Building on these results, lithium extraction from a simulated saline solution was investigated using 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([EMIM][NTf₂]) as the ionic liquid (IL) extractant. The IL effectively extracted Li⁺ from the aqueous phase in bulk extraction, demonstrating its potential for lithium recovery. However, the low interfacial tension between the IL and aqueous phases posed challenges for phase separation. The application of SAM and Zdol coatings effectively mitigated this issue. Overall, integrating tailored surface chemistry with the slope-plate MPSE design shows great promise as an efficient and scalable platform for studying liquid–liquid separation and optimizing ionic liquid–based extraction processes for lithium recovery from saline sources.

3D printing↗
DOE OSTI2025

Chemistry and Interfacial Structure Promoting Quasi-van der Waals Epitaxial Growth of WS 2 Nanosheets on Sapphire for Prospective Application in Field-Effect Transistors

How do chemical and structural modifications to the supporting crystal surface affect the subsequent van der Waals (vdW) or quasi(Q)-vdW epitaxial growth of 2D nanocrystals? Developing an atomic-scale picture of such an interfacial system is crucial for understanding its impact on the physical and chemical properties of the supported 2D materials. The elucidation of the interfacial structure and chemistry needed to promote the Q-vdW epitaxial growth of 2D tungsten disulfide (WS 2 ) nanocrystals contributes to the growth mechanism understanding, thus pushing forward the integration of such atomically thin semiconductors toward real field-effect transistor applications. In addition to an atomic-force microscopy top view, we showcase a combination of X-ray techniques for a top-to-bottom investigation of the complexities of the buried interface structures. Furthermore, this approach uses X-ray photoelectron spectroscopy, X-ray standing wave excited X-ray fluorescence, and crystal truncation rod scattering to produce a highly resolved chemical-state-specific 3D atomic map for the extended interface structure of WS 2 /α-Al 2 O 3 (001). Employing these detailed analysis methods, along with density functional theory to further refine the picoscale structure, we demonstrate how two different types of interface engineering during the pregrowth stage lead to significant differences in the chemical and structural modifications to the terminal surface of c-face sapphire, which in turn leads to substantial differences in the submonolayer growth of supported WS 2 2D nanocrystals in terms of lateral domain sizes, epitaxial registry, vdW gaps, and stability.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Integrated edge-to-exascale workflow for real-time steering in neutron scattering experiments

We introduce a computational framework that integrates artificial intelligence (AI), machine learning, and high-performance computing to enable real-time steering of neutron scattering experiments using an edge-to-exascale workflow. Focusing on time-of-flight neutron event data at the Spallation Neutron Source, our approach combines temporal processing of four-dimensional neutron event data with predictive modeling for multidimensional crystallography. At the core of this workflow is the Temporal Fusion Transformer model, which provides voxel-level precision in predicting 3D neutron scattering patterns. The system incorporates edge computing for rapid data preprocessing and exascale computing via the Frontier supercomputer for large-scale AI model training, enabling adaptive, data-driven decisions during experiments. This framework optimizes neutron beam time, improves experimental accuracy, and lays the foundation for automation in neutron scattering. Although real-time experiment steering is still in the proof-of-concept stage, the demonstrated potential of this system offers a substantial reduction in data processing time from hours to minutes via distributed training, and significant improvements in model accuracy, setting the stage for widespread adoption across neutron scattering facilities and more efficient exploration of complex material systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Status of the OECD/NEA Watts Bar unit 1 benchmark and calculations of local reactor core power data of the Zero Power physics tests

The paper aims to provide an update on the status of the Organization for Economic Cooperation and Development (OECD) / the Nuclear Energy Agency (NEA) Tennessee Valley Authority (TVA) Watts Bar 1 (WB1) Multi-Physics Multi-Cycle depletion benchmark and the results associated with Exercise 1. The benchmark relies on the set of benchmark progression problems developed by the Department of Energy (DOE) Consortium for the Advanced Simulation of Light Water Reactors (CASL) for the Virtual Environment for Reactor (VERA), which are based on real plant design and operational data. The OECD/NEA TVA WB1 benchmark is designed for validation of both traditional and novel high-fidelity multi-physics codes to analyze Pressurized Water Reactors (PWR) depletion cycles. The activities are conducted under the Expert Group on Reactor Systems Multi-Physics (EGMUP) at NEA/OECD. In this work, we analyze Exercise 1 of the benchmark: stand-alone Three-Dimensional (3D) neutronics at Start-up Zero Power Physics Test (ZPPT) at Hot Zero Power Conditions (HZP) and the resulting power maps. The code used to model the exercise is the continuous energy Monte Carlo code Serpent 2.1.31 along with the nuclear data library ENDF/B-VII.1. Serpent results are compared with the high-fidelity deterministic code MPACT, which is part of VERA. The paper presents three selected results from power calculations required output. The results compared are the normalized axially integrated radial core power maps, the normalized axial averaged core power shapes, and the normalized core hottest and coldest assemblies radial power maps. The Root Mean Square Deviation (RMSD) between Serpent and MPACT is 0.51 % for the axially integrated radial core power map, 1.41 % for the axial averaged core power shape, 2.42 % for the hottest assembly axial power shape, 0.88 % for the coldest assembly axial integrated hottest assembly radial pin power map, and 0.17 % for the axially integrated coldest assembly radial pin power map. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

SMATool: Strength of materials analysis toolkit

The study of the strength of materials is a cornerstone in material science and engineering, playing a critical role in shaping the progress and application of materials in diverse industrial sectors. The strength of a material is meticulously examined to understand the behavior of the material under different stress conditions and environments, thereby guiding material selection and structural design. Herein, we introduce the SMATool, a computational toolkit for the efficient calculation and analysis of material strength at both zero and finite temperatures for 3D, 2D, 1D, and tubular 2D-based nanostructures and nanotubes, as well as 1D nanoribbons. The toolkit is capable of calculating tensile, shear, ultimate, yield, and indentation (Vickers' hardness) strengths in various dimensions, as well as the energy storage capacity. We conducted several calculations both at zero and finite temperatures to validate the accuracy and reliability of the developed software. Here, the results show that the SMATool package provides accurate predictions that align with existing data on material strength. SMATool integrates seamlessly with widely used electronic structure codes like VASP and Quantum Espresso, providing a user-friendly interface catering to academic researchers and industry professionals.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Interpretative 3D MHD modelling of deuterium SPI into a JET H-mode plasma

Abstract The pre-thermal quench (pre-TQ) dynamics of a pure deuterium ( D 2 ) shattered pellet injection (SPI) into a 3 MA / 7 MJ JET H-mode plasma is studied via 3D non-linear MHD modelling with the JOREK code. The interpretative modelling captures the overall evolution of the measured density and radiated power. The simulations also identify the importance of the drifts of ablation plasmoids towards the tokamak low field side (LFS) and the impurities in the background plasma in fragment penetration, assimilation, radiative cooling and MHD activity in D 2 SPI experiments. It is found that plasmoid drifts lead to an about 70% reduction of the central line-integrated density (compared to a simulation without drifts) in the JET D 2 SPI discharge considered. Impurities that pre-exist before the SPI as well as those from possible impurity influxes related to the SPI are shown to dominate the radiation in the considered discharge. With inputs from JOREK simulations, modelling with the Lagrangian particle-based pellet code PELOTON reproduces the deviation of the SPI fragments in the direction of the major radius as observed by the fast camera. This confirms the role of rocket effects and plasmoid drifts in the considered discharge and reinforces the validity of the JOREK modelling. The limited core density rise due to plasmoid drifts and the strong radiative cooling and MHD activity with impurities (depending on their species and concentration) could limit the effectiveness of LFS D 2 SPI in runaway electron avoidance and are worth considering in the design of the ITER disruption mitigation system.

Physics↗
DOE OSTI2021

Time-Lapse Integration at FWU: Fluids, Rock Physics, Numerical Model Integration, and Field Data Comparison

We present the current status of time-lapse seismic integration at the Farnsworth (FWU) CO2 WAG (water-alternating-gas) EOR (Enhanced Oil Recovery) project at Ochiltree County, northwest Texas. As a potential carbon sequestration mechanism, CO2 WAG projects will be subject to some degree of monitoring and verification, either as a regulatory requirement or to qualify for economic incentives. In order to evaluate the viability of time-lapse seismic as a monitoring method the Southwest Partnership (SWP) has conducted time-lapse seismic monitoring at FWU using the 3D Vertical Seismic Profiling (VSP) method. The efficacy of seismic time-lapse depends on a number of key factors, which vary widely from one application to another. Most important among these are the thermophysical properties of the original fluid in place and the displacing fluid, followed by the petrophysical properties of the rock matrix, which together determine the effective elastic properties of the rock fluid system. We present systematic analysis of fluid thermodynamics and resulting thermophysical properties, petrophysics and rock frame elastic properties, and elastic property modeling through fluid substitution using data collected at FWU. These analyses will be framed in realistic scenarios presented by the FWU CO2 WAG development. The resulting fluid/rock physics models will be applied to output from the calibrated FWU compositional reservoir simulation model to forward model the time-lapse seismic response. Modeled results are compared with field time-lapse seismic measurements and strategies for numerical model feedback/update are discussed. While mechanical effects are neglected in the work presented here, complementary parallel studies are underway in which laboratory measurements are introduced to introduce stress dependence of matrix elastic moduli.

4D↗
DOE OSTI2020

A MEMS dynamic mechanical analyzer for in situ viscoelastic characterization of 3D printed nanostructures

Cellular metamaterial structures with sub-micron features have shown the ability to become excellent energy absorbing materials for impact mitigation due to the enhanced mechanical properties of materials at the nanoscale. However, in order to optimize the design of these energy absorbing metamaterial structures we need to be able to measure the dynamic properties of the sub-micron features such as storage and loss moduli and the loss factor. Therefore, at scale testing is required to capture the scale, temperature, and strain rate dependent material properties of these nanoscale materials. This study presents the design, fabrication, and calibration of a MEMS-based dynamic mechanical analyzer (DMA) that can be directly integrated with the two photon lithography (TPL) process commonly used to fabricate metamaterial structures with nanoscale features. The MEMS-based DMA consists of a chevron style thermal actuator used to generate a tensile load on the structure and two differential capacitive sensors on each side of the structure used to measure load and displacement. This design demonstrated 1.5 ± 0.75 nm displacement resolution and 104 ± 52 nN load resolution, respectively. Dynamic mechanical analysis was successfully conducted on a single nanowire feature printed between the load and displacement stages of the MEMS device with testing frequencies ranging between 0.01–10 Hz and testing temperatures ranging between 22 °C–47 °C. These initial tests on an exemplar TPL part demonstrate that the printed nanowire behaves as a viscoelastic material wherein the transition from glassy to viscous behavior has already set in at the room temperature.

42 ENGINEERING↗
DOE OSTI2024

High–Performance NiCo 2 O 4 /Graphene Quantum Dots for Asymmetric and Symmetric Supercapacitors with Enhanced Energy Efficiency

For the sustainable growth of future generations, energy storage technologies like supercapacitors and batteries are becoming more and more common. However, reliable and high-performance materials’ design and development is the key for the widespread adoption of batteries and supercapacitors. Quantum dots with fascinating and unusual properties are expected to revolutionize future technologies. However, while the recent discovery of quantum dots honored with a Nobel prize in Chemistry, their benefits for the tenacious problem of energy are not realized yet. In this context, herein, chemical-composition tuning enabled exceptional performance of NiCo 2 O 4 (NCO)/graphene quantum dots (GQDs) is reported, which outperform the existing similar materials, in supercapacitors. A comprehensive study is performed on the synthesis, characterization, and electrochemical performance evaluation of highly functional NCO/GQDs in supercapacitors delivering enhanced energy efficiency. The high-performance, functional NCO/GQDs electrode materials are synthesized by the incorporation of GQDs into NCO. The effect of variable amount of GQDs on the energy performance characteristics of NCO/GQDs in supercapacitors is studied systematically. In-depth structural and chemical bonding analyses using X-ray diffraction (XRD) and Raman spectroscopic studies indicate that all the NCO/GQDs composites crystallize in the spinel cubic phase of NiCo 2 O 4 while graphene integration evident in all the NCO/GQDs. The scanning electron microscopy imaging analysis reveals homogeneously distributed spherical particles with a size distribution of 5–9 nm validating the formation of QDs. The high-resolution transmission electron microscopy analyses reveal that the NCOQDs are anchored on graphene sheets, which provide a high surface area of 42.27 m 2 g –1 and high mesoporosity for the composition of NCO/GQDs-10%. In addition to establishing reliable electrical connection to graphene sheets, the NCOQDs provide reliable 3D-conductive channels for rapid transport throughout the electrode as well as synergistic effects. Chemical-composition tuning, and optimization yields NCO/GQDs-10% to deliver the best specific capacitance of 3940 Fg –1 at 0.5 Ag –1 , where the electrodes retain ≈98% capacitance after 5000 cycles. The NCO/GQD-10%//AC asymmetric supercapacitor device demonstrates outstanding energy density and power density values of 118.04 Wh kg –1 and 798.76 W kg –1 , respectively. The NCO/GQDs-10%//NCO/GQDs-10% symmetric supercapacitor device delivers excellent energy and power density of 24.30 Wh kg –1 and 500 W kg –1 , respectively. These results demonstrate and conclude that NCO/GQDs are exceptional and prospective candidates for developing next-generation high-performance and sustainable energy storage devices.

25 ENERGY STORAGE↗
DOE OSTI2022

Damage of prismatic lithium‐ion cells subject to bending: Test, model, and detection

Abstract The mechanically induced internal short circuit (ISC) is one of the major safety concerns of lithium‐ion batteries. Mechanical abuse tests are often performed to evaluate the integrity and safety of lithium‐ion batteries under mechanical loadings. Except for the widely explored compression‐dominated indentation tests, bending is another typical real‐world loading condition that is tension‐dominated. To investigate the mechanical damage and ISC behavior of batteries under bending, we carried out controlled three‐point bending tests in four progressive steps on prismatic battery cells with maximum deflections ranging from 38% to 76% of the cell thickness. None of the tested cells experienced an ISC. We then conducted 3D X‐ray computed tomography (CT) scanning on the bent cells after unloading. X‐ray CT images showed three out of the four tested cells have extensive cracking in the electrode layers at the bottom side (opposite to the loading head). This indicates that cracking does not necessarily lead to an ISC under bending. Electrochemical impedance spectroscopy was also measured on the bent cells and substantial changes were observed. Both the bulk resistance and charge‐transfer resistance increased significantly after bending, which could influence the battery performance and lifespan. We then developed a detailed finite (FE) element model to further investigate the mechanical deformation and failure mechanisms. The FE model successfully predicts the load–displacement response and reproduces the deformation patterns. The findings and the FE model developed in the present study provide useful insights and tools for the battery structure and crash safety design. image

25 ENERGY STORAGE↗
DOE OSTI2024

Finite elements for Matérn-type random fields: Uncertainty in computational mechanics and design optimization

This work highlights an approach for incorporating realistic uncertainties into scientific computing workflows based on finite elements, focusing on prevalent applications in computational mechanics and design optimization. We leverage Matérn-type Gaussian random fields (GRFs) generated using the SPDE method to model aleatoric uncertainties, including environmental influences, variating material properties, and geometric ambiguities. Our focus lies on delivering practical GRF realizations that accurately capture imperfections and variations and understanding how they impact the predictions of computational models as well as the shape and topology of optimized designs. Here we describe a numerical algorithm based on solving a generalized SPDE to sample GRFs on arbitrary meshed domains. The algorithm leverages established techniques and integrates seamlessly with the open-source finite element library MFEM and associated scientific computing workflows, like those found in industrial and national laboratory settings. Our solver scales efficiently for large-scale problems and supports various domain types, including surfaces and embedded manifolds. We showcase its versatility through biomechanics and topology optimization applications, emphasizing the potential to influence these domains. The flexibility and efficiency of SPDE-based GRF generation empowers us to run large-scale optimization problems on 2D and 3D domains, including finding optimized designs on embedded surfaces, and to generate design features and topologies beyond the reach of conventional techniques. Moreover, these capabilities allow us to model and quantify geometric uncertainties on reconstructed submanifolds, such as the interpolated surfaces of cerebral aneurysms provided by postprocessing CT scans. In addition to offering benefits in these specific domains, the proposed techniques transcend specific applications and generalize to arbitrary forward and backward problems in uncertainty quantification involving finite elements.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Comparison of Shape Optimization Methods for Heat Exchanger Fins Using Computational Fluid Dynamics

Inverse design techniques are one way to leverage advances in 3D printing, artificial intelligence, and computational resources to achieve increased performance of heat exchangers. Two optimization techniques (genetic algorithm and particle swarm) and three geometry representations (binary level set, composite Bézier, and free form deformation) are used to increase the heat transfer and reduce the pressure drop of a heat exchanger fin. After running 210,810 OpenFOAM simulations, results indicate that a significant performance increase of the fin can be realized in less than 48 hrs, allowing for such a process to be integrated in traditional design processes. The best design increased the performance of the objective function, compared to the baseline rectangular geometry by 75%. A custom distributed infrastructure was built, allowing for all methods to reach 95% of the final objective values in a little over 4 hrs, handling 1674 OpenFoam simulations per hour.

42 ENGINEERING↗
DOE OSTI2021

Computational Fluid Dynamics (CFD) Simulations of Taylor Bubbles in Vertical and Inclined Pipes with Upward and Downward Liquid Flow

Summary Two-phase flow is a common occurrence in pipes of oil and gas developments. Current predictive tools are based on the mechanistic two-fluid model, which requires the use of closure relations to predict integral flow parameters such as liquid holdup (or void fraction) and pressure gradient. However, these closure relations carry the highest uncertainties in the model. In particular, significant discrepancies have been found between experimental data and closure relations for the Taylor bubble velocity in slug flow, which has been determined to strongly affect the mechanistic model predictions (Lizarraga-García 2016). In this work, we study the behavior of Taylor bubbles in vertical and inclined pipes with upward and downward flow using a validated 3D computational fluid dynamics (CFD) approach with level set method implemented in a commercial code. A total of 56 cases are simulated, covering a wide range of fluid properties, pipe diameters, and inclination angles: Eo ∈ [10, 700]; Mo ∈ [1×10–6, 5×103]; ReSL ∈ [–40, 10]; θ ∈ [5°, 90°]. For bubbles in vertical upward flows, the simulated distribution parameter, C0, is successfully compared with an existing model. However, the C0 values of downward and inclined slug flows where the bubble becomes asymmetric are shown to be significantly different from their respective vertical upward flow values, and no current model exists for the fluids simulated here. The main contributions of this work are (1) the relatively large 3D numerical database generated for this type of flow, (2) the study of the asymmetric nature of inclined and some vertical downward slug flows, and (3) the analysis of its impact on the distribution parameter, C0.

Engineering↗
DOE OSTI2024

Additive Manufacturing of Thermal Energy Storage Composites with Microencapsulated Phase Change Materials Supported in a Multi-Polymer Matrix

Advanced manufacturing techniques, such as additive manufacturing (AM), that can directly integrate phase change materials (PCMs) have garnered interest in recent years due to their potential for development of highly efficient thermal energy storage architectures. Complex, high surface area geometries embedded with PCMs that are only feasible with AM can improve thermal management with reduced material waste. Our work focuses on developing composite filaments with microencapsulated phase change materials (MEPCM) bound within a single or dual polymer matrix that can be processed through standard filament extruders and additively manufactured using off-the-shelf 3D printers. Polymer powders, rather than polymer pellets, were key to homogenously mixed filaments achieving high MEPCM loadings with no deterioration in thermal energy storage (TES) capability during extrusion. Composite filaments contain upwards of 60 wt% MEPCM and were printed without loss in feature resolution, print speed, or layer adhesion. Storage enthalpies of printed composites range from 100 - 130 kJ/kg, which were within 5% of the theoretical enthalpy based on weight fraction of MEPCM and maintained enthalpies within 1% over 500 thermal cycles. We can reliably manufacture low density, high surface area structures like 15% gyroid infill, along with dense, compact pucks at a 100% concentric infill. Prints were also scalable to a 900 cm3 honeycomb infill heat exchanger model that has an estimated energy storage capacity of 9 Wh.

3D printing↗