Spectral computation of highly oscillatory integral equations in laser theory

2019 ◽  
Vol 395 ◽  
pp. 351-381
Author(s):  
Jing Gao ◽  
Marissa Condon ◽  
Arieh Iserles
Author(s):  
M. Condon ◽  
A. Iserles ◽  
S. P. Nørsett

The concern of this paper is in expanding and computing initial-value problems of the form y ′= f ( y )+ h ω ( t ), where the function h ω oscillates rapidly for ω ≫1. Asymptotic expansions for such equations are well understood in the case of modulated Fourier oscillators and they can be used as an organizing principle for very accurate and affordable numerical solvers. However, there is no similar theory for more general oscillators, and there are sound reasons to believe that approximations of this kind are unsuitable in that setting. We follow in this paper an alternative route, demonstrating that, for a much more general family of oscillators, e.g. linear combinations of functions of the form e i ωg k ( t ) , it is possible to expand y ( t ) in a different manner. Each r th term in the expansion is for some ς >0 and it can be represented as an r -dimensional highly oscillatory integral. Because computation of multivariate highly oscillatory integrals is fairly well understood, this provides a powerful method for an effective discretization of a numerical solution for a large family of highly oscillatory ordinary differential equations.


Author(s):  
Jae-bok Lee ◽  
Jun Zou ◽  
Benliang Li ◽  
Munno Ju

Purpose – The per-unit-length earth return mutual impedance of the overhead conductors plays an important role for analyzing electromagnetic transients or couplings of multi-conductor systems. It is impossible to have a closed-form expression to evaluate this kind of impedance. The purpose of this paper is to propose an efficient numerical approach to evaluate the earth return mutual impedance of the overhead conductors above horizontally multi-layered soils. Design/methodology/approach – The expression of the earth return mutual impedance, which contains a complex highly oscillatory semi-infinite integral, is divided into two parts intentionally, i.e. the definite and the tail integral, respectively. The definite integral is calculated using the proposed moment functions after fitting the integrand into the piecewise cubic spline functions, and the tail integral is replaced by exponential integrals with newly developed asymptotic integrands. Findings – The numerical examples show the proposed approach has a satisfactory accuracy for different parameter combinations. Compared to the direct quadrature approach, the computational time of the proposed approach is very competitive, especially, for the large horizontal distance and the low height of the conductors. Originality/value – The advantage of the proposed approach is that the calculation of the highly oscillatory integral is completely avoided due to the fact that the moment function can be evaluated analytically. The contribution of the tail integral is well included by means of the exponential integral, though in an asymptotic way. The proposed approach is completely general, and can be applied to calculate the earth return mutual impedance of overhead conductors above a soil structure with an arbitrary number of horizontal layers.


Symmetry ◽  
2019 ◽  
Vol 11 (6) ◽  
pp. 728 ◽  
Author(s):  
SAIRA ◽  
Shuhuang Xiang

In this paper, a fast and accurate numerical Clenshaw-Curtis quadrature is proposed for the approximation of highly oscillatory integrals with Cauchy and logarithmic singularities, ⨍ − 1 1 f ( x ) log ( x − α ) e i k x x − t d x , t ∉ ( − 1 , 1 ) , α ∈ [ − 1 , 1 ] for a smooth function f ( x ) . This method consists of evaluation of the modified moments by stable recurrence relation and Cauchy kernel is solved by steepest descent method that transforms the oscillatory integral into the sum of line integrals. Later theoretical analysis and high accuracy of the method is illustrated by some examples.


Sign in / Sign up

Export Citation Format

Share Document