1. |
- Appelo, Daniel, et al.
(författare)
-
Waveholtz : Iterative solution of the helmholtz equation via the wave equation
- 2020
-
Ingår i: SIAM Journal on Scientific Computing. - : Society for Industrial & Applied Mathematics (SIAM). - 1064-8275 .- 1095-7197. ; 42:4, s. A1950-A1983
-
Tidskriftsartikel (refereegranskat)abstract
- A new iterative method, the WaveHoltz iteration, for solution of the Helmholtz equation is presented. WaveHoltz is a fixed point iteration that filters the solution to the solution of a wave equation with time periodic forcing and boundary data. The WaveHoltz iteration corresponds to a linear and coercive operator which, after discretization, can be recast as a positive definite linear system of equations. The solution to this system of equations approximates the Helmholtz solution and can be accelerated by Krylov subspace techniques. Analysis of the continuous and discrete cases is presented, as are numerical experiments.
|
|
2. |
- Appelö, D., et al.
(författare)
-
El-WaveHoltz : A time-domain iterative solver for time-harmonic elastic waves
- 2022
-
Ingår i: Computer Methods in Applied Mechanics and Engineering. - : Elsevier BV. - 0045-7825 .- 1879-2138. ; 401
-
Tidskriftsartikel (refereegranskat)abstract
- We consider the application of the WaveHoltz iteration to time-harmonic elastic wave equations with energy conserving boundary conditions. The original WaveHoltz iteration for acoustic Helmholtz problems is a fixed-point iteration that filters the solution of the wave equation with time-harmonic forcing and boundary data. As in the original WaveHoltz method, we reformulate the fixed point iteration as a positive definite linear system of equations that is iteratively solved by a Krylov method. We present two time-stepping schemes, one explicit and one (novel) implicit, which completely remove time discretization error from the WaveHoltz solution by performing a simple modification of the initial data and time-stepping scheme. Numerical experiments indicate an iteration scaling similar to that of the original WaveHoltz method, and that the convergence rate is dictated by the shortest (shear) wave speed of the problem. We additionally show that the implicit scheme can be advantageous in practice for meshes with disparate element sizes.
|
|
3. |
- Appelö, Daniel, et al.
(författare)
-
WaveHoltz : Parallel and scalable solution of the Helmholtz equation via wave equation iteration
- 2020
-
Ingår i: SEG International Exposition and Annual Meeting 2019. - : Society of Exploration Geophysicists. ; , s. 1541-1545
-
Konferensbidrag (refereegranskat)abstract
- We introduce a novel idea, the WaveHoltz iteration, for solving the Helmholtz equation inspired by recent work on exact controllability (EC) methods. As in EC methods our method make use of time domain methods for wave equations to design frequency domain Helmholtz solvers but unlike EC methods we do not require adjoint solves. We show that the WaveHoltz iteration we propose is symmetric and positive definite in the continuous setting. We also present numerical examples, using various discretization techniques, that show that our method can be used to solve problems with rather high wave numbers.
|
|
4. |
- Arjmand, Doghonay, et al.
(författare)
-
A time dependent approach for removing the cell boundary error in elliptic homogenization problems
- 2016
-
Ingår i: Journal of Computational Physics. - : Elsevier. - 0021-9991 .- 1090-2716. ; 314, s. 206-227
-
Tidskriftsartikel (refereegranskat)abstract
- This paper concerns the cell-boundary error present in multiscale algorithms for elliptic homogenization problems. Typical multiscale methods have two essential components: a macro and a micro model. The micro model is used to upscale parameter values which are missing in the macro model. To solve the micro model, boundary conditions are required on the boundary of the microscopic domain. Imposing a naive boundary condition leads to O(epsilon/eta) error in the computation, where epsilon is the size of the microscopic variations in the media and eta is the size of the micro-domain. The removal of this error in modern multiscale algorithms still remains an important open problem. In this paper, we present a time-dependent approach which is general in terms of dimension. We provide a theorem which shows that we have arbitrarily high order convergence rates in terms of epsilon/eta in the periodic setting. Additionally, we present numerical evidence showing that the method improves the O(epsilon/eta) error to O(epsilon) in general non-periodic media.
|
|
5. |
- Arjmand, Doghonay, 1987-
(författare)
-
Analysis and Applications of Heterogeneous Multiscale Methods for Multiscale Partial Differential Equations
- 2015
-
Doktorsavhandling (övrigt vetenskapligt/konstnärligt)abstract
- This thesis centers on the development and analysis of numerical multiscale methods for multiscale problems arising in steady heat conduction, heat transfer and wave propagation in heterogeneous media. In a multiscale problem several scales interact with each other to form a system which has variations over a wide range of scales. A direct numerical simulation of such problems requires resolving the small scales over a computational domain, typically much larger than the microscopic scales. This demands a tremendous computational cost. We develop and analyse multiscale methods based on the heterogeneous multiscale methods (HMM) framework, which captures the macroscopic variations in the solution at a cost much lower than traditional numerical recipes. HMM assumes that there is a macro and a micro model which describes the problem. The micro model is accurate but computationally expensive to solve. The macro model is inexpensive but incomplete as it lacks certain parameter values. These are upscaled by solving the micro model locally in small parts of the domain. The accuracy of the method is then linked to how accurately this upscaling procedure captures the right macroscopic effects. In this thesis we analyse the upscaling error of existing multiscale methods and also propose a micro model which significantly reduces the upscaling error invarious settings. In papers I and IV we give an analysis of a finite difference HMM (FD-HMM) for approximating the effective solutions of multiscale wave equations over long time scales. In particular, we consider time scales T^ε = O(ε−k ), k =1, 2, where ε represents the size of the microstructures in the medium. In this setting, waves exhibit non-trivial behaviour which do not appear over short time scales. We use new analytical tools to prove that the FD-HMM accurately captures the long time effects. We first, in Paper I, consider T^ε =O(ε−2 ) and analyze the accuracy of FD-HMM in a one-dimensional periodicsetting. The core analytical ideas are quasi-polynomial solutions of periodic problems and local time averages of solutions of periodic wave equations.The analysis naturally reveals the role of consistency in HMM for high order approximation of effective quantities over long time scales. Next, in paperIV, we consider T^ε = O(ε−1 ) and use the tools in a multi-dimensional settingto analyze the accuracy of the FD-HMM in locally-periodic media where fast and slow variations are allowed at the same time. Moreover, in papers II and III we propose new multiscale methods which substantially improve the upscaling error in multiscale elliptic, parabolic and hyperbolic partial differential equations. In paper II we first propose a FD-HMM for solving elliptic homogenization problems. The strategy is to use the wave equation as the micro model even if the macro problem is of elliptic type. Next in paper III, we use this idea in a finite element HMM setting and generalize the approach to parabolic and hyperbolic problems. In a spatially fully discrete a priori error analysis we prove that the upscaling error can be made arbitrarily small for periodic media, even if we do not know the exact period of the oscillations in the media.
|
|
6. |
- Arjmand, Doghonay
(författare)
-
Analysis and Applications of the Heterogeneous Multiscale Methods for Multiscale Elliptic and Hyperbolic Partial Differential Equations
- 2013
-
Licentiatavhandling (övrigt vetenskapligt/konstnärligt)abstract
- This thesis concerns the applications and analysis of the Heterogeneous Multiscale methods (HMM) for Multiscale Elliptic and Hyperbolic Partial Differential Equations. We have gathered the main contributions in two papers.The first paper deals with the cell-boundary error which is present in multi-scale algorithms for elliptic homogenization problems. Typical multi-scale methods have two essential components: a macro and a micro model. The micro model is used to upscale parameter values which are missing in the macro model. Solving the micro model requires, on the other hand, imposing boundary conditions on the boundary of the microscopic domain. Imposing a naive boundary condition leads to $O(\varepsilon/\eta)$ error in the computation, where $\varepsilon$ is the size of the microscopic variations in the media and $\eta$ is the size of the micro-domain. Until now, strategies were proposed to improve the convergence rate up to fourth-order in $\varepsilon/\eta$ at best. However, the removal of this error in multi-scale algorithms still remains an important open problem. In this paper, we present an approach with a time-dependent model which is general in terms of dimension. With this approach we are able to obtain $O((\varepsilon/\eta)^q)$ and $O((\varepsilon/\eta)^q + \eta^p)$ convergence rates in periodic and locally-periodic media respectively, where $p,q$ can be chosen arbitrarily large. In the second paper, we analyze a multi-scale method developed under the Heterogeneous Multi-Scale Methods (HMM) framework for numerical approximation of wave propagation problems in periodic media. In particular, we are interested in the long time $O(\varepsilon^{-2})$ wave propagation. In the method, the microscopic model uses the macro solutions as initial data. In short-time wave propagation problems a linear interpolant of the macro variables can be used as the initial data for the micro-model. However, in long-time multi-scale wave problems the linear data does not suffice and one has to use a third-degree interpolant of the coarse data to capture the $O(1)$ dispersive effects apperaing in the long time. In this paper, we prove that through using an initial data consistent with the current macro state, HMM captures this dispersive effects up to any desired order of accuracy in terms of $\varepsilon/\eta$. We use two new ideas, namely quasi-polynomial solutions of periodic problems and local time averages of solutions of periodic hyperbolic PDEs. As a byproduct, these ideas naturally reveal the role of consistency for high accuracy approximation of homogenized quantities.
|
|
7. |
- Arjmand, Doghonay, 1987-, et al.
(författare)
-
Analysis of heterogeneous multiscale methods for long time wave propagation problems
- 2014
-
Ingår i: Multiscale Modeling & simulation. - : Society for Industrial & Applied Mathematics (SIAM). - 1540-3459 .- 1540-3467. ; 12:3, s. 1135-1166
-
Tidskriftsartikel (refereegranskat)abstract
- In this paper, we analyze a multiscale method developed under the heterogeneous multiscale method (HMM) framework for numerical approximation of multiscale wave propagation problems in periodic media. In particular, we are interested in the long time O(epsilon(-2)) wave propagation, where e represents the size of the microscopic variations in the media. In large time scales, the solutions of multiscale wave equations exhibit O(1) dispersive effects which are not observed in short time scales. A typical HMM has two main components: a macromodel and a micromodel. The macromodel is incomplete and lacks a set of local data. In the setting of multiscale PDEs, one has to solve for the full oscillatory problem over local microscopic domains of size eta = O(epsilon) to upscale the parameter values which are missing in the macroscopic model. In this paper, we prove that if the microproblems are consistent with the macroscopic solutions, the HMM approximates the unknown parameter values in the macromodel up to any desired order of accuracy in terms of epsilon/eta..
|
|
8. |
- Arjmand, Doghonay, 1987-, et al.
(författare)
-
Analysis of HMM for Long Time Multiscale Wave Propagation Problems in Locally-Periodic Media
-
Annan publikation (övrigt vetenskapligt/konstnärligt)abstract
- Multiscale wave propagation problems are difficult to solve numerically due to the interaction of different scales inherent in the problem. Extracting information about the average behaviour of the system requires resolving small scales in the problem. This leads to a tremendous computational burden if the size of microscopic variations are much smaller than the size of scales of interest. Heterogeneous multiscale methods (HMM) is a tool to avoid resolving the small scales everywhere. Nevertheless, it approximates the average part of the solution by upscaling the microscopic information on a small part of the domain. This leads to a substantial improvement in the computational cost. In this article, we analyze an HMM-based numerical method which approximates the long time behaviour of multiscale wave equations. In particular, we consider theoretically challenging case of locally-periodic media where fast and slow variations are allowed at the same time. We are interested in the long time regime (T=O(e^{-1})), where e represents the wavelength of the fast variations in themedia. We first use asymptotic expansions to derive effective equations describing the long time effects of the multiscale waves in multi-dimensional locally-periodic media. We then show that HMM captures these non-trivial long time eects. All the theoretical statements are general in terms of dimension. Two dimensional numericale xamples are considered to support our theoretical arguments
|
|
9. |
- Arjmand, Doghonay, et al.
(författare)
-
Estimates for the upscaling error in heterogeneous multiscale methods for wave propagation problems in locally periodic media
- 2017
-
Ingår i: Multiscale Modeling & simulation. - : Society for Industrial and Applied Mathematics. - 1540-3459 .- 1540-3467. ; 15:2, s. 948-976
-
Tidskriftsartikel (refereegranskat)abstract
- This paper concerns the analysis of a multiscale method for wave propagation problems in microscopically nonhomogeneous media. A direct numerical approximation of such problems is prohibitively expensive as it requires resolving the microscopic variations over a much larger physical domain of interest. The heterogeneous multiscale method (HMM) is an efficient framework to approximate the solutions of multiscale problems. In the HMM, one assumes an incomplete macroscopic model which is coupled to a known but expensive microscopic model. The micromodel is solved only locally to upscale the parameter values which are missing in the macro model. The resulting macroscopic model can then be solved at a cost independent of the small scales in the problem. In general, the accuracy of the HMM is related to how good the upscaling step approximates the right macroscopic quantities. The analysis of the method that we consider here was previously addressed only in purely periodic media, although the method itself is numerically shown to be applicable to more general settings. In the present study, we consider a more realistic setting by assuming a locally periodic medium where slow and fast variations are allowed at the same time. We then prove that the HMM captures the right macroscopic effects. The generality of the tools and ideas in the analysis allows us to establish convergence rates in a multidimensional setting. The theoretical findings here imply an improved convergence rate in one dimension, which also justifies the numerical observations from our earlier study.
|
|
10. |
- Benamou, J. D., et al.
(författare)
-
Numerical microlocal analysis of harmonic wavefields
- 2004
-
Ingår i: Journal of Computational Physics. - : Elsevier BV. - 0021-9991 .- 1090-2716. ; 199:2, s. 717-741
-
Tidskriftsartikel (refereegranskat)abstract
- We present and test a numerical method which, given an analytical or numerical solution of the Helmholtz equation in a neighborhood of a fixed observation point and assuming that the geometrical optics approximation is relevant, determines at this point the number of crossing rays and computes their directions and associated complex amplitudes.
|
|