Training stiff neural ordinary differential equations with explicit rational Taylor series methods.
Colby Fronk1, Linda Petzold2,3
1Department of Chemical Engineering, University of California, Santa Barbara, Santa Barbara, California 93106, USA.
Chaos (Woodbury, N.Y.)
|July 18, 2025
Summary
We developed new explicit rational Taylor methods for training stiff neural ordinary differential equations. These methods offer high efficiency and stability, reducing computational costs for complex dynamics modeling.
Area of Science:
- Numerical Analysis
- Machine Learning
- Dynamical Systems
Background:
- Stiff neural ordinary differential equations (NODEs) present significant computational challenges.
- Traditional implicit methods for stiff NODEs are often computationally intensive.
- Efficient and stable training of stiff dynamics is crucial for data-driven simulations.
Purpose of the Study:
- To introduce novel explicit rational Taylor series methods for directly training stiff NODEs.
- To enhance the efficiency and numerical stability of learning stiff dynamical systems.
- To provide a computationally reduced alternative to traditional implicit methods.
Main Methods:
- Development of second and third-order explicit rational Taylor series schemes.
- Demonstration of A-stability for the proposed explicit schemes.
- Application of these methods to train stiff systems, including the van der Pol oscillator.
Main Results:
- The explicit schemes achieve high efficiency with a single linear solve per time step.
- Proposed methods exhibit strong numerical stability, even at large step sizes.
- Effective learning of stiff dynamics was demonstrated without stability issues common to implicit schemes.
- Significantly reduced computational cost compared to traditional implicit methods.
Conclusions:
- Explicit rational Taylor methods provide an efficient and stable approach for training stiff NODEs.
- These methods expand the capabilities of data-driven simulation for complex dynamics.
- The findings support applications in mesh-based simulation and physics-informed neural networks.
Related Concept Videos
Linear Approximation in Frequency Domain
136
Linear systems are characterized by two main properties: superposition and homogeneity. Superposition allows the response to multiple inputs to be the sum of the responses to each individual input. Homogeneity ensures that scaling an input by a scalar results in the response being scaled by the same scalar.
In contrast, nonlinear systems do not inherently possess these properties. However, for small deviations around an operating point, a nonlinear system can often be approximated as linear....
In contrast, nonlinear systems do not inherently possess these properties. However, for small deviations around an operating point, a nonlinear system can often be approximated as linear....
136
Linear Approximation in Time Domain
125
Nonlinear systems often require sophisticated approaches for accurate modeling and analysis, with state-space representation being particularly effective. This method is especially useful for systems where variables and parameters vary with time or operating conditions, such as in a simple pendulum or a translational mechanical system with nonlinear springs.
For a simple pendulum with a mass evenly distributed along its length and the center of mass located at half the pendulum's length,...
For a simple pendulum with a mass evenly distributed along its length and the center of mass located at half the pendulum's length,...
125
Second Order systems II
173
In an underdamped second-order system, where the damping ratio ζ is between 0 and 1, a unit-step input results in a transfer function that, when transformed using the inverse Laplace method, reveals the output response. The output exhibits a damped sinusoidal oscillation, and the difference between the input and output is termed the error signal. This error signal also demonstrates damped oscillatory behavior. Eventually, as the system reaches a steady state, the error diminishes to zero.
173
Mechanistic Models: Compartment Models in Algorithms for Numerical Problem Solving
101
Mechanistic models play a crucial role in algorithms for numerical problem-solving, particularly in nonlinear mixed effects modeling (NMEM). These models aim to minimize specific objective functions by evaluating various parameter estimates, leading to the development of systematic algorithms. In some cases, linearization techniques approximate the model using linear equations.
In individual population analyses, different algorithms are employed, such as Cauchy's method, which uses a...
In individual population analyses, different algorithms are employed, such as Cauchy's method, which uses a...
101
Transmission-Line Differential Equations
405
Transmission lines are essential components of electrical power systems. They are characterized by the distributed nature of resistance (R), inductance (L), and capacitance (C) per unit length. To analyze these lines, differential equations are employed to model the variations in voltage and current along the line.
Line Section Model
A circuit representing a line section of length Δx helps in understanding the transmission line parameters. The voltage V(x) and current i(x) are measured...
Line Section Model
A circuit representing a line section of length Δx helps in understanding the transmission line parameters. The voltage V(x) and current i(x) are measured...
405
Discrete-Time Fourier Series
368
The Discrete-Time Fourier Series (DTFS) is a fundamental concept in signal processing, serving as the discrete-time counterpart to the continuous-time Fourier series. It allows for the representation and analysis of discrete-time periodic signals in terms of their frequency components. Unlike its continuous counterpart, which utilizes integrals, the calculation of DTFS expansion coefficients involves summations due to the discrete nature of the signal.
For a discrete-time periodic signal x[n]...
For a discrete-time periodic signal x[n]...
368


