TY - INPR U1 - Preprint A1 - Thiele, Michael A1 - Tscherpel, Tabea T1 - Asymptotic preserving discretisation and error estimates for compressible gas mixture models N2 - We consider a compressible mixture model for binary fluids with two velocities in one space dimension in the low-Mach, high-friction regime. Building on the framework in [Egger, Giesselmann 2023], we establish stability with respect to the data. Furthermore, we analyse a high-order, energy-consistent numerical discretisation. In particular, we prove asymptotic-preserving error estimates for a first-order time discretisation combined with a spatial discretisation of arbitrary order, under suitable regularity assumptions. The scheme combines a modified continuous Petrov-Galerkin discretisation in time, introduced in [Giesselmann, Karsai, Tscherpel 2025] for a general class of nonlinear port-Hamiltonian systems, with a mixed finite element discretisation in space, both of arbitrary order. By developing and using an approximation operator adapted to the nonlinear structure of the model, we obtain optimal-order error estimates in the spatial discretisation parameter. Our results also apply to a one-component fluid with friction, extending existing schemes and error estimates to higher-order spatial discretisations. AB - We consider a compressible mixture model for binary fluids with two velocities in one space dimension in the low-Mach, high-friction regime. Building on the framework in [Egger, Giesselmann 2023], we establish stability with respect to the data. Furthermore, we analyse a high-order, energy-consistent numerical discretisation. In particular, we prove asymptotic-preserving error estimates for a first-order time discretisation combined with a spatial discretisation of arbitrary order, under suitable regularity assumptions. The scheme combines a modified continuous Petrov-Galerkin discretisation in time, introduced in [Giesselmann, Karsai, Tscherpel 2025] for a general class of nonlinear port-Hamiltonian systems, with a mixed finite element discretisation in space, both of arbitrary order. By developing and using an approximation operator adapted to the nonlinear structure of the model, we obtain optimal-order error estimates in the spatial discretisation parameter. Our results also apply to a one-component fluid with friction, extending existing schemes and error estimates to higher-order spatial discretisations. KW - continuous Petrov-Galerkin KW - a priori error estimates KW - port-Hamiltonian KW - mixtures KW - energy-consistency Y1 - 2026 ER - TY - INPR U1 - Preprint A1 - Hante, Falk A1 - Hernandez, Martin A1 - Leugering, Günter A1 - Schmidt, Martin A1 - Zuazua, Enrique T1 - Domain Decomposition for Gas Network Control N2 - We present domain decomposition techniques for efficient and scalable simulation and optimization of partial differential equations on network domains. Motivated by transient gas network applications, the state-of-the-art of space-and-time decompositions for systems derived from Euler's equations are presented. Moreover, combinatorial aspects of switching valves are addressed in a decomposition framework for mixed-integer and ODE-constrained problems and a convergence analysis is presented. Motivated by stochastic approximation, we also propose and explore randomized decomposition approaches. Some of these advanced solution strategies can be combined and regarded as variants of certain penalty alternating direction methods. We present the potential applications of these techniques by summarizing numerical results from the literature for transient GasLib instances modeling a realistic gas transport network. AB - We present domain decomposition techniques for efficient and scalable simulation and optimization of partial differential equations on network domains. Motivated by transient gas network applications, the state-of-the-art of space-and-time decompositions for systems derived from Euler's equations are presented. Moreover, combinatorial aspects of switching valves are addressed in a decomposition framework for mixed-integer and ODE-constrained problems and a convergence analysis is presented. Motivated by stochastic approximation, we also propose and explore randomized decomposition approaches. Some of these advanced solution strategies can be combined and regarded as variants of certain penalty alternating direction methods. We present the potential applications of these techniques by summarizing numerical results from the literature for transient GasLib instances modeling a realistic gas transport network. Y1 - 2026 ER - TY - INPR U1 - Preprint A1 - Breitkopf, Jannik A1 - Ulbrich, Stefan T1 - An Adjoint Calculus for Optimal Control of Entropy Solutions of the Generalized Riemann Problem for Hyperbolic Systems of Conservation Laws with Shocks N2 - We derive an adjoint gradient representation for the objective functional of an optimal control problem of the unique entropy solution to the Generalized Riemann Problem (GRP) for strictly hyperbolic systems of conservation or balance laws. The GRP is an initial value problem where the initial state is piecewise C^1 with exactly one discontinuity. In a recent work, differentiability properties of the solution operator of the GRP were derived. In particular, the differentiability of a class of tracking-type functionals was shown. We build on these results to derive an adjoint gradient representation for such functionals. The adjoint problem is a system of linear transport equations with a discontinuous coefficient and discontinuous terminal state. We prove the existence and uniqueness of a piecewise Lipschitz continuous solution to the adjoint problem which satisfies suitable interior boundary conditions along all shock curves. The adjoint state generally contains infinitely many discontinuities but its total variation remains uniformly bounded in time. AB - We derive an adjoint gradient representation for the objective functional of an optimal control problem of the unique entropy solution to the Generalized Riemann Problem (GRP) for strictly hyperbolic systems of conservation or balance laws. The GRP is an initial value problem where the initial state is piecewise C^1 with exactly one discontinuity. In a recent work, differentiability properties of the solution operator of the GRP were derived. In particular, the differentiability of a class of tracking-type functionals was shown. We build on these results to derive an adjoint gradient representation for such functionals. The adjoint problem is a system of linear transport equations with a discontinuous coefficient and discontinuous terminal state. We prove the existence and uniqueness of a piecewise Lipschitz continuous solution to the adjoint problem which satisfies suitable interior boundary conditions along all shock curves. The adjoint state generally contains infinitely many discontinuities but its total variation remains uniformly bounded in time. KW - optimal control KW - adjoint calculus KW - hyperbolic systems of conservation laws Y1 - 2026 SP - 50 ER - TY - JOUR U1 - Wissenschaftlicher Artikel A1 - Karsai, Attila A1 - Schulze, Philipp T1 - A discrete gradient scheme for preserving QSR-dissipativity N2 - The notion of dissipative dynamical systems provides a formal description of processes that cannot generate energy internally. For these systems, changes in energy can only occur due to an external energy supply or dissipation effects. Unfortunately, dissipative properties tend to deteriorate in numerical computations, especially in nonlinear systems. Discrete gradient methods can help mitigate this problem. In this paper, we present a class of structure-preserving time discretization schemes based on discrete gradients for a special class of systems that are dissipative with respect to a quadratic supply rate. AB - The notion of dissipative dynamical systems provides a formal description of processes that cannot generate energy internally. For these systems, changes in energy can only occur due to an external energy supply or dissipation effects. Unfortunately, dissipative properties tend to deteriorate in numerical computations, especially in nonlinear systems. Discrete gradient methods can help mitigate this problem. In this paper, we present a class of structure-preserving time discretization schemes based on discrete gradients for a special class of systems that are dissipative with respect to a quadratic supply rate. Y1 - 2026 U6 - https://doi.org/10.48550/arXiv.2602.15445 DO - https://doi.org/10.48550/arXiv.2602.15445 ER - TY - JOUR U1 - Wissenschaftlicher Artikel A1 - Altmann, Robert A1 - Karsai, Attila A1 - Schulze, Philipp T1 - Structure-Preserving Discretization and Model Reduction for Energy-Based Models N2 - We investigate discretization strategies for a recently introduced class of energy-based models. The model class encompasses classical port-Hamiltonian systems, generalized gradient flows, and certain systems with algebraic constraints. Our framework combines existing ideas from the literature and systematically addresses temporal discretization, spatial discretization, and model order reduction, ensuring that all resulting schemes are dissipation-preserving in the sense of a discrete dissipation inequality. For this, we use a Petrov-Galerkin ansatz together with appropriate projections. Numerical results for a nonlinear circuit model and the Cahn-Hilliard equation illustrate the effectiveness of the approach. AB - We investigate discretization strategies for a recently introduced class of energy-based models. The model class encompasses classical port-Hamiltonian systems, generalized gradient flows, and certain systems with algebraic constraints. Our framework combines existing ideas from the literature and systematically addresses temporal discretization, spatial discretization, and model order reduction, ensuring that all resulting schemes are dissipation-preserving in the sense of a discrete dissipation inequality. For this, we use a Petrov-Galerkin ansatz together with appropriate projections. Numerical results for a nonlinear circuit model and the Cahn-Hilliard equation illustrate the effectiveness of the approach. Y1 - 2026 U6 - https://doi.org/10.48550/arXiv.2507.21552 DO - https://doi.org/10.48550/arXiv.2507.21552 ER - TY - THES U1 - Dissertation oder Habilitation A1 - Karsai, Attila T1 - Nonlinear energy-based systems: modeling, control and numerical realization N2 - Dynamical systems consist of ordinary and partial differential equations and are among the most prominent approaches to model physical processes. They describe the evolution of the system in terms of the current system state and external control inputs. To improve model accuracy, it can be beneficial to explicitly include properties such as energy conservation or energy dissipation, which are usually present in the real-world phenomena, in the mathematical problem description. Ideally, these properties should be kept in mind whenever one interacts with the model. The focus of this thesis is on two kinds of interactions: the discretization of such models, and their control. Discretization techniques are necessary whenever the system evolution is to be approximated using computational methods. Similarly fundamental is the numerical realization of control inputs that lead to desired outcomes. Both aspects require special attention to retain an energy-based perspective. For this, the first step is usually to encode the energy properties in the algebraic description of the model. This description must strike a balance between general applicability and its corresponding benefits. The next step is to leverage the algebraic description in further analysis. While the field is well-developed for linear systems, the nonlinear case often poses additional difficulties. First, there seems to be no clear consensus about what model class to use to describe nonlinear physical phenomena. Second, many discretization methods that preserve the energy-based viewpoint in the linear case are not trivial to generalize to nonlinear systems. Third, studying the behavior of energy-optimal controls and finding control laws that can be realized via energy-based models is usually more difficult. In this thesis, these points are addressed. We give an overview of energy-based model classes used for nonlinear phenomena in both finite and infinite dimensions and investigate their relationship. Moreover, using a modified Petrov--Galerkin method and discrete gradients, we present multiple structure-preserving discretization schemes. Furthermore, we show that similar to the linear case, energy-optimal controls steer the associated trajectories to the submanifold of the state space where no dissipation is present. Finally, we combine an optimal feedback law characterized by the Hamilton--Jacobi--Bellman equation with output feedback to state an energy-based feedback controller. Our theoretical results are illustrated using numerical experiments. AB - Dynamical systems consist of ordinary and partial differential equations and are among the most prominent approaches to model physical processes. They describe the evolution of the system in terms of the current system state and external control inputs. To improve model accuracy, it can be beneficial to explicitly include properties such as energy conservation or energy dissipation, which are usually present in the real-world phenomena, in the mathematical problem description. Ideally, these properties should be kept in mind whenever one interacts with the model. The focus of this thesis is on two kinds of interactions: the discretization of such models, and their control. Discretization techniques are necessary whenever the system evolution is to be approximated using computational methods. Similarly fundamental is the numerical realization of control inputs that lead to desired outcomes. Both aspects require special attention to retain an energy-based perspective. For this, the first step is usually to encode the energy properties in the algebraic description of the model. This description must strike a balance between general applicability and its corresponding benefits. The next step is to leverage the algebraic description in further analysis. While the field is well-developed for linear systems, the nonlinear case often poses additional difficulties. First, there seems to be no clear consensus about what model class to use to describe nonlinear physical phenomena. Second, many discretization methods that preserve the energy-based viewpoint in the linear case are not trivial to generalize to nonlinear systems. Third, studying the behavior of energy-optimal controls and finding control laws that can be realized via energy-based models is usually more difficult. In this thesis, these points are addressed. We give an overview of energy-based model classes used for nonlinear phenomena in both finite and infinite dimensions and investigate their relationship. Moreover, using a modified Petrov--Galerkin method and discrete gradients, we present multiple structure-preserving discretization schemes. Furthermore, we show that similar to the linear case, energy-optimal controls steer the associated trajectories to the submanifold of the state space where no dissipation is present. Finally, we combine an optimal feedback law characterized by the Hamilton--Jacobi--Bellman equation with output feedback to state an energy-based feedback controller. Our theoretical results are illustrated using numerical experiments. Y1 - 2026 U6 - https://doi.org/10.14279/depositonce-25690 DO - https://doi.org/10.14279/depositonce-25690 ER - TY - INPR U1 - Preprint A1 - Geiselmann, Zoe A1 - Joswig, Michael A1 - Kastner, Lars A1 - Mundinger, Konrad A1 - Pokutta, Sebastian A1 - Spiegel, Christoph A1 - Wack, Marcel A1 - Zimmer, Max T1 - 121 Patchworked Curves of Degree Seven N2 - The 121 real schemes, i.e., ambient isotopy classes, of smooth real plane algebraic curves of degree seven were classified by Viro (1984). By constructing one patchwork of the dilated triangle 7⋅Δ2 for each real scheme, we provide an explicit method for constructing polynomials realizing each real scheme. In particular, every real scheme of degree seven can be realized as a T-curve; this settles a question raised by Itenberg and Viro (1996). AB - The 121 real schemes, i.e., ambient isotopy classes, of smooth real plane algebraic curves of degree seven were classified by Viro (1984). By constructing one patchwork of the dilated triangle 7⋅Δ2 for each real scheme, we provide an explicit method for constructing polynomials realizing each real scheme. In particular, every real scheme of degree seven can be realized as a T-curve; this settles a question raised by Itenberg and Viro (1996). Y1 - 2026 ER - TY - INPR U1 - Preprint A1 - Geiselmann, Zoe A1 - Joswig, Michael A1 - Kastner, Lars A1 - Mundinger, Konrad A1 - Pokutta, Sebastian A1 - Spiegel, Christoph A1 - Wack, Marcel A1 - Zimmer, Max T1 - Fast Isotopy Computation for T-Curves N2 - A T-curve of degree d is given by a regular unimodular triangulation of d⋅Δ2 together with a sign distribution on its lattice points. By Viro's Patchworking Theorem, this determines the ambient isotopy type (a.k.a. real scheme) of a smooth real plane projective algebraic curve of the same degree. We present a near-quadratic time algorithm for extracting that isotopy type from the triangulation and the signs. Through a GPU-accelerated implementation, this allows one to compute billions of real schemes per second, enabling exhaustive enumeration at scale. This algorithm was essential for our recent construction of all 121 real schemes of degree seven by T-curves. AB - A T-curve of degree d is given by a regular unimodular triangulation of d⋅Δ2 together with a sign distribution on its lattice points. By Viro's Patchworking Theorem, this determines the ambient isotopy type (a.k.a. real scheme) of a smooth real plane projective algebraic curve of the same degree. We present a near-quadratic time algorithm for extracting that isotopy type from the triangulation and the signs. Through a GPU-accelerated implementation, this allows one to compute billions of real schemes per second, enabling exhaustive enumeration at scale. This algorithm was essential for our recent construction of all 121 real schemes of degree seven by T-curves. Y1 - 2026 ER - TY - INPR U1 - Preprint A1 - Göß, Adrian T1 - Clash of MINLP Relaxations: Piecewise Linear vs. Global Parabolic N2 - Solving mixed-integer nonlinear programs (MINLPs) typically relies on constructing relaxations that are easier to tackle than the original problem. Recently, global parabolic (PARA) relaxations were introduced, featuring separable quadratic functions – paraboloids – as global under- or overestimators of general nonlinear constraint functions. So far, the paraboloids are all computed at once by solving a mixed-integer linear program (MIP). For small tolerances or wide function domains, the corresponding MIP grows in size and is eventually intractable, preventing a meaningful comparison with established relaxation techniques. We therefore propose a novel iterative method to compute PARA approximations that succeeds on all tolerance-domain combinations where the original one has failed. The computational study is preceded by a thorough theoretical explanation and analysis. Finally, the improved method enables a computational comparison with piecewise linear (PWL) relaxations in terms of runtime on general MINLP instances. The results show that the modern solver SCIP can solve PWL relaxations faster when the olerance is high, shifting strongly in favor of PARA for tighter tolerances. We attribute the effect to the difference in the corresponding problem size: PWL relaxations introduce binary variables to identify the active linear piece and their number grows with decreasing tolerance. PARA, on the other hand, does not require additional variables such that the dimension is maintained. For problems with at least one (co)sine constraint, the effect significantly amplifies. Thereby, for medium tolerances, PARA relaxations outperform SCIP stand-alone. Applied problems like alternating current optimal power flow (AC-OPF) feature such constraint types, leaving PARA a viable relaxation strategy. AB - Solving mixed-integer nonlinear programs (MINLPs) typically relies on constructing relaxations that are easier to tackle than the original problem. Recently, global parabolic (PARA) relaxations were introduced, featuring separable quadratic functions – paraboloids – as global under- or overestimators of general nonlinear constraint functions. So far, the paraboloids are all computed at once by solving a mixed-integer linear program (MIP). For small tolerances or wide function domains, the corresponding MIP grows in size and is eventually intractable, preventing a meaningful comparison with established relaxation techniques. We therefore propose a novel iterative method to compute PARA approximations that succeeds on all tolerance-domain combinations where the original one has failed. The computational study is preceded by a thorough theoretical explanation and analysis. Finally, the improved method enables a computational comparison with piecewise linear (PWL) relaxations in terms of runtime on general MINLP instances. The results show that the modern solver SCIP can solve PWL relaxations faster when the olerance is high, shifting strongly in favor of PARA for tighter tolerances. We attribute the effect to the difference in the corresponding problem size: PWL relaxations introduce binary variables to identify the active linear piece and their number grows with decreasing tolerance. PARA, on the other hand, does not require additional variables such that the dimension is maintained. For problems with at least one (co)sine constraint, the effect significantly amplifies. Thereby, for medium tolerances, PARA relaxations outperform SCIP stand-alone. Applied problems like alternating current optimal power flow (AC-OPF) feature such constraint types, leaving PARA a viable relaxation strategy. Y1 - 2026 ER - TY - INPR U1 - Preprint A1 - Klimm, Max A1 - Pfetsch, Marc E. A1 - Skutella, Martin A1 - Strubberg, Lea T1 - Approximating the Network Design Problem for Potential-Based Flows N2 - We develop efficient algorithms for a fundamental network design problem arising in potential-based flow models, which are central to many energy transport networks (e.g., hydrogen and electricity). In contrast to classical network flow problems, the nonlinearities inherent in potential-based networks introduce significant new challenges. We address these challenges through intricate reductions to classical combinatorial optimization problems, such as (constrained) shortest path problems, enabling the application of well-established algorithmic techniques to compute exact and approximate solutions efficiently. Finally, we complement these algorithmic results with matching complexity results concerning the hardness and non-approximability of the considered problem variants. AB - We develop efficient algorithms for a fundamental network design problem arising in potential-based flow models, which are central to many energy transport networks (e.g., hydrogen and electricity). In contrast to classical network flow problems, the nonlinearities inherent in potential-based networks introduce significant new challenges. We address these challenges through intricate reductions to classical combinatorial optimization problems, such as (constrained) shortest path problems, enabling the application of well-established algorithmic techniques to compute exact and approximate solutions efficiently. Finally, we complement these algorithmic results with matching complexity results concerning the hardness and non-approximability of the considered problem variants. Y1 - 2026 ER -