Search

The modeling of electric machines and power transformers typically involves systems of nonlinear magnetostatics or -quasistatics, and their efficient and accurate simulation is required for the reliable design, control, and optimization of such devices. We study the numerical solution of the vector potential formulation of nonlinear magnetostatics by means of higher-order finite element methods. Numerical quadrature is used for the efficient handling of the nonlinearities and domain mappings are employed for the consideration of curved boundaries. The existence of a unique solution is proven on the continuous and discrete level and a full convergence analysis of the resulting finite element schemes is presented indicating order optimal convergence rates under appropriate smoothness assumptions. For the solution of the nonlinear discretized problems, we consider a Newton method with line search for which we establish global linear convergence with convergence rates that are independent of the discretization parameters. We further prove local quadratic convergence in a mesh-dependent neighborhood of the solution which becomes effective when high accuracy of the nonlinear solver is demanded. The assumptions required for our analysis cover inhomogeneous, nonlinear, and anisotropic materials, which may arise in typical applications, including the presence of permanent magnets. The theoretical results are illustrated by numerical tests for some typical benchmark problems.

We consider systematic numerical approximation of a viscoelastic phase separation model that describes the demixing of a polymer solvent mixture. An unconditionally stable discretisation method is proposed based on a finite element approximation in space and a variational time discretization strategy. The proposed method preserves the energy-dissipation structure of the underlying system exactly and allows to establish a fully discrete nonlinear stability estimate in natural norms based on the concept of relative energy. These estimates are used to derive order optimal error estimates for the method under minimal smoothness assumptions on the problem data, despite the presence of various strong nonlinearities in the equations. The theoretical results and main properties of the method are illustrated by numerical simulations which also demonstrate the capability to reproduce the relevant physical effects observed in experiments., Comment: 31 pages; 1 figure

The numerical solution of problems in nonlinear magnetostatics is typically based on a variational formulation in terms of magnetic potentials, the discretization by finite elements, and iterative solvers like the Newton method. The vector potential approach aims at minimizing a certain energy functional and, in three dimensions, requires the use of edge elements and appropriate gauging conditions. The scalar potential approach, on the other hand, seeks to maximize the negative coenergy and can be realized by standard Lagrange finite elements, thus reducing the number of degrees of freedom and simplifying the implementation. The number of Newton iterations required to solve the governing nonlinear system, however, has been observed to be usually higher than for the vector potential formulation. In this paper, we propose a method that combines the advantages of both approaches, i.e., it requires as few Newton iterations as the vector potential formulation while involving the magnetic scalar potential as the primary unknown. We discuss the variational background of the method, its well-posedness, and its efficient implementation. Numerical examples are presented for illustration of the accuracy and the gain in efficiency compared to other approaches.

A novel strategy is proposed for the coupling of field and circuit equations when modeling power devices in the low-frequency regime. The resulting systems of differential-algebraic equations have a particular geometric structure which explicitly encodes the energy storage, dissipation, and transfer mechanisms. This implies a power balance on the continuous level which can be preserved under appropriate discretization in space and time. The models and main results are presented in detail for linear constitutive models, but the extension to nonlinear elements and more general coupling mechanisms is possible. The theoretical findings are demonstrated by numerical results.

The magnetostatic field distribution in a nonlinear medium amounts to the unique minimizer of the magnetic coenergy over all fields that can be generated by the same current. This is a nonlinear saddlepoint problem whose numerical solution can in principle be achieved by mixed finite element methods and appropriate nonlinear solvers. The saddlepoint structure, however, makes the solution cumbersome. A remedy is to split the magnetic field into a known source field and the gradient of a scalar potential which is governed by a convex minimization problem. The penalty approach avoids the use of artificial potentials and Lagrange multipliers and leads to an unconstrained convex minimization problem involving a large parameter. We provide a rigorous justification of the penalty approach by deriving error estimates for the approximation due to penalization. We further highlight the close connections to the Lagrange-multiplier and scalar potential approach. The theoretical results are illustrated by numerical tests for a typical benchmark problem

We consider the convergence of iterative solvers for problems of nonlinear magnetostatics. Using the equivalence to an underlying minimization problem, we can establish global linear convergence of a large class of methods, including the damped Newton-method, fixed-point iteration, and the Kacanov iteration, which can all be interpreted as generalized gradient descent methods. Armijo backtracking isconsidered for an adaptive choice of the stepsize. The general assumptions required for our analysis cover inhomogeneous, nonlinear, and anisotropic materials, as well as permanent magnets. The main results are proven on the continuous level, but they carry over almost verbatim to various approximation schemes, including finite elements and isogeometric analysis, leading to bounds on the iteration numbers, which are independent of the particular discretization. The theoretical results are illustrated by numerical tests for a typical benchmark problem.

We consider systems of nonlinear magnetostatics and quasistatics that typically arise in the modeling and simulation of electric machines. The nonlinear problems, eventually obtained after time discretization, are usually solved by employing a vector potential formulation. In the relevant two-dimensional setting, a discretization can be obtained by H1-conforming finite elements. We here consider an alternative formulation based on the H-field which leads to a nonlinear saddlepoint problem. After commenting on the unique solvability, we study the numerical approximation by H(curl)-conforming finite elements and present the main convergence results. A particular focus is put on the efficient solution of the linearized systems arising in every step of the nonlinear Newton solver. Via hybridization, the linearized saddlepoint systems can be transformed into linear elliptic problems, which can be solved with similar computational complexity as those arising in the vector or scalar potential formulation. In summary, we can thus claim that the mixed finite element approach based on the $H$-field can be considered a competitive alternative to the standard vector or scalar potential formulations for the solution of problems in nonlinear magneto-quasistatics.

Implicit models for magnetic coenergy have been proposed by Pera et al. to describe the anisotropic nonlinear material behavior of electrical steel sheets. This approach aims at predicting magnetic response for any direction of excitation by interpolating measured of B--H curves in the rolling and transverse directions. In an analogous manner, an implicit model for magnetic energy is proposed. We highlight some mathematical properties of these implicit models and discuss their numerical realization, outline the computation of magnetic material laws via implicit differentiation, and discuss the potential use for finite element analysis in the context of nonlinear magnetostatics.

We study the numerical solution of a Cahn-Hilliard/Allen-Cahn system with strong coupling through state and gradient dependent non-diagonal mobility matrices. A fully discrete approximation scheme in space and time is proposed which preserves the underlying gradient flow structure and leads to dissipation of the free-energy on the discrete level. Existence and uniqueness of the discrete solution is established and relative energy estimates are used to prove optimal convergence rates in space and time under minimal smoothness assumptions. Numerical tests are presented for illustration of the theoretical results and to demonstrate the viability of the proposed methods., Comment: 23 pages, 14 figures, 1 table. arXiv admin note: text overlap with arXiv:2209.03849

A novel finite element scheme is studied for solving the time-dependent Maxwell's equations on unstructured grids efficiently. Similar to the traditional Yee scheme, the method has one degree of freedom for most edges and a sparse inverse mass matrix. This allows for an efficient realization by explicit time-stepping without solving linear systems. The method is constructed by algebraic reduction of another underlying finite element scheme which involves two degrees of freedom for every edge. Mass-lumping and additional modifications are used in the construction of this method to allow for the mentioned algebraic reduction in the presence of source terms and lossy media later on. A full error analysis of the underlying method is developed which by construction also carries over to the reduced scheme and allows to prove convergence rates for the latter. The efficiency and accuracy of both methods are illustrated by numerical tests. The proposed schemes and their analysis can be extended to structured grids and in special cases the reduced method turns out to be algebraically equivalent to the Yee scheme. The analysis of this paper highlights possible difficulties in extensions of the Yee scheme to non-orthogonal or unstructured grids, discontinuous material parameters, and non-smooth source terms, and also offers potential remedies.

We study a doubly nonlinear parabolic problem arising in the modeling of gas transport in pipelines. Using convexity arguments and relative entropy estimates we show uniform bounds and exponential stability of discrete approximations obtained by a finite element method and implicit time stepping. Due to convergence of the approximations to weak solutions of the problem, our results also imply regularity, uniqueness, and long time stability of weak solutions of the continuous problem.

The propagation of charged particles through a scattering medium in the presence of a magnetic field can be described by a Fokker-Planck equation with Lorentz force. This model is studied both, from a theoretical and a numerical point of view. A particular trace estimate is derived for the relevant function spaces to clarify the meaning of boundary values. Existence of a weak solution is then proven by the Rothe method. In the second step of our investigations, a fully practicable discretization scheme is proposed based on implicit time-stepping through the energy levels and a spherical-harmonics finite-element discretization with respect to the remaining variables. A full error analysis of the resulting scheme is given, and numerical results are presented to illustrate the theoretical results and the performance of the proposed method.

A system of phase-field equations with strong-coupling through state and gradient dependent non-diagonal mobility matrices is studied. Existence of weak solutions is established by the Galerkin approximation and a-priori estimates in strong norms. Relative energy estimates are used to derive a general nonlinear stability estimate. As a consequence, a weak-strong uniqueness principle is obtained and stability with respect to model parameters is investigated., Comment: 15 pages, 1 figure

In this work, we present a variational and quantitative phase-field model for non-isothermal sintering processes. The model is derived via an extended non-diagonal phase-field model. The model evolution equations have naturally cross-coupling terms between the conserved kinetics (i.e., mass and thermal transfer) and the non-conserved one (grain growth). These terms are shown via asymptotic analysis to be instrumental in ensuring the elimination of interface artefacts, while also examined to not modify the thermodynamic equilibrium condition (characterized by dihedral angle). Moreover, we demonstrate that the trapping effects and existence of surface diffusion in conservation laws are direction-dependent. An anisotropic interpolation scheme of the kinetic mobilities which differentiates the normal and the tangential directions along the interface is discussed. Numerically, we demonstrate the importance of the cross-couplings and the anisotropic interpolation via presenting thermal-microstructural evolutions.

A novel finite element method for the approximation of Maxwell's equations over hybrid two-dimensional grids is studied. The choice of appropriate basis functions and numerical quadrature leads to diagonal mass matrices which allow for efficient time integration by explicit methods. On purely rectangular grids, the proposed schemes coincide with well-established FIT and FDTD methods. Additional internal degrees of freedom introduced on triangles allow for mass-lumping without the usual constraints on the shape of these elements. A full error analysis of the method is developed and numerical tests are presented for illustration.

We study the numerical approximation of singularly perturbed convection-diffusion problems on one-dimensional pipe networks. In the vanishing diffusion limit, the number and type of boundary conditions and coupling conditions at network junctions changes, which gives rise to singular layers at the outflow boundaries of the pipes. A hybrid discontinuous Galerkin method is proposed, which provides a natural upwind mechanism for the convection-dominated case. Moreover, the method automatically handles the variable coupling and boundary conditions in the vanishing diffusion limit, leading to an asymptotic-preserving scheme. A detailed analysis of the singularities of the solution and the discretization error is presented, and an adaptive strategy is proposed, leading to order optimal error estimates that hold uniformly in the singular perturbation limit. The theoretical results are confirmed by numerical tests.

We propose and analyze a structure-preserving space-time variational discretization method for the Cahn-Hilliard-Navier-Stokes system. Uniqueness and stability for the discrete problem is established in the presence of concentration dependent mobility and viscosity parameters by means of the relative energy estimates and order optimal convergence rates are established for all variables using balanced approximation spaces and relaxed regularity conditions on the solution. Numerical tests are presented to demonstrate the proposed method is fully practical and yields the predicted convergence rates. The discrete stability estimates developed in this paper may also be used to analyse other discretization schemes, which is briefly outlined in the discussion., Comment: 18 pages + 14 appendix; 2 figures, 1 table

We consider the identifiability and stable numerical estimation of multiple parameters in a Cahn-Hilliard model for phase separation. Spatially resolved measurements of the phase fraction are assumed to be accessible, with which the identifiability of single and multiple parameters up to certain scaling invariances is established. A regularized equation error approach is proposed for the stable numerical solution of the parameter identification problems, and convergence of the regularized approximations is proven under reasonable assumptions on the data noise. The viability of the theoretical results and the proposed methods is demonstrated in numerical tests.

We present two strategies for designing passivity preserving higher order discretization methods for Maxwell's equations in nonlinear Kerr-type media. Both approaches are based on variational approximation schemes in space and time. This allows to rigorously prove energy conservation or dissipation, and thus passivity, on the fully discrete level. For linear media, the proposed methods coincide with certain combinations of mixed finite element and implicit Runge-Kutta schemes. The order optimal convergence rates, which can thus be expected for linear problems, are also observed for nonlinear problems in the numerical tests.

We consider the dynamics of an elastic continuum under large deformation but small strain. Such systems can be described by the equations of geometrically nonlinear elastodynamics in combination with the St. Venant-Kirchhoff material law. The velocity-stress formulation of the problem turns out to have a formal port-Hamiltonian structure. In contrast to the linear case, the operators of the problem are modulated by the displacement field which can be handled as a passive variable and integrated along with the velocities. A weak formulation of the problem is derived and essential boundary conditions are incorporated via Lagrange multipliers. This variational formulation explicitly encodes the transfer between kinetic and potential energy in the interior as well as across the boundary, thus leading to a global power balance and ensuring passivity of the system. The particular geometric structure of the weak formulation can be preserved under Galerkin approximation via appropriate mixed finite elements. In addition, a fully discrete power balance can be obtained by appropriate time discretization. The main properties of the system and its discretization are shown theoretically and demonstrated by numerical tests., Comment: The new version of the manuscript is a result of discussions with our new author. In the first version, we focused on the generation of explicit state space models, our attention is now on the consistent structure-preserving discretization in space and time for geometrically nonlinear elastodynamics. (15 pages, submitted to Mathematical and Computer Modelling of Dynamical Systems)

The use of trigonometric polynomials as Lagrange multipliers in the harmonic mortar method enables an efficient and elegant treatment of relative motion in the stator-rotor coupling of electric machine simulation. Explicit formulas for the torque computation are derived by energetic considerations, and their realization by harmonic mortar finite element and isogemetric analysis discretizations is discussed. Numerical tests are presented to illustrate the theoretical results and demonstrate the potential of harmonic mortar methods for the evaluation of torque ripples.

The modified nodal analysis (MNA) is probably the most widely used formulation for the modeling and simulation of electric circuits. Its conventional form uses electric node potentials and currents across inductors and voltage sources as unknowns, thus taking an electric viewpoint. In this paper, we propose a magnetic oriented nodal analysis (MONA) for electric circuits, which is based on magnetic node potentials and charges across capacitors and voltage sources as the primary degrees of freedom, thus giving direct access to these quantities. The resulting system has the structure of a generalized gradient system which immediately ensures passivity in the absence of sources. A complete index analysis is presented showing regularity of the magnetic oriented formulation under standard topological conditions on the network interconnection. In comparison to conventional MNA, the differential-algebraic index is reduced by one in most cases which facilitates the numerical solution. Some preliminary numerical experiments are presented for illustration of the feasibility and stability of the new approach.

The numerical solution of dynamical systems with memory requires the efficient evaluation of Volterra integral operators in an evolutionary manner. After appropriate discretisation, the basic problem can be represented as a matrix-vector product with a lower diagonal but densely populated matrix. For typical applications, like fractional diffusion or large scale dynamical systems with delay, the memory cost for storing the matrix approximations and complete history of the data then would become prohibitive for an accurate numerical approximation. For Volterra-integral operators of convolution type, the \emph{fast and oblivious convolution quadrature} method of Sch\"adle, Lopez-Fernandez, and Lubich allows to compute the discretized valuation with $N$ time steps in $O(N \log N)$ complexity and only requiring $O(\log N)$ active memory to store a compressed version of the complete history of the data. We will show that this algorithm can be interpreted as an $\mathcal{H}$-matrix approximation of the underlying integral operator and, consequently, a further improvement can be achieved, in principle, by resorting to $\mathcal{H}^2$-matrix compression techniques. We formulate a variant of the $\mathcal{H}^2$-matrix vector product for discretized Volterra integral operators that can be performed in an evolutionary and oblivious manner and requires only $O(N)$ operations and $O(\log N)$ active memory. In addition to the acceleration, more general asymptotically smooth kernels can be treated and the algorithm does not require a-priori knowledge of the number of time steps. The efficiency of the proposed method is demonstrated by application to some typical test problems.

Based on relative energy estimates, we study the stability of solutions to the Cahn-Hilliard equation with concentration dependent mobility with respect to perturbations. As a by-product of our analysis, we obtain a weak-strong uniqueness principle on the continuous level under realistic regularity assumptions on strong solutions. We then show that the stability estimates can be further inherited almost verbatim by appropriate Galerkin approximations in space and time. This allows us to derive sharp bounds for the discretization error in terms of certain projection errors and to establish order-optimal a-priori error estimates for semi- and fully discrete approximation schemes., Comment: 27 pages, 6 figures

(abridged) We present a detailed derivation of a simple hydrodynamic two-fluid model, which aims at the description of the phase separation of non-entangled polymer solutions, where viscoelastic effects play a role. It is directly based upon the coarse-graining of a well-defined molecular model, such that all degrees of freedom have a clear and unambiguous molecular interpretation. The considerations are based upon a free-energy functional, and the dynamics is split into a conservative and a dissipative part, where the latter satisfies the Onsager relations and the Second Law of thermodynamics. The model is therefore fully consistent with both equilibrium and non-equilibrium thermodynamics. The derivation proceeds in two steps: Firstly, we derive an extended model comprising two scalar and four vector fields, such that inertial dynamics of the macromolecules and of the relative motion of the two fluids is taken into account. In the second step, we eliminate these inertial contributions and, as a replacement, introduce phenomenological dissipative terms, which can be modeled easily by taking into account the principles of non-equilibrium thermodynamics. The final simplified model comprises the momentum conservation equation, which includes both interfacial and elastic stresses, a convection-diffusion equation where interfacial and elastic contributions occur as well, and a suitably convected relaxation equation for the end-to-end vector field. In contrast to the traditional two-scale description that is used to derive rheological equations of motion, we here treat the hydrodynamic and the macromolecular degrees of freedom on the same basis. Nevertheless, the resulting model is fairly similar, though not fully identical, to models that have been discussed previously. Notably, we find a rheological constitutive equation that differs from the standard Oldroyd-B model. ..., Comment: Submitted to Journal of Physics Condensed Matter

We consider the transport of gas in long pipes and pipeline networks for which the dynamics are dominated by friction at the pipe walls. The governing equations can be formulated as an abstract dissipative Hamiltonian system which allows us to derive perturbation bounds by means of relative energy estimates. As particular consequences, we obtain stability with respect to initial conditions and model parameters and quantitative estimates in the high friction limit. Our results are established in detail for the flow in a single pipe and through the energy-based modelling they naturally generalize also to pipe networks.

A new model for viscoelastic phase separation is proposed, based on a systematically derived conservative two-fluid model. Dissipative effects are included by phenomenological viscoelastic terms. By construction, the model is consistent with the second law of thermodynamics, and we study well-posedness of the model, i.e., existence of weak solutions, a weak-strong uniqueness principle, and stability with respect to perturbations, which are proven by means of relative energy estimates. A good qualitative agreement with mesoscopic simulations is observed in numerical tests., Comment: 10 pages, 8 figures

Wall shear stress (WSS) has been suggested as a potential biomarker in various cardiovascular diseases and it can be estimated from phase-contrast Magnetic Resonance Imaging (PC-MRI) velocity measurements. We present a new parametric sequential method for MRI-based WSS quantification consisting of a geometry identification and a subsequent approximation of the velocity field. The new method was validated on well controlled, high resolution in vitro measurements of turbulent stationary flows and physiological pulsatile flows in phantoms. Basic applicability to in vivo was demonstrated using 2D PC-MRI data of the ascending aorta of three volunteers.

Harmonic stator-rotor coupling offers a promising approach for the interconnection of rotating subsystems in the simulation of electric machines. This paper studies the stability of discretization schemes based on harmonic coupling in the framework of mortar methods for Poisson-like problems. A general criterion is derived that allows to ensure the relevant inf-sup stability condition for a variety of specific discretization approaches, including finite-element methods and isogeometric analysis with harmonic mortar coupling. The validity and sharpness of the theoretical results is demonstrated by numerical tests.

We study the efficient numerical solution of linear inverse problems with operator valued data which arise, e.g., in seismic exploration, inverse scattering, or tomographic imaging. The high-dimensionality of the data space implies extremely high computational cost already for the evaluation of the forward operator which makes a numerical solution of the inverse problem, e.g., by iterative regularization methods, practically infeasible. To overcome this obstacle, we take advantage of the underlying tensor product structure of the problem and propose a strategy for constructing low-dimensional certified reduced order models of quasi-optimal rank for the forward operator which can be computed much more efficiently than the truncated singular value decomposition. A complete analysis of the proposed model reduction approach is given in a functional analytic setting and the efficient numerical construction of the reduced order models as well as of their application for the numerical solution of the inverse problem is discussed. In summary, the setup of a low-rank approximation can be achieved in an offline stage at essentially the same cost as a single evaluation of the forward operator, while the actual solution of the inverse problem in the online phase can be done with extremely high efficiency. The theoretical results are illustrated by application to a typical model problem in fluorescence optical tomography.

We consider singularly perturbed convection-diffusion equations on one-dimensional networks (metric graphs) as well as the transport problems arising in the vanishing diffusion limit. Suitable coupling condition at inner vertices are derived that guarantee conservation of mass as well as dissipation of a mathematical energy which allows us to prove stability and well-posedness. For single intervals and appropriately specified initial conditions, it is well-known that the solutions of the convection-diffusion problem converge to that of the transport problem with order $O(\sqrt{\epsilon})$ in the $L^\infty(L^2)$-norm with diffusion $\epsilon \to 0$. In this paper, we prove a corresponding result for problems on one-dimensional networks. The main difficulty in the analysis is that the number and type of coupling conditions changes in the singular limit which gives rise to additional boundary layers at the interior vertices of the network. Since the values of the solution at these network junctions are not known a-priori, the asymptotic analysis requires a delicate choice of boundary layer functions that allows to handle these interior layers.

We study the systematic numerical approximation of Maxwell's equations in dispersive media. Two discretization strategies are considered, one based on a traditional leapfrog time integration method and the other based on convolution quadrature. The two schemes are proven to be equivalent and to preserve the underlying energy-dissipation structure of the problem. The second approach, however, is independent of the number of internal states and allows to handle rather general dispersive materials. Using ideas of fast-and-oblivious convolution quadrature, the method can be implemented efficiently.

We consider the numerical approximation of Maxwell's equations in time domain by a second order $H(curl)$ conforming finite element approximation. In order to enable the efficient application of explicit time stepping schemes, we utilize a mass-lumping strategy resulting from numerical integration in conjunction with the finite element spaces introduced by Elmkies and Joly. We prove that this method is second order accurate if the true solution is divergence free, but the order of accuracy reduces to one in the general case. We then propose a modification of the finite element space, which yields second order accuracy in the general case.

We discuss the mathematical modeling and numerical discretization of transport problems on one-dimensional networks. Suitable coupling conditions are derived that guarantee conservation of mass across network junctions and dissipation of a mathematical energy which allows to prove existence of unique solutions. We then consider the space discretization by a hybrid discontinuous Galerkin method which provides a suitable upwind mechanism to handle the transport problem and allows to incorporate the coupling conditions in a natural manner. In addition, the method inherits mass conservation and stability of the continuous problem. Order optimal convergence rates are established and illustrated by numerical tests.

We consider the systematic numerical approximation of Biot's quasistatic model for the consolidation of a poroelastic medium. Various discretization schemes have been analysed for this problem and inf-sup stable finite elements have been found suitable to avoid spurios pressure oscillations in the initial phase of the evolution. In this paper, we first clarify the role of the inf-sup condition for the well-posedness of the continuous problem and discuss the choice of appropriate initial conditions. We then develop an abstract error analysis that allows us to analyse some approximation schemes discussed in the literature in a unified manner. In addition, we propose and analyse the high-order time discretization by a scheme that can be interpreted as a variant of continuous-Galerkin or particular Runge-Kutta methods applied to a modified system. The scheme is designed to preserve both, the underlying differential-algebraic structure and energy-dissipation property of the problem. In summary, we obtain high-order Galerkin approximations with respect to space and time and derive order-optimal convergence rates. The numerical analysis is carried out in detail for the discretization of the two-field formulation by Taylor-Hood elements and a variant of a Runge-Kutta time discretization. Our arguments can however be extended to three- and four field formulations and other time discretization strategies.

We consider the efficient numerical approximation of acoustic wave propagation in time domain by a finite element method with mass lumping. In the presence of internal damping, the problem can be reduced to a second order formulation in time for the velocity field alone. For the spatial approximation we consider $H(\mathrm{div})$--conforming finite elements of second order. In order to allow for an efficient time integration, we propose a mass-lumping strategy based on approximation of the $L^2$-scalar product by inexact numerical integration which leads to a block-diagonal mass matrix. A careful error analysis allows to show that second order accuracy is not reduced by the quadrature errors which is illustrated also by numerical tests.

This paper addresses the numerical solution of the Westervelt equation, which arises as one of the model equations in nonlinear acoustics. The problem is rewritten in a canonical form that allows the systematic discretization by Galerkin approximation in space and time. Exact energy preserving methods of formally arbitrary order are obtained and their efficient realization as well as the relation to other frequently used methods is discussed.

We consider the stable reconstruction of flow geometry and wall shear stress from measurements obtained by magnetic resonance imaging. As noted in a review article by Petersson, most approaches considered so far in the literature seem not be satisfactory. We therefore propose a systematic reconstruction procedure that allows to obtain stable estimates of flow geometry and wall shear stress and we are able to quantify the reconstruction errors in terms of bounds for the measurement errors under reasonable smoothness assumptions. A full analysis of the approach is given in the framework of regularization methods. In addition, we discuss the efficient implementation of our method and we demonstrate its viability, accuracy, and regularizing properties for experimental data.

A general framework for the numerical approximation of evolution problems is presented that allows to preserve exactly an underlying Hamiltonian- or gradient structure. The approach relies on rewriting the evolution problem in a particular form that complies with the underlying geometric structure. The Galerkin approximation of a corresponding variational formulation in space then automatically preserves this structure which allows to deduce important properties for appropriate discretization schemes including projection based model order reduction. We further show that the underlying structure is preserved also under time discretization by a Petrov-Galerkin approach. The presented framework is rather general and allows the numerical approximation of a wide range of applications, including nonlinear partial differential equations and port-Hamiltonian systems. Some examples will be discussed for illustration of our theoretical results and connections to other discretization approaches will be revealed.

We consider the numerical approximation of single phase flow in porous media by a mixed finite element method with mass lumping. Our work extends previous results of Wheeler and Yotov, who showed that mass lumping together with an appropriate choice of basis allows to eliminate the flux variables locally and to reduced the mixed problem in this way to a finite volume discretization for the pressure only. Here we construct second order approximations for hybrid meshes in two and three space dimensions which, similar to the method of Wheeler and Yotov, allows the local elimination of the flux variables. A full convergence analysis of the method is given for which new arguments and, in part, also new quadrature rules and finite elements are required. Computational tests are presented for illustration of the theoretical results.

We consider the discretization of a semilinear damped wave equation arising, for instance, in the modeling of gas transport in pipeline networks. For time invariant boundary data, the solutions of the problem are shown to converge exponentially fast to steady states. We further prove that this decay behavior is inherited uniformly by a class of Galerkin approximations, including finite element, spectral, and structure preserving model reduction methods. These theoretical findings are illustrated by numerical tests.

We consider the efficient numerical solution of coupled dynamical systems, consisting of a small nonlinear part and a large linear time invariant part, possibly stemming from spatial discretization of an underlying partial differential equation. The linear subsystem can be eliminated in frequency domain and for the numerical solution of the resulting integro-differential algebraic equations, we propose a a combination of Runge-Kutta or multistep time stepping methods with appropriate convolution quadrature to handle the integral terms. The resulting methods are shown to be algebraically equivalent to a Runge-Kutta or multistep solution of the coupled system and thus automatically inherit the corresponding stability and accuracy properties. After a computationally expensive pre-processing step, the online simulation can, however, be performed at essentially the same cost as solving only the small nonlinear subsystem. The proposed method is, therefore, particularly attractive, if repeated simulation of the coupled dynamical system is required.

This paper studies the discretization of gas transport in pipeline networks by an inexact Petrov-Galerkin method. A full convergence analysis is presented for single pipes under the assumption of a linear friction law and the possible extension to pipe networks is discussed. The generalization to nonlinear gas transport models and the efficient implementation by hybridization is investigated numerically.

A novel mass-lumping strategy for a mixed finite element approximation of Maxwell's equations is proposed. On structured orthogonal grids the resulting method coincides with the spatial discretization of the Yee scheme. The proposed method, however, generalizes naturally to unstructured grids and anisotropic materials and thus yields a variational extension of the Yee scheme for these situations.