Problem of Determining the Anisotropic Conductivity in Electrodynamic Equations

For a system of electrodynamic equations, the inverse problem of determining an anisotropic conductivity is considered. It is supposed that the conductivity is described by a diagonal matrix σ(x) = \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\text{diag}}({{\sigma }_{1}}(x),{{\sigma }_{2}}(x)$$\end{document}, σ3(x)) with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma (x) = 0$$\end{document} outside of the domain Ω = \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\{ x \in {{\mathbb{R}}^{3}}|\left| x \right| < R\} $$\end{document}, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R > 0$$\end{document}, and the permittivity ε and the permeability μ of the medium are positive constants everywhere in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\mathbb{R}}^{3}}$$\end{document}. Plane waves coming from infinity and impinging on an inhomogeneity localized in Ω are considered. For the determination of the unknown functions \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\sigma }_{1}}(x)$$\end{document}, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\sigma }_{2}}(x)$$\end{document}, and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\sigma }_{3}}(x)$$\end{document}, information related to the vector of electric intensity is given on the boundary S of the domain Ω. It is shown that this information reduces the inverse problem to three identical problems of X-ray tomography.

Consider the nonstationary Maxwell equations (1) Here, and are the electric and magnetic field strengths, σ(x) = , σ 2 (x), σ 3 (x)) is a positive semidefinite diagonal matrix, and ε and μ are positive constants. Assume that outside the domain Ω = , where R > 0.
Let denote the velocity of propagation of electromagnetic waves. Let , , and j be the unit vector orthogonal to ν, i.e., . For Maxwell's equations (1) in a homogeneous medium ( ), there exist solutions of the form (2) = ε + σ , = −μ , = . curl where and f(t) is an arbitrary generalized function. Each such solution represents a plane wave propagating in the direction of the vector ν and is a weak solution of Maxwell's equations for a homogeneous medium.
Consider the Cauchy problem for an anisotropic medium: (3) where and are defined by formulas (2) and is a smooth function such that for and . Thus, and for all and t < 0. Let be the plane corresponding to the front of a plane wave at the time t = 0, when this front touches the domain Ω.

Let
be the boundary of and be its shadow part with respect to light propagating in the direction ν.
Below, problem (3) is considered for three different vectors j k , , and for corresponding orthogonal vectors depending on the angular parameter ϕ, namely, , given for all , , and , where T k (x, ϕ) = -t 0 + δ 0 and is an arbitrary number (possibly small). In other words, the task is to find σ(x) from the given functions (5) For stationary electrodynamic equations, inverse problems of determining conductivity representing a one-variable function were studied by Tikhonov [1][2][3][4] and Cagniard [5]. For nonstationary equations, the theory of inverse electrodynamic problems based on the full system of Maxwell equations was developed in [6][7][8]. The problem of determining the permittivity of an anisotropic medium was considered in [9]. Additionally, phaseless inverse problems of determining permittivity from the magnitude of the electric or magnetic component of a stationary electromagnetic field were studied (see [10] and the review article in [11]). The following result holds for the inverse problem stated above. Theorem 1. Suppose that the matrix σ(x) belongs to and vanishes outside , while the function has the form , where is a smooth function of class , , and is the Heaviside step function, i.e., for and for t < 0. Then all elements of σ(x) in are uniquely determined by data (5).
The study of the inverse problem is based on analyzing the structure of the solution to problem (3). In this case, it is convenient to use an integro-differential equation for the vector . To derive it, we apply the curl operator to the second equation in (3) and use the first equation to eliminate the emerging term . Then we obtain the equation Computing with the help of the first equation in (3) yields (7) = , , , ν ϕ = , ϕ, ϕ , ϕ ∈ ,π, = , , , ν ϕ = ϕ, , ϕ , ϕ∈ ,π, = , , , ν ϕ = ϕ, ϕ, , ϕ∈ ,π.
It follows from (6), (7), and (3) that the function E is a solution of the Cauchy problem (8) The following result holds for problem (8). Theorem 2. Suppose that the matrix σ(x) and the function satisfy the conditions of Theorem 1. Then the function for can be represented in the form (9) where , the function is a solution of the Cauchy problem (10) and is a bounded function of x and t for for any T > 0.
Equation (10) is derived by substituting representation (9) into Eq. (8) and equating the coefficients of to zero, while the initial data for the function follow from the initial data for . The function is the amplitude of at the electromagnetic wave front, i.e., at . Equation (10) No. 1 2021 ROMANOV These formulas imply that the integrals (11) are known for all and . Thus, for each , the right-hand side of (11) is known along any straight line intersecting Ω in the direction . By varying ϕ, we conclude that, in each section of Ω by the plane , the integrals along all possible straight lines lying in this plane are known. As a result, we obtain an X-ray tomography problem for determining , . It is well known that this problem is uniquely solvable (see [12][13][14]). This implies Theorem 1 on the uniqueness of a solution to the inverse problem and an algorithm for its solution.