Startseite Stability for the electromagnetic inverse source problem in inhomogeneous media
Artikel
Lizenziert
Nicht lizenziert Erfordert eine Authentifizierung

Stability for the electromagnetic inverse source problem in inhomogeneous media

  • Yue Zhao EMAIL logo
Veröffentlicht/Copyright: 31. Mai 2022

Abstract

This paper is concerned with the stability of the inverse source problem for Maxwell’s equations in an inhomogeneous background medium. We show that the stability estimate consists of the Lipschitz-type data discrepancy and the high frequency tail of the source function, where the latter decreases as the upper bound of the frequency increases. The analysis employs scattering theory to obtain the holomorphic domain and an upper bound for the resolvent of the elliptic operator.

MSC 2010: 35R30; 78A46

Award Identifier / Grant number: 12001222

Funding statement: The research of the author is supported in part by NSFC (No. 12001222).

A Proof of Proposition 2.3

Define

𝐄 ( t ) = sin t 𝒫 𝒫 .

By the spectral theorem, we have the following mapping properties for 𝐄 ( t ) :

𝐄 ( t ) : 𝒟 α 𝒟 α + 1 2 ,
𝐄 ( t ) : C ( t ; ( 𝒟 α , 𝒟 α + 1 2 ) ) .

It can be verified that 𝐄 ( t ) solves the equations

𝐄 ( t ) = 0 , 𝐄 ( 0 ) = 0 , t 𝐄 ( 0 ) = 𝐈 ,

where := t 2 + 𝒫 . Set 0 = t 2 + × ( × ) . We use the notation [ A , B ] = A B - B A for the commutator.

Choose a large enough such that supp χ B a , and take χ a C c ( B a ) , χ a 1 , on supp χ . Assume that ψ a C ( t × 3 ) satisfies

supp ψ a ( t × 3 B R ) { ( t , x ) : | x | + T a t | x | + T a + 1 }

and

ψ a | t × B R + ϵ = ψ a 0 ( t ) ,

where ψ a 0 C c ( T a + R , T a + R + 1 ) is a function of t only.

First, we show that

(A.1) ψ a 𝐄 ( t ) χ a C ( t ; ( , 𝒟 ) ) .

In fact, let χ 0 , χ 1 C 0 ( B a ) such that χ 0 1 on supp χ a and χ 1 1 on supp χ 0 . The non-trapping condition in Assumption 2.2 implies

(A.2) χ 0 𝐄 ( t ) χ a | t > T a C ( ( T a , ) ; ( , 𝒟 ) ) .

Let 0 := t t - × ( × ) . We have

0 ( ( 1 - χ 0 ) 𝐄 ( t ) χ a ) = - [ 0 , χ 0 ] ( χ 1 𝐄 ( t ) χ a ) = : 𝐆 a ( t ) ,
𝐆 a ( t ) | t > T a C ( ( T a , ) ; ( , H comp 1 ( B a B R ¯ ) 3 ) ) .

It is clear that ( 1 - χ 0 ) 𝐄 ( t ) χ a satisfies the initial conditions

( 1 - χ 0 ) 𝐄 ( 0 ) χ a = 0 , ( 1 - χ 0 ) t 𝐄 ( 0 ) χ a = 0 .

By the Duhamel principle and the support property of the propagator 𝐄 0 of free-space electromagnetic wave equations in odd dimensions

supp 𝐄 0 ( t , ) = { ( x , y ) : | x - y | = | t | } ,

we have

(A.3) ψ a ( 1 - χ 0 ) 𝐄 ( t ) χ a = ψ a T a t 𝐄 0 ( t - s ) 𝐆 a ( s ) d s C ( ( T a , ) ; ( , 𝒟 ) ) .

Since ψ 0 on { t < T a + R } , relations (A.2) and (A.3) imply (A.1).

Take a function ζ a C such that

ζ a ( x , t ) | | x | > R = { 1 , t | x | + T a , 0 , t | x | + T a + 1 ,
ζ a ( x , t ) | | x | < R + ϵ = ζ a 0 ( t ) = { 1 , t R + T a + ϵ , 0 , t R + T a + 1 .

It follows from the above claim that

𝐅 a ( t ) = [ , ζ a ] ( 𝐄 ( t ) χ a ) C ( ; ( , 𝒟 1 2 ) ) .

Notice that ζ a ( x , t ) 1 for t 0 . Thus, [ , ζ a ] = 0 on { t 0 } , and consequently 𝐅 a ( t ) 0 for t 0 .

Choose χ b C c ( B a ) such that χ b = 1 near B R and χ a 1 on supp χ b . Let 𝐖 a ( t ) solve

0 𝐖 a ( t ) = - ( 1 - χ b ) 𝐅 a ( t ) , 𝐖 a ( t ) 0 , t 0 .

It is clear that 𝐖 a ( t ) g C ( t ; 𝒟 1 / 2 ) .

In what follows, we show that for g ,

(A.4) supp ( 𝐖 a ( t ) g ) ( x ) { ( t , x ) : | x | + T a t | x | + C a } .

Here C a = 2 a + T a + 1 . We write

- 𝐖 a ( t ) = ( 1 - χ b ) ζ a H ( t ) 𝐄 ( t ) χ a - H ( t ) 𝐄 0 ( t ) ( 1 - χ b ) χ a + 𝐐 a ( t ) ,

where H ( t ) is the Heaviside function. It is easy to see that the first term has a support in

{ ( t , x ) : 0 t | x | + C a } .

The last term 𝐐 a ( t ) solves

0 𝐐 a ( t ) = 0 𝐖 a ( t ) - ( 1 - χ b ) 0 ( ζ a H ( t ) 𝐄 ( t ) χ a ) + [ 0 , χ b ] ζ a H ( t ) 𝐄 ( t ) χ a - 0 ( H ( t ) 𝐄 0 ( t ) ( 1 - χ b ) χ a )
= - ( 1 - χ b ) 𝐅 a ( t ) + ( 1 - χ b ) ( χ a + 𝐅 a ( t ) ) + [ 0 , χ b ] ζ a H ( t ) 𝐄 ( t ) χ a ( 1 - χ b ) χ a - χ a ( 1 - χ b ) χ a
= - [ × ( × ) , χ b ] ζ a H ( t ) 𝐄 ( t ) χ a .

Observe that

[ × ( × ) , χ b ] ζ a H ( t ) 𝐄 ( t ) χ a L comp 2 ( B a B R ) 3

vanishes on { t a + T a + 1 } . Therefore, it follows from the sharp Huygen principle that

supp ( 𝐐 a ( t ) g ) ( x ) { ( t , x ) : 0 t | x | + C a } .

Notice also that ( 1 - χ b ) 𝐅 a ( t ) g = 0 on { t | x | + T a } . By the finite speed of propagation, we have

supp ( 𝐖 a ( t ) g ) ( x ) { ( t , x ) : | x | + T a t } .

Choose χ c C c ( B a ) such that χ c = 1 near B R and χ b = 1 on supp χ c . Let

R a # ( λ ) = t λ ( ζ a H ( t ) 𝐄 ( t ) χ a + ( 1 - χ c ) 𝐖 a ( t ) ) .

A straightforward calculation yields

( 𝒫 - λ 2 ) R a # ( λ ) = t λ ( ( ζ a H ( t ) 𝐄 ( t ) χ a ) + ( ( 1 - χ c ) 𝐖 a ( t ) ) )
= t λ ( [ , ζ a ] ( H ( t ) 𝐄 ( t ) χ a ) + ζ a ( H ( t ) 𝐄 ( t ) χ a ) + 0 ( ( 1 - χ c ) 𝐖 a ( t ) ) )
= t λ ( 𝐅 a ( t ) + δ ( t ) χ a + [ 0 , 1 - χ c ] 𝐖 a ( t ) + ( 1 - χ c ) 0 𝐖 a ( t ) )
= t λ ( 𝐅 a ( t ) + δ ( t ) χ a + [ × ( × ) , χ c ] 𝐖 a ( t ) - ( 1 - χ c ) ( 1 - χ b ) 𝐅 a ( t ) )
= t λ ( χ b 𝐅 a ( t ) + δ ( t ) χ a + [ × ( × ) , χ c ] 𝐖 a ( t ) )
= χ a ( + 𝒦 a ( λ ) ) ,

where

𝒦 a ( λ ) = t λ ( χ b 𝐅 a ( t ) + [ × ( × ) , χ c ] 𝐖 a ( t ) ) .

We prove (2.1) by establishing the following estimates:

(A.5) χ R a # ( λ ) = 𝒪 ( e C ( Im λ ) - λ - 1 ) ,
(A.6) K a ( λ ) = 𝒪 ( e C ( Im λ ) - λ - N )

for any N. Then we have

(A.7) R ( λ ) χ a = R a # ( λ ) ( + 𝒦 a ( λ ) ) - 1 .

The invertibility of + 𝒦 a ( λ ) for λ Ω M is guaranteed by estimate (A.6).

First, we prove (A.5). Since 𝒫 is positive, we have

𝐄 ( t ) , t 𝐄 ( t ) = 𝒪 ( 1 ) .

Thus,

χ t λ ( ζ a H ( t ) 𝐄 ( t ) χ a ) = 𝒪 ( e C ( Im λ ) - )

and

i λ χ t λ ( ζ a H ( t ) 𝐄 ( t ) χ a ) = - t λ ( ζ a H ( t ) t 𝐄 ( t ) χ a ) = 𝒪 ( e C ( Im λ ) - ) .

Consequently, we have

χ t λ ( ζ a H ( t ) 𝐄 ( t ) χ a ) = 𝒪 ( e C ( Im λ ) - λ - 1 ) .

In view of (A.4), we have

t λ ( χ 𝐖 a ( t ) ) = 𝒪 ( e C ( Im λ ) - λ - N )

for any N.

Next, we prove (A.6). We only need to notice that χ b 𝐅 a ( t ) and [ × ( × ) , χ c ] 𝐖 a ( t ) both vanish outside a compact set in t. Hence, (2.1) is proved for α = 0 .

Now, we prove (2.1) for α = 1 2 . It suffices to prove

(A.8) χ R a # ( λ ) = 𝒪 ( e C ( Im λ ) - ) 𝒟 1 2

and use relation (A.7). Since

𝐄 ( t ) = sin t 𝒫 𝒫 ,

we have

𝐄 ( t ) 𝒟 1 2 C 𝐄 ( t ) + C 𝒫 𝐄 ( t ) = 𝒪 ( | t | + 1 ) .

Similarly, we obtain

χ t λ ( ζ a H ( t ) 𝐄 ( t ) χ a ) = 𝒪 ( e C ( Im λ ) - ) 𝒟 1 2 .

By

𝐖 a ( t ) 𝒟 1 2 = 𝒪 ( 1 )

and relation (A.4), we get

t λ ( χ 𝐖 a ( t ) ) = 𝒪 ( e C ( Im λ ) - λ - N ) 𝒟 1 2 .

Then (A.8) is proved.

For α = 1 , we take χ 1 C c ( 3 ) satisfying χ 1 = 1 near supp χ . A simple calculations gives

χ R ( λ ) χ 𝒟 C 𝒫 χ R ( λ ) χ + C χ R ( λ ) χ
C χ 𝒫 R ( λ ) χ + C [ 𝒫 , χ ] ( χ 1 R ( λ ) χ 1 ) χ + C χ R ( λ ) χ
| λ | 2 χ R ( λ ) χ + C χ 1 R ( λ ) χ 1 𝒟 1 2 + C χ R ( λ ) χ ,

which shows (2.1) for α = 1 .

References

[1] S. Acosta, S. Chow, J. Taylor and V. Villamizar, On the multi-frequency inverse source problem in heterogeneous media, Inverse Problems 28 (2012), no. 7, Article ID 075013. 10.1088/0266-5611/28/7/075013Suche in Google Scholar

[2] R. Albanese and P. B. Monk, The inverse source problem for Maxwell’s equations, Inverse Problems 22 (2006), no. 3, 1023–1035. 10.21236/ADA459256Suche in Google Scholar

[3] G. Alessandrini and A. Favaron, Sharp stability estimates of harmonic continuation along lines, Math. Methods Appl. Sci. 23 (2000), no. 12, 1037–1056. 10.1002/1099-1476(200008)23:12<1037::AID-MMA150>3.0.CO;2-BSuche in Google Scholar

[4] H. Ammari, G. Bao and J. L. Fleming, An inverse source problem for Maxwell’s equations in magnetoencephalography, SIAM J. Appl. Math. 62 (2002), no. 4, 1369–1382. 10.1137/S0036139900373927Suche in Google Scholar

[5] H. Ammari, J. Garnier and K. Sø lna, Resolution and stability analysis in full-aperture, linearized conductivity and wave imaging, Proc. Amer. Math. Soc. 141 (2013), no. 10, 3431–3446. 10.1090/S0002-9939-2013-11590-XSuche in Google Scholar

[6] G. Bao, C. Chen and P. Li, Inverse random source scattering problems in several dimensions, SIAM/ASA J. Uncertain. Quantif. 4 (2016), no. 1, 1263–1287. 10.1137/16M1067470Suche in Google Scholar

[7] G. Bao, C. Chen and P. Li, Inverse random source scattering for elastic waves, SIAM J. Numer. Anal. 55 (2017), no. 6, 2616–2643. 10.1137/16M1088922Suche in Google Scholar

[8] G. Bao, S.-N. Chow, P. Li and H. Zhou, An inverse random source problem for the Helmholtz equation, Math. Comp. 83 (2014), no. 285, 215–233. 10.1090/S0025-5718-2013-02730-5Suche in Google Scholar

[9] G. Bao and P. Li, Inverse medium scattering problems for electromagnetic waves, SIAM J. Appl. Math. 65 (2005), no. 6, 2049–2066. 10.1137/040607435Suche in Google Scholar

[10] G. Bao, P. Li, J. Lin and F. Triki, Inverse scattering problems with multi-frequencies, Inverse Problems 31 (2015), no. 9, Article ID 093001. 10.1088/0266-5611/31/9/093001Suche in Google Scholar

[11] G. Bao, P. Li and Y. Zhao, Stability for the inverse source problems in elastic and electromagnetic waves, J. Math. Pures Appl. (9) 134 (2020), 122–178. 10.1016/j.matpur.2019.06.006Suche in Google Scholar

[12] G. Bao, J. Lin and F. Triki, A multi-frequency inverse source problem, J. Differential Equations 249 (2010), no. 12, 3443–3465. 10.1016/j.jde.2010.08.013Suche in Google Scholar

[13] G. Bao, S. Lu, W. Rundell and B. Xu, A recursive algorithm for multifrequency acoustic inverse source problems, SIAM J. Numer. Anal. 53 (2015), no. 3, 1608–1628. 10.1137/140993648Suche in Google Scholar

[14] G. Bao and F. Triki, Stability for the multifrequency inverse medium problem, J. Differential Equations 269 (2020), no. 9, 7106–7128. 10.1016/j.jde.2020.05.021Suche in Google Scholar

[15] N. Bleistein and J. K. Cohen, Nonuniqueness in the inverse source problem in acoustics and electromagnetics, J. Math. Phys. 18 (1977), no. 2, 194–201. 10.1063/1.523256Suche in Google Scholar

[16] J. Cheng, V. Isakov and S. Lu, Increasing stability in the inverse source problem with many frequencies, J. Differential Equations 260 (2016), no. 5, 4786–4804. 10.1016/j.jde.2015.11.030Suche in Google Scholar

[17] A. J. Devaney, E. A. Marengo and M. Li, Inverse source problem in nonhomogeneous background media, SIAM J. Appl. Math. 67 (2007), no. 5, 1353–1378. 10.1137/060658618Suche in Google Scholar

[18] A. J. Devaney and G. C. Sherman, Nonuniqueness in inverse source and scattering problems, IEEE Trans. Antennas Propagation 30 (1982), no. 5, 1034–1042. 10.1109/TAP.1982.1142902Suche in Google Scholar

[19] S. Dyatlov and M. Zworski, Mathematical Theory of Scattering Resonances, Grad. Stud. Math. 200, American Mathematical Society, Providence, 2019. 10.1090/gsm/200Suche in Google Scholar

[20] A. El Badia and T. Nara, An inverse source problem for Helmholtz’s equation from the Cauchy data with a single wave number, Inverse Problems 27 (2011), no. 10, Article ID 105001. 10.1088/0266-5611/27/10/105001Suche in Google Scholar

[21] A. El Badia and T. Nara, Inverse dipole source problem for time-harmonic Maxwell equations: Algebraic algorithm and Hölder stability, Inverse Problems 29 (2013), no. 1, Article ID 15007. 10.1088/0266-5611/29/1/015007Suche in Google Scholar

[22] M. Eller and N. P. Valdivia, Acoustic source identification using multiple frequency information, Inverse Problems 25 (2009), no. 11, Article ID 115005. 10.1088/0266-5611/25/11/115005Suche in Google Scholar

[23] M. Entekhabi and V. Isakov, On increasing stability in the two dimensional inverse source scattering problem with many frequencies, Inverse Problems 34 (2018), no. 5, Article ID 055005. 10.1088/1361-6420/aab465Suche in Google Scholar

[24] M. Entekhabi and V. Isakov, Increasing stability in acoustic and elastic inverse source problems, SIAM J. Math. Anal. 52 (2020), no. 5, 5232–5256. 10.1137/19M1279885Suche in Google Scholar

[25] A. S. Fokas, Y. Kurylev and V. Marinakis, The unique determination of neuronal currents in the brain via magnetoencephalography, Inverse Problems 20 (2004), no. 4, 1067–1082. 10.1088/0266-5611/20/4/005Suche in Google Scholar

[26] K.-H. Hauer, L. Kühn and R. Potthast, On uniqueness and non-uniqueness for current reconstruction from magnetic fields, Inverse Problems 21 (2005), no. 3, 955–967. 10.1088/0266-5611/21/3/010Suche in Google Scholar

[27] V. Isakov and S. Kindermann, Subspaces of stability in the Cauchy problem for the Helmholtz equation, Methods Appl. Anal. 18 (2011), no. 1, 1–29. 10.4310/MAA.2011.v18.n1.a1Suche in Google Scholar

[28] A. Kirsch and F. Hettlich, The Mathematical Theory of Time-Harmonic Maxwell’s Equations Expansion-, Integral-, and Variational Methods, Appl. Math. Sci. 190, Springer, Cham, 2015. 10.1007/978-3-319-11086-8Suche in Google Scholar

[29] M. Li, C. Chen and P. Li, Inverse random source scattering for the Helmholtz equation in inhomogeneous media, Inverse Problems 34 (2018), no. 1, Article ID 015003. 10.1088/1361-6420/aa99d2Suche in Google Scholar

[30] P. Li, An inverse random source scattering problem in inhomogeneous media, Inverse Problems 27 (2011), no. 3, Article ID 035004. 10.1088/0266-5611/27/3/035004Suche in Google Scholar

[31] P. Li and G. Yuan, Increasing stability for the inverse source scattering problem with multi-frequencies, Inverse Probl. Imaging 11 (2017), no. 4, 745–759. 10.3934/ipi.2017035Suche in Google Scholar

[32] P. Li, J. Zhai and Y. Zhao, Stability for the acoustic inverse source problem in inhomogeneous media, SIAM J. Appl. Math. 80 (2020), no. 6, 2547–2559. 10.1137/20M1334267Suche in Google Scholar

[33] L. H. Nguyen, Q. Li and M. V. Klibanov, A convergent numerical method for a multi-frequency inverse source problem in inhomogenous media, Inverse Probl. Imaging 13 (2019), no. 5, 1067–1094. 10.3934/ipi.2019048Suche in Google Scholar

[34] G. Teschl, Mathematical Methods in Quantum Mechanics with Applications to Schrödinger Operators, 2nd ed., Grad. Stud. Math. 157, American Mathematical Society, Providence, 2014. Suche in Google Scholar

[35] D. Zhang and Y. Guo, Fourier method for solving the multi-frequency inverse source problem for the Helmholtz equation, Inverse Problems 31 (2015), no. 3, Article ID 035007. 10.1088/0266-5611/31/3/035007Suche in Google Scholar

Received: 2020-06-28
Revised: 2021-11-21
Accepted: 2022-03-07
Published Online: 2022-05-31
Published in Print: 2023-02-01

© 2022 Walter de Gruyter GmbH, Berlin/Boston

Heruntergeladen am 14.11.2025 von https://www.degruyterbrill.com/document/doi/10.1515/jiip-2020-0075/html
Button zum nach oben scrollen