When is a Stokes line not a Stokes line?
|
|
- Tracy Harmon
- 6 years ago
- Views:
Transcription
1 T H E U N I V E R S I T Y O H F E D I N B U R G When is a Stokes line not a Stokes line? C.J.Howls University of Southampton, UK A.B. Olde Daalhuis University of Edinburgh, UK Work supported by EPSRC, CWI, RIMS
2 Collaboration Chris Howls, Phil Langman (Southampton) Adri Olde Daalhuis (Edinburgh) Jon Chapman, David Mortimer (Oxford) John King, Richard Tew, Giles Body (Nottingham) EPSRC-funded research.
3 So far... We ve seen what exponential asymptotics is for. Extends range of validity of expansions. Increases numerical accuracy. Provides a systematic way to calculate Stokes constants.
4 This talk... Recall some basic ideas: Stokes lines and caustics. Introduce you to an apparent asymptotic paradox. Explain why our understanding of asymptotics is missing a key tool of analytic continuation.
5 Stokes phenomenon and caustics
6 Stokes phenomenon and caustics z(λ)
7 Stokes phenomenon and caustics z(λ) 2 Main Bow Supernumerary Fringes
8 !! Stokes phenomenon and caustics Ai(z)!!!! z
9 !! Stokes phenomenon and caustics Ai(z)!!!! z Convergent Series: Airy 1838
10 !! Stokes phenomenon and caustics Ai(z)!!!! z Divergent Asymptotic Series: Stokes 1864 Convergent Series: Airy 1838
11 !! Stokes phenomenon and caustics Ai(z)!!!! z Ai(z) sin(2 z 3/2 /3 + π/4) πz 1/4 Ai(z) e 2z3/2 /3 2 πz 1/4
12 ! Stokes phenomenon and caustics Divergence of asymptotics Ai(z)!!!! z Ai(z) sin(2 z 3/2 /3 + π/4) πz 1/4 Ai(z) e 2z3/2 /3 2 πz 1/4
13 ! Stokes phenomenon and caustics Divergence of asymptotics Ai(z)!!!! z Ai(z) sin(2 z 3/2 /3 + π/4) πz 1/4 Ai(z) e 2z3/2 /3 2 πz 1/4 Coalescence of two exponents as z 0
14 ! Stokes phenomenon and caustics Divergence of asymptotics Ai(z)!!!! z Full expansion Ai(z) sin(2 z 3/2 /3 + π/4) πz 1/4 Ai(z) e 2z3/2 /3 2 πz 1/4 a r (z) r=0 Coalescence of two exponents as z 0
15 ! Stokes phenomenon and caustics a r = Divergence of every term in asymptotics Ai(z) Γ(r) (Difference in exponents) Γ(r) ( 4 3 z3/2), r!!!! z Full expansion Ai(z) sin(2 z 3/2 /3 + π/4) πz 1/4 Ai(z) e 2z3/2 /3 2 πz 1/4 a r (z) r=0 Coalescence of two exponents as z 0
16 Stokes phenomenon and caustics z C 0 Second exponential switches on in complex z- plane, where first maximally dominates it in size: arg(z)= 2π/3.
17 Stokes phenomenon and caustics z C 0 e 2z3/2 /3 2 πz 1/4 Second exponential switches on in complex z- plane, where first maximally dominates it in size: arg(z)= 2π/3.
18 Stokes phenomenon and caustics z C e 2z3/2 /3 2 πz 1/4 0 Second exponential switches on in complex z- plane, where first maximally dominates it in size: arg(z)= 2π/3.
19 Stokes phenomenon and caustics z C e 2z3/2 /3 e2z3/2/3 2 + i πz1/4 2 πz 1/4 0 Second exponential switches on in complex z- plane, where first maximally dominates it in size: arg(z)= 2π/3.
20 Stokes phenomenon and caustics z C e 2z3/2 /3 e2z3/2/3 2 + i πz1/4 2 πz 1/4 0 Stokes constant Second exponential switches on in complex z- plane, where first maximally dominates it in size: arg(z)= 2π/3.
21 Stokes phenomenon and caustics z C e 2z3/2 /3 e2z3/2/3 2 + i πz1/4 2 πz 1/4 0 Stokes lines Second exponential switches on in complex z- plane, where first maximally dominates it in size: arg(z)= 2π/3.
22 Stokes phenomenon and caustics z C anti-stokes lines 0 Anti-Stokes lines, where real parts of exponentials are identical. Vital in nonlinear problems or for numerical integration.
23 Stokes phenomenon and caustics z C turning point/ caustic 0 Stokes-lines/anti-Stokes lines sprout from turning point/caustic and extend to infinity, or to another turning point or singularity, e.g...
24 Gravitational waves in Reissner-Nordström black hole Turning points c anti-stokes lines a b t 1 t 2 t 3 t 4 f e' r- r+ t 6 t 5 d e pole c' Complex r plane Stokes lines
25 Pearcey function ( ) = exp i t 4 + yt 2 + xt Ψ 2 x,y + ( ( )) A highly oscillatory problem...
26 Pearcey function
27 Pearcey function Stokes line Stokes line Number and type of contributing saddles caustics/turning points
28 Summary Stokes lines, anti-stokes lines, turning points/ caustics exist and play key role in asymptotic continuation of approximations. Stokes line, anti-stokes lines sprout from turning point or singularities Stokes line, anti-stokes lines run between turning points or singularities or run to infinity.
29 Summary End of talk? Until recently, pretty much yes... But... The activity of the Stokes lines themselves can be dynamically switched on and off!
30 New phenomenon Activity of Stokes line plays a significant role in wide classes of asymptotics problems involving a (non-asymptotic) parameter: Integrals, ODEs and PDE time-evolution problems (linear and nonlinear) Effect on highly oscillatory problems. We show how an exponential approach is vital to the explanation of activity of Stokes lines. Explained by introducing a new concept: Higher Order Stokes Phenomenon
31 Activity of Stokes lines Sudden change in activity of Stokes line noticed by Berk, Nevins and Roberts (1982) in third order WKB approach to: i d3 y dy + 3i dx3 dx + xy = 0 In complex x plane they noticed crossing of Stokes lines. From monodromy of solutions they noticed the need for a third Stokes line at each crossing point.
32 Activity of Stokes lines x-plane i d3 y dy + 3i dx3 dx + xy = 0 Im x Re x
33 Activity of Stokes lines x-plane i d3 y dy + 3i dx3 dx + xy = 0 y exp ( x ) i k j (ξ)dξ x T P Im x Re x
34 Activity of Stokes lines x-plane i d3 y dy + 3i dx3 dx + xy = 0 ( x ) y exp i k j (ξ)dξ x T P Im x k 3 j 3k j + x = 0 Re x
35 Activity of Stokes lines x-plane { } Turning points i d3 y dy + 3i dx3 dx + xy = 0 ( x ) y exp i k j (ξ)dξ x T P k 3 j 3k j + x = 0 k i (x T P ) = k j (x T P )
36 Activity of Stokes lines x-plane Stokes lines i d3 y dy + 3i dx3 dx + xy = 0 ( x ) y exp( i k j (ξ)dξ ) x T P k 3 j 3k j + x = 0 Im { x } i (k i k j ) dx x T P = 0
37 Activity of Stokes lines x-plane i d3 y dy + 3i dx3 dx + xy = 0 ( x ) y exp( i k j (ξ)dξ ) x T P k 3 j 3k j + x = 0 Im { x } i (k i k j ) dx x T P = 0
38 Activity of Stokes lines x-plane i d3 y dy + 3i dx3 dx + xy = 0 ( x ) y exp( i k j (ξ)dξ ) x T P k 3 j 3k j + x = 0 Im { x } i (k i k j ) dx x T P = 0
39 Activity of Stokes lines x-plane i d3 y dy + 3i dx3 dx + xy = 0 ( x ) y exp( i k j (ξ)dξ ) x T P k 3 j 3k j + x = 0 Im { x } i (k i k j ) dx x T P = 0
40 Activity of Stokes lines x-plane i d3 y dy + 3i dx3 dx + xy = 0 ( x ) y exp( i k j (ξ)dξ ) x T P k 3 j 3k j + x = 0 Im { x } i (k i k j ) dx x T P = 0 New Stokes line
41 Activity of Stokes lines x-plane i d3 y dy + 3i dx3 dx + xy = 0 ( x ) y exp( i k j (ξ)dξ ) x T P k 3 j 3k j + x = 0 Im New Stokes line j = 1, 2, 3 { x } i (k i k j ) dx x T P = 0
42 Activity of Stokes lines x-plane i d3 y dy + 3i dx3 dx + xy = 0 ( x ) y exp( i k j (ξ)dξ ) x T P k 3 j 3k j + x = 0 2 > 3 1 > 2 Im { x } i (k i k j ) dx x T P = 0 New Stokes line j = 1, 2, 3
43 Activity of Stokes lines x-plane i d3 y dy + 3i dx3 dx + xy = 0 ( x ) y exp( i k j (ξ)dξ ) x T P k 3 j 3k j + x = 0 2 > 3 1 > 3 1 > 2 Im { x } i (k i k j ) dx x T P = 0 New Stokes line j = 1, 2, 3
44 Activity of Stokes lines x-plane i d3 y dy + 3i dx3 dx + xy = 0 ( x ) y exp( i k j (ξ)dξ ) x T P k 3 j 3k j + x = 0 Im { x } i (k i k j ) dx x T P Why change in activity of Stokes line inside crossing point? = 0
45 Activity of Stokes lines x-plane i d3 y dy + 3i dx3 dx + xy = 0 ( x ) y exp( i k j (ξ)dξ ) x T P 1 > 3 k 3 j 3k j + x = 0 Im { x } i (k i k j ) dx x T P Why change in activity of Stokes line inside crossing point? = 0
46 Activity of Stokes lines x-plane i d3 y dy + 3i dx3 dx + xy = 0 ( x ) y exp( i k j (ξ)dξ ) x T P k 3 j 3k j + x = 0 Im { x } i (k i k j ) dx x T P Virtual turning point: 1 & 3 coalesce, but regular point. = 0
47 Activity of Stokes lines BNR explained in terms of somewhat complicated integral representations of the specific equation. Here we show that change in activity of a Stokes line is a general feature of asymptotic expansions involving 3 or more asymptotic contributions with parameters. Not confined to complex plane: real effects possible. Start explanation with an integral.
48 An integral with 3 saddles ( ) = exp f ( z;a) ε I ε;a V 2 V 1 ( ) where V 1 = e 3πi 8, V 2 = e πi 8 f ( z) = i( 1 4 z z2 + az) ε <<1
49 An integral with 3 saddles where Not quite Pearcey, but found in uniform asymptotics ( ) = exp f ( z;a) ε I ε;a V 2 V 1 ( ) V 1 = e 3πi 8, V 2 = e πi 8 f ( z) = i( 1 4 z z2 + az) ε <<1 Control Parameter
50 An integral with 3 saddles ( ) = exp f ( z;a) ε I ε;a V 2 V 1 ( ) Saddles: z 0, z 1, z 2 z n 3 + zn + a = 0, n = 0,1, 2
51 An integral with 3 saddles ( ) = exp f ( z;a) ε I ε;a V 2 V 1 ( ) Saddles: z 0, z 1, z 2 z n 3 + zn + a = 0, n = 0,1, 2 Expansion: ( ) ~ exp f n ( a) ε I ε;a ( ) subset a n ( ) T r f n ( a) f ( z n ;a) = 1 4 iz n( z n + 3a ) r ( ) ε r
52 An integral with 3 saddles ( ) = exp f ( z;a) ε I ε;a V 2 V 1 ( ) Saddles: z 0, z 1, z 2 z n 3 + zn + a = 0, n = 0,1, 2 Expansion: ( ) ~ exp f n ( a) ε I ε;a ( ) subset a n ( ) T r f n ( a) f ( z n ;a) = 1 4 iz n( z n + 3a ) r ( ) ε r Subset changes when: ( ) f m ( a) { a : ( f n a ) ε > 0} Stokes phenomenon
53 Integration plane a-plane Borel plane
54 Stokes phenomenon Integration plane a-plane Borel plane
55 Stokes phenomenon Stokes line Integration plane a-plane Borel plane
56 Stokes phenomenon Saddles connect Integration plane a-plane Borel plane
57 Stokes phenomenon Integration plane a-plane Corresponding singularities line up horizontally Borel plane
58 Integration plane 3 a-plane Borel plane
59 Integration plane 3 a-plane Borel plane
60 Stokes phenomenon a 5 Integration plane a 5 3 a-plane Borel plane
61 a 6 Integration plane a 6 3 a-plane Borel plane
62 Stokes phenomenon a 7 Integration plane a 7 3 a-plane Borel plane
63 Stokes phenomenon a 7 Integration plane a 7 3 a-plane Borel plane
64 No Stokes phenomenon a 7 Integration plane a 7 3 a-plane Borel plane
65 Connection is between saddles that are are currently not contributing a 7 Integration plane a 7 3 a-plane Borel plane
66 Connection is between saddles that are are currently not contributing Irrelevant Stokes line a 7 Integration plane a 7 3 a-plane Borel plane
67 a-plane a 7 3
68 Stokes crossing point: a regular point a-plane
69 Stokes crossing point: a regular point Three in a row a 7 3 a-plane Borel plane
70 Stokes crossing point: a regular point a 7 3 Stokes line a-plane Borel plane
71 Stokes crossing point: a regular point a 7 3 No Stokes line Stokes line a-plane Borel plane
72 Stokes crossing point: a regular point even though singularity 2 is crossing cut from singularity 0: signature of Stokes phenomenon! a 7 3 No Stokes line Stokes line a-plane Borel plane
73 Stokes crossing point: a regular point Stokes line stops at regular point!?! a 7 3 No Stokes line Stokes line a-plane Borel plane
74 Eh?
75 Eh? How can we explain this?
76 Eh? How can we explain this? Higher order Stokes phenomenon!
77 Higher order Stokes line Integration plane a-plane Borel plane
78 Higher order Stokes line Integration plane a-plane Collinearity of singularities Borel plane
79 Higher order Stokes line Nothing happening to saddles! Integration plane a-plane Collinearity of singularities Borel plane
80 What does happen at a higher order Stokes phenomenon? Ordinary Stokes phenomenon: Exponentially small contribution is switched on or off Stokes curves: Higher order Stokes phenomenon ( ) > 0 { a :ε 1 f 1 ( a) f 0 ( a) } Possibility of a Stokes phenomenon is switched on/off A Stokes multipler changes value Higher order Stokes curves: ( f 1 ( a) f 0 ( a) ) a : f 2 ( a) f 1 ( a) ( ) > 0
81 What does happen at a higher order Stokes phenomenon? Ordinary Stokes phenomenon: Exponentially small contribution is switched on or off Stokes curves: Higher order Stokes phenomenon ( ) > 0 { a :ε 1 f 1 ( a) f 0 ( a) } Possibility of a Stokes phenomenon is switched on/off A Stokes multipler changes value Higher order Stokes curves: ( f 1 ( a) f 0 ( a) ) a : f 2 ( a) f 1 ( a) ( ) > 0
82 What does happen at a higher order Stokes phenomenon? Ordinary Stokes phenomenon: Exponentially small contribution is switched on or off Stokes curves: Higher order Stokes phenomenon ( ) > 0 { a :ε 1 f 1 ( a) f 0 ( a) } Possibility of a Stokes phenomenon is switched on/off A Stokes multipler itself changes value Higher order Stokes curves: ( f 1 ( a) f 0 ( a) ) a : f 2 ( a) f 1 ( a) Condition for collinearity ( ) > 0
83 What does happen at a Stokes phenomenon? Subdominant singularity crosses integration contour from dominant singularity.
84 What does happen at a Stokes phenomenon? Subdominant singularity crosses integration contour from dominant singularity. Integration contour loops around subdominant singularity.
85 What does happen at a Stokes phenomenon? Subdominant singularity crosses integration contour from dominant singularity. Integration contour loops around subdominant singularity. SUBDOMINANT CONTRIBUTION ACQUIRED.
86 What does happen at a higher order Stokes phenomenon?
87 What does happen at a higher order Stokes phenomenon? 0 can t see 2: Stokes constant =0
88 What does happen at a higher order Stokes phenomenon? 0 can t see 2: Stokes constant =0 0, 1 & 2 collinear
89 What does happen at a higher order Stokes phenomenon? 0 can t see 2: Stokes constant =0 0, 1 & 2 collinear 0 sees 2: Stokes constant 0
90 What does happen at a higher order Stokes phenomenon? 0 can t see 2: Stokes constant =0 0, 1 & 2 colinear 0 sees 2: Stokes constant 0 NB: NO CHANGE IN FORM OF ORDINARY ASYMPTOTICS!
91 What does happen at a higher order Stokes phenomenon? Singularity generating dominant contribution
92 What does happen at a higher order Stokes phenomenon? Singularity generating dominant contribution is on different sheet to subdominant singularity NB: NO CHANGE IN FORM OF ORDINARY ASYMPTOTICS!
93 Frequency of Higher Order Stokes Stokes line: 2 exponents line up horizontally in Borel plane. f 1 f 2 f 1 F 12 >0 Higher order Stokes line: 3 exponents line up in any direction in Borel plane. Same number of constraints (1): f 2 f 3 F 12 /F 23 >0 Same codimension Just as common
94 Higher Order Stokes Phenomenon Higher order Stokes phenomenon can occur whenever 3 or more different asymptotic contributions/independent Borel plane singularities exist. Inhomogeneous second order linear odes Higher order linear homogeneous odes Second order nonlinear odes Single/multiple integrals PDEs...
95 Ordinary vs. Higher order Stokes ε real { a :ε 1 ( f 1 ( a) f 0 ( a) ) > 0} ε complex
96 Ordinary vs. Higher order Stokes ε real { a :ε 1 ( f 1 ( a) f 0 ( a) ) > 0} ε complex
97 Ordinary vs. Higher order Stokes ε real ε complex Ordinary Stokes dependent on asymptotic parameter.
98 Ordinary vs. Higher order Stokes a : ( f 1 ( a) f 0 ( a) ) ε real f 2 ( a) f 1 ( a) ε complex Higher order Stokes independent of asymptotic parameter! ( ) > 0
99 Link to Hyperasymptotics Hyperasymptotics previously used to find Stokes constants. If Stokes constants are changing... Hyper and higher order Stokes phenomenon intimately linked. Examples shows how HoSP can be crucial for deducing larger time behaviour in evolution equations.
100 PDE illustrative example Derive large time behaviour of: ( ) θ t θ x = ε 2 θ xxx x 2 x (, + ), t 0, ( ) = arctan x [ ) θ x,0 θ xx,θ x,θ bounded as x
101 PDE illustrative example Dominant balances θ x ~1 ( 1 + x 2 ) θ( x) ~ arctan x
102 PDE illustrative example Dominant balances θ x ~1 ( 1 + x 2 ) θ( x) ~ arctan x Satisfies initial condition
103 PDE illustrative example Dominant balances θ x ~1 ( 1 + x 2 ) θ( x) ~ arctan x Satisfies initial condition θ x ~ ε 2 θ xxx θ ~ Asin ( x ε) + Bcos( x ε ) + C
104 PDE illustrative example Dominant balances θ x ~1 ( 1 + x 2 ) θ( x) ~ arctan x Satisfies initial condition θ x ~ ε 2 θ xxx θ ~ Asin ( x ε) + Bcos( x ε ) + C Doesn t satisfy initial condition, not obviously present in solution.
105 PDE illustrative example Can solve PDE by Fourier integral (cheating!) u(x, t; ε) = I 1 (x, t; ε) + I 2 (x, t; ε) + I 3 (x; ε) + I 4 (x; ε) I 1 (x, t; ε) = I 3 (x; ε) = 0 0 iπp 2 1 p 2 e (p(1+ix)+ip(1 p )t)/ε dp iπ p(1 p 2 ) e p(1+ix)/ε dp Endpoint, saddle point, pole asymptotics. BUT can derive same results directly from expansion coefficients.
106 PDE illustrative example Base (hyper)asymptotics on
107 PDE illustrative example Base (hyper)asymptotics on ( ) ~ arctan x + a n x,t θ x,t r=1 n ( )ε n
108 PDE illustrative example Base (hyper)asymptotics on ( ) ~ arctan x + a n x,t θ x,t r=1 ( )ε n Large time behaviour evolves from sum, via hyperasymptotics: 3 different types of contribution. n
109 PDE illustrative example θ s θ e ( x,t) = arctan x + 4πε 2 x + t 1 + x + t Always contributing ( x,t) = 2 Re ( 3t )1 4 π 3 2 i + t + x ( i + 2t x) θ p Sometimes contributing ( x,t ) = 2π 2 e 1 ε cos x ε ( ) ( ( ) 2 ) 2 x ( ) 1 4 exp ( 1 + x 2 ) 2 + O ε 3 2i ( i + t + x ) tε ( ) i π 4 Sub-sub dominant near t=0, but dominates large t behaviour!
110 PDE illustrative example θ s θ e ( x,t) = arctan x + 4πε 2 x + t 1 + x + t Always contributing ( x,t) = 2 Re ( 3t )1 4 π 3 2 i + t + x ( i + 2t x) θ p Sometimes contributing ( x,t ) = 2π 2 e 1 ε cos x ε ( ) ( ( ) 2 ) 2 x ( ) 1 4 exp ( 1 + x 2 ) 2 + O ε 3 2i ( i + t + x ) tε ( ) i π 4 Sub-sub dominant near t=0, but dominates large t behaviour!
111 PDE illustrative example θ s θ e ( x,t) = arctan x + 4πε 2 x + t 1 + x + t Always contributing ( x,t) = 2 Re ( 3t )1 4 π 3 2 i + t + x ( i + 2t x) θ p Sometimes contributing ( x,t ) = 2π 2 e 1 ε cos x ε ( ) ( ( ) 2 ) 2 x ( ) 1 4 exp ( 1 + x 2 ) 2 + O ε 3 2i ( i + t + x ) tε ( ) i π 4 Sub-sub dominant near t=0, but dominates large t behaviour!
112 Higher Order Stokes Phenomenon e+p+s e+s? e+p+s Contradiction: no active Stokes line on x = 0 for t > A e+p e+p A e e+s
113 Higher Order Stokes Phenomenon Switches off Stokes line ep Higher order Stokes line
114 Higher Order Stokes Phenomenon: Hyper explanation Ordinary Stokes phenomenon
115 Higher Order Stokes Phenomenon: Hyper explanation Ordinary Stokes phenomenon Birth of exponentially small terms across curve (surface) in (x,t) plane: switches on/off terms. Change in number of contributions to full asymptotic expansion.
116 Higher Order Stokes Phenomenon: Hyper explanation Ordinary Stokes phenomenon Birth of exponentially small terms across curve (surface) in (x,t) plane: switches on/off terms. Change in number of contributions to full asymptotic expansion. Higher order Stokes phenomenon
117 Higher Order Stokes Phenomenon: Hyper explanation Ordinary Stokes phenomenon Birth of exponentially small terms across curve (surface) in (x,t) plane: switches on/off terms. Change in number of contributions to full asymptotic expansion. Higher order Stokes phenomenon Change in value of Stokes constants across curve (surface) in (x,t) plane: switches on/off Stokes lines. Change in number of remainder terms.
118 Higher Order Stokes Phenomenon: Hyper explanation ( ( N y (1) (ɛ) e f1/ɛ ɛ µ (0) 1 n +ɛ N (0) 1 r=0 K 12 2πi K 13 2πi K 14 2πi a (1) r ɛ r N (1) 2 1 a (2) s=0 N (1) 3 1 a (3) s=0 N (1) 4 1 a (4) s=0 ( s F (1) ɛ; ) + K 21 2πi K 23 2πi K 24 2πi ( s F (1) ɛ; ; ) { + ( s F (1) ɛ; ; ) { + N (2) 1 1 a (1) s=0 N (2) 3 1 a (3) s=0 N (2) 4 1 a (4) s=0 ( s F (2) ɛ; ; ) + ( s F (2) ɛ; ; ) + ( s F (2) ɛ; ; ) + Tree structure: any function possessing a Borel transform
119 Higher Order Stokes Phenomenon: Hyper explanation θ e N 1 e ( x,t) = a r r=0 ( ) ε r p= p 1,p 2 Level 1 Hyper re-expansion ε N g p F ep N +1 du e u u N uε F ep 0 ( ) m=s 1,s 2 p= p 1, p 2 dv e v v N M 1 2 ( 1 vε F em ) 0 ε N g p 2πiF em N M+1 2 F mp M+1 2 du e u u M 1/2 ( 1 uf em vf mp ) +third order hyperterms 0 F ep Responsible For Stokes line ep, when F ep >0 ( x,t), F mp ( x,t), F em ( x,t)
120 θ e Higher Order Stokes Phenomenon: Hyper explanation N 1 e ( x,t) = a r r=0 m=s 1,s 2 ( ) ε r p= p 1, p 2 p= p 1,p 2 dv e v v N M 1 2 ( 1 vε F em ) 0 ε N g p F ep N +1 du e u u N uε F ep 0 ( ) ε N g p 2πiF em N M+1 2 F mp M+1 2 du e u u M 1/2 ( 1 uf em vf mp ) +third order hyperterms 0 Level 1 Hyper re-expansion Residue cancels F ep term exactly Level 2 Hyper re-expansion Pole when F em /F mp > 0
121 Higher Order Stokes Phenomenon: Hyper explanation θ e N 1 e ( x,t) = a r r=0 ( ) ε r m=s 1,s 2 p= p 1, p 2 dv e v v N M 1 2 ( 1 vε F em ) 0 ε N g p 2πiF em N M+1 2 F mp M+1 2 du e u u M 1/2 ( 1 uf em vf mp ) +third order hyperterms 0 Residue cancels F ep term exactly Pole when F em /F mp > 0 Condition for colinearity of singularities in the Borel plane
122 Higher Order Stokes Phenomenon: Hyper explanation θ e N 1 e ( x,t) = a r r=0 m=s 1,s 2 ( ) ε r p= p 1, p 2 dv e v v N M 1 2 ( 1 vε F em ) 0 ε N g p 2πiF em N M+1 2 F mp M+1 2 du e u u M 1/2 ( 1 uf em vf mp ) 0 No further possibility of Stokes line between e and p: Higher order Stokes Phenomenon +third order hyperterms
123 PDE illustrative example t x Exact ε=0.5
124 PDE illustrative example NB t ε=0.5 x Asymptotics neglecting initial sub-subdominant term
125 PDE illustrative example NB t x Exact ε=0.5
126 PDE illustrative example NB t ε=0.5 x Asymptotics including initially sub-sub dominant term and taking account of higher order Stokes phenomenon
127 PDE illustrative example Exact θ x
128 PDE illustrative example Exact With HoSP θ x
129 PDE illustrative example Exact With HoSP HoSP ignored θ x
130 PDE illustrative example Exact With HoSP HoSP ignored Exact vs HoSP θ x
131 PDE illustrative example Exact With HoSP HoSP ignored Exact vs HoSP θ x Neighborhood of complex caustics
132 Hyper-tree and higher order Stokes ( y (1) (ɛ) e f 1/ɛ ɛ µ n +ɛ N (0) 1 ( N (0) 1 r=0 K 12 2πi K 13 2πi K 14 2πi a (1) r ɛ r N (1) 2 1 a (2) s=0 N (1) 3 1 a (3) s=0 N (1) 4 1 a (4) s=0 ( s F (1) ɛ; ) + K 21 2πi K 23 2πi K 24 2πi ( s F (1) ɛ; ; ) { + ( s F (1) ɛ; ; ) { + N (2) 1 1 a (1) s=0 N (2) 3 1 a (3) s=0 N (2) 4 1 a (4) s=0 ( s F (2) ɛ; ; ) + ( s F (2) ɛ; ; ) + ( s F (2) ɛ; ; ) + Tree structure: any function possessing a Borel transform
133 Hyper-tree and higher order Stokes ( y (1) (ɛ) e f 1/ɛ ɛ µ n +ɛ N (0) 1 ( N (0) 1 r=0 K 12 2πi K 13 2πi K 14 2πi a (1) r ɛ r N (1) 2 1 a (2) s=0 N (1) 3 1 a (3) s=0 N (1) 4 1 a (4) s=0 ( s F (1) ɛ; ) + K 21 2πi K 23 2πi K 24 2πi ( s F (1) ɛ; ; ) { + ( s F (1) ɛ; ; ) { + N (2) 1 1 a (1) s=0 N (2) 3 1 a (3) s=0 N (2) 4 1 a (4) s=0 ( s F (2) ɛ; ; ) + ( s F (2) ɛ; ; ) + ( s F (2) ɛ; ; ) + Higher order twig on one branch...
134 Hyper-tree and higher order Stokes ( y (1) (ɛ) e f 1/ɛ ɛ µ n +ɛ N (0) 1 ( N (0) 1 r=0 K 12 2πi K 13 2πi K 14 2πi a (1) r ɛ r N (1) 2 1 a (2) s=0 N (1) 3 1 a (3) s=0 N (1) 4 1 a (4) s=0 ( s F (1) ɛ; ) + K 21 2πi K 23 2πi K 24 2πi ( s F (1) ɛ; ; ) { + ( s F (1) ɛ; ; ) { + N (2) 1 1 a (1) s=0 N (2) 3 1 a (3) s=0 N (2) 4 1 a (4) s=0 ( s F (2) ɛ; ; ) + ( s F (2) ɛ; ; ) + ( s F (2) ɛ; ; ) +...can switch off other branches
135 Hyper-tree and higher order Stokes ( y (1) (ɛ) e f 1/ɛ ɛ µ n +ɛ N (0) 1 ( N (0) 1 r=0 K 12 2πi K 13 2πi K 14 2πi a (1) r ɛ r N (1) 2 1 a (2) s=0 N (1) 3 1 a (3) s=0 N (1) 4 1 a (4) s=0 ( s F (1) ɛ; ) + K 21 2πi K 23 2πi K 24 2πi ( s F (1) ɛ; ; ) { + ( s F (1) ɛ; ; ) { + N (2) 1 1 a (1) s=0 N (2) 3 1 a (3) s=0 N (2) 4 1 a (4) s=0 ( s F (2) ɛ; ; ) + ( s F (2) ɛ; ; ) + ( s F (2) ɛ; ; ) +...can switch off other branches
136 Hyper-tree and higher order Stokes ( y (1) (ɛ) e f 1/ɛ ɛ µ n +ɛ N (0) 1 ( N (0) 1 r=0 K 12 2πi K 13 2πi K 14 2πi a (1) r ɛ r N (1) 2 1 a (2) s=0 N (1) 3 1 a (3) s=0 N (1) 4 1 a (4) s=0 ( s F (1) ɛ; ) + K 21 2πi K 23 2πi K 24 2πi ( s F (1) ɛ; ; ) { + ( s F (1) ɛ; ; ) { + N (2) 1 1 a (1) s=0 N (2) 3 1 a (3) s=0 N (2) 4 1 a (4) s=0 ( s F (2) ɛ; ; ) + ( s F (2) ɛ; ; ) + ( s F (2) ɛ; ; ) + Equivalent to Stokes constant changing, e.g., vanishing.
137 Extensions to nonlinear PDEs Higher order Stokes phenomenon plays a significant role in smoothed-shock formation. Without it the smoothed-shock would be a caustic! Example: Burgers
138 Burgers Equation u t + uu x = εu xx, x C, t 0, ε 0 +. u(x, 0) = 1, and u 0 as x. 1 + x2 tion of the form
139 Burgers Equation ( ) x = x j + a 0 (x j )t, j = 0, 1, 2, rves are shown in Figure Inviscid multivaluedness
140 Burgers Equation x = x j + a 0 (x j )t, j = 0, 1, 2, 10 Inviscid shock
141 Burgers Equation u t + uu x = εu xx, x C, t 0, ε 0 +. u(x, 0) = 1, and u 0 as x. 1 + x2 tion of the form u(x, t; ε) u (0) ( ) (x, t; ε) + K( 01u n (n,1) (x, t; ε), ) Transseries n=1 ( ( ) ( ) ( ) a (n,1) 0 = a (1,1) n 0 2 f ) 1 n 1. x a (1,1) 0 (x 0, x 1 ) = (a 0 (x 1 ) a 0 (x 0 )) we find = = 0, so that u(x, t; ε) = a 0 (x, t) + 2K 01a (1,1) 0 (x, t) f 1 a 0 (x 1 ) a 0 (x 0 ) a 0 (x 0)(x 1 x 0 ) a 0 (x 1 ) a 0 (x 0 ) a 0 (x 1)(x 1 x 0 ). x e f 1/ε 2 f 1 x + K 01a (1,1) 0 (x, t)e f 1/ε + O(ε), result (38) is valid everywhere in the region where Re 0 and
142 Complex x-plane at t = 0.25, =0.05 ɛ 2 AS 01 Pre-shock S 0>1' C 01' 10 1 S 0>1 8 6 C 01 S 0> S 0>2' Active Stokes Anti-Stokes Higher-Order-Stokes Inactive Stokes Branch Cut Irrelevant Stokes curves not shown S 0>2 C 02 S 0>2-2 C 02' AS
143 Complex x-plane at t = 0.25, =0.05 ɛ 2 AS 01 Pre-shock S 0>1' C 01' 1 S 0>1 0-1 S 0>2' Active Stokes Anti-Stokes Higher-Order-Stokes Inactive Stokes Branch Cut Irrelevant Stokes curves not shown S 0>2 C 01 C 02 S 0>1 S 0>2 Inactive Stokes line: -2 C 02' AS
144 Complex x-plane at t = 0.25, =0.05 ɛ 2 AS 01 Pre-shock S 0>1' C 01' 1 S 0>1 C 01 S 0>1 0 Active Stokes Anti-Stokes Higher-Order-Stokes Inactive Stokes Branch Cut Irrelevant Stokes curves not shown C 02 S 0>2 Inactive Stokes line: switched off by Higher order Stokes line -1 S 0>2' S 0>2 C 02' -2 AS
145 2 Complex x-plane at t = 0.25, ɛ=0.05 AS 01 Pre-shock 1 S 0>1' C 01' S 0> C 01 S 0> C 02 S 0>2-1 S 0>2' S 0>2-2 C 02' AS
146 2 Complex x-plane at t = 0.25, ɛ=0.05 AS 01 Pre-shock S 0>1' C 01' 1 S 0>1 0 Without HOS a Stokes phenomenon would take place across real space even before shock formation. C 01 C 02 S 0>1 S 0>2-1 S 0>2' S 0>2-2 C 02' AS
147 Complex x-plane at t = 5, ɛ= C 01' AS 01 Post-shock t S 0>1 S C R C R x s C R S 0>1 S 0>2 x S 0> C 02' AS 02 Poles given by summed transseries
148 Complex x-plane at t = 5, ɛ=0.05 AS C 01' S 0>1 S 021 Shock=virtual turning point, not a caustic 0 C R x s C R S 0>1 S 0>2 S 0> C 02' AS 02
149 Complex x-plane at t = 5, ɛ=0.05 AS C 01' S 0>1 S 021 Shock=virtual turning point, not a caustic C 02' S 0>2 C R x s C R S 0>1 S 0>2 AS 02 HOSP forces asymptotic contributions onto different Riemann sheets, preventing coalescence and hence a caustic/ turning point
150 Generalisations u t + u 2 u x = εu xx, u t + F (u)u x = εu xx, etc...
151 Hokusai
152 Higher order Stokes curve
153 Virtual turning point
154 Conclusions Stokes lines switch off at regular points! Higher order Stokes phenomenon can explain activity of Stokes lines. Has crucial role role to play in asymptotic evolution of system, especially PDEs. Should include sub-sub(-sub-sub...) dominant terms! Simple hyperasymptotics can cope with this
155 KdV u t + uu x + ɛ 2 u xxx = 0 u(x, 0) = x 2 ɛ = 0.02
156 KdV t =5,ε =0.02 "kdv4fout.dat" u x Work in progress...
Exponentially accurate solution tracking for nonlinear ODEs, the higher order Stokes phenomenon and double transseries resummation
Exponentially accurate solution tracking for nonlinear ODEs, the higher order Stokes phenomenon and double transseries resummation C.J. Howls and A. B. Olde Daalhuis 2 School of Mathematics, University
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 informationPower series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) a n z n. n=0
Lecture 22 Power series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) Recall a few facts about power series: a n z n This series in z is centered at z 0. Here z can
More informationOn positivity of Fourier transforms
On positivity of Fourier transforms by E.O. Tuck Applied Mathematics The University of Adelaide AUSTRALIA 55 April 1, 26 Abstract This note concerns Fourier transforms on the real positive line. In particular,
More informationComplex Variables. Instructions Solve any eight of the following ten problems. Explain your reasoning in complete sentences to maximize credit.
Instructions Solve any eight of the following ten problems. Explain your reasoning in complete sentences to maximize credit. 1. The TI-89 calculator says, reasonably enough, that x 1) 1/3 1 ) 3 = 8. lim
More informationMath Homework 2
Math 73 Homework Due: September 8, 6 Suppose that f is holomorphic in a region Ω, ie an open connected set Prove that in any of the following cases (a) R(f) is constant; (b) I(f) is constant; (c) f is
More information2 Write down the range of values of α (real) or β (complex) for which the following integrals converge. (i) e z2 dz where {γ : z = se iα, < s < }
Mathematical Tripos Part II Michaelmas term 2007 Further Complex Methods, Examples sheet Dr S.T.C. Siklos Comments and corrections: e-mail to stcs@cam. Sheet with commentary available for supervisors.
More informationPartial Differential Equations
M3M3 Partial Differential Equations Solutions to problem sheet 3/4 1* (i) Show that the second order linear differential operators L and M, defined in some domain Ω R n, and given by Mφ = Lφ = j=1 j=1
More informationSolutions to practice problems for the final
Solutions to practice problems for the final Holomorphicity, Cauchy-Riemann equations, and Cauchy-Goursat theorem 1. (a) Show that there is a holomorphic function on Ω = {z z > 2} whose derivative is z
More informationFINAL EXAM MATH303 Theory of Ordinary Differential Equations. Spring dx dt = x + 3y dy dt = x y.
FINAL EXAM MATH0 Theory of Ordinary Differential Equations There are 5 problems on 2 pages. Spring 2009. 25 points Consider the linear plane autonomous system x + y x y. Find a fundamental matrix of the
More information{ } is an asymptotic sequence.
AMS B Perturbation Methods Lecture 3 Copyright by Hongyun Wang, UCSC Recap Iterative method for finding asymptotic series requirement on the iteration formula to make it work Singular perturbation use
More informationHamiltonian partial differential equations and Painlevé transcendents
The 6th TIMS-OCAMI-WASEDA Joint International Workshop on Integrable Systems and Mathematical Physics March 22-26, 2014 Hamiltonian partial differential equations and Painlevé transcendents Boris DUBROVIN
More informationDivergent series: taming the tails
Divergent series: taming the tails M. V. Berry & C. J. Howls 2 H. H. Wills Physics Laboratory, Tyndall Avenue, Bristol BS8 TL, UK 2 Mathematical Sciences, University of Southampton, Southampton, SO7 BJ,
More informationMath 4263 Homework Set 1
Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that
More informationSynopsis of Complex Analysis. Ryan D. Reece
Synopsis of Complex Analysis Ryan D. Reece December 7, 2006 Chapter Complex Numbers. The Parts of a Complex Number A complex number, z, is an ordered pair of real numbers similar to the points in the real
More informationBessel s Equation. MATH 365 Ordinary Differential Equations. J. Robert Buchanan. Fall Department of Mathematics
Bessel s Equation MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Fall 2018 Background Bessel s equation of order ν has the form where ν is a constant. x 2 y + xy
More informationBorel Summability in PDE initial value problems
Borel Summability in PDE initial value problems Saleh Tanveer (Ohio State University) Collaborator Ovidiu Costin & Guo Luo Research supported in part by Institute for Math Sciences (IC), EPSRC & NSF. Main
More information221A Lecture Notes Steepest Descent Method
Gamma Function A Lecture Notes Steepest Descent Method The best way to introduce the steepest descent method is to see an example. The Stirling s formula for the behavior of the factorial n! for large
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 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 informationMATH 452. SAMPLE 3 SOLUTIONS May 3, (10 pts) Let f(x + iy) = u(x, y) + iv(x, y) be an analytic function. Show that u(x, y) is harmonic.
MATH 45 SAMPLE 3 SOLUTIONS May 3, 06. (0 pts) Let f(x + iy) = u(x, y) + iv(x, y) be an analytic function. Show that u(x, y) is harmonic. Because f is holomorphic, u and v satisfy the Cauchy-Riemann equations:
More informationSelected Solutions To Problems in Complex Analysis
Selected Solutions To Problems in Complex Analysis E. Chernysh November 3, 6 Contents Page 8 Problem................................... Problem 4................................... Problem 5...................................
More informationA REVIEW OF RESIDUES AND INTEGRATION A PROCEDURAL APPROACH
A REVIEW OF RESIDUES AND INTEGRATION A PROEDURAL APPROAH ANDREW ARHIBALD 1. Introduction When working with complex functions, it is best to understand exactly how they work. Of course, complex functions
More informationSuggested Homework Solutions
Suggested Homework Solutions Chapter Fourteen Section #9: Real and Imaginary parts of /z: z = x + iy = x + iy x iy ( ) x iy = x #9: Real and Imaginary parts of ln z: + i ( y ) ln z = ln(re iθ ) = ln r
More informationMath Final Exam.
Math 106 - Final Exam. This is a closed book exam. No calculators are allowed. The exam consists of 8 questions worth 100 points. Good luck! Name: Acknowledgment and acceptance of honor code: Signature:
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 information1 What s the big deal?
This note is written for a talk given at the graduate student seminar, titled how to solve quantum mechanics with x 4 potential. What s the big deal? The subject of interest is quantum mechanics in an
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 informationLocal smoothing and Strichartz estimates for manifolds with degenerate hyperbolic trapping
Local smoothing and Strichartz estimates for manifolds with degenerate hyperbolic trapping H. Christianson partly joint work with J. Wunsch (Northwestern) Department of Mathematics University of North
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 informationMATH 425, FINAL EXAM SOLUTIONS
MATH 425, FINAL EXAM SOLUTIONS Each exercise is worth 50 points. Exercise. a The operator L is defined on smooth functions of (x, y by: Is the operator L linear? Prove your answer. L (u := arctan(xy u
More informationApplied Asymptotic Analysis
Applied Asymptotic Analysis Peter D. Miller Graduate Studies in Mathematics Volume 75 American Mathematical Society Providence, Rhode Island Preface xiii Part 1. Fundamentals Chapter 0. Themes of Asymptotic
More informationPhysics 307. Mathematical Physics. Luis Anchordoqui. Wednesday, August 31, 16
Physics 307 Mathematical Physics Luis Anchordoqui 1 Bibliography L. A. Anchordoqui and T. C. Paul, ``Mathematical Models of Physics Problems (Nova Publishers, 2013) G. F. D. Duff and D. Naylor, ``Differential
More informationLECTURE-13 : GENERALIZED CAUCHY S THEOREM
LECTURE-3 : GENERALIZED CAUCHY S THEOREM VED V. DATAR The aim of this lecture to prove a general form of Cauchy s theorem applicable to multiply connected domains. We end with computations of some real
More informationSalmon: Lectures on partial differential equations
4 Burger s equation In Lecture 2 we remarked that if the coefficients in u x, y,! "! "x + v x,y,! "! "y = 0 depend not only on x,y but also on!, then the characteristics may cross and the solutions become
More informationGreen s Theorem, Cauchy s Theorem, Cauchy s Formula
Math 425 Spring 2003 Green s Theorem, Cauchy s Theorem, Cauchy s Formula These notes supplement the discussion of real line integrals and Green s Theorem presented in.6 of our text, and they discuss applications
More informationExercises for Part 1
MATH200 Complex Analysis. Exercises for Part Exercises for Part The following exercises are provided for you to revise complex numbers. Exercise. Write the following expressions in the form x + iy, x,y
More informationEE2 Mathematics : Complex Variables
EE Mathematics : omplex Variables J. D. Gibbon (Professor J. D Gibbon 1, Dept of Mathematics) j.d.gibbon@ic.ac.uk http://www.imperial.ac.uk/ jdg These notes are not identical word-for-word with my lectures
More informationGeneral Technical Remarks on PDE s and Boundary Conditions Kurt Bryan MA 436
General Technical Remarks on PDE s and Boundary Conditions Kurt Bryan MA 436 1 Introduction You may have noticed that when we analyzed the heat equation on a bar of length 1 and I talked about the equation
More informationIntroduction and some preliminaries
1 Partial differential equations Introduction and some preliminaries A partial differential equation (PDE) is a relationship among partial derivatives of a function (or functions) of more than one variable.
More informationOn Hamiltonian perturbations of hyperbolic PDEs
Bologna, September 24, 2004 On Hamiltonian perturbations of hyperbolic PDEs Boris DUBROVIN SISSA (Trieste) Class of 1+1 evolutionary systems w i t +Ai j (w)wj x +ε (B i j (w)wj xx + 1 2 Ci jk (w)wj x wk
More informationMAE294B/SIOC203B: Methods in Applied Mechanics Winter Quarter sgls/mae294b Solution IV
MAE9B/SIOC3B: Methods in Applied Mechanics Winter Quarter 8 http://webengucsdedu/ sgls/mae9b 8 Solution IV (i The equation becomes in T Applying standard WKB gives ɛ y TT ɛte T y T + y = φ T Te T φ T +
More informationRepresentations of solutions of general nonlinear ODEs in singular regions
Representations of solutions of general nonlinear ODEs in singular regions Rodica D. Costin The Ohio State University O. Costin, RDC: Invent. 2001, Kyoto 2000, [...] 1 Most analytic differential equations
More informationStokes-multiplier expansion in an inhomogeneous differential. equation with a small parameter
Stokes-multiplier expansion in an inhomogeneous differential equation with a small parameter E. I. Ólafsdóttir, A. B. Olde Daalhuis and J. Vanneste School of Mathematics, University of Edinburgh, Edinburgh
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 informationFunctions of a Complex Variable and Integral Transforms
Functions of a Complex Variable and Integral Transforms Department of Mathematics Zhou Lingjun Textbook Functions of Complex Analysis with Applications to Engineering and Science, 3rd Edition. A. D. Snider
More informationLeplace s Equations. Analyzing the Analyticity of Analytic Analysis DEPARTMENT OF ELECTRICAL AND COMPUTER ENGINEERING. Engineering Math EECE
Leplace s Analyzing the Analyticity of Analytic Analysis Engineering Math EECE 3640 1 The Laplace equations are built on the Cauchy- Riemann equations. They are used in many branches of physics such as
More information2-Motivations From the Theory of ODE s MATH 22C
2-Motivations From the Theory of ODE s MATH 22C 1. Ordinary Differential Equations (ODE) and the Fundamental Role of the Derivative in the Sciences Recall that a real valued function of a real variable
More information13 Definite integrals
3 Definite integrals Read: Boas h. 4. 3. Laurent series: Def.: Laurent series (LS). If f() is analytic in a region R, then converges in R, with a n = πi f() = a n ( ) n + n= n= f() ; b ( ) n+ n = πi b
More informationFirst-Order Ordinary Differntial Equations: Classification and Linear Equations. David Levermore Department of Mathematics University of Maryland
First-Order Ordinary Differntial Equations: Classification and Linear Equations David Levermore Department of Mathematics University of Maryland 1 February 2009 These notes cover some of the material that
More informationPropagation of discontinuities in solutions of First Order Partial Differential Equations
Propagation of discontinuities in solutions of First Order Partial Differential Equations Phoolan Prasad Department of Mathematics Indian Institute of Science, Bangalore 560 012 E-mail: prasad@math.iisc.ernet.in
More informationAdditional material: Linear Differential Equations
Chapter 5 Additional material: Linear Differential Equations 5.1 Introduction The material in this chapter is not formally part of the LTCC course. It is included for completeness as it contains proofs
More informationComplex Analysis Math 205A, Winter 2014 Final: Solutions
Part I: Short Questions Complex Analysis Math 205A, Winter 2014 Final: Solutions I.1 [5%] State the Cauchy-Riemann equations for a holomorphic function f(z) = u(x,y)+iv(x,y). The Cauchy-Riemann equations
More information2. Complex Analytic Functions
2. Complex Analytic Functions John Douglas Moore July 6, 2011 Recall that if A and B are sets, a function f : A B is a rule which assigns to each element a A a unique element f(a) B. In this course, we
More informationPartial Differential Equations Separation of Variables. 1 Partial Differential Equations and Operators
PDE-SEP-HEAT-1 Partial Differential Equations Separation of Variables 1 Partial Differential Equations and Operators et C = C(R 2 ) be the collection of infinitely differentiable functions from the plane
More informationAsymptotic Expansions
Asymptotic Expansions V 1.5.4 217/1/3 Suppose we have a function f(x) of single real parameter x and we are interested in an approximation to f(x) for x close to x. The approximation might depend on whether
More informationSolution for Final Review Problems 1
Solution for Final Review Problems Final time and location: Dec. Gymnasium, Rows 23, 25 5, 2, Wednesday, 9-2am, Main ) Let fz) be the principal branch of z i. a) Find f + i). b) Show that fz )fz 2 ) λfz
More information2.20 Fall 2018 Math Review
2.20 Fall 2018 Math Review September 10, 2018 These notes are to help you through the math used in this class. This is just a refresher, so if you never learned one of these topics you should look more
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 informationComplex functions, single and multivalued.
Complex functions, single and multivalued. D. Craig 2007 03 27 1 Single-valued functions All of the following functions can be defined via power series which converge for the entire complex plane, and
More informationAsymptotic Quasinormal Frequencies for d Dimensional Black Holes
Asymptotic Quasinormal Frequencies for d Dimensional Black Holes José Natário (Instituto Superior Técnico, Lisbon) Based on hep-th/0411267 with Ricardo Schiappa Oxford, February 2009 Outline What exactly
More informationUniformly accurate averaging numerical schemes for oscillatory evolution equations
Uniformly accurate averaging numerical schemes for oscillatory evolution equations Philippe Chartier University of Rennes, INRIA Joint work with M. Lemou (University of Rennes-CNRS), F. Méhats (University
More informationQuantum Physics Lecture 8
Quantum Physics ecture 8 Steady state Schroedinger Equation (SSSE): eigenvalue & eigenfunction particle in a box re-visited Wavefunctions and energy states normalisation probability density Expectation
More information1 Holomorphic functions
Robert Oeckl CA NOTES 1 15/09/2009 1 1 Holomorphic functions 11 The complex derivative The basic objects of complex analysis are the holomorphic functions These are functions that posses a complex derivative
More informationWave propagation in an inhomogeneous plasma
DRAFT Wave propagation in an inhomogeneous plasma Felix I. Parra Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Oxford OX NP, UK This version is of 7 February 208. Introduction In
More informationDifferential equations, comprehensive exam topics and sample questions
Differential equations, comprehensive exam topics and sample questions Topics covered ODE s: Chapters -5, 7, from Elementary Differential Equations by Edwards and Penney, 6th edition.. Exact solutions
More information(1) Let f(z) be the principal branch of z 4i. (a) Find f(i). Solution. f(i) = exp(4i Log(i)) = exp(4i(π/2)) = e 2π. (b) Show that
Let fz be the principal branch of z 4i. a Find fi. Solution. fi = exp4i Logi = exp4iπ/2 = e 2π. b Show that fz fz 2 fz z 2 fz fz 2 = λfz z 2 for all z, z 2 0, where λ =, e 8π or e 8π. Proof. We have =
More informationGlobal reconstruction of analytic functions from local expansions and a new general method of converting sums into integrals
Global reconstruction of analytic functions from local expansions and a new general method of converting sums into integrals O. Costin 1 Math Tower, 231 West 18th Avenue, Columbus, OH 4321 e-mail: costin@math.ohio-state.edu;
More informationMath 312 Fall 2013 Final Exam Solutions (2 + i)(i + 1) = (i 1)(i + 1) = 2i i2 + i. i 2 1
. (a) We have 2 + i i Math 32 Fall 203 Final Exam Solutions (2 + i)(i + ) (i )(i + ) 2i + 2 + i2 + i i 2 3i + 2 2 3 2 i.. (b) Note that + i 2e iπ/4 so that Arg( + i) π/4. This implies 2 log 2 + π 4 i..
More informationComplex Differentials and the Stokes, Goursat and Cauchy Theorems
Complex Differentials and the Stokes, Goursat and Cauchy Theorems Benjamin McKay June 21, 2001 1 Stokes theorem Theorem 1 (Stokes) f(x, y) dx + g(x, y) dy = U ( g y f ) dx dy x where U is a region of the
More informationLecture Introduction
Lecture 1 1.1 Introduction The theory of Partial Differential Equations (PDEs) is central to mathematics, both pure and applied. The main difference between the theory of PDEs and the theory of Ordinary
More informationThe first order quasi-linear PDEs
Chapter 2 The first order quasi-linear PDEs The first order quasi-linear PDEs have the following general form: F (x, u, Du) = 0, (2.1) where x = (x 1, x 2,, x 3 ) R n, u = u(x), Du is the gradient of u.
More informationDepartment of Mathematics, University of California, Berkeley. GRADUATE PRELIMINARY EXAMINATION, Part A Fall Semester 2016
Department of Mathematics, University of California, Berkeley YOUR 1 OR 2 DIGIT EXAM NUMBER GRADUATE PRELIMINARY EXAMINATION, Part A Fall Semester 2016 1. Please write your 1- or 2-digit exam number on
More informationLecture17: Generalized Solitary Waves
Lecture17: Generalized Solitary Waves Lecturer: Roger Grimshaw. Write-up: Andrew Stewart and Yiping Ma June 24, 2009 We have seen that solitary waves, either with a pulse -like profile or as the envelope
More informationMAT665:ANALYTIC FUNCTION THEORY
MAT665:ANALYTIC FUNCTION THEORY DR. RITU AGARWAL MALAVIYA NATIONAL INSTITUTE OF TECHNOLOGY JAIPUR Contents 1. About 2 2. Complex Numbers 2 3. Fundamental inequalities 2 4. Continuously differentiable functions
More informationPartial Differential Equations
Partial Differential Equations Partial differential equations (PDEs) are equations involving functions of more than one variable and their partial derivatives with respect to those variables. Most (but
More informationMATH 311: COMPLEX ANALYSIS CONTOUR INTEGRALS LECTURE
MATH 3: COMPLEX ANALYSIS CONTOUR INTEGRALS LECTURE Recall the Residue Theorem: Let be a simple closed loop, traversed counterclockwise. Let f be a function that is analytic on and meromorphic inside. Then
More informationWe start with a discussion of the Sampling Theorem.
Math 625 Lecture Notes - March 10, 2014 We start with a discussion of the Sampling Theorem. Discussion 0.1 Assume that the Fourier transform f of a function f satisfies f(w) = 0 for w > Ω, that is, f is
More information1. The COMPLEX PLANE AND ELEMENTARY FUNCTIONS: Complex numbers; stereographic projection; simple and multiple connectivity, elementary functions.
Complex Analysis Qualifying Examination 1 The COMPLEX PLANE AND ELEMENTARY FUNCTIONS: Complex numbers; stereographic projection; simple and multiple connectivity, elementary functions 2 ANALYTIC FUNCTIONS:
More informationSample Final Questions: Solutions Math 21B, Winter y ( y 1)(1 + y)) = A y + B
Sample Final Questions: Solutions Math 2B, Winter 23. Evaluate the following integrals: tan a) y y dy; b) x dx; c) 3 x 2 + x dx. a) We use partial fractions: y y 3 = y y ) + y)) = A y + B y + C y +. Putting
More informationOn the singularities of non-linear ODEs
On the singularities of non-linear ODEs Galina Filipuk Institute of Mathematics University of Warsaw G.Filipuk@mimuw.edu.pl Collaborators: R. Halburd (London), R. Vidunas (Tokyo), R. Kycia (Kraków) 1 Plan
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 informationMOMENTS OF HYPERGEOMETRIC HURWITZ ZETA FUNCTIONS
MOMENTS OF HYPERGEOMETRIC HURWITZ ZETA FUNCTIONS ABDUL HASSEN AND HIEU D. NGUYEN Abstract. This paper investigates a generalization the classical Hurwitz zeta function. It is shown that many of the properties
More informationMeromorphic connections and the Stokes groupoids
Meromorphic connections and the Stokes groupoids Brent Pym (Oxford) Based on arxiv:1305.7288 (Crelle 2015) with Marco Gualtieri and Songhao Li 1 / 34 Warmup Exercise Find the flat sections of the connection
More informationLeplace s Equations. Analyzing the Analyticity of Analytic Analysis DEPARTMENT OF ELECTRICAL AND COMPUTER ENGINEERING. Engineering Math 16.
Leplace s Analyzing the Analyticity of Analytic Analysis Engineering Math 16.364 1 The Laplace equations are built on the Cauchy- Riemann equations. They are used in many branches of physics such as heat
More informationUNIVERSITY OF SURREY FACULTY OF ENGINEERING AND PHYSICAL SCIENCES DEPARTMENT OF PHYSICS. BSc and MPhys Undergraduate Programmes in Physics LEVEL HE2
Phys/Level /1/9/Semester, 009-10 (1 handout) UNIVERSITY OF SURREY FACULTY OF ENGINEERING AND PHYSICAL SCIENCES DEPARTMENT OF PHYSICS BSc and MPhys Undergraduate Programmes in Physics LEVEL HE PAPER 1 MATHEMATICAL,
More informationComplex WKB analysis of energy-level degeneracies of non-hermitian Hamiltonians
INSTITUTE OF PHYSICS PUBLISHING JOURNAL OF PHYSICS A: MATHEMATICAL AND GENERAL J. Phys. A: Math. Gen. 4 (001 L1 L6 www.iop.org/journals/ja PII: S005-4470(01077-7 LETTER TO THE EDITOR Complex WKB analysis
More informationComplex varibles:contour integration examples
omple varibles:ontour integration eamples 1 Problem 1 onsider the problem d 2 + 1 If we take the substitution = tan θ then d = sec 2 θdθ, which leads to dθ = π sec 2 θ tan 2 θ + 1 dθ Net we consider the
More informationSecond Midterm Exam Name: Practice Problems March 10, 2015
Math 160 1. Treibergs Second Midterm Exam Name: Practice Problems March 10, 015 1. Determine the singular points of the function and state why the function is analytic everywhere else: z 1 fz) = z + 1)z
More informationSolutions for Problem Set #5 due October 17, 2003 Dustin Cartwright and Dylan Thurston
Solutions for Problem Set #5 due October 17, 23 Dustin Cartwright and Dylan Thurston 1 (B&N 6.5) Suppose an analytic function f agrees with tan x, x 1. Show that f(z) = i has no solution. Could f be entire?
More informationLecture Notes Complex Analysis. Complex Variables and Applications 7th Edition Brown and Churchhill
Lecture Notes omplex Analysis based on omplex Variables and Applications 7th Edition Brown and hurchhill Yvette Fajardo-Lim, Ph.D. Department of Mathematics De La Salle University - Manila 2 ontents THE
More informationFinal: Solutions Math 118A, Fall 2013
Final: Solutions Math 118A, Fall 2013 1. [20 pts] For each of the following PDEs for u(x, y), give their order and say if they are nonlinear or linear. If they are linear, say if they are homogeneous or
More informationlaplace s method for ordinary differential equations
Physics 24 Spring 217 laplace s method for ordinary differential equations lecture notes, spring semester 217 http://www.phys.uconn.edu/ rozman/ourses/p24_17s/ Last modified: May 19, 217 It is possible
More informationAdvanced Mathematical Methods for Scientists and Engineers I
Carl M. Bender Steven A. Orszag Advanced Mathematical Methods for Scientists and Engineers I Asymptotic Methods and Perturbation Theory With 148 Figures Springer CONTENTS! Preface xiii PART I FUNDAMENTALS
More informationMat104 Fall 2002, Improper Integrals From Old Exams
Mat4 Fall 22, Improper Integrals From Old Eams For the following integrals, state whether they are convergent or divergent, and give your reasons. () (2) (3) (4) (5) converges. Break it up as 3 + 2 3 +
More informationPower series and Taylor series
Power series and Taylor series D. DeTurck University of Pennsylvania March 29, 2018 D. DeTurck Math 104 002 2018A: Series 1 / 42 Series First... a review of what we have done so far: 1 We examined series
More informationMATH 205C: STATIONARY PHASE LEMMA
MATH 205C: STATIONARY PHASE LEMMA For ω, consider an integral of the form I(ω) = e iωf(x) u(x) dx, where u Cc (R n ) complex valued, with support in a compact set K, and f C (R n ) real valued. Thus, I(ω)
More informationAnalyzability in the context of PDEs and applications
Analyzability in the context of PDEs and applications O. Costin Math Department, Rutgers University S. Tanveer Math Department, The Ohio State University January 30, 2004 Abstract We discuss the notions
More information( ), λ. ( ) =u z. ( )+ iv z
AMS 212B Perturbation Metho Lecture 18 Copyright by Hongyun Wang, UCSC Method of steepest descent Consider the integral I = f exp( λ g )d, g C =u + iv, λ We consider the situation where ƒ() and g() = u()
More informationApplied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.
Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R
More information