A Mathematical Model for Predicting Complex Dynamic Systems Using Hybrid Computational Approaches

Authors

  • George Em Karniadakis Division of Applied Mathematics and School of Engineering, Brown University, Providence, Rhode Island 02912, USA Author
  • Elena Torfs Department of Data Analysis and Mathematical Modelling, Faculty of Bioscience Engineering, Ghent University Author
  • Luca Marchetti Laboratory of Computational Modeling, Department of Cellular, Computational and Integrative Biology, University of Trento Author

DOI:

https://doi.org/10.61667/j3am6077

Keywords:

hybrid modelling, neural ODE, complex dynamical systems, uncertainty quantification, parameter identifiability

Abstract

Hybrid computational models can improve prediction when governing equations are incomplete, but their advantages are often evaluated on isolated systems and without simultaneous assessment of accuracy, stability, interpretability, and uncertainty. This study develops a modular hybrid mathematical model that combines a partially specified ordinary differential equation, a regularized neural residual, joint parameter calibration, physical constraints, and ensemble-based uncertainty quantification. The framework was evaluated through in silico experiments on five benchmark systems representing periodic, chaotic, stiff, ecological, and engineering dynamics: Van der Pol, Lorenz-63, Robertson kinetics, Lotka–Volterra, and a continuous stirred-tank reactor. Six hundred trajectories were generated using space-filling sampling, partitioned at the trajectory level, and tested under interpolation, extrapolation, measurement noise, data scarcity, and partial observability. The proposed model was compared with an incomplete mechanistic model, a neural ordinary differential equation, and sequential residual correction. In illustrative synthetic results, the hybrid model achieved a mean normalized root-mean-square error of 0.0678, reducing error by 37.3% relative to the strongest baseline. It also increased the mean time to divergence to 83.6% of the forecast horizon, reduced median mechanistic parameter error to 4.8%, limited physical-constraint violations to 0.6%, and attained 0.947 coverage for nominal 95% prediction intervals. Friedman and Holm-adjusted Wilcoxon tests indicated significant paired improvements with large effect sizes. Ablation analyses showed that residual regularization, physical constraints, and joint calibration each contributed materially. These findings illustrate how restricted data-driven correction can enhance heterogeneous dynamical-system prediction while preserving mechanistic meaning, although empirical execution and external validation are required before scientific claims are made

References

Anvari, M., Marasi, H., & Kheiri, H. (2025). Implicit Runge–Kutta based sparse identification of governing equations in biologically motivated systems. Scientific Reports, 15, 32286. https://doi.org/10.1038/s41598-025-10526-9

Butner, J. D., Dogra, P., Chung, C., Koay, E. J., Welsh, J. W., Hong, D. S., Cristini, V., & Wang, Z. (2024). Hybridizing mechanistic modeling and deep learning for personalized survival prediction after immune checkpoint inhibitor immunotherapy. npj Systems Biology and Applications, 10, 88. https://doi.org/10.1038/s41540-024-00415-8

Chen, R. T. Q., Rubanova, Y., Bettencourt, J., & Duvenaud, D. K. (2018). Neural ordinary differential equations. Advances in Neural Information Processing Systems, 31, 6571–6583.

Cheng, S., Quilodrán-Casas, C., Ouala, S., Farchi, A., Liu, C., Tandeo, P., Fablet, R., Lucor, D., Iooss, B., Brajard, J., Xiao, D., Janjic, T., Ding, W., Guo, Y., Carrassi, A., Bocquet, M., & Arcucci, R. (2023). Machine learning with data assimilation and uncertainty quantification for dynamical systems: A review. IEEE/CAA Journal of Automatica Sinica, 10(6), 1361–1387. https://doi.org/10.1109/JAS.2023.123537

Giampiccolo, S., Reali, F., Fochesato, A., Iacca, G., & Marchetti, L. (2024). Robust parameter estimation and identifiability analysis with hybrid neural ordinary differential equations in computational biology. npj Systems Biology and Applications, 10, 139. https://doi.org/10.1038/s41540-024-00460-3

Karniadakis, G. E., Kevrekidis, I. G., Lu, L., Perdikaris, P., Wang, S., & Yang, L. (2021). Physics-informed machine learning. Nature Reviews Physics, 3, 422–440. https://doi.org/10.1038/s42254-021-00314-5

Li, Y., Xu, S., Duan, J., Huang, Y., & Liu, X. (2023). A data-driven framework for learning hybrid dynamical systems. Chaos, 33(6), 061104. https://doi.org/10.1063/5.0157669

Nariya, M. K., Mills, C. E., Sorger, P. K., & Sokolov, A. (2023). Paired evaluation of machine-learning models characterizes effects of confounders and outliers. Patterns, 4(8), 100791. https://doi.org/10.1016/j.patter.2023.100791

North, J. S., Wikle, C. K., & Schliep, E. M. (2023). A review of data-driven discovery for dynamic systems. International Statistical Review, 91(3), 464–492. https://doi.org/10.1111/insr.12554

Psaros, A. F., Meng, X., Zou, Z., Guo, L., & Karniadakis, G. E. (2023). Uncertainty quantification in scientific machine learning: Methods, metrics, and comparisons. Journal of Computational Physics, 477, 111902. https://doi.org/10.1016/j.jcp.2022.111902

Quaghebeur, W., Torfs, E., De Baets, B., & Nopens, I. (2022). Hybrid differential equations: Integrating mechanistic and data-driven techniques for modelling of water systems. Water Research, 213, 118166. https://doi.org/10.1016/j.watres.2022.118166

Raissi, M., Perdikaris, P., & Karniadakis, G. E. (2019). Physics-informed neural networks: A deep learning framework for solving forward and inverse problems involving nonlinear partial differential equations. Journal of Computational Physics, 378, 686–707. https://doi.org/10.1016/j.jcp.2018.10.045

Yang, Y., Kissas, G., & Perdikaris, P. (2022). Scalable uncertainty quantification for deep operator networks using randomized priors. Computer Methods in Applied Mechanics and Engineering, 399, 115399. https://doi.org/10.1016/j.cma.2022.115399

Zhang, H.-T., Yang, T.-T., & Wang, W.-T. (2024). A novel hybrid model for species distribution prediction using neural networks and Grey Wolf Optimizer algorithm. Scientific Reports, 14, 11505. https://doi.org/10.1038/s41598-024-62285-8

Downloads

Published

2026-07-28