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

Axel Målqvist

Axel Målqvist appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NANumerical Analysis

Catalog footprint

What is connected

16works
2topics
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

16 published item(s)

preprint2022arXiv

Iterative solution of spatial network models by subspace decomposition

We present and analyze a preconditioned conjugate gradient method (PCG) for solving spatial network problems. Primarily, we consider diffusion and structural mechanics simulations for fiber based materials, but the methodology can be applied to a wide range of models, fulfilling a set of abstract assumptions. The proposed method builds on a classical subspace decomposition into a coarse subspace, realized as the restriction of a finite element space to the nodes of the spatial network, and localized subspaces with support on mesh stars. The main contribution of this work is the convergence analysis of the proposed method. The analysis translates results from finite element theory, including interpolation bounds, to the spatial network setting. A convergence rate of the PCG algorithm, only depending on global bounds of the operator and homogeneity, connectivity and locality constants of the network, is established. The theoretical results are confirmed by several numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Numerical homogenization of spatial network models

We present and analyze a methodology for numerical homogenization of spatial networks, modelling e.g. diffusion processes and deformation of mechanical structures. The aim is to construct an accurate coarse model of the network. By solving decoupled problems on local subgraphs we construct a low dimensional subspace of the solution space with good approximation properties. The coarse model of the network is expressed by a Galerkin formulation and can be used to perform simulations with different source and boundary data at a low computational cost. We prove optimal convergence of the proposed method under mild assumptions on the homogeneity, connectivity, and locality of the network on the coarse scale. The theoretical findings are numerically confirmed for both scalar-valued (heat conduction) and vector-valued (structural) models.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Well-Posedness and Finite Element Approximation of Mixed Dimensional Partial Differential Equations

We consider a mixed dimensional elliptic partial differential equation posed in a bulk domain with a large number of embedded interfaces. In particular, we study well-posedness of the problem and regularity of the solution. We also propose a fitted finite element approximation and prove an a priori error bound. For the solution of the arising linear system we propose and analyze an iterative method based on subspace decomposition. Finally, we present numerical experiments and achieve rapid convergence using the proposed preconditioner, confirming our theoretical findings.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A numerical multiscale method for fiber networks

Fiber network modeling can be used for studying mechanical properties of paper. The individual fibers and the bonds in-between constitute a detailed representation of the material. However, detailed microscale fiber network models must be resolved with efficient numerical methods. In this work, a numerical multiscale method for discrete network models is proposed that is based on the localized orthogonal decomposition method [arXiv:1810.05059]. The method is ideal for these network problems, because it reduces the maximum size of the problem, it is suitable for parallelization, and it can effectively solve fracture propagation. The problem analyzed in this work is the nodal displacement of a fiber network given an applied load. This problem is formulated as a linear system that is solved by using the aforementioned multiscale method. To solve the linear system, the multiscale method constructs a low-dimensional solution space with good approximation properties. The method is observed to work well for unstructured fiber networks, with optimal rates of convergence obtainable for highly localized configurations of the method.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Numerical upscaling for heterogeneous materials in fractured domains

We consider numerical solution of elliptic problems with heterogeneous diffusion coefficients containing thin highly conductive structures. Such problems arise e.g. in fractured porous media, reinforced materials, and electric circuits. The main computational challenge is the high resolution needed to resolve the data variation. We propose a multiscale method that models the thin structures as interfaces and incorporate heterogeneities in corrected shape functions. The construction results in an accurate upscaled representation of the system that can be used to solve for several forcing functions or to simulate evolution problems in an efficient way. By introducing a novel interpolation operator, defining the fine scale of the problem, we prove exponential decay of the shape functions which allows for a sparse approximation of the upscaled representation. An a priori error bound is also derived for the proposed method together with numerical examples that verify the theoretical findings. Finally we present a numerical example to show how the technique can be applied to evolution problems.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A generalized finite element method for linear thermoelasticity

We propose and analyze a generalized finite element method designed for linear quasistatic thermoelastic systems with spatial multiscale coefficients. The method is based on the local orthogonal decomposition technique introduced by Målqvist and Peterseim (Math. Comp., 83(290): 2583--2603, 2014). We prove convergence of optimal order, independent of the derivatives of the coefficients, in the spatial $H^1$-norm. The theoretical results are confirmed by numerical examples.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Multiscale methods for problems with complex geometry

We propose a multiscale method for elliptic problems on complex domains, e.g. domains with cracks or complicated boundary. For local singularities this paper also offers a discrete alternative to enrichment techniques such as XFEM. We construct corrected coarse test and trail spaces which takes the fine scale features of the computational domain into account. The corrections only need to be computed in regions surrounding fine scale geometric features. We achieve linear convergence rate in energy norm for the multiscale solution. Moreover, the conditioning of the resulting matrices is not affected by the way the domain boundary cuts the coarse elements in the background mesh. The analytical findings are verified in a series of numerical experiments.

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

Generalized finite element methods for quadratic eigenvalue problems

We consider a large-scale quadratic eigenvalue problem (QEP), formulated using P1 finite elements on a fine scale reference mesh. This model describes damped vibrations in a structural mechanical system. In particular we focus on problems with rapid material data variation, e.g., composite materials. We construct a low dimensional generalized finite element (GFE) space based on the localized orthogonal decomposition (LOD) technique. The construction involves the (parallel) solution of independent localized linear Poisson-type problems. The GFE space is then used to compress the large-scale algebraic QEP to a much smaller one with a similar modeling accuracy. The small scale QEP can then be solved by standard techniques at a significantly reduced computational cost. We prove convergence with rate for the proposed method and numerical experiments confirm our theoretical findings.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Multiscale techniques for parabolic equations

We use the local orthogonal decomposition technique to derive a generalized finite element method for linear and semilinear parabolic equations with spatial multiscale diffusion coefficient. We consider nonsmooth initial data and a backward Euler scheme for the temporal discretization. Optimal order convergence rate, depending only on the contrast, but not on the variations in the diffusion coefficient, is proven in the $L_\infty(L_2)$-norm. We present numerical examples, which confirm our theoretical findings.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A multilevel Monte Carlo method for computing failure probabilities

We propose and analyze a method for computing failure probabilities of systems modeled as numerical deterministic models (e.g., PDEs) with uncertain input data. A failure occurs when a functional of the solution to the model is below (or above) some critical value. By combining recent results on quantile estimation and the multilevel Monte Carlo method we develop a method which reduces computational cost without loss of accuracy. We show how the computational cost of the method relates to error tolerance of the failure probability. For a wide and common class of problems, the computational cost is asymptotically proportional to solving a single accurate realization of the numerical model, i.e., independent of the number of samples. Significant reductions in computational cost are also observed in numerical experiments.

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

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
preprint2012arXiv

Convergence of a discontinuous Galerkin multiscale method

A convergence result for a discontinuous Galerkin multiscale method for a second order elliptic problem is presented. We consider a heterogeneous and highly varying diffusion coefficient in $L^\infty(Ω,\mathbb{R}^{d\times d}_{sym})$ with uniform spectral bounds and without any assumption on scale separation or periodicity. The multiscale method uses a corrected basis that is computed on patches/subdomains. The error, due to truncation of corrected basis, decreases exponentially with the size of the patches. Hence, to achieve an algebraic convergence rate of the multiscale solution on a uniform mesh with mesh size $H$ to a reference solution, it is sufficient to choose the patch sizes as $\mathcal{O}(H|\log(H^{-1})|)$. We also discuss a way to further localize the corrected basis to element-wise support leading to a slight increase of the dimension of the space. Improved convergence rate can be achieved depending on the piecewise regularity of the forcing function. Linear convergence in energy norm and quadratic convergence in $L^2$-norm is obtained independently of the forcing function. A series of numerical experiments confirms the theoretical rates of convergence.

0 citations0 reviewsNo code linkedOpen access