skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on September 17, 2026

Title: Learning the governing partial differential equation of solid diffusion from atomistic simulations
Governing partial differential equations (PDEs) play a critical role in materials research and applications, as they describe essential physics underlying materials behaviour. Traditionally, these equations are developed through phenomenological modelling of experimental results or first principle analysis based on conservation laws. In addition, molecular dynamics (MD) simulations capture atomistic-scale behaviour with detailed physics. However, translating atomistic insights into continuum-scale governing equations remains a significant challenge. Empowered by recent advances in data-driven modelling, we develop a computational framework to learn governing PDEs directly from atomistic simulation data. The framework integrates numerical differentiation of MD data with the identification of constitutive relationships. It proves effective and efficient in learning governing PDEs from noisy and limited MD datasets, without requiring prior knowledge of the final PDEs. Using this framework, we identify a nonlinear PDE governing solid-state diffusion in nickel–hydrogen alloys. This PDE reveals a highly concentration-dependent diffusivity that varies over an order of magnitude. Our data-driven computational framework paves the way for cross-scale constitutive modelling.  more » « less
Award ID(s):
2138431
PAR ID:
10636345
Author(s) / Creator(s):
;
Publisher / Repository:
The Royal Soceity
Date Published:
Journal Name:
Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences
Volume:
481
Issue:
2322
ISSN:
1364-5021
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Traditional data-driven deep learning models often struggle with high training costs, error accumulation, and poor generalizability in complex physical processes. Physics-informed deep learning (PiDL) addresses these challenges by incorporating physical principles into the model. Most PiDL approaches regularize training by embedding governing equations into the loss function, yet this depends heavily on extensive hyperparameter tuning to weigh each loss term. To this end, we propose to leverage physics prior knowledge by “baking” the discretized governing equations into the neural network architecture via the connection between the partial differential equations (PDE) operators and network structures, resulting in a PDE-preserved neural network (PPNN). This method, embedding discretized PDEs through convolutional residual networks in a multi-resolution setting, largely improves the generalizability and long-term prediction accuracy, outperforming conventional black-box models. The effectiveness and merit of the proposed methods have been demonstrated across various spatiotemporal dynamical systems governed by spatiotemporal PDEs, including reaction-diffusion, Burgers’, and Navier-Stokes equations. 
    more » « less
  2. Due to the complex behaviour of amorphous shape memory polymers (SMPs), traditional constitutive models often struggle with material-specific limitations, challenging curve-fitting, history-dependent stress calculations and error accumulation from stepwise calculation for governing equations. In this study, we propose a physics-informed artificial neural network (PIANN) that integrates a conventional neural network with a strain-based phase transition framework to predict the constitutive behaviour of amorphous SMPs. The model is validated using five temperature–stress datasets and four temperature–strain datasets, including experimental data from four types of SMPs and simulation results from a widely accepted model. PIANN predicts four key shape memory behaviours: stress evolution during hot programming, stress recovery following both cold and hot programming and free strain recovery during heating branch. Notably, it predicts recovery strain during heating without using any heating data for training. Comparisons with experimental data show excellent agreement in both programming (cooling) and recovery (heating) branches. Remarkably, the model achieves this performance with as few as two temperature–stress curves in the training set. Overall, PIANN addresses common challenges in SMP modelling by eliminating history dependence, improving curve-fitting accuracy and significantly enhancing computational efficiency. This work represents a substantial step forward in developing generalizable models for SMPs. 
    more » « less
  3. Abstract Melt pool dynamics in metal additive manufacturing (AM) is critical to process stability, microstructure formation, and final properties of the printed materials. Physics-based simulation, including computational fluid dynamics (CFD), is the dominant approach to predict melt pool dynamics. However, the physics-based simulation approaches suffer from the inherent issue of very high computational cost. This paper provides a physics-informed machine learning method by integrating the conventional neural networks with the governing physical laws to predict the melt pool dynamics, such as temperature, velocity, and pressure, without using any training data on velocity and pressure. This approach avoids solving the nonlinear Navier–Stokes equation numerically, which significantly reduces the computational cost (if including the cost of velocity data generation). The difficult-to-determine parameters' values of the governing equations can also be inferred through data-driven discovery. In addition, the physics-informed neural network (PINN) architecture has been optimized for efficient model training. The data-efficient PINN model is attributed to the extra penalty by incorporating governing PDEs, initial conditions, and boundary conditions in the PINN model. 
    more » « less
  4. Mean-field games (MFGs) provide a statistical physics-inspired modelling framework for decision-making in large populations of strategic, non-cooperative agents. Mathematically, these systems consist of a forwards–backwards in time-system of two coupled nonlinear partial differential equations (PDEs), namely, the Fokker–Plank (FP) and the Hamilton–Jacobi–Bellman (HJB) equations, governing the agent state and control distribution, respectively. In this work, we study a finite-time MFG with a rich global bifurcation structure using a reduced-order model (ROM). The ROM is a four-dimensional (4D) two-point boundary value problem (BVP) obtained by restricting the controlled dynamics to the first two moments of the agent state distribution, i.e. the mean and the variance. Phase space analysis of the ROM reveals that the invariant manifolds of periodic orbits around the so-called ‘ergodic MFG equilibrium’ play a crucial role in determining the bifurcation diagram and imparting a topological signature to various solution branches. We show a qualitative agreement of these results with numerical solutions of the full-order MFG PDE system. 
    more » « less
  5. Included in this dataset are: 1. Monte Carlo (MC) simulation data used for PDE identification 2. Molecular dynamics (MD) simulation data used for PDE identification 3. LAMMPS input file for PDE validation, Case 1: linear initial C profile 4. LAMMPS input file for PDE validation, Case 2: exponential initial C profile They are used in our manuscript, Learning the Governing PDE of Solid Diffusion from Atomistic Simulations, by Wongelemengist Nadew and Haoran Wang 
    more » « less