Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Displacement”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Complexity of Fault Rupture and Fluid Leakage in Shale: Insights From a Controlled Fault Activation Experiment

Abstract We observed rupture growth caused by controlled fluid injections at 340‐m depth within a fault zone in the low‐permeability Opalinus Clay in the Mont Terri Underground Research Laboratory (Switzerland). The rupture mechanisms were evaluated using measurements of the three‐component borehole wall displacements and fluid pressure in two sections of the fault zone and located horizontally 3 m apart from each other. One section was set across a secondary segment of the fault and used for stepwise fluid injection intended to trigger rupture growth. The other section was set across the principal shear zone of the fault for monitoring. After stepwise pressure increase up to 5.95 MPa at injection, rupture initiated as slip activation, followed by an overall opening of the fault planes connected to the injection. After 19 s of continued injection, displacements arrived at the monitoring point on the principal shear zone. These displacements are about 2.4 times larger than in the secondary fault segment. Overall, the displacements corresponded to a normal fault activation. About 9 s after the displacement front arrived, a strong pressure increase of 4.17 MPa was measured at the monitoring point, indicating a hydraulic connection had formed along the initially very low permeability fault planes between the injection and the monitoring points. Our analyses highlight that the fault activation is consistent with the state of stress but that injection pressure must be close to the normal stress acting on the fault for permeability to be generated and for fluid leakage to occur.

Guglielmi, Yves↗
DOE OSTI2024

Validation of MARS-F modeling of plasma response to RMPs using internal measurements on DIII-D

The linear resistive plasma response model is validated against the plasma internal measurement data from DIII-D edge-localized mode (ELM) control experiments with applied resonant magnetic perturbation (RMP). Considered are DIII-D discharges where the n = 1, 2, and 3 (n is the toroidal mode number) RMP field was applied. Experimental data for the plasma boundary displacement, as well as the three-dimensional (3D) pressure perturbation in the edge pedestal region, are deduced from the vertical Thomson scattering (TS) system and the horizontal charge exchange recombination (CER) system on DIII-D. The linear response model produces results that are in reasonable quantitative agreement with the DIII-D internal measurements. The plasma boundary displacement of up to 15 mm is modeled, with the pedestal pressure perturbation reaching 3 kPa. As an important insight, the larger plasma displacement measured by the vertical TS system, as compared to that measured by the horizontal CER system, is due to the contribution from the tangential component of the plasma displacement to the former. This mixing of displacement components is also found to influence the sensitivity of the CER measurement comparisons. The results of this study provide further confidence in the linear resistive plasma response model for analyzing ELM control experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Orbital-Dependent Coulomb Drag in Electron-Hole Bilayer Graphene Heterostructures

Here, we report Coulomb drag studies in an electron-hole bilayer graphene heterostructure in a magnetic field, where the orbital, spin, and valley degrees of freedom are lifted by the combined effects of exchange interaction, Zeeman energy, and a vertical displacement field. Our device enables the application of a large vertical displacement field across both layers. In addition to the well-established strong Coulomb drag between the Landau levels with an orbital quantum number 𝑁 = 0, we observe a Coulomb drag signal between the 𝑁 = 1 Landau levels under a suitable vertical displacement field. As the vertical displacement field increases further, the Coulomb drag signal between 𝑁 = 1 Landau levels weakens, and a Coulomb drag signal emerges between the 𝑁 = 0 and 𝑁 = 1 Landau levels. These findings suggest the important roles of the orbital index and the vertical displacement field in interlayer Coulomb interaction within the quantum Hall regime of coupled bilayer systems.

Landau levels↗
DOE OSTI2025

Chiral microwave nonreciprocity demonstrated via Rayleigh and Sezawa modes supported in an Al 0.58 ⁢Sc 0.42 ⁢N/4⁢H–Si⁢C platform

Chirality plays a crucial role in the helicity mismatch between surface acoustic waves and magnetic spin waves, leading to nonreciprocal transmission of acoustic power for coupled magnetoacoustic modes. Acoustic modes with both longitudinal and shear strain exhibit elliptical particle displacements, making them chiral, and different acoustic modes can exhibit different helicities of this elliptical particle displacement. Here, we study chiral acoustic modes with different helicities supported on the same piezoelectric platform and their interaction with magnetic spin waves. Our study demonstrates that the nonreciprocal transmission of acoustic power is driven by the helicity mismatch effect and, specifically, that the handedness of the nonreciprocity is based on whether the surface acoustic wave has retrograde (Rayleigh mode) or prograde (Sezawa mode) elliptical particle displacement with respect to the propagation direction. We found the transmission nonreciprocity to be significant, with 7.3 dB/mm for the retrograde particle displacement (Rayleigh mode at 2.358 GHz) and 3.3 dB/mm for prograde particle displacement (Sezawa mode at 3.112 GHz). Furthermore, this work highlights that piezoelectric platforms can be engineered to support acoustic modes with opposite helicities to enable frequency-selective nonreciprocal radiofrequency and microwave components, such as isolators and circulators, through coupled acoustic spin wave interactions.

42 ENGINEERING↗
DOE OSTI2025

Microstructural and rheological training and memory of nanocolloidal soft glasses under cyclic shear

An intrinsic feature of disordered and out-of-equilibrium materials, such as glasses, is the dependence of their properties on their history. An important example is rheological memory, in which disordered solids obtain properties based on their deformation history. Here, in this study, we employ x-ray photon correlation spectroscopy with in situ rheometry to characterize memory formation in a nanocolloidal soft glass due to cyclic shear. During a cycle, particles undergo irreversible displacements composed of a combination of shear-induced diffusion and heterogeneous, residual strain fields. At lower shear amplitudes, the displacements resemble a random walk in which the directions in each cycle are independent of those in preceding cycles, while at high amplitude, the irreversible displacements in consecutive cycles become correlated. The magnitudes of the displacements decrease with each cycle before reaching a steady state where the microstructure has been trained to achieve enhanced reversibility even at shear amplitudes well above yielding and despite the presence of thermal fluctuations. At amplitudes below and near yielding, these decreases are monotonic, while well above yielding, they are nonmonotonic, suggesting evidence of shear banding. Accompanying this microstructural training are corresponding decreases in the dissipation during each cycle and the magnitude of the residual stress toward steady-state values. Memory of the training is revealed by measurements in which the amplitude of the shear is changed after steady state is reached. The magnitude of the particle displacements, as well as the dissipation and the change in residual stress, vary nonmonotonically with the new shear amplitude, having minima near the training amplitude, thereby revealing correlated microscopic and macroscopic signatures of memory.

Chen, Yihao [Johns Hopkins Univ., Baltimore, MD (U↗
DOE OSTI2024

Recovery of Oxpure 612C-50 Coconut Carbon Particles Using Bump Arrays

Cyanide is used to leach gold from crushed ore (solid matrix) into a gold cyanide solution. The gold is extracted from the cyanide solution by adsorption onto activated carbon. The gold extraction process occurs when the activated coconut carbon is placed into tanks that contain the gold cyanide solution either in a batch process or into a continuous flow circuit. The coconut carbon is then removed from the solution for gold recovery. The use of the mesofluidic separation system (a deterministic lateral displacement system) can significantly reduce costs and waste generated from the process of removing the coconut carbon from the solution. In the coconut carbon recovery process used at many gold mines in Nevada, the gold is attached to carbon particles that are still in the cyanide liquor. The mesofluidic system can rapidly remove the gold bearing carbon particles from the cyanide liquor quickly and with no operating costs. The benefits include: • Reduce operational costs related to filtration and hydro-cyclones • Reduced acid usage for the elution gold stripping phase as only 25% of the fluids will follow the carbon particles to the final express lane when two mesofluidic systems are used in series This alternative particle removal techniques would be advantageous. A promising technique for removing larger particles from slurries is mesofluidic separation, similar to deterministic lateral displacement arrays or “bump” arrays [1] but operates at much larger flow rates. As described by Pease, et al. [2], “Bump arrays in deterministic lateral displacement devices separate large particles from small particles using arrays of staggered posts. Large particles, defined as those with radii larger than the distance between the edge of a post and the stagnation streamline from the next downstream post, must bump toward one side of the device, whereas particles smaller than this distance slalom from entrance to exit without net lateral displacement.” Unlike filters or sieves, the posts in the arrays that cause separation do not block or occlude particle flow but work because particles go around the posts. Unlike hydrocyclones separation, separation is not driven by particle density, because gravitational forces are unimportant to mesofluidic separation. Burns, et al., and Pease, et al., have shown that particles may be separated from complex suspensions, under turbulent flow conditions, and at industrially important flow rates [3-5]. They have also recently shown that mesofluidic devices may be arranged in series to increase separation performance [6]. In this paper, we evaluate the mesofluidic system for the separation of commercial OxPure GR 612 Charred Coconut shell particles as a proof of principle test for the rapid separation using the 1500 micron cut mesofluidic separator. We first describe the experimental system and conditions. We then present the experimental results. [1] Huang, L.R., E.C. Cox, R.H. Austin, J.C. Sturm. 2004. Continuous particle separation through deterministic lateral displacement. Science, 304, 987–990. https://doi.org/10.1126/science.1094567. [2] Pease L.F., J.A. Bamberger, C.A. Burns, and M.J. Minette. 2021. Large Particle Separation from Non-Newtonian Slurries using Bump Arrays. In Proceedings of the ASME 2021 Fluids Engineering Division Summer Meeting FEDSM2021-65904, V003T08A023. New York, New York: ASME. doi:10.1115/FEDSM2021-65904 [3] Burns C.A., T.G. Veldman, J. Serkowski, R.C. Daniel, X.-Y. Yu, M.J. Minette, L.F. Pease. 2021. Mesofluidic separation versus dead-end filtration. Separation and Purification Technology, 254, 117256. https://doi.org/10.1016/j.seppur.2020.117256. [4] Pease L.F., J.E. Serkowski, T.G. Veldman, J. Williams, X.-Y. Yu, M.J. Minette, J.A. Bamberger, C.A. Burns. 2021. Can Bump Arrays Separate Particles from Turbulent Flows?. In Proceedings of the ASME 2021 Fluids Engineering Division Summer Meeting FEDSM2021-67696, V003T08A024. New York, New York: ASME. doi:10.11

mesofluidic separation, slurry, bump array↗
DOE OSTI2026

Estimating IDP Origins Using ACLED Political Violence Events: Validation with Lebanon IDP Flow Data

A key input in modeling population distribution and flows following conflict events is the inclusion of Internally Displaced Person (IDP) flows between administrative units within a country. These flows are critical for capturing population movement and redistribution driven by current events, particularly conflict. In some cases, IDP destination data are available while origin data is incomplete or unavailable. This creates a gap in understanding where displacement is occurring, limiting the ability to model population redistribution accurately. Without origin data, it is not possible to reallocate population flows or accurately represent where displacement is occurring within the country. This report evaluates whether Armed Conflict Location Event Data (ACLED) political violence event data can be used to estimate IDP origin distributions when direct origin data are unavailable. The approach is validated using historical IDP flow data from Lebanon, where both origins and destinations are observed. Results show that ACLED event distributions strongly correspond to observed IDP-origin patterns, particularly when using cumulative 60-day event windows. The method is most reliable for identifying major origin districts and approximating proportional origin shares. However, it is not intended to reconstruct exact individual displacement flows, but rather to provide a probabilistic spatial allocation of displacement origins.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2023

Patterns and rates of soil movement and shallow failures across several small watersheds on the Seward Peninsula, Alaska

Abstract. Thawing permafrost can alter topography, ecosystems, and sediment and carbon fluxes, but predicting landscape evolution of permafrost-influenced watersheds in response to warming and/or hydrological changes remains an unsolved challenge. Sediment flux and slope instability in sloping saturated soils have been commonly predicted from topographic metrics (e.g., slope, drainage area). In addition to topographic factors, cohesion imparted by soil and vegetation and melting ground ice may also control spatial trends in slope stability, but the distribution of ground ice is poorly constrained and hard to predict. To address whether slope stability and surface displacements follow topographic-based predictions, we document recent drivers of permafrost sediment flux present on a landscape in western Alaska that include creep, solifluction, gullying, and catastrophic hillslope failures ranging in size from a few meters to tens of meters, and we find evidence of rapid and substantial landscape change on an annual timescale. We quantify the timing and rate of surface movements using a multi-pronged, multi-scalar dataset including aerial surveys, interannual GPS surveys, synthetic aperture radar interferometry (InSAR), and climate data. Despite clear visual evidence of downslope soil transport of solifluction lobes, we find that the interannual downslope surface displacement of these features does not outpace downslope displacement of soil in locations where lobes are absent (downslope movement means: 7 cm yr−1 for lobes over 2 years vs. 10 cm yr−1 in landscape positions without lobes over 1 year). Annual displacements do not appear related to slope, drainage area, or modeled total solar radiation but are likely related to soil thickness, and volumetric sediment fluxes are high compared to temperate landscapes of comparable bedrock lithology. Time series of InSAR displacements show accelerated movement in late summer, associated with intense rainfall and/or deep thaw. While mapped slope failures do cluster at slope–area thresholds, a simple slope stability model driven with hydraulic conductivities representative of throughflow in mineral and organic soil drastically overpredicts the occurrence of slope failures. This mismatch implies permafrost hillslopes have unaccounted-for cohesion and/or throughflow pathways, perhaps modulated by vegetation, which stabilize slopes against high rainfall. Our results highlight the breadth and complexity of soil transport processes in Arctic landscapes and demonstrate the utility of using a range of synergistic data collection methods to observe multiple scales of landscape change, which can aid in predicting periglacial landscape evolution.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

In situ mechanical characterization of a single cell-cell adhesion interface under large strain

A method of measuring a stress-strain curve in a cell-cell adhesion interface, the method including: providing a structure including a first movable island supported by a first beam, a second movable island supported by a second beam, and a gap therebetween connected by a pair of cells forming a junction, the pair of cells comprising a cell-cell adhesion interface having an initial length defined by a distance between nuclei of the pair of cells; moving the second movable island with a defined displacement; determining a displacement of the first movable island based on moving the second movable island; calculating a difference between the displacement of the first movable island and the defined displacement of the second movable island based on moving the second movable island; determining an applied strain in the cell-cell adhesion interface between the pair of cells based on the difference divided by the initial length of the cell-cell adhesion interface; calculating a force between the cell-cell adhesion interface of the pair of cells based on the displacement of the first movable island; calculating a stress in the cell-cell adhesion interface between the pair of cells based on the force; and determining the stress-strain curve of the cell-cell adhesion interface between the pair of cells by plotting the calculated stress against the applied strain.

Yang, Ruiguo↗
DOE OSTI2026

HDG-1 Experiment Irradiation Monitoring Data Qualification Final Report

SUMMARY The U.S. Department of Energy (DOE) Advanced Reactor Technologies (ART) Graphite Research and Development (GRD) Program is conducting a series of six experiments to quantify the effects of irradiation on nuclear-grade graphite. This report documents the qualification of irradiation monitoring data for the fifth experiment, High Dose Graphite-1 (HDG-1). Qualified monitoring data are required by the ART program to support the design and licensing of the first high-temperature reactor (HTR) nuclear plant. Data are classified as Qualified if they meet the usage requirements described in the experiment planning and quality assurance (QA) documents, Failed if they do not meet those requirements and provide no usable information, or Trend if they do not fully meet all requirements but still provide useful information subject to an assessment of how any deficiencies may affect a particular use of the data. HDG-1 irradiation began with Advanced Test Reactor (ATR) Cycle 168B on August 24, 2020, and concluded after Cycle 173C on January 27, 2025. The HDG-1 capsule was removed from the reactor core twice—during core internal change (CIC) Cycle 170A and powered axial locator mechanism (PALM) Cycle 172A—to prevent overheating of the graphite specimens during high-power PALM cycles. The capsule was therefore irradiated during a total of seven normal ATR cycles: 168B, 169A, 171A, 171B, 173A, 173B, and 173C. Irradiation monitoring data evaluated in this report include thermocouple (TC) temperature, gas flow rate, gas moisture, gas pressure, specimen load, and graphite stack displacement. Temperature. A total of 14,508,065 TC temperature records were captured. Of these, 13,901,785 (95.8%) are Qualified and 606,280 (4.2%) are Failed. The principal source of failed temperature data was the instrument failure of TC-9 (Zone 2) on June 24, 2024, and TC-10 (Zone 1) on July 5, 2024, near the end of Cycle 173A, which resulted in 595,554 Failed readings. An additional 379 missing values and 10,347 slightly negative values from TC-13 during ATR outages are also Failed. Neither TC-9 nor TC-10 was used as a temperature-control TC, and their failures did not compromise capsule condition monitoring. Correlation analysis of all 13 TCs found no evidence of virtual junction formation. Control chart analysis revealed clear downward drift of approximately 80°C for TC-6 (Zone 3) relative to other stable TCs, and possible downward drift of approximately 60°C for TC-13 relative to the Zone 5 control TC (TC-1), though TC-13 remained consistent with the Zone 2 control TC (TC-12). Gas flow. A total of 20,088,090 gas flow rate records were captured. Of these, 19,941,463 (99.3%) are Qualified and 146,627 (0.7%) are Failed due to missing values. All argon, helium, and total gas flow data were within expected ranges throughout the irradiation. Gas moisture. A total of 1,116,005 outlet gas moisture values were captured. Of these, 1,101,421 (98.7%) are Qualified and 14,584 (1.3%) are Failed, comprising 14,556 out-of-range values and 28 missing values. The out-of-range moisture values exceeded 22,000 ppmv for approximately 1 week at the beginning of Cycle 173A, when accumulated moisture evaporated after the capsule was retrieved from water storage during PALM Cycle 172A and reinserted into the east flux trap. Moisture levels returned to below 25 ppmv for the remaining three cycles, and the transient high-moisture event did not affect the integrity of specimen irradiation. Gas pressure. A total of 7,812,035 gas pressure values were captured. Of these, 6,642,048 (85.0%) are Qualified and 1,169,987 (15.0%) outlet pressure values are Failed, comprising 718,537 zero outlet pressure values due to sensor failure from Cycle 168B through Cycle 171B, 54,550 missing values, and 396,900 too-low outlet pressure values, ranging from 1.1 to 1.6 psia after sensor replacement during Cycle 173A. Load. A total of 6,696,030 load values were captured. Of these, 6,694,580 (99.98%) are Qualified and 1,450 (0.02%) are Failed due to missing values. Applied loads to the six specimen stacks were stable throughout the irradiation. Stack displacement. A total of 6,696,030 displacement values were captured. Of these, 5,713,297 (85.32%) are Qualified and 3,781 (0.06%) are Failed due to missing values. Stack displacement increased consistently throughout the irradiation, reaching approximately 3.08 in. for Channels 5 and 6 by the end of irradiation. 978,952 (14.62%) substantially elevated displacements observed for Channel 6 beginning in Cycle 171A and for Channel 5 beginning in Cycle 173A are assigned Trend status. Raising pressure. A total of 1,115,999 raising pressure values were captured. Of these, 1,115,430 (99.95%) are Qualified and 569 (0.05%) are Failed due to missing values. Ram pressure. A total of 6,696,030 ram pressure values were captured. Of these, 6,692,249 (99.95%) are Qualified and 3,484 (0.05%) are Failed due to missing values. Stack raising was perf

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Utah FORGE: Fault Shear Reactivation Experimental Data for Fluid Injection-Rate Controls on Seismic Moment

Included are experimental data recorded from shear experiments that specifically explore the link between fluid-injection rate and seismic moment resulting from shear reactivation of laboratory faults. Raw mechanical data from three experiments are included alongside corresponding MATLAB scripts that import and plot the data, as well as use it to calculate shear and normal stress. Experiments are performed on 2.5-3 inch long granitoid cores from the Utah FORGE EGS demonstration site, containing a single inclined fracture with small-scale roughness added to the fracture surface. The raw data included here were recorded from an aluminum triaxial pressure vessel (TEMCO) configured with three independent servo-controlled pumps, with distilled water used as the working fluid. The pumps control confining pressure, upstream pore pressure, and axial pressure, with each pump connected to a LabView interface to record applied pressures, cumulative injected water volumes, and pump flow rates. The downstream outlet from the fracture is closed to allow pressurization, which is measured by an external pressure transducer. A linear variable differential transformer (LVDT) attached to the axial piston measures axial displacement, from which we calculate shear displacement along the fracture. Additionally, P-wave transducers are used to record acoustic signatures, where acoustic emission events and maximum amplitudes are compared against seismic moment and shear slip velocity. Fluid injection rates range between 0.05 mL/min, 0.25 mL/min, and 0.75 mL/min for each experiment. Along-fault pressure distributions are progressively less uniform as injection rates increase, representing a switch from steady-state to transient conditions. Triggered shear displacement is used as a proxy for seismic moment and is indexed against cumulative injection volume and rate. Each experiment is performed under constant shear stress conditions, and the sample is fully saturated with DI water. Axial and confining stresses are applied to 3 MPa through pressure-stepping in 500 kPa increments. The pore pressure is held constant at 200 kPa prior to initiating the experiment, and initial axial displacement is recorded. The axial stress is then increased to initiate shear mobilization during the loading phase (run-in) until a peak steady state is achieved. The initial shear stress is reduced to approximately 80% of the peak shear stress by decreasing the axial stress, then held constant for the duration of each experiment.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Effect of Stoichiometry on the Structure and Polarization of BaTiO 3

Barium titanate (BaTiO 3 ) is a material of interest for photonic device applications due to its strong optical non-linearity. However, BaTiO 3 -based devices have not found widespread adoption, in part due to the challenges associated with synthesizing high quality thin-films. Here, high-resolution scanning transmission electron microscope (STEM) imaging is used to investigate the atomic structure of both on- and off-stoichiometric BaTiO 3 synthesized by molecular beam epitaxy (MBE). Here, this investigation reveals an asymmetry in the way the BaTiO 3 atomic lattice accommodates off-stoichiometry growth and unveils features beyond what is expected from diffraction or surface characterization techniques. Excess titanium incorporates into the BaTiO 3 lattice to form pervasive defects despite titanium-rich films having a low surface roughness and high-quality appearance in diffraction. Excess barium forms a rough, water-soluble surface layer but does not significantly impact the quality of the BaTiO 3 lattice. STEM is used to map titanium atom displacement in real-space. The average displacement distance is 30–60 pm in the strained thin-films, higher than the <20 pm displacement in bulk BaTiO 3 . Additionally, the titanium atom displacement direction deviates from the c-axis of the unit cell, which may have implications for the material's electro-optic tensor and thus for electro-optic device design.

Cavanagh, Ashley E. [Harvard Univ., Cambridge, MA ↗
DOE OSTI2020

Interlaboratory Study of Digital Volume Correlation Error Due to X-Ray Computed Tomography Equipment and Scan Parameters: an Update from the DVC Challenge

Background: The quality of Digital Volume Correlation (DVC) full-field displacement measurements depends directly on the characteristics of the X-ray Computed Tomography (XCT) equipment, and scan procedures used to acquire the tomographic images. Objective: In this work, we seek to experimentally study the effects of XCT equipment and tomographic scan procedures on the quality of these images for DVC analysis, and to survey the level of DVC error that may be achieved using standard XCT operating procedures. Methods: Six participants in an interlaboratory study acquired high-quality XCT scans of a syntactic foam before and after rigid body motion. The resulting images were correlated using commercial DVC software to quantify error sources due to random image noise, reconstruction artifacts, as well as systematic spatial or temporal distortion. Results: In the absence of rigid body motion, the standard deviation of the displacement measurements ranged from 0.012 to 0.043 voxels using a moderate subvolume size, indicating that subvoxel measurement resolution could readily be achieved with a variety of XCT equipment and scan recipes. Comparison of consecutive scans without rigid body motion showed transient dilatational displacement gradients due to self-heating of the X-ray source and/or thermal expansion of the foam. Evaluation of the scans after rigid body motion showed significant, machine-specific spatial distortion in the displacement fields of up to 0.5 voxels; new approaches to remove this error need to be developed. Conclusions: Analysis of the scan protocols used in the interlaboratory study, as well as a complementary parametric sensitivity study, showed that the DVC error was strongly influenced by the XCT equipment, but could be mitigated by adjusting the total scan duration.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Propagation of Noise Uncertainty Through Virtual Strain Gauge Formulations for 2D Digital Image Correlation

The effect of displacement uncertainty is examined on 2-dimensional strain, calculated using linear surfaces fitted to the displacement field. A classical engineering error propagation method is used to calculate uncertainty in Green-Lagrangian strain calculations. The derived uncertainty is compared to a Monte Carlo simulation and discrepancies under 2% are seen between these two methods. The effect of virtual strain gauge size, displacement uncertainty, and boundaries on the region of interest on the strain uncertainty are considered. Here, an exponential decay relationship is observed between strain uncertainty and virtual strain gauge size, while a linear relationship is seen between strain and displacement uncertainty. For boundaries in the region of interest, strain uncertainty is affected by the reduced number of points available to perform the regression.

42 ENGINEERING↗
DOE OSTI2020

A coupled multipoint stress–multipoint flux mixed finite element method for the Biot system of poroelasticity

In this work, we present a mixed finite element method for a five-field formulation of the Biot system of poroelasticity that reduces to a cell-centered pressure–displacement system on simplicial and quadrilateral grids. A mixed stress–displacement–rotation formulation for elasticity with weak stress symmetry is coupled with a mixed velocity–pressure Darcy formulation. The spatial discretization is based on combining the multipoint stress mixed finite element (MSMFE) method for elasticity and the multipoint flux mixed finite element (MFMFE) method for Darcy flow. It uses the lowest order Brezzi–Douglas–Marini mixed finite element spaces for the poroelastic stress and Darcy velocity, piecewise constant displacement and pressure, and continuous piecewise linear or bilinear rotation. A vertex quadrature rule is applied to the velocity, stress, and stress–rotation bilinear forms, which block-diagonalizes the corresponding matrices and allows for local velocity, stress, and rotation elimination. This leads to a cell-centered positive-definite system for pressure and displacement at each time step. We perform error analysis for the semidiscrete and fully discrete formulations, establishing first order convergence for all variables in their natural norms. The numerical tests confirm the theoretical convergence rates and illustrate the locking-free property of the method.

42 ENGINEERING↗
DOE OSTI2024

Stress-hybrid virtual element method on six-noded triangular meshes for compressible and nearly-incompressible linear elasticity

In this paper, we present a first-order Stress-Hybrid Virtual Element Method (SH-VEM) on six-noded triangular meshes for linear plane elasticity. Here, we adopt the Hellinger–Reissner variational principle to construct a weak equilibrium condition and a stress based projection operator. In each element, the stress projection operator is expressed in terms of the nodal displacements, which leads to a displacement based formulation. This stress-hybrid approach assumes a globally continuous displacement field while the stress field is discontinuous across each element. The stress field is initially represented by divergence-free tensor polynomials based on Airy stress functions, but we also present a formulation that uses a penalty term to enforce the element equilibrium conditions, referred to as the Penalty Stress-Hybrid Virtual Element Method (PSH-VEM). Numerical results are presented for PSH-VEM and SH-VEM, and we compare their convergence to the composite triangle FEM and B-bar VEM on benchmark problems in linear elasticity. The SH-VEM converges optimally in the L 2 norm of the displacement, energy seminorm, and the L 2 norm of hydrostatic stress. Furthermore, the results reveal that PSH-VEM converges in most cases at a faster rate than the expected optimal rate, but it requires the selection of a suitably chosen penalty parameter.

42 ENGINEERING↗
DOE OSTI2020

Time-Lapse Gravity Monitoring of CO2 Migration Based on Numerical Modeling of a Faulted Storage Complex

In this study, the performance of both surface and borehole time-lapse gravity monitoring to detect CO2 leakage from a carbon storage site is evaluated. Several hypothetical scenarios of CO2 migration in a leaky fault, and thief zones at different depths at the Kimberlina site (California, USA) constitute the basis of the approach. The CO2 displacement is simulated using the TOUGH2 simulator applied to a detailed geological model of the site. The gravity responses to these CO2 plumes are simulated using forward modeling with sensors at ground surface and in vertical boreholes. Results of inversion on one scenario are also presented. The surface-based gravity responses obtained for the different leakage scenarios demonstrate that leakage can be detected at the surface in all the scenarios but the time to detection is highly variable (10 to 40 years) and dependent on the detection threshold considered. Borehole measurements of the vertical component of gravity provide excellent constrains in depth when they are located in proximity of the density anomaly associated with the presence of CO2, thus discriminating multiple leaks in different thief zones. Joint inversion of surface and borehole data can bring valuable information of the occurrence of leakages and their importance by providing a reasonable estimate of mass of displaced fluids. This study demonstrates the importance of combining multiphase flow simulations with gravity modeling in order to define if and when gravity monitoring would be applicable at a given storage site.

Time-lapse gravity monitoring, leak detection, CCS↗
DOE OSTI2024

A machine-learning-aided data recovery approach for predicting multi-material thermal behaviors in advanced test reactor capsules

Instrumented experiments conducted at test reactors are essential to the deployment of new advanced reactor systems. Designing new experiments and generating data on specific reactor conditions require significant investments in terms of both time and cost. Finite element analysis software can be used to create high-fidelity models of experiment environments in order to support the actual experiments, but computation time remains a concern in terms of applying outcomes to real-time usage of data (e.g., a digital twin [DT]). Here, the present research proposes a machine-learning (ML) aided approach to making temperature and displacement predictions based on the thickness of the outer gas gap on the experimental capsule used for in-pile demonstration of a novel new thermal conductivity probe in the Advanced Test Reactor (ATR). This capsule consisted of U10Zr fuel, a rodlet, sodium, and inner and outer capsules. Gas gaps existed between the fuel and the rodlet, and between the inner and the outer capsule. The learning data pertained to an experimental capsule's radial distributions of temperature and displacement, as obtained based on Abaqus and the physical features. For the first step of ML sequence, the temperature was predicted using three positional parameters. Next, the displacement was predicted using seven additional parameters. Each physical feature was normalized in order to be both nondimensional and standardized. The temperature and displacement predictions showed good agreement with the simulation results in all cases involving interpolation and extrapolation. Furthermore, data similarity enhancement increased the similarity between the training and the target data, thereby increasing the predictive accuracy of the ML models. In certain extrapolation cases involving limited original ML model accuracy, data similarity enhancement and data recovery was able to somewhat improve this accuracy.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗