Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “two dimensional materials”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Effect of Thermal Aging on Microstructure and Stress Corrosion Cracking Behavior of an Alloy 152 1st Layer Butter Weldment

Nickel-based Alloy 690 and the associated weld Alloys 52 and 152 are typically used for nozzle penetrations in replacement heads for pressurized water reactor (PWR) vessels, because of their excellent overall resistance to general corrosion and environmental degradation, primarily stress corrosion cracking (SCC). However, many of the existing PWRs are expected to operate for 40- 80 years. Likewise, water-cooled small modular reactors (SMRs) will use Ni-Cr alloys and are expected to receive initial operating licenses for 60 years. Hence, the thermal stability of Ni-Cr alloys is critical for the long-term performance of both existing and advanced nuclear power plants, and possibly spent fuel storage containers. The objective of this research is to understand the microstructural changes occurring in high-Cr, Ni-based Alloy 152 weldments during long time exposure to the reactor operating temperatures, and the effect of these changes on the service performance. One area of particular concern is the potential for long range ordering (LRO), i.e. formation of the intermetallic Ni2Cr phase under prolonged exposure to reactor temperatures and/or irradiation, which can increase strength, decrease ductility, and cause dimensional changes or lead to in-service embrittlement of components made with these alloys. Hence, this research focused on the microstructural evolution and the SCC response of Alloy 152 under accelerated thermal aging. The materials studied involved three heats of Alloy 152 used to produce a dissimilar metal weld (DMW) joining an Alloy 690 plate to an Alloy 533 low alloy steel (LAS) plate, thermally aged at three different temperatures (370°C, 400°C and 450°C) for different durations up to 75,000h (equivalent to 60 years of reactor service). The microstructural characterization by means of synchrotron X-ray conducted in small, 0.2 mm - step line scans in the high-deformation regions of the weld root – covering areas spanning from the weld heat affected zone (HAZ) in Alloy 690 to the weld and weld butter on LAS - did not show evidence of LRO in any of the three Alloy 152 heats aged at 370°C and 450°C to an equivalent of 60 years of service. However, the first weld butter layer has high levels of deformation and is highly susceptible to SCC even in its non-aged condition. Nanohardness testing confirmed the extreme hardening with aging (ΔHV ≅ 100) at two locations within this weldment. In absence of LRO, hardening is suspected to be due to thermally-induced Cr carbide precipitation and coarsening. Testing in a primary water environment of the 1st layer of Alloy 152 weld butter aged at 370°C to a 60-year service equivalent revealed a fatigue and corrosion fatigue crack growth responses similar to those measured on the un-aged alloys. Similarly, the SCC CGR response of the aged weld butter does not appear to show a deterioration in performance, however, the difficult-to-test weldment geometry may affect the test results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Angle-Resolved Polarized Raman Study of Layered Cr 2 Se 3

The polarization-resolved Raman spectra of two-dimensional Cr 2 Se 3 synthesized via chemical vapor deposition (CVD) and chemical vapor transport (CVT) techniques were investigated in detail. The samples were characterized using X-ray diffraction (XRD), transmission electron microscopy (TEM), and energy-dispersive X-ray spectroscopy (EDS). A distinct polarization dependence was observed in the Raman intensity of all the Cr-Cr, Cr-Se, and Se-Se modes in both samples. The observed angle-dependent Raman intensities of each peak could be related to the crystal structure-specific Raman tensor. XRD results of the bulk Cr 2 Se 3 sample synthesized via CVT confirm its trigonal crystal structure, and the Raman peaks can be fitted using the Raman tensors for the A g and E g modes for both the parallel and crossed polarizations. However, for the Cr 2 Se 3 samples directly grown on Si/SiO 2 substrates by CVD, it was necessary to assume the triclinic crystal structure in order to explain the polarized Raman dependence of all the peaks in both parallel and crossed polarization directions. Furthermore, this is the first experimental result suggesting the existence of triclinic Cr 2 Se 3 crystal structure, which has been theoretically predicted in the Materials Project database.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Control of core–shell nanoparticles properties through plasma synthesis: a computational study

The improved properties of core–shell nanoparticles (CSNPs) over homogeneous nanoparticles (NPs) have expanded and diversified the applications of these nanomaterials. However, controlling the properties of CSNPs can be a challenging task. Low temperature plasmas have proven to be an effective method of producing NPs with uniform size and morphology, and high yield. That said, NP transport and growth dynamics are sensitive to LTP properties. We report on a computational investigation of the evolution of Ge–Si CSNP properties as a function of operating conditions through the modeling of a flowing, two-zone inductively coupled plasma (ICP) reactor. Ar/GeH 4 and Ar/SiH 4 gas mixtures were supplied to separate plasma zones at a pressure of 1 Torr to promote growth of Ge cores and Si shells. The negatively charged CSNPs are trapped electrostatically in the vicinity of the antennas where the plasma is generated and where the majority of particle growth occurs. Particles that grow to a critical size are then de-trapped by fluid drag due to neutral gas flow. A two-dimensional hybrid plasma model coupled with a three-dimensional kinetic NP transport model were utilized to resolve plasma chemistry and NP growth processes that take place on distinct timescales. The trends in CSNP properties and trapping mechanisms associated with flow rate, applied ICP power and inlet precursor fraction are discussed. While the spatial distribution of plasma produced radical species can have significant impact on the NP growth process, the NP transport dynamics are what ultimately dictates the growth environment that is unique to each particle and so determines their final dimension and composition. The key to optimizing reactor conditions involves controlling the spatial density of growth species and plasma profile as a means to tailor particle trapping dynamics suitable to produce CSNPs for a specific application.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Development, Verification, and Validation of an OpenFOAM-Based Solver for Modeling Inertial Fusion Energy Chambers

Our work seeks to introduce a computational tool tailored to the physics of inertial fusion energy chambers, in particular, those concepts based on thick liquid walls. In this approach, the structural materials are protected by several neutron mean-free-paths of renewable liquid and thus will be able to survive much longer than un-shielded walls, with virtually all structures lasting for the life of the plant and enabling the use of commercially available and qualified materials. The OpenFOAM-based solver named rhoCentralFoam has been used as a starting point. rhoCentralFoam belongs to the standard OpenFOAM solver toolset. It is a high-speed, explicit compressible flow solver with shock-capturing capability. While the main features have been retained, the solver had to be restructured to make use of tabular data for equations of states, a necessary addition to model the complex thermo-physical properties of ionized gasses. This entailed the need to change the independent state variables used by the solver, resulting in a new thermodynamic library and slightly different solution algorithm. Moreover, a radiation heat transfer model based on the P-1 approximation was added to the solver. The solver is verified against an analytical solution from the Sedov-Taylor-Neumann test problem to showcase the ability of the hydrodynamic solvers to handle strong shocks, whereas the P-1 model was verified using a simple one-dimensional problem with an analytical solution. Additionally, a validation case involving shock-wave propagation through jet array is presented, and the results are compared with experimental data from the open literature. Lastly, in order to showcase the utility of the solver for practical cases, we applied the refined solver to two representative scenarios: gas venting within the HYLIFE-II chamber and the compression of the gas following the partial ablation of the liquid wall.

Chamber dynamics↗
DOE OSTI2025

Bayesian inference of anisotropic 2D small-angle scattering from sparse measurement

Here, we present a Bayesian inference framework for reconstructing anisotropic two-dimensional small-angle scattering (2D SAS) patterns from sparse, noisy, or partially missing data. The method combines a symmetry-aware angular basis with radial Gaussian process priors to enable accurate, training-free interpolation and denoising. Computational benchmarks demonstrate reliable recovery of both isotropic and high-order anisotropic features under severe data reduction. Experimental validations on stretched polymers, sheared wormlike micelles, and carbon fibers show improved fidelity and resolution compared to raw measurements, achieving comparable accuracy with up to 50-fold fewer detected neutrons. This approach enables quantitative structural analysis under low-flux, time-limited, or single-shot conditions, extending the applicability of 2D SAS techniques to compact neutron sources and mechanically driven soft matter systems undergoing transient structural changes.

Tung, Chi-Huan [Oak Ridge National Laboratory (ORN↗
DOE OSTI2025

Niobium coaxial cavities with internal quality factors exceeding 1.4 × 10 9 for circuit quantum electrodynamics

Group-V materials such as niobium and tantalum have become popular choices for extending the performance of circuit quantum electrodynamics (cQED) platforms, allowing for quantum processors and memories with reduced error rates and more modes. The complex surface chemistry of niobium, however, makes identifying the main modes of decoherence difficult at millikelvin temperatures and single-photon powers. We use niobium coaxial quarter-wave cavities to study the impact of etch chemistry, prolonged atmospheric exposure, and the significance of cavity conditions prior to and during cooldown—in particular, niobium hydride evolution—on single-photon coherence. We demonstrate cavities with quality factors Q int ≳ 1.4 × 10 9 in the single-photon regime, a 15-fold improvement over aluminum cavities of the same geometry. We rigorously quantify the sensitivity of our fabrication process to various loss mechanisms and demonstrate a two- to fourfold reduction in the two-level system loss tangent and a three- to fivefold improvement in the residual resistance over traditional buffered chemical polishing etching techniques. Finally, we demonstrate transmon integration and coherent cavity control while maintaining a cavity coherence of 11.3 ms. The accessibility of our method, which can be easily replicated in academic laboratory settings, together with the demonstration of its performance, mark an advancement in three-dimensional cQED.

Oriani, Andrew E. [Chicago U., EFI; Chicago U.] (O↗
DOE OSTI2024

Development of a Griffin model of the advanced test reactor

In the pursuit of a higher fidelity deterministic simulation capability of the Advanced Test Reactor, it is important to have a fast yet accurate deterministic neutronics model. Here, to achieve this, we employed an advanced two-step method. The first step involves generating homogenized cross sections using OpenMC, a cutting-edge Monte Carlo neutron transport code. OpenMC offers excellent modular capabilities, allowing for easy component integration and flexibility in incorporating new designs into the model. The second step involves deterministic transport calculations, which are performed using Griffin, a reactor physics application based on the Multiphysics Object-Oriented Simulation Environment (MOOSE). To ensure the accurate spatial resolution and assignment of material cross sections, a Cubit-generated mesh for the Advanced Test Reactor is utilized as an intermediate step between the OpenMC and Griffin models; Griffin utilizes the mesh for its finite element solution, while OpenMC material identifications are written to the mesh file to be used in Griffin material assignments. Additionally, a Python-based script converts the cross sections generated by OpenMC into the ISOXML format required by Griffin. Initial comparisons using the Griffin diffusion solver indicated good agreement between the neutron multiplication factors obtained from the standalone OpenMC model and the Griffin model, with differences of less than 10 pcm in the 2D geometry configuration; it was later determined that this agreement was likely due to compensating effect and was more likely on the order of –700 pcm relative to the OpenMC solution. However, in three-dimensional calculations, an unacceptably large error (almost 8,000 pcm) was found in the Griffin solution with the diffusion solver. Subsequent calculations using Griffin’s discrete ordinates solver demonstrated substantially improved agreement, within 116 pcm of the OpenMC solution used to generate the cross sections for Griffin. Building on this capability, future work will seek to perform more detailed validation calculations. The ultimate goal is to evaluate both transient and multiphysics simulations of the reactor.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2026

Pixelated plastic scintillator array manufacturing using fast-, photo-curable resin

Pixelated plastic scintillator arrays can serve as high efficiency and high resolution neutron imaging detectors. Manufacturing these arrays is intensive in both time and labor. This article presents a fabrication method based on additive manufacturing for two-dimensional plastic organic scintillator arrays using a custom-built automated assembly machine and a custom photocurable resin that has significant non-aromatic acrylate oligomer content. The process involves two main stages: fully autonomous production of one-dimensional layered arrays, followed by semi-autonomous cutting and stacking to form two-dimensional pixel arrays. One-dimensional arrays were manufactured at a rate of around 4 layers per hour with minimal defects and tight dimensional tolerances, while two-dimensional arrays up to 7 x 7 pixels and 70 mm in length were completed in approximately 3.5 hours. Final arrays exhibited dimensional deviations of less than 0.5 mm. Two-dimensional arrays read out by a multi-anode photomultiplier tube demonstrated per-pixel position resolution and pulse-shape discrimination, enabling gamma–neutron interaction separation in mixed radiation environments.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Enhanced Interlayer Coupling and Excitons in Twin-Stacked Two-Dimensional Magnetic CrSBr Bilayers

The degree of electronic coupling between individual layers in van der Waals heterostructures offers a route to engineer their magnetic, electronic, and optical functionalities. Using state-of-the-art first-principles calculations, in this work, we demonstrate that the electronic coupling between two monolayers of CrSBr─an anisotropic two-dimensional magnetic semiconductor─is highly nonlinear and nonmonotonic with respect to their relative twist angle, exhibiting a pronounced maximum at the twin-stacking configuration. The coupling strength scales with both the degree of overlap of Br orbitals adjacent to the van der Waals gap and the cosine of half of the interlayer spin angle. This enhanced interlayer electronic coupling leads to excitons delocalized across the two layers, with a polarization dependence that reflects the interlayer spin alignment. Our results reveal a sensitive interplay among twist angle, magnetism, and excitonic properties in twin-stacked CrSBr bilayers, suggesting twin stacking as an effective means for engineering interlayer coupling.

2D materials↗
DOE OSTI2025

Strain-induced topological transitions in epitaxial films of cadmium arsenide

Strain in epitaxial thin films can have a profound influence on their electronic bands. Here, we investigate the influence of epitaxial strains on the electronic states of (001)-oriented cadmium arsenide (Cd 3 As 2 ) thin films, which are two-dimensional topological insulators. Cd 3 As 2 films were grown coherently on Al 1−x In x Sb buffer layers with varying degrees of lattice mismatch to Cd 3 As 2 , producing in-plane stresses that ranged from compressive to tensile. Using magneto-transport measurements of the Landau levels of the lowest subbands, we show that films under compression remain in the two-dimensional topological insulator state, while a sufficiently large tensile stress induces a transition to a topologically trivial phase. These results qualitatively agree with expectations from theoretical models for these films. At the topological transition, the gap avoids completely closing, and we discuss possible origins.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Separate Surface and Bulk Topological Anderson Localization Transitions in Disordered Axion Insulators

In topological phases of matter for which the bulk and boundary support distinct electronic gaps, there exists the possibility of decoupled mobility gaps in the presence of disorder. This is in analogy with the well-studied problem of realizing separate or concomitant bulk-boundary criticality in conventional Landau theory. Using a three-dimensional axion insulator having clean, gapped surfaces with 𝑒 2 /2⁢ℎ quantized Hall conductance, we show that the bulk and surface mobility gap evolve differently in the presence of disorder. The decoupling of the bulk and surface topology yields a regime that realizes a two-dimensional, unquantized anomalous Hall metal in the Gaussian unitary ensemble on each surface, which shares some spectral and response properties akin to the surface states of a conventional 3D topological insulator. The generality of these results, as well as extensions to other insulators and superconductors, is discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Spin-dependent electron transfer in electrochemically transparent van der Waals heterostructures for oxygen evolution reaction

Spin selective catalysis is an emerging approach for improving the thermodynamics and kinetics of reactions. The role of electron spins has been scarcely studied in catalytic reactions. One exception is the oxygen evolution reaction (OER) where strongly correlated metals and oxides are used as catalysts. In OER, spin alignment facilitates the transition of singlet state of the reactant to the triplet state of O 2 . However, the influence of strong correlations on spin exchange mechanism and spin selective thermodynamics of most catalytic reactions remain unclear. Here we decouple the strongly correlated catalyst from the electrolyte to study spin exchange in two-dimensional (2D) magnetic iron germanium telluride (FGT) heterostructure. We demonstrate that transmission of spin and electrochemical information between the catalyst and the reactant can occur through quantum exchange interaction despite the catalyst of FGT being completely encapsulated by graphene or hexagonal boron nitride (hBN). The strong correlations in FGT that lead to enhanced spin exchange in OER are observed in graphene or hBN layers with thicknesses of up to 6 nm. We demonstrate that spin alignment in FGT leads to a lowering of thermodynamic barrier for adsorption of hydroxide ion and electron transfer to the catalyst. This results in up to fivefold enhancement in OER performance and improved kinetics. Our results provide clear evidence that transmission of both quantum mechanical and electrochemical information through quantum spin exchange interaction in FGT leads to an enhancement in catalytic performance.

2D materials↗
DOE OSTI2025

Synthetic Pathways to gamma-Graphyne and Related Allotropes of Carbon

Graphynes, two-dimensional carbon lattices combining sp 1 and sp 2 hybridized atoms, were predicted theoretically more than three decades ago, but few structures have been realized to date. These carbons are believed to possess remarkable mechanical and electronic properties, including high charge carrier mobilities comparable to those in graphene (10 4 to 10 5 cm 2 V -1 s -1 ). Unlike graphene, certain graphynes are predicted to be intrinsic semiconductors. Among these intriguing structures, γ-graphyne stands out as the structurally simplest and most symmetric sp 1 /sp 2 lattice. γ-Graphyne was first theorized in 1987. In contrast with graphene, γ-graphyne will be a semiconductor with a small band gap suitable for fabrication of electronic devices. This solves one of the fundamental problems of carbon-based electronics, the necessity for inducing a band gap in graphene. γ-Graphyne has the potential to form the basis for the next generation of carbon-based electronics operating at speeds unattainable by traditional silicon chips. Unlike silicon, γ-graphyne is a direct band gap semiconductor, and it will feature exceptional strength comparable to that of other 2D carbon allotropes. Such combination of properties may enable a new generation of highly efficient, ultra-light and flexible solar cells. Despite being a potentially “magical” material, γ-graphyne remained synthetically elusive for over three decades. The primary goals of this project were: (1) Synthesis of bulk γ-graphyne phases through solution-phase 2D polymerizations; (2) Experimental exploration of the physical and chemical properties of γ-graphyne; and (3) Mechanistic and theoretical studies of the novel chemical transformations developed in Goal 1. Common pyrolytic and vapor-deposition methodologies used for the synthesis of graphitic allotropes are unsuitable for graphyne and other sp 1 -contaning structures, as acetylenes readily convert to graphene and amorphous carbon at high temperatures. In contrast, this proposal is based on solution-based 2D polymerization. The major advantages of this approach over the traditional high temperature techniques are the potential to adjust the structure of the material with atomic precision, and the possibility of using structurally complex and relatively fragile repeat units. The outcomes of this research can revolutionize carbon nanotechnology, expanding the field’s structural toolbox beyond primarily graphitic and benzenoid structures. Understanding the chemistry of sp 1 carbon allotropes can lead to entirely new classes of structures with unique properties, including graphyne ribbons, nanotubes, quantum dots, and heterostructures with other 2D materials. Furthermore, the development of reliable and robust synthetic pathways towards periodic covalent molecular sheets with atomically precise structures shall have a profound impact on chemistry and materials science.

2D polymerization↗
DOE OSTI2025

Nanoscale control of LaAlO 3 /SrTiO 3 conductance beyond the polar catastrophe

Conductive atomic force microscope (c-AFM) lithography can be utilized to create a wide range of LaAlO 3 /SrTiO 3 (LAO/STO)-based nanoelectronic devices in a reconfigurable manner. Experiments were generally performed with intrinsically insulating LAO/STO heterostructures, with LAO thickness less than the critical value at which a polar catastrophe takes place [<4 unit cell (u.c.)]. Here, we use inductively coupled plasma reactive ion etching (ICPRIE) to fabricate c-AFM “canvases” on intrinsically conducting LAO/STO samples with ≥4 u.c. LAO. We observe that its interfacial two-dimensional electron gas (2DEG) can be pinched off and then switched back on by c-AFM lithography. Nanowires created with initially conductive LAO/STO interfaces have an order-of-magnitude longer lifetime in ambient conditions, when compared to an identically created 3.4 u.c. LAO/STO nanowire. We also demonstrate key nanoscale properties such as ballistic transport in a quasi-one-dimensional electron waveguide at a 5 u.c. LAO/STO interface. This approach frees c-AFM-written nanodevice designs from time constraints in air associated with <4 u.c. LAO/STO heterostructures.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Intern Poster

Digital Image Correlation (DIC) is an optical technique that combines image registration and tracking methods for accurate two-dimensional and three-dimensional changes in images. DIC software can be used to track the contour, deformation, and strain of a sample. In the Advanced Test Reactor (ATR) at INL (Idaho National Laboratory) there exists a small working window of samples that can become irradiated. Hundreds of graphite disks called piggybacks have undergone this irradiation as part of the Advanced Reactor Technologies (ART) program. After irradiation, it is desirable to understand the change in tensile strength as a function of dose. Due to the limited space in the ATR, typical dog bone tensile tests reduce the number of graphite samples from hundreds to tens. However, there does exist an ASTM standard, D8289, which uses disc compression of graphite to estimate the tensile strength of the specimen with the Brazilian Disk test fixture. While only used as an estimate, which is typically off by a third, it is the purpose of this study to identify how to amend D8289 to remove the word "estimate" with the help of DIC.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Magnetic Phase Diagram of 2D Magnet Fe 5– x GeTe 2 Revealed by Cryogenic Lorentz 4D-STEM

A two-dimensional van der Waals ferromagnet Fe 5–x GeTe 2 (F5GT) exhibits room-temperature magnetic transition and magnetic anisotropy. Studies have identified diverse magnetic states, including stripe domains, skyrmionic type-I and topologically trivial type-II bubbles, depending on the magnetic and thermal history. Yet, the underlying micromagnetic energetics driving these states remains unclear. Here, we establish the magnetic phase diagram of bulk F5GT using cryogenic Lorentz four-dimensional scanning transmission electron microscopy. We reveal that stripe domains spontaneously form upon zero-field cooling, while metastable magnetic bubbles are stabilized by external fields. Transitions from type-I to type-II bubbles are driven by an oblique external field. Micromagnetic simulations confirm that these transitions arise from the interplay of uniaxial anisotropy, dipolar interactions, and external fields, without requiring significant Dzyaloshinskii–Moriya Interaction (DMI). In conclusion, these findings clarify the micromagnetic origin of spin textures of F5GT and establish it as a tunable platform for nanoscale topological magnetism.

2D Magnets↗
DOE OSTI2026

Atomic-Scale Imaging Reveals Polar-π Interactions in Two-Dimensional Molecular Superlattices

Controlling coassembly of synthetic oligomers into binary superlattices at the atomic level is challenging. Here, we report a strategy for programming polar-π interactions in oligomeric peptoids, a class of sequence-defined peptidomimetics, facilitating the formation of homogeneous two-dimensional (2D) superlattices. N-2-phenylethyl and N-(2-perfluorophenyl)ethyl side chains, similar in size, but with contrasting electrostatic characteristics, were introduced at defined sequence positions to generate favorable dipolar aromatic interactions. The resulting nanosheets exhibit different crystal motifs depending on the side chain interactions: systems containing only one type of aromatic side chain form a parallel V-shaped motif driven by π-π interactions, whereas a combination of both types of aromatic side chains, either within one backbone or through the coassembly of two distinct peptoids, adopt an antiparallel V-shaped superlattice with higher thermal stability, driven by polar-π interactions. Cryogenic transmission electron microscopy directly resolved the packing arrangement of perfluorophenyl and phenyl rings in individual nanosheet superlattices, confirming that intermolecular polar-π interaction dominates the superlattice motifs and increases lattice stability. Molecular dynamics simulations and density functional theory calculations further substantiate the energetic favorability of polar-π interactions over π-π interactions, rationalizing the formation of homogeneous superlattices with enhanced thermal stability. Our discoveries establish a design principle for binary coassembly using sequence-defined oligomers, which enables control over unit cell geometry, lattice stability, and molecular registration through aromatic side chain polarization and sequence control. This ability to program atomic-scale binary superlattices opens new avenues for designing functional 2D soft materials.

Lee, Yen Jea [Lawrence Berkeley National Laborator↗
DOE OSTI2026

A Data-Driven Method for Modeling Creep-Fatigue Stress- Strain Behavior Using Neural ODEs

In this paper, we introduce a data-driven machine learning approach for modeling one-dimensional stress–strain behavior under cyclic loading, utilizing experimental data from the nickel-based Alloy 617. The study employs uniaxial creep–fatigue test data acquired under various loading histories and compares two distinct neural network-based ODE models. The first model, known as the black-box model, comprehensively describes the strain–stress relationship using a Neural ODE equation. To interpret this black-box model, we apply the Sparse Identification of Nonlinear Dynamical Systems (SINDy) technique, transforming the black-box model into an equation-based model using symbolic regression. The second model, the Neural flow rule model, incorporates Hooke’s Law for the linear elastic component, with the nonlinear part characterized by a Neural ODE. Both models are trained with experimental data to accurately reflect the observed stress–strain behavior. We conduct a detailed comparison with the standard Chaboche model, which includes three back stresses. Our results demonstrate that the neural network-based ODE models precisely capture the experimental creep–fatigue mechanical behavior, exceeding the standard Chaboche model’s accuracy. Furthermore, an interpretable model derived from the black-box neural ODE model through symbolic regression achieves accuracy comparable to the Chaboche model, enhancing its interpretability. The results highlight the potential of neural network-based ODE models to depict complex creep–fatigue behavior, eliminating the necessity for experts to define a specific, material-focused model form.

creep-fatigue↗