Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “generalizability”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A Nonaqueous Redox‐Matched Flow Battery with Charge Storage in Insoluble Polymer Beads**

Abstract We describe the nonaqueous redox‐matched flow battery (RMFB), where charge is stored on redox‐active moieties covalently tethered to non‐circulating, insoluble polymer beads and charge is transferred between the electrodes and the beads via soluble mediators with redox potentials matched to the active moieties on the beads. The RMFB reported herein uses ferrocene and viologen derivatives bound to crosslinked polystyrene beads. Charge storage in the beads leads to a high (approximately 1.0–1.7 M) effective concentration of active material in the reservoirs while preventing crossover of that material. The relatively low concentration of soluble mediators (15 mM) eliminates the need for high‐solubility molecules to create high energy density batteries. Nernstian redox exchange between the beads and redox‐matched mediators was fast relative to the cycle time of the RMFB. This approach is generalizable to many different redox‐active moieties via attachment to the versatile Merrifield resin.

25 ENERGY STORAGE↗
DOE OSTI2020

Application of multidimensional structural characterization to detect and describe moderate forest disturbance

The study of vegetation community and structural change has been central to ecology for over a century, yet the ways in which disturbances reshape the physical structure of forest canopies remain relatively unknown. Moderate severity disturbances affect different canopy strata and plant species, resulting in variable structural outcomes and ecological consequences. Terrestrial lidar (light detection and ranging) offers an unprecedented view of the interior arrangement and distribution of canopy elements, permitting the derivation of multidimensional measures of canopy structure that describe several canopy structural traits (CSTs) with known links to ecosystem function. We used lidar-derived CSTs within a machine learning framework to detect and describe the structural changes that result from various disturbance agents, including moderate severity fire, ice storm damage, age-related senescence, hemlock woolly adelgid, beech bark disease, and chronic acidification. We found that fire and ice storms primarily affected the amount and position of vegetation within canopies, while acidification, senescence, pathogen, and insect infestation altered canopy arrangement and complexity. Only two of the six disturbance agents significantly reduced leaf area, counter to common assumptions regarding many moderate severity disturbances. While findings are limited in their generalizability due to lack of replication among disturbances, they do suggest that the current limitations of standard disturbance detection methods—such as optical-based remote sensing platforms, which are often above-canopy perspectives—limit our ability to understand the full ecological and structural impacts of disturbance, and to evaluate the consistency of structural patterns within and among disturbance agents. A more broadly inclusive definition of ecological disturbance that incorporates multiple aspects of canopy structural change may potentially improve the modeling, detection, and prediction of functional implications of moderate severity disturbance as well as broaden our understanding of the ecological impacts of disturbance.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Clarifying the trophic state concept to advance macroscale freshwater science and management

For over a century, ecologists have used the concept of trophic state (TS) to characterize an aquatic ecosystem's biological productivity. However, multiple TS classification schemes, each relying on a variety of measurable parameters as proxies for productivity, have emerged to meet use‐specific needs. Frequently, chlorophyll a, phosphorus, and Secchi depth are used to classify TS based on autotrophic production, whereas phosphorus, dissolved organic carbon, and true color are used to classify TS based on both autotrophic and heterotrophic production. Both classification approaches aim to characterize an ecosystem's function broadly, but with varying degrees of autotrophic and heterotrophic processes considered in those characterizations. Moreover, differing classification schemes can create inconsistent interpretations of ecosystem integrity. For example, the US Clean Water Act focuses exclusively on algal threats to water quality, framed in terms of eutrophication in response to nutrient loading. This usage lacks information about non‐algal threats to water quality, such as dystrophication in response to dissolved organic carbon loading. Consequently, the TS classification schemes used to identify eutrophication and dystrophication may refer to ecosystems similarly (e.g., oligotrophic and eutrophic), yet these categories are derived from different proxies. These inconsistencies in TS classification schemes may be compounded when interdisciplinary projects employ varied TS frameworks. Even with these shortcomings, TS can still be used to distill information on complex aquatic ecosystem function into a set of generalizable expectations. The usefulness of distilling complex information into a TS index is substantial such that usage inconsistencies should be explicitly addressed and resolved. To emphasize the consequences of diverging TS classification schemes, we present three case studies for which an improved understanding of the TS concept advances freshwater research, management efforts, and interdisciplinary collaboration. To increase clarity in TS, the aquatic sciences could benefit from including information about the proxy variables, ecosystem type, as well as the spatiotemporal domains used to classify TS. As the field of aquatic sciences expands and climatic irregularity increases, we highlight the importance of re‐evaluating fundamental concepts, such as TS, to ensure their compatibility with evolving science.

classification↗
DOE OSTI2024

Brain‐age prediction: Systematic evaluation of site effects, and sample age range and size

Abstract Structural neuroimaging data have been used to compute an estimate of the biological age of the brain (brain‐age) which has been associated with other biologically and behaviorally meaningful measures of brain development and aging. The ongoing research interest in brain‐age has highlighted the need for robust and publicly available brain‐age models pre‐trained on data from large samples of healthy individuals. To address this need we have previously released a developmental brain‐age model. Here we expand this work to develop, empirically validate, and disseminate a pre‐trained brain‐age model to cover most of the human lifespan. To achieve this, we selected the best‐performing model after systematically examining the impact of seven site harmonization strategies, age range, and sample size on brain‐age prediction in a discovery sample of brain morphometric measures from 35,683 healthy individuals (age range: 5–90 years; 53.59% female). The pre‐trained models were tested for cross‐dataset generalizability in an independent sample comprising 2101 healthy individuals (age range: 8–80 years; 55.35% female) and for longitudinal consistency in a further sample comprising 377 healthy individuals (age range: 9–25 years; 49.87% female). This empirical examination yielded the following findings: (1) the accuracy of age prediction from morphometry data was higher when no site harmonization was applied; (2) dividing the discovery sample into two age‐bins (5–40 and 40–90 years) provided a better balance between model accuracy and explained age variance than other alternatives; (3) model accuracy for brain‐age prediction plateaued at a sample size exceeding 1600 participants. These findings have been incorporated into CentileBrain ( https://centilebrain.org/#/brainAGE2 ), an open‐science, web‐based platform for individualized neuroimaging metrics.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Emerging technologies and radical collaboration to advance predictive understanding of watershed hydrobiogeochemistry

Increasing population and resource-intensive lifestyles are driving enhanced demands for clean water, food, and energy. In parallel, land-use change, climate change, and perturbations—including drought, floods, fires, and early snowmelt—are significantly reshaping interactions within watersheds throughout the world. While watersheds are the Earth's key functional unit for assessing and managing water resources, hydrological processes in watersheds also mediate biogeochemical interactions that support terrestrial life on Earth (Kaushal, Gold, Bernal, & Tank, 2018; National Research Council, 2012). Although society is dependent upon clean water availability, tractable prediction of watershed hydrobiogeochemical behavior, including watershed response to perturbations, remains a challenge. Central to the challenge are complex, multiscale interactions between plants, microorganisms, organic matter, minerals, dissolved constituents, and migrating fluids, which occur within and across bedrock-to-canopy compartments and along extensive lateral gradients of a watershed. Several recent community reports have synthesized formidable challenges associated with watershed science and technology (AGU, 2018; Blöschl et al., 2019). Here, we discuss emerging technologies and collaboration modes that are critical for developing generalizable insights about and predictive understanding of complex watershed hydrobiogeochemical behavior, which are important for underpinning optimized natural resource management.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Plant‐mediated methane transport in emergent and floating‐leaved species of a temperate freshwater mineral‐soil wetland

Abstract Methane flux from freshwater mineral‐soil (FWMS) wetlands and its variability among sites is largely modulated by plant‐mediated transport. However, plant‐mediated transport processes are rarely resolved in land surface models and are poorly parametrized for plants commonly found in FWMS wetlands. Here, relationships between methane flux and CO 2 uptake, as well as plant conductance of methane were evaluated for three plant species and two characteristic functional types: emergent (narrow‐leaved cattail) and floating‐leaved (American lotus and water lily). We found significant but contrasting correlations between methane flux and CO 2 uptake in cattails ( r 2 = 0.51, slope = −0.16, during morning) and water lily ( r 2 = 0.32, slope = 0.064, after midday). This relationship was not significant in American lotus, showing that stomata regulation of methane fluxes is species‐specific and not generalizable across the floating‐leaved plant functional type. Conductance of methane per leaf area showed distinct seasonal dynamics across species. Conductance was similar among the floating‐leaved species (6.2 × 10 −3 m d −1 in lotus and 7.2 × 10 −3 m d −1 in water lily) and higher than conductance in the emergent species (2.7 × 10 −3 m d −1 ). Our results provide direct observations of plant conductance rates and identify the vegetation parameters (leaf area, stomatal conductance) that modify them. Our results further suggest that models of methane emissions from FWMS should parameterize plant‐mediated transport in different plant functional types, scaled by leaf area and with variable seasonal phenological dynamics, and consider possible species‐specific mechanisms that control methane transport through plants.

Villa, Jorge A.↗
DOE OSTI2026

Data‐Driven Engineering of Thermostable Collagen‐Mimetic Peptoid Triple Helices

Collagen-mimetic peptides (CMPs) are engineered molecules designed to replicate the triple-helical structure of natural collagen. A repeating x–y-Gly sequence is the defining motif of CMPs and is critical to their triple-helical structure and stability. Substitutions to the residues occupying the x and y positions present a means to modulate the CMP structure and properties. Peptoid residues—N-substituted glycine derivatives—present an attractive potential substitution due to their thermal stability, proteolytic resistance, biocompatibility, and diverse palette of non-natural side chains, but also tend to introduce a high degree of backbone flexibility that can diminish the stability of the triple helix. In this work, we report a computational active learning cycle comprising molecular dynamics simulation, Gaussian process regression, and Bayesian optimization to computationally identify a number of promising peptoid substitutions predicted to stabilize the desired quaternary structure through side chain interactions and produce stable peptoid-based collagen-like triple helices. To experimentally test the computational predictions, a top candidate identified by the screen was synthesized and imaged using scanning electron microscopy to resolve fibril-like bundles consistent with collagen-like triple helices. This work predicts a number of CMP peptoid substitutions capable of forming stable triple-helical structures, presents a generalizable design strategy for engineering desired peptoid structures, and opens new avenues for the design of peptoid-based biomimetic materials.

active learning↗
DOE OSTI2020

Voxel carving–based 3D reconstruction of sorghum identifies genetic determinants of light interception efficiency

Changes in canopy architecture traits have been shown to contribute to yield increases. Optimizing both light interception and light interception efficiency of agricultural crop canopies will be essential to meeting the growing food needs. Canopy architecture is inherently three-dimensional (3D), but many approaches to measuring canopy architecture component traits treat the canopy as a two-dimensional (2D) structure to make large scale measurement, selective breeding, and gene identification logistically feasible. We develop a high throughput voxel carving strategy to reconstruct 3D representations of sorghum from a small number of RGB photos. Our approach builds on the voxel carving algorithm to allow for fully automatic reconstruction of hundreds of plants. It was employed to generate 3D reconstructions of individual plants within a sorghum association population at the late vegetative stage of development. Light interception parameters estimated from these reconstructions enabled the identification of known and previously unreported loci controlling light interception efficiency in sorghum. The approach is generalizable and scalable, and it enables 3D reconstructions from existing plant high throughput phenotyping datasets. We also propose a set of best practices to increase 3D reconstructions’ accuracy.

3D plant reconstruction↗
DOE OSTI2022

Position-specific carbon isotope analysis of serine by gas chromatography/Orbitrap mass spectrometry, and an application to plant metabolism

Position-specific 13 C/ 12 C ratios within amino acids remain largely unexplored in environmental samples due to methodological limitations. We hypothesized that natural-abundance isotope patterns in serine may serve as a proxy for plant metabolic fluxes including photorespiration. Here we describe an Orbitrap method optimized for the position-specific carbon isotope analysis of serine to test our hypothesis and discuss the generalizability of this method to other amino acids. Position-specific carbon isotope ratios of serine were measured using a Thermo Scientific™ Q Exactive™ GC Orbitrap™. Amino acids were hydrolyzed from Arabidopsis biomass, purified from potential matrix interferences, and derivatized alongside standards. Derivatized serine (N,O-bis(trifluoroacetyl)methyl ester) was isolated using gas chromatography, trapped in a reservoir, and purged into the electron ionization source over tens of minutes, producing fragment ions containing different combinations of atoms from the serine-derivative molecule. The 13 C/ 12 C ratios of fragments with monoisotopic masses of 110.0217, 138.0166, and 165.0037 Da were monitored in the mass analyzer and used to calculate position-specific δ 13 C values relative to a working standard. This methodology constrains position-specific δ 13 C values for nanomole amounts of serine isolated from chemically complex mixtures. The δ 13 C values of fragment ions of serine were characterized with ≤1‰ precisions, leading to propagated standard errors of 0.7–5‰ for each carbon position. Position-specific δ 13 C values differed by up to ca 28 ± 5‰ between serine molecules hydrolyzed from plants grown under contrasting pCO 2 , selected to promote different fluxes through photosynthesis and photorespiration. The method was validated using pure serine standards characterized offline. Here this study presents the first Orbitrap-based measurements of natural-abundance, position-specific carbon isotope variation in an amino acid isolated from a biological matrix. We present a method for the precise characterization of isotope ratios in serine and propose applications probing metabolism in plants. We discuss the potential for extending these approaches to other amino acids, paving the way for novel applications.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Use of Longitudinal Serum Analysis and Machine Learning to Develop a Classifier for Cancer Early Detection

Early detection of solid tumors through a simple screening process, such as the proteomic analysis of biofluids, has the potential to significantly alter the management and outcomes of cancers. The application of advanced targeted proteomics measurements and data analysis strategies to uniformly collected serum or plasma samples would enable longitudinal studies of cancer risk, progression, and response to therapy that have the potential to significantly reduce cancer burden in general. In this article, we describe a generalizable workflow combining robust, multiplexed targeted proteomics measurements applied to longitudinal samples from the Department of Defense Serum Repository with a Random Forest machine learning method for developing and initially evaluating the performance of candidate biomarker panels for early detection of cancers. The effectiveness of this approach was demonstrated in a cohort of 175 head and neck squamous cell carcinoma patients. The outlined protocols include methods for sample preparation, instrument analysis, and data analysis and interpretation using this workflow.

Longitudinal analysis, machine learning, cancer, e↗
DOE OSTI2021

Energy and Emission Prediction for Mixed-Vehicle Transit Fleets Using Multi-task and Inductive Transfer Learning

Public transit agencies are focused on making their fixed-line bus systems more energy efficient by introducing electric (EV) and hybrid (HV) vehicles to their fleets. However, because of the high upfront cost of these vehicles, most agencies are tasked with managing a mixed-fleet of internal combustion vehicles (ICEVs), EVs, and HVs. In managing mixed-fleets, agencies require accurate predictions of energy use for optimizing the assignment of vehicles to transit routes, scheduling charging, and ensuring that emission standards are met. The current state-of-the-art is to develop separate neural network models to predict energy consumption for each vehicle class. Although different vehicle classes’ energy consumption depends on a varied set of covariates, we hypothesize that there are broader generalizable patterns that govern energy consumption and emissions. In this paper, we seek to extract these patterns to aid learning to address two problems faced by transit agencies. First, in the case of a transit agency which operates many ICEVs, HVs, and EVs, we use multi-task learning (MTL) to improve accuracy of forecasting energy consumption. Second, in the case where there is a significant variation in vehicles in each category, we use inductive transfer learning (ITL) to improve predictive accuracy for vehicle class models with insufficient data. As this work is to be deployed by our partner agency, we also provide an online pipeline for joining the various sensor streams for fixed-line transit energy prediction. Here, we find that our approach outperforms vehicle-specific baselines in both the MTL and ITL settings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Multi-task Parallelism for Robust Pre-training of Graph Foundation Models on Multi-source, Multi-fidelity Atomistic Modeling Data

Graph foundation models using graph neural networks promise sustainable, efficient atomistic modeling. To tackle challenges of processing multi-source, multi-fidelity data during pre-training, recent studies employ multi-task learning, in which shared message passing layers initially process input atomistic structures regardless of source, then route them to multiple decoding heads that predict data-specific outputs. This approach stabilizes pre-training and enhances a model’s transferability to unexplored chemical regions. Preliminary results on approximately four million structures are encouraging, yet questions remain about generalizability to larger, more diverse datasets and scalability on supercomputers. We propose a multi-task parallelism method that distributes each head across computing resources with GPU acceleration. Implemented in the open-source HydraGNN architecture, our method was trained on over 24 million structures from five datasets and tested on the Perlmutter, Aurora, and Frontier supercomputers, demonstrating efficient scaling on all three highly heterogeneous super-computing architectures.

Lupo Pasini, Massimiliano [ORNL] (ORCID:0000000249↗
DOE OSTI2023

Fast and accurate reduced-order modeling of a MOOSE-based additive manufacturing model with operator learning

One predominant challenge in additive manufacturing (AM) is to achieve specific material properties by manipulating manufacturing process parameters during the runtime. Such manipulation tends to increase the computational load imposed on existing simulation tools employed in AM. The goal of the present work is to construct a fast and accurate reduced-order model (ROM) for an AM model developed within the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework, ultimately reducing the time/cost of AM control and optimization processes. Our adoption of the operator learning (OL) approach enabled us to learn a family of differential equations produced by altering process variables in the laser’s Gaussian point heat source. More specifically, we used the Fourier neural operator (FNO) and deep operator network (DeepONet) to develop ROMs for time-dependent responses. Furthermore, we benchmarked the performance of these OL methods against a conventional deep neural network (DNN)-based ROM. Ultimately, we found that OL methods offer comparable performance and, in terms of accuracy and generalizability, even outperform DNN at predicting scalar model responses. The DNN-based ROM afforded the fastest training time. Furthermore, all the ROMs were faster than the original MOOSE model yet still provided accurate predictions. FNO had a smaller mean prediction error than DeepONet, with a larger variance for time-dependent responses. Unlike DNN, both FNO and DeepONet were able to simulate time series data without the need for dimensionality reduction techniques. Finally, the present work can help facilitate the AM optimization process by enabling faster execution of simulation tools while still preserving evaluation accuracy.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Chatter detection in simulated machining data: a simple refined approach to vibration data

Vibration monitoring is a critical aspect of assessing the health and performance of machinery and industrial processes. This study explores the application of machine learning techniques, specifically the Random Forest (RF) classification model, to predict and classify chatter—a detrimental self-excited vibration phenomenon—during machining operations. While sophisticated methods have been employed to address chatter, this research investigates the efficacy of a novel approach to an RF model. The study leverages simulated vibration data, bypassing resource-intensive real-world data collection, to develop a versatile chatter detection model applicable across diverse machining configurations. The feature extraction process combines time-series features and Fast Fourier Transform (FFT) data features, streamlining the model while addressing challenges posed by feature selection. By focusing on the RF model’s simplicity and efficiency, this research advances chatter detection techniques, offering a practical tool with improved generalizability, computational efficiency, and ease of interpretation. The study demonstrates that innovation can reside in simplicity, opening avenues for wider applicability and accelerated progress in the machining industry.

42 ENGINEERING↗
DOE OSTI2026

Rerouting reductant flux via protein tethering enhances biohydrogen production in Thermococcus kodakarensis

Microbes that generate copious amounts of hydrogen (H 2 ) via dark fermentation are a promising means to evolve and improve renewable biofuels. Many anaerobic hyperthermophilic archaea, such as the fast-growing, genetically tractable, heterotroph Thermococcus kodakarensis, produce generous quantities of H2 and provide an idealized platform to further optimize naturally high levels of biohydrogen reduction. Precise genetic manipulations and modifications to growth conditions have already resulted in substantial increases to H2 output but additional improvements are desired. An unexamined and potentially valuable route towards increased H 2 production is to tether select electron donor and acceptor proteins together to reroute and maximize the flow of electrons towards H 2 production. Such strategies have shown promise in Bacteria and Eukarya but have not yet been investigated in thermophilic Archaea. Here, we generate and evaluate twelve novel T. kodakarensis strains wherein a proteinaceous electron carrier (a ferredoxin, Fd) is physically tethered to the membrane-bound-hydrogenase (MBH), the sole H 2 producing enzyme, to direct electron flux towards biohydrogen generation. Growth assessments and H 2 output measurements demonstrate that strains encoding protein-fusions evolve up to ~ 40% more H2 per cell than the host strain. Eliminating H 2 consumption and alternative routes of electron sinks in concert with protein tethering further increased H2 output per cell for a maximum increase of ~ 66% over the host strain. Our results demonstrate that rerouting electron flux via protein tethering coupled with the elimination of reductant sinks is a promising means towards improved biohydrogen production in T. kodakarensis. KEY POINTS: Protein tethering between redox proteins can reroute electron flux in vivo. Enforced protein proximity results in ~ 40% increases in H2 production per cell. Protein-tethering provides a generalizable framework to redirect redox metabolism.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Integrated modeling for assessing climate change impacts on water resources and hydropower potential in the Himalayas

Regional hydroclimatic variability and change can affect water resources and hydropower generation. It is essential to assess hydropower potential under current and future climatic conditions to inform the design and operation of hydropower infrastructures. Here, we employ an integrated modeling framework to assess the impact of projected hydroclimatic conditions on water resource systems and hydropower generation. The integrated framework samples climate model outputs under different scenarios to force a hydrologic model and produces streamflow projections. The projected streamflows are inputs for the future hydropower potential assessment. We implement the framework in the central Himalayan river basin. Our results demonstrate substantial spatiotemporal variability in different water balance components (precipitation, evapotranspiration, and water yield) under current and future climatic conditions. For the Himalayan Tila river basin, the annual average energy production is expected to increase under future hydroclimatic conditions (up to 39% in Tila-2 hydropower project, suggested by ensemble mean). Furthermore, this increase in energy is driven mainly by the increased streamflow projections, particularly during the dry season and in the late century. Our results highlight the impacts of hydroclimatic variability in hydropower productions and are of practical use to provide decision-relevant information for designing and operating hydropower infrastructures. The integrated modeling framework presented here is region-specific; however, the approach is reproducible, and the overall insights are generalizable across the Himalayan region.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Integrated energy-water-land nexus planning in the Colorado River Basin (Argentina)

Integrated energy-water-land (EWL) planning promotes synergies and avoids conflicts in ways that sector-specific planning approaches cannot. Many important decisions that influence emerging EWL nexus issues are implemented at regional (e.g., large river basin, electricity grid) and sub-regional (e.g., small river basin, irrigation district) scales. However, actual implementation of integrated planning at these scales has been limited. Simply collecting and visualizing data and interconnections across multiple sectors and sub-regions in a single modeling platform is a unique endeavor in many regions. This study introduces and applies a novel approach to linking together multiple sub-regions in a single platform to characterize and visualize EWL resource use, EWL system linkages within and among sub-regions, and the EWL nexus implications of future policies and investments. This integrated planning methodology is applied in the water-stressed Colorado River Basin in Argentina, which is facing increasing demands for agricultural and fossil fuel commodities. Guided by stakeholders, this study seeks to inform basin planning activities by characterizing and visualizing (1) the basin’s current state of EWL resources, (2) the linkages between sectors within and among basin sub-regions, and (3) the EWL nexus implications of planned future agricultural development activities. Results show that water scarcity, driven in part by human demands that have historically reached 60% of total surface water supply, poses a substantial constraint to economic development in the basin. The Colorado basin has the potential to serve as a testbed for crafting novel and generalizable sub-regional EWL planning approaches capable of informing the EWL planning dialogue globally.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Towards a unified nonlocal, peridynamics framework for the coarse-graining of molecular dynamics data with fractures

Molecular dynamics (MD) has served as a powerful tool for designing materials with reduced reliance on laboratory testing. However, the use of MD directly to treat the deformation and failure of materials at the mesoscale is still largely beyond reach. In this work, we propose a learning framework to extract a peridynamics model as a mesoscale continuum surrogate from MD simulated material fracture data sets. Firstly, we develop a novel coarse-graining method, to automatically handle the material fracture and its corresponding discontinuities in the MD displacement data sets. Inspired by the weighted essentially non-oscillatory (WENO) scheme, the key idea lies at an adaptive procedure to automatically choose the locally smoothest stencil, then reconstruct the coarse-grained material displacement field as the piecewise smooth solutions containing discontinuities. Then, based on the coarse-grained MD data, a two-phase optimization-based learning approach is proposed to infer the optimal peridynamics model with damage criterion. In the first phase, we identify the optimal nonlocal kernel function from the data sets without material damage to capture the material stiffness properties. Then, in the second phase, the material damage criterion is learnt as a smoothed step function from the data with fractures. As a result, a peridynamics surrogate is obtained. As a continuum model, our peridynamics surrogate model can be employed in further prediction tasks with different grid resolutions from training, and hence allows for substantial reductions in computational cost compared with MD. We illustrate the efficacy of the proposed approach with several numerical tests for the dynamic crack propagation problem in a single-layer graphene. Our tests show that the proposed data-driven model is robust and generalizable, in the sense that it is capable of modeling the initialization and growth of fractures under discretization and loading settings that are different from the ones used during training.

97 MATHEMATICS AND COMPUTING↗