C04
Refine
Language
- English (19)
Has Fulltext
- yes (19)
Keywords
- First-Discretize-Then-Optimize (1)
- Galerkin projection (1)
- Mixed-Integer Nonlinear Programming (1)
- Power-to-Gas (1)
- Storage Capacity Maximization (1)
- Transient Gas Transport Optimization (1)
- asymptotic limits (1)
- gas transport on networks (1)
- hyperbolic balance laws (1)
- partial differential-algebraic equations (1)
We consider the identification of nonlinear diffusion coefficients of the form a(t,u) or a(u) in quasi-linear parabolic and elliptic equations. Uniqueness for this inverse problem is established under very general assumptions using partial knowledge of the Dirichlet-to-Neumann map. The proof of our main result relies on the construction of a series of appropriate Dirichlet data and test functions with a particular singular behavior at the boundary. This allows us to localize the analysis and to separate the principal part of the equation from the remaining terms. We therefore do not require specific knowledge of lower order terms or initial data which allows to apply our results to a variety of applications. This is illustrated by discussing some typical examples in detail.
Novel experimental modalities acquire spatially resolved velocity measurements for steady state and transient flows which are of interest for engineering and biological applications. One of the drawbacks of such high resolution velocity data is their susceptibility to measurement errors. In this paper, we propose a novel filtering strategy that allows enhancement of noisy measurements to obtain reconstruction of smooth divergence free velocity and corresponding pressure fields, which together approximately comply to a prescribed flow model. The main step in our approach consists of the appropriate use of the velocity measurements in the design of a linearized flow model which can be shown to be well-posed and consistent with the true velocity and pressure fields up to measurement and modeling errors. The reconstruction procedure is formulated
as a linear quadratic optimal control problem and the resulting filter has analyzable smoothing and approximation properties. We also discuss briefly the discretization of our approach by finite element methods and comment on the efficient solution of the linear optimality system by iterative solvers. The capability of the proposed method to significantly reduce data noise is demonstrated by numerical tests in which we also compare to other methods like smoothing and solenoidal filtering.
We consider a parameter dependent family of damped hyperbolic equations with interesting limit behavior: the system approaches steady states exponentially fast and for parameter to zero the solutions converge to that of a parabolic limit problem. We establish sharp estimates and elaborate their dependence on the model parameters. For the numerical approximation we then consider a mixed finite element method in space together with a Runge-Kutta method in time. Due to the variational and dissipative nature of this approximation, the limit behavior of the infinite dimensional level is inherited almost automatically by the discrete problems. The resulting numerical method thus is asymptotic preserving in the parabolic limit and uniformly exponentially stable. These results are further shown to be independent of the discretization parameters. Numerical tests are presented for a simple model problem which illustrate that the derived estimates are sharp in general.
The numerical solution of time-dependent radiative transfer problems is
challenging, both, due to the high dimension as well as the anisotropic structure of the underlying integro-partial differential equation. In this paper we propose a general framework for designing numerical methods for time-dependent radiative transfer based on a Galerkin discretization in space and angle combined with appropriate time stepping schemes. This allows us to systematically incorporate boundary conditions and to preserve basic properties like exponential stability and decay to equilibrium also on the discrete level. We present the basic a-priori error analysis and provide abstract error estimates that cover a wide class of methods. The starting point for our considerations is to rewrite the radiative transfer problem as a system of evolution equations which has a similar structure like first order hyperbolic systems in acoustics or electrodynamics. This analogy allows us to generalize the main arguments of the numerical analysis for such applications to the radiative transfer problem under investigation. We also discuss a particular discretization scheme based on a truncated spherical harmonic expansion in angle, a finite element discretization in space, and the implicit Euler method in time. The performance of the resulting mixed PN-finite element time stepping scheme
is demonstrated by computational results.
We consider the discretization and subsequent model reduction of a system of partial differential-algebraic equations describing the propagation of pressure waves in a pipeline network. Important properties like conservation of mass, dissipation of energy, passivity, existence of steady states, and exponential stability can be preserved by an appropriate semi-
discretization in space via a mixed finite element method and also during the further dimension reduction by structure preserving Galerkin projection which is the main focus of this paper. Krylov subspace methods are employed for the construction of the reduced models and we discuss modifications needed to satisfy certain algebraic compatibility conditions; these are required to ensure the well-posedness of the reduced models and the preservation of the key properties. Our analysis is based on the underlying infinite dimensional problem and its Galerkin approximations. The proposed algorithms therefore have a direct interpretation in function spaces; in principle, they are even applicable directly to the original system of partial differential-algebraic
equations while the intermediate discretization by finite elements is only required for the actual
computations. The performance of the proposed methods is illustrated with numerical tests and the necessity for the compatibility conditions is demonstrated by examples.
We consider the identification of a nonlinear friction law in a one-dimensional
damped wave equation from additional boundary measurements. Well-posedness of the
governing semilinear hyperbolic system is established via semigroup theory and con-
traction arguments. We then investigte the inverse problem of recovering the unknown
nonlinear damping law from additional boundary measurements of the pressure drop
along the pipe. This coefficient inverse problem is shown to be ill-posed and a varia-
tional regularization method is considered for its stable solution. We prove existence of
minimizers for the Tikhonov functional and discuss the convergence of the regularized so-
lutions under an approximate source condition. The meaning of this condition and some
arguments for its validity are discussed in detail and numerical results are presented for
illustration of the theoretical findings
We consider the numerical approximation of compressible flow in a pipe net-
work. Appropriate coupling conditions are formulated that allow us to derive a variational
characterization of solutions and to prove global balance laws for the conservation of mass
and energy on the whole network. This variational principle, which is the basis of our fur-
ther investigations, is amenable to a conforming Galerkin approximation by mixed finite
elements. The resulting semi-discrete problems are well-posed and automatically inherit the
global conservation laws for mass and energy from the continuous level. We also consider the
subsequent discretization in time by a problem adapted implicit time stepping scheme which
leads to conservation of mass and a slight dissipation of energy of the full discretization.
The well-posedness of the fully discrete scheme is established and a fixed-point iteration is
proposed for the solution of the nonlinear systems arising in every single time step. Some
computational results are presented for illustration of our theoretical findings and for demon-
stration of the robustness and accuracy of the new method
We consider the numerical approximation of linear damped wave systems
by Galerkin approximations in space and appropriate time-stepping schemes. Based on
a dissipation estimate for a modified energy, we prove exponential decay of the physical
energy on the continuous level provided that the damping is effective everywhere in the
domain. The methods of proof allow us to analyze also a class of Galerkin approximations
based on a mixed variational formulation of the problem. Uniform exponential stabil-
ity can be guaranteed for these approximations under a general compatibility condition
on the discretization spaces. As a particular example, we discuss the discretization by
mixed finite element methods for which we obtain convergence and uniform error esti-
mates under minimal regularity assumptions. We also prove unconditional and uniform
exponential stability for the time discretization by certain one-step methods. The valid-
ity of the theoretical results as well as the necessity of some of the conditions required
for our analysis are demonstrated in numerical tests
We consider a damped linear hyperbolic system modelling the propagation
of pressure waves in a network of pipes. Well-posedness is established via semi-group
theory and the existence of a unique steady state is proven in the absence of driving
forces. Under mild assumptions on the network topology and the model parameters,
we show exponential stability and convergence to equilibrium. This generalizes related
results for single pipes and multi-dimensional domains to the network context. Our proof
of the exponential stability estimate is based on a variational formulation of the problem,
some graph theoretic results, and appropriate energy estimates. The main arguments
are rather generic and can be applied also for the analysis of Galerkin approximations.
Uniform exponential stability can be guaranteed for the resulting semi-discretizations
under mild compatibility conditions on the approximation spaces. A particular realiza-
tion by mixed finite elements is discussed and the theoretical results are illustrated by
numerical tests in which also bounds for the decay rate are investigated.