arxiv: v1 [math.ap] 11 Jan 2014
|
|
- Ashlyn Wiggins
- 5 years ago
- Views:
Transcription
1 THE UNIFIED TRANSFORM FOR THE MODIFIED HELMHOLTZ EQUATION IN THE EXTERIOR OF A SQUARE A. S. FOKAS AND J. LENELLS arxiv:4.252v [math.ap] Jan 24 Abstract. The Unified Transform provides a novel method for analyzing boundary value problems for linear and for integrable nonlinear PDEs. The numerical implementation of this method to linear elliptic PDEs formulated in the interior of a polygon has been investigated by several authors (see the article by Iserles, Smitheman, and one of the authors in this book). Here, we show that the Unified Transform also yields a novel numerical technique for computing the solution of linear elliptic PDEs in the exterior of a polygon. One of the advantages of this new technique is that it actually yields directly the scattering amplitude. Details are presented for the modified Helmholtz equation in the exterior of a square. AMS Subject Classification (2): 35J5, 35J25, 65N35. Keywords: linear elliptic PDE, boundary value problem, exterior domain, scattering theory.. Introduction The so-called Unified Transform, introduced by the first author in the late nineties [3, 6], has been implemented to a variety of boundary value problems for linear elliptic PDEs. In particular, it has led to the emergence of a new numerical technique for the Laplace, the Helmholtz, and the modified Helmholtz equations formulated in the interior of a polygon [, 2, 4, 7, 8,, 2, 4 8]. This technique is illustrated for the particular case of a square in this book []. A novel numerical implementation of the Unified Transform to the Helmholtz and modified Helmholtz equations formulated in the exterior of a polygon is presented in [9]. Here, we illustrate this implementation in the particular case of a square. 2. The Unified Transform The starting point of the Unified Transform is rewriting the given PDE as the condition that a one-parameter family of differential forms is closed. For the modified Helmholtz equation 2 u x u y 2 4β2 u =, β >, (x, y) D R 2, (2.) the relevant differential form is given by [( ) W (z, z, ) = e iβ(z z u ) z + iβu where dz ( u z + βu ) ] d z, (2.2) i z = x + iy, z = x iy, z D. (2.3) Indeed, recall that the differential of the one-form W, W (z, z) = A(z, z)dz + B(z, z)d z,
2 2 A. S. FOKAS AND J. LENELLS is given by dw = ( B z A ) dz d z. (2.4) z Using the above expression for the differential form W defined by (2.2) we find ( dw = 2e iβ(z z ) 2 ) u z z β2 u dz d z, z D. (2.5) However, employing the change of variables defined by (2.3), we find that the modified Helmholtz equation (2.) can be written in the form thus (2.5) implies 2 u z z β2 u =, z D, (2.6) dw =, z D. (2.7) 3. The Global Relation for Exterior Domains Suppose that the domain D is exterior to some bounded domain D B. Then equation (2.7) and the complex form of Green s theorem imply D B W (z, z, ) = I(), C, (3.) where I() is defined by I() = lim W (z, z, ), C. (3.2) r z =r Assume that u satisfies the Sommerfeld radiation condition, i.e. ( ) lim r r r + 2β u =. (3.3) It is well known (see [3]) that as r, u admits the asymptotic expansion ( u πiβr ei(2iβr π 4 ) f (ϕ) + f (ϕ) + f ) 2(ϕ) r r 2 +, r. (3.4) Employing the above representation in (3.2), it is shown in [9] that I() is given by I() = 4if (i ln( i)), C. (3.5) It is shown in [9] that f (ϕ) admits an analytic continuation to complex values of ϕ and hence the right-hand side of (3.5) is well defined. We recall that if the modified Helmholtz equation (2.) is defined in the interior domain D B, then W satisfies the equation D B W (z, z, ) =, C. (3.6) It turns out that if D B is the interior of a polygon, then equation (3.6) provides an efficient way for computing the generalized Dirichlet to Neumann map (see []). Indeed, for a well posed problem, one specifies a relationship between the Dirichlet and the Neumann boundary values. Thus, if the relevant polygon consists of n sides, then equation (3.6) contains n unknowns. For example, in the case of the Dirichlet problem, equation (3.6) contains the n unknown Neumann boundary values. However, equation (3.6) is valid for all complex values of. Thus, by expanding the unknown boundary
3 THE UNIFIED TRANSFORM FOR THE MODIFIED HELMHOLTZ 3 values in an appropriate basis, and by evaluating equation (3.6) at a sufficiently large number of collocation points, we can obtain the unknown coefficients appearing in the above expansions. Fornberg and colleagues have shown that by choosing Legendre polynomials as base functions, and by overdetermining the relevant linear system (i.e. by choosing the number of collocation points much larger than the number of unknown coefficients), the above technique yields a linear system with a small condition number. Furthermore, it is shown in [8], that by choosing the collocation points to be on certain curves, it is possible to obtain a condition number which is independent of β. Equation (3.) differs from equation (3.6) only by the presence of the explicit term defined in (3.5). Thus, we can still apply the techniques developed for analyzing equation (3.6), but we must supplement equation (3.) with an additional equation in order to compensate for the existence of the additional unknown function f. 4. A Supplement to the Global Relation Suppose that D is the exterior of the n-gon with vertices at {z j } n. It is shown in [9] and [9] that in this case the following relation is valid: n j= e iβ(z z l j ) û j () =, z D B, (4.) where D B denotes the interior of the polygon, the rays l j oriented from the origin to infinity are defined by and the functions {û j ()} n l j = { C arg = arg(z j+ z j )}, j =,..., n, (4.2) û j () = zj+ z j are defined by ( ) ] e iβ(z z [iu ) d z N + iβ dt + dz u dt, (4.3) dt where u N denotes the Neumann boundary value. Let us parametrize the side [z j, z j+ ] by where z(t) = m j + th j, t, (4.4) m j = z j + z j+ 2, h j = z j+ z j. (4.5) 2 Then, the expression û j () defined in (4.3) can be written in the form û j () = û D j () + û N j (), j =,..., n, (4.6) where ( ) û D j () = iβ j + h j e h iβ(m j m j ) e iβt(h j h j ) udt, (4.7) û N j () = ie iβ(m j m j ) e iβt(h j h j ) u N dt. (4.8)
4 4 A. S. FOKAS AND J. LENELLS 5. The Dirichlet Problem For the Dirichlet problem the functions {û D j ()}n are known (see (4.7)). In order to compute û N j (), we approximate the Neumann boundary value on the side j as follows: M u N j (t) c (j) m P m (t), j =,..., n, (5.) m= where P m is the Legendre polynomial of degree m, and the unknown constants c (j) m are to be determined. Equation (4.8) yields û N j () = ie iβ(m j m j ) M m= c (j) m e iβt(h j h j ) P m (t)dt. (5.2) Thus, using the expression [8] m e iλx A mp e iλ + B mp e iλ P m (x)dx = Λ p+, m, Λ C \ {}, (5.3) where equation (5.2) becomes p= (m + p)!( )m+p A mp = 2 p p!(m p)!i p+, B (m + p)! mp = 2 p, (5.4) p!(m p)!ip+ û N j () = ie iβ(m j m j ) Thus, equation (4.) yields where u D (z) + i n M c (j) m j= m= m p= M m= u D (z) = c (j) m m A mp e iλ + B mp e iλ Λ p+. (5.5) Λ=β(hj h j ) p= A mp Q + pj (z) + B mpq pj (z) β p+ =, z D B, (5.6) n j= and the functions Q ± pj (z) are defined by Q ± pj (z) = z iβ(z e ˆlj e iβ(z z m j+ m j ±(h j h j )) ˆlj (h j h j )p+ ) û D j (), (5.7), (5.8) with ˆl j denoting the contour l j with small indentations inserted appropriately so that ˆl j passes around the points = ± exp( i arg h j ) (these points are removable singularities, thus we deform l j to ˆl j before splitting the relevant integrals). Suppose we take the limit as z approaches the j th side. Then, using for this side the parametrization z = m j + t h j, t, (5.9)
5 THE UNIFIED TRANSFORM FOR THE MODIFIED HELMHOLTZ 5 we find Q σ pj (z ) = Q σ pj(z ) = ˆlj e iβ(t +σ)[h j h j ] (h j h j )p+ ˆlj e iβ[(m j m j )+(t h j +σh j ), σ = ±, (5.) (h j h j )p+ m j m j t h j +σ h j ], j j, σ = ±. By employing the above limiting procedure, we can obtain n additional equations. Let D B be the square with corners at 6. The Exterior of a Square z = + i, z 2 = + i, z 3 = i, z 4 = i. Then, the four midpoints of the sides are m = i, m 2 =, m 3 = i, m 4 =, with the corresponding tangential vectors Also, h =, h 2 = i, h 3 =, h 4 = i. l = (, ), l 2 = (, i ), l 3 = (, ), l 4 = (, i ). Thus, if j = so that z = i t, t, belongs to the top side, then Q σ p(z ) = Q σ p2(z ) = Q σ p3(z ) = Q σ p4(z ) = ˆl e iβ(t +σ)[ ] ( + )p+ ˆl2 e iβ[i( σ)(+ )+( t )( )] ( i i )p+, σ = ±,, σ = ±, ˆl3 e iβ[2i(+ )+(σ t )( )] ( )p+, σ = ±, ˆl4 e iβ[i(+σ)(+ ) (+t )( )] (i + i )p+, σ = ±. (5.) We choose the indented contours ˆl j, j =,..., 4, displayed in Figure. The exponential e iβ(t +σ)[ ] has decay for Im < if σ = and for Im > if σ =. Using that e iβ(t +σ)[ ] Res = ( + )p+ = d p ( ) p+ p! p e iβ(t +σ)[ ] =, we find Q σ p(z ) = e 3iπ 4 + 2πi p! e iβ(t +σ)[ ] ( + )p+ d p ( p = ) p+ e iβ(t +σ)[ ], σ =,
6 6 A. S. FOKAS AND J. LENELLS ˆl2 i ˆl3 ˆl i ˆl4 Figure. The indented contours ˆl j, j =,..., 4, utilized for the exterior of the square. and Similarly, we find Q σ p(z ) = e 3iπ 4 e iβ(t +σ)[ ] ( + )p+, σ =. Q σ p2(z ) = Q σ p3(z ) = iπ e 4 iπ e 4 e iβ[i( σ)(+ )+( t )( )] ( i i )p+, σ = ±, e iβ[2i(+ )+(σ t )( )] ( )p+, σ =, e iπ 4 e iβ[2i(+ )+(σ t )( )] Q σ p3(z ) = ( )p+ 2πi d p ( ) p+ p! p e iβ[2i(+ )+(σ t )( )] = +, σ =, Q σ p4(z ) = e iπ 4 e iβ[ (+t )( )+i(+σ)(+ )] (i + i )p+, σ = ±. The above integrals have exponential decay and are suitable for numerical evaluation. Similar considerations lead to analogous expressions for Q ± pj (z ) for z belonging to the left, bottom, and right sides of the square. Substituting the above expressions into relation (5.6) evaluated at z = z where z belongs to one of the sides of the square, we find relations for the c (j) m s which can be used to supplement the exterior global relation. Taking into consideration that the global relation now involves one additional unknown (the scattering amplitude f ), one would expect that it is sufficient to supplement the global relation with the equation obtained by varying z on any one of the sides of the square. However, it appears that in order to obtain accurate results it is necessary to consider all four sides.
7 THE UNIFIED TRANSFORM FOR THE MODIFIED HELMHOLTZ 7 Re f (ϕ) ϕ Im f (ϕ) ϕ Figure 2. The real and imaginary parts of the scattering amplitude f (ϕ) for Example. Numerical values (solid) compared with true values (dashed). Re u y (i t) t Im u y (i t) t Figure 3. The real and imaginary parts of the Neumann value on the top side for Example. It was noted earlier that the term u D defined by (5.7) is known. However, the question of computing this term numerically in the most efficient way remains open. In the examples below we compute the known expression (5.7) by approximating the given Dirichlet data in terms of Legendre polynomials. 6.. Examples. Let D be the domain exterior to the square with vertices at ± ± i. The function u(r, ϕ) = H () n (2iβr)(Ae inϕ + Be inϕ ), (6.) satisfies the modified Helmholtz equation (2.) in D for any n and any complex constants A and B Example. For the first example, we let u(r, ϕ) = H () (2iβr). This solution satisfies ( u(r, ϕ) = πiβr ei(2iβr π 4 ) ) 6βr β 2 r 2 + O(r 3 ), r. Hence equation (3.4) is satisfied with f (ϕ). We use the values of u on the four sides of the square as the given Dirichlet data. Letting β =, we apply the above technique with M = 8 and with 54 collocation points. The numerical scheme yields the scattering amplitude f (ϕ) displayed in Figure 2. The numerically obtained and
8 8 A. S. FOKAS AND J. LENELLS Re f (ϕ) Im f (ϕ) ϕ ϕ Figure 4. The real and imaginary parts of the scattering amplitude f (ϕ) for Example 2. Re u y (i t) Im u y (i t) t t Figure 5. The real and imaginary parts of the Neumann value on the top side for Example 2. true Neumann values on the top side are shown in Figure 3; the results for the other three sides are similar. The numerics gives the correct profiles with a maximal error of approximately Example 2. For the second example, we let u(r, ϕ) = H () (2iβr)eiϕ. This solution satisfies ) u(r, ϕ) = πiβr ei(2iβr π 4 ( ie ) iϕ + 3eiϕ 6iβr + 5ieiϕ 52β 2 r 2 + O(r 3 ), r. Hence equation (3.4) is satisfied with f (ϕ) = ie iϕ. Letting β =, the numerical scheme with M = 8 and with 54 collocation points yields the scattering amplitude f (ϕ) displayed in Figure 4. The numerically obtained and true Neumann values on the top side are shown in Figure 5. The numerics gives the correct profiles with a maximal error of approximately Conclusions We have illustrated a novel numerical technique for linear elliptic PDEs formulated in the exterior of a polygon. This technique, which is based on the Unified Transform, uses two main ingredients: (a) the global relation; in the case of the modified Helmholtz equation, the global relation is given by (3.), where I() is defined by equation (3.5). This equation can be analyzed numerically in the same way that
9 THE UNIFIED TRANSFORM FOR THE MODIFIED HELMHOLTZ 9 the global relation is analyzed in the case of linear elliptic PDEs formulated in the interior of a polygon. However, now the global relation, in addition to the unknown boundary values, also involves the unknown function f (the scattering amplitude). Hence, it must be supplemented with an additional equation. (b) Additional equations can be obtained by taking the limit of a certain equation which is valid for z in the interior of the given polygon, as z approaches the boundary; in the case of the modified Helmholtz equation, this equation is given by (5.6). Actually, by taking the limit as z approaches each side of an n-gon, one obtains n additional equations. These equations provide the formulation in the complex Fourier plane (complex extension of ) of the classical boundary integral method. By exploiting the analyticity structure in the complex Fourier plane, it is possible to perform contour deformations and to obtain integrals which decay exponentially for large. This yields a most efficient numerical evaluation of these integrals. Detailed comparisons will be presented elsewhere. Acknowledgement The authors acknowledge support from the EPSRC, UK. References [] A. C. L. Ashton, On the rigorous foundations of the Fokas method for linear elliptic partial differential equations, Proc. R. Soc. A 468 (22), [2] A. C. L. Ashton, The spectral Dirichlet-Neumann map for the Laplace equation in a convex polygon, SIAM J. Math. Anal. 45 (23), [3] A. S. Fokas, A unified transform method for solving linear and certain nonlinear PDEs, Proc. Roy. Soc. Lond. A 453 (997), [4] A. S. Fokas, Two-dimensional linear partial differential equations in a convex polygon, R. Soc. Lond. Proc. Ser. A Math. Phys. Eng. Sci. 457 (2), [5] A. S. Fokas, Integrable nonlinear evolution equations on the half-line, Comm. Math. Phys. 23 (22), 39. [6] A. S. Fokas, A unified approach to boundary value problems, CBMS-NSF regional conference series in applied mathematics, SIAM (28). [7] A. S. Fokas, N. Flyer, S. A. Smitheman, and E. A. Spence, A semi-analytical numerical method for solving evolution and elliptic partial differential equations, J. Comput. Appl. Math. 227 (29), [8] A. S. Fokas, A. C. L. Ashton, and S. A. Smitheman, A numerical technique for the interior of a polygonal domain, preprint. [9] A. S. Fokas and J. Lenells, The Unified Transform for the Helmholtz and modified Helmholtz equations in the exterior of a polygon, in preparation. [] B. Fornberg and N. Flyer, A numerical implementation of Fokas boundary integral approach: Laplace s equation on a polygonal domain, Proc. R. Soc. Lond. Ser. A Math. Phys. Eng. Sci. 467 (2), [] A. S. Fokas, A. Iserles, and S. A. Smitheman, The unified transform in polygonal domains via the explicit Fourier transform of Legendre polynomials, in this volume. [2] S. R. Fulton, A. S. Fokas, and C. A. Xenophontos, An analytical method for linear elliptic PDEs and its numerical implementation, J. Comput. Appl. Math. 67 (24), [3] S. N. Karp, A convergent farfield expansion for two-dimensional radiation functions, Comm. Pure Appl. Math. 4 (96), [4] Y. G. Saridakis, A. G. Sifalakis, and E. P. Papadopoulou, Efficient numerical solution of the generalized Dirichlet-Neumann map for linear elliptic PDEs in regular polygon domains, J. Comput. Appl. Math. 236 (22), [5] A. G. Sifalakis, A. S. Fokas, S. R. Fulton, and Y. G. Saridakis, The generalized Dirichlet- Neumann map for linear elliptic PDEs and its numerical implementation, J. Comput. Appl. Math. 29 (28), 9 34.
10 A. S. FOKAS AND J. LENELLS [6] A. G. Sifalakis, S. R. Fulton, E. P. Papadopoulou, and Y. G. Saridakis, Direct and iterative solution of the generalized Dirichlet-Neumann map for elliptic PDEs on square domains, J. Comput. Appl. Math. 227 (29), [7] A. G. Sifalakis, E. P. Papadopoulou, and Y. G. Saridakis, Numerical study of iterative methods for the solution of the Dirichlet-Neumann map for linear elliptic PDEs on regular polygon domains, Int. J. Appl. Math. Comput. Sci 4 (27), [8] S. A. Smitheman, E. A. Spence, and A. S. Fokas, A spectral collocation method for the Laplace and modified Helmholtz equations in a convex polygon, IMA J. Numer. Anal. 3 (2), [9] E. A. Spence, Boundary Value Problems for Linear Elliptic PDEs, Ph.D. thesis, University of Cambridge, Cambridge, UK, 2. A.S.F.: Department of Applied Mathematics and Theoretical Physics, University of Cambridge, Cambridge CB3 WA, United Kingdom, and Research Center of Mathematics, Academy of Athens, 527, Greece. address: T.Fokas@damtp.cam.ac.uk J.L.: Department of Mathematics and Center for Astrophysics, Space Physics & Engineering Research, Baylor University, One Bear Place #97328, Waco, TX 76798, USA. address: Jonatan Lenells@baylor.edu
arxiv: v1 [math.na] 21 Nov 2012
The Fourier Transforms of the Chebyshev and Legendre Polynomials arxiv:111.4943v1 [math.na] 1 Nov 1 A S Fokas School of Engineering and Applied Sciences, Harvard University, Cambridge, MA 138, USA Permanent
More informationBoundary value problems and medical imaging
Journal of Physics: Conference Series OPEN ACCESS Boundary value problems and medical imaging To cite this article: Athanasios S Fokas and George A Kastis 214 J. Phys.: Conf. Ser. 49 1217 View the article
More informationThe elliptic sinh-gordon equation in the half plane
Available online at www.tjnsa.com J. Nonlinear Sci. Appl. 8 25), 63 73 Research Article The elliptic sinh-gordon equation in the half plane Guenbo Hwang Department of Mathematics, Daegu University, Gyeongsan
More informationDirect and Iterative Solution of the Generalized Dirichlet-Neumann Map for Elliptic PDEs on Square Domains
Direct and Iterative Solution of the Generalized Dirichlet-Neumann Map for Elliptic PDEs on Square Domains A.G.Sifalakis 1, S.R.Fulton, E.P. Papadopoulou 1 and Y.G.Saridakis 1, 1 Applied Mathematics and
More informationIBVPs for linear and integrable nonlinear evolution PDEs
IBVPs for linear and integrable nonlinear evolution PDEs Dionyssis Mantzavinos Department of Applied Mathematics and Theoretical Physics, University of Cambridge. Edinburgh, May 31, 212. Dionyssis Mantzavinos
More informationWell-posedness of boundary-value problems for the linear Benjamin-Bona-Mahony equation
Well-posedness of boundary-value problems for the linear Benjamin-Bona-Mahony equation Vishal Vasan and Bernard Deconinck Department of Applied Mathematics University of Washington Seattle WA 98195 July
More informationBoundary value problems for the elliptic sine Gordon equation in a semi strip
Boundary value problems for the elliptic sine Gordon equation in a semi strip Article Accepted Version Fokas, A. S., Lenells, J. and Pelloni, B. 3 Boundary value problems for the elliptic sine Gordon equation
More informationA Generalization of both the Method of Images and of the Classical Integral Transforms
A Generalization of both the Method of Images and of the Classical Integral Transforms Athanassios S. Fokas 1 and Daniel ben-avraham 2 1 Institute for Nonlinear Studies, Clarkson University, Potsdam, NY
More informationc 2007 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol. 9, No. 3, pp. 07 6 c 007 Society for Industrial and Applied Mathematics NUMERICAL QUADRATURES FOR SINGULAR AND HYPERSINGULAR INTEGRALS IN BOUNDARY ELEMENT METHODS MICHAEL CARLEY
More informationA new method for the solution of scattering problems
A new method for the solution of scattering problems Thorsten Hohage, Frank Schmidt and Lin Zschiedrich Konrad-Zuse-Zentrum Berlin, hohage@zibde * after February 22: University of Göttingen Abstract We
More informationWell posed two point initial boundary value problems with arbitrary boundary conditions
Well posed two point initial boundary value problems with arbitrary boundary conditions Article Accepted Version Smith, D. A. (2012) Well posed two point initial boundary value problems with arbitrary
More informationProceedings of the 5th International Conference on Inverse Problems in Engineering: Theory and Practice, Cambridge, UK, 11-15th July 2005
Proceedings of the 5th International Conference on Inverse Problems in Engineering: Theory and Practice, Cambridge, UK, 11-15th July 2005 SOME INVERSE SCATTERING PROBLEMS FOR TWO-DIMENSIONAL SCHRÖDINGER
More informationThe Pole Condition: A Padé Approximation of the Dirichlet to Neumann Operator
The Pole Condition: A Padé Approximation of the Dirichlet to Neumann Operator Martin J. Gander and Achim Schädle Mathematics Section, University of Geneva, CH-, Geneva, Switzerland, Martin.gander@unige.ch
More informationRational Chebyshev pseudospectral method for long-short wave equations
Journal of Physics: Conference Series PAPER OPE ACCESS Rational Chebyshev pseudospectral method for long-short wave equations To cite this article: Zeting Liu and Shujuan Lv 07 J. Phys.: Conf. Ser. 84
More informationOn spherical-wave scattering by a spherical scatterer and related near-field inverse problems
IMA Journal of Applied Mathematics (2001) 66, 539 549 On spherical-wave scattering by a spherical scatterer and related near-field inverse problems C. ATHANASIADIS Department of Mathematics, University
More informationThe Watson transformation revisited
The Watson transformation revisited E. A. Spence June 17, 214 Abstract The Watson transformation is essentially the idea that a series can be converted into an integral in the complex plane via an inverse
More informationThe Factorization Method for Inverse Scattering Problems Part I
The Factorization Method for Inverse Scattering Problems Part I Andreas Kirsch Madrid 2011 Department of Mathematics KIT University of the State of Baden-Württemberg and National Large-scale Research Center
More informationNEW ANALYTICAL SOLUTION FOR SOLVING STEADY-STATE HEAT CONDUCTION PROBLEMS WITH SINGULARITIES
THERMAL SCIENCE: Year 3, Vol. 7, No. 3, pp. 665-67 665 NEW ANALYTICAL SOLUTION FOR SOLVING STEADY-STATE HEAT CONDUCTION PROBLEMS WITH SINGULARITIES by Najib LARAQI a* and Eric MONIER-VINARD b a Paris West
More informationProblems for MATH-6300 Complex Analysis
Problems for MATH-63 Complex Analysis Gregor Kovačič December, 7 This list will change as the semester goes on. Please make sure you always have the newest version of it.. Prove the following Theorem For
More informationCONSTRUCTION OF THE HALF-LINE POTENTIAL FROM THE JOST FUNCTION
CONSTRUCTION OF THE HALF-LINE POTENTIAL FROM THE JOST FUNCTION Tuncay Aktosun Department of Mathematics and Statistics Mississippi State University Mississippi State, MS 39762 Abstract: For the one-dimensional
More informationLinear independence of exponentials on the real line
Linear independence of exponentials on the real line Alexandre Eremenko August 25, 2013 The following question was asked on Math Overflow. Let (λ n ) be a sequence of complex numbers tending to infinity.
More informationarxiv: v1 [math.na] 15 Nov 2008
arxiv:0811.2481v1 [math.na] 15 Nov 2008 A Phase-Fitted Runge-Kutta-Nyström method for the Numerical Solution of Initial Value Problems with Oscillating Solutions Abstract D. F. Papadopoulos a, Z. A. Anastassi
More informationConsidering our result for the sum and product of analytic functions, this means that for (a 0, a 1,..., a N ) C N+1, the polynomial.
Lecture 3 Usual complex functions MATH-GA 245.00 Complex Variables Polynomials. Construction f : z z is analytic on all of C since its real and imaginary parts satisfy the Cauchy-Riemann relations and
More informationTHE ZERO-DISTRIBUTION AND THE ASYMPTOTIC BEHAVIOR OF A FOURIER INTEGRAL. Haseo Ki and Young One Kim
THE ZERO-DISTRIBUTION AND THE ASYMPTOTIC BEHAVIOR OF A FOURIER INTEGRAL Haseo Ki Young One Kim Abstract. The zero-distribution of the Fourier integral Q(u)eP (u)+izu du, where P is a polynomial with leading
More informationOn the decay of elements of inverse triangular Toeplitz matrix
On the decay of elements of inverse triangular Toeplitz matrix Neville Ford, D. V. Savostyanov, N. L. Zamarashkin August 03, 203 arxiv:308.0724v [math.na] 3 Aug 203 Abstract We consider half-infinite triangular
More informationAdditive Schwarz method for scattering problems using the PML method at interfaces
Additive Schwarz method for scattering problems using the PML method at interfaces Achim Schädle 1 and Lin Zschiedrich 2 1 Zuse Institute, Takustr. 7, 14195 Berlin, Germany schaedle@zib.de 2 Zuse Institute,
More informationThe Residue Theorem. Integration Methods over Closed Curves for Functions with Singularities
The Residue Theorem Integration Methods over losed urves for Functions with Singularities We have shown that if f(z) is analytic inside and on a closed curve, then f(z)dz = 0. We have also seen examples
More informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
More informationDispersive and soliton perturbations of finite-genus solutions of the KdV equation: computational results
Dispersive and soliton perturbations of finite-genus solutions of the KdV equation: computational results Thomas Trogdon and Bernard Deconinck Department of Applied Mathematics University of Washington
More informationGeometry-Fitted Fourier-Mellin Transform Pairs
Geometry-Fitted Fourier-Mellin Transform Pairs Darren Crowdy Abstract The construction of novel Fourier/Mellin-type transform pairs that are tailor-made for given planar regions within the special class
More informationMath 411, Complex Analysis Definitions, Formulas and Theorems Winter y = sinα
Math 411, Complex Analysis Definitions, Formulas and Theorems Winter 014 Trigonometric Functions of Special Angles α, degrees α, radians sin α cos α tan α 0 0 0 1 0 30 π 6 45 π 4 1 3 1 3 1 y = sinα π 90,
More informationThe Linear KdV Equation with an Interface
The Linear KdV Equation with an Interface Bernard Deconinc, Natalie E. Sheils, and David A. Smith 3. Department of Applied Mathematics, University of Washington, Seattle, WA 9895-395. School of Mathematics,
More informationComplex geometrical optics solutions for Lipschitz conductivities
Rev. Mat. Iberoamericana 19 (2003), 57 72 Complex geometrical optics solutions for Lipschitz conductivities Lassi Päivärinta, Alexander Panchenko and Gunther Uhlmann Abstract We prove the existence of
More informationPSI Lectures on Complex Analysis
PSI Lectures on Complex Analysis Tibra Ali August 14, 14 Lecture 4 1 Evaluating integrals using the residue theorem ecall the residue theorem. If f (z) has singularities at z 1, z,..., z k which are enclosed
More informationAnalytical solution for determination the control parameter in the inverse parabolic equation using HAM
Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 12, Issue 2 (December 2017, pp. 1072 1087 Applications and Applied Mathematics: An International Journal (AAM Analytical solution
More informationTechnische Universität Graz
Technische Universität Graz A note on the stable coupling of finite and boundary elements O. Steinbach Berichte aus dem Institut für Numerische Mathematik Bericht 2009/4 Technische Universität Graz A
More informationAsymptotic of Enumerative Invariants in CP 2
Peking Mathematical Journal https://doi.org/.7/s4543-8-4-4 ORIGINAL ARTICLE Asymptotic of Enumerative Invariants in CP Gang Tian Dongyi Wei Received: 8 March 8 / Revised: 3 July 8 / Accepted: 5 July 8
More informationThe Method of Fundamental Solutions applied to the numerical calculation of eigenfrequencies and eigenmodes for 3D simply connected domains
ECCOMAS Thematic Conference on Meshless Methods 2005 C34.1 The Method of Fundamental Solutions applied to the numerical calculation of eigenfrequencies and eigenmodes for 3D simply connected domains Carlos
More informationFINITE ELEMENT APPROXIMATION OF ELLIPTIC DIRICHLET OPTIMAL CONTROL PROBLEMS
Numerical Functional Analysis and Optimization, 28(7 8):957 973, 2007 Copyright Taylor & Francis Group, LLC ISSN: 0163-0563 print/1532-2467 online DOI: 10.1080/01630560701493305 FINITE ELEMENT APPROXIMATION
More informationExplicit kernel-split panel-based Nyström schemes for planar or axisymmetric Helmholtz problems
z Explicit kernel-split panel-based Nyström schemes for planar or axisymmetric Helmholtz problems Johan Helsing Lund University Talk at Integral equation methods: fast algorithms and applications, Banff,
More informationarxiv: v1 [math.cv] 18 Aug 2015
arxiv:508.04376v [math.cv] 8 Aug 205 Saddle-point integration of C bump functions Steven G. Johnson, MIT Applied Mathematics Created November 23, 2006; updated August 9, 205 Abstract This technical note
More informationEvaluation of integrals
Evaluation of certain contour integrals: Type I Type I: Integrals of the form 2π F (cos θ, sin θ) dθ If we take z = e iθ, then cos θ = 1 (z + 1 ), sin θ = 1 (z 1 dz ) and dθ = 2 z 2i z iz. Substituting
More informationThe Pearcey integral in the highly oscillatory region
The Pearcey integral in the highly oscillatory region José L. López 1 and Pedro Pagola arxiv:1511.0533v1 [math.ca] 17 Nov 015 1 Dpto. de Ingeniería Matemática e Informática and INAMAT, Universidad Pública
More informationA FAST SOLVER FOR ELLIPTIC EQUATIONS WITH HARMONIC COEFFICIENT APPROXIMATIONS
Proceedings of ALGORITMY 2005 pp. 222 229 A FAST SOLVER FOR ELLIPTIC EQUATIONS WITH HARMONIC COEFFICIENT APPROXIMATIONS ELENA BRAVERMAN, MOSHE ISRAELI, AND ALEXANDER SHERMAN Abstract. Based on a fast subtractional
More informationInteraction of lumps with a line soliton for the DSII equation
Physica D 152 153 (2001) 189 198 Interaction of lumps with a line soliton for the DSII equation A.S. Fokas a,b, D.E. Pelinovsky c,, C. Sulem c a Department of Mathematics, Imperial College, London SW7
More informationQ You mentioned that in complex analysis we study analytic functions, or, in another name, holomorphic functions. Pray tell me, what are they?
COMPLEX ANALYSIS PART 2: ANALYTIC FUNCTIONS Q You mentioned that in complex analysis we study analytic functions, or, in another name, holomorphic functions. Pray tell me, what are they? A There are many
More informationSurface Tension Effect on a Two Dimensional. Channel Flow against an Inclined Wall
Applied Mathematical Sciences, Vol. 1, 007, no. 7, 313-36 Surface Tension Effect on a Two Dimensional Channel Flow against an Inclined Wall A. Merzougui *, H. Mekias ** and F. Guechi ** * Département de
More informationDmitry Pelinovsky Department of Mathematics, McMaster University, Canada
Spectrum of the linearized NLS problem Dmitry Pelinovsky Department of Mathematics, McMaster University, Canada Collaboration: Marina Chugunova (McMaster, Canada) Scipio Cuccagna (Modena and Reggio Emilia,
More informationExponentially Accurate Semiclassical Tunneling Wave Functions in One Dimension
Exponentially Accurate Semiclassical Tunneling Wave Functions in One Dimension Vasile Gradinaru Seminar for Applied Mathematics ETH Zürich CH 8092 Zürich, Switzerland, George A. Hagedorn Department of
More informationDunkl operators and Clifford algebras II
Translation operator for the Clifford Research Group Department of Mathematical Analysis Ghent University Hong Kong, March, 2011 Translation operator for the Hermite polynomials Translation operator for
More informationlim when the limit on the right exists, the improper integral is said to converge to that limit.
hapter 7 Applications of esidues - evaluation of definite and improper integrals occurring in real analysis and applied math - finding inverse Laplace transform by the methods of summing residues. 6. Evaluation
More informationMath 185 Fall 2015, Sample Final Exam Solutions
Math 185 Fall 2015, Sample Final Exam Solutions Nikhil Srivastava December 12, 2015 1. True or false: (a) If f is analytic in the annulus A = {z : 1 < z < 2} then there exist functions g and h such that
More information1 Discussion on multi-valued functions
Week 3 notes, Math 7651 1 Discussion on multi-valued functions Log function : Note that if z is written in its polar representation: z = r e iθ, where r = z and θ = arg z, then log z log r + i θ + 2inπ
More informationThe Calculus of Residues
hapter 7 The alculus of Residues If fz) has a pole of order m at z = z, it can be written as Eq. 6.7), or fz) = φz) = a z z ) + a z z ) +... + a m z z ) m, 7.) where φz) is analytic in the neighborhood
More informationT u
WANG LI-LIAN Assistant Professor Division of Mathematical Sciences, SPMS Nanyang Technological University Singapore, 637616 T 65-6513-7669 u 65-6316-6984 k lilian@ntu.edu.sg http://www.ntu.edu.sg/home/lilian?
More informationHere are brief notes about topics covered in class on complex numbers, focusing on what is not covered in the textbook.
Phys374, Spring 2008, Prof. Ted Jacobson Department of Physics, University of Maryland Complex numbers version 5/21/08 Here are brief notes about topics covered in class on complex numbers, focusing on
More informationc 2012 Society for Industrial and Applied Mathematics
SIAM J. NUMER. ANAL. Vol. 50, No. 4, pp. 849 860 c 0 Society for Industrial and Applied Mathematics TWO RESULTS CONCERNING THE STABILITY OF STAGGERED MULTISTEP METHODS MICHELLE GHRIST AND BENGT FORNBERG
More information18.04 Practice problems exam 1, Spring 2018 Solutions
8.4 Practice problems exam, Spring 8 Solutions Problem. omplex arithmetic (a) Find the real and imaginary part of z + z. (b) Solve z 4 i =. (c) Find all possible values of i. (d) Express cos(4x) in terms
More informationEvaluation of Chebyshev pseudospectral methods for third order differential equations
Numerical Algorithms 16 (1997) 255 281 255 Evaluation of Chebyshev pseudospectral methods for third order differential equations Rosemary Renaut a and Yi Su b a Department of Mathematics, Arizona State
More informationPublications: Charles Fulton. Papers On Sturm-Liouville Theory
Publications: Charles Fulton Papers On Sturm-Liouville Theory References [1] C. Fulton, Parametrizations of Titchmarsh s m(λ)- Functions in the limit circle case, Trans. Amer. Math. Soc. 229, (1977), 51-63.
More informationResearch Statement. James Bremer Department of Mathematics, University of California, Davis
Research Statement James Bremer Department of Mathematics, University of California, Davis Email: bremer@math.ucdavis.edu Webpage: https.math.ucdavis.edu/ bremer I work in the field of numerical analysis,
More informationFree-surface potential flow of an ideal fluid due to a singular sink
Journal of Physics: Conference Series PAPER OPEN ACCESS Free-surface potential flow of an ideal fluid due to a singular sink To cite this article: A A Mestnikova and V N Starovoitov 216 J. Phys.: Conf.
More informationNotes by Zvi Rosen. Thanks to Alyssa Palfreyman for supplements.
Lecture: Hélène Barcelo Analytic Combinatorics ECCO 202, Bogotá Notes by Zvi Rosen. Thanks to Alyssa Palfreyman for supplements.. Tuesday, June 2, 202 Combinatorics is the study of finite structures that
More informationNonlinear Integral Equation Formulation of Orthogonal Polynomials
Nonlinear Integral Equation Formulation of Orthogonal Polynomials Eli Ben-Naim Theory Division, Los Alamos National Laboratory with: Carl Bender (Washington University, St. Louis) C.M. Bender and E. Ben-Naim,
More informationQualitative Theory of Differential Equations and Dynamics of Quadratic Rational Functions
Nonl. Analysis and Differential Equations, Vol. 2, 2014, no. 1, 45-59 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/nade.2014.3819 Qualitative Theory of Differential Equations and Dynamics of
More informationTHE MAGNUS METHOD FOR SOLVING OSCILLATORY LIE-TYPE ORDINARY DIFFERENTIAL EQUATIONS
THE MAGNUS METHOD FOR SOLVING OSCILLATORY LIE-TYPE ORDINARY DIFFERENTIAL EQUATIONS MARIANNA KHANAMIRYAN Department of Applied Mathematics and Theoretical Physics, University of Cambridge, Wilberforce Road,
More informationChapter II. Complex Variables
hapter II. omplex Variables Dates: October 2, 4, 7, 2002. These three lectures will cover the following sections of the text book by Keener. 6.1. omplex valued functions and branch cuts; 6.2.1. Differentiation
More informationNumerical Simulations on Two Nonlinear Biharmonic Evolution Equations
Numerical Simulations on Two Nonlinear Biharmonic Evolution Equations Ming-Jun Lai, Chun Liu, and Paul Wenston Abstract We numerically simulate the following two nonlinear evolution equations with a fourth
More informationAn Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations
An Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations Moshe Israeli Computer Science Department, Technion-Israel Institute of Technology, Technion city, Haifa 32000, ISRAEL Alexander
More informationWeakly non-linear completely resonant hamiltonian PDEs and the problem of weak turbulence
Sergei Kuksin Weakly non-linear completely resonant hamiltonian PDEs and the problem of weak turbulence (Toronto, 10 January, 2014 ) 1 1 Introduction: weak turbulence (WT) (One of) origines: Rudolf Peierls,
More informationGABRIELA JARAMILLO EMPLOYMENT & EDUCATION. Assistant Professor, University of Houston, TX
GABRIELA JARAMILLO ADDRESS Department of Mathematics University of Houston 3551 Cullen Blvd. Room 641 Philip Guthrie Hoffman Hall Houston, TX 77204 CONTACT INFORMATION email: gabriela@math.uh.edu website:
More informationFourier transforms, Generalised functions and Greens functions
Fourier transforms, Generalised functions and Greens functions T. Johnson 2015-01-23 Electromagnetic Processes In Dispersive Media, Lecture 2 - T. Johnson 1 Motivation A big part of this course concerns
More informationDiagonalization of the Coupled-Mode System.
Diagonalization of the Coupled-Mode System. Marina Chugunova joint work with Dmitry Pelinovsky Department of Mathematics, McMaster University, Canada Collaborators: Mason A. Porter, California Institute
More informationTOPICS IN NONLINEAR ANALYSIS AND APPLICATIONS. Dipartimento di Matematica e Applicazioni Università di Milano Bicocca March 15-16, 2017
TOPICS IN NONLINEAR ANALYSIS AND APPLICATIONS Dipartimento di Matematica e Applicazioni Università di Milano Bicocca March 15-16, 2017 Abstracts of the talks Spectral stability under removal of small capacity
More informationarxiv:math/ v1 [math.ap] 28 Oct 2005
arxiv:math/050643v [math.ap] 28 Oct 2005 A remark on asymptotic completeness for the critical nonlinear Klein-Gordon equation Hans Lindblad and Avy Soffer University of California at San Diego and Rutgers
More informationUltraspherical moments on a set of disjoint intervals
Ultraspherical moments on a set of disjoint intervals arxiv:90.049v [math.ca] 4 Jan 09 Hashem Alsabi Université des Sciences et Technologies, Lille, France hashem.alsabi@gmail.com James Griffin Department
More informationSplitting methods with boundary corrections
Splitting methods with boundary corrections Alexander Ostermann University of Innsbruck, Austria Joint work with Lukas Einkemmer Verona, April/May 2017 Strang s paper, SIAM J. Numer. Anal., 1968 S (5)
More informationInitial Boundary Value Problems for Scalar and Vector Burgers Equations
Initial Boundary Value Problems for Scalar and Vector Burgers Equations By K. T. Joseph and P. L. Sachdev In this article we stu Burgers equation and vector Burgers equation with initial and boundary conditions.
More informationMidterm Examination #2
Anthony Austin MATH 47 April, 9 AMDG Midterm Examination #. The integrand has poles of order whenever 4, which occurs when is equal to i, i,, or. Since a R, a >, the only one of these poles that lies inside
More informationInitial-Boundary Value Problem for Two-Component Gerdjikov Ivanov Equation with 3 3 Lax Pair on Half-Line
Commun. Theor. Phys. 68 27 425 438 Vol. 68, No. 4, October, 27 Initial-Boundary Value Problem for Two-Component Gerdjiov Ivanov Equation with 3 3 Lax Pair on Half-Line Qiao-Zhen Zhu 朱巧珍,, En-Gui Fan 范恩贵,,
More informationNon-radial solutions to a bi-harmonic equation with negative exponent
Non-radial solutions to a bi-harmonic equation with negative exponent Ali Hyder Department of Mathematics, University of British Columbia, Vancouver BC V6TZ2, Canada ali.hyder@math.ubc.ca Juncheng Wei
More informationPerformance Evaluation of Generalized Polynomial Chaos
Performance Evaluation of Generalized Polynomial Chaos Dongbin Xiu, Didier Lucor, C.-H. Su, and George Em Karniadakis 1 Division of Applied Mathematics, Brown University, Providence, RI 02912, USA, gk@dam.brown.edu
More informationNumerical Analysis and Methods for PDE I
Numerical Analysis and Methods for PDE I A. J. Meir Department of Mathematics and Statistics Auburn University US-Africa Advanced Study Institute on Analysis, Dynamical Systems, and Mathematical Modeling
More informationON FACTORIZATIONS OF ENTIRE FUNCTIONS OF BOUNDED TYPE
Annales Academiæ Scientiarum Fennicæ Mathematica Volumen 29, 2004, 345 356 ON FACTORIZATIONS OF ENTIRE FUNCTIONS OF BOUNDED TYPE Liang-Wen Liao and Chung-Chun Yang Nanjing University, Department of Mathematics
More informationOn Bank-Laine functions
Computational Methods and Function Theory Volume 00 0000), No. 0, 000 000 XXYYYZZ On Bank-Laine functions Alastair Fletcher Keywords. Bank-Laine functions, zeros. 2000 MSC. 30D35, 34M05. Abstract. In this
More informationSpectral Functions for Regular Sturm-Liouville Problems
Spectral Functions for Regular Sturm-Liouville Problems Guglielmo Fucci Department of Mathematics East Carolina University May 15, 13 Regular One-dimensional Sturm-Liouville Problems Let I = [, 1 R, and
More informationA modification of the factorization method for scatterers with different physical properties
A modification of the factorization method for scatterers with different physical properties Takashi FURUYA arxiv:1802.05404v2 [math.ap] 25 Oct 2018 Abstract We study an inverse acoustic scattering problem
More informationCoercivity of high-frequency scattering problems
Coercivity of high-frequency scattering problems Valery Smyshlyaev Department of Mathematics, University College London Joint work with: Euan Spence (Bath), Ilia Kamotski (UCL); Comm Pure Appl Math 2015.
More informationChapter 30 MSMYP1 Further Complex Variable Theory
Chapter 30 MSMYP Further Complex Variable Theory (30.) Multifunctions A multifunction is a function that may take many values at the same point. Clearly such functions are problematic for an analytic study,
More information1 Complex Numbers. 1.1 Sums and Products
1 Complex Numbers 1.1 Sums Products Definition: The complex plane, denoted C is the set of all ordered pairs (x, y) with x, y R, where Re z = x is called the real part Imz = y is called the imaginary part.
More informationSpectra of the Gurtin-Pipkin type equations
Spectra of the Gurtin-Pipkin type equations A. Eremenko and S. Ivanov September 28, 2 Abstract We study the spectra of certain integro-differential equations arising in applications. Under some conditions
More informationDiscontinuous Galerkin methods for fractional diffusion problems
Discontinuous Galerkin methods for fractional diffusion problems Bill McLean Kassem Mustapha School of Maths and Stats, University of NSW KFUPM, Dhahran Leipzig, 7 October, 2010 Outline Sub-diffusion Equation
More informationZeta Functions and Regularized Determinants for Elliptic Operators. Elmar Schrohe Institut für Analysis
Zeta Functions and Regularized Determinants for Elliptic Operators Elmar Schrohe Institut für Analysis PDE: The Sound of Drums How Things Started If you heard, in a dark room, two drums playing, a large
More informationHarmonic Mappings and Shear Construction. References. Introduction - Definitions
Harmonic Mappings and Shear Construction KAUS 212 Stockholm, Sweden Tri Quach Department of Mathematics and Systems Analysis Aalto University School of Science Joint work with S. Ponnusamy and A. Rasila
More informationNew series expansions of the Gauss hypergeometric function
New series expansions of the Gauss hypergeometric function José L. López and Nico. M. Temme 2 Departamento de Matemática e Informática, Universidad Pública de Navarra, 36-Pamplona, Spain. e-mail: jl.lopez@unavarra.es.
More informationA Direct Method for reconstructing inclusions from Electrostatic Data
A Direct Method for reconstructing inclusions from Electrostatic Data Isaac Harris Texas A&M University, Department of Mathematics College Station, Texas 77843-3368 iharris@math.tamu.edu Joint work with:
More informationDetermination of thin elastic inclusions from boundary measurements.
Determination of thin elastic inclusions from boundary measurements. Elena Beretta in collaboration with E. Francini, S. Vessella, E. Kim and J. Lee September 7, 2010 E. Beretta (Università di Roma La
More informationFeedback control of transient energy growth in subcritical plane Poiseuille flow
Feedback control of transient energy growth in subcritical plane Poiseuille flow Fulvio Martinelli, Maurizio Quadrio, John McKernan and James F. Whidborne Abstract Subcritical flows may experience large
More informationComplex Inversion Formula for Stieltjes and Widder Transforms with Applications
Int. J. Contemp. Math. Sciences, Vol. 3, 8, no. 16, 761-77 Complex Inversion Formula for Stieltjes and Widder Transforms with Applications A. Aghili and A. Ansari Department of Mathematics, Faculty of
More informationTrefftz type method for 2D problems of electromagnetic scattering from inhomogeneous bodies.
Trefftz type method for 2D problems of electromagnetic scattering from inhomogeneous bodies. S. Yu. Reutsiy Magnetohydrodynamic Laboratory, P. O. Box 136, Mosovsi av.,199, 61037, Kharov, Uraine. e mail:
More information