Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Failure Rate”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Step-loaded creep testing of Zircaloy-4 cladding at higher temperatures in the α-phase

A refined understanding of zirconium-based cladding thermomechanical performance during rapid transients is essential for enhancing the safety and operation of light-water reactors. Traditional models for zirconium alloys under accident conditions generally assume that creep dominates fuel cladding performance. Here, these historic models have largely remained unchanged and serve as the basis for safety criteria development. As the U.S. nuclear industry pursues higher burnup levels, the increased release of fission gases during transients raises the risk of cladding failure in the low-temperature hcp α-phase, making the fidelity of these models of greater importance. Creep testing was conducted from 550–700°C with 25–120 MPa applied hoop stresses to investigate Zircaloy-4 deformation at accident-relevant temperatures in the α-phase. Step-loading was employed to capture creep behavior across a wide stress range from a single sample. The stress-strain rate data at higher temperatures (650 and 700°C) were well-described by isotropic versions of the Erbacher and Kaddour models, while the lower temperature data (550 and 600°C) were underpredicted by both anisotropic and isotropic model variants. Greater strain rates during the initial loading step at 650 and 700°C were attributed to recrystallization and grain growth of sub-micron crystallites. Yet, texture analysis revealed the basal split texture remained after testing. These observations produced results suggesting Zircaloy-4 claddings experience higher creep rates across the α-phase than previously thought, possibly related to dynamic anisotropy due to temperature dependent activation of deformation mechanisms, effects of biaxial loading, and variation in material condition between the current testing used in previous model development.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Switching from batch to continuous granulation: A case study of metoprolol succinate ER tablets

Continuous manufacturing (CM) has been used to produce several immediate release drug products. No extended-release (ER) product manufactured employing CM technology has been approved yet. Herein this study investigated the critical aspects of switching from the batch mode of high shear granulation to the continuous operation of twin-screw granulation for extended-release tablets. Metoprolol succinate ER tablets was used as a model ER formulation for this purpose. A central composite design (CCD) was employed to determine the effects of high shear granulator (HSG) parameters, namely impeller speed, granulation time, and binder liquid feeding rate, on the critical granulation characteristics important for product performance. These critical granulation characteristics served as a guide for switching from the batch processing to the continuous operation for achieving the same breaking strength and dissolution for this ER metoprolol tablets. The granulation time was the most critical factor affecting the bulk properties of granules which contributed to tablet dissolution. The higher density and lower compressibility of granules were attained at the longest granulation time of 5.4 min with the fastest liquid feeding rate of 75 g/min. The granules’ density was the primary factor negatively affecting the dissolution of metoprolol tablets. However, the breaking strength of tablets confounded the effect of granules density on metoprolol dissolution. Switching the processing parameters of high shear granulation to twin-screw granulation achieved similar dissolution profiles (F2 greater than 50). The screw speed was not found to affect bulk properties of granules. The root cause of granulation failures in twin-screw granulation, such as premature consolidation, excessive swelling, poor cohesion, inconsistent shearing effects, and formation of deformed agglomerates, were identified. In conclusion, the use of critical granulation characteristics through a performance-based approach of ER tablets facilitated the switching of manufacturing of an ER formulation form batch to continuous operation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Physical Interpretation of Early Battery Life Prediction Models

Early battery life prediction models are most useful for R&D if they help us understand the early changes in battery electrochemical response that correspond with long-term degradation and failure. Linear regression models such as Fused lasso and Partial Least Squares can fit coefficients directly to high-dimensional electrochemical data like capacity-voltage and ΔV–state-of-charge, i.e., Q(V) and ΔV(SOC) curves, learning coefficients that can be physically interpreted. We leverage the ISU-ILCC battery aging data set to learn high-dimensional coefficients for early battery life prediction from traditional slow-rate capacity check data, demonstrating learning on Q(V), d Q· d V −1 , and ΔV(SOC) curves. A thorough study on the dependence of coefficient values on train/test size and data preprocessing methods is made, demonstrating the reliability of high-dimensional regression approaches unless very small amounts of data are used for model training. For this data set, coefficients from Q(V) and d Q· d V −1 models highlight changes in electrode stoichiometry due to lithium loss, while ΔV(SOC) coefficients highlight changes in positive electrode diffusivity due to particle cracking as well as electrode stoichiometry shifts. By directly interpreting the coefficients of a regression model, we make physical insights into battery degradation mechanisms without requiring the assumptions of traditional battery data analysis methods.

25 ENERGY STORAGE↗
DOE OSTI2026

CDRL: Certification-Driven Reinforcement Learning for Neutrino Flavor Model Discovery

Many scientific discovery problems require searching combinatorial hypothesis spaces under complex domain constraints. Reinforcement learning (RL) offers a promising approach, but existing methods rely on scalar rewards that provide limited information about why candidate solutions fail, leading agents to repeatedly explore invalid regions. We introduce Certification-Driven Reinforcement Learning (CDRL), a framework that leverages structured feedback from symbolic reasoning tools. When a candidate violates domain constraints, these tools produce certificates identifying the actions responsible for failure. CDRL converts these certificates into reusable constraints that eliminate classes of invalid solutions and guide exploration toward valid regions. We evaluate CDRL on neutrino flavor model discovery in theoretical particle physics, where the hypothesis space exceeds $10^{26}$ possible models, and compare it with the state-of-the-art RL approach previously used for this task. Across three theory spaces, CDRL achieves up to 1.95$\times$ higher valid model rates and up to 6.33$\times$ higher neutrino model rates while evaluating up to 4$\times$ fewer candidates. We further extract 40 interpretable rules from search trajectories using a post-hoc decision-tree framework and show that reusing them as soft constraints yields gains of up to 2$\times$ in valid model rates and 3$\times$ in neutrino model discovery across all three theory spaces. These results suggest that CDRL uncovers reusable structure in combinatorial search spaces and provides a general framework for scientific model discovery.

Jha, Piyush [Georgia Tech., Atlanta; Georgia Tech]↗
DOE OSTI2023

Silver-carbon interlayers in anode-free solid-state lithium metal batteries: Current development, interfacial issues, and instability challenges

As an interlayer between the anode and the electrolyte of the all-solid-state lithium metal batteries (ASSLMBs), the silver-carbon (Ag-C) nanocomposite has been reported to significantly increase the energy density and cycle rate of solid-state lithium metal batteries. Ag-C interlayers serve as mixed ionic-electronic conductor that conducts both Li+ ions and electrons and lithium storage capacity. Unfortunately, it was unclear how the Ag-C interlayer regulated lithium plating and stripping. Moreover, the structural and chemical instabilities between the interlayer and the electrolyte, within the interlayer, or beneath the interlayer on lithium substrate are likely to cause cell failure. In this review, we discuss interfacial issues and summarize recent progress in solution strategies for ASSLMBs, with a specific focus on the use of a silver-carbon (Ag-C) nanocomposite interlayer in anode-free setups. Based on the Li transport kinetics among the Ag-C interlayers, the interfacial configurations of Ag-C interlayers are classified as either exterior or internal. Furthermore, the review concludes with a discussion of the perspectives and future prospects, allowing for the improvement of interlayer techniques for solid-state batteries.

25 ENERGY STORAGE↗
DOE OSTI2021

Optimization of the microstructure and mechanical properties of electron beam welded high-strength medium-entropy alloy (NiCoCr) 94 Al 3 Ti 3

For this work, the weldability and effects of subsequent heat treatment on the microstructure and mechanical properties of the medium-entropy alloy (MEA) (NiCoCr) 94 Al 3 Ti 3 have been explored using electron beam welding (EBW). It was found that EBW joints consist of a supersaturated f.c.c. matrix with coarse, asymmetric columnar grains (17.7–24 μm). The fraction of twins and Σ3 n CSL boundaries in weld joints was only 10% of that in the base metal, a value which increased to 30% after aging at 800 °C. The aging also led to improvements in the hardness, tensile strength, and work hardening rate, a result of the precipitation of a high volume fraction (0.29–0.33) of ultrafine (15–17.4 nm dia.) coherent spherical L1 2 -structured particles. The contributions of the L1 2 nanoparticles and Hall-Petch strengthening to the yield strength were estimated to be 537–573 MPa and 105–123 MPa, respectively. The MEA shows excellent weldability and, after aging, good strength (YS~1098 MPa, UTS~1368 MPa) and elongation to failure (~11%) by ductile fracture.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Characterization of a Small-Scale Prototype Detector With Wide Dynamic Range for Time-Resolved High-Energy X-Ray Applications

We present characterization measurements of a fast-framing, wide-dynamic-range x-ray area detector intended for high-energy applications (≥20 keV photons). The MM-PAD-2.1 combines an integrating pixel front-end with a charge-removal mechanism to extend the maximum measurable signal to > 107 20-keV ph/pixel/frame. The charge-removal mechanism is dead-time-less (i.e., incoming signal continues to be integrated by the front-end while charge removal is taking place) up to an incoming photon rate of > 109 20-keV ph/pix/s. A 16×16 pixel prototype was fabricated and characterized using both laboratory and synchrotron x-ray sources, allowing for measurements from single-photon to high-flux conditions. High-flux measurements show linear performance up to ~1010 ph/s in a spot spanning approximately 2×1.5 pixels (FWHM). The onset of temporary radiation damage in the detector was observed during these measurements, manifesting as failure of the digital readout from the detector. These effects were seen to reverse after annealing at 30° C for several hours.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Spatiotemporal Variability in Wind Turbine Blade Leading Edge Erosion

Wind turbine blade leading edge erosion (LEE) reduces energy production and increases wind energy operation and maintenance costs. Degradation of the blade coating and ultimately damage to the underlying blade structure are caused by collisions of falling hydrometeors with rotating blades. The selection of optimal methods to mitigate/reduce LEE are critically dependent on the rates of coating fatigue accumulation at a given location and the time variance in the accumulation of material stresses. However, no such assessment currently exists for the United States of America (USA). To address this research gap, blade coating lifetimes at 883 sites across the USA are generated based on high-frequency (5-min) estimates of material fatigue derived using a mechanistic model and robust meteorological measurements. Results indicate blade coating failure at some sites in as few as 4 years, and that the frequency and intensity of material stresses are both highly episodic and spatially varying. Time series analyses indicate that up to one-third of blade coating lifetime is exhausted in just 360 5-min periods in the Southern Great Plains (SGP). Conversely, sites in the Pacific Northwest (PNW) exhibit the same level of coating lifetime depletion in over three times as many time periods. Thus, it may be more cost-effective to use wind turbine deregulation (erosion-safe mode) for damage reduction and blade lifetime extension in the SGP, while the application of blade leading edge protective measures may be more appropriate in the PNW. Annual total precipitation and mean wind speed are shown to be poor predictors of blade coating lifetime, re-emphasizing the need for detailed modeling studies such as that presented herein.

Pryor, Sara C. (ORCID:0000000348473440)↗
DOE OSTI2023

Laboratory investigation of hydraulic fracturing in granitic rocks using active and passive seismic monitoring

SUMMARY Knowledge of the fracturing processes can be important for the optimization of pressurized fluid injection operations in the deep underground rock mass. Active and passive seismic monitoring techniques have been used in the field for tracking or mapping the propagating hydraulic fracture. Although both these monitoring techniques provide valuable information about the generated fracture network, it is difficult for either technique to comprehensibly identify the different processes associated with hydraulic fracturing. The combined active and passive monitoring has the potential for better characterization of the complex hydraulic fracturing phenomena. In this study, laboratory hydraulic fracturing experiments with combined active and passive seismic monitoring were conducted on true triaxially loaded Barre granite cubes with different fluid injection rates. The seismic inelastic fracturing was detected by 16 passive acoustic emission sensors, where 3678 and 2370 seismic source events were detected for the high and low injection rate experiments, respectively. For active monitoring, strong variations in the attributes of signals were observed which were transmitted through four source–receiver pairs, placed both perpendicular and parallel to the generated hydraulic fracture. Positive velocity changes were observed for active sensor pairs with ray paths passing through the generated hydraulic fracture indicating fluid permeation, whereas isolated dry deformation was characterized by a slight but permanent velocity decrease. Compared to velocity, the energy of the active signals was 1–2 orders of magnitude more sensitive to different hydraulic fracturing processes. However, the sensitivity and signatures of the active signal attributes were found to be dependent on the frequency range and direction of ray path with respect to the location of the generated fracture network. Using the coupled evaluation of the active and passive signals we were able to systematically identify various hydraulic fracturing processes including: (1) aseismic deformation, (2) fracture initiation and fluid permeation, (3) pressure build-up, (4) fracture propagation and (5) pressure release and leak-off. The results of this study showed that combining the respective advantages of active and passive seismic techniques and using both of them to monitor the failure processes can facilitate a more comprehensive understanding and better control of the hydraulic stimulations in subsurface operations.

Geochemistry & Geophysics↗
DOE OSTI2025

Self-Sensing Composites via an Embedded 3D-Printed PVDF-MoS 2 Nanosensor for Structural Health Monitoring

Carbon fiber (CF)-reinforced epoxy composites are widely used in vehicle applications, where early damage detection is crucial for reliability and safety. To address this need, we developed a self-sensing epoxy/CF composite by embedding a PVDF-MoS 2 nanosensor via an embedded 3D printing method. By harnessing the intrinsic curing kinetics of epoxy, we tailored its rheological properties to optimize the embedded printing process, enabling precise and reliable support for sensor filaments without compromising the composite’s structural and functional integrity. Through comprehensive rheological and kinetic analysis, we established a quantitative relationship among curing temperature, conversion rate, and resulting yield modulus─defining a narrow processing window essential for successful sensor integration. Specifically, we identified that an epoxy yield modulus range of 180–294 Pa and a conversion rate below 10% are critical to support the PVDF-MoS 2 filament architecture. Here, this embedded 3D printing method produces complex and multimaterial PVDF-MoS 2 sensors within an epoxy matrix with minimal deformation and reduced postprocessing, which is scalable and adaptable for industrial applications. Under cyclic loading, the embedded sensors exhibited stable signals under constant loads and increased voltage signals in response to crack formation (17–35% higher) and catastrophic failure (1 order of magnitude higher), effectively capturing structural changes in real time. This study demonstrates the potential of PVDF-MoS 2 nanocomposite sensor materials for real-time structural health monitoring in epoxy–CF composite systems, enabling early detection of defects and stress anomalies, significantly reducing the risk of unexpected failures, and enhancing structural reliability.

PVDF-MoS2 sensor↗
DOE OSTI2022

The effect of differential mineral shrinkage on crack formation and network geometry

Rock, concrete, and other engineered materials are often composed of several minerals that change volumetrically in response to variations in the moisture content of the local environment. Such differential shrinkage is caused by varying shrinkage rates between mineral compositions during dehydration. Using both 3D X-ray imaging of geo-architected samples and peridynamic (PD) numerical simulations, we show that the spatial distribution of the clay affects the crack network geometry with distributed clay particles yielding the most complex crack networks and percent damage (99.56%), along with a 60% reduction in material strength. We also demonstrate that crack formation, growth, coalescence, and distribution during dehydration, are controlled by the differential shrinkage rates between a highly shrinkable clay and a homogeneous mortar matrix. Sensitivity tests performed with the PD models show a clay shrinkage parameter of 0.4 yields considerable damage, and reductions in the parameter can result in a significant reduction in fracturing and an increase in material strength. Additionally, isolated clay inclusions induced localized fracturing predominantly due to debonding between the clay and matrix. These insights indicate differential shrinkage is a source of potential failure in natural and engineered barriers used to sequester anthropogenic waste.

58 GEOSCIENCES↗
DOE OSTI2022

Foam Fracturing Study for Stimulation Development of Enhanced Geothermal Systems

The large thermal gradients and high subsurface temperatures of the western region of the U.S. hold great potential for the implementation of enhanced geothermal systems (EGS). The development of these potential EGS resources requires stimulation of the reservoir to enhance permeability and it has been widely reported that a substantial amount of water will be required should conventional hydraulic stimulation be used. This presents a huge challenge and a high risk to the geothermal development because the water stress1 in these areas is already high or extremely high. The use of foam, a gas/liquid mixture predominantly composed of gas, in fracturing is considered and explored in this project as a potential approach to address water concerns with hydraulic stimulation in the development of EGS. This project, led by Oak Ridge National Laboratory (ORNL) in collaboration with Temple University, was awarded in an open lab call in 2018, and was part of the DOE GTO waterless stimulation initiative. The goal of the project was to demonstrate the feasibility of foam fracturing for EGS development through two primary tasks: Task 1: Laboratory study of the effectiveness of foam fracturing for representative geological materials, including cyclic pressurization using foam (led by ORNL) and Task 2: High temperature foam material selection and characterization (led by Temple University). In FY19, ORNL finished the critical review on serval issues associated with foam fracturing and the implementation of the proposed tasks in a lab study (Wang, et al., 2019), and completed the foam fracturing testing using cement as a model material (Wang, et al, 2020a). The work at ORNL was geared up to develop a brand-new foam testing system in FY20. The purchase of main components for the new system was finished in the first half of the FY20. The assembly of the foam testing system and foam fracturing testing were completed in the second half of the FY20 (Wang, et al., 2021a). Task 1 required the development of a test system which can be used to perform hydraulic fracturing of geological specimens with both water and foamed liquids at pressure up to 6,000 psi (41.4 MPa). The system possesses several capabilities that conventional injection systems lack for hydraulic fracturing. In addition to its ability to generate foam with controlled quality, it is capable of cycling pressure levels between specified values up to frequencies of 50 Hz. The latter capability was developed to evaluate the hypothesis that cyclic loading of samples would produce enhanced fracturing. The system consists of two sections: one for foam generation and another for foam injection. The foam is generated through separate control and pressurization of liquid and gas phases with controlled flow rates. The injection section is equipped with a low-flow Coriolis flow that monitors the density of foam to ensure the injection is in the range of target foam quality2. Experimental results of foam fracturing are reported for cylindrical granite specimens using water and aqueous N2 foam as the fracturing fluids. All experiments were performed for unconfined conditions. The effects of injection mode (i.e., monotonic vs cyclic pressurization) on breakdown pressure and failure response sample were investigated using water alone as a fracturing fluid and foams with a range of compositions. It was found that in the case of monotonic injection, the breakdown pressure of granite specimens tended to be slightly higher when fracturing with foam. Additionally, with a foam quality of 90%, the water use can be reduced by 50 to 84%, depending on hole size. On the other hand, it was observed that the breakdown pressure can be brought down to 70% of the monotonic breakdown pressure by using low cycle fatigue. Finally, discussions are presented regarding injectivity and water use reduction.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

DROP DURABILITY ASSESSMENT OF ELECTRONIC ASSEMBLIES UNDER OFF-AXIS LOADING WITH SKEWED FIXTURES

This thesis studies drop durability of electronic assemblies when the acceleration vector is oriented at 45° to the out-of-plane direction of the circuit card. The off-axis drop tests are accomplished with a skewed fixture and are conducted as a proxy for multiaxial drop testing. Advanced shock testing and vibration test methods have been developed over the last few decades to better represent real-world field environments during ground-based laboratory testing. However, many of these test methods require expensive and specialized equipment not available in most laboratories. An alternative approach for approximating simultaneous loading along multiple axes on conventional equipment utilizes skewed fixtures which have seen use in off-axis random vibration and drop impact testing. These methods generally rely on the conversion of a uniaxial input load from the test equipment (using a uniaxial drop tower or shaker) into a multiaxial load when resolved in the reference frame of the test article (mounted on a skewed fixture). Skewed fixture design is presented and recommendations for conducting skewed angle drop testing are introduced based on local measurements along the skewed face of the fixture to accurately monitor the impact event. Characterization tests were performed with a skewed fixture, at simultaneous acceleration loads from 500 to 3,000 g in two (in-plane and out-of-plane) directions, while meeting standard time domain tolerances. Upon experimental characterization, drop shock durability tests were conducted on a printed circuit assembly (PCA). Mean drops-to-failure were measured and quantified with Weibull statistics. Dominant solder joint failure modes were identified via failure analysis. Prior work on inclined angle impact testing is limited, and the majority of solder joint interconnect level fatigue studies are conducted considering perpendicular loading normal the circuit card. Low-cycle fatigue curves are generated based on plastic strain and plastic work density within the solder joint. A multiscale nonlinear finite element model is used to relate board-level flexure to solder joint interconnect level plastic strain. A high strain rate solder constitutive model allows for accurate modeling of solder plasticity resulting from high-impact drop shock. Fatigue parameters are computed from the Coffin-Manson relation and Palmgren-Miner damage accumulation. This work serves to apply established low-cycle fatigue methods for conventional drop shock loading (impact normal to circuit card) to non-perpendicular loading with a skewed fixture.

Hower, Jonathan [Kansas City National Security Cam↗
DOE OSTI2025

The Ductility of 49Fe-49Co-2V Soft Magnetic Alloy Bar: Surface Effects and Test Methods

The tensile ductility of 49Fe-49Co-2V (Hiperco® 50A) bar was investigated in both as-received and heat-treated conditions. The as-received/machined specimens exhibit very low ductility compared to samples where heat treatment was the final step prior to testing. Microstructural characterization showed that internal residual strain from bar processing and, most importantly, surface machining damage, cause lower elongation in the as-received material. Because fracture of this intermetallic alloy initiates at the surface, it is particularly susceptible to surface machining damage, i.e., the near-surface region has already exhausted most of its ability to accumulate tensile strain. During heat treatment, the internal residual strain and near-surface machining damage are eliminated and ductility is improved, despite a higher degree of crystallographic ordering in the heat-treated condition (which typically lowers ductility). Furthermore, if machining is again performed after heat treatment, the material again exhibits brittle behavior, even with only light touch-up machining passes. Here, in this work, methods of tensile strain measurement were investigated, namely conventional knife-edge extensometry and noncontact digital image correlation (DIC) on heat-treated material. For clip-on knife-edge extensometry, the range of failure strain was 2.5-5.5% for heat-treated Hiperco. For noncontact methods, ductility up to 7% was observed. The results highlight the tendency for the alloy to fail at surface imperfections, even those produced by application of the extensometer itself. Noncontact laser extensometry is recommended for determining the intrinsic ductility of the alloy. A method of laser surface modification was developed which increased ductility by ~ 100% compared to unmodified samples. The high cooling rates achieved during laser surface processing can bypass the ordering reaction and produce a ductile disordered structure at the surface that exhibits ductile fracture characteristics.

EBSD↗
DOE OSTI2022

Microstructure-sensitive modeling of high temperature creep in grade-91 alloy

Predicting the effects of microstructure on high-temperature creep responses of steel components is critical to minimize the risks of failure and maximize economic viability in the energy sector. In this research, a recently developed advanced mechanistic constitutive model is employed to study the effect of microstructure on the creep responses and to rationalize experimentally reported variability in the performances of grade-91 alloy. Hundreds of experimental creep tests from the literature are used to assess the predictability of the model. The model proposed is shown to accurately predict the steady-state creep rates and also the less commonly considered, yet important, primary-to-secondary transients for a wide range of creep conditions. Further, the model is capable of extrapolating the creep response of grade 91 alloy even outside of the calibration regime. Using this model, the roles of initial microstructure described in terms of grain size, dislocation density, and precipitate content, on the creep responses are investigated. The effect of initial microstructure on the steady-state creep rates and creep-strain is found to vary with stress and temperature. Grain size plays a significant role in the low-stress creep, where the diffusional creep is dominant. On the other hand, the effects of dislocation densities and precipitate content are significant in the dislocation-plasticity-dominated high-stress and temperature regimes. Furthermore, by comparing model predictions against a large experimental creep database, it is found that variability in the creep responses can be rationalized on the basis of differences in the initial microstructure. Overall, this work provides a robust pathway for microstructure-aware engineering design, verification, and validation of metallic components for energy applications.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A fourth-order phase-field fracture model: Formulation and numerical solution using a continuous/discontinuous Galerkin method

Modeling crack initiation and propagation in brittle materials is of great importance to be able to predict sudden loss of load-carrying capacity and prevent catastrophic failure under severe dynamic loading conditions. Second-order phase-field fracture models have gained wide adoption given their ability to capture the formation of complex fracture patterns, e.g. via crack merging and branching, and their suitability for implementation within the context of the conventional finite element method. Higher-order phase-field models have also been proposed to increase the regularity of the exact solution and thus increase the spatial convergence rate of its numerical approximation. However, they require special numerical techniques to enforce the necessary continuity of the phase field solution. In this paper, we derive a fourth-order phase-field model of fracture in two independent ways; namely, from Hamilton’s principle and from a higher-order micromechanics-based approach. The latter approach is novel, and provides a physical interpretation of the higher-order terms in the model. In addition, we propose a continuous/discontinuous Galerkin (C/DG) method for use in computing the approximate phase-field solution. This method employs Lagrange polynomial shape functions to guarantee -continuity of the solution at inter-element boundaries, and enforces the required regularity with the aid of additional variational and interior penalty terms in the weak form. Finally, the phase-field equation is coupled with the momentum balance equation to model dynamic fracture problems in hyper-elastic materials. Two benchmark problems are presented to compare the numerical behavior of the C/DG method with mixed finite element methods.

42 ENGINEERING↗
DOE OSTI2022

Asynchronous Reciprocal Coupling of Martini 2.2 Coarse-Grained and CHARMM36 All-Atom Simulations in an Automated Multiscale Framework

The appeal of multiscale modeling approaches is predicated on the promise of combinatorial synergy. However, this promise can only be realized when distinct scales are combined with reciprocal consistency. Here, we consider multiscale molecular dynamics (MD) simulations that combine the accuracy and macromolecular flexibility accessible to fixed-charge all-atom (AA) representations with the sampling speed accessible to reductive, coarse-grained (CG) representations. AA-to-CG conversions are relatively straightforward because deterministic routines with unique outcomes are achievable. Conversely, CG-to-AA conversions have many solutions due to a surge in the number of degrees of freedom. While automated tools for biomolecular CG-to-AA transformation exist, we find that one popular option, called Backward, is prone to stochastic failure and the AA models that it does generate frequently have compromised protein structure and incorrect stereochemistry. Although these shortcomings can likely be circumvented by human intervention in isolated instances, automated multiscale coupling requires reliable and robust scale conversion. In this work, we detail an extension to Multiscale Machine-learned Modeling Infrastructure (MuMMI), including an improved CG-to-AA conversion tool called sinceCG. This tool is reliable (~98% weakly correlated repeat success rate), automatable (no unrecoverable hangs), and yields AA models that generally preserve protein secondary structure and maintain correct stereochemistry. We describe how the MuMMI framework identifies CG system configurations of interest, converts them to AA representations, and simulates them at the AA scale while on-the-fly analyses provide feedback to update CG parameters. Application to systems containing the peripheral membrane protein RAS and proximal components of RAF kinase on complex eight-component lipid bilayers with ~1.5 million atoms is discussed in the context of MuMMI.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Algorithm 1049: The Delaunay Density Diagnostic

Accurate approximation of a real-valued function depends on two aspects of the available data: the density of inputs within the domain of interest and the variation of the outputs over that domain. There are few methods for assessing whether the density of inputs is sufficient to identify the relevant variations in outputs—i.e., the “geometric scale” of the function—despite the fact that sampling density is closely tied to the success or failure of an approximation method. In this article, we introduce a general purpose, computational approach to detecting the geometric scale of real-valued functions over a fixed domain using a deterministic interpolation technique from computational geometry. The algorithm is intended to work on scalar data in moderate dimensions (2–10). Our algorithm is based on the observation that a sequence of piecewise linear interpolants will converge to a continuous function at a quadratic rate (in L 2 norm) if and only if the data are sampled densely enough to distinguish the feature from noise (assuming sufficiently regular sampling). We present numerical experiments demonstrating how our method can identify feature scale, estimate uncertainty in feature scale, and assess the sampling density for fixed (i.e., static) datasets of input–output pairs. Finally, we include analytical results in support of our numerical findings and have released lightweight code that can be adapted for use in a variety of data science settings.

97 MATHEMATICS AND COMPUTING↗