The problem of recovering the permittivity coefficient from the modulus of the scattered electromagnetic field

Under consideration is the stationary system of equations of electrodynamics relating to a nonmagnetic nonconducting medium. We study the problem of recovering the permittivity coefficient ε from given vectors of electric or magnetic intensities of the electromagnetic field. It is assumed that the field is generated by a point impulsive dipole located at some point y. It is also assumed that the permittivity differs from a given constant ε₀ only inside some compact domain Ω ? R³ with smooth boundary S. To recover ε inside Ω, we use the information on a solution to the corresponding direct problem for the system of equations of electrodynamics on the whole boundary of Ω for all frequencies from some fixed frequency ω ₀ on and for all y ∈ S. The asymptotics of a solution to the direct problem for large frequencies is studied and it is demonstrated that this information allows us to reduce the initial problem to the well-known inverse kinematic problem of recovering the refraction index inside Ω with given travel times of electromagnetic waves between two arbitrary points on the boundary of Ω. This allows us to state uniqueness theorem for solutions to the problem in question and opens up a way of its constructive solution.     

A two-dimensional inverse problem for an integro-differential equation of electrodynamics

An integro-differential equation corresponding to a two-dimensional problem of electrodynamics with dispersion is considered. It is assumed that the electrodynamic properties of a nonconducting medium with a constant magnetic permeability and the external current are independent of the x ₃ coordinate. In this case, the third component of the electric field vector satisfies a second-order scalar integro-differential equation with a variable permittivity of the medium. For this equation, we study the problem of finding the spatial part of the kernel entering the integral term. This corresponds to finding the part of the permittivity that depends on the electromagnetic frequency. It is assumed that the permittivity support is contained in some compact domain Ω ? ?². To find this coefficient inside Ω, we use information on the solution of the corresponding direct problem on the boundary of Ω on a finite time interval. An estimate for the conditional stability of the solution of the inverse problem is established under the assumption that the time interval is sufficiently large.     

Convergence of a family of solutions to a Fujita-type equation in domains with cavities

The Dirichlet problem for a Fujita-type equation, i.e., a second-order quasilinear uniformly elliptic equation is considered in domains Ω_ε with spherical or cylindrical cavities of characteristic size ε. The form of the function in the condition on the cavities’ boundaries depends on ε. For ε tending to zero and the number of cavities increasing simultaneously, sufficient conditions are established for the convergence of the family of solutions {u _ε(x)} of this problem to the solution u(х) of a similar problem in the domain Ω with no cavities with the same boundary conditions imposed on the common part of the boundaries ?Ω and ?Ω_ε. Convergence rate estimates are given.     

A stability estimate for a solution to an inverse problem of electrodynamics

We consider the integrodifferential system of equations of electrodynamics which corresponds to a dispersive nonmagnetic medium. For this system we study the problem of determining the spatial part of the kernel of the integral term. This corresponds to finding the part of dielectric permittivity depending nonlinearly on the frequency of the electromagnetic wave. We assume that the support of dielectric permittivity lies in some compact domain Ω ⊂ ℝ³. In order to find it inside Ω we start with known data about the solution to the corresponding direct problem for the equations of electrodynamics on the whole boundary of Ω for some finite time interval. On assuming that the time interval is sufficiently large we estimate the conditional stability of the solution to this inverse problem.     

Estimate of solution stability in a two-dimensional inverse problem for elasticity equations

The problem of determining the density of the medium and one of its elasticity moduli is considered. Properties of the elastic medium and external forces are assumed to be independent of the coordinate x ₃. In this case, the third component of the displacement vector satisfies a scalar equation of the second order, which contains the density ρ of the medium and elasticity modulus μ as coefficients. The parameters ρ and μ are known to be positive and constant everywhere outside some compact domain D ? ?², but they are unknown inside D. The problem of determining these coefficients in D via information, given on the boundary of the domain D for some finite time interval, about a solution of two direct problems is considered. An estimate of the conditional stability of a solution of the inverse problem under consideration is established.     

On the problem of determining the structure of a layered medium and the shape of an impulse source

For the equation of wave propagation in the half-space ? ₊ ² + = {(x, y) ∈ ?² | y > 0} we consider the problem of determining the speed of wave propagation that depends only on the variable y and the shape of a point impulse source on the boundary of the half-space. We show that, under some assumptions on the shape of the source and the structure of the medium, both unknown functions of one variable are uniquely determined by the displacements of boundary points of the medium. We estimate stability of a solution to the problem.     

Method of Integral Equations for the Three-Dimensional Problem of Wave Reflection from an Irregular Surface

The problem on the reflection of the field of a plane H-polarized three-dimensional electromagnetic wave from a perfectly conducting interface between media which contains a local perfectly conducting inhomogeneity is considered. To construct a numerical algorithm, the boundary value problem for the system of Maxwell equations in an infinite domain with irregular boundary is reduced to a system of singular integral equations, which is solved by the approximation–collocation method. The elements of the resulting complex matrix are calculated by a specially developed algorithm. The solution of the system of singular integral equations is used to obtain an integral representation for the reflected electromagnetic field and computational formulas for the directional diagram of the reflected electromagnetic field in the far region.     

Asymptotics of the Velocity Potential of an Ideal Fluid Flowing around a Thin Body

We consider the Neumann problem outside a small neighborhood of a planar disk in the three-dimensional space. The surface of this neighborhood is assumed to be smooth, and its thickness is characterized by a small parameter ε. A uniform asymptotic expansion of the solution of this problem with respect to ε is constructed by the matching method. Since the problem turned out to be bisingular, an additional inner asymptotic expansion in the so-called stretched variables is constructed near the edge of the disk. A physical interpretation of the solution of this boundary value problem is the velocity potential of a laminar flow of an ideal fluid around a thin body, which is the neighborhood of the disk. It is assumed that this flow has unit velocity at a large distance from the disk, which is equivalent to the following condition for the potential: u(x₁, x₂, x₃, ε) = x₃+O(r^?2) as r → ∞, where r is the distance to the origin. The boundary condition of this problem is the impermeability of the surface of the body: ?u/?n = 0 at the boundary. After subtracting x₃ from the solution u(x₁, x₂, x₃, ε), we get a boundary value problem for the potential ?(x₁, x₂, x₃, ε) of the perturbed motion. Since the integral of the function ??/?n over the surface of the body is zero, we have ?(x₁, x₂, x₃, ε) = O(r^?2) as r → ∞. Hence, all the coefficients of the outer asymptotic expansion with respect to ε have the same behavior at infinity. However, these coefficients have growing singularities at the approach to the edge of the disk, which implies the bisingularity of the problem.     

Uniform grid approximation of nonsmooth solutions to the mixed boundary value problem for a singularly perturbed reaction-diffusion equation in a rectangle

A mixed boundary value problem for a singularly perturbed reaction-diffusion equation in a square is considered. A Neumann condition is specified on one side of the square, and a Dirichlet condition is set on the other three. It is assumed that the coefficient of the equation, its right-hand side, and the boundary values of the desired solution or its normal derivative on the sides of the square are smooth enough to ensure the required smoothness of the solution in a closed domain outside the neighborhoods of the corner points. No compatibility conditions are assumed to hold at the corner points. Under these assumptions, the desired solution in the entire closed domain is of limited smoothness: it belongs only to the Hölder class C ^μ, where μ ∈ (0, 1) is arbitrary. In the domain, a nonuniform rectangular mesh is introduced that is refined in the boundary domain and depends on a small parameter. The numerical solution to the problem is based on the classical five-point approximation of the equation and a four-point approximation of the Neumann boundary condition. A mesh refinement rule is described under which the approximate solution converges to the exact one uniformly with respect to the small parameter in the L _∞ ^h norm. The convergence rate is O(N ^?2ln² N), where N is the number of mesh nodes in each coordinate direction. The parameter-uniform convergence of difference schemes for mixed problems without compatibility conditions at corner points was not previously analyzed.     

Two-Parameter Asymptotics in a Bisingular Cauchy Problem for a Parabolic Equation

The Cauchy problem for a quasilinear parabolic equation with a small parameter ε at the highest derivative is considered. The initial function, which has the form of a smoothed step, depends on a “stretched” variable x/ρ, where ρ is another small parameter. This problem statement is of interest for applications as a model of propagation of nonlinear waves in physical systems in the presence of small dissipation. In the case corresponding to a compression wave, asymptotic solutions of the problem are constructed in the parameters ε and ρ independently tending to zero. It is assumed that ε/ρ → 0. Far from the line of discontinuity of the limit solution, asymptotic solutions are constructed in the form of series in powers of ε and ρ. In a small domain of linear approximation, an asymptotic solution is constructed in the form of a series in powers of the ratio ρ/ε. The coefficients of the inner expansion are determined from a recursive chain of initial value problems. The asymptotics of these coefficients at infinity is studied. The time of reconstruction of the scale of the internal space variable is determined.     

Singularly Perturbed Parabolic Problems with Multidimensional Boundary Layers

The first boundary value problem for a multidimensional parabolic differential equation with a small parameter ε multiplying all derivatives is studied. A complete (i.e., of any order with respect to the parameter) regularized asymptotics of the solution is constructed, which contains a multidimensional boundary layer function that is bounded for x = (x¹, x₂) = 0 and tends to zero as ε → +0 for x ≠ 0. In addition, it contains corner boundary layer functions described by the product of a boundary layer function of the exponential type by a multidimensional parabolic boundary layer function.     

The Problem of Finding the One-Dimensional Kernel of the Thermoviscoelasticity Equation

The problem of determining the kernel h(t), t ∈ [0, T], appearing in the system of integro-differential thermoviscoelasticity equations is considered. It is assumed that the coefficients of the equations depend only on one space variable. The inverse problem is replaced by the equivalent system of integral equations for unknown functions. The contraction mapping principle with weighted norms is applied to this system in the space of continuous functions. A global unique solvability theorem is proved and an estimate of the stability of the solution of the inverse problem is obtained.     

Difference scheme for an initial–boundary value problem for a singularly perturbed transport equation

An initial–boundary value problem for a singularly perturbed transport equation with a perturbation parameter ε multiplying the spatial derivative is considered on the set ? = G ∪ S, where ? = D? × [0 ≤ t ≤ T], D? = {0 ≤ x ≤ d}, S = S ^l ∪ S, and S ^l and S₀ are the lateral and lower boundaries. The parameter ε takes arbitrary values from the half-open interval (0,1]. In contrast to the well-known problem for the regular transport equation, for small values of ε, this problem involves a boundary layer of width O(ε) appearing in the neighborhood of S ^l ; in the layer, the solution of the problem varies by a finite value. For this singularly perturbed problem, the solution of a standard difference scheme on a uniform grid does not converge ε-uniformly in the maximum norm. Convergence occurs only if h=dN^-1 ? ε and N₀^-1 ? 1, where N and N₀ are the numbers of grid intervals in x and t, respectively, and h is the mesh size in x. The solution of the considered problem is decomposed into the sum of regular and singular components. With the behavior of the singular component taken into account, a special difference scheme is constructed on a Shishkin mesh, i.e., on a mesh that is piecewise uniform in x and uniform in t. On such a grid, a monotone difference scheme for the initial–boundary value problem for the singularly perturbed transport equation converges ε-uniformly in the maximum norm at an ?(N^?1 + N₀^?1) rate.     

Homogenization of the stationary periodic Maxwell system in the case of constant permeability

The homogenization problem in the small period limit for the stationary periodic Maxwell system in ℝ³ is considered. It is assumed that the permittivity η^ε(x)=η(ε⁻x), ε > 0, is a rapidly oscillating positive matrix function and the permeability μ₀ is a constant positive matrix. For all four physical fields (the electric and magnetic field intensities, the electric displacement field, and the magnetic flux density), we obtain uniform approximations in the L ₂(ℝ³)-norm with order-sharp remainder estimates. __________ Translated from Funktsional’nyi Analiz i Ego Prilozheniya, Vol. 41, No. 2, pp. 3–23, 2007 Original Russian Text Copyright ? by M. Sh. Birman and T. A. Suslina Dedicated to the memory of the great mathematician Mark Grigor’evich Krein Supported by RFBR grants no. 05-01-01076-a, 05-01-02944-YaF-a.     

An overfilled cavity problem for Maxwell's equations

This paper is concerned with the mathematical analysis of the scattering of a time‐harmonic electromagnetic plane wave by an open and overfilled cavity that is embedded in a perfect electrically conducting infinite ground plane, where the electromagnetic wave propagation is governed by the Maxwell equations. Above the flat ground surface and the open aperture of the cavity, the space is assumed to be filled with a homogeneous medium with a constant permittivity and permeability, whereas the interior of the cavity is filled with some inhomogeneous medium with a variable permittivity and permeability. The scattering problem is modeled as a boundary value problem over a bounded domain, with transparent boundary condition proposed on the hemisphere enclosing the inhomogeneity represented by the cavity. The existence and uniqueness of the weak solution for the model problem are established by using a variational approach. The perfectly matched layer (PML) method is investigated to truncate the unbounded electromagnetic cavity scattering problem. It is shown that the truncated PML problem attains a unique solution. An explicit error estimate is given between the solution of the original scattering problem and that of the truncated PML problem. The error estimate implies that the PML solution converges exponentially to the original cavity scattering problem by increasing either the PML medium parameter or the PML layer thickness. The convergence result is expected to be useful for determining the PML medium parameter in the computational electromagnetic scattering problem. Copyright © 2012 John Wiley & Sons, Ltd.     

A spectral method for elliptic equations: the Dirichlet problem

Let Ω be an open, simply connected, and bounded region in ? ^d , d?≥?2, and assume its boundary \(\partial\Omega\) is smooth. Consider solving an elliptic partial differential equation Lu?=?f over Ω with zero Dirichlet boundary values. The problem is converted to an equivalent elliptic problem over the unit ball B; and then a spectral Galerkin method is used to create a convergent sequence of multivariate polynomials u _n of degree ≤?n that is convergent to u. The transformation from Ω to B requires a special analytical calculation for its implementation. With sufficiently smooth problem parameters, the method is shown to be rapidly convergent. For \(u\in C^{\infty}( \overline{\Omega})\) and assuming \(\partial\Omega\) is a C ^?∞? boundary, the convergence of \(\left\Vert u-u_{n}\right\Vert _{H^{1}}\) to zero is faster than any power of 1/n. Numerical examples in ?² and ?³ show experimentally an exponential rate of convergence.     

Computer difference scheme for a singularly perturbed elliptic convection–diffusion equation in the presence of perturbations

A grid approximation of a boundary value problem for a singularly perturbed elliptic convection–diffusion equation with a perturbation parameter ε, ε ∈ (0,1], multiplying the highest order derivatives is considered on a rectangle. The stability of a standard difference scheme based on monotone approximations of the problem on a uniform grid is analyzed, and the behavior of discrete solutions in the presence of perturbations is examined. With an increase in the number of grid nodes, this scheme does not converge -uniformly in the maximum norm, but only conditional convergence takes place. When the solution of the difference scheme converges, which occurs if N ₁ ^-1 N ₂ ^-1 ? ε, where N ₁ and N ₂ are the numbers of grid intervals in x and y, respectively, the scheme is not -uniformly well-conditioned or ε-uniformly stable to data perturbations in the grid problem and to computer perturbations. For the standard difference scheme in the presence of data perturbations in the grid problem and/or computer perturbations, conditions imposed on the “parameters” of the difference scheme and of the computer (namely, on ε, N ₁,N ₂, admissible data perturbations in the grid problem, and admissible computer perturbations) are obtained that ensure the convergence of the perturbed solutions as N ₁,N ₂ → ∞, ε ∈ (0,1]. The difference schemes constructed in the presence of the indicated perturbations that converges as N ₁,N ₂ → ∞ for fixed ε, ε ∈ (0,1, is called a computer difference scheme. Schemes converging ε-uniformly and conditionally converging computer schemes are referred to as reliable schemes. Conditions on the data perturbations in the standard difference scheme and on computer perturbations are also obtained under which the convergence rate of the solution to the computer difference scheme has the same order as the solution of the standard difference scheme in the absence of perturbations. Due to this property of its solutions, the computer difference scheme can be effectively used in practical computations.     

Equivalence of two definitions of a generalized <Emphasis Type="Italic">L</Emphasis><Subscript><Emphasis Type="Italic">p</Emphasis></Subscript> solution to the initial-boundary value problem for the wave equation

In our previous papers, we introduced the notion of a generalized solution to the initial-boundary value problem for the wave equation with a boundary function µ(t) such that the integral ∫ ₀ ^T (T ? t)|µ(t)| ^p dt exists. Here we prove that this solution is a unique solution to the problem in L _p that satisfies the corresponding integral identity.     

Mixed problem for the wave equation with integrable potential in the case of two-point boundary conditions of distinct orders

We study a mixed problem for the wave equation with integrable potential and with two-point boundary conditions of distinct orders for the case in which the corresponding spectral problem may have multiple spectrum. Based on the resolvent approach in the Fourier method and the Krylov convergence acceleration trick for Fourier series, we obtain a classical solution u(x, t) of this problem under minimal constraints on the initial condition u(x, 0) = ?(x). We use the Carleson–Hunt theorem to prove the convergence almost everywhere of the formal solution series in the limit case of ?(x) ∈ L _p[0, 1], p > 1, and show that the formal solution is a generalized solution of the problem.     

Finite difference schemes for the singularly perturbed reaction-diffusion equation in the case of spherical symmetry

The boundary value problem for the singularly perturbed reaction-diffusion parabolic equation in a ball in the case of spherical symmetry is considered. The derivatives with respect to the radial variable appearing in the equation are written in divergent form. The third kind boundary condition, which admits the Dirichlet and Neumann conditions, is specified on the boundary of the domain. The Laplace operator in the differential equation involves a perturbation parameter ?², where ? takes arbitrary values in the half-open interval (0, 1]. When ? → 0, the solution of such a problem has a parabolic boundary layer in a neighborhood of the boundary. Using the integro-interpolational method and the condensing grid technique, conservative finite difference schemes on flux grids are constructed that converge ?-uniformly at a rate of O(N ^?2ln² N + N ₀ ^?1 ), where N + 1 and N ₀ + 1 are the numbers of the mesh points in the radial and time variables, respectively.