pith. sign in

arxiv: 2509.21280 · v2 · submitted 2025-09-25 · 🧮 math.NA · cs.NA

Model reduction of parametric ordinary differential equations via autoencoders: representation properties and convergence analysis

Pith reviewed 2026-05-18 13:49 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords model reductionautoencodersparametric ODEsreduced-order modelingconvergence analysisnonlinear dynamicsneural networksdimensionality reduction
0
0 comments X

The pith

Autoencoders with exact representation properties yield convergent reduced-order models for parametric nonlinear ODEs.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a reduced-order modeling method for nonlinear ordinary differential equations that depend on parameters by using autoencoders to compress the high-dimensional state into a low-dimensional latent space. The low-dimensional dynamics are then integrated with standard time-stepping schemes and the result is lifted back to the original dimension through the decoder. The authors focus on neural-network architectures that guarantee exact reconstruction of the solution manifold so that a rigorous convergence analysis between the reduced and full models can be carried out. Numerical tests in several nonlinear regimes confirm that the approach remains accurate while also revealing how the reduction step modifies the stability of the reconstructed trajectories. The work matters because it supplies a practical route to cheaper parametric studies of stiff or high-dimensional dynamical systems without sacrificing the essential nonlinear behavior.

Core claim

By training autoencoders that possess exact representation capabilities for the input manifold, the high-dimensional parametric ODE can be replaced by a low-dimensional ODE whose solutions, when reconstructed, converge to the solutions of the original system; the convergence holds under standard assumptions on the network approximation error and the time integrator.

What carries the argument

Autoencoder neural networks engineered for exact representation of the solution manifold, which map the original ODE to a low-dimensional surrogate while preserving the structure needed for convergence proofs.

If this is right

  • Standard time integrators applied to the low-dimensional ODE produce trajectories whose lift to the original space approximates the high-fidelity solution.
  • The stability properties of the reconstructed solution are directly inherited from or modified in a quantifiable way by the reduction step.
  • The method remains accurate for strongly nonlinear and parametric regimes where linear projection techniques typically degrade.
  • The same autoencoder-based reduction framework can be reused across multiple parameter values once the network is trained.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same representation-property approach could be tested on systems whose solution manifolds are known to be low-dimensional, such as certain chemical kinetics models.
  • If the exact-representation condition can be relaxed to controlled approximation error, the framework would apply to noisy or incomplete data sets.
  • Extending the convergence analysis to include discretization error in both space and time would connect the method to existing finite-element reduced-basis theory.

Load-bearing premise

The autoencoder can be chosen and trained so that it reconstructs the input manifold with exact representation capabilities.

What would settle it

A family of test problems in which the reconstructed high-dimensional solution diverges from the full-order reference even when the autoencoder reconstruction error is driven to zero would disprove the convergence claim.

Figures

Figures reproduced from arXiv: 2509.21280 by Alessio Fumagalli, Anna Scotti, Enrico Ballini, Luca Formaggia, Marco Gambarini, Paolo Zunino.

Figure 1
Figure 1. Figure 1: Example of a 1-manifold in R 3 that can be reproduced exactly with linear transformations. Black point (overlapped by red points): original manifold. Red points: output of autoencoder. Here, n = 1. 3.2.2 Linear encoder - Non-linear decoder Here, we consider a linear encoder and a non-linear decoder, with the same objective as before: having a null representation error. Proposition 4. (A linear–non-linear a… view at source ↗
Figure 2
Figure 2. Figure 2: (a) Manifold defined as graph of a function. Left panel, black dots: encoder inputs; colored dots: decoder output. The shaded surface is reconstructed with a Delaunay triangulation. The color scale is related to the z￾coordinate and it is used just to ease the graphical interpretation. Right panel: flat surface obtained from ϕ ◦ Ψ′ . (b) Coil geometry. Black dots: encoder input. Red dots: decoder output. O… view at source ↗
Figure 3
Figure 3. Figure 3: Noisy coil. A fully non-linear autoencoder in necessary to reproduce Id [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Lyapunov stability. Bounds (black lines) for the given [PITH_FULL_IMAGE:figures/full_fig_p012_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Linear case. Three components, one for each panel, of the reconstructed solutions in the test dataset. [PITH_FULL_IMAGE:figures/full_fig_p022_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Two perspectives, (a) and (b), of the curved manifold M. Black cones indicate the data used to train the autoencoder, while the colored points represent the reconstructed Ψ(un). The color scale (from blue to yellow) is related to the time, solely for the purpose ease the readability. It is worth noting that, in this test case, despite the extrapolation in time, the reconstructed solution still converges to… view at source ↗
Figure 7
Figure 7. Figure 7: 23 [PITH_FULL_IMAGE:figures/full_fig_p023_7.png] view at source ↗
Figure 7
Figure 7. Figure 7: Non linear case. Components of Ψ(un), one for each panel. First row of panels corresponds to N = 3, second and third row corresponds to N = 100, while the last two rows corresponds to N = 500. For visualization purposes, only the first six components of Ψ(un) are displayed. un is computed using FE with ∆t = 0.99∆tmax. The vertical blue-dashed line is located at the end of the training time. 24 [PITH_FULL_… view at source ↗
Figure 8
Figure 8. Figure 8: Time convergence for different strategies and integration schemes is shown. The observed orders of [PITH_FULL_IMAGE:figures/full_fig_p026_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Time evolution of the components of uN . Solid black lines represent the reference solution obtained by accurately solving the N-ODE; dashed green lines show the reconstructed solution NΨ(Un). The two solutions exhibit good agreement visually. The training data, and hence the learned manifold, are restricted to the interval t ≤ 10, which is marked by a vertical yellow dashed line. The black dashed lines in… view at source ↗
Figure 10
Figure 10. Figure 10: Set of reduced solutions, U k n , at timesteps k = 0, . . . , K, with µ sampled from the training dataset, along with the extrapolated solutions for unseen future times. 8 Conclusions In this work, we have presented a data-driven reduced-order modeling method entirely based on neural networks. The methodology can be applied to generic ODEs, such as those resulting from the spatial discretization of PDEs. … view at source ↗
read the original abstract

We propose a reduced-order modeling approach for nonlinear, parameter-dependent ordinary differential equations (ODE). Dimensionality reduction is achieved using nonlinear maps represented by autoencoders. The resulting low-dimensional ODE is then solved using standard integration in time schemes, and the high-dimensional solution is reconstructed from the low-dimensional one. We investigate the architecture of neural networks for constructing effective autoencoders that hold necessary properties to reconstruct the input manifold with exact representation capabilities. We study the convergence of the reduced-order model to the high-fidelity one. Numerical experiments show the robustness and accuracy of our approach in different scenarios, highlighting its effectiveness in highly complex and nonlinear settings without sacrificing accuracy. Moreover, we examine how the reduction influences the stability properties of the reconstructed high-dimensional solution.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper proposes a reduced-order modeling method for nonlinear parametric ODEs that uses autoencoders to perform nonlinear dimensionality reduction. The resulting low-dimensional ODE is integrated with standard time-stepping schemes and the high-dimensional solution is reconstructed via the decoder. The authors investigate autoencoder architectures that achieve exact representation of the solution manifold, derive a convergence analysis of the reduced model to the full-order parametric ODE, and present numerical experiments demonstrating accuracy and robustness in complex nonlinear regimes. They additionally examine how the reduction affects stability properties of the reconstructed trajectories.

Significance. If the convergence result can be extended to account for finite training and approximation errors, the work would supply a theoretically grounded nonlinear alternative to classical projection-based reduction for parametric ODEs. The combination of representation-property analysis, convergence statements, and stability examination is a constructive contribution; the numerical examples already illustrate practical utility in highly nonlinear settings.

major comments (2)
  1. [Convergence analysis] Convergence analysis (following the architecture discussion): the argument requires the autoencoder to possess exact representation properties so that the reduced ODE converges to the high-fidelity parametric system. No explicit error-propagation estimate is supplied that bounds the effect of residual training error, finite network width/depth, or parametric sampling on the Lipschitz constants or stability margins of the reconstructed high-dimensional trajectory. This assumption is load-bearing for the central convergence claim.
  2. [Architecture investigation] Architecture section: while the paper studies network designs that enable exact manifold reconstruction, it does not detail how the parametric dependence of the ODE is incorporated into the training loss or sampling strategy to guarantee that the representation properties hold uniformly over the parameter domain. Without such uniformity the subsequent convergence and stability statements may not transfer directly to the parametric setting.
minor comments (2)
  1. Define the precise norms and function spaces used for the error bounds and stability analysis to make the theoretical statements unambiguous.
  2. Clarify whether the reported numerical experiments employ the same autoencoder weights for all parameter values or retrain per parameter; this affects interpretation of the robustness claims.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive overall assessment and the detailed, constructive comments. We address each major comment below, indicating the revisions we will make.

read point-by-point responses
  1. Referee: [Convergence analysis] Convergence analysis (following the architecture discussion): the argument requires the autoencoder to possess exact representation properties so that the reduced ODE converges to the high-fidelity parametric system. No explicit error-propagation estimate is supplied that bounds the effect of residual training error, finite network width/depth, or parametric sampling on the Lipschitz constants or stability margins of the reconstructed high-dimensional trajectory. This assumption is load-bearing for the central convergence claim.

    Authors: We appreciate the referee highlighting this point. The convergence theorem is derived under the exact-representation assumption established in the preceding architecture analysis; this permits a clean statement that the reduced ODE converges to the full-order parametric system in the limit of increasing latent dimension. We agree that the absence of an explicit error-propagation bound for residual training error, finite network capacity, or sampling density constitutes a limitation of the current analysis. Deriving such bounds would require additional technical machinery from approximation theory and would lengthen the paper considerably. In the revised manuscript we will add a short subsection that (i) explicitly states the ideal-case nature of the result, (ii) discusses the practical implications of approximate autoencoders, and (iii) sketches a possible route toward quantitative error estimates under standard Lipschitz and boundedness assumptions on the network. This constitutes a partial revision. revision: partial

  2. Referee: [Architecture investigation] Architecture section: while the paper studies network designs that enable exact manifold reconstruction, it does not detail how the parametric dependence of the ODE is incorporated into the training loss or sampling strategy to guarantee that the representation properties hold uniformly over the parameter domain. Without such uniformity the subsequent convergence and stability statements may not transfer directly to the parametric setting.

    Authors: We thank the referee for this observation. In the current manuscript the training data consist of solution snapshots generated for a dense, uniform grid of parameter values drawn from the admissible parameter domain; the autoencoder is trained with the standard mean-squared reconstruction loss evaluated on these parametric snapshots. This procedure is intended to promote uniform representation properties across the parameter range. We acknowledge that the manuscript does not spell out the sampling density or the precise manner in which uniformity is enforced. In the revised version we will expand the architecture section with an explicit description of the parameter-sampling strategy, the construction of the training set, and a brief argument why the resulting representation properties transfer to the full parameter domain, thereby supporting the subsequent convergence and stability claims. revision: yes

Circularity Check

0 steps flagged

No significant circularity; convergence analysis rests on explicit assumption of exact autoencoder reconstruction.

full rationale

The paper states its convergence results under the assumption that the autoencoder is trained to possess exact representation properties for the solution manifold, then analyzes the reduced ODE under that hypothesis using standard numerical integration. This is a forward assumption rather than a self-referential definition or fitted quantity renamed as prediction. No load-bearing step reduces by construction to the same data or prior self-citation chain; the architecture investigation and stability examination remain independent of the target convergence claim. The numerical experiments provide external validation in nonlinear regimes, confirming the derivation chain is self-contained.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

Limited information available from abstract only. The central claim rests on the existence of autoencoder architectures that achieve exact manifold reconstruction and on standard assumptions of numerical ODE theory.

axioms (2)
  • domain assumption Autoencoders can be constructed to possess exact representation capabilities for the solution manifold of the parametric ODE.
    Invoked in the investigation of neural network architectures for reconstruction properties.
  • standard math Standard time-integration schemes applied to the reduced ODE produce solutions that converge to the high-fidelity solution under suitable conditions.
    Basis for the convergence analysis stated in the abstract.

pith-pipeline@v0.9.0 · 5672 in / 1286 out tokens · 30249 ms · 2026-05-18T13:49:10.292669+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

64 extracted references · 64 canonical work pages · 4 internal anchors

  1. [1]

    In Proceedings of the Conference Algoritmy, pages 1–12, 2016

    Reduced basis methods: success, limitations and future challenges. In Proceedings of the Conference Algoritmy, pages 1–12, 2016

  2. [2]

    Ansel, E

    J. Ansel, E. Yang, H. He, et al. PyTorch 2: Faster Machine Learning Through Dynamic Python Bytecode Transformation and Graph Compilation. In 29th ACM International Conference on Architectural Support for Programming Languages and Operating Systems, Volume 2 (ASPLOS ’24). ACM, 2024

  3. [3]

    Ballini, A

    E. Ballini, A. Cominelli, L. Dovera, A. Forello, L. Formaggia, A. Fumagalli, S. Nardean, A. Scotti, and P. Zunino. Enhancing computational efficiency of numerical simulation for subsurface fluid-induced deformation using deep learning reduced order models. In SPE Reservoir Simulation Conference, 25RSC. SPE, Mar. 2025

  4. [4]

    J. C. Butcher. Numerical Methods for Ordinary Differential Equations. Wiley, July 2016

  5. [5]

    Y. Cai. Achieve the minimum width of neural networks for universal approximation, 2022. arXiv:2209.11395

  6. [6]

    P. Y. Chen, J. Xiang, D. H. Cho, Y. Chang, G. A. Pershing, H. T. Maia, M. M. Chiaramonte, K. T. Carlberg, and E. Grinspun. CROM: Continuous reduced-order modeling of PDEs using implicit neural representations. In The Eleventh International Conference on Learning Representations, 2023

  7. [7]

    R. T. Q. Chen, B. Amos, and M. Nickel. Learning neural event functions for ordinary differential equations

  8. [8]

    arXiv:2011.03902. 29

  9. [9]

    R. T. Q. Chen, Y. Rubanova, J. Bettencourt, and D. K. Duvenaud. Neural ordinary differential equations. In S. Bengio, H. Wallach, H. Larochelle, K. Grauman, N. Cesa-Bianchi, and R. Garnett, editors, Advances in Neural Information Processing Systems, volume 31. Curran Associates, Inc., 2018

  10. [10]

    Conti, M

    P. Conti, M. Guo, A. Manzoni, A. Frangi, S. L. Brunton, and J. Nathan Kutz. Multi-fidelity reduced-order surrogate modelling. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, 480(2283), 2024

  11. [11]

    Q. Du, Y. Gu, H. Yang, and C. Zhou. The discovery of dynamics via linear multistep methods and deep learning: Error estimation. SIAM Journal on Numerical Analysis, 60(4):2014–2045, 2022

  12. [12]

    D’Ambrosio and S

    R. D’Ambrosio and S. D. Giovacchino. Mean-square contractivity of stochastic theta-methods.Communications in Nonlinear Science and Numerical Simulation, 96:105671, May 2021

  13. [13]

    Farenga, S

    N. Farenga, S. Fresca, S. Brivio, and A. Manzoni. On latent dynamics learning in nonlinear reduced order modeling. Neural Networks, 185:107146, May 2025

  14. [14]

    Floryan and M

    D. Floryan and M. D. Graham. Data-driven discovery of intrinsic dynamics. Nature Machine Intelligence, 4(12):1113–1120, Dec. 2022

  15. [15]

    N. R. Franco, S. Fresca, F. Tombari, and A. Manzoni. Deep learning-based surrogate models for parametrized pdes: Handling geometric variability through graph neural networks. Chaos: An Interdisciplinary Journal of Nonlinear Science, 33(12), 2023

  16. [16]

    N. R. Franco, A. Manzoni, and P. Zunino. A deep learning approach to reduced order modelling of parameter dependent partial differential equations. Mathematics of Computation, 92(340):483–524, 2022

  17. [17]

    N. R. Franco, A. Manzoni, and P. Zunino. Mesh-informed neural networks for operator learning in finite element spaces. Journal of Scientific Computing, 97(2), 2023

  18. [18]

    Fresca, L

    S. Fresca, L. Dede’, and A. Manzoni. A comprehensive deep learning-based approach to reduced order modeling of nonlinear time-dependent parametrized PDEs. Journal of Scientific Computing, 87(2), 2021

  19. [19]

    Fresca, F

    S. Fresca, F. Fatone, and A. Manzoni. Long-time prediction of nonlinear parametrized dynamical systems by deep learning-based reduced order models. Mathematics in Engineering, 5(6):1–36, 2023

  20. [20]

    Fresca and A

    S. Fresca and A. Manzoni. POD-DL-ROM: Enhancing deep learning-based reduced order models for nonlin- ear parametrized PDEs by proper orthogonal decomposition. Computer Methods in Applied Mechanics and Engineering, 388:114181, 2022

  21. [21]

    R. Fu, D. Xiao, I. Navon, F. Fang, L. Yang, C. Wang, and S. Cheng. A non-linear non-intrusive reduced order model of fluid flow by auto-encoder and self-attention deep learning methods. International Journal for Numerical Methods in Engineering, 124(13):3087–3111, 2023

  22. [22]

    Fumagalli and A

    A. Fumagalli and A. Scotti. A mathematical model for thermal single-phase flow and reactive transport in fractured porous media. Journal of Computational Physics, 434:110205, 2021

  23. [23]

    F. J. Gonzalez and M. Balajewicz. Deep convolutional recurrent autoencoders for learning low-dimensional feature dynamics of fluid systems. arXiv:1808.01346, 2018

  24. [24]

    Goodfellow, Y

    I. Goodfellow, Y. Bengio, and A. Courville. Deep Learning. MIT Press, 2017

  25. [25]

    Gupta and R

    R. Gupta and R. Jaiman. Three-dimensional deep learning-based reduced order model for unsteady flow dynamics with variable reynolds number. Physics of Fluids, 34(3), 2022

  26. [26]

    G¨ uhring and M

    I. G¨ uhring and M. Raslan. Approximation rates for neural networks with encodable weights in smoothness spaces. Neural Networks, 134:107–130, Feb. 2021

  27. [27]

    Hairer and G

    E. Hairer and G. Wanner. Solving Ordinary Differential Equations II. Springer Berlin Heidelberg, 1996. 30

  28. [28]

    Approximating Continuous Functions by ReLU Nets of Minimal Width

    B. Hanin and M. Sellke. Approximating continuous functions by ReLU nets of minimal width.arXiv:1710.11278, 2017

  29. [29]

    Hasegawa, K

    K. Hasegawa, K. Fukami, T. Murata, and K. Fukagata. Machine-learning-based reduced-order modeling for un- steady flows around bluff bodies of various shapes. Theoretical and Computational Fluid Dynamics, 34(4):367– 383, 2020

  30. [30]

    K. He, X. Zhang, S. Ren, and J. Sun. Deep residual learning for image recognition. In 2016 IEEE Conference on Computer Vision and Pattern Recognition (CVPR), pages 770–778, 2016

  31. [31]

    J. S. Hesthaven, G. Rozza, and B. Stamm. Certified Reduced Basis Methods for Parametrized Partial Differential Equations. Springer International Publishing, 2016

  32. [32]

    Kadeethum, F

    T. Kadeethum, F. Ballarin, Y. Choi, D. O’Malley, H. Yoon, and N. Bouklas. Non-intrusive reduced order modeling of natural convection in porous media using convolutional autoencoders: Comparison with linear subspace techniques. Advances in Water Resources, 160:104098, 2022

  33. [33]

    P. Kidger. On Neural Differential Equations. Ph.D. thesis, University of Oxford, 2021

  34. [34]

    D. P. Kingma and J. Ba. Adam: A method for stochastic optimization. arXiv:1412.6980, 2014

  35. [35]

    Kumar, R

    A. Kumar, R. Hu, and S. D. C. Walsh. Development of reduced order hydro-mechanical models of fractured media. Rock Mechanics and Rock Engineering, 55(1):235–248, 2021

  36. [36]

    J. D. Lambert. Numerical methods for ordinary differential systems. Wiley, 1999

  37. [37]

    J. M. Lee. Introduction to Topological Manifolds. Springer New York, 2000

  38. [38]

    Lee and K

    K. Lee and K. T. Carlberg. Model reduction of dynamical systems on nonlinear manifolds using deep convo- lutional autoencoders. Journal of Computational Physics, 404:108973, 2020

  39. [39]

    Legaard, T

    C. Legaard, T. Schranz, G. Schweiger, J. Drgoˇ na, B. Falay, C. Gomes, A. Iosifidis, M. Abkar, and P. Larsen. Constructing neural network based models for simulating dynamical systems. ACM Computing Surveys, 55(11):1–34, Feb. 2023

  40. [40]

    F. Li, Y. Zhang, and C. Xiao. Surrogate-based hydraulic fracture propagation prediction using deep neural network proxy. In 58th U.S. Rock Mechanics/Geomechanics Symposium, ARMA24. ARMA, June 2024

  41. [41]

    L. Li, Y. Duan, G. Ji, and Y. Cai. Minimum width of leaky-relu neural networks for uniform universal approximation, 2023. arXiv:2305.18460

  42. [42]

    Z. Li, N. Kovachki, K. Azizzadenesheli, B. Liu, K. Bhattacharya, A. Stuart, and A. Anandkumar. Fourier neural operator for parametric partial differential equations. arXiv:2010.08895, 2020

  43. [43]

    Z. Li, S. Patil, F. Ogoke, D. Shu, W. Zhen, M. Schneier, J. R. Buchanan, and A. Barati Farimani. Latent neural pde solver: A reduced-order modeling framework for partial differential equations. Journal of Computational Physics, 524:113705, Mar. 2025

  44. [44]

    J. Lu, Z. Shen, H. Yang, and S. Zhang. Deep network approximation for smooth functions. SIAM Journal on Mathematical Analysis, 53(5):5465–5506, Jan. 2021

  45. [45]

    L. Lu, P. Jin, G. Pang, Z. Zhang, and G. E. Karniadakis. Learning nonlinear operators via deeponet based on the universal approximation theorem of operators. Nature Machine Intelligence, 3(3):218–229, Mar. 2021

  46. [46]

    Y. Lu, A. Zhong, Q. Li, and B. Dong. Beyond finite layer neural networks: Bridging deep architectures and numerical differential equations. In J. Dy and A. Krause, editors, Proceedings of the 35th International Conference on Machine Learning, volume 80 of Proceedings of Machine Learning Research, pages 3276–3285. PMLR, 10–15 Jul 2018

  47. [47]

    J. E. Marsden. Mathematical Foundations of Elasticity. Dover Civil and Mechanical Engineering. Dover Publications, Newburyport, 2012. Description based upon print version of record. 31

  48. [48]

    J. Munkres. Topology. Featured Titles for Topology. Prentice Hall, Incorporated, 2000

  49. [49]

    Murata, K

    T. Murata, K. Fukami, and K. Fukagata. Nonlinear mode decomposition with convolutional neural networks for fluid dynamics. Journal of Fluid Mechanics, 882, 2019

  50. [50]

    J. M. Nordbotten and M. A. Celia. Geological storage of CO2 modeling approaches for large-scale simulation. Wiley, 2012

  51. [51]

    P. Pant, R. Doshi, P. Bahl, and A. Barati Farimani. Deep learning for reduced order modelling and efficient temporal evolution of fluid simulations. Physics of Fluids, 33(10), 2021

  52. [52]

    S. Park, C. Yun, J. Lee, and J. Shin. Minimum width for universal approximation. arXiv:2006.08859, 2020

  53. [53]

    Quarteroni, A

    A. Quarteroni, A. Manzoni, and F. Negri. Reduced Basis Methods for Partial Differential Equations. Springer International Publishing, 2016

  54. [54]

    Quarteroni, R

    A. Quarteroni, R. Sacco, and F. Saleri. Numerical Mathematics. Number 37 in Texts in Applied Mathematics. Springer, Berlin, second edition edition, 2007. Literaturverzeichnis: Seite 635-644

  55. [55]

    Raissi, P

    M. Raissi, P. Perdikaris, and G. Karniadakis. Physics-informed neural networks: A deep learning frame- work for solving forward and inverse problems involving nonlinear partial differential equations. Journal of Computational Physics, 378:686–707, 2019

  56. [56]

    Regazzoni, S

    F. Regazzoni, S. Pagani, M. Salvador, L. Dede’, and A. Quarteroni. Learning the intrinsic dynamics of spatio- temporal processes through latent dynamics networks. Nature Communications, 15(1), Feb. 2024

  57. [57]

    Ronneberger, P

    O. Ronneberger, P. Fischer, and T. Brox. U-Net: Convolutional Networks for Biomedical Image Segmentation, pages 234–241. Springer International Publishing, 2015

  58. [58]

    Rubanova, R

    Y. Rubanova, R. T. Q. Chen, and D. Duvenaud. Latent ODEs for irregularly-sampled time series. Curran Associates Inc., Red Hook, NY, USA, 2019

  59. [59]

    Sabzi Shahrebabaki, S

    A. Sabzi Shahrebabaki, S. Fouladi, E. Holtar, and L. Vynnytska. Autoencoder-based generation of subsurface models. In Second EAGE Digitalization Conference and Exhibition, pages 1–5. European Association of Geoscientists & Engineers, 2022

  60. [60]

    C. I. Steefel and A. C. Lasaga. A coupled model for transport of multiple chemical species and kinetic precipitation/dissolution reactions with application to reactive flow in single phase hydrothermal systems. American Journal of Science, 294(5):529–592, May 1994

  61. [61]

    P. A. Thompson. Compressible-Fluid Dynamics. Advanced engineering series. McGraw-Hill, Inc., 1988

  62. [62]

    Tripura and S

    T. Tripura and S. Chakraborty. Wavelet neural operator for solving parametric partial differential equations in computational mechanics problems. Computer Methods in Applied Mechanics and Engineering, 404:115783, Feb. 2023

  63. [63]

    Virtanen, R

    P. Virtanen, R. Gommers, T. E. Oliphant, et al. SciPy 1.0: fundamental algorithms for scientific computing in Python. Nature Methods, 17(3):261–272, 2020

  64. [64]

    Linear(Affine)x×y

    G. Ziarelli, N. Parolini, and M. Verani. Learning epidemic trajectories through kernel operator learning: From modelling to optimal control. Numerical Mathematics: Theory, Methods and Applications, Apr. 2025. 32 Appendix A Lipschitz constants and stability Lipschitz constantL Fn.LetM y(x) = supx ∥y(x)∥andL y the Lipschitz constant of the generic functiony...