Please wait ... |

Refine search result

CiteExportLink to result list
http://kth.diva-portal.org/smash/resultList.jsf?query=&language=en&searchType=SIMPLE&noOfRows=50&sortOrder=author_sort_asc&sortOrder2=title_sort_asc&onlyFullText=false&sf=all&aq=%5B%5B%7B%22personId%22%3A%22u1na5bn9%22%7D%5D%5D&aqe=%5B%5D&aq2=%5B%5B%5D%5D&af=%5B%5D $(function(){PrimeFaces.cw("InputTextarea","widget_formSmash_upper_j_idt483_recordPermLink",{id:"formSmash:upper:j_idt483:recordPermLink",widgetVar:"widget_formSmash_upper_j_idt483_recordPermLink",autoResize:true});}); $(function(){PrimeFaces.cw("OverlayPanel","widget_formSmash_upper_j_idt483_j_idt485",{id:"formSmash:upper:j_idt483:j_idt485",widgetVar:"widget_formSmash_upper_j_idt483_j_idt485",target:"formSmash:upper:j_idt483:permLink",showEffect:"blind",hideEffect:"fade",my:"right top",at:"right bottom",showCloseIcon:true});});

Permanent link

Cite

Citation styleapa ieee modern-language-association-8th-edition vancouver Other style $(function(){PrimeFaces.cw("SelectOneMenu","widget_formSmash_upper_j_idt501",{id:"formSmash:upper:j_idt501",widgetVar:"widget_formSmash_upper_j_idt501",behaviors:{change:function(ext) {PrimeFaces.ab({s:"formSmash:upper:j_idt501",e:"change",f:"formSmash",p:"formSmash:upper:j_idt501",u:"formSmash:upper:otherStyle"},ext);}}});});

- apa
- ieee
- modern-language-association-8th-edition
- vancouver
- Other style

Languagede-DE en-GB en-US fi-FI nn-NO nn-NB sv-SE Other locale $(function(){PrimeFaces.cw("SelectOneMenu","widget_formSmash_upper_j_idt512",{id:"formSmash:upper:j_idt512",widgetVar:"widget_formSmash_upper_j_idt512",behaviors:{change:function(ext) {PrimeFaces.ab({s:"formSmash:upper:j_idt512",e:"change",f:"formSmash",p:"formSmash:upper:j_idt512",u:"formSmash:upper:otherLanguage"},ext);}}});});

- de-DE
- en-GB
- en-US
- fi-FI
- nn-NO
- nn-NB
- sv-SE
- Other locale

Output formathtml text asciidoc rtf $(function(){PrimeFaces.cw("SelectOneMenu","widget_formSmash_upper_j_idt522",{id:"formSmash:upper:j_idt522",widgetVar:"widget_formSmash_upper_j_idt522"});});

- html
- text
- asciidoc
- rtf

Rows per page

- 5
- 10
- 20
- 50
- 100
- 250

Sort

- Standard (Relevance)
- Author A-Ö
- Author Ö-A
- Title A-Ö
- Title Ö-A
- Publication type A-Ö
- Publication type Ö-A
- Issued (Oldest first)
- Issued (Newest first)
- Created (Oldest first)
- Created (Newest first)
- Last updated (Oldest first)
- Last updated (Newest first)
- Disputation date (earliest first)
- Disputation date (latest first)

- Standard (Relevance)
- Author A-Ö
- Author Ö-A
- Title A-Ö
- Title Ö-A
- Publication type A-Ö
- Publication type Ö-A
- Issued (Oldest first)
- Issued (Newest first)
- Created (Oldest first)
- Created (Newest first)
- Last updated (Oldest first)
- Last updated (Newest first)
- Disputation date (earliest first)
- Disputation date (latest first)

Select

The maximal number of hits you can export is 250. When you want to export more records please use the Create feeds function.

1. Arjmand, Doghonay et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_0_j_idt588",{id:"formSmash:items:resultList:0:j_idt588",widgetVar:"widget_formSmash_items_resultList_0_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:0:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Runborg, OlofKTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:0:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); A time dependent approach for removing the cell boundary error in elliptic homogenization problems2016In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 314, p. 206-227Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_0_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:0:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_0_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); 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.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:0:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 2. Arjmand, Doghonay PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_1_j_idt585",{id:"formSmash:items:resultList:1:j_idt585",widgetVar:"widget_formSmash_items_resultList_1_j_idt585",onLabel:"Arjmand, Doghonay ",offLabel:"Arjmand, Doghonay ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_1_j_idt588",{id:"formSmash:items:resultList:1:j_idt588",widgetVar:"widget_formSmash_items_resultList_1_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA. KTH, Centres, SeRC - Swedish e-Science Research Centre.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:1:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Runborg, OlofKTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA. KTH, Centres, SeRC - Swedish e-Science Research Centre.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:1:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Analysis of heterogeneous multiscale methods for long time wave propagation problems2014In: Multiscale Modeling & simulation, ISSN 1540-3459, E-ISSN 1540-3467, Vol. 12, no 3, p. 1135-1166Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_1_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:1:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_1_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); 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..

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:1:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 3. Arjmand, Doghonay PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_2_j_idt585",{id:"formSmash:items:resultList:2:j_idt585",widgetVar:"widget_formSmash_items_resultList_2_j_idt585",onLabel:"Arjmand, Doghonay ",offLabel:"Arjmand, Doghonay ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_2_j_idt588",{id:"formSmash:items:resultList:2:j_idt588",widgetVar:"widget_formSmash_items_resultList_2_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:2:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Runborg, OlofKTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:2:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Analysis of HMM for Long Time Multiscale Wave Propagation Problems in Locally-Periodic MediaManuscript (preprint) (Other academic)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_2_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:2:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_2_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); 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

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:2:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 4. Arjmand, Doghonay et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_3_j_idt588",{id:"formSmash:items:resultList:3:j_idt588",widgetVar:"widget_formSmash_items_resultList_3_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:3:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Runborg, OlofKTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA. KTH, Centres, SeRC - Swedish e-Science Research Centre.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:3:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Estimates for the upscaling error in heterogeneous multiscale methods for wave propagation problems in locally periodic media2017In: Multiscale Modeling & simulation, ISSN 1540-3459, E-ISSN 1540-3467, Vol. 15, no 2, p. 948-976Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_3_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:3:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_3_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); 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.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:3:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 5. Benamou, J. D. et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_4_j_idt588",{id:"formSmash:items:resultList:4:j_idt588",widgetVar:"widget_formSmash_items_resultList_4_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:4:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Collino, F.Runborg, OlofKTH, Superseded Departments, Numerical Analysis and Computer Science, NADA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:4:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Numerical microlocal analysis of harmonic wavefields2004In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 199, no 2, p. 717-741Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_4_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:4:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_4_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); 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.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:4:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 6. Castella, F. et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_5_j_idt588",{id:"formSmash:items:resultList:5:j_idt588",widgetVar:"widget_formSmash_items_resultList_5_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:5:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Perthame, B.Runborg, OlofPrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:5:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); High frequency limit of the Helmholtz equation. II. Source on a general smooth manifold2002In: Communications in Partial Differential Equations, ISSN 0360-5302, E-ISSN 1532-4133, Vol. 27, no 04-mar, p. 607-651Article in journal (Refereed)7. Daubechies, I. et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_6_j_idt588",{id:"formSmash:items:resultList:6:j_idt588",widgetVar:"widget_formSmash_items_resultList_6_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:6:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Runborg, OlofKTH, Superseded Departments, Numerical Analysis and Computer Science, NADA.Sweldens, W.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:6:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Normal multiresolution approximation of curves2004In: Constructive approximation, ISSN 0176-4276, E-ISSN 1432-0940, Vol. 20, no 3, p. 399-463Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_6_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:6:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_6_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); A multiresolution analysis of a curve is normal if each wavelet detail vector with respect to a certain subdivision scheme lies in the local normal direction. In this paper we study properties such as regularity, convergence, and stability of a normal multiresolution analysis. In particular, we show that these properties critically depend on the underlying subdivision scheme and that, in general, the convergence of normal multiresolution approximations equals the convergence of the underlying subdivision scheme.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:6:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 8. Daubechies, Ingrid et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_7_j_idt588",{id:"formSmash:items:resultList:7:j_idt588",widgetVar:"widget_formSmash_items_resultList_7_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:7:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Runborg, OlofKTH, School of Computer Science and Communication (CSC), Numerical Analysis and Computer Science, NADA.Zou, JingPrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:7:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); A sparse spectral method for homogenization multiscale problems2007In: Multiscale Modeling & simulation, ISSN 1540-3459, E-ISSN 1540-3467, Vol. 6, no 3, p. 711-740Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_7_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:7:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_7_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); We develop a new sparse spectral method, in which the fast Fourier transform (FFT) is replaced by RAlSFA (randomized algorithm of sparse Fourier analysis); this is a sublinear randomized algorithm that takes time O(B log N) to recover a B-term Fourier representation for a signal of length N, where we assume B << N. To illustrate its potential, we consider the parabolic homogenization problem with a characteristic. ne scale size epsilon. For fixed tolerance the sparse method has a computational cost of O(vertical bar log epsilon vertical bar) per time step, whereas standard methods cost at least O(epsilon(-1)). We present a theoretical analysis as well as numerical results; they show the advantage of the new method in speed over the traditional spectral methods when epsilon is very small. We also show some ways to extend the methods to hyperbolic and elliptic problems.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:7:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 9. Di, Yana et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_8_j_idt588",{id:"formSmash:items:resultList:8:j_idt588",widgetVar:"widget_formSmash_items_resultList_8_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:8:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Popovic, JelenaKTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.Runborg, OlofKTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA. KTH, Centres, SeRC - Swedish e-Science Research Centre.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:8:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); AN ADAPTIVE FAST INTERFACE TRACKING METHOD2015In: Journal of Computational Mathematics, ISSN 0254-9409, E-ISSN 1991-7139, Vol. 33, no 6, p. 576-586Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_8_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:8:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_8_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); An adaptive numerical scheme is developed for the propagation of an interface in a velocity field based on the fast interface tracking method proposed in [2]. A multiresolution stategy to represent the interface instead of point values, allows local grid refinement while controlling the approximation error on the interface. For time integration, we use an explicit Runge-Kutta scheme of second-order with a multiscale time step, which takes longer time steps for finer spatial scales. The implementation of the algorithm uses a dynamic tree data structure to represent data in the computer memory. We briefly review first the main algorithm, describe the essential data structures, highlight the adaptive scheme, and illustrate the computational efficiency by some numerical examples.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:8:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 10. Engquist, Björn PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_9_j_idt585",{id:"formSmash:items:resultList:9:j_idt585",widgetVar:"widget_formSmash_items_resultList_9_j_idt585",onLabel:"Engquist, Björn ",offLabel:"Engquist, Björn ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_9_j_idt588",{id:"formSmash:items:resultList:9:j_idt588",widgetVar:"widget_formSmash_items_resultList_9_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA (closed 2012-06-30).PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:9:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Holst, HenrikKTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA (closed 2012-06-30).Runborg, OlofKTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA (closed 2012-06-30).PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:9:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Analysis of HMM for One Dimensional Wave Propagation Problems Over Long Time2011Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_9_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:9:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_9_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); Multiscale problems are computationally costly to solve by direct simulation because the smallest scales must be represented over a domain determined by the largest scales of the problem. We have developed and analyzed new numerical methods for multiscale wave propagation following the framework of the heterogeneous multiscale method. The numerical methods couple simulations on macro- and microscales for problems with rapidly fluctuating material coefficients. The computational complexity of the new method is significantly lower than that of traditional techniques. We focus on HMM approximation applied to long time integration of one-dimensional wave propagation problems in both periodic and non-periodic medium and show that the dispersive effect that appear after long time is fully captured.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:9:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 11. Engquist, Björn PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_10_j_idt585",{id:"formSmash:items:resultList:10:j_idt585",widgetVar:"widget_formSmash_items_resultList_10_j_idt585",onLabel:"Engquist, Björn ",offLabel:"Engquist, Björn ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_10_j_idt588",{id:"formSmash:items:resultList:10:j_idt588",widgetVar:"widget_formSmash_items_resultList_10_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); University of Texas Austin.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:10:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Holst, HenrikKTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA.Runborg, OlofKTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:10:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Multiscale Methods for One Dimensional Wave Propagation with High Frequency Initial Data2011Report (Other academic)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_10_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:10:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_10_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); High frequency wave propagation problems are computationally costly to solve by traditional techniques because the short wavelength must be well represented over a domain determined by the largest scales of the problem. We have developed and analyzed a new numerical method for high frequency wave propagation in the framework of heterogeneous multiscale methods, closely related to the analytical method of geometrical optics. The numerical method couples simulations on macro- and micro-scales for problems with highly oscillatory initial data. The method has a computational complexity essentially independent of the wavelength. We give one numerical example with a sharp but regular jump in velocity on the microscopic scale for which geometrical optics fails but our HMM gives correct results. We briefly discuss how the method can be extended to higher dimensional problems.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:10:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 12. Engquist, Björn PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_11_j_idt585",{id:"formSmash:items:resultList:11:j_idt585",widgetVar:"widget_formSmash_items_resultList_11_j_idt585",onLabel:"Engquist, Björn ",offLabel:"Engquist, Björn ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_11_j_idt588",{id:"formSmash:items:resultList:11:j_idt588",widgetVar:"widget_formSmash_items_resultList_11_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:11:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Holst, HenrikKTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA.Runborg, OlofKTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:11:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Multiscale methods for the wave equation2007In: PAMM · Proc. Appl. Math. Mech. 7, 2007, p. 1140903-1140904Conference paper (Other academic)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_11_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:11:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_11_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); We consider the wave equation in a medium with a rapidly varying speed of propagation. We construct a multiscale schemebased on the heterogeneous multiscale method, which can compute the correct coarse behavior of wave pulses traveling in themedium, at a computational cost essentially independent of the size of the small scale variations. This is verified by theoreticalresults and numerical examples.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:11:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 13. Engquist, Björn PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_12_j_idt585",{id:"formSmash:items:resultList:12:j_idt585",widgetVar:"widget_formSmash_items_resultList_12_j_idt585",onLabel:"Engquist, Björn ",offLabel:"Engquist, Björn ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_12_j_idt588",{id:"formSmash:items:resultList:12:j_idt588",widgetVar:"widget_formSmash_items_resultList_12_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:12:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Holst, HenrikKTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA.Runborg, OlofKTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:12:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Multi-scale methods for wave propagation in heterogeneous media2011In: Communications in Mathematical Sciences, ISSN 1539-6746, E-ISSN 1945-0796, Vol. 9, no 1, p. 33-56Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_12_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:12:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_12_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); Multi-scale wave propagation problems are computationally costly to solve by traditional techniques because the smallest scales must be represented over a domain determined by the largest scales of the problem. We have developed and analyzed new numerical methods for multi-scale wave propagation in the framework of heterogeneous multi-scale method. The numerical methods couple simulations on macro-and micro-scales for problems with rapidly oscillating coefficients. We show that the complexity of the new method is significantly lower than that of traditional techniques with a computational cost that is essentially independent of the micro-scale. A convergence proof is given and numerical results are presented for periodic problems in one, two, and three dimensions. The method is also successfully applied to non-periodic problems and for long time integration where dispersive effects occur.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:12:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 14. Engquist, Björn PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_13_j_idt585",{id:"formSmash:items:resultList:13:j_idt585",widgetVar:"widget_formSmash_items_resultList_13_j_idt585",onLabel:"Engquist, Björn ",offLabel:"Engquist, Björn ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_13_j_idt588",{id:"formSmash:items:resultList:13:j_idt588",widgetVar:"widget_formSmash_items_resultList_13_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA (closed 2012-06-30).PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:13:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Holst, HenrikKTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA (closed 2012-06-30).Runborg, OlofKTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA (closed 2012-06-30).PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:13:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Multiscale Methods for Wave Propagation in Heterogeneous Media Over Long Time2012In: Numerical Analysis of Multiscale Computations: Proceedings Workshop on Numerical Analysis and Multiscale Computations, Banff Int Res Stn, Banff, CANADA, DEC 06-11, 2009 / [ed] Björn Engquist, Olof Runborg, Yen-Hsi R. Tsai, Springer Verlag , 2012, p. 167-186Chapter in book (Other academic)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_13_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:13:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_13_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); Multiscale wave propagation problems are computationally costly to solve by traditional techniques because the smallest scales must be represented over a domain determined by the largest scales of the problem. We have developed and analyzed new numerical methods for multiscale wave propagation in the framework of the heterogeneous multiscale method (HMM). The numerical methods couple simulations on macro- and microscales for problems with rapidly oscillating coefficients. The complexity of the new method is significantly lower than that of traditional techniques with a computational cost that is essentially independent of the smallest scale, when computing solutions at a fixed time and accuracy. We show numerical examples of the HMM applied to long time integration of wave propagation problems in both periodic and non-periodic medium. In both cases our HMM accurately captures the dispersive effects that occur. We also give a stability proof for the HMM, when it is applied to long time wave propagation problems.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:13:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 15. Engquist, Björn PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_14_j_idt585",{id:"formSmash:items:resultList:14:j_idt585",widgetVar:"widget_formSmash_items_resultList_14_j_idt585",onLabel:"Engquist, Björn ",offLabel:"Engquist, Björn ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_14_j_idt588",{id:"formSmash:items:resultList:14:j_idt588",widgetVar:"widget_formSmash_items_resultList_14_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); Department of Mathematics and Institute for Computational Engineering and Sciences, The University of Texas at Austin, USA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:14:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Häggblad, JonKTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA.Runborg, OlofKTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:14:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Numerical subgrid scale models for the Yee schemeManuscript (preprint) (Other academic)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_14_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:14:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_14_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); The Yee scheme is a very common and practical algorithm for the simulation of wave propagation on uniform grids. We develop numerical subgrid scale models in order to incorporate effects of obstacles and holes that are smaller than the grid spacing. The models are based on pre-computing at the microscale, and are thus including the effect of the detailed small scale shape. Numerical examples in 1D, 2D and 3D are given.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:14:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 16. Engquist, Björn PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_15_j_idt585",{id:"formSmash:items:resultList:15:j_idt585",widgetVar:"widget_formSmash_items_resultList_15_j_idt585",onLabel:"Engquist, Björn ",offLabel:"Engquist, Björn ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_15_j_idt588",{id:"formSmash:items:resultList:15:j_idt588",widgetVar:"widget_formSmash_items_resultList_15_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); Univ Texas Austin.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:15:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Häggblad, JonKTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA (closed 2012-06-30).Runborg, OlofKTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA (closed 2012-06-30).PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:15:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); On Energy Preserving Consistent Boundary Conditions for the Yee Scheme in 2D2012In: BIT Numerical Mathematics, ISSN 0006-3835, E-ISSN 1572-9125, Vol. 52, no 3, p. 615-637Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_15_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:15:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_15_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); The Yee scheme is one of the most popular methods for electromagnetic wave propagation. A main advantage is the structured staggered grid, making it simple and efficient on modern computer architectures. A downside to this is the difficulty in approximating oblique boundaries, having to resort to staircase approximations. In this paper we present a method to improve the boundary treatment in two dimensions by, starting from a staircase approximation, modifying the coefficients of the update stencil so that we can obtain a consistent approximation while preserving the energy conservation, structure and the optimal CFL-condition of the original Yee scheme. We prove this in L_2 and verify it by numerical experiments.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:15:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 17. Engquist, Björn PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_16_j_idt585",{id:"formSmash:items:resultList:16:j_idt585",widgetVar:"widget_formSmash_items_resultList_16_j_idt585",onLabel:"Engquist, Björn ",offLabel:"Engquist, Björn ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_16_j_idt588",{id:"formSmash:items:resultList:16:j_idt588",widgetVar:"widget_formSmash_items_resultList_16_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); Department of Mathematics and Institute for Computational Engineering and Sciences, The University of Texas at Austin, USA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:16:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Häggblad, JonKTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA.Runborg, OlofKTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA.Tornberg, Anna-KarinKTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:16:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); On Consistent Boundary Conditions for the Yee Scheme in 3DManuscript (preprint) (Other academic)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_16_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:16:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_16_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); The standard staircase approximation of curved boundaries in the Yee scheme is inconsistent. Consistency can however be achieved by modifying the algorithm close to the boundary. We consider a technique to consistently model curved boundaries where the coefficients of the update stencil is modified, thus preserving the Yee structure. The method has previously been successfully applied to acoustics in two and three dimension, as well as electromagnetics in two dimensions. In this paper we generalize to electromagnetics in three dimensions. Unlike in previous cases there is a non-zero divergence growth along the boundary that needs to be projected away. We study the convergence and provide numerical examples that demonstrates the improved accuracy.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:16:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 18. Engquist, Björn PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_17_j_idt585",{id:"formSmash:items:resultList:17:j_idt585",widgetVar:"widget_formSmash_items_resultList_17_j_idt585",onLabel:"Engquist, Björn ",offLabel:"Engquist, Björn ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_17_j_idt588",{id:"formSmash:items:resultList:17:j_idt588",widgetVar:"widget_formSmash_items_resultList_17_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:17:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Lotstedt, P.Runborg, O.KTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:17:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Multiscale Modeling and Simulation in Science2009Book (Refereed)19. Engquist, Björn PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_18_j_idt585",{id:"formSmash:items:resultList:18:j_idt585",widgetVar:"widget_formSmash_items_resultList_18_j_idt585",onLabel:"Engquist, Björn ",offLabel:"Engquist, Björn ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_18_j_idt588",{id:"formSmash:items:resultList:18:j_idt588",widgetVar:"widget_formSmash_items_resultList_18_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA (closed 2012-06-30).PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:18:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Lotstedt, P.Runborg, OlofKTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA (closed 2012-06-30).PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:18:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Multiscale methods in science and engineering2005Book (Refereed)20. Engquist, Björn PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_19_j_idt585",{id:"formSmash:items:resultList:19:j_idt585",widgetVar:"widget_formSmash_items_resultList_19_j_idt585",onLabel:"Engquist, Björn ",offLabel:"Engquist, Björn ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_19_j_idt588",{id:"formSmash:items:resultList:19:j_idt588",widgetVar:"widget_formSmash_items_resultList_19_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); Department of Mathematics, The University of Texas at Austin, Austin, USA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:19:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Runborg, OlofKTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:19:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Wavelet based numerical homogenization2009In: Highly Oscillatory Problems / [ed] Bjorn Engquist, Athanasios Fokas, Ernst Hairer, Arieh Iserles, Cambridge University Press, 2009, p. 98-126Chapter in book (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_19_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:19:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_19_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); We consider multiscale differential equations in which the operator varies rapidly over fine scales. Direct numerical simulation methods need to resolve the small scales and they therefore become very expensive for such problems when the computational domain is large. Inspired by classical homogenization theory, we describe a numerical procedure for homogenization, which starts from a fine discretization of a multiscale differential equation, and computes a discrete coarse grid operator which incorporates the influence of finer scales. In this procedure the discrete operator is represented in a wavelet space, projected onto a coarser subspace and approximated by a banded or block-banded matrix. This wavelet homogenization applies to a wider class of problems than classical homogenization. The projection procedure is general and we give a presentation of a framework in Hilbert spaces, which also applies to the differential equation directly. We show numerical results when the wavelet based homogenization technique is applied to discretizations of elliptic and hyperbolic equations, using different approximation strategies for the coarse grid operator.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:19:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 21. Engquist, Björn PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_20_j_idt585",{id:"formSmash:items:resultList:20:j_idt585",widgetVar:"widget_formSmash_items_resultList_20_j_idt585",onLabel:"Engquist, Björn ",offLabel:"Engquist, Björn ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_20_j_idt588",{id:"formSmash:items:resultList:20:j_idt588",widgetVar:"widget_formSmash_items_resultList_20_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:20:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Runborg, OlofKTH, Superseded Departments, Numerical Analysis and Computer Science, NADA.Tornberg, Anna KarinPrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:20:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); High-frequency wave propagation by the segment projection method2002In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 178, no 2, p. 373-390Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_20_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:20:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_20_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); Geometrical optics is a standard technique used for the approximation of high-frequency wave propagation. Computational methods based on partial differential equations instead of the traditional ray tracing have recently been applied to geometrical optics. These new methods have a number of advantages but typically exhibit difficulties with linear superposition of waves. In this paper we introduce a new partial differential technique based on the segment projection method in phase space. The superposition problem is perfectly resolved and so is the problem of computing amplitudes in the neighborhood of caustics. The computational complexity is of the same order as that of ray tracing. The new algorithm is described and a number of computational examples are given. including a simulation of waveguides.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:20:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 22. Engquist, Björn et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_21_j_idt588",{id:"formSmash:items:resultList:21:j_idt588",widgetVar:"widget_formSmash_items_resultList_21_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:21:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Runborg, OlofKTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA. KTH, Centres, SeRC - Swedish e-Science Research Centre.Tsai, Y. -HR.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:21:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Preface2012In: Workshop on Numerical Analysis and Multiscale Computations, 2009, Springer Verlag , 2012Conference paper (Refereed)23. Gosse, L. et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_22_j_idt588",{id:"formSmash:items:resultList:22:j_idt588",widgetVar:"widget_formSmash_items_resultList_22_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:22:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Runborg, OlofKTH, School of Computer Science and Communication (CSC), Numerical Analysis and Computer Science, NADA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:22:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Finite moment problems and applications to multiphase computations in geometric optics2005In: Communications in Mathematical Sciences, ISSN 1539-6746, E-ISSN 1945-0796, Vol. 3, no 3, p. 373-392Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_22_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:22:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_22_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); Recovering a function out of a finite number of moments is generally an ill-posed inverse problem. We focus on two special cases arising from applications to multiphase geometric optics computations where this problem can be carried out in a restricted class of given densities. More precisely, we present a simple algorithm for the inversion of Markov's moment problem which appears in the treatment of Brenier and Corrias' K-multibranch solutions and study Stieltje's algorithm in order to process moment systems arising from a Wigner analysis. Numerical results are provided for moderately intricate wave-fields.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:22:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 24. Gosse, L. et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_23_j_idt588",{id:"formSmash:items:resultList:23:j_idt588",widgetVar:"widget_formSmash_items_resultList_23_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:23:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Runborg, OlofKTH, School of Computer Science and Communication (CSC), Numerical Analysis and Computer Science, NADA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:23:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Resolution of the finite Markov moment problem2005In: Comptes rendus. Mathematique, ISSN 1631-073X, E-ISSN 1778-3569, Vol. 341, no 12, p. 775-780Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_23_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:23:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_23_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); We expose in full detail a constructive procedure to invert the so-called 'finite Markov moment problem'. The proofs rely on the general theory of Toeplitz matrices together with the classical Newton's relations.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:23:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 25. Gosse, Laurent et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_24_j_idt588",{id:"formSmash:items:resultList:24:j_idt588",widgetVar:"widget_formSmash_items_resultList_24_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:24:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Runborg, OlofKTH, School of Computer Science and Communication (CSC), Numerical Analysis and Computer Science, NADA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:24:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Existence, uniqueness, and a constructive solution algorithm for a class of finite Markov moment problems2007In: SIAM Journal on Applied Mathematics, ISSN 0036-1399, E-ISSN 1095-712X, Vol. 68, no 6, p. 1618-1640Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_24_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:24:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_24_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); We consider a class of finite Markov moment problems with an arbitrary number of positive and negative branches. We show criteria for the existence and uniqueness of solutions, and we characterize in detail the nonunique solution families. Moreover, we present a constructive algorithm to solve the moment problems numerically and prove that the algorithm computes the right solution.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:24:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 26. Häggblad, Jon PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_25_j_idt585",{id:"formSmash:items:resultList:25:j_idt585",widgetVar:"widget_formSmash_items_resultList_25_j_idt585",onLabel:"Häggblad, Jon ",offLabel:"Häggblad, Jon ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_25_j_idt588",{id:"formSmash:items:resultList:25:j_idt588",widgetVar:"widget_formSmash_items_resultList_25_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA (closed 2012-06-30).PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:25:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Runborg, OlofKTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA (closed 2012-06-30). KTH, Centres, SeRC - Swedish e-Science Research Centre.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:25:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Accuracy of staircase approximations in finite-difference methods for wave propagation2014In: Numerische Mathematik, ISSN 0029-599X, E-ISSN 0945-3245, Vol. 128, no 4, p. 741-771Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_25_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:25:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_25_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); While a number of increasingly sophisticated numerical methods have been developed for time-dependent problems in electromagnetics, the Yee scheme is still widely used in the applied fields, mainly due to its simplicity and computational efficiency. A fundamental drawback of the method is the use of staircase boundary approximations, giving inconsistent results. Usually experience of numerical experiments provides guidance of the impact of these errors on the final simulation result. In this paper, we derive exact discrete solutions to the Yee scheme close to the staircase approximated boundary, enabling a detailed theoretical study of the amplitude, phase and frequency errors created. Furthermore, we show how evanescent waves of amplitude occur along the boundary. These characterize the inconsistencies observed in electromagnetic simulations and the locality of the waves explain why, in practice, the Yee scheme works as well as it does. The analysis is supported by detailed proofs and numerical examples.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:25:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 27. Jarlebring, Elias PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_26_j_idt585",{id:"formSmash:items:resultList:26:j_idt585",widgetVar:"widget_formSmash_items_resultList_26_j_idt585",onLabel:"Jarlebring, Elias ",offLabel:"Jarlebring, Elias ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_26_j_idt588",{id:"formSmash:items:resultList:26:j_idt588",widgetVar:"widget_formSmash_items_resultList_26_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA. KTH, Centres, SeRC - Swedish e-Science Research Centre.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:26:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Mele, GiampaoloKTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA. KTH, Centres, SeRC - Swedish e-Science Research Centre.Runborg, OlofKTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA. KTH, Centres, SeRC - Swedish e-Science Research Centre.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:26:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); The waveguide eigenvalue problem and the tensor infinite Arnoldi method2017In: SIAM Journal on Scientific Computing, ISSN 1064-8275, E-ISSN 1095-7197, Vol. 39, no 3, p. A1062-A1088Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_26_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:26:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_26_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); We present a new computational approach for a class of large-scale nonlinear eigenvalue problems (NEPs) that are nonlinear in the eigenvalue. The contribution of this paper is two fold. We derive a new iterative algorithm for NEPs, the tensor infinite Arnoldi method (TIAR), which is applicable to a general class of NEPs, and we show how to specialize the algorithm to a specific NEP: the waveguide eigenvalue problem. The waveguide eigenvalue problem arises from a finite-element discretization of a partial differential equation used in the study waves propagating in a periodic medium. The algorithm is successfully applied to accurately solve benchmark problems as well as complicated waveguides. We study the complexity of the specialized algorithm with respect to the number of iterations "m" and the size of the problem "n", both from a theoretical perspective and in practice. For the waveguide eigenvalue problem, we establish that the computationally dominating part of the algorithm has complexity O(nm^2+sqrt(n)m^3). Hence, the asymptotic complexity of TIAR applied to the waveguide eigenvalue problem, for n→ ∞, is the same as for Arnoldi’s method for standard eigenvalue problems.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:26:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 28. Kieri, Emil et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_27_j_idt588",{id:"formSmash:items:resultList:27:j_idt588",widgetVar:"widget_formSmash_items_resultList_27_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:27:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Kreiss, GunillaRunborg, OlofKTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA. KTH, Centres, SeRC - Swedish e-Science Research Centre.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:27:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Coupling of Gaussian Beam and Finite Difference Solvers for Semiclassical Schrodinger Equations2015In: Advances in Applied Mathematics and Mechanics, ISSN 2070-0733, E-ISSN 2075-1354, Vol. 7, no 6, p. 687-714Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_27_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:27:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_27_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); In the semiclassical regime, solutions to the time-dependent Schrodinger equation for molecular dynamics are highly oscillatory. The number of grid points required for resolving the oscillations may become very large even for simple model problems, making solution on a grid intractable. Asymptotic methods like Gaussian beams can resolve the oscillations with little effort and yield good approximations when the atomic nuclei are heavy and the potential is smooth. However, when the potential has variations on a small length-scale, quantum phenomena become important. Then asymptotic methods are less accurate. The two classes of methods perform well in different parameter regimes. This opens for hybrid methods, using Gaussian beams where we can and finite differences where we have to. We propose a new method for treating the coupling between the finite difference method and Gaussian beams. The new method reduces the needed amount of overlap regions considerably compared to previous methods, which improves the efficiency.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:27:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 29. Krishnan, J. et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_28_j_idt588",{id:"formSmash:items:resultList:28:j_idt588",widgetVar:"widget_formSmash_items_resultList_28_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:28:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Runborg, OlofKevrekidis, I. G.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:28:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Bifurcation analysis of nonlinear reaction-diffusion problems using wavelet-based reduction techniques2004In: Computers and Chemical Engineering, ISSN 0098-1354, E-ISSN 1873-4375, Vol. 28, no 4, p. 557-574Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_28_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:28:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_28_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); Using a computational method for numerical homogenization, we perform the coarse-scale bifurcation analysis of nonlinear reaction-diffusion problems in both uniform and spatially varying media. The method is based on wavelet decomposition and projection of the differential equation on coarse scale wavelet spaces. The approach is capable of capturing turning points and pitchfork bifurcations of sharp, front-like solutions at the coarse level.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:28:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 30. Liu, H. et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_29_j_idt588",{id:"formSmash:items:resultList:29:j_idt588",widgetVar:"widget_formSmash_items_resultList_29_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:29:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Runborg, OlofKTH, School of Engineering Sciences (SCI), Mathematics (Dept.). KTH, Centres, SeRC - Swedish e-Science Research Centre.Tanushev, N. M.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:29:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Sobolev and max norm error estimates for Gaussian beam superpositions2016In: Communications in Mathematical Sciences, ISSN 1539-6746, E-ISSN 1945-0796, Vol. 14, no 7, p. 2037-2072Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_29_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:29:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_29_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); This work is concerned with the accuracy of Gaussian beam superpositions, which are asymptotically valid high frequency solutions to linear hyperbolic partial differential equations and the Schrödinger equation. We derive Sobolev and max norms estimates for the difference between an exact solution and the corresponding Gaussian beam approximation, in terms of the short wavelength e. The estimates are performed for the scalar wave equation and the Schrödinger equation. Our result demonstrates that a Gaussian beam superposition with kth order beams converges to the exact solution as O(εk/2-s) in order s Sobolev norms. This result is valid in any number of spatial dimensions and it is unaffected by the presence of caustics in the solution. In max norm, we show that away from caustics the convergence rate is O(ε⌈k/2⌉) and away from the essential support of the solution, the convergence is spectral in ε. However, in the neighborhood of a caustic point we are only able to show the slower, and dimensional dependent, rate O(ε(k-n)/2) in n spatial dimensions.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:29:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 31. Liu, Hailiang et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_30_j_idt588",{id:"formSmash:items:resultList:30:j_idt588",widgetVar:"widget_formSmash_items_resultList_30_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:30:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Ralston, JamesRunborg, OlofKTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA. KTH, Centres, SeRC - Swedish e-Science Research Centre.Tanushev, Nicolay M.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:30:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Gaussian beam methods for the helmholtz equation2014In: SIAM Journal on Applied Mathematics, ISSN 0036-1399, E-ISSN 1095-712X, Vol. 74, no 3, p. 771-793Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_30_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:30:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_30_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); In this work we construct Gaussian beam approximations to solutions of the high frequency Helmholtz equation with a localized source. Under the assumption of nontrapping rays we show error estimates between the exact outgoing solution and Gaussian beams in terms of the wave number k, both for single beams and superposition of beams. The main result is that the relative local L-2 error in the beam approximations decay as k(-N/2) independent of dimension and presence of caustics for Nth order beams.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:30:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 32. Liu, Hailiang et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_31_j_idt588",{id:"formSmash:items:resultList:31:j_idt588",widgetVar:"widget_formSmash_items_resultList_31_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:31:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Runborg, OlofKTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA. KTH, Centres, SeRC - Swedish e-Science Research Centre.Tanushev, Nicolay M.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:31:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Error estimates for Gaussian beam superpositions2013In: Mathematics of Computation, ISSN 0025-5718, E-ISSN 1088-6842, Vol. 82, no 282, p. 919-952Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_31_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:31:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_31_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); Gaussian beams are asymptotically valid high frequency solutions to hyperbolic partial differential equations, concentrated on a single curve through the physical domain. They can also be extended to some dispersive wave equations, such as the Schrodinger equation. Superpositions of Gaussian beams provide a powerful tool to generate more general high frequency solutions that are not necessarily concentrated on a single curve. This work is concerned with the accuracy of Gaussian beam superpositions in terms of the wavelength epsilon. We present a systematic construction of Gaussian beam superpositions for all strictly hyperbolic and Schrodinger equations subject to highly oscillatory initial data of the form Ae(i Phi/) (epsilon). Through a careful estimate of an oscillatory integral operator, we prove that the k-th order Gaussian beam superposition converges to the original wave field at a rate proportional to epsilon(k/2) in the appropriate norm dictated by the well-posedness estimate. In particular, we prove that the Gaussian beam superposition converges at this rate for the acoustic wave equation in the standard, epsilon-scaled, energy norm and for the Schrodinger equation in the L-2 norm. The obtained results are valid for any number of spatial dimensions and are unaffected by the presence of caustics. We present a numerical study of convergence for the constant coefficient acoustic wave equation in R-2 to analyze the sharpness of the theoretical results.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:31:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 33. Malenova, Gabriela PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_32_j_idt585",{id:"formSmash:items:resultList:32:j_idt585",widgetVar:"widget_formSmash_items_resultList_32_j_idt585",onLabel:"Malenova, Gabriela ",offLabel:"Malenova, Gabriela ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_32_j_idt588",{id:"formSmash:items:resultList:32:j_idt588",widgetVar:"widget_formSmash_items_resultList_32_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:32:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Motamed, M.Runborg, OlofKTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.Tempone, R.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:32:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); A Sparse Stochastic Collocation Technique for High-Frequency Wave Propagation with Uncertainty2016In: SIAM/ASA Journal on Uncertainty Quantification, ISSN 1560-7526, E-ISSN 2166-2525, Vol. 4, no 1, p. 1084-1110Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_32_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:32:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_32_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); We consider the wave equation with highly oscillatory initial data, where there is uncertainty in the wave speed, initial phase, and/or initial amplitude. To estimate quantities of interest related to the solution and their statistics, we combine a high-frequency method based on Gaussian beams with sparse stochastic collocation. Although the wave solution, u(epsilon), is highly oscillatory in both physical and stochastic spaces, we provide theoretical arguments for simplified problems and numerical evidence that quantities of interest based on local averages of |u(epsilon)|(2) are smooth, with derivatives in the stochastic space uniformly bounded in epsilon, where epsilon denotes the short wavelength. This observable related regularity makes the sparse stochastic collocation approach more efficient than Monte Carlo methods. We present numerical tests that demonstrate this advantage.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:32:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 34. Malenova, Gabriela PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_33_j_idt585",{id:"formSmash:items:resultList:33:j_idt585",widgetVar:"widget_formSmash_items_resultList_33_j_idt585",onLabel:"Malenova, Gabriela ",offLabel:"Malenova, Gabriela ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_33_j_idt588",{id:"formSmash:items:resultList:33:j_idt588",widgetVar:"widget_formSmash_items_resultList_33_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA. KTH, Centres, SeRC - Swedish e-Science Research Centre.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:33:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Motamed, MohammadDepartment of Mathematics and Statistics, The University of New Mexico, Albuquerque, NM 87131, USA.Runborg, OlofKTH, Centres, SeRC - Swedish e-Science Research Centre. KTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:33:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Regularity of stochastic observables in Gaussian beam superposition of high-frequency waves2017In: Research in mathematical sciences, ISSN 2197-9847Article in journal (Other academic)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_33_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:33:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_33_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); We consider high-frequency waves satisfying the scalar wave equation with highly oscillatory initial data. The wave speed, and the phase and amplitude of the initial data are assumed to be uncertain, described by a finite number of random variables with known probability distributions. We define quantities of interest (QoIs), or observables, as local averages of the squared modulus of the wave solution. We aim to quantify the regularity of these QoIs in terms of the input random parameters, and the wave length, i.e., to estimate the size of their derivatives. The regularity is important for uncertainty quantification methods based on interpolation in the stochastic space. In particular, the size of the derivatives should be bounded independently of the wave length. In this paper, we are able to show that when these QoIs are approximated by Gaussian beam superpositions, they indeed have this property, despite the highly oscillatory character of the waves.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:33:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 35. Malenova, Gabriela PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_34_j_idt585",{id:"formSmash:items:resultList:34:j_idt585",widgetVar:"widget_formSmash_items_resultList_34_j_idt585",onLabel:"Malenova, Gabriela ",offLabel:"Malenova, Gabriela ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_34_j_idt588",{id:"formSmash:items:resultList:34:j_idt588",widgetVar:"widget_formSmash_items_resultList_34_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:34:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Runborg, OlofKTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:34:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Fast evaluation of high frequency observables by by Gaussian beams and the numerical steepest descent methodManuscript (preprint) (Other academic)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_34_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:34:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_34_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); We consider observables (quantities of interest) of the wave equation where the initial data has a wavelength which is very small compared to the distance traveled by the wave. Computation by direct methods becomes expensive owing to the short wavelength. Therefore, two fast asymptotic methods are proposed to address the difficulties: the Gaussian beam method to approximate the wave solution, and the numerical steepest descent method to compute the quantity of interest. In this work, we demonstrate the use of the latter method and present the corresponding error estimates. We show that the latter method has a lower computational cost compared to the previously considered methods and present the corresponding error estimates.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:34:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 36. Malenova, Gabriela PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_35_j_idt585",{id:"formSmash:items:resultList:35:j_idt585",widgetVar:"widget_formSmash_items_resultList_35_j_idt585",onLabel:"Malenova, Gabriela ",offLabel:"Malenova, Gabriela ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_35_j_idt588",{id:"formSmash:items:resultList:35:j_idt588",widgetVar:"widget_formSmash_items_resultList_35_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:35:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Runborg, OlofKTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA. KTH, Centres, SeRC - Swedish e-Science Research Centre.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:35:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Stochastic regularity of general quadratic observables of high frequency wavesManuscript (preprint) (Other academic)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_35_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:35:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_35_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); We consider solutions to the wave equation with uncertain initial data and medium, whose wavelength is short compared to thedistance traveled by the wave. We are interested in the statistics of the observables, i.e. functionals of the wave solution. Computation by direct methods gets very expensive or outright non-feasible as the wavelength decreases. To address the difficulties, we proposed a method consisting of the Gaussian beam method to treat the high frequencies and the sparse stochastic collocation method to remedy the curse of dimensionality in the stochastic space. For the latter method to converge, we need the observables to satisfy certain stochastic regularity conditions. The main contribution of this work is to show this regularity for a set of quadratic observables obtained by the Gaussian beam approximation of the wave solution.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:35:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 37. Marin, Oana PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_36_j_idt585",{id:"formSmash:items:resultList:36:j_idt585",widgetVar:"widget_formSmash_items_resultList_36_j_idt585",onLabel:"Marin, Oana ",offLabel:"Marin, Oana ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_36_j_idt588",{id:"formSmash:items:resultList:36:j_idt588",widgetVar:"widget_formSmash_items_resultList_36_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:36:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Runborg, OlofKTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA. KTH, Centres, SeRC - Swedish e-Science Research Centre.Tornberg, Anna-KarinKTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA. KTH, Centres, SeRC - Swedish e-Science Research Centre.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:36:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Corrected trapezoidal rules for a class of singular functions2014In: IMA Journal of Numerical Analysis, ISSN 0272-4979, E-ISSN 1464-3642, Vol. 34, no 4, p. 1509-1540Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_36_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:36:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_36_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); A set of accurate quadrature rules applicable to a class of integrable functions with isolated singularities is designed and analysed theoretically in one and two dimensions. These quadrature rules are based on the trapezoidal rule with corrected quadrature weights for points in the vicinity of the singularity. To compute the correction weights, small-size ill-conditioned systems have to be solved. The convergence of the correction weights is accelerated by the use of compactly supported functions that annihilate boundary errors. Convergence proofs with error estimates for the resulting quadrature rules are given in both one and two dimensions. The tabulated weights are specific for the singularities under consideration, but the methodology extends to a large class of functions with integrable isolated singularities. Furthermore, in one dimension we have obtained a closed form expression based on which the modified weights can be computed directly.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:36:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 38. Motamed, M. et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_37_j_idt588",{id:"formSmash:items:resultList:37:j_idt588",widgetVar:"widget_formSmash_items_resultList_37_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:37:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Runborg, OlofKTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA. KTH, Centres, SeRC - Swedish e-Science Research Centre.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:37:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); A wavefront-based Gaussian beam method for computing high frequency wave propagation problems2015In: Computers and Mathematics with Applications, ISSN 0898-1221, E-ISSN 1873-7668, Vol. 69, no 9, p. 949-963Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_37_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:37:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_37_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); We present a novel wavefront method based on Gaussian beams for computing high frequency wave propagation problems. Unlike standard geometrical optics, Gaussian beams compute the correct solution of the wave field also at caustics. The method tracks a front of two canonical beams with two particular initial values for width and curvature. In a fast post-processing step, from the canonical solutions we recreate any other Gaussian beam with arbitrary initial data on the initial front. This provides a simple mechanism to include a variety of optimization processes, including error minimization and beam width minimization, for a posteriori selection of optimal beam initial parameters. The performance of the method is illustrated with two numerical examples.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:37:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 39. Motamed, Mohammad PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_38_j_idt585",{id:"formSmash:items:resultList:38:j_idt585",widgetVar:"widget_formSmash_items_resultList_38_j_idt585",onLabel:"Motamed, Mohammad ",offLabel:"Motamed, Mohammad ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_38_j_idt588",{id:"formSmash:items:resultList:38:j_idt588",widgetVar:"widget_formSmash_items_resultList_38_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Computer Science and Communication (CSC), Numerical Analysis and Computer Science, NADA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:38:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Runborg, OlofKTH, School of Computer Science and Communication (CSC), Numerical Analysis and Computer Science, NADA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:38:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); A fast method for the creeping ray contribution to scattering problems2006In: Mathematical Modeling of Wave Phenomena / [ed] Nilsson, B; Fishman, L, MELVILLE, NY: AMER INST PHYSICS , 2006, Vol. 834, p. 56-64Conference paper (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_38_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:38:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_38_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); Creeping rays can give an important contribution to the solution of medium to high frequency scattering problems. They are generated at the shadow lines of the illuminated scatterer by grazing incident rays and propagate along geodesics on the scatterer surface, continuously shedding diffracted rays in their tangential direction. In this paper we show how the ray propagation problem can be formulated as a partial differential equation (PDE) in a three-dimensional phase space. To solve the PDE we use a fast marching method. The PDE solution contains information about all possible creeping rays. This information includes the phase and amplitude of the field, which are extracted by a fast postprocessing. Computationally the cost of solving the PDE is less than tracing all rays individually by solving a system of ordinary differential equations(ODE). We consider an application to monostatic radar cross section problems where creeping rays from all illumination angles must be computed.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:38:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 40. Motamed, Mohammad PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_39_j_idt585",{id:"formSmash:items:resultList:39:j_idt585",widgetVar:"widget_formSmash_items_resultList_39_j_idt585",onLabel:"Motamed, Mohammad ",offLabel:"Motamed, Mohammad ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_39_j_idt588",{id:"formSmash:items:resultList:39:j_idt588",widgetVar:"widget_formSmash_items_resultList_39_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Computer Science and Communication (CSC), Numerical Analysis and Computer Science, NADA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:39:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Runborg, OlofKTH, School of Computer Science and Communication (CSC), Numerical Analysis and Computer Science, NADA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:39:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); A fast phase space method for computing creeping rays2006In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 219, no 1, p. 276-295Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_39_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:39:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_39_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); Creeping rays can give an important contribution to the solution of medium to high frequency scattering problems. They are generated at the shadow lines of the illuminated scatterer by grazing incident rays and propagate along geodesics on the scatterer surface, continuously shedding diffracted rays in their tangential direction. In this paper, we show how the ray propagation problem can be formulated as a partial differential equation (PDE) in a three-dimensional phase space. To solve the PDE we use a fast marching method. The PDE solution contains information about all possible creeping rays. This information includes the phase and amplitude of the field, which are extracted by a fast post-processing. Computationally, the cost of solving the PDE is less than tracing all rays individually by solving a system of ordinary differential equations. We consider an application to mono-static radar cross section problems where creeping rays from all illumination angles must be computed. The numerical results of the fast phase space method and a comparison with the results of ray tracing are presented.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:39:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 41. Motamed, Mohammad PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_40_j_idt585",{id:"formSmash:items:resultList:40:j_idt585",widgetVar:"widget_formSmash_items_resultList_40_j_idt585",onLabel:"Motamed, Mohammad ",offLabel:"Motamed, Mohammad ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_40_j_idt588",{id:"formSmash:items:resultList:40:j_idt588",widgetVar:"widget_formSmash_items_resultList_40_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Computer Science and Communication (CSC), Numerical Analysis and Computer Science, NADA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:40:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Runborg, OlofKTH, School of Computer Science and Communication (CSC), Numerical Analysis and Computer Science, NADA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:40:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); A multiple-patch phase space method for computing trajectories on manifolds with applications to wave propagation problems2007In: Communications in Mathematical Sciences, ISSN 1539-6746, E-ISSN 1945-0796, Vol. 5, no 3, p. 617-648Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_40_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:40:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_40_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); We present a multiple-patch phase space method for computing trajectories on two-dimensional manifolds possibly embedded in a higher-dimensional space. The dynamics of trajectories are given by systems of ordinary differential equations (ODEs). We split the manifold into multiple patches where each patch has a well-defined regular parameterization. The ODEs are formulated as escape equations, which are hyperbolic partial differential equations (PDEs) in a three- dimensional phase space. The escape equations are solved in each patch, individually. The solutions of individual patches are then connected using suitable inter-patch boundary conditions. Properties for particular families of trajectories are obtained through a fast post-processing. We apply the method to two different problems : the creeping ray contribution to mono-static radar cross section computations and the multivalued travel-time of seismic waves in multi-layered media. We present numerical examples to illustrate the accuracy and efficiency of the method.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:40:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 42. Motamed, Mohammad PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_41_j_idt585",{id:"formSmash:items:resultList:41:j_idt585",widgetVar:"widget_formSmash_items_resultList_41_j_idt585",onLabel:"Motamed, Mohammad ",offLabel:"Motamed, Mohammad ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_41_j_idt588",{id:"formSmash:items:resultList:41:j_idt588",widgetVar:"widget_formSmash_items_resultList_41_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:41:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Runborg, OlofKTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:41:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); A Wave Front-based Gaussian Beam Method for Computing High Frequency WavesManuscript (Other academic)43. Motamed, Mohammad PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_42_j_idt585",{id:"formSmash:items:resultList:42:j_idt585",widgetVar:"widget_formSmash_items_resultList_42_j_idt585",onLabel:"Motamed, Mohammad ",offLabel:"Motamed, Mohammad ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_42_j_idt588",{id:"formSmash:items:resultList:42:j_idt588",widgetVar:"widget_formSmash_items_resultList_42_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); Department of Mathematics, Simon Fraser University.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:42:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Runborg, OlofKTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:42:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Taylor expansion and discretization errors in Gaussian beam superposition2010In: Wave motion, ISSN 0165-2125, E-ISSN 1878-433X, Vol. 47, no 7, p. 421-439Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_42_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:42:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_42_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); The Gaussian beam superposition method is an asymptotic method for computing high frequency wave fields in smoothly varying inhomogeneous media. In this paper we study the accuracy of the Gaussian beam superposition method and derive error estimates related to the discretization of the superposition integral and the Taylor expansion of the phase and amplitude off the center of the beam. We show that in the case of odd order beams, the error is smaller than a simple analysis would indicate because of error cancellation effects between the beams. Since the cancellation happens only when odd order beams are used, there is no remarkable gain in using even order beams. Moreover, applying the error estimate to the problem with constant speed of propagation, we show that in this case the local beam width is not a good indicator of accuracy, and there is no direct relation between the error and the beam width. We present numerical examples to verify the error estimates.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:42:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 44. Motamed, Mohammad PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_43_j_idt585",{id:"formSmash:items:resultList:43:j_idt585",widgetVar:"widget_formSmash_items_resultList_43_j_idt585",onLabel:"Motamed, Mohammad ",offLabel:"Motamed, Mohammad ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_43_j_idt588",{id:"formSmash:items:resultList:43:j_idt588",widgetVar:"widget_formSmash_items_resultList_43_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:43:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Runborg, OlofKTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:43:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Taylor Expansion Errors in Gaussian Beam SummationManuscript (Other academic)45. Möller, Joakim PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_44_j_idt585",{id:"formSmash:items:resultList:44:j_idt585",widgetVar:"widget_formSmash_items_resultList_44_j_idt585",onLabel:"Möller, Joakim ",offLabel:"Möller, Joakim ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_44_j_idt588",{id:"formSmash:items:resultList:44:j_idt588",widgetVar:"widget_formSmash_items_resultList_44_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:44:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Runborg, OlofKTH, School of Computer Science and Communication (CSC), Numerical Analysis and Computer Science, NADA.Kevrekidis, P.G.Lust, K.Kevrekidis, I.G.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:44:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Equation-free, effective computation for discrete systems: a time stepper based approach2005In: International Journal of Bifurcation and Chaos in Applied Sciences and Engineering, ISSN 0218-1274, Vol. 15, no 3, p. 975-996Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_44_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:44:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_44_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); We propose a computer-assisted approach to studying the effective continuum behavior of spatially discrete evolution equations. The advantage of the approach is that the "coarse model" (the continuum, effective equation) need not be explicitly constructed. The method only uses a time-integration code for the discrete problem and judicious choices of initial data and integration times; our bifurcation computations are based on the so-called Recursive Projection Method (RPM) with arc-length continuation [Shroff & Keller, 1993]. The technique is used to monitor features of the genuinely discrete problem such as the pinning of coherent structures and its results are compared to quasi-continuum approaches such as the ones based on Pade approximations.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:44:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 46. Persson, P. O. et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_45_j_idt588",{id:"formSmash:items:resultList:45:j_idt588",widgetVar:"widget_formSmash_items_resultList_45_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:45:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Runborg, OlofPrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:45:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Simulation of a waveguide filter using wavelet-based numerical homogenization2001In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 166, no 2, p. 361-382Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_45_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:45:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_45_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); We apply wavelet-based numerical homogenization to the simulation of an optical waveguide filter. We use the method to derive approximate one-dimensional models and subgrid models of the filter. Numerical examples of the technique are presented, and the computational gains are investigated.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:45:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 47. Popovic, Jelena PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_46_j_idt585",{id:"formSmash:items:resultList:46:j_idt585",widgetVar:"widget_formSmash_items_resultList_46_j_idt585",onLabel:"Popovic, Jelena ",offLabel:"Popovic, Jelena ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_46_j_idt588",{id:"formSmash:items:resultList:46:j_idt588",widgetVar:"widget_formSmash_items_resultList_46_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:46:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Runborg, OlofKTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:46:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); A Fast Method for Solving the Helmholtz Equation Based on Wave-Splitting2009In: WAVES 2009 / [ed] Barucq, H.; Bonnet-Bendhia, A.-S.; Cohen, G.; Diaz, J.; Ezziani, A.; Joly, P., 2009, p. 220-221Conference paper (Refereed)48. Popovic, Jelena et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_47_j_idt588",{id:"formSmash:items:resultList:47:j_idt588",widgetVar:"widget_formSmash_items_resultList_47_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:47:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Runborg, OlofKTH, Centres, SeRC - Swedish e-Science Research Centre. KTH, School of Engineering Sciences (SCI), Mathematics (Dept.).PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:47:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Adaptive fast interface tracking methods2017In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 337, p. 42-61Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_47_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:47:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_47_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); In this paper, we present a fast time adaptive numerical method for interface tracking. The method uses an explicit multiresolution description of the interface, which is represented by wavelet vectors that correspond to the details of the interface on different scale levels. The complexity of standard numerical methods for interface tracking, where the interface is described by N marker points, is 0 (N/Delta t),when a time step At is used. The methods that we propose in this paper have 0 (TOL-1/P log N + N log N) computational cost, at least for uniformly smooth problems, where TOL is some given tolerance and p is the order of the time stepping method that is used for time advection of the interface. The adaptive method is robust in the sense that it can handle problems with both smooth and piecewise smooth interfaces (e.g. interfaces with corners) while keeping a low computational cost. We show numerical examples that verify these properties.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:47:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 49. Popovic, Jelena PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_48_j_idt585",{id:"formSmash:items:resultList:48:j_idt585",widgetVar:"widget_formSmash_items_resultList_48_j_idt585",onLabel:"Popovic, Jelena ",offLabel:"Popovic, Jelena ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_48_j_idt588",{id:"formSmash:items:resultList:48:j_idt588",widgetVar:"widget_formSmash_items_resultList_48_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:48:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Runborg, OlofKTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:48:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Adaptive Fast Interface Tracking Methods: Part I: Time AdaptivityManuscript (preprint) (Other academic)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_48_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:48:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_48_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); In this paper, we present a fast adaptive numerical method for interface tracking that uses an explicit multiresolution description of the interface. The interface is represented by wavelet vectors that correspond to the details of the interface on different scale levels.The complexity of standard numerical methods for interface tracking, where the interface is described by marker points, is O(N/dt), where N is the number of points on the interface and dt is the time step. The methods that we propose in this paper have O(tol^(-1/p)log N) computational cost, where tol is some given tolerance and p is the order of the numerical method for ordinary differential equations that is used for time advection of the interface. The adaptivity makes methods robust in the sense that they can handle problems with both smooth and non-smooth interfaces (i.e. interfaces with corners) while keeping low computational cost.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:48:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 50. Popovic, Jelena PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_49_j_idt585",{id:"formSmash:items:resultList:49:j_idt585",widgetVar:"widget_formSmash_items_resultList_49_j_idt585",onLabel:"Popovic, Jelena ",offLabel:"Popovic, Jelena ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_49_j_idt588",{id:"formSmash:items:resultList:49:j_idt588",widgetVar:"widget_formSmash_items_resultList_49_j_idt588",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:49:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Runborg, OlofKTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:49:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Adaptive Fast Interface Tracking Methods: Part II: Spatial AdaptivityManuscript (preprint) (Other academic)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_49_j_idt623_0_j_idt624",{id:"formSmash:items:resultList:49:j_idt623:0:j_idt624",widgetVar:"widget_formSmash_items_resultList_49_j_idt623_0_j_idt624",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); In this paper, we present a fast space-time adaptive numerical method for interface propagation in a time varying velocity field based on a multiresolution description of the interface. The interface is represented by wavelet vectors that correspond to the details of the interface on different scale levels.The method is an extension of the method proposed in "J. Popovic and O. Runborg, Adaptive fast interface tracking methods: Part I, preprint (2012)", which is only time adaptive and it is thus not suitable for problems with expanding interfaces. The method that we propose in this paper, remedies that disadvantage of the time adaptive method in an efficient way.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:49:j_idt623:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500});

CiteExportLink to result list
http://kth.diva-portal.org/smash/resultList.jsf?query=&language=en&searchType=SIMPLE&noOfRows=50&sortOrder=author_sort_asc&sortOrder2=title_sort_asc&onlyFullText=false&sf=all&aq=%5B%5B%7B%22personId%22%3A%22u1na5bn9%22%7D%5D%5D&aqe=%5B%5D&aq2=%5B%5B%5D%5D&af=%5B%5D $(function(){PrimeFaces.cw("InputTextarea","widget_formSmash_lower_j_idt903_recordPermLink",{id:"formSmash:lower:j_idt903:recordPermLink",widgetVar:"widget_formSmash_lower_j_idt903_recordPermLink",autoResize:true});}); $(function(){PrimeFaces.cw("OverlayPanel","widget_formSmash_lower_j_idt903_j_idt905",{id:"formSmash:lower:j_idt903:j_idt905",widgetVar:"widget_formSmash_lower_j_idt903_j_idt905",target:"formSmash:lower:j_idt903:permLink",showEffect:"blind",hideEffect:"fade",my:"right top",at:"right bottom",showCloseIcon:true});});

Permanent link

Cite

Citation styleapa ieee modern-language-association-8th-edition vancouver Other style $(function(){PrimeFaces.cw("SelectOneMenu","widget_formSmash_lower_j_idt921",{id:"formSmash:lower:j_idt921",widgetVar:"widget_formSmash_lower_j_idt921",behaviors:{change:function(ext) {PrimeFaces.ab({s:"formSmash:lower:j_idt921",e:"change",f:"formSmash",p:"formSmash:lower:j_idt921",u:"formSmash:lower:otherStyle"},ext);}}});});

- apa
- ieee
- modern-language-association-8th-edition
- vancouver
- Other style

Languagede-DE en-GB en-US fi-FI nn-NO nn-NB sv-SE Other locale $(function(){PrimeFaces.cw("SelectOneMenu","widget_formSmash_lower_j_idt932",{id:"formSmash:lower:j_idt932",widgetVar:"widget_formSmash_lower_j_idt932",behaviors:{change:function(ext) {PrimeFaces.ab({s:"formSmash:lower:j_idt932",e:"change",f:"formSmash",p:"formSmash:lower:j_idt932",u:"formSmash:lower:otherLanguage"},ext);}}});});

- de-DE
- en-GB
- en-US
- fi-FI
- nn-NO
- nn-NB
- sv-SE
- Other locale

Output formathtml text asciidoc rtf $(function(){PrimeFaces.cw("SelectOneMenu","widget_formSmash_lower_j_idt942",{id:"formSmash:lower:j_idt942",widgetVar:"widget_formSmash_lower_j_idt942"});});

- html
- text
- asciidoc
- rtf