arxiv: v1 [math.ap] 2 Nov 2012
|
|
- Miles Glenn
- 5 years ago
- Views:
Transcription
1 PARAXIAL COUPLING OF PROPAGAING MODES IN HREE-DIMENSIONAL WAVEGUIDES WIH RANDOM BOUNDARIES LILIANA BORCEA AND JOSSELIN GARNIER arxiv:.468v math.ap Nov Abstract. We analyze long range wave propagation in three-dimensional random waveguides. he waves are trapped by top and bottom boundaries, but the medium is unbounded in the two remaining directions. We consider scalar waves, and motivated by applications in underwater acoustics, we take a pressure release boundary condition at the top surface and a rigid bottom boundary. he wave speed in the waveguide is known and smooth, but the top boundary has small random fluctuations that cause significant cumulative scattering of the waves over long distances of propagation. o quantify the scattering effects, we study the evolution of the random amplitudes of the waveguide modes. We obtain that in the long range limit they satisfy a system of paraxial equations driven by a Brownian field. We use this system to estimate three important mode-dependent scales: the scattering mean free path, the cross-range decoherence length and the decoherence frequency. Understanding these scales is important in imaging and communication problems, because they encode the cumulative scattering effects in the wave field measured by remote sensors. As an application of the theory, we analyze time reversal and coherent interferometric imaging in strong cumulative scattering regimes. Key words. Waveguides, random media, asymptotic analysis. AMS subject classifications. 76B5, 35Q99, 6F5.. Introduction. We study long range scalar acoustic wave propagation in a three-dimensional waveguide. he setup is illustrated in Figure., and it is motivated by problems in underwater acoustics. We denote by z R the range, the main direction of propagation of the waves. he medium is unbounded in the cross-range direction x R, but it is confined in depth y by two boundaries which trap the waves, thus creating the waveguide effect. he acoustic pressure field is denoted by pt, x, y, z, and it satisfies the wave equation x + y + z c y t pt, x, y, z = ft, x, yδz, y, x, z, x, z R, t >,. in a medium with wave speed cy. he excitation is due to a source located in the plane z =, emitting the pulse ft. he medium is quiescent before the source excitation, he bottom of the waveguide is assumed rigid pt, x, y, z =, t.. y pt, x, y =, z =,.3 and we take a pressure release boundary condition at the perturbed top boundary pt, x, y = x, z, z =..4 Perturbed means that the boundary y = x, z has small fluctuations around the mean depth D, x, z D D..5 We choose this setup for simplicity. he results extend readily to other boundary conditions and to fluctuating bottoms. Such boundaries were considered recently in,, in two-dimensional waveguides. Extensions to media with small x, y, z-dependent fluctuations of the wave speed can also be made using the techniques developed in, 6, 8, 9,. he goal of our study is to quantify the effect of scattering at the surface. Because in applications it is not feasible to know the boundary fluctuations in detail, we model them with a random process. he solution pt, x, y, z of equations.-.4 is therefore a random field, and we describe in detail its statistics at long ranges, where cumulative scattering is significant. We use the results for two applications: time reversal and sensor array imaging. Computational and Applied Mathematics, Rice University, Houston, X 775. borcea@rice.edu Laboratoire de Probabilités et Modèles Aléatoires & Laboratoire Jacques-Louis Lions, Université Paris VII, Site Chevaleret, 755 Paris Cedex 3, France. garnier@math.jussieu.fr
2 y y = D x y = Fig... Schematic of the problem setup. he system of coordinates has range origin z = at the source. he rigid bottom boundary y = is assumed flat and the pressure release top boundary has fluctuations around the value y = D. he cross-range x and the range z are unbounded, that is x, z R. z Our method of solution uses a change of coordinates to straighten the boundary. he transformed problem has a simple geometry but a randomly perturbed differential operator. Its solution is given by a superposition of propagating and evanescent waveguide modes, with random amplitudes. We show that in the long range limit these amplitudes satisfy a system of paraxial equations that are driven by a Brownian field. he detailed characterization of the statistics of pt, x, y, z follows from this system. It involves the calculation of the mode-dependent scattering mean free path, which is the distance over which the modes lose coherence; the mode-dependent decoherence length, which is the cross-range offset over which the mode amplitudes decorrelate; and the mode-dependent decoherence frequency, which is the frequency offset over which the mode amplitudes decorrelate. hese scales are important in studies of time reversal and imaging, because they dictate the resolution of focusing and the robustness statistical stability of the results with respect to realizations of the random fluctuations of the boundary. he paper is organized as follows: We begin in section with the description of the reference pressure field p o t, x, y, z in ideal waveguides with planar boundaries. he random field pt, x, y, z derived in section 3 may be viewed as a perturbation of p o t, x, y, z, in the sense that it is decomposed in the same waveguide modes. However, the amplitudes of the modes are random and coupled. Because the fluctuations of the boundary are small, we consider in section 4 a long range limit, so that we can observe significant cumulative scattering. he statistics of the wave field at such long ranges is described in section 5. he results are summarized in section 6 and are used in sections 7 and 8 for analyzing time reversal and imaging with sensor arrays. We end with a summary in section 9.. Wave propagation in ideal waveguides. he pressure field in ideal waveguides, with planar boundaries, is given by p o t, x, y, z = dω π p oω, x, y, ze iωt,. with Fourier coefficients satisfying a separable problem for the Helmholtz equation x + y + z + ω c p o ω, x, y, z = y fω, x, yδz, ω ω B, x, z R, y, D,. with boundary conditions y p o ω, x, y =, z = p o ω, x, y = D, z =,.3 and outgoing radiation conditions at x + z. he Fourier transform of the source fω, x, y = dt ft, x, ye iωt,.4 is compactly supported in ω B/, ω + B/, for any x and y. Here ω is the central frequency and B is the bandwidth.
3 .. Propagating and evanescent modes. he solution of the Helmholtz equation. is a superposition of Nω propagating modes, and infinitely many evanescent ones, p o ω, x, y, z = Nω φ j ω, yû j,o ω, x, z + j=nω+ φ j ω, y v j,o ω, x, z..5 he decomposition is in the L, D orthonormal basis of the eigenfunctions φ j ω, y of the self-adjoint differential operator in y, y + ω c φ j ω, y = λ j ωφ j ω, y, y φ j ω, D = y φ j ω, =, j =,,..., with eigenvalues λ j ω that are simple 6. o simplify the analysis, we assume in this paper that the wave speed is homogeneous cy = c o..6 he results for variable wave speeds are similar in all the essential aspects. he simplification brought by.6 amounts to having explicit expressions of the eigenfunctions, which are independent of the frequency φ j y = D cos π j y..7 D he eigenvalues are π kd λ j ω = D π j,.8 where k = ω/c o is the wavenumber, and only the first Nω of them are non-negative kd Nω = π +..9 he notation stands for the integer part. We suppose for simplicity that Nω remains constant in the bandwidth ω B/, ω + B/, and write from now on Nω = N. he propagating components in.5 satisfy the two-dimensional Helmholtz equation x + z + β j ω û j,o ω, x, z = F j ω, xδz, j =,..., N,. with outgoing, radiation conditions at x + z. he evanescent components solve x + z β j ω v j,o ω, x, z = F j ω, xδz, j > N,. with decay condition v j,o ω, x, z at x + z. Here we introduced the coefficients of the source profile in the basis of the eigenfunctions and the mode wavenumbers β j ω = F j ω, x = D λ j ω = π D kd π dy φ j y fω, x, y, j,. j, j..3 We assume that none of the β j ω vanishes in the bandwidth, so that there are no standing waves. hat is to say, kd π = N + αω, αω, for all ω ω B/, ω + B/..4 3
4 .. he paraxial regime. We now introduce the paraxial scaling for the ideal waveguide, with the source emitting a beam that propagates along the z-axis. As we show below, this happens when the crossrange profile of the source is larger than the wavelength. he source generates a quasi-plane wave, with slowly varying envelope satisfying a Schrödinger-like equation. Explicitly, we assume that the source is of the form f ε t, x, y = ft, εx, y.5 where ε is a small dimensionless parameter defined as the ratio of the central wavelength λ and the transverse width r of the source. Standard diffraction theory gives that the Rayleigh length for a beam with initial width r = λ /ε is of the order of r /λ = λ /ε. he Rayleigh length is the distance along the z axis from the beam waist to the place where the beam area is doubled by diffraction. herefore, we look at the wavefield at Oε cross-range scales, similar to r, and at Oε range scale, similar to the Rayleigh length. We rename the field in this scaling as p ε ot, X, y, Z = p o t, X ε, y, Z ε..6 he Fourier coefficients of.6 are given by the scaled version of.5 p ε oω, X, y, Z = Nω φ j yû ε j,oω, X, Z + j=nω+ φ j y v ε j,oω, X, Z,.7 with propagating mode amplitudes û ε j,o F j ω, εx = X. hey can be written as satisfying the scaled equation., with the source replaced by û ε j,oω, X, Z = ε in terms of the outgoing Green s function Ĝ o βj ω, x, z = i 4 H dx Fj ω, X Ĝo β j ω, X X, Z, ε ε β j ω x + z. Here H is the Hankel function of the first kind, and because ε, we can use its asymptotic form for a scaled range Z > i X X 4 H β j ω ε + Z ε 4 4 ε i X X πβ j ω ε i πβ j ωz exp he propagating components of the wave field become + Z ε 4 { iβ j ω û ε j,oω, X, Z a j,o ω, X, Z exp / exp iβ j ω Z ε + X X }. Z iβ j ω Zε, X X ε + Z ε 4 with a j,o ω, X, Z = i πβ j ωz dx iβj ωx X exp F j ω, X,.8 Z 4
5 for j =,..., N. he evanescent components are obtained similarly from. v j,oω, ε X, Z e j,o ω, X, Z exp β j ω Zε, and e j,o ω, X, Z = πβ j ωz dx exp β jωx X F j ω, X,.9 Z for j N +. hese modes are exponentially damped and can be neglected. In summary, the paraxial approximation of the wave field is given by p ε oω, X, y, Z φ j ya j,o ω, X, Ze iβjω Z ε. It is a superposition of forward going modes with complex valued amplitudes a j,o given by.8, and solving the paraxial equations iβj ω Z + X aj,o ω, X, Z =, j =,..., N,. with initial conditions a j,o ω, X, Z = = a j,ini ω, X := iβ j ω F j ω, X, j =,..., N.. 3. Wave propagation in random waveguides. In this section we consider a waveguide with fluctuating boundary and analyze the wave field under the following scaling assumptions:. he transverse width r of the source and the central wavelength λ satisfy r = ε λ, 3. as in the previous section. he correlation length l ε of the boundary fluctuations is similar to r, l ε = ε l r, 3. so that there is a non-trivial interaction between the boundary fluctuations and the wavefield. Here l is the scaled order-one correlation length defined below.. he scale L ε of the propagation distance is much larger than λ. More precisely, L ε /λ = Oε. 3.3 Recall that the Rayleigh length for a beam with initial width r and central wavelength λ is of the order of r /λ ε λ in absence of random fluctuations. he high-frequency scaling assumption 3.3 ensures that the propagation distance is similar to the Rayleigh length. 3. he amplitude of the boundary fluctuations is small, of the order of ε 3/ λ. As we will show, this scaling is precisely the one that gives a cumulative scattering effect of order one after the propagation distance L ε. We use the hyperbolicity of the problem to truncate mathematically the boundary fluctuations to the range interval, L/ε. he bound L/ε is the maximum range of the fluctuations that can affect the waves up to the observation time ε of order ε. he lower bound in the range interval coincides with the location of the source. It is motivated by two facts: First, we observe the waves at positive ranges. Second, the backscattered field is negligible in the scaling regime defined above, as we show later in section 4.3. he boundary fluctuations are modeled with a random process µ ε x, z = D + ε 3/ µ εx, εz, z, L/ε
6 he process µ is bounded, zero-mean, stationary and mixing, meaning in particular that its covariance is integrable. Because our method of solution flattens the boundary by changing coordinates, we require that µ is twice differentiable, with almost surely bounded derivatives. Its covariance function is given by and we denote by R o ξ its integral over ζ, R ξ, ζ = E µξ + ξ, ζ + ζµξ, ζ, 3.5 R o ξ = dζ Rξ, ζ. 3.6 Our assumption on the differentiability of µ implies that R o is four times differentiable. Note that ξ = is the maximum of the integrated covariance R o ξ, so we have R o =. 3.7 We define the scaled square amplitude σ and correlation length l of the boundary fluctuations through the equations R o = σ l, R o R o = l Change of coordinates. We introduce the change of coordinates from x, y, z to x, η, z, with η = yd ε x, z. 3.9 It straightens the boundary y = ε x, z to η = D, for any x R and z, L/ε. he pressure field in the new coordinates is denoted by P ω, x, η, z = p ω, x, η ε x, z, z. 3. D It satisfies the simple boundary conditions P ω, x, D, z = η P ω, x,, z =, 3. and the partial differential equation D x + z + ε + ε η ε η η ε ε η + η ε ε η ε ε η + k P = f ε ω, x, ηδz, 3. derived from. and 3. using the chain rule. Here and are the gradient and Laplacian operators in x, z and f ε is the source of the form.5. When substituting the model 3.4 in 3., we obtain that P satisfies a randomly perturbed problem x + z + ε 3/ µεx, εz η + k +... P ω, x, η, z = fω, εx, ηδz. 3.3 he higher-order terms denoted by the dots are ε 5/ + Oε 3/ η µ η P + 3ε 3 µ + Oε η P ε 7/ µ + Oε 3/ η η P. hey come from the expansions in ε of the coefficients in 3., and are negligible in the limit ε considered in section 4. More precisely, µ is a ϕ-mixing process, with ϕ L / R +, as stated in 3,
7 3.. Wave decomposition. Equation 3.3 is not separable, but we can still write its solution in the L, D basis of the eigenfunctions.7. he expansion is similar to.5 P ω, x, η, z = φ j ηû j ω, x, z + j>n We define the forward and backward going wave mode amplitudes a j and b j by φ j η v j ω, x, z. 3.4 a j ω, x, z = ûjω, x, z + iβ j ω zû j ω, x, z e iβjωz, b j ω, x, z = ûjω, x, z iβ j ω zû j ω, x, z e iβjωz, 3.5 so that the complex valued amplitudes of the propagating modes can be written as Definition 3.5 implies that û j ω, x, z = a j ω, x, ze iβjωz + b j ω, x, ze iβjωz. z a j ω, x, ze iβjωz + z b j ω, x, ze iβjωz =, j =,..., N. 3.6 his equation is needed to specify uniquely the propagating mode amplitudes, because they each satisfy a single boundary condition in the range, L/ε of the fluctuations. o derive these boundary conditions, let us observe that a j and b j must be constant in z, and in z L/ε,, because the boundary is flat outside, L/ε. Moreover, the radiation conditions lim a jω, x, z =, z lim b jω, x, z =, z imply that the mode amplitudes satisfy a j ω, x, z = =, 3.7 b j ω, x, z = L/ε =. 3.8 he last equation is the boundary condition for b j. he boundary value a j ω, x, z = + follows from the jump conditions across the plane z = of the source in equation 3.3. We have with F j defined by.. his gives û j + =, zû j + = F j ω, εx, â j + b j + =, iβ jâ j b j + = F j ω, εx, and therefore a j ω, x, + = iβ j ω F j ω, εx. 3.9 Substituting 3.4 in 3.3, and using the orthonormality of the eigenfunctions φ j, we find that the wave mode amplitudes solve paraxial equations coupled by the random fluctuations in z, L/ε, N iβj z + x aj + e iβjz xb j ε 3/ µεx, εze iβjz q jl al e iβlz + b l e iβ lz + q jl v l, 3. l= l>n N iβj z + x bj + e iβjz xa j ε 3/ µεx, εze iβjz q jl al e iβlz + b l e iβ lz + q jl v l. 3. l>n 7 l=
8 We dropped the higher-order terms that do not play a role in the limit ε, and replaced the equality with the approximate sign. o simplify our notation, we omit henceforth all the arguments in the equations, except those of µ. he arguments will be spelled out only in definitions. he coupling matrix in is given by q jl = D π dη φ j ηφ l η = j δjl. 3. D It takes this simple diagonal form because we assumed a homogeneous background speed c o. If we had a variable speed cy, the matrix {q jl } would not be diagonal, and the modes with j l would be coupled. However, the results of the asymptotic analysis below would still hold, because the coupling would become negligible in the limit ε considered in section 4, due to rapid phases arising in the right hand sides of 3., 3.. he equations for the evanescent components are obtained similarly, z + x βj vj ε 3/ µεx, εz q jj v j, 3.3 and they are augmented with the decay conditions v j ω, x, z as x + z, for all j N he limit process. We characterize next the wave field in the asymptotic limit ε. We begin with the paraxial long range scaling that gives significant net scattering, and then take the limit. he scaling has already been described at the beginning of section Asymptotic scaling. We obtain from that the propagating mode amplitudes satisfy the block diagonal system of partial differential equations iβj z + x e iβjz x aj e iβjz x iβ j z + x ε 3/ e iβ jz aj q b jj µεx, εz j e iβjz, 4. b j for j =,..., N. Again, the approximate sign means equal to leading order. Because the right hand side in 4. is small, of order ε 3/, and has zero statistical expectation, it follows from 8, Chapter 6 that there is no net scattering effect until we reach ranges of order ε. hus, we let z = Z/ε, 4. with scaled range Z independent of ε. he source directivity in the range direction suggests observing the wavefield on a cross-range scale that is smaller than that in range. We choose it as x = X/ε, 4.3 with scaled cross-range X independent of ε, to balance the two terms in the paraxial operators in 4.. Our goal is to characterize the ε limit of the mode amplitudes in the paraxial long range scaling regime We denote them by a ε jω, X, Z = a j ω, X ε, Z ε and b ε jω, X, Z = b j ω, X ε, Z ε, 4.4 and obtain from that they satisfy the scaled system iβ j Z + X e iβjz/ε X a ε j e iβjz/ε X iβ j Z + X b ε X, j ε µ Z e q iβjz/ε / jj ε e iβjz/ε a ε j b ε j, 4.5 for j =,..., N, with initial conditions and end conditions a ε jω, X, = a j,ini ω, X := iβ j ω F j ω, X, 4.6 b ε jω, X, L =
9 4.. he random propagator. Let us rewrite 4.5 in terms of the random propagator matrix P ε ω, X, X, Z C N N, the solution of the initial value problem Z P ε ω, X, X, Z = X, ε µ Z H ω, X, Z / ε ε + G ω, X, Z ε P ε ω, X, X, Z, Z >, P ε ω, X, X, = δx X I. 4.8 Here I is the N N identity matrix, δx is the Dirac delta distribution in X, and G and H are matrices with entries given by partial differential operators in X, with deterministic coefficients. We can define them from 4.5 once we note that the solution a ε ω, X, Z = a ε ω, X, Z a ε N. ω, X, Z, b ε ω, X, Z = b ε ω, X, Z b ε N. ω, X, Z 4.9 follows from a ε ω, X, Z b ε ω, X, Z = dx P ε ω, X, X a, Z ε ω, X, b ε ω, X,. 4. Here b ε ω, X, is the vector of backward going amplitudes at the beginning of the randomly perturbed section of the waveguide, and it can be eliminated using the boundary identity a ε ω, X, L = dx P ε ω, X, X a, L ε ω, X, b ε ω, X. 4., he initial conditions a ε ω, X, are given in 4.6. We obtain from 4.5 that H and G have the block form H a H b H = H b H a G a G b, G = G b G a where the bar denotes complex conjugation. he blocks are diagonal, with entries, 4. H a jl = i δ jl q jj, H b jl = i δ jl q jj β j β j e iβjz/ε 4.3 and G a jl = i δ jl β j X, G b jl = i δ jl e iβjz/ε β X, j, l =,..., N. 4.4 j he entries of the diagonal blocks depend only on the mode indices and the frequency, via β j ω. he entries of the off-diagonal blocks are rapidly oscillating, due to the large phases proportional to Z/ε. he symmetry relations satisfied by the blocks in H and G imply that the propagator has the form P ε ω, X, X, Z = with N N complex, diagonal blocks ε and R ε. ε ω, X, X, Z R ε ω, X, X, Z R ε ω, X, X, Z ε ω, X, X, Z, he diffusion limit. he limit of P ε as ε is a multi-dimensional Markov diffusion process, with entries satisfying a system of Itô-Schrödinger equations. his follows from the diffusion approximation theorem 4, 5, see also 8, Chapter 6, applied to system 4.8. When computing the generator of the limit process, we obtain that due to the fast phases in the offdiagonal blocks of H and G, the forward and backward going amplitudes decouple as ε. his implies 9
10 that there is no backscattered field in the limit, because the backward going amplitudes b ε are set to zero at Z = L. Equation 4. simplifies as a ε ω, X, Z = dx ε ω, X, X, Za ε ω, X,, 4.6 where the initial conditions a ε ω, X, are given in 4.6. We call the complex diagonal matrix ε ω, X, X, Z = diag ε ω, X, X, Z,..., ε Nω, X, X, Z 4.7 the transfer process, because it gives the amplitudes of the forward going modes at positive ranges Z, in terms of the initial conditions at Z =. he limit transfer process is described in the next proposition. It follows straight from 4, 5. Proposition 4.. As ε, ε ω, X, X, Z converges weakly and in distribution to the diffusion Markov process ω, X, X, Z. his process is complex and diagonal matrix valued, with entries j ω, X, X, Z solving the Itô-Schrödinger equations d j ω, X, X i, Z = β j ω X q jj R o 8βj ω j ω, X, X, ZdZ + i q jj β j ω jω, X, X, ZdBX, Z, 4.8 for Z >, and initial conditions j ω, X, X, = δx X, j =,..., N. 4.9 Equations 4.8 are uncoupled, but they are driven by the same Brownian field BX, Z, satisfying E BX, Z =, E BX, ZBX, Z = min{z, Z }R o X X, 4. with R o defined in 3.6. hus, the transfer coefficients j are statistically correlated. he weak convergence in distribution means that we can calculate the limit ε of statistical moments of ε, smoothed by integration over X against the initial conditions, using the Markov diffusion defined by In applications we have a fixed ε, and we use Proposition 4. to approximate the statistical moments of the amplitudes of the forward going waveguide modes. When comparing the Itô-Schrödinger equations 4.8 to the deterministic Schrödinger equations. satisfied by the amplitudes in the ideal waveguides, we see that the random boundary scattering effect amounts to a net diffusion, as described by the last two terms in 4.8. We show next how this leads to loss of coherence of the waves, that is to exponential decay in range of the mean field. We also study the propagation of energy of the modes and quantify the decorrelation properties of the random fluctuations of their amplitudes. 5. Statistics of the wave field. We begin in section 5. with the analysis of the coherent field. Explicitly, we estimate the mean forward going mode amplitudes in the paraxial long range regime. raditional imaging methods rely on these being large with respect to their random fluctuations. However, this is not the case, because E a ε j ω, X, Z decay exponentially with Z, at rates that increase monotonically with mode indices j. he second moments of the amplitudes do not decay, but there is decorrelation over the modes and the frequency and cross-range offsets, as shown in sections 5.3 and 5.4. Understanding these decorrelations is key to designing time reversal and imaging methods that are robust at low SNR. Robust means that wave focusing in time reversal or imaging is essentially independent of the realization of the random boundary fluctuations, it is statistically stable. Low SNR means that the coherent mean field, the signal, is faint with respect to its random fluctuations, the noise. 5.. he coherent field. he mean modal amplitudes are E a ε jω, X, Z dx E j ω, X, X, Z a j,ini ω, X, 5.
11 with mean transfer matrix satisfying the partial differential equation Z E j ω, X, X i, Z = β j ω X E j ω, X, X, Z, S j ω Z >, 5. with mode-dependent damping coefficients S j ω = 8β j ω qjj R o = D σ π l N + αω / j / j / 4, 5.3 with units of length. Here we used definitions.3,.4 and 3.8, and obtained equation 5.3 by taking expectations in 4.8. Its solution is given by E j ω, X, X βj ω, Z = πiz exp Z S j ω + iβ jωx X, 5.4 Z and the mean modal amplitudes are obtained from equations 5. and 4.6 E a ε jω, X, Z i dx Fj ω, X exp Z πβ j ωz S j ω + iβ jωx X Z = a j,o ω, X, Z exp Z, 5.5 S j ω with a j,o the solution of the paraxial wave equation.-. in the ideal waveguide. he mean wave field follows from 3.4, after neglecting the evanescent part, E P ω, X ε, η, Z ε It is different than the field in the ideal waveguides p o ω, X ε, η, Z ε φ j ηa j,o ω, X, Z exp Z S j ω + iβ jω Z ε. 5.6 φ j ηa j,o ω, X, Z exp iβ j ω Zε, 5.7 because of the exponential decay of the mean mode amplitudes, on range scales S j ω. 5.. High-frequency and low-snr regime. We call the length scales S j ω the mode-dependent scattering mean free paths, because they give the range over which the modes become essentially incoherent, with low SNR, E a ε j SNR j,ω = ω, X, Z E a εj ω, X, Z E a εj ω, X, Z exp Z, if Z S j ω. 5.8 S j ω he second moments E a ε j ω, X, Z are calculated in the next section, and they do not decay with range. his is why equation 5.8 holds. he scattering mean free paths decrease monotonically with mode indices j, as shown in 5.3. he first mode encounters less often the random boundary, and has the longest scattering mean free path S ω = 3D N + αω / σ π /4 3D N l σ π l. 5.9 he highest indexed mode scatters most frequently at the boundary, and its scattering mean free path S N ω = D αω N + αω σ π l N / 4 αω 8 S ω N 5 5.
12 is much smaller than S ω, when N is large. o be complete, we also have and S j ω S ω s4 s 4, if j = sn, s,, N 4 S j ω S ω, if j = on. j 4 Our analysis of time reversal and imaging is carried in a high-frequency regime, with waveguide depth D much larger than the central wavelength λ o or, equivalently, with N. We also assume a low-snr regime, with scaled range Z exceeding the scattering mean free path of all the modes, so that none of the amplitudes a j are coherent. his is the most challenging case for sensor array imaging, because the wave field measured at the sensors is essentially just noise. We model the low-snr regime using the dimensionless large parameter and observe from 5.3 that γ = Z, 5. S ω Z γ, for all j =,..., N. 5. S j ω 5.3. he second moments. he quantification of SNR and the analysis of time reversal and imaging involves the second moments of the mode amplitudes. Recall that a ε jω, X, Z dx ε j ω, X, X, Z a j,ini ω, X 5.3 with j ε the entries of the diagonal transfer matrix ε. o calculate the second moments, we need to estimate E j ε l ε. he equations for ε j ω, X, X, Zl εω, X, X, Z follow from the forward scattering approximation of 4.8, for Z >, with initial condition Z j ε l ε i i β j ω X i qjj µx, Z/ε β j ω ε / β l ω X ε j ε l q ll µx, Z/ε β l ω ε j ε l, 5.4 ε j ω, X, X, ε l ω, X, X, = δx X δx X. 5.5 heir statistical distribution is characterized in the limit ε by the diffusion approximation theorem 4, 5, see also 8, Chapter 6. It is the distribution of j ω, X, Z l ω, X, Z, with j the limit transfer coefficients in Proposition 4.. his gives the approximate relation E a ε jω, X, Za ε l ω, X, Z dx dx a j,ini ω, X a l,ini ω, X E j ω, X, X, Z l ω, X, X, Z. 5.6 he calculation of E j l is given in appendix A. We summarize the results in Propositions
13 5.3.. he single mode and frequency moments. It is easier to calculate the diagonal moments, with j = l, and the same frequency ω = ω = ω. We have the following result proved in appendix A. Proposition 5.. For all j =,..., N, and all the frequencies ω ω πb, ω + πb, E j ω, X, X, Z j ω, X, X, Z = β { jω πz exp iβj ωx X X X Z with kernel C o defined by Z S j ω ds C o X X s + X X s }, 5.7 C o X = R ox R o. 5.8 he general second moment formula does not have an explicit form in arbitrary regimes. But it can be approximated in the low-snr regime 5.. he expression 5.7 also simplifies in that regime, as stated in the following proposition, which we prove below. Proposition 5.. In the low-snr regime 5., and under the assumption X = X = X, the right hand side in 5.7 is essentially zero, unless X X 3 S j ω, 5.9 l γ S ω and the moment formula simplifies to E j ω, X, X, Z j ω, X, X, Z β j πz exp iβ j X X X X X X Z Xd,j ω, 5. with 3 Sj ω 3 S j ω X d,j ω = l = l l. 5. Z γ S ω If the initial points X and X are different, but still close enough to satisfy X X l, 5. the moment formula becomes E j ω, X, X, Z j ω, X, X, Z β j πz exp iβj X X X X Z exp X X + X X + X X X X Xd,j ω. 5.3 Proof. We see from definitions 3.6 and 5.8 that C o X for X l. herefore, ds C o X X s if X X l, and the right hand side in 5.7 becomes negligible, of order O e Z/Sj. In the case X X l we obtain similarly that the damping term is of order Z, and the right hand side in 5.7 is exponentially 3
14 small. It is only when X X l that the moment does not vanish. hen, we can approximate the kernel C o in the integral with its first nonzero term in the aylor expansion around zero, using the relations that follow from We have C o =, C o =, and C o = R o R o = l, 5.4 ds C o X X s X X 6l, if X X l, and the right hand side in 5.7 is of the order exp X X Z 3l S j. his gives the condition 5.9, and the simpler moment formula 5. follows. Essentially the same proof applies in the case X X, because we can still expand the integrand in 5.7 by assumption he two mode and frequency moments. he general second moment formula is derived in appendix A, in the low-snr regime 5.. It has a complicated expression that we do not repeat here, but it simplifies for nearby frequencies, as stated below. Proposition 5.3. he modes decorrelate under the low-snr assumption 5. E j ω, X, X, Z l ω, X, X, Z if j l, 5.5 for any two frequencies ω, ω and cross-ranges X, X. he modes also decorrelate for frequency offsets that exceed Ω d,j ω = S jωβ j ωl Z β j ω = β jω β j ω S j ωβ j ωl γ S ω, 5.6 where β j ω is the derivative of β jω with respect to ω. For much smaller frequency offsets satisfying ω ω Ω d,j ω, ω = ω + ω, 5.7 the moment formula is { E j ω, X, X, Z j ω, X, X, Z β jω i πz exp βj ω X X β j ω X X Z X X + X X + X X X X X d,j ω } Decorrelation properties. We already stated the decorrelation of the modes in Proposition 5.3. But even for a single mode, we have decorrelation over cross-range and frequency offsets. he decoherence length of mode j is denoted by X d,j ω, and it is defined in 5.. It is the length scale over which the second moment at frequency ω decays with cross-range. It follows from 5. that X d,j is much smaller than the correlation length, for all the modes, and that it decreases monotonically with j. he first mode has the largest decorrelation length 3 X d, ω = l γ. 5.9 because it scatters less often at the boundary. he decoherence length of the highest mode is much smaller in high-frequency regimes with N, X d,n ω = X d, ω S N ω S ω l 3αω N 5/. 8 γ 5.3 4
15 he decorrelation frequency is derived in appendix A.. It is given by 5.6 or, more explicitly, by Ω d,j ω ωσ π 3 64γ l λ 3 N + αω j 5/ N 9 j / 4, 5.3 it is much smaller than ω for all the modes, and it decreases monotonically with j, starting from Ω d, ω ωσ π 3 l/λ 3 4γ N he forward model. Let us gather the results and summarize them in the following model of the pressure field P ω, X ε, η, Z ε φ j η Z iβ j ω eiβjω ε dx j ω, X, X, Z D dη φ j η fω, X, η, 6. where the symbol stands for approximate, in distribution. hat is to say, the statistical moments of the random pressure field P are approximately equal to those of the right handside. he first and second moments follow from Propositions In our analysis of time reversal and imaging we take small frequency offsets, satisfying ω Ω d,j ω, so that we can use the simpler moment formula 5.8. he computation of the fourth moments of the transfer coefficients is quite involved. We estimate in appendix B some of them, for a particular combination of the mode indices and arguments. hese moments are used in the next sections to show the statistical stability of the time reversal and coherent interferometric imaging functions. We analyze next time reversal and imaging in the low SNR regime, and assume for convenience that the source.5 has the separable form ft, X, η = ϕt X X ρ, η η, 6. θ X θ η θ X θ η meaning that the same pulse ϕt is emitted from all the points in the support of the non-negative source density ρ. We scale this support with the dimensionless parameters θ X and θ η, and normalize the source by dx θ X dη X X ρ, η η =. 6.3 θ η θ X θ η he coefficients F j ω, X = ϕω θ X θ η D X X dη φ j ηρ, η η, 6.4 θ X θ η are proportional to the Fourier coefficients ϕω of the pulse, and are thus supported in the frequency interval ω o B/, ω o + B/. he bandwidth B is small enough so that we can freeze the number of propagating modes to that at the central frequency, as explained in section.. he width of the pulse ϕt is inverse proportional to B, and we distinguish two regimes: he broadband regime with B ε ω o, and the narrowband regime with B ε ω o. he comparison with ε is because the source is at range Z A /ε from the array, and the modes arrive at time intervals of order /ε. Broadband pulses have smaller support than these travel times, meaning that we can observe the different arrivals of the modes, at least in the ideal waveguides. o analyze the resolution of time reversal and imaging, we study in detail the case of a source density localized around the point X, η,. We say that we study the point spread time reversal and imaging functions, because the source has small support. Note however that it is not a point source. Its support is quantified by the positive parameters θ X and θ η which are small, but independent of ε. 5
16 7. ime reversal. Let us denote by Dt, X, η the pressure field measured in a time window ψt/ ε at an array A, with aperture modeled by the indicator function A X, η = AX X Aη η, 7. at range z A = Z A /ε. Here X is the scaled cross-range in the array, related to the cross-range x by x = X/ε, and A X R and A η, D are intervals in X and η. he window ψ is a function of dimensionless arguments, of support of order one, and ε denotes the length of time of the measurements. Because the waves travel distances of order ε, we scale ε as ε = /ε, with of order one. In time reversal, the array takes the recorded field Dt, X, η, time reverses it and emits D ε t, X, η back in the medium. We study in this section the resolution of the refocusing of the waves at the source, in the high-frequency and low-snr regime described in section 5.. Because we have a random waveguide, the resolution analysis includes that of statistical stability, given in section Mathematical model of time reversal. We have in our notation t Dt, X, η = A X, ηψ ε P t, x = X ε, η, z A = Z A ε, 7. with mathematical model following from 6., t Dt, X, η A X, ηψ he time reversed field has Fourier transform with ε N φ j η Dω, X, η A X, η Z A dω eiβjω ε iωt iβ j ω dx Fj ω, X j ω, X, X, Z A. 7.3 D R t, X, η = D ε t, X, η 7.4 D R ω, X, η = e iω ε Dω, X, η 7.5 du φ j η π u/ Z A eiβjω ε ε ψu iβ j ω ε u/ dx Fj ω ε u, X j ω ε u, X, X, Z A. 7.6 he small frequency shifts ε u/ are due to the time scaling, and we can neglect them in the source terms F j and in the amplitude factor /β j. he model of the observed wave field at search locations x s, η s, z s = X s ε, ηs, Zs ε is given by Ot, X s, η s, Z s = dx dω exp φ j η s π iβ j ω Z A Z s iωt ε iβ j ω dη φ j η D R ω, X, η j ω, X, X s, Z A Z s, 7.7 using reciprocity. Note the similarity with equation 6., except that the source is now at the array, which we approximate in 7.7 as a continuum, instead of a discrete collection of sensors. his approximation is convenient for the analysis, because sums over the sensors are replaced by integrals over the X and η apertures, of lengths A X and A η. Using 7.5 in 7.7 and letting Γ jl = D dη Aη η φ j ηφ l η, 7.8 6
17 we obtain dx Ot, X s, η s, Z s dη θ X θ η du π ψu We define the time reversal function by X X ρ, η η N θ X θ η j,l= dω Γ ε jl π eiω t ϕω φ jη s φ l η 4β j ωβ l ω dx AX X j ω, X, X s, Z A Z s l ω ε u exp iβ j ω Z A Z s ε iβ l ω ε u, X, X, Z A ZA ε. 7.9 J R X s, η s = Ot = ε, X s, η s, Z s =. 7. It models the wave field observed at the time instant t = ε, at the source range Z A. his is when and where the refocusing occurs. In the case of a source density that is tightly supported around X, η, we may approximate J R by J R X s, η s with frequency-dependent kernel point spread function M R ω, X s, η s j,l= Here we used the source normalization 6.3. dω π ϕω MR ω, X s, η s, 7. φ j η s φ l η du Γ jl 4β j ωβ l ω π ψu j ω, X, X s, Z A l ω ε u exp iβ j ω Z A ε iβ l dx AX X, X, X, Z A ω ε u ZA ε Resolution analysis. If the time reversal process is statistically stable, then we can estimate its refocusing resolution by studying the mean of 7.. We refer to the next section for the analysis of the statistical stability of J R. he mean time reversal function follows from 5.8 and E J R X s, η s = dω π ϕω E MR ω, X s, η s, 7.3 with E M R ω, X s, η s A X 8πZ A φ j η s φ j η β j ωz A Γ jj ψ β j ω exp iβ jω Z A dx A X X A X exp Xs X Xd,j ω X Xs + X X s X. 7.4 Moreover, letting A X = A X, A X, 7.5 7
18 we obtain after integrating in X that Note that E M R ω, X s, η s A X 8πZ A sinc φ j η s φ j η β j ωz A Γ jj ψ β j j ω βj ω A X X s X exp Z A exp Xs X Xd,j ω { iβj ω Z A X s X }. 7.6 τ j = β jωz A 7.7 are the scaled travel times of the modes, so only those modes that arrive within the support of the window ψ contribute in Cross-range resolution. We observe in 7.6 that modes contribute differently to the focusing in cross-range X, with resolution { } X s X πz A X,j ω := min X d,j ω,. 7.8 β j ω A X Recall from 5. and 5.9 that X d,j decreases monotonically with j X d,j ω X d,ω N + αω / j / / 3 4j /, X d, ω = l N γ, 7.9 whereas πz A β j ω A X Z AD N + αω / j / /, 7. A X increases with j. hus, in the high-frequency regime with N, the cross-range resolution for the highorder modes is determined by the decorrelation length, even for large apertures. he cross-range resolution of the first modes may be determined by the aperture, but only if it is large enough, A X Z AD l γ N It may appear at this point that the time reversal process can give good results even for small apertures A X. However, we will see in section 7.3 that large apertures are needed for statistical stability. he modes with higher indices give the best cross-range resolution, but they travel at smaller speed. hus, the focusing improves when we increase the recording time, because the array can capture the late arrivals of the high-order modes see Figure Depth resolution. o study the focusing in η, we evaluate the point spread function at crossrange X s = X. We have N M R ω, X, η s, X, η φ j η s φ j η Γ jj, 7. β j ω where N is the number of modes with arrival times in the recording window, he coefficients Γ jj are given by Γ jj = D dη Aη ηφ jη = A η D + η D sinc τ j <, for j =,,..., N N. 7.3 πj / η D 8 η D sinc πj / η, 7.4 D
19 depth profile a.u..5 N =5 N = N =9 cross range profile a.u N =5 N = N = η.5.5 x Fig. 7.. Depth profile left and cross range profile right of the mean point spread function for the time reversal functional. Here Z A =, l =, σ =.5, k = 6, D = so that N = 9. he array diameter A X is supposed to be smaller than the critical value 7. which is about. N is the cut-off number modes smaller than N are recorded and reemitted. Note that the high modes play an important role. he larger N is, the better the resolution. for an array in the set A η = η, η, D. hey satisfy Γ jj = in the full aperture case A η =, D. he sum in 7. is maximum at η s = η, because all the terms are positive. he point spread function is smaller at other depths, because of cancellations in the sum of the oscillatory terms. We can make this more explicit in the high-frequency regime, with N, if we write D πn k, φ jη D cos j and interpret 7. as a Riemann sum, which we then approximate with an integral. Consider for simplicity the full aperture case, where E M R ω, X, η s A X 8πZ A N A X 8π Z A N Dβ j cos N A X 8π Z A Λ N /N j kη s cos j N cos j / N kη η / j / N kη s η α, Λ α x = kη, 7.5 N kη N ds cossx s. 7.6 he function Λ α becomes proportional to the Bessel function of first kind J as α, more explicitly, we have Λ x = π/j x so that E M R ω, X, η s A X 6πZ A J kη s η, if N N. 7.7 We can then estimate the depth resolution as the distance between the peak of J, that occurs when η s = η, and its first zero, that occurs when k η s η.4 first zero of J. herefore, the depth resolution of time reversal with full aperture is equal to the diffraction limit η s η η ω =.4 k, 7.8 if the array records the waves long enough to capture almost all the propagating modes. he resolution deteriorates if N is much smaller than N. Indeed for small α we have Λ α x α/π sincαx and therefore the depth resolution is η ω πn/kn see Figure 7.. 9
20 7.3. Statistical stability. We now show that the time reversal function is statistically stable, meaning that the refocusing of the wave at the original source location does not depend on the the realization of the random medium but only on its statistical distribution, and the point spread function is approximately equal to its expectation. We restrict the analysis of statistical stability to the case of full aperture, where the calculations are simpler because the coupling matrix Γ jl becomes the identity. he point spread function follows from 7. M R ω, X s, η s = φ j η s φ j η ψ β j ω β j ωz A dx A X X j ω, X, X s, Z A j ω, X, X A X, Z A, and its variance at the source location is Var M R ω, X, η = φ j η φ J η β J ω ψ β j ωz A ψ β J ωz A dxdy A X X AX Y A X β j ω j,j= { E j ω, X, X, Z A J ω, Y, X, Z A E j ω, X, X, Z A E J ω, Y, X, Z A }. From appendix B first case we find that the variance is much smaller than the square expectation when A X l, and therefore the point spread function is equal to its mean approximately. he results contained in appendix B second case also show that if A X < l, the variance of the point spread function is large, and therefore that time reversal refocusing may be unstable in this case. here is however another mechanism that can ensure statistical stability of the focal spot if the array is small. Indeed, if the bandwidth of ϕ is larger than the decorrelation frequency, then the variance dω dω Var J R X, η = π π ϕω, X ϕω Cov M R ω, X, η, M R ω, X, η is small because the covariance of the point spread function at two frequencies becomes approximately zero if the frequency gap is large enough. herefore, if the pulse has large bandwidth, then the time-reversal focal spot is statistically stable even for small arrays. 8. Imaging. he sharp and stable focusing of the time reversal process in the random waveguide is due the backpropagation of the time reversed field D R in exactly the same waveguide. ime reversal is a physical experiment, where the waves can be observed in the vicinity of the source, as they refocus. In imaging we only have access to the data measured at the array, and the backpropagation to the search points is synthetic. Because we cannot know the fluctuations of the boundary, we simply ignore them in the synthetic backpropagation and obtain the so-called reversed time migration imaging function. We analyze it in section 8. and show that it does not give useful results in the low-snr regime. In particular, we show that the images are not statistically stable with respect to realizations of the fluctuations. Stability can be achieved by imaging with local cross-correlations of the array measurements. Local means that we recall the decorrelation properties of the random mode amplitudes described in section 5.4, and cross-correlate the measurements over receivers located at nearby cross-ranges X, and projected on the same eigenfunctions. he resulting coherent interferometric imaging method is analyzed in section Reverse time migration. he reverse time migration function is given by the time reversed data D R propagated migrated in the ideal waveguide to the search points x s, η s, z s = X s ε, ηs, Zs ε. Its mathematical expression follows from., with amplitudes.8 replaced by a j,o ω, X, Z dx j,o ω, X, X, Z iβ j ω D dη φ j η D R ω, X, η. 8. he ideal transfer coefficients j,o ω, X, X βj ω, Z = πiz exp iβj ωx X Z 8.
21 are defined by the Green s functions of the paraxial operator in.. We obtain dω exp iβ j ω Z A Z s J M X s, η s, Z s = φ j η s ε iωt dx dη φ j η π iβ j ω D R ω, X, η j,o ω, X, X s, Z A Z s, 8.3 t= ε with the right hand side evaluated at the same time t = ε as in time reversal. We assume again a tightly supported source density normalized by 6.3 and substitute the model 7.5 R of D in 8.3 to obtain dω J M X s, η s, Z s = π ϕω MM ω, X s, η s, Z s, 8.4 with frequency-dependent kernel point spread function M M ω, X s, η s, Z s j,l= φ j η s φ l η du Γ jl 4β j ωβ l ω π ψu dx AX X j,o ω, X, X s, Z A Z s l ω ε u, X, X, Z A exp iβ j ω Z A Z s ε iβ l ω ε u ZA ε he mean imaging function. Let us take for simplicity the case of full aperture in depth, where the coupling matrix Γ jl given by 7.8 becomes the identity. We obtain from 8.5 and the moment formula 5.4 that E M M ω, X s, η s, Z s = A X 8πZ A φ j η s φ j η ψ β j ω dx A X X A X β j ωz A exp iβ jω Z A exp Z A S j ω iβ jω Zs ε X Xs + X Moreover, assuming the aperture A X defined in 7.5, and integrating in X, we get E M M ω, X s, η s, Z s = A X 8πZ A φ j η s φ j η β j ωz A ψ β j ω βj ω A X sinc X s X exp Z A exp Z A X s X. 8.6 S j ω iβ jω Zs ε }. 8.7 { iβj ω Z A X s X his result is almost the same as in the ideal waveguide, except for the damping coefficients exp Z A /S j. he sinc kernel in the mean point spread function gives the focusing in cross-range, with mode-dependent resolution X S X X,j ω = πz A β j ω A X. 8.8 he best resolution is for the first mode, that has the largest wavenumber β ω πn/d k, and gives the Rayleigh cross-range resolution X, ω πz A k A X. 8.9
22 he focusing of the point spread function M M in range can only be due to the summation of the rapidly oscillating terms exp iβ j Z s /ε. But these terms are weighted by exp Z A /S j, which decay fast in j. he first term dominates in E M M ω, X, η, Z s = A X 8πZ A φ j η β j ω ψ β j ωz A exp Z A S j ω iβ jω Zs ε, 8. so the mode diversity does not lead to focusing in range, as is the case in ideal waveguides. Nevertheless, the mean reverse time migration function peaks at Z s = because of the integral over the bandwidth in dω E J M X, η, Z s = π ϕωe MM ω, X, η, Z s, 8. and the range resolution is of the order ε /β ω o B. When we evaluate the point spread function at Z s = and X s = X, we obtain E M M ω, X s = X, η s, Z s = = A X 8πZ A φ j η s φ j η ψ β j ω β j ωz A his is a sum of the oscillatory functions φ j η s φ j η = { cos π j η s η + cos π j η s + η } D D D e Z A/S jω. 8. multiplied by positive weights, which are small and decay fast in j. he first term dominates in 8. and there is no depth resolution at all. We show next that these small weights also indicate the lack of statistical stability of the reverse time migration function Stability analysis. o assess the stability of the reverse time migration, we calculate its variance at the source location Var J M X, η, = E J M X, η, E J M X, η,. We have from the results above that E J M X, η, A X 8πZ A dω π ϕω N φ j η β j ω ψ β j ωz A he second moment of J M is E J M X, η, dω dω π π ϕω φ j ϕω η φ l η 6β j,l= j ω βl ω iβ j ω u β l exp ω u Z A dx l,o ω, X, X, Z A E j ω ε u, X, X, Z A l e Z A/S jω. 8.3 du π du π ψu ψu dx AX X AX X j,o ω, X, X, Z A ω ε u, X, X, Z A, 8.4 and we recall from Proposition 5.3 that only the diagonal terms j = l contribute to the expectation. We also assume a small bandwidth B Ω d,j, for all the modes j, so that we can use the simpler moment formula 5.8. We obtain E J M X, η, ϕ φ 4 j η β 8πZ A β j ω o ψ j ω o Z A dx dx AX X AX X exp X X Xd,j ω, 8.5 o
23 after approximating the modal wavenumbers by their value at the central frequency. his expression can be approximated further, after integrating in X and X, and supposing that the decoherence lengths X d,j are much smaller than the array aperture, E J M X, η, ϕ A X π 8πZ A X d,j ω o φ 4 j η β βj ω o ψ j ω o Z A. 8.6 he second moment 8.6 is clearly much larger than the square of the mean 8.3, which is exponentially small in range. Although the mean of the imaging function is focused at the source, it cannot be observed because it is dominated by its random fluctuations. he reverse time migration lacks statistical stability with respect to the realizations of the random fluctuations of the boundary of the waveguide. he calculations above are for a small bandwidth, satisfying B Ω d,j for all the modes captured in the recording window. he calculations are more complicated for a larger bandwidth, but the conclusion remains that reverse time migration is not stable with respect to different realization of the random boundary fluctuations. 8.. Coherent interferometric imaging. he main idea of the coherent interferometric CIN imaging approach is to backpropagate synthetically to the imaging points the local cross-correlations of the array measurements, instead of the measurements themselves. By local we mean that because of the statistical decorrelation properties of the random mode amplitudes described in section 5.4, we cross-correlate the data Dω, X, η at nearby frequencies and cross-ranges X, after projecting it on the subspace of one eigenfunction φ j at a time. he projection gives the coefficients D j ω, X = D dη φ j η Dω, X, η, 8.7 which are directly proportional to the coefficients F j of the source only in the case of an array spanning the entire depth of the waveguide. We assume this case here, because it simplifies the analysis of the focusing and stability of the CIN function. We also take a small source, meaning that we essentially compute the CIN point spread function. he model of the coefficients 8.7 is D j ω, X AX X ϕωφ jη e iβjωz A/ε du iβ j ω π iβ ψue j ωuza/ j ω ε u, X, X, Z A, 8.8 and we cross-correlate them at cross-ranges satisfying X X X d,j ω, and at frequency offsets ω ω Ω Ω d,j. 8.9 We take such small Ω to simplify the second moment formulas. he CIN image is formed by backpropagating the cross-correlations to the imaging point, using the Green s function in the ideal waveguide. We first define the CIN image in the X, Z-domain: with J CIN X s, Z s = J CIN j X s, Z s 8. Jj CIN X s, Z s dω dω = π π Ωω ω e iβjω βjω Z s Z A ε dx dx Xd,j X X D j ω, X D j ω, X j,o ω, X, X s, Z A Z s j,o ω, X, X s, Z A Z s, 8. where Xd,j are indicator functions of the cross-range interval X d,j ω, X d,j ω calculated at the central frequency ω = ω + ω /. Similarly, Ω is the indicator function of the frequency interval Ω, Ω. 3
WAVE PROPAGATION AND IMAGING IN RANDOM WAVEGUIDES. Liliana Borcea
WAVE PROPAGATION AND IMAGING IN RANDOM WAVEGUIDES by Liliana Borcea Contents 1. Introduction..................................................... 2 2. Set-up...........................................................
More informationPulse propagation in random waveguides with turning points and application to imaging
Pulse propagation in random waveguides with turning points and application to imaging Liliana Borcea Mathematics, University of Michigan, Ann Arbor and Josselin Garnier Centre de Mathématiques Appliquées,
More informationSOURCE ESTIMATION WITH INCOHERENT WAVES IN RANDOM WAVEGUIDES
SOURCE ESTIMATION WITH INCOHERENT WAVES IN RANDOM WAVEGUIDES SEBASTIAN ACOSTA, RICARDO ALONSO AND LILIANA BORCEA Abstract. We study an inverse source problem for the acoustic wave equation in a random
More informationarxiv: v1 [math-ph] 15 Jul 2016
IMAGING IN RANDOM MEDIA WITH CONVEX OPTIMIZATION LILIANA BORCEA AND ILKER KOCYIGIT arxiv:67.465v math-ph 5 Jul 6 Abstract. We study an inverse problem for the wave equation where localized wave sources
More informationGhost Imaging. Josselin Garnier (Université Paris Diderot)
Grenoble December, 014 Ghost Imaging Josselin Garnier Université Paris Diderot http:/www.josselin-garnier.org General topic: correlation-based imaging with noise sources. Particular application: Ghost
More informationInterferometric array imaging in clutter
INSTITUTE OF PHYSICS PUBLISHING Inverse Problems 1 005 1419 1460 INVERSE PROBLEMS doi:10.1088/066-5611/1/4/015 Interferometric array imaging in clutter Liliana Borcea 1, George Papanicolaou and Chrysoula
More informationA Mode Coupling Theory for Random Waveguides with Turning Points
A Mode Coupling Theory for Random Waveguides with Turning Points by Derek MacLean Wood A dissertation submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy Mathematics
More informationSéminaire Laurent Schwartz
Séminaire Laurent Schwartz EDP et applications Année 2013-2014 Josselin Garnier Multiscale analysis of wave propagation in random media. Application to correlation-based imaging Séminaire Laurent Schwartz
More information2D Imaging in Random Media with Active Sensor Array
2D Imaging in Random Media with Active Sensor Array Paper Presentation: Imaging and Time Reversal in Random Media Xianyi Zeng rhodusz@stanford.edu icme, Stanford University Math 221 Presentation May 18
More informationPASSIVE SENSOR IMAGING USING CROSS CORRELATIONS OF NOISY SIGNALS IN A SCATTERING MEDIUM
PASSIVE SENSOR IMAGING USING CROSS CORRELATIONS OF NOISY SIGNALS IN A SCATTERING MEDIUM JOSSELIN GARNIER AND GEORGE PAPANICOLAOU Abstract. It is well known that the travel time or even the full Green s
More informationImaging with Ambient Noise II
Imaging with Ambient Noise II George Papanicolaou Stanford University Michigan State University, Department of Mathematics Richard E. Phillips Lecture Series April 21, 2009 With J. Garnier, University
More informationSchrödinger equation with random potential
Lille January, 213 Schrödinger equation with random potential Josselin Garnier Université Paris VII http://www.proba.jussieu.fr/~garnier Lille January, 213 Physical derivation of the paraxial wave equation
More informationRESOLUTION ENHANCEMENT FROM SCATTERING IN PASSIVE SENSOR IMAGING WITH CROSS CORRELATIONS. Josselin Garnier. George Papanicolaou
Inverse Problems and Imaging Volume 8, No. 3, 214, 645 683 doi:1.3934/ipi.214.8.645 RESOLUTION ENHANCEMENT FROM SCATTERING IN PASSIVE SENSOR IMAGING WITH CROSS CORRELATIONS Josselin Garnier Laboratoire
More informationAdaptive interferometric imaging in clutter and optimal illumination
Adaptive interferometric imaging in clutter and optimal illumination Liliana Borcea and George Papanicolaou and Chrysoula Tsogka Computational and Applied Mathematics, MS 134, Rice University, 61 Main
More information1 Coherent-Mode Representation of Optical Fields and Sources
1 Coherent-Mode Representation of Optical Fields and Sources 1.1 Introduction In the 1980s, E. Wolf proposed a new theory of partial coherence formulated in the space-frequency domain. 1,2 The fundamental
More informationUncertainty Principle Applied to Focused Fields and the Angular Spectrum Representation
Uncertainty Principle Applied to Focused Fields and the Angular Spectrum Representation Manuel Guizar, Chris Todd Abstract There are several forms by which the transverse spot size and angular spread of
More informationImaging and time reversal in random media
Imaging and time reversal in random media Liliana Borcea George Papanicolaou Chrysoula Tsogka James Berryman June 4, 2002 Abstract We present a general method for estimating the location of small, well-separated
More informationACMAC s PrePrint Repository
ACMAC s rerint Repository Filtering deterministic layer effects in imaging Liliana Borcea and Fernando Gonzalez del Cueto and George apanicolaou and Chrysoula Tsogka Original Citation: Borcea, Liliana
More informationc 2013 International Press
COMMUN. MATH. SCI. Vol., No., pp. 33 67 c 3 International Press WAVE PROPAGATION IN WAVEGUIDES WITH RANDOM BOUNDARIES RICARDO ALONSO, LILIANA BORCEA, AND JOSSELIN GARNIER Abstract. We give a detailed analysis
More informationDerivation of the General Propagation Equation
Derivation of the General Propagation Equation Phys 477/577: Ultrafast and Nonlinear Optics, F. Ö. Ilday, Bilkent University February 25, 26 1 1 Derivation of the Wave Equation from Maxwell s Equations
More informationEquations of linear stellar oscillations
Chapter 4 Equations of linear stellar oscillations In the present chapter the equations governing small oscillations around a spherical equilibrium state are derived. The general equations were presented
More informationChapter 4. Adjoint-state methods
Chapter 4 Adjoint-state methods As explained in section (3.4), the adjoint F of the linearized forward (modeling) operator F plays an important role in the formula of the functional gradient δj of the
More informationImaging and time reversal in random media
0 1 2 Imaging and time reversal in random media Liliana Borcea George Papanicolaou Chrysoula Tsogka James Berryman December 15, 2002 Abstract We present a general method for estimating the location of
More information1.1 A Scattering Experiment
1 Transfer Matrix In this chapter we introduce and discuss a mathematical method for the analysis of the wave propagation in one-dimensional systems. The method uses the transfer matrix and is commonly
More informationDecoherence and the Classical Limit
Chapter 26 Decoherence and the Classical Limit 26.1 Introduction Classical mechanics deals with objects which have a precise location and move in a deterministic way as a function of time. By contrast,
More informationWave equation techniques for attenuating multiple reflections
Wave equation techniques for attenuating multiple reflections Fons ten Kroode a.tenkroode@shell.com Shell Research, Rijswijk, The Netherlands Wave equation techniques for attenuating multiple reflections
More informationA universal filter for enhanced imaging with small arrays
A universal filter for enhanced imaging with small arrays Liliana Borcea, Thomas Callaghan 2, Josselin Garnier 3, and George Papanicolaou 4 Computational and Applied Mathematics, MS34, Rice University,
More informationResolution Analysis of Passive Synthetic Aperture Imaging of Fast Moving Objects
SIAM J. IMAGING SCIENCES Vol. 10, No. 2, pp. 665 710 c 2017 Society for Industrial and Applied Mathematics Resolution Analysis of Passive Synthetic Aperture Imaging of Fast Moving Objects L. Borcea, J.
More information2. As we shall see, we choose to write in terms of σ x because ( X ) 2 = σ 2 x.
Section 5.1 Simple One-Dimensional Problems: The Free Particle Page 9 The Free Particle Gaussian Wave Packets The Gaussian wave packet initial state is one of the few states for which both the { x } and
More informationLecture 2 Supplementary Notes: Derivation of the Phase Equation
Lecture 2 Supplementary Notes: Derivation of the Phase Equation Michael Cross, 25 Derivation from Amplitude Equation Near threshold the phase reduces to the phase of the complex amplitude, and the phase
More informationTRANSPORT OF POWER IN RANDOM WAVEGUIDES WITH TURNING POINTS
COMMUN. MATH. SCI. Vol. 15, No. 8, pp. 37 371 c 17 International Press TRANSPORT OF POWER IN RANDOM WAVEGUIDES WITH TURNING POINTS LILIANA BORCEA, JOSSELIN GARNIER, AND DEREK WOOD Abstract. We present
More informationThe Generation of Ultrashort Laser Pulses II
The Generation of Ultrashort Laser Pulses II The phase condition Trains of pulses the Shah function Laser modes and mode locking 1 There are 3 conditions for steady-state laser operation. Amplitude condition
More informationTheory of Ship Waves (Wave-Body Interaction Theory) Quiz No. 2, April 25, 2018
Quiz No. 2, April 25, 2018 (1) viscous effects (2) shear stress (3) normal pressure (4) pursue (5) bear in mind (6) be denoted by (7) variation (8) adjacent surfaces (9) be composed of (10) integrand (11)
More informationkg meter ii) Note the dimensions of ρ τ are kg 2 velocity 2 meter = 1 sec 2 We will interpret this velocity in upcoming slides.
II. Generalizing the 1-dimensional wave equation First generalize the notation. i) "q" has meant transverse deflection of the string. Replace q Ψ, where Ψ may indicate other properties of the medium that
More informationMathematical Introduction
Chapter 1 Mathematical Introduction HW #1: 164, 165, 166, 181, 182, 183, 1811, 1812, 114 11 Linear Vector Spaces: Basics 111 Field A collection F of elements a,b etc (also called numbers or scalars) with
More informationPage 404. Lecture 22: Simple Harmonic Oscillator: Energy Basis Date Given: 2008/11/19 Date Revised: 2008/11/19
Page 404 Lecture : Simple Harmonic Oscillator: Energy Basis Date Given: 008/11/19 Date Revised: 008/11/19 Coordinate Basis Section 6. The One-Dimensional Simple Harmonic Oscillator: Coordinate Basis Page
More informationOPTIMAL PERTURBATION OF UNCERTAIN SYSTEMS
Stochastics and Dynamics c World Scientific Publishing Company OPTIMAL PERTURBATION OF UNCERTAIN SYSTEMS BRIAN F. FARRELL Division of Engineering and Applied Sciences, Harvard University Pierce Hall, 29
More informationNondiffracting Waves in 2D and 3D
Nondiffracting Waves in 2D and 3D A thesis submitted in partial fulfillment of the requirements for the degree of Bachelor of Science in Physics from the College of William and Mary by Matthew Stephen
More information1. Reflection and Refraction of Spherical Waves
1. Reflection and Refraction of Spherical Waves Our previous book [1.1] was completely focused on the problem of plane and quasi-plane waves in layered media. In the theory of acoustic wave propagation,
More informationReview of Linear Time-Invariant Network Analysis
D1 APPENDIX D Review of Linear Time-Invariant Network Analysis Consider a network with input x(t) and output y(t) as shown in Figure D-1. If an input x 1 (t) produces an output y 1 (t), and an input x
More informationRESOLUTION ANALYSIS OF PASSIVE SYNTHETIC APERTURE IMAGING OF FAST MOVING OBJECTS
RESOLUTION ANALYSIS OF PASSIVE SYNTHETIC APERTURE IMAGING OF FAST MOVING OBJECTS L. BORCEA, J. GARNIER, G. PAPANICOLAOU, K. SOLNA, AND C. TSOGKA Abstract. In this paper we introduce and analyze a passive
More information6. LIGHT SCATTERING 6.1 The first Born approximation
6. LIGHT SCATTERING 6.1 The first Born approximation In many situations, light interacts with inhomogeneous systems, in which case the generic light-matter interaction process is referred to as scattering
More informationTRANSPORT IN POROUS MEDIA
1 TRANSPORT IN POROUS MEDIA G. ALLAIRE CMAP, Ecole Polytechnique 1. Introduction 2. Main result in an unbounded domain 3. Asymptotic expansions with drift 4. Two-scale convergence with drift 5. The case
More informationStochastic Volatility and Correction to the Heat Equation
Stochastic Volatility and Correction to the Heat Equation Jean-Pierre Fouque, George Papanicolaou and Ronnie Sircar Abstract. From a probabilist s point of view the Twentieth Century has been a century
More informationLecture Notes 1: Vector spaces
Optimization-based data analysis Fall 2017 Lecture Notes 1: Vector spaces In this chapter we review certain basic concepts of linear algebra, highlighting their application to signal processing. 1 Vector
More informationOPTIMAL PERTURBATION OF UNCERTAIN SYSTEMS
Stochastics and Dynamics, Vol. 2, No. 3 (22 395 42 c World Scientific Publishing Company OPTIMAL PERTURBATION OF UNCERTAIN SYSTEMS Stoch. Dyn. 22.2:395-42. Downloaded from www.worldscientific.com by HARVARD
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 informationLecture 12: Detailed balance and Eigenfunction methods
Lecture 12: Detailed balance and Eigenfunction methods Readings Recommended: Pavliotis [2014] 4.5-4.7 (eigenfunction methods and reversibility), 4.2-4.4 (explicit examples of eigenfunction methods) Gardiner
More informationENGI 9420 Lecture Notes 8 - PDEs Page 8.01
ENGI 940 ecture Notes 8 - PDEs Page 8.0 8. Partial Differential Equations Partial differential equations (PDEs) are equations involving functions of more than one variable and their partial derivatives
More information4 Classical Coherence Theory
This chapter is based largely on Wolf, Introduction to the theory of coherence and polarization of light [? ]. Until now, we have not been concerned with the nature of the light field itself. Instead,
More informationStochastic Spectral Approaches to Bayesian Inference
Stochastic Spectral Approaches to Bayesian Inference Prof. Nathan L. Gibson Department of Mathematics Applied Mathematics and Computation Seminar March 4, 2011 Prof. Gibson (OSU) Spectral Approaches to
More informationWe finally come to the determination of the CMB anisotropy power spectrum. This set of lectures will be divided into five parts:
Primary CMB anisotropies We finally come to the determination of the CMB anisotropy power spectrum. This set of lectures will be divided into five parts: CMB power spectrum formalism. Radiative transfer:
More informationSynthetic Aperture Radar Imaging and Motion Estimation via Robust Principal Component Analysis
SIAM J. IMAGING SCIENCES Vol. 6, No. 3, pp. 5 76 c 3 Society for Industrial and Applied Mathematics Synthetic Aperture Radar Imaging and Motion Estimation via Robust Principal Component Analysis Liliana
More informationPhysics 506 Winter 2004
Physics 506 Winter 004 G. Raithel January 6, 004 Disclaimer: The purpose of these notes is to provide you with a general list of topics that were covered in class. The notes are not a substitute for reading
More informationABOUT TRAPPED MODES IN OPEN WAVEGUIDES
ABOUT TRAPPED MODES IN OPEN WAVEGUIDES Christophe Hazard POEMS (Propagation d Ondes: Etude Mathématique et Simulation) CNRS / ENSTA / INRIA, Paris MATHmONDES Reading, July 2012 Introduction CONTEXT: time-harmonic
More informationLast Update: April 7, 201 0
M ath E S W inter Last Update: April 7, Introduction to Partial Differential Equations Disclaimer: his lecture note tries to provide an alternative approach to the material in Sections.. 5 in the textbook.
More informationFILTERING RANDOM LAYERING EFFECTS IN IMAGING
FILTERING RANDOM LAYERING EFFECTS IN IMAGING L. BORCEA, F. GONZÁLEZ DEL CUETO, G. PAPANICOLAOU, AND C. TSOGKA Abstract. Objects that are buried deep in heterogeneous media produce faint echoes which are
More informationStochastic Structural Dynamics Prof. Dr. C. S. Manohar Department of Civil Engineering Indian Institute of Science, Bangalore
Stochastic Structural Dynamics Prof. Dr. C. S. Manohar Department of Civil Engineering Indian Institute of Science, Bangalore Lecture No. # 33 Probabilistic methods in earthquake engineering-2 So, we have
More informationAddition of Angular Momenta
Addition of Angular Momenta What we have so far considered to be an exact solution for the many electron problem, should really be called exact non-relativistic solution. A relativistic treatment is needed
More informationReview of Fundamental Equations Supplementary notes on Section 1.2 and 1.3
Review of Fundamental Equations Supplementary notes on Section. and.3 Introduction of the velocity potential: irrotational motion: ω = u = identity in the vector analysis: ϕ u = ϕ Basic conservation principles:
More informationCoherent imaging without phases
Coherent imaging without phases Miguel Moscoso Joint work with Alexei Novikov Chrysoula Tsogka and George Papanicolaou Waves and Imaging in Random Media, September 2017 Outline 1 The phase retrieval problem
More information. (70.1) r r. / r. Substituting, we have the following equation for f:
7 Spherical waves Let us consider a sound wave in which the distribution of densit velocit etc, depends only on the distance from some point, ie, is spherically symmetrical Such a wave is called a spherical
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 informationA system that is both linear and time-invariant is called linear time-invariant (LTI).
The Cooper Union Department of Electrical Engineering ECE111 Signal Processing & Systems Analysis Lecture Notes: Time, Frequency & Transform Domains February 28, 2012 Signals & Systems Signals are mapped
More informationSeparation of Variables in Linear PDE: One-Dimensional Problems
Separation of Variables in Linear PDE: One-Dimensional Problems Now we apply the theory of Hilbert spaces to linear differential equations with partial derivatives (PDE). We start with a particular example,
More information4. Complex Oscillations
4. Complex Oscillations The most common use of complex numbers in physics is for analyzing oscillations and waves. We will illustrate this with a simple but crucially important model, the damped harmonic
More informationAnswers to Problem Set Number MIT (Fall 2005).
Answers to Problem Set Number 6. 18.305 MIT (Fall 2005). D. Margetis and R. Rosales (MIT, Math. Dept., Cambridge, MA 02139). December 12, 2005. Course TA: Nikos Savva, MIT, Dept. of Mathematics, Cambridge,
More informationGraceful exit from inflation for minimally coupled Bianchi A scalar field models
Graceful exit from inflation for minimally coupled Bianchi A scalar field models Florian Beyer Reference: F.B. and Leon Escobar (2013), CQG, 30(19), p.195020. University of Otago, Dunedin, New Zealand
More informationJ10M.1 - Rod on a Rail (M93M.2)
Part I - Mechanics J10M.1 - Rod on a Rail (M93M.2) J10M.1 - Rod on a Rail (M93M.2) s α l θ g z x A uniform rod of length l and mass m moves in the x-z plane. One end of the rod is suspended from a straight
More informationAn inverse scattering problem in random media
An inverse scattering problem in random media Pedro Caro Joint work with: Tapio Helin & Matti Lassas Computational and Analytic Problems in Spectral Theory June 8, 2016 Outline Introduction and motivation
More informationDamped Harmonic Oscillator
Damped Harmonic Oscillator Wednesday, 23 October 213 A simple harmonic oscillator subject to linear damping may oscillate with exponential decay, or it may decay biexponentially without oscillating, or
More informationCorrelation based imaging
Correlation based imaging George Papanicolaou Stanford University International Conference on Applied Mathematics Heraklion, Crete September 17, 2013 G. Papanicolaou, ACMAC-Crete Correlation based imaging
More informationUltrasound experiment by M. Fink
Ultrasound experiment by M. Fink cf. A. Tourin, M. Fink, and A. Derode, Multiple scattering of sound, Waves Random Media (), R3-R6. Experimental set-up for a time-reversal experiment through a multiple-scattering
More informationIntegrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows
Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows Alexander Chesnokov Lavrentyev Institute of Hydrodynamics Novosibirsk, Russia chesnokov@hydro.nsc.ru July 14,
More informationFiltering Deterministic Layer Effects in Imaging
SIAM REVIEW Vol. 54, No. 4, pp. 757 798 c 01 Society for Industrial and Applied Mathematics Filtering Deterministic Layer Effects in Imaging L. Borcea F. González del Cueto G. Papanicolaou C. Tsogka Abstract.
More informationMTH Linear Algebra. Study Guide. Dr. Tony Yee Department of Mathematics and Information Technology The Hong Kong Institute of Education
MTH 3 Linear Algebra Study Guide Dr. Tony Yee Department of Mathematics and Information Technology The Hong Kong Institute of Education June 3, ii Contents Table of Contents iii Matrix Algebra. Real Life
More informationarxiv: v1 [physics.comp-ph] 28 Jan 2008
A new method for studying the vibration of non-homogeneous membranes arxiv:0801.4369v1 [physics.comp-ph] 28 Jan 2008 Abstract Paolo Amore Facultad de Ciencias, Universidad de Colima, Bernal Díaz del Castillo
More informationIdentification Methods for Structural Systems. Prof. Dr. Eleni Chatzi Lecture March, 2016
Prof. Dr. Eleni Chatzi Lecture 4-09. March, 2016 Fundamentals Overview Multiple DOF Systems State-space Formulation Eigenvalue Analysis The Mode Superposition Method The effect of Damping on Structural
More informationElectromagnetic fields and waves
Electromagnetic fields and waves Maxwell s rainbow Outline Maxwell s equations Plane waves Pulses and group velocity Polarization of light Transmission and reflection at an interface Macroscopic Maxwell
More informationLecture 5 Notes, Electromagnetic Theory II Dr. Christopher S. Baird, faculty.uml.edu/cbaird University of Massachusetts Lowell
Lecture 5 Notes, Electromagnetic Theory II Dr. Christopher S. Baird, faculty.uml.edu/cbaird University of Massachusetts Lowell 1. Waveguides Continued - In the previous lecture we made the assumption that
More informationAppendix C: Recapitulation of Numerical schemes
Appendix C: Recapitulation of Numerical schemes August 31, 2009) SUMMARY: Certain numerical schemes of general use are regrouped here in order to facilitate implementations of simple models C1 The tridiagonal
More informationStochastic nonlinear Schrödinger equations and modulation of solitary waves
Stochastic nonlinear Schrödinger equations and modulation of solitary waves A. de Bouard CMAP, Ecole Polytechnique, France joint work with R. Fukuizumi (Sendai, Japan) Deterministic and stochastic front
More informationMath Ordinary Differential Equations
Math 411 - Ordinary Differential Equations Review Notes - 1 1 - Basic Theory A first order ordinary differential equation has the form x = f(t, x) (11) Here x = dx/dt Given an initial data x(t 0 ) = x
More informationWeek 9 Generators, duality, change of measure
Week 9 Generators, duality, change of measure Jonathan Goodman November 18, 013 1 Generators This section describes a common abstract way to describe many of the differential equations related to Markov
More informationThis is a Gaussian probability centered around m = 0 (the most probable and mean position is the origin) and the mean square displacement m 2 = n,or
Physics 7b: Statistical Mechanics Brownian Motion Brownian motion is the motion of a particle due to the buffeting by the molecules in a gas or liquid. The particle must be small enough that the effects
More informationStep index planar waveguide
N. Dubreuil S. Lebrun Exam without document Pocket calculator permitted Duration of the exam: 2 hours The exam takes the form of a multiple choice test. Annexes are given at the end of the text. **********************************************************************************
More informationLecture 4 Fiber Optical Communication Lecture 4, Slide 1
ecture 4 Dispersion in single-mode fibers Material dispersion Waveguide dispersion imitations from dispersion Propagation equations Gaussian pulse broadening Bit-rate limitations Fiber losses Fiber Optical
More informationLinear Algebra. Min Yan
Linear Algebra Min Yan January 2, 2018 2 Contents 1 Vector Space 7 1.1 Definition................................. 7 1.1.1 Axioms of Vector Space..................... 7 1.1.2 Consequence of Axiom......................
More informationLecture 12: Detailed balance and Eigenfunction methods
Miranda Holmes-Cerfon Applied Stochastic Analysis, Spring 2015 Lecture 12: Detailed balance and Eigenfunction methods Readings Recommended: Pavliotis [2014] 4.5-4.7 (eigenfunction methods and reversibility),
More informationFiber Optics. Equivalently θ < θ max = cos 1 (n 0 /n 1 ). This is geometrical optics. Needs λ a. Two kinds of fibers:
Waves can be guided not only by conductors, but by dielectrics. Fiber optics cable of silica has nr varying with radius. Simplest: core radius a with n = n 1, surrounded radius b with n = n 0 < n 1. Total
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 informationEXISTENCE OF GUIDED MODES ON PERIODIC SLABS
SUBMITTED FOR: PROCEEDINGS OF THE FIFTH INTERNATIONAL CONFERENCE ON DYNAMICAL SYSTEMS AND DIFFERENTIAL EQUATIONS June 16 19, 2004, Pomona, CA, USA pp. 1 8 EXISTENCE OF GUIDED MODES ON PERIODIC SLABS Stephen
More informationSpectrum and Exact Controllability of a Hybrid System of Elasticity.
Spectrum and Exact Controllability of a Hybrid System of Elasticity. D. Mercier, January 16, 28 Abstract We consider the exact controllability of a hybrid system consisting of an elastic beam, clamped
More information13.42 READING 6: SPECTRUM OF A RANDOM PROCESS 1. STATIONARY AND ERGODIC RANDOM PROCESSES
13.42 READING 6: SPECTRUM OF A RANDOM PROCESS SPRING 24 c A. H. TECHET & M.S. TRIANTAFYLLOU 1. STATIONARY AND ERGODIC RANDOM PROCESSES Given the random process y(ζ, t) we assume that the expected value
More informationKirchhoff, Fresnel, Fraunhofer, Born approximation and more
Kirchhoff, Fresnel, Fraunhofer, Born approximation and more Oberseminar, May 2008 Maxwell equations Or: X-ray wave fields X-rays are electromagnetic waves with wave length from 10 nm to 1 pm, i.e., 10
More informationBackground velocity estimation with cross correlations of incoherent waves in the parabolic scaling
IOP PUBLISHING 34pp INVERSE PROBLEMS doi:.88/66-56/5/4/455 Background velocity estimation with cross correlations of incoherent waves in the parabolic scaling Josselin Garnier and Knut Sølna Laboratoire
More informationLow Emittance Machines
Advanced Accelerator Physics Course RHUL, Egham, UK September 2017 Low Emittance Machines Part 1: Beam Dynamics with Synchrotron Radiation Andy Wolski The Cockcroft Institute, and the University of Liverpool,
More informationSummary of free theory: one particle state: vacuum state is annihilated by all a s: then, one particle state has normalization:
The LSZ reduction formula based on S-5 In order to describe scattering experiments we need to construct appropriate initial and final states and calculate scattering amplitude. Summary of free theory:
More informationPhysics 202 Laboratory 5. Linear Algebra 1. Laboratory 5. Physics 202 Laboratory
Physics 202 Laboratory 5 Linear Algebra Laboratory 5 Physics 202 Laboratory We close our whirlwind tour of numerical methods by advertising some elements of (numerical) linear algebra. There are three
More information5. LIGHT MICROSCOPY Abbe s theory of imaging
5. LIGHT MICROSCOPY. We use Fourier optics to describe coherent image formation, imaging obtained by illuminating the specimen with spatially coherent light. We define resolution, contrast, and phase-sensitive
More information