Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model deformation”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A Fusion of Geothermal and InSAR Data with Machine Learning for Enhanced Deformation Forecasting at the Geysers

The Geysers geothermal field in California is experiencing land subsidence due to the seismic and geothermal activities taking place. This poses a risk not only to the underlying infrastructure but also to the groundwater level which would reduce the water availability for the local community. Because of this, it is crucial to monitor and assess the surface deformation occurring and adjust geothermal operations accordingly. In this study, we examine the correlation between the geothermal injection and production rates as well as the seismic activity in the area, and we show the high correlation between the injection rate and the number of earthquakes. This motivates the use of this data in a machine learning model that would predict future deformation maps. First, we build a model that uses interferometric synthetic aperture radar (InSAR) images that have been processed and turned into a deformation time series using LiCSBAS, an open-source InSAR time series package, and evaluate the performance against a linear baseline model. The model includes both convolutional neural network (CNN) layers as well as long short-term memory (LSTM) layers and is able to improve upon the baseline model based on a mean squared error metric. Then, after getting preprocessed, we incorporate the geothermal data by adding them as additional inputs to the model. This new model was able to outperform both the baseline and the previous version of the model that uses only InSAR data, motivating the use of machine learning models as well as geothermal data in assessing and predicting future deformation at The Geysers as part of hazard mitigation models which would then be used as fundamental tools for informed decision making when it comes to adjusting geothermal operations.

Yazbeck, Joe (ORCID:0000000302235260)↗
DOE OSTI2023

Thermomechanical Modeling of Copper

The plastic deformation of metals is a dissipative process. Some fraction of the plastic work is converted to heat which, given the temperature dependent response of metals, produces a thermal-mechanical coupling. In various cases, for instance when the loading is dynamic, this interaction can impact the resulting response of a material and/or system. Thus, appropriately capturing the heat generation from plastic work is necessary for various solid mechanics analysis. Determination of the fraction of work converted to heat has been long studied. Recent developments have demonstrated that the fraction is not constant but depends on various state variables. Resolving these features requires combined modeling and experimental studies. To this end, 304L stainless steel – a poor thermal conductor – was recently subjected to such an investigation. Advanced modeling capabilities were deployed to assess novel thermomechanically coupled experiments. As a complement to that study, in the current work a similar investigation is performed on copper – a good thermal conductor – to assess performance on the opposite end of the spectrum. The current document discusses these modeling efforts.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Three-dimensional modeling and simulations of single-crystal and bi-crystal titanium for high-strain-rate loading conditions

High purity single crystal titanium (Ti) under shock wave loading is modeled under both one-dimensional and three-dimensional cylindrical conditions. Cylinder sizes of 10 μm and 20 μm radius are both considered in order to assess influence of boundary conditions. A thermodynamically consistent single crystal model for application to shock conditions is presented. The model accounts for the coupled non-linear elastic, dislocation slip, deformation twinning, and structural phase transformation response of the titanium material. Plate impact experiment results using a copper flyer are used to compare against the simulations for crystals oriented in [0001] and [10$\overline1$1] crystallographic directions. The one-dimensional and threedimensional simulations of the two differently oriented single crystals indicate differences between the one-dimensional and three-dimensional representation, especially for the [10$\overline1$1] oriented single crystal. This orientation breaks the relative orientation symmetry between the crystal and cylinder which otherwise exists for the [0001] oriented single crystal. A significant amount of heterogeneity in the field response of the[10$\overline1$1] oriented simulation was demonstrated due to the highly coupled nature of the deformation. A bi-crystal model composed of both the [0001] and [10$\overline1$1] orientations with the boundary between the two along the axis of the cylinder is also considered for a cylinder model size of 10 μm. The results indicate a strong interaction between the two grains that affects the ω phase volume fraction achieved relative to the single crystal calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Density matrix renormalization group description of the island of inversion isotopes F 28 – 33

Recent experiments have confirmed that the neutron-rich isotopes 28,29 F belong to the so-called island of inversion (IOI), a region of the nuclear chart around Z = 10 and N = 20 where nuclear structure deviates from the standard shell model predictions due to deformation and continuum effects. However, while the general principles leading to the IOI are relatively well understood, the details of the low-lying structure of the exotic fluorine isotopes 28–33 F are basically unknown. In this study, we perform large-scale shell model calculations including continuum states to investigate the properties of the neutron-rich isotopes 25–33 F, from a core of 24 O and using an effective two-body interaction with a small number of adjustable parameters in the central and tensor channels. We develop two models adjusted on experimentally confirmed states in 25,26 O and 25–27 F based on different assumptions concerning the positions of the neutron 0d 3/2 and 1p 3/2 shells, and solve the many-body problem using the density matrix renormalization group (DMRG) method for open quantum systems in an sd–fp model space. We obtain the first detailed spectroscopy of 25–33F in the continuum and show how the interplay between continuum effects and deformation explains the recent data on 28,29 F. Several deformed one- and two-neutron halo states are predicted in 29,31 F, and we provide some information about the possible structure of the heaviest fluorine isotopes. We also suggest several experimental studies of interest to constraint models and test the present predictions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Landscape of pear-shaped even-even nuclei

The phenomenon of reflection-asymmetric nuclear shapes is relevant to nuclear stability, nuclear spectroscopy, nuclear decays and fission, and the search for new physics beyond the standard model. Global surveys of ground-state octupole deformation, performed with a limited number of models, suggest that the number of pear-shaped isotopes is fairly limited across the nuclear landscape. Here, we carry out a global analysis of ground-state octupole deformations for particle-bound even-even nuclei with Z ≤ 110 and N ≤ 210 using nuclear density functional theory (DFT) with several nonrelativistic and covariant energy density functionals. In this way, we can identify the best candidates for reflection-asymmetric shapes. The calculations are performed in the frameworks of axial reflection-asymmetric Hartree-Fock-Bogoliubov theory and relativistic Hartree-Bogoliubov theory using DFT solvers employing harmonic oscillator basis expansion. We consider five Skyrme and four covariant energy density functionals. We predict several regions of ground-state octupole deformation. In addition to the “traditional” regions of neutron-deficient actinide nuclei around 224 Ra and neutron-rich lanthanides around 146 Ba, we identified vast regions of reflection-asymmetric shapes in very neutron-rich nuclei around 200 Gd and 288 Pu, as well as in several nuclei around 112 Ba. Our analysis suggests several promising candidates with stable ground-state octupole deformation, primarily in the neutron-deficient actinide region, that can be reached experimentally. Detailed comparison between Skyrme and covariant models is performed. Octupole shapes predicted in this study are consistent with the current experimental information. This work can serve as the starting point of a systematic search for parity doublets in odd-mass and odd-odd nuclei, which will be of interest in the context of new physics searches.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Phase-field modeling of rate-dependent fluid-driven fracture initiation and propagation

The rate-dependent behavior associated with deformation and fracturing of materials, such as natural rocks, poses significant challenges for modeling. In addition to the complications of the viscoelastic response, the speed of fracture propagation reflects micromechanical mechanisms in the fracture process zone(FPZ). In order to represent these complicated behaviors, a thermodynamically consistent, rate-dependent fracture model is required. Based on rigorous thermodynamic principles, we derive a rate-dependent phase-field mechanical model coupled with single-phase fluid flow in both the matrix and the fracture. The model is guaranteed to satisfy energy conservation during fracture propagation. Here, the system of equations is solved using the introduced solution procedure and a novel preconditioner that accounts for the complex fluid-structure interaction. The proposed phase-field model is tested against several benchmark problems on solid-fluid coupling, fluid-driven fracture propagation and rate-dependent viscoelastic deformation. The model serves as a strong basis for investigating rate-dependent fracturing experiments and for making predictions of material behaviors under new conditions.

42 ENGINEERING↗
DOE OSTI2021

Seismic anisotropy, dominant slip systems and phase transitions in the lowermost mantle

SUMMARY The presence of seismic anisotropy at the base of the Earth's mantle is well established, but there is no consensus on the deformation mechanisms in lower mantle minerals that could explain it. Strong anisotropy in magnesium post-perovskite (pPv) has been invoked, but different studies disagree on the dominant slip systems at play. Here, we aim to further constrain this by implementing the most recent results from atomistic models and high-pressure deformation experiments, coupled with a realistic composition and a 3-D geodynamic model, to compare the resulting deformation-induced anisotropy with seismic observations of the lowermost mantle. We account for forward and reverse phase transitions from bridgmanite (Pv) to pPv. We find that pPv with either dominant (001) or (010) slip can both explain the seismically observed anisotropy in colder regions where downwellings turn to horizontal flow, but only a model with dominant (001) slip matches seismic observations at the root of hotter large-scale upwellings. Allowing for partial melt does not change these conclusions, while it significantly increases the strength of anisotropy and reduces shear and compressional velocities at the base of upwellings.

58 GEOSCIENCES↗
DOE OSTI2024

Life at the Landau Pole

If a quantum field theory has a Landau pole, the theory is usually called ‘sick’ and dismissed as a candidate for an interacting UV-complete theory. In a recent study on the interacting 4d O(N) model at large N, it was shown that at the Landau pole, observables remain well-defined and finite. In this work, I investigate both relevant and irrelevant deformations of the said model at the Landau pole, finding that physical observables remain unaffected. Apparently, the Landau pole in this theory is benign. As a phenomenological application, I compare the O(N) model to QCD by identifying $Λ_{\overline{MS}}$ with the Landau pole in the O(N) model.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Micro-structuration effects on local magneto-transport in [Co/Pd]IrMn thin films

We measured the local magneto-transport (MT) signal with an out-of-plane magnetic field, including magneto-resistance (MR) and Extraordinary Hall effect (EHE), in exchange-biased [Co/Pd]IrMn thin multilayers that are micro-structured with a 100 μm window. We found that when measured locally around the window, the MT signal deviate from the expected behavior. We studied possible causes, including film micro-structuration, electrical contact geometry as well as magnetic field angular tilt. We found that tilting the magnetic field direction with respect to the normal direction does not significantly affect the MT signal, whereas the positioning and geometry of the contacts seem to highly affect the MT signal. For comparison purposes, we carried these MT measurements using the Van-der-Pauw method on a set of four microscopic contacts directly surrounding the window, and on another set of micro-contacts located outside the window, as well as a set of four contacts positioned several millimeters away of each other at the corners of the wafer. If the contacts are sufficiently far apart, the EHE and MR signals have the expected shape and are not significantly affected by the presence of the window. If, on the other hand, the contacts are micro-positioned, the shape of the EHE signal is drastically deformed, and may be modeled as a mix of the standard EHE and MR signals measured on the outer contacts. Furthermore, if the micro-contacts are located directly around the window, the deformation is amplified, and the weight of the MR signal in the mix is further increased by about 40 %. This suggests that the electron path in the Hall geometry is disturbed by both the proximity of the electrodes and by the presence of the window, which both contribute to the deformation for about two-third and one third, respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A multiscale cohesive law for carbon fiber networks

Better predictive models of mechanical failure in low-weight heat shield composites would aid material certification for missions with aggressive atmospheric entry conditions. In this study, we develop such a model for the rapid engineering analysis of the failure limits of phenolic impregnated carbon ablator (PICA) - a leading heat shield material whose structural component is a carbon fiber network. We hypothesize inelastic deformation failure mechanisms and model their behavior using molecular dynamics simulations to calculate the binding energy. We then upscale this binding energy to the macroscale using a renormalization argument. The approach delivers insightful and reasonably accurate macroscale predictions that compare favorably to experiments. In application, the model is validated for a particular variety of PICA by comparison to experiment and would then be used to study design scenarios in different entry conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Towards a multiscale approach for understanding irradiation induced swelling and creep in 316 stainless steels - A coupled cluster dynamics and crystal plasticity approach

Structural materials undergo mechanical degradation, in part due to irradiation-induced swelling and creep, under nuclear reactor conditions. While swelling results from the migration and clustering of irradiation-induced atomic-scale mobile defects, the interaction of mesoscale dislocations with these defects causes creep deformation. A coupled crystal plasticity (CP) and mean-field cluster dynamics (CD) approach is presented to investigate the effect of irradiation on the long-term mechanical behavior of 316 stainless steel, which are under consideration for use in nuclear reactors. The temporal evolution of Frenkel pairs and extended defect population, under a chosen irradiation flux and temperature, is predicted using the CD model. The impact of the irradiation defects on the stress state, and the resulting dislocation-mediated inelastic deformation, is modeled concurrently with the CP model. The inelastic deformation is irradiation flux dependent, and early-stage defect evolution determines the later-stage mechanical behavior in 316 stainless steel.

36 - MATERIALS SCIENCE↗
DOE OSTI2022

Maximum strength and dislocation patterning in multi–principal element alloys

Multi–principal element alloys (MPEAs) containing three or more components in high concentrations render a tunable chemical short-range order (SRO). Leveraging large-scale atomistic simulations, we probe the limit of Hall-Petch strengthening and deformation mechanisms in a model CrCoNi alloy and unravel chemical ordering effects. The presence of SRO appreciably increases the maximum strength and lowers the propensity for faulting and structure transformation, accompanied by intensification of planar slip and strain localization. Deformation grains exhibit notably different microstructures and dislocation patterns that prominently depend on their crystallographic orientation and the number of active slip planes. Grain of single-planar slip attains the highest volume fraction of deformation-induced structure transformation, and grain with double-slip planes develops the densest dislocation network. These results advancing the fundamental understanding of deformation mechanisms and dislocation patterning in MPEAs suggest a mechanistic strategy for tuning mechanical behavior through simultaneously tailoring grain texture and local chemical order.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Stress evolution and creep deformation in solid-oxide electrolysis cell systems – Dynamic modeling and multi-objective optimization to maximize stack life and efficiency

Here, this study develops a thermal stress model of solid-oxide electrolysis cells (SOECs) including a model for creep strain and failure probability that is integrated with a dynamic plant-wide model of a hydrogen production process. Uncertainties in key material properties of the cell are quantified to assess their impact on stress profile variability. The oxygen electrode is found to have about 10 times higher failure probability compared to the fuel electrode. The study shows that if the stack operation is not optimized, cycling operation would lead to stress build-up eventually leading to catastrophic failure. A dynamic optimization problem is set up for obtaining the optimal operational profile considering a variable hydrogen production rate. Due to the tradeoff between the efficiency and stress build-up, the dynamic optimization problem is multi-objective. It is observed that the optimizer can considerably reduce the stress build-up (i.e., can increase the stack life) albeit at the cost of a lower efficiency thus exhibiting strong tradeoffs between capital and operating costs. For example, if the stack would be replaced in 0.5 yr, specific energy requirement would be 48.5 kWh/kg H 2 while for a stack replacement time of about 6 yr, the specific energy requirement rises by about 4.2 %.

SOEC↗
DOE OSTI2022

Analysis of an Elasto-Hydrodynamic Seal by Using the Reynolds Equation

This paper reports numerical studies of an Elasto-Hydrodynamic (EHD) seal, which is being developed for supercritical CO2 (sCO2) turbomachinery applications. Current sCO2 turbomachinery suffers from high leakage rates, which is creating a major roadblock to the full realization of sCO2 power technology. The high leakage rates not only penalize the efficiencies but also create environmental concerns due to greenhouse effects caused by the increased CO2 discharge to the atmosphere. The proposed EHD seal needs to work at elevated pressures (10–35 MPa) and temperatures (350–700 °C) with low leakage and minimal wear. The unique mechanism of the EHD seal provides a self-regulated constriction effect to restrict the flow without substantial material contact, thereby minimizing leakage and wear. This work utilizes a physics-based modeling approach. The flow through the gradually narrowing seal clearance is modeled by the well-known Reynolds equation in EHD lubrication theory, while the deformation of the seal is modeled by using the governing equations of three-dimensional solid mechanics. As for the solution methodology, COMSOL’s Thin-Film Flow and Solid Mechanics modules were employed with their powerful capabilities. The numerical results were presented and discussed. It was observed that the Reynolds equation fully coupled with the surface deformation was able to successfully capture the constriction effect. The maximum and minimum leakages were calculated to be 2.25 g/s and 0.1 g/s at P = 5.5 MPa and P = 11 MPa for the design seal, respectively. It was interesting to observe that the seal leakage followed a quadratic trend with increasing pressure differential, which can become advantageous for high-pressure applications such as sCO2 power generation technology.

42 ENGINEERING↗
DOE OSTI2024

A macro-micro approach for identifying crystal plasticity parameters for necking and failure in nickel-based alloy haynes 282

Here, this work develops a two-scales macro-micro approach to address the challenge in calibrating crystal plasticity microstructural models when samples undergo necking prior to fracture. The crystal plasticity models are crucial for predicting the materials’ plastic deformation and failure at the microstructure level, identifying the materials’ intrinsic properties as well as investigating the microstructure-properties relationships. However, after necking occurs, the experimentally measured stress-strain curves fail to reflect the materials ‘true’ stress-strain behavior and cannot be directly fitted into crystal plasticity models. The proposed macro-micro approach employs a top-down strategy to address this challenge, which has been studied with experimental tests on precipitation-strengthened Ni-based superalloy Haynes® 282®. In this approach, a macro rate-dependent anisotropic plasticity model with Voce-type hardening and Rice-Tracey damage law is first utilized to model the deformation and failure of the tensile bar, and calibrated by matching the stress-strain curves, necking strain, and reduction of area. Especially, to match the testing results under different applied strain rates, the rate-sensitivity parameter m and saturation stress in the elasticity model are modified to incorporate dependence on the local strain rate. Then, the ‘true’ stress-strain behaviors are extracted from the necking zone of the macro-model, which are used to calibrate a micro-model with explicit microstructures and governed by an extended crystal plasticity law. The consistency between the micro-model and macro-model are enforced during calibration. The calibration outcomes from the crystal plasticity model elucidate the materials intrinsic properties for slip, hardening, and failure, which is vital for further investigations into the microstructure-properties relationship and for accurate prediction of the material behavior under various test and service conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Localization of deformation in fluid-saturated sandstone

Herein we discuss a laboratory study of inelastic deformation and localization was conducted to evaluate microseismicity and fracture in fluid-saturated rock. Plane strain compression experiments were performed on oil- and water-saturated Berea sandstone under dry, drained, undrained, and quasi-unjacketed conditions. Parameters associated with the elastic and inelastic material response were accurately measured and used in a constitutive model to predict localization of deformation under the drained and undrained conditions. Acoustic emission (AE) activity was recorded and the onset of inelastic response in fluid-saturated rock, similar to dry rock, coincided with an increase in AE rate. Fracture initiation and propagation were successfully monitored by AE locations. Clustering of AE events and the start of inhomogeneous deformation inferred from displacement measurements are shown to be well correlated with the model prediction of the onset of localized deformation under various rock-fluid interaction conditions.

42 ENGINEERING↗
DOE OSTI2023

A mechanistic model for creep and thermal aging in Alloy 709

This report describes a physics-based model for creep and thermal aging in Alloy 709. Alloy 709 is an advanced austenitic alloy, targeted for use in future Sodium Fast Reactors (SFRs) and other advanced reactors. The material has superior high temperature properties compared to currently qualified 316 and 304 stainless steels. However, the available creep and thermal aging test database for Alloy 709 is significantly more limited compared to the historical materials. The physics-based model developed here is one way to accelerate the qualification of the material by providing more accurate long-term predictions for creep properties and thermal aging, compared to current empirical time-extrapolate techniques. The crystal plasticity finite element model is used to predict the deformation and failure of alloy 709. The same setup for the CPFE model is used in both the baseline model calibration process and the simulation campaigns for parameter inference. Specific constitutive choices are made for Alloy 709 to capture the primary deformation mechanisms. The dislocation creep formulation developed by Hu and Cocks is extended to account for coupled precipitation formation and the grain boundary cavitation model developed by Sham, Needleman, et al. is used to model grain boundary cavitation-induced failure. A novel update algorithm is proposed to render the semi-discrete constitutive update for the Sham-Needleman model unconditionally stable. A progressive calibration approach is adopted based on the observations that several types of material responses can be effectively decoupled. A surrogate model is trained based on full-fledged CPFE simulations to accelerate the forward model evaluations, and stochastic variational inference (SVI) is used to calibrate the unknown microstructural model parameters. The calibrated mechanistic model is used to predict the long-term creep life of Alloy 709, and the predictions are compared against classical empirical approaches.

36 MATERIALS SCIENCE↗