Abstract:In this paper, we consider the use of discrete gradients for differential-algebraic equations (DAEs) with a conservation/dissipation law. As one of the most popular numerical methods for conservative/dissipative ordinary differential equations, the framework of discrete gradient methods has been intensively developed over recent decades. Although discrete gradients have been applied to several specific conservative/dissipative DAEs, no unified framework for DAEs has yet been constructed. In this paper, we move toward the establishment of such a framework, and introduce concepts including an appropriate linear gradient structure for DAEs. Then, we reveal that the simple use of discrete gradients does not imply the discrete conservation/dissipation laws. Fortunately, however, we can successfully construct a new discrete gradient method for the case of index-1 DAEs. We believe this first attempt provides an indispensable basis for constructing a unified framework of discrete gradient methods for DAEs.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to establish a unified discrete - gradient method framework for differential - algebraic equations (DAEs) in order to deal with DAEs with conservation or dissipation laws. Specifically, the author hopes to overcome some limitations in existing methods and provide a new numerical integration method for this type of equations. ### Main Problem Description 1. **Lack of a Unified Framework**: Although the discrete - gradient method has been widely used in ordinary differential equations (ODEs), there is no unified framework for applying these methods to DAEs. 2. **Treatment of Conservative and Dissipative Quantities**: When dealing with DAEs with conservation or dissipation laws, the simple discrete - gradient method cannot guarantee that the discretized equations still maintain these properties. 3. **Complexity of Linear Conservation Quantities**: Unlike ODEs, in DAEs, even linear conservation quantities may be difficult to maintain in numerical methods, which increases the complexity of the problem. ### Solutions To address these problems, the author has taken the following measures: 1. **Introducing an Appropriate Linear Gradient Structure**: The author introduced a linear gradient structure suitable for DAEs, making it more convenient to analyze and design discrete - gradient methods. 2. **Proposing the Concept of "Appropriate Function"**: By defining "appropriate functions", the author can better handle the conservative and dissipative quantities in DAEs and provide simple criteria for these quantities. 3. **Constructing a New Discrete - Gradient Method**: For DAEs with an index of 1, the author successfully constructed a new discrete - gradient method that can maintain the required conserved quantities while maintaining the constraint conditions. ### Conclusions The author believes that these preliminary results lay the foundation for constructing a unified discrete - gradient method framework and are helpful for further research and development of efficient numerical methods applicable to DAEs. This contribution will have an important impact on the field of geometric numerical integration. ### Formula Summary - Linear Gradient Form: \[ \dot{z} = S(z)\nabla V(z) \] where \(S(z)\) is a skew - symmetric matrix (for conservative systems) or a negative semi - definite matrix (for dissipative systems). - Discrete - Gradient Method: \[ \frac{z^{(m + 1)} - z^{(m)}}{\Delta t}=S(z^{(m + 1)}, z^{(m)})\nabla V(z^{(m + 1)}, z^{(m)}) \] - Definition of Appropriate Function: \[ \nabla V(z)\in \text{car}(A)\quad \forall z\in M \] where \(\text{car}(A)=\text{null}(A)^{\perp}\). Through these formulas and methods, the author hopes to provide more reliable and effective tools for the numerical integration of DAEs.

Linear gradient structures and discrete gradient methods for conservative/dissipative differential-algebraic equations

A Linearly-Fitted Conservative (dissipative) Scheme for Efficiently Solving Conservative (dissipative) Nonlinear Wave PDEs

Improved structural methods for nonlinear differential-algebraic equations via combinatorial relaxation

Higher order discrete gradient method by the discontinuous Galerkin time-stepping method

Fully discrete energy-dissipative and conservative discrete gradient particle methods for a class of continuity equations

Discrete gradients in short-range molecular dynamics simulations

On the equivalence between SOR-type methods for linear systems and discrete gradient methods for gradient systems

Analysis for Implicit and Implicit-Explicit ADER and DeC Methods for Ordinary Differential Equations, Advection-Diffusion and Advection-Dispersion Equations

High-order energy stable discrete variational derivative schemes for gradient flows

A structure preserving discretization for the Derrida-Lebowitz-Speer-Spohn equation based on diffusive transport

A geometric integration approach to smooth optimisation: Foundations of the discrete gradient method

A Path-Conservative ADER Discontinuous Galerkin Method for Non-Conservative Hyperbolic Systems: Applications to Shallow Water Equations

Computational enhancement of discrete gradient method

On microscopic origins of generalized gradient structures

Energy-dissipative spectral renormalization exponential integrator method for gradient flow problems

An optimization-based approach for high-order accurate discretization of conservation laws with discontinuous solutions

Scalable Gradients for Stochastic Differential Equations

Exponential Collocation Methods for Conservative or Dissipative Systems.

Fully-discrete provably Lyapunov consistent discretizations for convection-diffusion-reaction PDE systems

Structure Preserving Primal Dual Methods for Gradient Flows with Nonlinear Mobility Transport Distances

Conservation Laws with Discontinuous Gradient-Dependent Flux: the Unstable Case