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

Patrick Henning

Patrick Henning 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-phmath.MPmath.SP

Catalog footprint

What is connected

19works
5topics
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

19 published item(s)

preprint2026arXiv

Nonlinear Inverse Iterations for Spin-Orbit Coupled Quantum Gases

This work concerns the computation of ground states of two-component spin-orbit coupled Bose-Einstein condensates (SO-coupled BECs), modelled by a coupled nonlinear eigenvalue problem of Gross-Pitaevskii type. Spin-orbit coupling gives rise to fascinating phenomena, including supersolid-like phases with spatially modulated densities. However, in such complex settings, conventional numerical approaches, such as generalized inverse iterations or gradient descent, often converge very slowly. To overcome this issue, we apply the concept of the J-method [E.~Jarlebring, S.~Kvaal, W.~Michiels. SIAM~J.~Sci.~Comput.~36-4,~2014] to construct a nonlinear inverse iteration scheme whose convergence can be accelerated through spectral shifting, analogous to techniques used for linear eigenproblems. For a fixed shift parameter, we establish local linear convergence rates determined by spectral gaps in the neighbourhood of each quasi-unique ground state. With adaptively chosen shifts, superlinear convergence is observed, which we verify through numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Shadow Lagrangian dynamics for superfluidity

Motivated by a similar approach for Born-Oppenheimer molecular dynamics, this paper proposes an extended "shadow" Lagrangian density for quantum states of superfluids. The extended Lagrangian contains an additional field variable that is forced to follow the wave function of the quantum state through a rapidly oscillating extended harmonic oscillator. By considering the adiabatic limit for large frequencies of the harmonic oscillator, we can derive the two equations of motions, a Schrödinger-type equation for the quantum state and a wave equation for the extended field variable. The equations are coupled in a nonlinear way, but each equation individually is linear with respect to the variable that it defines. The computational advantage of this new system is that it can be easily discretized using linear time stepping methods, where we propose to use a Crank-Nicolson-type approach for the Schrödinger equation and an extended leapfrog scheme for the wave equation. Furthermore, the difference between the quantum state and the extended field variable defines a consistency error that should go to zero if the frequency tends to infinity. By coupling the time-step size in our discretization to the frequency of the harmonic oscillator we can extract an easily computable consistency error indicator that can be used to estimate the numerical error without any additional costs. The findings are illustrated in numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A note on optimal $H^1$-error estimates for Crank-Nicolson approximations to the nonlinear Schrödinger equation

In this paper we consider a mass- and energy--conserving Crank-Nicolson time discretization for a general class of nonlinear Schrödinger equations. This scheme, which enjoys popularity in the physics community due to its conservation properties, was already subject to several analytical and numerical studies. However, a proof of optimal $L^{\infty}(H^1)$-error estimates is still open, both in the semi-discrete Hilbert space setting, as well as in fully-discrete finite element settings. This paper aims at closing this gap in the literature. We also suggest a fixed point iteration to solve the arising nonlinear system of equations that makes the method easy to implement and efficient. This is illustrated by numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Computational homogenization of time-harmonic Maxwell's equations

In this paper we consider a numerical homogenization technique for curl-curl-problems that is based on the framework of the Localized Orthogonal Decomposition and which was proposed in [D. Gallistl, P. Henning, B. Verfürth. SIAM J. Numer. Anal. 56-3:1570-1596, 2018] for problems with essential boundary conditions. The findings of the aforementioned work establish quantitative homogenization results for the time-harmonic Maxwell's equations that hold beyond assumptions of periodicity, however, a practical realization of the approach was left open. In this paper, we transfer the findings from essential boundary conditions to natural boundary conditions and we demonstrate that the approach yields a computable numerical method. We also investigate how boundary values of the source term can effect the computational complexity and accuracy. Our findings will be supported by various numerical experiments, both in $2D$ and $3D$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Quantitative Anderson localization of Schrödinger eigenstates under disorder potentials

This paper concerns spectral properties of linear Schrödinger operators under oscillatory high-amplitude potentials on bounded domains. Depending on the degree of disorder, we prove the existence of spectral gaps amongst the lowermost eigenvalues and the emergence of exponentially localized states. We quantify the rate of decay in terms of geometric parameters that characterize the potential. The proofs are based on the convergence theory of iterative solvers for eigenvalue problems and their optimal local preconditioning by domain decomposition.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Sobolev gradient flow for the Gross-Pitaevskii eigenvalue problem: global convergence and computational efficiency

We propose a new normalized Sobolev gradient flow for the Gross-Pitaevskii eigenvalue problem based on an energy inner product that depends on time through the density of the flow itself. The gradient flow is well-defined and converges to an eigenfunction. For ground states we can quantify the convergence speed as exponentially fast where the rate depends on spectral gaps of a linearized operator. The forward Euler time discretization of the flow yields a numerical method which generalizes the inverse iteration for the nonlinear eigenvalue problem. For sufficiently small time steps, the method reduces the energy in every step and converges globally in $H^1$ to an eigenfunction. In particular, for any nonnegative starting value, the ground state is obtained. A series of numerical experiments demonstrates the computational efficiency of the method and its competitiveness with established discretizations arising from other gradient flows for this problem.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A multiscale method for linear elasticity reducing Poisson locking

We propose a generalized finite element method for linear elasticity equations with highly varying and oscillating coefficients. The method is formulated in the framework of localized orthogonal decomposition techniques introduced by Målqvist and Peterseim (Math. Comp., 83(290): 2583--2603, 2014). Assuming only $L_\infty$-coefficients we prove linear convergence in the $H^1$-norm, also for materials with large Lamé parameter $λ$. The theoretical a priori error estimate is confirmed by numerical examples.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Multiscale methods for wave problems in heterogeneous media

In this paper we give a survey on various multiscale methods for the numerical solution of second order hyperbolic equations in highly heterogeneous media. We concentrate on the wave equation and distinguish between two classes of applications. First we discuss numerical methods for the wave equation in heterogeneous media without scale separation. Such a setting is for instance encountered in the geosciences, where natural structures often exhibit a continuum of different scales, that all need to be resolved numerically to get meaningful approximations. Approaches tailored for these settings typically involve the construction of generalized finite element spaces, where the basis functions incorporate information about the data variations. In the second part of the paper, we discuss numerical methods for the case of structured media with scale separation. This setting is for instance encountered in engineering sciences, where materials are often artificially designed. If this is the case, the structure and the scale separation can be explicitly exploited to compute appropriate homogenized/upscaled wave models that only exhibit a single coarse scale and that can be hence solved at significantly reduced computational costs.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Multiscale mixed finite elements

In this work, we propose a mixed finite element method for solving elliptic multiscale problems based on a localized orthogonal decomposition (LOD) of Raviart-Thomas finite element spaces. It requires to solve local problems in small patches around the elements of a coarse grid. These computations can be perfectly parallelized and are cheap to perform. Using the results of these patch problems, we construct a low dimensional multiscale mixed finite element space with very high approximation properties. This space can be used for solving the original saddle point problem in an efficient way. We prove convergence of our approach, independent of structural assumptions or scale separation. Finally, we demonstrate the applicability of our method by presenting a variety of numerical experiments, including a comparison with an MsFEM approach.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Finite Element Method for the time-dependent Gross-Pitaevskii equation with angular momentum rotation

We consider the time-dependent Gross-Pitaevskii equation describing the dynamics of rotating Bose-Einstein condensates and its discretization with the finite element method. We analyze a mass conserving Crank-Nicolson-type discretization and prove corresponding a priori error estimates with respect to the maximum norm in time and the $L^2$- and energy-norm in space. The estimates show that we obtain optimal convergence rates under the assumption of additional regularity for the solution to the Gross-Pitaevskii equation. We demonstrate the performance of the method in numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A new Heterogeneous Multiscale Method for time-harmonic Maxwell's equations based on divergence-regularization

In this paper, we suggest a new heterogeneous multiscale method (HMM) for the time-harmonic Maxwell equations in locally periodic media. The method is constructed by using a divergence-regularization in one of the cell problems. This allows us to introduce fine-scale correctors that are not subject to a cumbersome divergence-free constraint and which can hence easily be implemented. To analyze the method, we first revisit classical homogenization theory for time-harmonic Maxwell equations and derive a new homogenization result that makes use of the divergence-regularization in the two-scale homogenized equation. We then show that the HMM is equivalent to a discretization of this equation. In particular, writing both problems in a fully coupled two-scale formulation is the crucial starting point for a corresponding numerical analysis of the method. With this approach we are able to prove rigorous a priori error estimates in the $\mathbf{H}(\mbox{curl})$- and the $H^{-1}$-norm and we derive reliable and efficient localized residual-based a posteriori error estimates.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A reduced basis localized orthogonal decomposition

In this work we combine the framework of the Reduced Basis method (RB) with the framework of the Localized Orthogonal Decomposition (LOD) in order to solve parametrized elliptic multiscale problems. The idea of the LOD is to split a high dimensional Finite Element space into a low dimensional space with comparably good approximation properties and a remainder space with negligible information. The low dimensional space is spanned by locally supported basis functions associated with the node of a coarse mesh obtained by solving decoupled local problems. However, for parameter dependent multiscale problems, the local basis has to be computed repeatedly for each choice of the parameter. To overcome this issue, we propose an RB approach to compute in an "offline" stage LOD for suitable representative parameters. The online solution of the multiscale problems can then be obtained in a coarse space (thanks to the LOD decomposition) and for an arbitrary value of the parameters (thanks to a suitable "interpolation" of the selected RB). The online RB-LOD has a basis with local support and leads to sparse systems. Applications of the strategy to both linear and nonlinear problems are given.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Localized orthogonal decomposition method for the wave equation with a continuum of scales

This paper is devoted to numerical approximations for the wave equation with a multiscale character. Our approach is formulated in the framework of the Localized Orthogonal Decomposition (LOD) interpreted as a numerical homogenization with an $L^2$-projection. We derive explicit convergence rates of the method in the $L^{\infty}(L^2)$-, $W^{1,\infty}(L^2)$- and $L^{\infty}(H^1)$-norms without any assumptions on higher order space regularity or scale-separation. The order of the convergence rates depends on further graded assumptions on the initial data. We also prove the convergence of the method in the framework of G-convergence without any structural assumptions on the initial data, i.e. without assuming that it is well-prepared. This rigorously justifies the method. Finally, the performance of the method is demonstrated in numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Adaptive Heterogeneous Multiscale Methods for immiscible two-phase flow in porous media

In this contribution we present the first formulation of a heterogeneous multiscale method for an incompressible immiscible two-phase flow system with degenerate permeabilities. The method is in a general formulation which includes oversampling. We do not specify the discretization of the derived macroscopic equation, but we give two examples of possible realizations, suggesting a finite element solver for the fine scale and a vertex centered finite volume method for the effective coarse scale equations. Assuming periodicity, we show that the method is equivalent to a discretization of the homogenized equation. We provide an a-posteriori estimate for the error between the homogenized solutions of the pressure and saturation equations and the corresponding HMM approximations. The error estimate is based on the results recently achieved in [C. Canc{è}s, I. S. Pop, and M. Vohral\'ık. An a posteriori error estimate for vertex-centered finite volume discretizations of immiscible incompressible two-phase flow. Math. Comp., 2014].

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Localized orthogonal decomposition techniques for boundary value problems

In this paper we propose a Local Orthogonal Decomposition method (LOD) for elliptic partial differential equations with inhomogeneous Dirichlet- and Neumann boundary conditions. For this purpose, we present new boundary correctors which preserve the common convergence rates of the LOD, even if the boundary condition has a rapidly oscillating fine scale structure. We prove a corresponding a-priori error estimate and present numerical experiments. We also demonstrate numerically that the method is reliable with respect to thin conductivity channels in the diffusion matrix. Accurate results are obtained without resolving these channels by the coarse grid and without using patches that contain the channels.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Multiscale Partition of Unity

We introduce a new Partition of Unity Method for the numerical homogenization of elliptic partial differential equations with arbitrarily rough coefficients. We do not restrict to a particular ansatz space or the existence of a finite element mesh. The method modifies a given partition of unity such that optimal convergence is achieved independent of oscillation or discontinuities of the diffusion coefficient. The modification is based on an orthogonal decomposition of the solution space while preserving the partition of unity property. This precomputation involves the solution of independent problems on local subdomains of selectable size. We deduce quantitative error estimates for the method that account for the chosen amount of localization. Numerical experiments illustrate the high approximation properties even for 'cheap' parameter choices.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On Multiscale Methods in Petrov-Galerkin formulation

In this work we investigate the advantages of multiscale methods in Petrov-Galerkin (PG) formulation in a general framework. The framework is based on a localized orthogonal decomposition of a high dimensional solution space into a low dimensional multiscale space with good approximation properties and a high dimensional remainder space{, which only contains negligible fine scale information}. The multiscale space can then be used to obtain accurate Galerkin approximations. As a model problem we consider the Poisson equation. We prove that a Petrov-Galerkin formulation does not suffer from a significant loss of accuracy, and still preserve the convergence order of the original multiscale method. We also prove inf-sup stability of a PG Continuous and a Discontinuous Galerkin Finite Element multiscale method. Furthermore, we demonstrate that the Petrov-Galerkin method can decrease the computational complexity significantly, allowing for more efficient solution algorithms. As another application of the framework, we show how the Petrov-Galerkin framework can be used to construct a locally mass conservative solver for two-phase flow simulation that employs the Buckley-Leverett equation. To achieve this, we couple a PG Discontinuous Galerkin Finite Element method with an upwind scheme for a hyperbolic conservation law.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Two-Level discretization techniques for ground state computations of Bose-Einstein condensates

This work presents a new methodology for computing ground states of Bose-Einstein condensates based on finite element discretizations on two different scales of numerical resolution. In a pre-processing step, a low-dimensional (coarse) generalized finite element space is constructed. It is based on a local orthogonal decomposition and exhibits high approximation properties. The non-linear eigenvalue problem that characterizes the ground state is solved by some suitable iterative solver exclusively in this low-dimensional space, without loss of accuracy when compared with the solution of the full fine scale problem. The pre-processing step is independent of the types and numbers of bosons. A post-processing step further improves the accuracy of the method. We present rigorous a priori error estimates that predict convergence rates H^3 for the ground state eigenfunction and H^4 for the corresponding eigenvalue without pre-asymptotic effects; H being the coarse scale discretization parameter. Numerical experiments indicate that these high rates may still be pessimistic.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Oversampling for the Multiscale Finite Element Method

This paper reviews standard oversampling strategies as performed in the Multiscale Finite Element Method (MsFEM). Common to those approaches is that the oversampling is performed in the full space restricted to a patch but including coarse finite element functions. We suggest, by contrast, to perform local computations with the additional constraint that trial and test functions are linear independent from coarse finite element functions. This approach re-interprets the Variational Multiscale Method in the context of computational homogenization. This connection gives rise to a general fully discrete error analysis for the proposed multiscale method with constrained oversampling without any resonance effects. In particular, we are able to give the first rigorous proof of convergence for a MsFEM with oversampling.

0 citations0 reviewsNo code linkedOpen access