| 1. |
- Altmann, Robert, et al.
(author)
-
A decoupling and linearizing discretization for poroelasticity with nonlinear permeability
- 2021
-
In: ArXiv Preprint 2104.10092.
-
Other publication (other academic/artistic)abstract
- We analyze a semi-explicit time discretization scheme of first order for poroelasticity with nonlinear permeability provided that the elasticity model and the flow equation are only weakly coupled. The approach leads to a decoupling of the equations and, at the same time, linearizes the nonlinearity without the need of further inner iteration steps. Hence, the computational speed-up is twofold without a loss in the convergence rate. We prove optimal first-order error estimates by considering a related delay system and investigate the method numerically for different examples with various types of nonlinear displacement-permeability relations.
|
|
| 2. |
- Altmann, R., et al.
(author)
-
Semi-explicit discretization schemes for weakly coupled elliptic-parabolic problems
- 2021
-
In: Mathematics of Computation. - : American Mathematical Society (AMS). - 0025-5718 .- 1088-6842. ; 90:329, s. 1089-1118
-
Journal article (peer-reviewed)abstract
- We prove first-order convergence of the semi-explicit Euler scheme combined with a finite element discretization in space for elliptic-parabolic problems which are weakly coupled. This setting includes poroelasticity, thermoelasticity, as well as multiple-network models used in medical applications. The semi-explicit approach decouples the system such that each time step requires the solution of two small and well-structured linear systems rather than the solution of one large system. The decoupling improves the computational efficiency without decreasing the convergence rates. The presented convergence proof is based on an interpretation of the scheme as an implicit method applied to a constrained partial differential equation with delay term. Here, the delay time equals the used step size. This connection also allows a deeper understanding of the weak coupling condition, which we accomplish to quantify explicitly.
|
|
| 3. |
- Caiazzo, A., et al.
(author)
-
Reconstruction of Quasi-Local Numerical Effective Models from Low-Resolution Measurements
- 2020
-
In: Journal of Scientific Computing. - : Springer Science and Business Media LLC. - 0885-7474 .- 1573-7691. ; 85:1
-
Journal article (peer-reviewed)abstract
- We consider the inverse problem of reconstructing an effective model for a prototypical diffusion process in strongly heterogeneous media based on coarse measurements. The approach is motivated by quasi-local numerical effective forward models that are provably reliable beyond periodicity assumptions and scale separation. The goal of this work is to show that an identification of the matrix representation related to these effective models is possible. On the one hand, this provides a reasonable surrogate in cases where a direct reconstruction is unfeasible due to a mismatch between the coarse data scale and the microscopic quantities to be reconstructed. On the other hand, the approach allows us to investigate the requirement for a certain non-locality in the context of numerical homogenization. Algorithmic aspects of the inversion procedure and its performance are illustrated in a series of numerical experiments.
|
|
| 4. |
- Geevers, Sjoerd, et al.
(author)
-
Fast mass lumped multiscale wave propagation modelling
- 2021
-
In: ArXiv Preprint 2104.08346.
-
Other publication (other academic/artistic)abstract
- In this paper, we investigate the use of a mass lumped fully explicit time stepping scheme for the discretisation of the wave equation with underlying material parameters that vary at arbitrarily fine scales. We combine the leapfrog scheme for the temporal discretisation with the multiscale technique known as Localized Orthogonal Decomposition for the spatial discretisation. To speed up the method and to make it fully explicit, a special mass lumping approach is introduced that relies on an appropriate interpolation operator. This operator is also employed in the construction of the Localized Orthogonal Decomposition and is a key feature of the approach. We prove that the method converges with second order in the energy norm, with a leading constant that does not depend on the scales at which the material parameters vary. We also illustrate the performance of the mass lumped method in a set of numerical experiments.
|
|
| 5. |
- Kröpfl, Fabian, et al.
(author)
-
Operator compression with deep neural networks
- 2021
-
In: ArXiv Preprint 2105.12080.
-
Other publication (other academic/artistic)abstract
- This paper studies the compression of partial differential operators using neural networks. We consider a family of operators, parameterized by a potentially high-dimensional space of coefficients that may vary on a large range of scales. Based on existing methods that compress such a multiscale operator to a finite-dimensional sparse surrogate model on a given target scale, we propose to directly approximate the coefficient-to-surrogate map with a neural network. We emulate local assembly structures of the surrogates and thus only require a moderately sized network that can be trained efficiently in an offline phase. This enables large compression ratios and the online computation of a surrogate based on simple forward passes through the network is substantially accelerated compared to classical numerical upscaling approaches. We apply the abstract framework to a family of prototypical second-order elliptic heterogeneous diffusion operators as a demonstrating example.
|
|
| 6. |
- Ljung, Per, 1994, et al.
(author)
-
A space-time multiscale method for parabolic problems
- 2022
-
In: Multiscale Modeling and Simulation. - 1540-3467 .- 1540-3459. ; 20:2, s. 714-740
-
Journal article (peer-reviewed)abstract
- We present a space-time multiscale method for a parabolic model problem with an underlying coefficient that may be highly oscillatory with respect to both the spatial and the temporal variables. The method is based on the framework of the Variational Multiscale Method in the context of a space-time formulation and computes a coarse-scale representation of the differential operator that is enriched by auxiliary space-time corrector functions. Once computed, the coarse-scale representation allows us to efficiently obtain well-approximating discrete solutions for multiple right-hand sides. We prove first-order convergence independently of the oscillation scales in the coefficient and illustrate how the space-time correctors decay exponentially in both space and time, making it possible to localize the corresponding computations. This localization allows us to define a practical and computationally efficient method in terms of complexity and memory, for which we provide a posteriori error estimates and present numerical examples.
|
|
| 7. |
- Ljung, Per, et al.
(author)
-
A space-time multiscale method for parabolic problems
- 2021
-
In: ArXiv Preprint 2109.06647.
-
Other publication (other academic/artistic)abstract
- We present a space-time multiscale method for a parabolic model problem with an underlying coefficient that may be highly oscillatory with respect to both the spatial and the temporal variables. The method is based on the framework of the Variational Multiscale Method in the context of a space-time formulation and computes a coarse-scale representation of the differential operator that is enriched by auxiliary space-time corrector functions. Once computed, the coarse-scale representation allows us to efficiently obtain well-approximating discrete solutions for multiple right-hand sides. We prove first-order convergence independently of the oscillation scales in the coefficient and illustrate how the space-time correctors decay exponentially in both space and time, making it possible to localize the corresponding computations. This localization allows us to define a practical and computationally efficient method in terms of complexity and memory, for which we provide a posteriori error estimates and present numerical examples.
|
|
| 8. |
- Maier, Roland, 1993
(author)
-
A high-order approach to elliptic multiscale problems with general unstructured coefficients
- 2021
-
In: SIAM Journal on Numerical Analysis. - 1095-7170 .- 0036-1429. ; 59:2, s. 1067-1089
-
Journal article (peer-reviewed)abstract
- We propose a multiscale approach for an elliptic multiscale setting with general unstructured diffusion coefficients that is able to achieve high-order convergence rates with respect to the mesh parameter and the polynomial degree. The method allows for suitable localization and does not rely on additional regularity assumptions on the domain, the diffusion coefficient, or the exact (weak) solution as typically required for high-order approaches. Rigorous a priori error estimates are presented with respect to the involved discretization parameters, and the interplay between these parameters as well as the performance of the method are studied numerically.
|
|
| 9. |
- Maier, Roland, 1993, et al.
(author)
-
Adaptive refinement for unstructured T-splines with linear complexity
- 2021
-
In: ArXiv Preprint 2109.00448.
-
Other publication (other academic/artistic)abstract
- We present an adaptive refinement algorithm for T-splines on unstructured 2D meshes. While for structured 2D meshes, one can refine elements alternatingly in horizontal and vertical direction, such an approach cannot be generalized directly to unstructured meshes, where no two unique global mesh directions can be assigned. To resolve this issue, we introduce the concept of direction indices, i.e., integers associated to each edge, which are inspired by theory on higher-dimensional structured T-splines. Together with refinement levels of edges, these indices essentially drive the refinement scheme. We combine these ideas with an edge subdivision routine that allows for I-nodes, yielding a very flexible refinement scheme that nicely distributes the T-nodes, preserving global linear independence, analysis-suitability (local linear independence) except in the vicinity of extraordinary nodes, sparsity of the system matrix, and shape regularity of the mesh elements. Further, we show that the refinement procedure has linear complexity in the sense of guaranteed upper bounds on a) the distance between marked and additionally refined elements, and on b) the ratio of the numbers of generated and marked mesh elements.
|
|
| 10. |
- Maier, Roland, 1993, et al.
(author)
-
Multiscale Scattering in Nonlinear Kerr-Type Media
- 2020
-
In: ArXiv Preprint 2011.09168.
-
Other publication (other academic/artistic)abstract
- We propose a multiscale approach for a nonlinear Helmholtz problem with possible oscillations in the Kerr coefficient, the refractive index, and the diffusion coefficient. The method does not rely on structural assumptions on the coefficients and combines the multiscale technique known as Localized Orthogonal Decomposition with an adaptive iterative approximation of the nonlinearity. We rigorously analyze the method in terms of well-posedness and convergence properties based on suitable assumptions on the initial data and the discretization parameters. Numerical examples illustrate the theoretical error estimates and underline the practicability of the approach.
|
|
