Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “full-field”

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 109 records · Page 6

Combined thermographic phosphor and digital image correlation (TP + DIC) for simultaneous temperature and strain measurements

We report thermographic phosphors (TP) are combined with stereo digital image correlation (DIC) in a novel diagnostic, TP + DIC, to measure full-field surface strains and temperatures simultaneously. The TP + DIC method is presented, including corrections for nonlinear CMOS camera detectors and generation of pixel-wise calibration curves to relate the known temperature to the ratio of pixel intensities between two distinct wavelength bands. Additionally, DIC is employed not only for strain measurements but also for accurate image registration between the two cameras for the two-colour ratio method approach of phosphoric thermography. TP + DIC is applied to characterize the thermo-mechanical response of 304L stainless steel dog bones during tensile testing at different strain rates. The dog bones are patterned for DIC with Mg 3 F 2 GeO 4 :Mn (MFG) via aerosol deposition through a shadow mask. Temperatures up to 425°K (150°C) and strains up to 1.0 mm/mm are measured in the localized necking region, with conservative noise levels of 10°K and 0.01 mm/mm or less. Finally, TP + DIC is compared to the more established method of combining infrared (IR) thermography with DIC (IR + DIC), with results agreeing favourably. Three topics of continued research are identified, including cracking of the aerosol-deposited phosphor DIC features, incomplete illumination for pixels on the border of the phosphor features, and phosphor emission evolution as a function of applied substrate strain. This work demonstrates the combination of phosphor thermography and DIC and lays the foundation for further development of TP + DIC for testing in combined thermo-mechancial environments.

304L stainless steel↗

Physics-informed neural networks for identification of material properties using standing waves

A metallic structure in its initial stage of failure involves plastic deformation or environmental degradation that changes the elastic modulus and density. This work presents the detection of change in wave velocity (a function of elastic modulus and density) as a system identification problem. A physics-informed neural network (PINN) is proposed to solve the system identification problem. The PINN takes the spatial coordinates of scanning locations and time as inputs and provides the displacement and wave velocity as outputs. The governing partial differential equation of standing waves in a rod is incorporated into the neural network as physics in the form of a loss function. The wave velocity vector is randomly initiated. During the training of the network, physics is used to determine and update the wave velocity target vector from the network’s displacement predictions. The measured data, comprising sparse displacement response on the rod structure, are used to train the PINN. The wave velocity at the sparse locations on the rod is learned from the predicted displacements during the training. Using the predictions of the trained network, the response of free vibration or material property variation can be reconstructed at unscanned locations on the structure to obtain high-resolution maps for full-field imaging to detect and localize the changes caused by plastic deformation. The PINN’s sparse scanning and simultaneous prediction capability during training can lead to high scanning and data-processing speeds. This capability yields a nondestructive evaluation system that can predict the presence of degraded material locations as the structural vibrations are scanned and processed in real time.

Rathod, Vivek↗

Discovery of spontaneous mesoscopic strain waves in nematic domains using dark-field x-ray microscopy

Electronic nematic order arises when correlated electrons spontaneously break the rotational symmetry of a crystal lattice. When electronic nematic order couples bilinearly to symmetry-breaking lattice strain, both appear together at a single ferroelastic phase transition, producing structural twin domains with distinct orientations of the nematic director. While the effects of applied strain on these domains are well established, the intrinsic behavior of spontaneous subdomain strain fields has remained unexplored. Here, we report the discovery of spontaneous mesoscopic strain waves within individual nematic domains of an iron-based superconductor, observed using dark-field x-ray microscopy (DFXM). Using this advanced full-field imaging technique, we visualize subdomain strain modulations emerging concurrently with nematic order. Elastic compatibility relations governing inhomogeneous strains provide a natural mechanism for the strain waves. Our findings reveal a broadly relevant form of strain self-organization and position DFXM as a powerful probe of the local interplay between lattice strain and electronic order.

36 MATERIALS SCIENCE↗

DiscoFluxM version 1.x

DiscoFluxM is software for modeling the deformation behavior of crystalline materials within the third-party finite-element based software, MOOSE. The software enables full-field crystal plasticity finite element simulations of deformation processes of metals and other crystalline materials. DiscoFluxM is an implementation of the DiscoFlux theoretical and algorithmic framework for coupling single-crystal plasticity with dislocation transport. The continuum mesoscale modeling framework is founded upon a nonlocal dislocation-density based crystal plasticity theory. The nonlocal theory couples continuum dislocation transport with an otherwise spatially-local single-crystal model employing nonlinear thermoelasticity and crystallographic plasticity. Dislocation transport is modeled by enforcing dislocation conservation at a slip-system level through the solution of advection-diffusion equations. The configuration of geometrically necessary dislocation density gives rise to a back-stress that inhibits or accentuates the flow of dislocations.

Luscher, Darby↗

High fidelity blade-resolved and actuator line data from a 16 turbine wind farm simulation using ExaWind

This data was generated with the ExaWind code suite (https://github.com/Exawind) as a demonstration of a large, 16 turbine wind farm simulation, calculated using two different levels of fidelity. The lower level of fidelity approach uses an actuator line approach to represent the turbines, and was simulated with AMR-Wind (https://github.com/Exawind/amr-wind/) as the background flow solver, coupled to OpenFAST (https://github.com/OpenFAST/openfast). The higher level of fidelity simulation uses a blade-resolved approach, and is done using AMR-Wind, Nalu-Wind (https://github.com/Exawind/nalu-wind), OpenFAST, and TIOGA (https://github.com/Exawind/tioga). In the blade-resolved simulation, ExaWind couples together a background flow solver, AMR-Wind, and a near-body solver, Nalu-Wind, through an overset technique from the TIOGA application. OpenFAST handles the structural dynamics of the turbine blades and towers, which informs the fluid-structure interaction of the wind turbines with the flow solvers. In the actuator line simulation, a mesh of 295M elements was used for a 5km x 5km domain, and it was simulated using 256 nodes (2048 GPU's) on the Oak Ridge Leadership Computing Facility Frontier supercomputer. For the blade-resolved simulation, 1.5B element mesh was used in the AMR-Wind background 5km x 5km domain, and 16M elements were used for each turbine in the Nalu-Wind domains, for a total of 1.7B elements. This was simulated using 384 nodes on Frontier, with each node using 56 cores for Nalu-Wind and 8 GPU cores. The data in this archive includes the turbine outputs from OpenFAST, 2D sampling planes from AMR-Wind, and full-field solution files from AMR-Wind and Nalu-Wind.

17 WIND ENERGY↗

Correcting for digital image correlation speckle inversion at high temperature using color cameras

Digital image correlation (DIC) is a popular, noncontacting technique to measure full-field deformation by using cameras to track the motion of an applied surface pattern. Because it is noncontacting, DIC can be performed for extreme temperature applications (e.g., hot-fire rocket testing of carbon composite rocket nozzles) under harsh conditions during which bonded gauges are damaged. Speckle pattern inversion is a phenomenon that sometimes occurs while performing high-temperature DIC. During speckle pattern inversion, portions of the surface pattern that were initially darker at room temperature (e.g., graphite) may emit more light due to blackbody radiation than the portions that were initially paler, thereby producing images in which the pattern appears inverted at high temperature relative to the initial pattern at room temperature. This phenomenon can prevent the correlation algorithm from being able to resolve the displacements between images. This work compares three methods to mitigate speckle pattern inversion: (A) the subtraction method, a recently-published technique in which two high-temperature images are subtracted to remove unwanted light; (B) the filtering method, a popular technique in which optical bandpass filters screen out unwanted light; and (C) the histogram rescaling method, a proposed new method that pairs a color camera with a blue light source and uses information from the green sensor of the camera to correct against inversion in the blue sensor through postprocessing. The histogram rescaling method is shown to successfully eliminate speckle pattern inversion and has the added advantages that it does not require quasi-static loading to be able to compensate for speckle pattern inversion, nor does it impose thick-glass distortions caused by the optical filter.

Rowley, Lindsey J.↗

Refractive axicon for X-ray microscopy applications: design, optimization, and experiment

In a full-field transmission X-ray microscopy (TXM) setup, a condenser X-ray optical element is used to illuminate the sample by condensing the X-ray beam delivered by the synchrotron storage ring. On-going and future upgrades of synchrotron facilities to diffraction-limited storage rings will pose new challenges to these TXM setups, such as much smaller X-ray beams on the condenser. Here, we demonstrate that a refractive axicon can be used as an X-ray beam shaper to match the ring-shaped aperture of the condenser. Aiming at more efficient use of the incoming X-ray intensity, we explore several axicon designs both analytically and with numerical simulations. The axicons were produced by two-photon polymerization 3D printing on thin silicon nitride membrane substrates. The first characterization of the axicon was carried out at the TOMCAT beamline of the Swiss Light Source (Switzerland).

47 OTHER INSTRUMENTATION↗

Evaluation of Desiccation Behavior in Basalt Microfiber–Reinforced Bentonite Clay for Geological Repositories of Nuclear Spent Fuel Using Digital Image Correlation

ABSTRACT Secure storage of nuclear spent fuel (NSF) is of great concern for protecting public health and safety. The preferred long-term solution is underground containment in geological repositories, where one or more engineered barrier materials (EBM) encapsulate the NSF and separate it from the natural rock. Bentonite clay is commonly used as an EBM due to its many advantageous properties including low hydraulic conductivity, which ensures limitation of water infiltration to the system and the subsequent risk of corrosion in NSF canisters. However, bentonite clay subjected to heating from nuclear decay may form desiccation cracking. This study conducted disk-shaped free shrinkage tests and ring-shaped restrained shrinkage tests of bentonite clay samples reinforced with basalt microfibers. Digital image correlation was used as a noncontact full-field displacement measurement to track the time-evolving shrinkage and desiccation cracking phenomena and make quantified comparisons between plain bentonite and bentonite with varying contents of basalt microfibers (i.e., 0.0, 0.5, 1.0, and 1.5 % wt.). Results indicate that plain bentonite and basalt microfiber-reinforced samples showed similar free shrinkage behavior, while desiccation cracking behavior was significantly altered by adding basalt microfibers. Microfiber reinforcement effectively reduced major cracks through a “crack-bridging” effect while causing minor cracks to initiate earlier and at higher moisture contents than plain bentonite. Results infer that reinforcing plain bentonite with inorganic microfibers can potentially control desiccation cracking, leading to safer and improved nuclear waste management.

Materials Science↗

Robust Carbon Dioxide Plume Imaging Using Joint Tomographic Inversion of Seismic Onset Time and Distributed Pressure and Temperature Measurements (Final Report)

We develop and demonstrate rapid and cost-effective methodologies for spatiotemporal tracking of CO2 plumes during geologic sequestration using joint inversion of seismic data and distributed pressure and temperature measurements. Key elements of our methodology are: (a) a computationally efficient approach to pressure and temperature propagation, (b) analysis of time lapse seismic data using a novel ‘seismic onset time’ approach to detect fluid front propagation, and (c) data assimilation and uncertainty assessment via joint inversion of pressure, temperature and time lapse seismic data, and (d) validating the numerical tomographic inversion using a CO2 injection demonstration projects, specifically data collected from the from the Petra Nova Parish Holdings CCUS project in the West Ranch Field, Texas and the Chester-16 reef CO2 injection site in Northern Michigan which is part of the DOE Midwestern Carbon Sequestration Project. The research team is led by Texas A&M University and includes Battelle as a subcontractor with support from Shell, Anadarko, Chevron and JX Nippon. A carbon dioxide (CO2) water-alternating-gas (WAG) pilot was conducted to gain insights into tertiary oil recovery potential via CO2 flood in the West Ranch Field as part of the Petra Nova project, the world’s largest post-combustion CO2 capture and utilization initiative. With a fluvial formation geology and large contrasts in permeability, this is a challenging and novel application of CO2 enhanced oil recovery (EOR). We build a predictive dynamic model of the subsurface that incorporates the multiphase and compositional data acquired during the pilot operation. The calibrated model is used for the carbon dioxide plume imaging. The study began with an initialization of the pilot sector model extracted from a calibrated full-field model. The pilot model calibration follows a two-step hierarchical workflow. First, we performed a large-scale update of the permeability distribution by integrating available bottomhole pressure and multiphase production data. In the second step, local permeability field is fine-tuned using a streamline-based method to match CO2 breakthrough times at the producers. The predictive capability of the calibrated model was verified through two blind validation tests: (1) the model showed good agreement with saturation logs acquired at two observation wells; and (2) the model reproduced the CO2 recovery as a fraction of the injected CO2. The use of seismic onset times has shown great promise for integrating near-continuous seismic surveys for updating geologic models. In this study, we analyze the impact of seismic survey frequency on the onset time approach aiming to extend the application of onset time to infrequent seismic surveys. In addition, we quantitatively examine the nonlinearity of the onset time method and compare it to the commonly used amplitude inversion method. We carry out a sensitivity analysis of seismic survey frequency based on the complete seismic survey data (over 175 surveys) of steam injection in a heavy oil reservoir (Peace River Unit) in Canada. Our results show that an adequate onset time map can be obtained from the infrequent seismic surveys by interpolation between seismic surveys as long as there is no change in the dominant underlying physics between the successive surveys. The study also shows that nonlinearity of the onset time method can be -smaller than that of the amplitude inversion method by several orders of magnitude. Application to the Brugge benchmark case shows that the onset time method obtains comparable permeability update as the traditional seismic amplitude inversion method with faster computation and improved convergence characteristics. We extend the streamline-based data integration approach to incorporate distributed temperature sensor (DTS) data using the concept of thermal tracer travel time. Then, a hierarchical workflow composed of evolutionary and streamline methods is employed to jointly history match the DTS and pressure data. Finally, CO2 saturation and streamline maps are used to visualize the CO2 plume movement during the sequestration process. The hierarchical workflow is applied to a carbon sequestration project in a carbonate reef reservoir within the Northern Niagaran Pinnacle Reef Trend in Michigan, USA. The monitoring data set consists of distributed temperature sensing (DTS) data acquired at the injection well and a monitoring well, flowing bottom-hole pressure data at the injection well, and time-lapse pressure measurements at several locations along the monitoring well. The history matching results indicate that the CO2 movement is mostly restricted to the intended zones of injection which is consistent with an independent warm-back analysis of the temperature data. In addition to employing simulation models and inverse methods for CO2 plume imaging, we also initialized a data-driven technology for detecting inter-well connectivity based on production and pressure data. Our machine-learning framework is built on the statistical recurrent unit (SRU) model and interprets well-based injection/production data into inter-well connectivity without relying on a geologic model. We test it on synthetic and field-scale CO2 EOR projects utilizing the water-alternating-gas (WAG) process. The validation of the proposed data-driven inter-well connectivity assessment is performed using synthetic data from simulation models where inter-well connectivity can be easily measured using the streamline-based flux allocation. The SRU model is shown to offer excellent prediction performance on the synthetic case. Despite significant measurement noise and frequent well shut-ins imposed in the field-scale case, the SRU model offers good prediction accuracy, the overall relative error of the phase production rates at most producers ranges from 10% to 30%. It is shown that the dominant connections identified by the data-driven method and streamline method are in close agreement. Texas A&M University, the lead organization in the project, was primarily responsible for the development of tomographic approaches for CO2 plume mapping in conjunction with distributed pressure, temperature and seismic onset time data. Battelle, as a subcontractor, was primarily responsible for the development of analytical and empirical methods for analyzing transient injection rate and pressure data from point/line sources such as injection and monitoring wells. An additional area of emphasis for Battelle was the use of machine learning for such tasks as inferring reservoir connectivity information from injection-production data, and identifying variable importance for machine learning-based proxy models developed from full-physics simulations. The two organizations also collaborated on the application of the tomographic inversion methodology for a field data set.

02 PETROLEUM↗

Benchmarking Microscale Ductility Measurements (Final Report of the Project DE-NE0008799)

Conventional macroscale experimentation is generally considered to be straightforward with few limitations. Conversely, micro/nanoscale experimentation presents numerous challenges in loading device design, sample preparation and handling, as well as accurate understanding of grain size and local texture effects on recorded measurements. Despite these challenges, nanopillar compression, MEMs based micro-tension, and nanoindentation approaches have been able to provide fundamental contributions to the understanding of material behavior at small lengthscales. However, the overarching shortcoming of these micro/nanoscale experimentation approaches, is the inability to directly translate measurements evaluated at the nm and µm length scales (e.g., hardness) to macroscale tensile material behavior (i.e., elastic modulus, yield strength, and ductility). The objectives of the proposed study are, 1) to establish best practices for obtaining tensile microscale ductility measurements, and 2) to validate methodologies to for comparing microscale ductility measurements to macroscale ductility measurements. In order to achieve these objectives, a multi-lengthscale, multi-temperature testing protocol and simulation framework are executed first on copper as a model material to validate the following approach, and second on reactor grade Zircaloy-2. Experiments are conducted on specimens extracted from the same test piece to ensure nominally identical grain size and texture from specimens to specimen. Motivated by the need to isolate the contribution of size-effects on obtained mechanical property measurements, specimens are manufactured with thicknesses at the micro- (1-10 µm), meso- (10-100's µm), and macroscales (sub-sized ASTM E8). In-situ full-field deformation techniques (scanning electron microscopy (SEM) grid methods and optical DIC) are incorporated into testing at each specimen length-scale to capture plasticity localization and evolution. Experimental testing for all specimens is conducted at both room temperature and elevated temperatures to probe the role of thermal activation on plastic deformation accommodation processes. Simulation efforts focus on examining the mechanical behavior of microscale specimens using a finite element approach with explicitly resolved grain morphologies, and an embedded crystal plasticity model. The cost-efficient implementation method allows for the modeling of a statistically significant number of both real (i.e., digital twin) and generated microstructures to obtain an understanding of the interrelationships between specimen microstructure and geometric variables (grain size, texture, specimen geometry, etc.) on microscale mechanical behavior.

36 MATERIALS SCIENCE↗

Solute effects upon dislocation motion and recovery in Mg alloys (Final Report)

The objective the research was to develop a firmer understanding of the interactions between substitutional solute atoms and dislocations within Mg alloys. These interactions govern the absolute and relative mobilities of various types of dislocations in Mg (e.g., basal < a>, non-basal < a>, and pyramidal < c+a>). Furthermore, they have an impact on dislocation recovery processes (e.g., cross-glide, climb, rearrangement, and annihilation). Ultimately, solute-dislocation interactions strongly impact a) strain hardening, b) strain rate sensitivity, c) plastic anisotropy of textured polycrystals, d) texture evolution, e) dislocation substructures evolution, and even f) recrystallization behavior. It is broadly known that Mg alloys exhibit poor low temperature formability. Since formability is largely governed by properties (a) – (c) in this list, it is critical to develop a better understanding of solute-dislocation interactions, if one hopes to improve the situation. While the theory of static solute strengthening is well developed, especially for alloys with face centered cubic (FCC) crystal structures, there are outstanding questions related to applications to hexagonal close packed (HCP) crystal structures and dynamic strain aging (DSA) in materials of various structures. DSA has far-reaching implications for metal formability and, in the case of Mg alloys, appears to correlate with the so-called rare-earth (RE) texture which has been shown to benefit formability. The Portevin-Le Chatelier (PLC) effect, and associated negative strain rate sensitivity, occur at higher temperatures (>100℃) in Mg alloys as compared with similar Al alloys, even though they have similar melting points and solute diffusivities. Our preliminary research has shown that modern, physics-based models of DSA can be tuned to describe the behavior of Mg alloys if a rather higher activation enthalpy is assumed for cross-core diffusivity. While this partially explains the delay in DSA to higher temperatures, it is also hypothesized that this delay is due in part to the intrinsically more thermally activated (rate sensitive) nature of non-basal < a> dislocation motion which is required for macroscopic flow of Mg alloys, whereas octahedral slip in many FCC metals like aluminum is essentially athermal at room temperature. It was originally proposed to employ a combination of in-situ diffraction-based experimental characterization to validate existing theory. It was envisioned to perform in-situ transmission electron microscopy (TEM) to assess individual dislocation behavior and in-situ high-energy X-ray diffraction (HEXRD) techniques which were showing great promise for elucidating collective dislocation behavior, 2 including recovery, especially if the contributions to various forms of diffraction peak broadening (𝜂𝜂,𝜔𝜔,and 2𝜃𝜃) can be effectively integrated. Finally, it was envisioned to perform discrete dislocation dynamics (DDD) modeling approaches to aide in the interpretation of both TEM and HEXRD experiments. In the end, mechanical tests were performed on more complex Mg alloys which exhibited evidence of dynamic strain aging, and this led to the establishment of another project. Mechanical test data obtained at McMaster University served as the basis of an assessment of the applicability of Bazinski’s “stress equivalence” theory of solute strengthening to polycrystalline alloys of Mg. Although we did not succeed in applying the approach to Mg alloys, we did develop expertise with the HEXRD approach using a BCC, β-Ti alloy and demonstrated numerous new capabilities that may be applied to any polycrystalline material in collaboration with researchers at CHESS and around the world: (1) assessment of details of the elastoplastic transition (yielding) using a combination of HEXRD and full-field polycrystal plasticity modeling, (2) the first-ever experimental observation of strong stress rotation within the individual grains of a polycrystalline material, and (3) a comprehensive analysis of the grain-level dislocation density evolution based upon diffraction peak broadening along 𝜂𝜂,𝜔𝜔,and 2𝜃𝜃 directions. This final aspect allowed us to confirm that dislocations were gliding on multiple plane types and not restricted to {110} type planes, and it also provided clues as to why some grains were unloading during straining, with surprising implications for our understanding of the effects of geometrically necessary dislocations (GNDs). Finally, graduate student, Mohammed Shabana, developed a MATLAB code which confirmed the conclusions of Prof. Catalin Picu (Rensallear Polytechnic Institute, RPI) regarding the effect of solute-trapped, forest dislocations on the breaking stress of Lomer lock junctions in FCC metal alloys. He applied the same anisotropic line-tension model to a variety of dislocation junction configurations and found an inconsistency in the widely cited results of Dupuy and Fivel regarding the Hirth Lock, and he outlined an approach to extend these finding to HCP Mg alloys that we are still pursuing with discretionary fundings at UVA.

36 MATERIALS SCIENCE↗

Enhanced convection for higher module and system efficiency (Final Technical Report)

The goal of this work is to develop new solar photovoltaic (PV) modules and solar system-scale designs that promote an increase of the convective heat transfer coefficient of at least 40%, reducing the operating temperature of the solar PV panels, and leading to a boost on the annual energy yield by at least 5%. By achieving this goal, we further foresee decreasing solar panel degradation by +0.3%/year, reducing the LCOE by 2.9-4.5¢/kWh. Currently, this question and/or approach is not considered though it is a relevant factor in performance/efficiency. The technology proposed will be at a mature stage by the end of the project and could be aligned to be proven on a full-field test. PV users, installers and manufacturers will be interested in the technology due to the aforementioned benefits. The work was conducted through a suite of wind tunnel experiments, numerical simulations and field experiments. These were then evaluated through a correlation encompassing the entire data set and implementing findings into NREL SAM. The work resulted in eight peer reviewed publications with three others in preparation. The key results of the project are 1.) Showing the capability of increasing the solar farm power production by increasing the heat transfer coefficient; 2.) Providing strategies for achieving this through large scale arrangement manipulation (heights, spacings, angles), flow deflectors and vortex generators; 3.) Demonstrating scalability of results by comparing numerical simulations, wind tunnel experiments and field experiments; 4.) Providing a correlation to be able to capture the various arrangements while exploiting the heat transfer coefficient; and 5.) Proving the feasibility of such strategies and their capacity through a comprehensive techno-economic analysis. For the heat transfer coefficients, gains of 40% were attained and increases in power in the excess of 5% were observed. For heat mitigation strategies, consideration at the PV array scale were inter-panel spacing, panel height and panel angles while at the module level, wind breaks such as solid walls and selected gaps were evaluated as well as vortex generators on the panel. All proposed strategies saw an increase in heat transfer coefficient. For scalability, results were compared against field experiments and a correlation for the heat transfer quantity was attained leading to a parametric understanding of the variations and impact in the heat transfer coefficient independent of tools implemented (wind tunnel/numerical simulations). The strategies were also further assessed by a technoeconomic analysis finding improvements based on LCOE with respect to location (several states with varies atmospheric conditions/geographic locations) yielded different results. Improvements on this end were also observed especially towards the colder climates and when considering bi-facial modules. The quantitative demonstrations of such conclusions are provided in the main portion of this document with complementary results and their narration. Milestones and go/no-go points were successfully met in the project.

14 SOLAR ENERGY↗

Chapter 4: Physically informed deep learning networks for simulating microstructure evolution of 3D polycrystals

As discussed in the previous chapter, high energy diffraction microscopy (HEDM) is used to study the micromechanical evolution of a material during in situ loading. HEDM experiments have been used to verify crystal plasticity (CP) simulations [119, 91, 90, 120], for experimental planning, material design, and to further analyze experimental results. However, Fast Fourier transform-based CP (CP-FFT) or finite element-based CP (CP-FE) methods are often too slow to be used in real-time during an experiment. CP-FFT is faster than CP-FE simulations due to the absence of meshing, but can still take hours to simulate the response of a single volume depending on the size and number of strain steps [127]. Reducing computation time would create a larger exploration space in planning and design, and enable faster analysis of experimental results and real-time feedback during an experiment. This research expands upon previous works to develop a workflow for predicting the full-field evolution of a 3D polycrystal. The workflow is simplified from previous works to predict only orientation and elastic strain tensors (from which stress tensors are calculated). The network is physically informed through loss functions and network architecture for a more robust model. The orientation predictions are informed about the cubic crystal symmetry of the material by incorporating disorientation and misorientation information into the network architecture and loss. The Von Mises stress is used to enforce the correct stress-strain trends in the strain tensor predictions. Additional total strain steps from the elastic and elastoplastic region are included to better capture the stress-strain evolution at smaller total strain steps. Material and hardening parameters are additional inputs into the networks to further inform the network and to study the network’s ability to predict different materials other than those used for training.

36 MATERIALS SCIENCE↗

Exploring the process-structure-property relationship of Aerosol Deposition to phosphor coatings for non-contact thermometry

Full-field, multi-measurand diagnostics provide rich validation data necessary to improve the product life cycle time of nuclear safety components. Thermophosphor digital image correlation (TP+DIC) is a method of simultaneously measuring strain and temperature fields using patterned phosphor coatings deposited with aerosol deposition (AD). While TP+DIC produces a functional diagnostic, the coating’s reproducibility and the effect of the patterned features on the inferred temperature remains uncharacterized. This NSR&D project provided the opportunity to study two areas: 1) the tunability and repeatability of aerosol deposition and 2) the robustness of aerosol deposition phosphor on deforming substrates. The first area explores the process-property relationship of parameters elucidating the significance of each on the coating. The second area explores the relationship between the features’ characteristics (namely thickness) and the phosphor emission and inferred temperature. Together, the results will lead to the improved accuracy and functionality of TP+DIC for qualification testing of nuclear safety components.

36 MATERIALS SCIENCE↗

Enabling Digital Image Correlation for Deployment in the Hot Fuels Examination Facility

Digital image correlation (DIC) is a non-contact image-based technique for displacement and strain measurement. It is also a full-field method that can provide excellent spatial resolution. DIC has been widely used for measurements in many styles of mechanical tests and is desirable in part because of this flexibility to measure nearly any material across a range of time and length scales. Its ability to measure localized strain for tests and in environments not accessible to other strain measurement techniques makes it an ideal option for several mechanical test cases that exist in hot cell facilities such as the Hot Fuel Examination Facility (HFEF) at the Idaho National Laboratory (INL). However, DIC has not been demonstrated at HFEF, and several barriers to implementation in the hot cell environment exist. This report outlines the possible paths for deploying DIC at HFEF, weighing the technical advantages and risks with each path. It also discusses the steps which must be taken to add this capability to the suite of post-irradiation examination (PIE) techniques available at HFEF. Finally, it concludes with recommendations for the best method to perform DIC at the facility.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Interlaced Characterization and Calibration (ICC) for Improved Computational Simulation Credibility

Accurate material characterization and model calibration are pivotal for simulations used for high-consequence engineering decisions. Current characterization and calibration methods (1) use simplified test specimen geometries and global data, (2) cannot guarantee that sufficient characterization data is collected for a specific model of interest, (3) provide only mean parameter values with no uncertainty quantification, and (4) are sequential, inflexible, and time-consuming. This work developed a new paradigm—coined Interlaced Characterization and Calibration (ICC)—which drives forward the state-of-the-art in model calibration by bringing together recent advancements into one improved workflow. The ICC paradigm (1) employs tools to efficiently use full-field data to calibrate high-fidelity material models, (2) aligns the data needed with the data collected by adopting an optimal experimental design protocol, (3) provides uncertainty metrics on the calibrated model parameters, and (4) incorporates these advances into a quasi real-time feedback loop. The ICC framework was validated synthetically with both low-fidelity and high-fidelity simulations paired with several different elastoplastic material models, and was also demonstrated experimentally with an aluminum 6061 cruciform exemplar specimen. Results showed that the ICC framework—in which Bayesian optimal experimental design actively guided the experiment— resulted in calibrations with similar or better accuracy than predetermined experiments based on subject matter expertise. Moreover, the ICC framework produced a complete model calibration— with quantified uncertainties on model parameters—in 1 week, a 5 - 10× increase in efficiency over traditional approaches. Thus, the ICC paradigm improves both the calibration process and quality, by (1) improving efficiency, which increases agility of solid mechanics modeling and enables utilization of computational simulation (CompSim) at earlier stages of the design cycle and (2) providing quantified, and in some cases reduced, parameter uncertainties, which increases confidence in model predictions and supports credible decision making.

97 MATHEMATICS AND COMPUTING↗

Introduction to Digital Image Correlation (DIC) with annotated bibliography

Digital Image Correlation (DIC) is “a non-contact means of measuring motion and deformation using digital images of the object of interest”, capable of full-field measurements over large areas. Getting started in DIC can be difficult, since there are many different techniques, and the body of literature on DIC is extensive. IEEE Xplore alone lists nearly 11,000 references, including magazines, conference proceedings, journal articles, books, etc. A Google search returns around 145 million results. The intent of this annotated bibliography is to provide an entry point for beginners and to collect references that may be useful to more advanced practitioners.

42 ENGINEERING↗

Can changes in deformation regimes be inferred from crystallographic preferred orientations in polar ice?

Creep due to ice flow is generally thought to be the main cause for the formation of crystallographic preferred orientations (CPOs) in polycrystalline anisotropic ice. However, linking the development of CPOs to the ice flow history requires a proper understanding of the ice aggregate's microstructural response to flow transitions. In this contribution the influence of ice deformation history on the CPO development is investigated by means of full-field numerical simulations at the microscale. We simulate the CPO evolution of polycrystalline ice under combinations of two consecutive deformation events up to high strain, using the code VPFFT (visco-plastic fast Fourier transform algorithm) within ELLE. A volume of ice is first deformed under coaxial boundary conditions, which results in a CPO. The sample is then subjected to different boundary conditions (coaxial or non-coaxial) in order to observe how the deformation regime switch impacts the CPO. The model results indicate that the second flow event tends to destroy the first, inherited fabric with a range of transitional fabrics. However, the transition is slow when crystallographic axes are critically oriented with respect to the second imposed regime. Therefore, interpretations of past deformation events from observed CPOs must be carried out with caution, particularly in areas with complex deformation histories.

54 ENVIRONMENTAL SCIENCES↗