Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mechanical discontinuity”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A generalized machine learning workflow to visualize mechanical discontinuity

Accurate detection and mapping of mechanical discontinuity in materials has widespread industrial and research applications. Herein, we developed a generalized machine-learning framework for visualizing single mechanical discontinuity embedded in material of any composition, velocity, density, porosity, and size with limited data. The proposed visualization of discontinuity requires accurate estimations of the length, location, and orientation of the embedded discontinuity by processing multipoint wave-transmission measurements. k-Wave simulator is used to create a large dataset of elastic waveforms recorded during multi-point wave-transmission measurements through materials containing single mechanical discontinuity. k-Wave simulator considers the wave attenuation, dispersion, and mode conversion in wave motion. Discrete wavelet transform (DWT) and statistical feature extraction are essential for data preprocessing prior to the data-driven model development. DWT also minimizes the effect of noise. Using hyper-parameter tuning and cross validation, gradient boosting regression can visualize the mechanical discontinuity with an accuracy of 0.85, in terms of coefficient of determination. A double-layered neural network-based regression has better performance with an accuracy of 0.95. Use of convolutional neural network converts the predictive task from a waveform processing to an image processing problem. Convolutional neural network achieved a generalization performance of 0.91. The proposed generalized workflow requires robust simulation of wave propagation, signal processing, feature engineering, and model evaluation. Sensors closest to the source and those located opposite the source are the most significant for the desired visualization. Notably, the sensors closest to the source capture the non-linear associations, whereas the sensor on the border opposite to the source capture the linear associations between the measured waveforms and the properties of the mechanical discontinuity.

42 ENGINEERING↗
DOE OSTI2022

Monitoring the propagation of mechanical discontinuity using data-driven causal discovery and supervised learning

Mechanical wave transmission through a material is influenced by the mechanical discontinuity in the material. The propagation of embedded discontinuities can be monitored by analyzing the wave-transmission measurements recorded by a multipoint sensor system placed on the surface of the material. The proposed workflow monitors the propagation of mechanical discontinuity through three stages, namely initial, intermediate, and final stages, by using supervised learning followed by data-driven causal discovery. To the end, the workflow processes the multipoint waveform measurements resulting from a single impulse source, while considering the effects of wave attenuation, dispersion and multiple wave-propagation modes due to the discontinuity and material boundaries. Among various feature reduction techniques ranging from decomposition methods to manifold approximation methods, the features derived based on statistical parameterizations of the measured waveforms lead to reliable monitoring that is robust to changes in precision, resolution, and signal-to-noise ratio of the multipoint sensor measurements. The numbers of zero-crossing, negative-turning, and positive turning in the waveforms are the strongest causal signatures of the propagation of mechanical discontinuity. Higher order moments of the waveforms, such as variance, skewness and kurtosis, are also strong causal signatures of the propagation. Finally, the newly discovered causal signatures confirm that the statistical correlations and conventional feature rankings are not always statistically significant indicators of causality.

42 ENGINEERING↗
DOE OSTI2021

Characterization of mechanical discontinuities based on data-driven classification of compressional-wave travel times

Wave propagation and diffusive transport phenomena are influenced by the mechanical discontinuities in material. Here this study shows that certain bulk properties of the network of low-velocity mechanical discontinuities (e.g. air-filled cracks) in a material can be characterized by processing compressional-wave travel times using traditional data-driven classification techniques. To that end, we perform three tasks in chronological order: (1) use the discrete fracture network (DFN) method to create two-dimensional (2D) numerical models of crack-bearing material embedded with various types of low-velocity mechanical discontinuities, (2) use the fast marching method (FMM) to simulate the propagation of the wave/diffusion front from a single source through the 2D crack-bearing material to multiple receivers placed on the boundary of the material, and (3) train 9 data-driven classifiers to characterize the crack-bearing materials (i.e. bulk properties of the network of mechanical discontinuities in the crack-bearing material) by learning from the simulations of travel times detected by multiple receivers placed around the crack-bearing material. The classifiers identified the orientation, spatial distribution, and dispersion of the low-velocity mechanical discontinuities. Voting classifier performs the best among the 9 classifiers. For the characterization of bulk dispersion and distribution of discontinuities, the sensors located on the adjacent boundaries are more important; whereas for the characterization of bulk orientation of discontinuities, the sensors located on the opposite side are more important.

42 ENGINEERING↗
DOE OSTI2024

Quantification of the Crack Evolution Process by Extracting Relevant Signal Components from Wave Propagation and Diffusive Transport Front Measurements

Wave propagation and diffusive transport phenomena in a geological rock sample undergoing crack evolution process are expected to interact with the mechanical discontinuities in the medium. The measurements of the signals associated with these phenomena can be used to assess and monitor the crack-driven micromechanical alterations in the rock. Different wave/diffusion phenomena, such as sonic propagation, pressure diffusion, and acoustic emission (AE), are sensitive to different elements of the mechanical discontinuities generated during the evolution of the crack clusters from initiation to coalescence. Sonic propagation, AE, and pressure diffusion monitoring have the potential to map the crack evolution because the transmitter-receiver arrays can be designed, arranged and tuned to (1) achieve maximum recovery of the scattered waveforms and travel times, (2) capture the later arrivals and multiple reflections, and (3) illuminate large rock volume. However, the structural/topological complexities of the mechanical discontinuities, complex distribution of the stress fields, complex mechanical alterations in media, and fluid redistribution in the crack system pose serious challenges for the detection and modeling of the crack evolution process (from here on, we will use the term ‘crack evolution process’ to mean that the crack evolution occurred under shallow crustal conditions). For purposes of accurately accounting such complexities and heterogeneities in the absence of reliable physical laws, simulation methods, and signal processing techniques, my early-career research proposal will develop and apply novel data-driven machine learning methods to: (1) extract signal components relevant to the various phases of crack evolution and (2) generate a 2D visual map of the crack evolution process.

58 GEOSCIENCES↗
DOE OSTI2026

Mechanically graded granular scaffolds for osteochondral tissue engineering

Engineered scaffolds designed to approximate the mechanical microenvironment of the osteochondral unit often address this complexity using discrete, two-phase architectures that introduce mechanical discontinuities and interfacial stress concentrations rather than a contiguous stiffness transition. To address this challenge, we created a photoannealed polyethylene glycol (PEG) granular scaffold with a spatially controlled stiffness gradient within a cell-permissive, macroporous architecture. Stiffness was dictated by photoannealing microgels using a photomask. We tuned void volume and available surface area by varying microgel diameter and tested how mesenchymal stromal cells (MSCs) interpret local mechanical environments. MSCs exhibited position-dependent differences in morphology, cytoskeletal structure, matrix deposition, and lineage-specific gene expression within the gradient scaffolds. Softer regions supported rounded cell morphology and deposition of a glycosaminoglycan-rich matrix, whereas stiffer regions promoted cell elongation, increased cytoskeletal tension, and expression of mineral-associated markers. Gradients formed from smaller microgels magnified these spatial responses by increasing cellular confinement and adhesion site availability. Disruption of actomyosin contractility eliminated these regional differences, demonstrating that MSCs rely on tension-dependent mechanotransduction to interpret the gradient. These findings reveal that coupling microgel architecture with continuous stiffness transitions provides a tractable platform to study multiscale mechanobiologic regulation and spatially guide osteochondral tissue formation.

Biological and medical sciences↗
DOE OSTI2024

Multistage Condensation Pathway Minimizes Hysteresis in Water Harvesting with Large-Pore Metal–Organic Frameworks

Metal-organic frameworks (MOFs) have emerged as promising materials for atmospheric water harvesting (AWH). Large-pore MOFs provide high water capacity, but their significant hysteresis between sorption and desorption makes them unsuitable for AWH. Co2Cl2(BTDD) is a noteworthy exception. This MOF has large, 2.2 nm diameter one-dimensional pores, and combines both record-high water capacity and minimal hysteresis, making it an excellent material for water capture in arid areas. Sorption reversibility in Co2Cl2(BTDD) has been attributed to continuous water uptake. However, the sharp adsorption/desorption in the isotherms supports a discontinuous first-order transition. Here we use molecular simulations to compute the water adsorption and desorption pathways and isotherms in a Co2Cl2(BTDD) model, to elucidate how does this MOF achieve reversibility despite its large pore size. The simulations reveal a multi-stage mechanism of discontinuous water uptake facilitated by spatial segregation of rows of hydrophilic metal sites bridged by ~1 nm hydrophobic ligands. The multi-stage mechanism breaks the barrier of capillary condensation into smaller, easier to surmount ones, resulting in a facile process despite the sharp density discontinuity between confined liquid and vapor. Our results explain why exchanging Co2+ for Ni2+ or Clfor F- in the MOF have minimal impact on the condensation and desorption pressures. On the other hand, we predict that a decrease in hydrophilicity of the MOF vertices would strongly increase the hysteresis. We expect that the relationships between spatial distribution of hydrophilic sites and hysteresis unraveled in this study assist the design of water harvesting materials with maximal capacity and reversibility.

Zaragoza, Alberto↗
DOE OSTI2026

Peridynamic modeling of cementitious materials for nuclear waste management

Radioactive and hazardous waste generated from fuel processing plants, nuclear reactors, and hospitals, requires effective management strategies. Cementitious materials are widely applied for these needs, serving as structural materials, reactive barriers, or waste forms. In these applications cracking poses a significant risk to performance, driven by inconsistent shrinkage of the components and the varying strength and permeability of their interfaces. Traditional modeling approaches face challenges in representing the complex fracture behavior of cementitious materials due to the reliance on spatial derivatives and difficulties with mesh generation. Here, to overcome these limitations, we employ peridynamics, a novel continuum mechanics formulation that uses integrals to describe mechanical equilibrium, avoiding discontinuities associated with traditional methods. Through incorporation of a bi-linear softening model and quasistatics, an experimentally validated model for Portland cement concrete samples was created. Mechanical parameters, including compressive strength and elastic modulus were validated and variation due to aggregate packing was evaluated. Additionally, sensitivity analysis of the peridynamic parameters to the bulk material properties was established. The results lay the groundwork for evaluating the impact of unique conditions of cementitious waste forms that can be assessed to improve the reliability of waste management strategies.

Aggregates↗
DOE OSTI2026

Segmentation method comparison for residual fiber length measurement across tiled microscopy images

Fiber length distribution (FLD), in part, governs mechanical properties in discontinuous fiber composites, yet manual measurement methods limit the high-throughput characterization needed for materials design optimization. This study compares deep learning segmentation approaches for automated FLD measurement in large-field microscopy, evaluating how method choice affects the microstructural descriptors used in structure-property-processing relationships. A critical challenge is that high-resolution microscopy images (10,000×10,000 pixels) must be tiled for deep learning analysis, fragmenting fibers at boundaries. We demonstrate that segmentation method proves crucial for measurement accuracy. For example, instance segmentation with Slicing Aided Hyper Inference (SAHI) preserves individual fiber integrity across tiles while semantic segmentation prioritizes speed. Comparing against manual measurement of extracted carbon fibers, YOLOv11-SAHI matched manual ground truth (238 μm weighted mean) with 40x speedup (4.5 vs 167 minutes per image). U-Net provides rapid quantification although it is at the cost of reduced accuracy due only reliably measuring stand-alone fibers. Our comparative analysis reveals that instance segmentation with SAHI better preserves length measurements while semantic segmentation prioritizes speed, providing empirical guidance for method selection. The characterization provides essential inputs for mechanical property prediction models and inverse design workflows, accelerating composite materials development cycles.

Additive manufacturing↗
DOE OSTI2024

Understanding interfacial crystallization dynamics on carbon fiber reinforced polypropylene composite manufacturing

Reinforcing polymers with discontinuous fibers improves mechanical properties, such as strength and stiffness, and in some cases achieve isotropic properties, rendering them suitable for various engineering applications. Matrix materials are generally highly engineered thermosets (e.g. crosslinked epoxies), bonded to the fiber periphery by proprietary surface and sizing chemistries. Semicrystalline thermoplastic matrices are less utilized due to poor fiber-matrix bonding resulting in inefficient interfacial load-transfer in reinforced composites. However, flexibility with melt-processing or molding conditions can be leveraged to promote non-covalent interfacial bonding between matrix and fiber via crystallization of the matrix onto fiber surface. In the present study, we utilize a co-mingle chopped carbon and isotactic polypropylene fibers to form isotropic composites, tailoring interfacial immobilized matrix or interphase morphology to optimize performance through precise control of thermal processing/molding windows. Calorimetry and optical microscopy were employed to investigate the impact of carbon fiber at various volume fractions (10, 20, and 30 %) on isotactic polypropylene crystallization and mechanical performance. Variations in mechanical properties correspond to the structural evolution of the interfacial region and are correlated to underlying microstructural attributes using wide-angle X-ray scattering, thermal analysis, and low-field nuclear magnetic resonance spectroscopy. These results provide a practical framework for the manufacturing of thermoplastic matrix composites. In conclusion, the results presented provide a guide for the strategic optimization of interphase design, showcasing tailorable tensile strengths which outperform any isotactic polypropylene carbon fiber composites previously reported in literature.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Neutronics Enhanced Meshing Operations (NEMO) Users Manual V.1.0

The is a Python framework developed at Idaho National Laboratory that uses the Cubit Python to build Cubit models and export them in EXODUS II format. These EXODUS mesh files can be used by any -based application, but NEMO is specialized for coupled neutronics calculations. Therefore, the primary purpose of NEMO is to prepare mesh files for reactor geometries and generate a homogenization overlay file for the Serpent Monte Carlo code. This homogenization overlay file ensures the consistency of the mesh with the generation of homogenized cross sections and reference fluxes for the Superhomogenization method. NEMO also supports the preparation of mesh files for thermal and thermo-mechanical models, particularly the discontinuous finite element gap heat transfer model in Griffin. This manual introduces the user to the basic capabilities in NEMO for prismatic geometries.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Heterogeneous microstructural evolution during hydrodynamic penetration of a high-velocity copper microparticle impacting copper

Microparticle hydrodynamic penetration (HDP) may be associated with the erosion regime in cold spray processing and other high-velocity impact events. Here, in an experimental approach where we can individually launch particles and study the impact sites, we explore copper microparticles impacted on copper substrates at velocities above 900 m/s where HDP begins. We lift cross-sectional lamellae from the impact sites with a focused-ion beam for further microstructural characterization using electron backscatter diffraction and scanning transmission electron microscopy. Due to the gradients of strain, strain rate, and temperature associated with HDP, heterogeneous microstructures result. The structural evolution processes observed include deformation twinning and multiple dislocation-mediated grain recrystallization mechanisms—geometric dynamic recrystallization (gDRX), discontinuous DRX (dDRX), and meta DRX (mDRX). The higher strains at the interface lead to the most significant structural changes and complex mechanisms. In contrast, there is a gradient to more conventional dislocation plasticity away from the interface (on either the particle or substrate side). Here, these microstructural observations are consistent with the deformation map for copper and extend the observations of impact-induced recrystallization across new regimes of behavior.

Cold spray process↗
DOE OSTI2026

3D seismic imaging of a fracture damage zone controlling reservoir compartmentalization at the Raft River EGS using multi-azimuth walkaway VSP

Accurate imaging of steeply dipping fracture zones in crystalline enhanced geothermal systems (EGS) is critical for constraining permeability architecture and guiding stimulation design. However, such structures remain poorly resolved by conventional surface seismic methods. We present a fully three-dimensional (3D) elastic-waveform inversion-migration workflow applied to multi-azimuth walkaway vertical seismic profiling (VSP) data acquired in a deviated borehole at the Raft River EGS. The workflow integrates first-arrival traveltime tomography, multi-scale elastic waveform inversion (EWI), and elastic least-squares reverse-time migration (ELSRTM) to recover high-resolution compressional-and shear-wave velocity models and to image structural discontinuities in the crystalline basement. The results reveal a laterally continuous low-velocity anomaly, with shear-wave velocity reductions of 25-30%, consistent with fractureinduced mechanical weakening. Two steeply dipping discontinuities bound a 50-80 m wide fracture damage zone. Independent constraints from microseismic clustering and geochemical compartmentalization corroborate the geometry and structural significance of this feature. Synthetic modeling further confirms that structures of this scale are resolvable (~30 m). These findings indicate that the Narrows structure is a distributed fracture damage zone rather than a discrete fault plane. By resolving fracture-zone geometry at the tens-of-meters scale using a single borehole, this workflow provides a practical and transferable approach for improving structural characterization, reducing uncertainty in permeability architecture, and supporting reservoir modeling and stimulation design in fractured crystalline EGS reservoirs.

58 GEOSCIENCES↗
DOE OSTI2022

Implicit shock tracking for unsteady flows by the method of lines

Here, a recently developed high-order implicit shock tracking (HOIST) framework for resolving discontinuous solutions of inviscid, steady conservation laws is extended to the unsteady case. Central to the framework is an optimization problem which simultaneously computes a discontinuity-aligned mesh and the corresponding high-order approximation to the flow, which provides nonlinear stabilization and a high-order approximation to the solution. This work extends the implicit shock tracking framework to the case of unsteady conservation laws using a method of lines discretization via a diagonally implicit Runge-Kutta method by “solving a steady problem at each timestep”. We formulate and solve an optimization problem that produces a feature-aligned mesh and solution at each Runge-Kutta stage of each timestep, and advance this solution in time by standard Runge-Kutta update formulas. A Rankine-Hugoniot based prediction of the shock location together with a high-order, untangling mesh smoothing procedure provides a high-quality initial guess for the optimization problem at each time, which results in rapid convergence of the sequential quadratic programing (SQP) optimization solver. This method is shown to deliver highly accurate solutions on coarse, high-order discretizations without nonlinear stabilization and recover the design accuracy of the Runge-Kutta scheme. We demonstrate this framework on a series of inviscid, unsteady conservation laws in both one- and two- dimensions. We also verify that our method is able to recover the design order of accuracy of our time integrator in the presence of a strong discontinuity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Peridynamics and surrogate modeling of pressure-driven well stimulation

In this work we use the peridynamics theory of solid mechanics to simulate fracture in an annular rock domain subject to an in-situ stress and create surrogate models that predict the area of the resulting cracks. Peridynamics is a non-local formulation of continuum mechanics that naturally accommodates material discontinuities. Furthermore, unlike other fracture modeling techniques there is no need to provide information about the crack path. We utilize the peridynamics code Peridigm and take a two-stage approach to fracture modeling. First an implicit solve is performed to compute the in-situ stress state. We then execute an explicit solve where a pressure loading designed to emulate fluid-driven hydraulic fracture is applied at the borehole and transmitted to the pre-stressed rock. We present results from polynomial and single and multi-level Gaussian process surrogate models constructed from a sampling study of the peridynamics model. The surrogates predict crack area given a measure of the in-situ stress anisotropy and rise time and amplitude of the pressure loading. These surrogates take a minuscule fraction of peridynamics model's running time to evaluate and are a step towards enabling advanced optimization and uncertainty quantification workflows that require many model evaluations.

42 ENGINEERING↗
DOE OSTI2023

Natural Fracture Pattern Development Using 3D Coupled Numerical Model: Mechanics and Diagenesis

In the earth's crust, while tectonic processes can mechanically create natural fractures, cementation can simultaneously grow on the surfaces of the fractures. This crystallization can affect the growth of the fractures. To more accurately address the modeling of natural fracture development, both mechanical deformation and chemical diagenesis should be carried out in an integrated manner. Here, we achieve this with a three-dimensional (3D) simulator that couples mechanical and chemical processes for natural fracture pattern development in layer confined situations. Tectonic loading over geologic time is mimicked by displacement driven boundary conditions in a 3D displacement discontinuity method (DDM) scheme. Quartz diagenesis competes with mechanical fracture opening, resulting in three different cementation conditions - open, bridging or fully infilled – which define fracture element stiffness that is used as feedback within the 3D-DDM code during propagation. An important milestone of this work is that this coupled structural diagenetic model can reproduce power-law aperture population statistics which cannot be generated using mechanics alone. More accurately estimated mechanical fracture apertures and their cementation condition is a critical step toward better permeability estimation and fluid path analysis for subsurface reservoirs.

58 GEOSCIENCES↗
DOE OSTI2023

Rapid grain refinement and compositional homogenization in a cast binary Cu50Ni alloy achieved by friction stir processing

Friction stir processing (FSP) has been increasingly adopted for joining and processing materials in automotive, aerospace and industrial construction. During FSP a dynamic competition between high-speed shear deformation and deformation-induced heating brings about a complex competition between multiple dynamic microstructural evolution mechanisms making it difficult to predict the microstructural evolution pathway. Hence, improved understanding of microstructural evolution mechanisms during FSP can be beneficial for continued growth in the adoption of FSP for demanding applications of future. Towards this goal, this study uses a model binary Cu – 50 at.% Ni alloy to clarify the effect of single and double pass FSP on the microstructural evolution of a coarse grained and compositionally heterogeneous cast microstructure. High energy synchrotron X-ray diffraction, electron backscatter diffraction, and nanoindentation are used to clarify the microstructural evolution due to FSP. The process of compositional homogenization of as-cast segregations is studied by energy dispersive spectroscopy and atom probe tomography. Our results show that a single fast FSP pass at 30 mm.s -1 produces a 100 µm deep layer of submicrometric and hall-petch hardened CuNi grains. The initial cast compositional heterogeneities in a micrometric scale is rapidly transformed to nano-sized domains, mainly confined at grain boundaries. Double pass FSP increases the penetration depth of the processed layer and leads to a 2.9 times grain growth relative to single pass FSP. Grain fragmentation, discontinuous dynamic recrystallization, grain growth, and twinning mechanisms are discussed. In conclusion, these results highlight the value of FSP for ultrafast grain refinement and compositional homogenization of cast alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A numerical extension of the spatially-filtered Euler equations for contact discontinuities

Solving the spatially-filtered Euler equations with a kinetic energy-preserving centered discretization is an efficient framework for simulating shock-dominated flows. Here, the present work describes a numerical treatment of the spatially-filtered Euler equations to minimize oscillations when simulating contact discontinuities. To counteract dispersive errors inherent to centered schemes, a WENO-like correction term is applied to the enthalpy transport in the energy equation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗