This paper is devoted to the construction of exponential integrators of first and second order for the time discretization of constrained parabolic systems. For this extend, we combine well-known exponential integrators for unconstrained systems with the solution of certain saddle point problems in order to meet the constraints throughout the integration process. The result is a novel class of semi-explicit time integration schemes. We prove the expected convergence rates and illustrate the performance on two numerical examples including a parabolic equation with nonlinear dynamic boundary conditions.

We consider partial differential equations on networks with a small parameter $\varepsilon$, which are hyperbolic for $\varepsilon>0$ and parabolic for $\varepsilon=0$. With a combination of an $\varepsilon$-expansion and Runge-Kutta schemes for constrained systems of parabolic type, we derive a new class of time discretization schemes for hyperbolic systems on networks, which are constrained due to interconnection conditions. For the analysis we consider the coupled system equations as partial differential-algebraic equations based on the variational formulation of the problem. We discuss well-posedness of the resulting systems and estimate the error caused by the $\varepsilon$-expansion.

Metriplectic systems are state space formulations that have become well-known under the acronym GENERIC. In this work we present a GENERIC based state space formulation in an operator setting that encodes a weak-formulation of the field equations describing the dynamics of a homogeneous mixture of compressible heat-conducting Newtonian fluids consisting of reactive constituents. We discuss the mathematical model of the fluid mixture formulated in the framework of continuum thermodynamics. The fluid mixture is considered an open thermodynamic system that moves free of external body forces. As closure relations we use the linear constitutive equations of the phenomenological theory known as Thermodynamics of Irreversible Processes (TIP). The phenomenological coefficients of these linear constitutive equations satisfy the Onsager-Casimir reciprocal relations. We present the state space representation of the fluid mixture, formulated in the extended GENERIC framework for open systems, specified by a symmetric, mixture related dissipation bracket and a mixture related Poisson-bracket for which we prove the Jacobi-identity.