BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Ulrich Langer

Ulrich Langer appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NANumerical Analysismath.OCastro-ph.IM

Catalog footprint

What is connected

26works
4topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

26 published item(s)

preprint2022arXiv

Matrix-free Monolithic Multigrid Methods for Stokes and Generalized Stokes Problems

We consider the widely used continuous $\mathcal{Q}_{k}$-$\mathcal{Q}_{k-1}$ quadrilateral or hexahedral Taylor-Hood elements for the finite element discretization of the Stokes and generalized Stokes systems in two and three spatial dimensions. For the fast solution of the corresponding symmetric, but indefinite system of finite element equations, we propose and analyze matrix-free monolithic geometric multigrid solvers that are based on appropriately scaled Chebyshev-Jacobi smoothers. The analysis is based on results by Schöberl and Zulehner (2003). We present and discuss several numerical results for typical benchmark problems.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Robust finite element discretization and solvers for distributed elliptic optimal control problems

We consider standard tracking-type, distributed elliptic optimal control problems with $L^2$ regularization, and their finite element discretization. We are investigating the $L^2$ error between the finite element approximation $u_{\varrho h}$ of the state $u_\varrho$ and the desired state (target) $\bar{u}$ in terms of the regularization parameter $\varrho$ and the mesh size $h$ that leads to the optimal choice $\varrho = h^4$. It turns out that, for this choice of the regularization parameter, we can devise simple Jacobi-like preconditioned MINRES or Bramble-Pasciak CG methods that allow us to solve the reduced discrete optimality system in asymptotically optimal complexity with respect to the arithmetical operations and memory demand. The theoretical results are confirmed by several benchmark problems with targets of various regularities including discontinuous targets.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Robust discretization and solvers for elliptic optimal control problems with energy regularization

We consider the finite element discretization and the iterative solution of singularly perturbed elliptic reaction-diffusion equations in three-dimensional computational domains. These equations arise from the optimality conditions for elliptic distributed optimal control problems with energy regularization that were recently studied by M.~Neumüller and O.~Steinbach (2020). We provide quasi-optimal a priori finite element error estimates which depend both on the mesh size $h$ and on the regularization parameter $\varrho$. The choice $\varrho = h^2$ ensures optimal convergence which only depends on the regularity of the target function. For the iterative solution, we employ an algebraic multigrid preconditioner and a balancing domain decomposition by constraints (BDDC) preconditioner. We numerically study robustness and efficiency of the proposed algebraic preconditioners with respect to the mesh size $h$, the regularization parameter $\varrho$, and the number of subdomains (cores) $p$. Furthermore, we investigate the parallel performance of the BDDC preconditioned conjugate gradient solver.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Adaptive space-time finite element methods for non-autonomous parabolic problems with distributional sources

We consider locally stabilized, conforming finite element schemes on completely unstructured simplicial space-time meshes for the numerical solution of parabolic initial-boundary value problems with variable, possibly discontinuous in space and time, coefficients. Distributional sources are also admitted. Discontinuous coefficients, non-smooth boundaries, changing boundary conditions, non-smooth or incompatible initial conditions, and non-smooth right-hand sides can lead to non-smooth solutions. We present new a priori and a posteriori error estimates for low-regularity solutions. In order to avoid reduced convergence rates appearing in the case of uniform mesh refinement, we also consider adaptive refinement procedures based on residual a posteriori error indicators and functional a posteriori error estimators. The huge system of space-time finite element equations is then solved by means of GMRES preconditioned by space-time algebraic multigrid. In particular, in the 4d space-time case that is 3d in space, simultaneous space-time parallelization can considerably reduce the computational time. We present and discuss numerical results for several examples possessing different regularity features.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Efficient Direct Space-Time Finite Element Solvers for Parabolic Initial-Boundary Value Problems in Anisotropic Sobolev Spaces

We consider a space-time variational formulation of parabolic initial-boundary value problems in anisotropic Sobolev spaces in combination with a Hilbert-type transformation. This variational setting is the starting point for the space-time Galerkin finite element discretization that leads to a large global linear system of algebraic equations. We propose and investigate new efficient direct solvers for this system. In particular, we use a tensor-product approach with piecewise polynomial, globally continuous ansatz and test functions. The developed solvers are based on the Bartels-Stewart method and on the Fast Diagonalization method, which result in solving a sequence of spatial subproblems. The solver based on the Fast Diagonalization method allows to solve these spatial subproblems in parallel leading to a full parallelization in time. We analyze the complexity of the proposed algorithms, and give numerical examples for a two-dimensional spatial domain, where sparse direct solvers for the spatial subproblems are used.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Parallel matrix-free higher-order finite element solvers for phase-field fracture problems

Phase-field fracture models lead to variational problems that can be written as a coupled variational equality and inequality system. Numerically, such problems can be treated with Galerkin finite elements and primal-dual active set methods. Specifically, low-order and high-order finite elements may be employed, where, for the latter, only few studies exist to date. The most time-consuming part in the discrete version of the primal-dual active set (semi-smooth Newton) algorithm consists in the solutions of changing linear systems arising at each semi-smooth Newton step. We propose a new parallel matrix-free monolithic multigrid preconditioner for these systems. We provide two numerical tests, and discuss the performance of the parallel solver proposed in the paper. Furthermore, we compare our new preconditioner with a block-AMG preconditioner available in the literature.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Reliability and efficiency of DWR-type a posteriori error estimates with smart sensitivity weight recovering

We derive efficient and reliable goal-oriented error estimations, and devise adaptive mesh procedures for the finite element method that are based on the localization of a posteriori estimates. In our previous work [SIAM J. Sci. Comput., 42(1), A371--A394, 2020], we showed efficiency and reliability for error estimators based on enriched finite element spaces. However, the solution of problems on a enriched finite element space is expensive. In the literature, it is well known that one can use some higher-order interpolation to overcome this bottleneck. Using a saturation assumption, we extend the proofs of efficiency and reliability to such higher-order interpolations. The results can be used to create a new family of algorithms, where one of them is tested on three numerical examples (Poisson problem, p-Laplace equation, Navier-Stokes benschmark), and is compared to our previous algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Space-time finite element discretization of parabolic optimal control problems with energy regularization

We analyze space-time finite element methods for the numerical solution of distributed parabolic optimal control problems with energy regularization in the Bochner space $L^2(0,T;H^{-1}(Ω))$. By duality, the related norm can be evaluated by means of the solution of an elliptic quasi-stationary boundary value problem. When eliminating the control, we end up with the reduced optimality system that is nothing but the variational formulation of the coupled forward-backward primal and adjoint equations. Using Babuška's theorem, we prove unique solvability in the continuous case. Furthermore, we establish the discrete inf-sup condition for any conforming space-time finite element discretization yielding quasi-optimal discretization error estimates. Various numerical examples confirm the theoretical findings. We emphasize that the energy regularization results in a more localized control with sharper contours for discontinuous target functions, which is demonstrated by a comparison with an $L^2$ regularization and with a sparse optimal control approach.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Unstructured space-time finite element methods for optimal control of parabolic equations

This work presents and analyzes space-time finite element methods on fully unstructured simplicial space-time meshes for the numerical solution of parabolic optimal control problems. Using Babuška's theorem, we show well-posedness of the first-order optimality systems for a typical model problem with linear state equations, but without control constraints. This is done for both continuous and discrete levels. Based on these results, we derive discretization error estimates. Then we consider a semilinear parabolic optimal control problem arising from the Schlögl model. The associated nonlinear optimality system is solved by Newton's method, where a linear system, that is similar to the first-order optimality systems considered for the linear model problems, has to be solved at each Newton step. We present various numerical experiments including results for adaptive space-time finite element discretizations based on residual-type error indicators. In the last two examples, we also consider semilinear parabolic optimal control problems with box constraints imposed on the control.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Unstructured Space-Time Finite Element Methods for Optimal Sparse Control of Parabolic Equations

We consider a space-time finite element method on fully unstructured simplicial meshes for optimal sparse control of semilinear parabolic equations. The objective is a combination of a standard quadratic tracking-type functional including a Tikhonov regularization term and of the $L^1$-norm of the control that accounts for its spatio-temporal sparsity. We use a space-time Petrov-Galerkin finite element discretization for the first-order necessary optimality system of the associated discrete optimal sparse control problem. The discretization is based on a variational formulation that employs piecewise linear finite elements simultaneously in space and time. Finally, the discrete nonlinear optimality system that consists of coupled forward-backward state and adjoint state equations is solved by a semismooth Newton method.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Multigoal-oriented optimal control problems with nonlinear PDE constraints

In this work, we consider an optimal control problem subject to a nonlinear PDE constraint and apply it to the regularized $p$-Laplace equation. To this end, a reduced unconstrained optimization problem in terms of the control variable is formulated. Based on the reduced approach, we then derive an a posteriori error representation and mesh adaptivity for multiple quantities of interest. All quantities are combined to one, and then the dual-weighted residual (DWR) method is applied to this combined functional. Furthermore, the estimator allows for balancing the discretization error and the nonlinear iteration error. These developments allow us to formulate an adaptive solution strategy, which is finally substantiated via several numerical examples.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Discontinuous Galerkin Isogeometric Analysis on Non-matching Segmentation: Error Estimates and Efficient Solvers

The Isogeometric Analysis (IgA) of boundary value problems in complex domains often requires a decomposition of the computational domain into patches such that each of which can be parametrized by the so-called geometrical mapping. In this paper, we develop discontinuous Galerkin (dG) IgA techniques for solving elliptic diffusion problems on decompositions that can include non-matching parametrizations of the interfaces, i.e., the interfaces of the adjacent patches may be not identical. The lack of the exact parametrization of the patches leads to the creation of gap and overlapping regions between the patches. This does not allow the immediate use of the classical numerical fluxes that are known in the literature. The unknown normal fluxes of the solution on the non-matching interfaces are approximated by Taylor expansions using the values of the solution computed on the boundary of the patches These approximations are used in order to build up the numerical fluxes of the final dG IgA scheme and to couple the local patch-wise discrete problems. The resulting linear systems are solved by using efficient domainecomposition methods based on the tearing and interconnecting technology. We present numerical results of a series of test problems that validate the theoretical estimates presented.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Discontinuous Galerkin Isogeometric Analysis of Elliptic Diffusion Problems on Segmentations with Gaps

We propose a new discontinuous Galerkin Isogeometric Analysis (IgA) technique for the numerical solution of elliptic diffusion problems in computational domains decomposed into volumetric patches with non-matching interfaces. Due to an incorrect segmentation procedure, it may happen that the interfaces of adjacent subdomains don't coincide. In this way, gap regions, which are not present in the original physical domain, are created. In this paper, the gap region is considered as a subdomain of the decomposition of the computational domain and the gap boundary is taken as an interface between the gap and the subdomains. We apply a multi-patch approach and derive a subdomain variational formulation which includes interface continuity conditions and is consistent with the original variational formulation of the problem. The last formulation is further modified by deriving interface conditions without the presence of the solution in the gap. Finally, the solution of this modified problem is approximated by a special discontinuous Galerkin IgA technique. The ideas are illustrated on a model diffusion problem with discontinuous diffusion coefficients. We develop a rigorous theoretical framework for the proposed method clarifying the influence of the gap size onto the convergence rate of the method. The theoretical estimates are supported by numerical examples in two- and three-dimensional computational domains.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Dual-Primal Isogeometric Tearing and Interconnecting Solvers for large-scale systems of multipatch continuous Galerkin IgA equations

The dual-primal isogeometric tearing and interconnecting (IETI-DP) method is the adaption of the dual-primal finite element tearing and interconnecting (FETI-DP) method to isogeometric analysis of scalar elliptic boundary value problems like, e.g., diffusion problems with heterogeneous diffusion coefficients. The purpose of this paper is to extent the already existing results on condition number estimates to multi-patch domains, which consist of different geometrical mappings for each patch. Another purpose is to prove a polylogarithmic condition number bound for the preconditioned system with stiffness scaling in case of $C^0$ smoothness across patch interfaces. Numerical experiments validate the presented theory.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Functional a posteriori error estimates for time-periodic parabolic optimal control problems

This paper is devoted to the a posteriori error analysis of multiharmonic finite element approximations to distributed optimal control problems with time-periodic state equations of parabolic type. We derive a posteriori estimates of functional type, which are easily computable and provide guaranteed upper bounds for the state and co-state errors as well as for the cost functional. These theoretical results are confirmed by several numerical tests that show high efficiency of the a posteriori error bounds.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Mesh Grading in Isogeometric Analysis

This paper is concerned with the construction of graded meshes for approximating so-called singular solutions of elliptic boundary value problems by means of multipatch discontinuous Galerkin Isogeometric Analysis schemes. Such solutions appear, for instance, in domains with re-entrant corners on the boundary of the computational domain, in problems with changing boundary conditions, in interface problems, or in problems with singular source terms. Making use of the analytic behavior of the solution, we construct the graded meshes in the neighborhoods of such singular points following a multipatch approach. We prove that appropriately graded meshes lead to the same convergence rates as in the case of smooth solutions with approximately the same number of degrees of freedom. Representative numerical examples are studied in order to confirm the theoretical convergence rates and to demonstrate the efficiency of the mesh grading technology in Isogeometric Analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Numerical simulation of parabolic moving and growing interface problems using small mesh deformation

In this work, we develop a cutting method for solving problems with moving and growing interfaces in 3D. This new method is able to resolve large displacement or deformation of immersed objects by combining the Arbitrary Lagrangian-Eulerian method with only small local mesh deformation defined on the reference domain, that is decomposed into the macro-elements. The linear system of algebraic equations arising after the temporal and spatial discretizations of a model parabolic interface heat-conduction-like problem with vector-valued functions is solved by either an all-at-once or a segregated algebraic multigrid method.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Space-Time Isogeometric Analysis of Parabolic Evolution Equations

We present and analyze a new stable space-time Isogeometric Analysis (IgA) method for the numerical solution of parabolic evolution equations in fixed and moving spatial computational domains. The discrete bilinear form is elliptic on the IgA space with respect to a discrete energy norm. This property together with a corresponding boundedness property, consistency and approximation results for the IgA spaces yields an a priori discretization error estimate with respect to the discrete norm. The theoretical results are confirmed by several numerical experiments with low- and high-order IgA spaces.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Mid-Infrared Instrument for JWST, II: Design and Build

The Mid-InfraRed Instrument (MIRI) on the James Webb Space Telescope (JWST) provides measurements over the wavelength range 5 to 28.5 microns. MIRI has, within a single 'package', four key scientific functions: photometric imaging, coronagraphy, single-source low-spectral resolving power (R ~ 100) spectroscopy, and medium-resolving power (R ~ 1500 to 3500) integral field spectroscopy. An associated cooler system maintains MIRI at its operating temperature of < 6.7 K. This paper describes the driving principles behind the design of MIRI, the primary design parameters, and their realization in terms of the 'as-built' instrument. It also describes the test program that led to delivery of the tested and calibrated Flight Model to NASA in 2012, and the confirmation after delivery of the key interface requirements.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Analysis of Multipatch Discontinuous Galerkin IgA Approximations to Elliptic Boundary Value Problems

In this work, we study the approximation properties of multi-patch dG-IgA methods, that apply the multipatch Isogeometric Analysis (IgA) discretization concept and the discontinuous Galerkin (dG) technique on the interfaces between the patches, for solving linear diffusion problems with diffusion coefficients that may be discontinuous across the patch interfaces. The computational domain is divided into non-overlapping sub-domains, called patches in IgA, where $B$-splines, or NURBS finite dimensional approximations spaces are constructed. The solution of the problem is approximated in every sub-domain without imposing any matching grid conditions and without any continuity requirements for the discrete solution across the interfaces. Numerical fluxes with interior penalty jump terms are applied in order to treat the discontinuities of the discrete solution on the interfaces. We provide a rigorous a priori discretization error analysis for problems set in 2d- and 3d- dimensional domains, with solutions belonging to $W^{l,p}, l\geq 2,{\ } p\in ({2d}/{(d+2(l-1))},2]$. In any case, we show optimal convergence rates of the discretization with respect to the dG - norm.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Discontinuous Galerkin Isogeometric Analysis of Elliptic PDEs on Surfaces

Isogeometric analysis uses the same class of basis functions for both, representing the geometry of the computational domain and approximating the solution. In practical applications, geometrical patches are used in order to get flexibility in the geometrical representation. This patch representation corresponds to a domain decomposition. In this paper, we present for the first time a Discontinuous Galerkin (DG) Method that allows for discontinuities only along the sub-domain (patch) boundaries. The required smoothness is obtained by the DG terms associated with the boundary of the sub-domains. The construction and corresponding discretization error analysis of such DG scheme is presented for Elliptic PDEs in a 2D as well as on open and closed surfaces. Furthermore, we present numerical results to confirm the theory presented.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Functional a posteriori error estimates for parabolic time-periodic boundary value problems

The paper is concerned with parabolic time-periodic boundary value problems which are of theoretical interest and arise in different practical applications. The multiharmonic finite element method is well adapted to this class of parabolic problems. We study properties of multiharmonic approximations and derive guaranteed and fully computable bounds of approximation errors. For this purpose, we use the functional a posteriori error estimation techniques earlier introduced by S. Repin. Numerical tests confirm the efficiency of the a posteriori error bounds derived.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Multipatch Discontinuous Galerkin Isogeometric Analysis

Isogeometric analysis (IgA) uses the same class of basis functions for both, representing the geometry of the computational domain and approximating the solution. In practical applications, geometrical patches are used in order to get flexibility in the geometrical representation. This multi-patch representation corresponds to a decomposition of the computational domain into non-overlapping subdomains also called patches in the geometrical framework. We will present discontinuous Galerkin (dG) methods that allow for discontinuities across the subdomain (patch) boundaries. The required interface conditions are weakly imposed by the dG terms associated with the boundary of the sub-domains. The construction and the corresponding discretization error analysis of such dG multi-patch IgA schemes will be given for heterogeneous diffusion model problems in volumetric 2d and 3d domains as well as on open and closed surfaces. The theoretical results are confirmed by numerous numerical experiments which have been performed in G+SMO. The concept and the main features of the IgA library G+SMO are also described.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Numerical Simulation of Fluid-Structure Interaction Problems with Hyperelastic Models: A Monolithic Approach

In this paper, we consider a monolithic approach to handle coupled fluid-structure interaction problems with different hyperelastic models in an all-at-once manner. We apply Newton's method in the outer iteration dealing with nonlinearities of the coupled system. We discuss preconditioned Krylov sub-space, algebraic multigrid and algebraic multilevel methods for solving the linearized algebraic equations. Finally, we compare the results of the monolithic approach with those of the corresponding partitioned approach that we studied in our previous work.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Numerical Simulation of Fluid-Structure Interaction Problems with Hyperelastic Models I: A Partitioned Approach

In this work, we consider fluid-structure interaction simulation with nonlinear hyperelastic models in the solid part. We use a partitioned approach to deal with the coupled nonlinear fluid-structure interaction problems. We focus on handling the nonlinearity of the fluid and structure sub-problems, the near-incompressibility of materials, the stabilization of employed finite element discretization, and the robustness and efficiency of Krylov subspace and algebraic multigrid methods for the linearized algebraic equations.

0 citations0 reviewsNo code linkedOpen access