Search NASA⌕ Search

SEARCH · Search NASA

Results for “strain simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Effect of an oxide layer on high velocity impact of tantalum particles characterized using molecular dynamics

Micro-cold spray, also commonly referred to as the aerosol deposition method or vacuum kinetic spraying, is a process for depositing nanostructured films of metals and ceramics by impacting solid nanoparticles at high velocity onto a substrate at room temperature and at low pressures. Tantalum films deposited by micro-cold spray are desirable because tantalum is non-reactive to certain molten metals and can therefore serve as a non-reactive barrier coating. However, unlike most other metals deposited by micro-cold spray, tantalum quickly forms an oxide layer when exposed to air. In the deposition of films using traditional cold spray with much larger particles, native oxides must fracture for particle–substrate bonding to occur. Furthermore, the fraction of the particle that is oxidized in micro-cold spray and the strain rates experienced by the particles upon impact are both much larger than in cold spray, suggesting that the role of the oxide on film deposition may be distinct. The effect of this oxide layer on particle deformation and adherence to a substrate under high strain rates is studied using molecular dynamics simulations. The effect of oxide layer thickness, particle diameter, and impact velocity are examined.

42 ENGINEERING↗
DOE OSTI2025

Finite domain solution of a hydraulic fracture in a permeable rock

In this work, we present a domain-based algorithm to simulate the propagation of a plane-strain hydraulic fracture in a zero-toughness permeable elastic medium. The algorithm utilizes a domain-based method to solve the elasticity equation and integrates a multi-scale tip asymptote, which is particular to hydraulic fractures, into this framework. This integration is key to accurately model the energy dissipation and the fluid leak-off in the fracture tip region. The algorithm combines a 2D finite volume method (FVM) for solving the elasticity equation with a 1D FVM for solving the nonlinear lubrication equation. Incorporating the far-field asymptotics and using a moving-mesh scheme reduces the computational burden while improving the accuracy of the scheme. The paper concludes with an analysis of the numerical results. This study demonstrates the potential of this domain-based approach for modeling hydraulic fractures in poroelastic media.

Domain-based method↗
DOE OSTI2025

Prediction of carbon nanostructure mechanical properties and the role of defects using machine learning

Graphene-based nanostructures hold immense potential as strong and lightweight materials, however, their mechanical properties such as modulus and strength are difficult to fully exploit due to challenges in atomic-scale engineering. This study presents a database of over 2,000 pristine and defective nanoscale CNT bundles and other graphitic assemblies, inspired by microscopy, with associated stress–strain curves from reactive molecular dynamics (MD) simulations using the reactive INTERFACE force field (IFF-R). These 3D structures, containing up to 80,000 atoms, enable detailed analyses of structure-stiffness-failure relationships. By leveraging the database and physics- and chemistry-informed machine learning (ML), accurate predictions of elastic moduli and tensile strength are demonstrated at speeds 1,000 to 10,000 times faster than efficient MD simulations. Hierarchical Graph Neural Networks with Spatial Information (HS-GNNs) are introduced, which integrate chemistry knowledge. HS-GNNs as well as extreme gradient boosted trees (XGBoost) achieve forecasts of mechanical properties of arbitrary carbon nanostructures with only 3 to 6% mean relative error. The reliability equals experimental accuracy and is up to 20 times higher than other ML methods. Predictions maintain 8 to 18% accuracy for large CNT bundles, CNT junctions, and carbon fiber cross-sections outside the training distribution. The physics- and chemistry-informed HS-GNN works remarkably well for data outside the training range while XGBoost works well with limited training data inside the training range. The carbon nanostructure database is designed for integration with multimodal experimental and simulation data, scalable beyond 100 nm size, and extendable to chemically similar compounds and broader property ranges. The ML approaches have potential for applications in structural materials, nanoelectronics, and carbon-based catalysts.

Winetrout, Jordan J.↗
DOE OSTI2021

Anisotropic strength behavior of single-crystal TATB

High-rate strength behavior plays an important role in the shock initiation of high explosives, with plastic deformation serving to localize heat into hot spots and as a mechanochemical means to enhance reactivity. Recent simulations predict that detonation-like shocks produce highly reactive nanoscale shear bands in the layered crystalline explosive TATB (1,3,5-triamino-2,4,6-trinitrobenzene), but the thresholds leading to this response are poorly understood. We utilize molecular dynamics to simulate the high-rate compressive stress–strain response of TATB, with a focus on understanding flow behavior. Here, the dependence of strength on pressure and loading axis (crystal orientation) is explored. The deformation mechanisms fall broadly into two categories, with compression along crystal layers activating a buckling/twinning mode and compression normal to the layers producing nanoscale shear bands. Despite the complexity of the underlying mechanisms, the crystal exhibits relatively straightforward stress–strain curves. Most of the crystal orientations studied show rapid strain softening following the onset of yielding, which settles to a steady flow state. Trajectories are analyzed using five metrics for local states and structural order, but most of these metrics yield similar distributions for these deformation mechanisms. On the other hand, a recently proposed measure of intramolecular strain energy is found to most cleanly distinguish between these mechanisms, while also providing a plausible connection with mechanochemically accelerated decomposition kinetics. Localization of intramolecular strain energy is found to depend strongly on crystal orientation and pressure.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Summary of Geometric Parameters and Their Effects on Performance of U-10Mo Fuel Plates

A monolithic plate-type fuel system has been under development to convert high-performance test reactors from highly enriched uranium to low-enrichment uranium fuels and is now moving into the qualification phase, a predecessor to the timely conversion of the target reactors. To qualify this fuel system, the plates must meet safety standards and perform well in a reactor. The plates must maintain mechanical integrity, exhibit geometric stability, and have stable and predictable in-reactor behavior. The requirement to maintain mechanical integrity under normal operating conditions is primarily demonstrated by successful testing. However, each high-performance reactor employs a distinct design, resulting in distinct plate geometries, with unique features, attributes, irregularities, and tolerances. Due to the abundance of such distinct geometric varieties, a single “generic” plate geometry capturing all extremes is not achievable. It is also impractical to test each of these proposed designs in a reactor. This limitation necessitates cautious evaluations since the thermomechanical response of a plate with a certain geometry may not be representative for a plate with a significantly different geometry. To address concerns related to in-reactor plate performance, large set of sensitivity studies were performed. These parametric studies aimed to better understand irradiation performance, while evaluating the sensitivity of results to various modeling inputs, including geometric, operational, and material parameters. This work studied selected geometric parameters based on provided fuel specifications and performed a series of parametric simulations. The resulting temperature, displacement, and stress strains were comparatively evaluated to determine the effects of various geometric parameters. This draft provides a “high-level summary” of our parametric sensitivity studies and summarizes the key findings.

42 ENGINEERING↗
DOE OSTI2026

Modeling the contributions to acoustic nonlinearity from complex dislocation networks using 3D dislocation dynamics

Nonlinear ultrasonic parameters are highly sensitive to microstructural features that affect macroscale material behavior, providing a nondestructive means to characterize their evolution. Although dislocations are known to be a strong source of acoustic nonlinearity, establishing quantitative links between the acoustic nonlinearity parameter (β), measured via Second Harmonic Generation, and dislocation morphology—such as dislocation length and density—remains an open challenge. This work advances the numerical modeling of dislocation–β relationships using 3D dislocation dynamics (DD) simulations in two approaches: a “static” method computing strain and stress fields from dislocation configurations in the absence of external loading, and a “quasi-static” method to estimate β from the curvature of dislocation lines under applied load. First, the static method is combined with finite element analysis to investigate a recent assertion that heterogeneous initial strain fields can induce higher harmonic generation in a linear elastic medium; the present results do not corroborate this outcome. Then, the quasi-static method is applied to multiple-dislocation scenarios through parametric studies, revealing behaviors not predicted by analytical models, such as the competing interactions of edge and screw dislocations and the significant influence of applied stress on β. Finally, the simulations are used to model SHG experimental results and validate the hypothesis that β can decrease during plastic deformation, despite increasing dislocation density. As the DD code used here is open-source, it provides a practical platform for future investigation into microstructure–β relationships important to the interpretation of SHG results.

Materials science↗
DOE OSTI2021

Anomalous strain-energy-driven macroscale translation of grains during nonisothermal annealing

We report a mode of grain growth, involving the macroscopic translation of grain centers during nonisothermal annealing. Through synchrotron high-energy x-ray diffraction microscopy, we find dissolution of semicoherent precipitates generates dislocations, thereby raising the stored strain energy within grains. Here, the subsequent evolution of grains shows unexpected grain translations over length scales of 10-100 mu m. Phase-field simulations reveal that such translations are not uncommon in strain-energy-driven grain growth, wherein different regions of a grain may grow and shrink simultaneously.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Distributed electric field sensing using fibre optics in borehole environments

In the past decade, rapid advances in distributed optical fibre sensing technologies have made it possible to record various geophysical data (e.g. strain, temperature and pressure) continuously in both time and space along the fibre, providing an unprecedented quantity and spatial density of data compared to traditional geophysical measurements as well as reducing data acquisition cost. To date, no distributed fibre-based electromagnetic field sensing system has been implemented although electromagnetic sensing could have a broad range of applications to geophysical imaging and monitoring in borehole environments. The goal of this paper is to provide a theoretical feasibility study regarding the design and use of an electromagnetic sensing optical fibre for geophysical applications. First, we present the sensitivity analysis of a ‘hypothetical’ optical fibre coated with polyvinylidene fluoride, a polymer that provides relatively high piezoelectric properties, yet unlike ceramics, is flexible. Next, using a two-dimensional electromagnetic modelling algorithm, we simulate the earth electric-field-to-fibre-strain transfer function and estimate the theoretical sensitivity of the optical fibre to electric fields. Given the state-of-the-art distributed acoustic sensing strain sensitivities in the picometres strain range, our numerical modelling analysis suggests that a perfectly coupled polyvinylidene fluoride–coated optical fibre can measure electric field values in the mV/m to V/m amplitude range. We then apply a cylindrically symmetric modelling algorithm to simulate numerical models demonstrating the applicability of such a fibre in an oilfield environment. Scenarios investigated employ an electric field source and suggest that the measurements can be used to distinguish the oil versus water ratio with a fibre mounted inside a producing steel cased oil well as well as distinguishing between brine and hydrocarbon filled reservoir zones with a fibre located outside of the casing.

58 GEOSCIENCES↗
DOE OSTI2025

The influence of cooling rates on strain phase diagrams and domain structures of ferroelectric thin films: A case study of PbTiO 3

Strain engineering has been established as an effective approach to control phase equilibria, domain configurations, and functional properties of ferroelectric thin films. Temperature-strain phase diagrams have been used as powerful tools for providing insights into strain engineering. However, almost all existing phase diagrams established using the phase-field approach assume quenching conditions without considering actual cooling rates during the post-deposition annealing process of ferroelectric thin films. Within this work, we systematically investigate the influence of cooling rates on domain structures and the strain-phase diagram of ferroelectric thin films using phase-field simulations, taking PbTiO 3 thin films as a model system. We found that both the position of phase boundaries in the strain phase diagrams and the domain morphology are significantly influenced by the cooling rates. It is revealed that while the paraelectric-ferroelectric phase boundary remains invariant, the phase boundaries between single-phase and multi-phase regions tend to shift toward the corresponding multi-phase region as the cool rate reduces. Slow cooling generally leads to more ordered domain structures with increased domain size. Using the obtained equilibrium domain structures, we calculated effective thermal conductivities and found significant variations that can be tuned by the cooling rates. In conclusion, this work reveals an underexplored yet critical impact of cooling rates on phase equilibria and domain structures in ferroelectric thin films, which may inspire further fine-tuning of domains and domain walls in low-dimensional ferroelectrics for multifunctional applications.

Domain structures↗
DOE OSTI2023

BISON analysis of FeCrAl and Zircaloy cladding deformation during simulated BWR cyclic dryout conditions

The BISON fuel performance analysis code has been used to study FeCrAl and Zircaloy cladding performance to better understand cladding deformation behavior in boiling water reactors during conditions caused by anticipated operational occurrences (AOOs) or by anticipated transients without scram (ATWS). BISON models of novel separate effects tests simulating ATWS events predicted mechanical responses to thermal cycling boundary conditions. Here, we assessed FeCrAl and Zircaloy deformation and burst models to determine those best suited for cyclic dryout event analysis. In addition, we compared diametral strain profiles of the post-experiment claddings to simulation results. The capability of the BISON code for modeling BWR AOO or ATWS events is promising.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Simulation Analysis of Geometry and Material Effects for Dropkinson Bar

The reported research is motivated by the need to address a key issue affecting the Dropkinson bar apparatus. This unresolved issue is the interference of the stress wave reflected from the bar-beam boundary with the measurement of the stress-strain response of a material tested in the apparatus. The purpose of the wave beam that is currently connected to the bar is to dissipate the stress wave, but the portion of the wave reflected from the bar-beam boundary is still significant. First, we focused on understanding which parameters affect the reflected wave's arrival time at a strain gauge. Specifically, we used finite-element numerical simulations with the Sierra/SM module to study the effects of various bar-beam connection fixities, alternative wave beam materials, and alternative geometries of the Dropkinson bar system based on a monolithic design. The conclusion of this study is that a partial reflection always occurs at the bar-beam boundary (or, for a monolithic design, at a point where the bar geometry changes). Therefore, given a fixed total length of the bar, it is impossible to increase the reflected wave's arrival time by any significant amount. After reaching this conclusion, we focused instead on trying to minimize the energy of the reflected stress wave circulating up and down through the bar over a relatively long period of time (10 ms). Once again, we used numerical simulations with the Sierra/SM module to investigate the effects of various bar-beam connection fixities, alternative wave beam materials, and parameters of an asymmetric monolithic design of the bar-and-beam system. This study demonstrated that various parameters can significantly affect the energy of the wave reflections, with the difference between best and worst configurations being about one order of magnitude in terms of energy. Based on the obtained results, we conclude with concrete takeaways for Dropkinson bar users and propose potential directions for future research and optimization.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Unveiling damage mechanisms of chromium-coated zirconium-based fuel claddings at LWR operating temperature by in-situ digital image correlation

Here, we investigate the coupled thermomechanical fracture mechanisms of coated nuclear fuel claddings at Light Water Reactor (LWR) operating temperatures. These coated claddings are a highly attractive, near-term solution, which addresses the demands for accident-tolerant fuel systems and provide greater oxidation resistance. However, the fracture mechanisms of these coatings, which may create channels for oxidation ingression, must be fully understood prior to implementation. Thus, high-temperature expanding plug experiments were conducted on coated cladding specimens at a temperature of 315 °C, consistent with the operating environment of LWRs. Furthermore, in-situ thermomechanical deformation was measured with stereo digital image correlation during heating and mechanical testing to separately resolve contributions of thermal and mechanical strain. Digital image correlation, supported by acoustic emissions (AE) detection, was also leveraged to track cracking activity during loading. Coating fracture was found to initiate at total hoop strains of 0.34%. The thermal deformation of the coated claddings was investigated via finite element simulations, revealing a bidirectional stress-state within the coating with axial and circumferential strains of 0.026 and 0.031%. This bidirectional stress-state was attributed with the generation of off-axis fracture pattern within the coating as identified via post-experiment scanning electron microscopy. Thus, this study unveiled critical, coupled thermomechanical mechanisms governing the coating fracture of coated claddings at LWR temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Additively manufactured strain sensors for in-pile applications

Accurate, real-time monitoring of strain in fuel, cladding, and structural components of nuclear reactors is critical to better understand radiation induced phenomena during reactor tests and operations. The data provided is crucial to verify physics-based, multiscale modeling and simulation efforts, which aim to shorten the timeline for the development of new nuclear materials. Resistive strain gauges have limited performance during in-pile experiments due to the harsh operating conditions and limited physical space between fuel and cladding components. In this work, aerosol jet printing using silver nanoparticle inks was used to fabricate interdigitated electrode capacitive strain gauges on aluminum alloy 6061 tensile specimens. Here to simulate the temperatures of a traditional light water reactor, the capacitive strain gauges were tested with a mechanical test frame up to 300 °C and compared to commercially available bondable resistive strain gauges. The printed capacitive strain gauges exhibited a gauge factor of 1.0 and showed higher reproducibility and predictability of strain sensing performance than the resistive strain gauge. The results demonstrate the potential of aerosol jet printing to fabricate strain sensors with predictable performance and reduced invasiveness for high-temperature applications with confined spacing.

36 MATERIALS SCIENCE↗
DOE OSTI2020

The effects of heat treatment on the response of WE43 Mg alloy: crystal plasticity finite element simulation and SEM-DIC experiment

The present report addressed the effects of heat treatment on the mechanical response of a WE43 Mg alloy using an integrated framework of SEM-DIC experiment and CPFE simulation. Both macroscopic responses including yield strength, ultimate strength, ductility, and microscopic responses, including local displacement and strain maps, were experimentally investigated. The focus of this work is to use the CPFE simulation as an integrated computational tool to study the effects of heat treatment. The CPFE framework was evaluated using the local fields of displacement and strain obtained from the SEM-DIC experiment rather than the conventional scheme of using macroscopic responses at different loading directions. Subsequently, the information which is available using CPFE, such as the critical resolved shear stress (CRSS) and relative slip activity, was used to study the effects of heat treatment on the response of WE43 Mg alloy. The contributions of different strengthening mechanisms on the CRSS were addressed. The results show that effects of heat treatment can be captured using the predominant mechanisms of the grain size effect and the influence of precipitates. Finally, it has been shown that classical Hall-Petch in which one constant can capture the size effects, should be modified. To do so, each deformation mode should have a unique Hall-Petch constant, which are calculated here for the WE43 Mg alloy.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Atomistic insights into the synergistic effects of tensile strain and hydroxyl group on increasing the thermal conductivity of monohydric alcohols as latent heat storage materials

The low thermal conductivity of linear monohydric alcohols significantly limits their performance when used as latent heat storage materials. It was found that imposing certain mechanical strains on nanomaterials with highly ordered structures can lead to appreciable thermal conductivity enhancement. To explore the potential of enhancing the intrinsic thermal conductivity of monohydric alcohols, a simplified ideal crystal model was employed for molecular arrangement in this work. The variation of atomistic heat transfer as a function of tensile strains was exploited by non-equilibrium molecular dynamics simulations. Furthermore, the influence of stretching on the vibration density of states and the corresponding molecular morphology were revealed. The results suggested that the thermal conductivity of monohydric alcohols with such simplified ideal crystalline structures can increase by about 170 % at the strain of 0.1. The distance between neighboring interfaces was found to determine the energy of inter-molecular interactions at the crystal grain boundaries. The hydrogen bonds energy has not reached the critical value when increasing tensile strains from 0 to 0.1. In contrast, the strongest van der Walls interactions for alkanes were obtained at the strain of 0.06 where the interfacial distances of alkanes have achieved their equilibrium distances. In addition, the interfacial heat transfer coefficient of ideal crystal monohydric alcohols at hydroxyl interfaces was found to be threefold than that between methyl interfaces, in consistent with the strength of inter-molecular interactions between functional groups at the interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Dynamic Strength of Copper at High Pressures Using Pressure Shear Plate Experiments

Measurement of strength at high pressures is critical to understanding the behavior of materials subjected to extreme loading conditions. Recent developments have extended the pressure shear plate impact (PSPI) technique to the high-pressure regime. Such experiments provide a unique opportunity to extract the complete stress-strain behavior of materials at relatively high pressures. The modified technique includes a new fiber-optic heterodyne transverse velocity interferometer system and new analysis methods to account for the inelastic response of the anvil plates. In this study, PSPI experiments are conducted on oxygen-free high conductivity (OFHC) copper at pressures ranging from 10 to 43 GPa and at strain rates of ~10 5 s -1 . Complete stress-strain curves of copper are obtained using a hybrid methodology that relies on numerical simulations and particle velocity records gathered in these experiments. Experiments are conducted at similar pressures and strain rates using two types of anvils that undergo different levels of plastic strains: tool steel and tungsten carbide, to check the robustness of the hybrid method. The behavior of copper at these high pressures and strain rates is compared with previous literature data, and the effect of pressure and strain rate on the strength behavior of copper are discussed. Here, it was observed that the rate of change in yield strength change in pressure was 2.5 times the rate of change of shear modulus with pressure. The strength models that predict the response of copper at these extreme conditions are also examined and validated.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Modeling Occupant Core Temperatures Across the Boston Building Stock to Advance Public Health

In recent history, extreme heat has been the cause of most deaths from a natural disaster. Exposure to extreme heat can aggravate preexisting conditions, increase hospitalization, and even cause death. Thus, modeling the thermal resilience of households across the United States will allow for a quantitative assessment of the health and safety risks posed by extreme heat. For the first time, we simulate occupant comfort in representative households across Boston by combining the granular results of the ResStock(TM) model with a two-node heat strain model. This model calculates occupants' core body temperature in each simulated household over a year. We compare the simulated core temperatures to two public health metrics: hyperthermia and heat stroke. Our results show for households in Boston that do not have or use aid conditioning, thousands potentially experience many dangerous heat events each summer and these events can last for more than a day at a time. These heat events peak during the late afternoon and evening, just as residents are coming home, cooking meals, and trying to go to sleep. We have found that multi-family buildings, renters, and low-income homes in more airtight and insulated homes are the most at risk for these events. These results demonstrate the scale and urgency of exposure to extreme heat.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Effect of grain structure and strain rate on dynamic recrystallization and deformation behavior: A phase field-crystal plasticity model

Herein, a phase field (PF)-crystal plasticity (CP) model has been developed to study the effects of initial grain structures, recrystallized grain orientation, and strain rates on dynamic recrystallization and deformation behavior under stress states typical of the Shear Assisted Processing and Extrusion (ShAPE) process. The model takes the geometrically necessary dislocation density as the nucleation criterion for recrystallization and minimizes the deformation energy as the driving force of recrystallized grain growth. Using this model, the effects of grain orientations and strain rates on dynamic recrystallization and deformation behavior were simulated. The results demonstrate that (1) a polycrystalline structure with the texture observed in ShAPE process has the lowest yield stress under the ShAPE stress conditions, (2) the recrystallized grains with the texture observed in ShAPE process largely soften the materials, and (3) higher shear strain rate or rotation speed results in larger magnitude of material softening.

36 MATERIALS SCIENCE↗