
The current manuscript focuses on the solution and analysis of space and time fractional coupled Schrödinger system that belongs to a class of evolution equations. These systems encounter in different fields related to plasma waves, optics, and quantum physics. The fractional He-Laplace approach is proposed for the series form solutions of fractional systems. This approach contains hybrid of Laplace transform and homotopy perturbation along with Caputo fractional derivative. The current study provide new results on time and space fractional coupled Schrödinger systems which are not captured in existing literature. Reliability of proposed algorithm in both time and space fractional scenarios is observed through residual error concept throughout fractional domains. The effect of fractional parameters on wave profiles are analyzed numerically and graphically as 2D and 3D illustrations. Analysis reveals that proposed algorithm is suitable for non-linear time-space fractional systems encountering in different fields of sciences.
Citation: Mubashir Qayyum, Efaza Ahmad, Hijaz Ahmad, Bandar Almohsen. New solutions of time-space fractional coupled Schrödinger systems[J]. AIMS Mathematics, 2023, 8(11): 27033-27051. doi: 10.3934/math.20231383
[1] | Wei Fan, Kangqun Zhang . Local well-posedness results for the nonlinear fractional diffusion equation involving a Erdélyi-Kober operator. AIMS Mathematics, 2024, 9(9): 25494-25512. doi: 10.3934/math.20241245 |
[2] | Min Jiang, Rengang Huang . Existence of solutions for q-fractional differential equations with nonlocal Erdélyi-Kober q-fractional integral condition. AIMS Mathematics, 2020, 5(6): 6537-6551. doi: 10.3934/math.2020421 |
[3] | Kangqun Zhang . Existence and uniqueness of positive solution of a nonlinear differential equation with higher order Erdélyi-Kober operators. AIMS Mathematics, 2024, 9(1): 1358-1372. doi: 10.3934/math.2024067 |
[4] | Dumitru Baleanu, S. Hemalatha, P. Duraisamy, P. Pandiyan, Subramanian Muthaiah . Existence results for coupled differential equations of non-integer order with Riemann-Liouville, Erdélyi-Kober integral conditions. AIMS Mathematics, 2021, 6(12): 13004-13023. doi: 10.3934/math.2021752 |
[5] | Kishor D. Kucche, Sagar T. Sutar, Kottakkaran Sooppy Nisar . Analysis of nonlinear implicit fractional differential equations with the Atangana-Baleanu derivative via measure of non-compactness. AIMS Mathematics, 2024, 9(10): 27058-27079. doi: 10.3934/math.20241316 |
[6] | Mohamed Jleli, Bessem Samet . Nonexistence for fractional differential inequalities and systems in the sense of Erdélyi-Kober. AIMS Mathematics, 2024, 9(8): 21686-21702. doi: 10.3934/math.20241055 |
[7] | Wedad Albalawi, Muhammad Imran Liaqat, Kottakkaran Sooppy Nisar, Abdel-Haleem Abdel-Aty . Qualitative study of Caputo Erdélyi-Kober stochastic fractional delay differential equations. AIMS Mathematics, 2025, 10(4): 8277-8305. doi: 10.3934/math.2025381 |
[8] | Choukri Derbazi, Hadda Hammouche . Caputo-Hadamard fractional differential equations with nonlocal fractional integro-differential boundary conditions via topological degree theory. AIMS Mathematics, 2020, 5(3): 2694-2709. doi: 10.3934/math.2020174 |
[9] | Bashir Ahmad, Manal Alnahdi, Sotiris K. Ntouyas, Ahmed Alsaedi . On a mixed nonlinear boundary value problem with the right Caputo fractional derivative and multipoint closed boundary conditions. AIMS Mathematics, 2023, 8(5): 11709-11726. doi: 10.3934/math.2023593 |
[10] | Manal Elzain Mohamed Abdalla, Hasanen A. Hammad . Solving functional integrodifferential equations with Liouville-Caputo fractional derivatives by fixed point techniques. AIMS Mathematics, 2025, 10(3): 6168-6194. doi: 10.3934/math.2025281 |
The current manuscript focuses on the solution and analysis of space and time fractional coupled Schrödinger system that belongs to a class of evolution equations. These systems encounter in different fields related to plasma waves, optics, and quantum physics. The fractional He-Laplace approach is proposed for the series form solutions of fractional systems. This approach contains hybrid of Laplace transform and homotopy perturbation along with Caputo fractional derivative. The current study provide new results on time and space fractional coupled Schrödinger systems which are not captured in existing literature. Reliability of proposed algorithm in both time and space fractional scenarios is observed through residual error concept throughout fractional domains. The effect of fractional parameters on wave profiles are analyzed numerically and graphically as 2D and 3D illustrations. Analysis reveals that proposed algorithm is suitable for non-linear time-space fractional systems encountering in different fields of sciences.
Generally, one of the purposes of the microwave imaging is to localize or detect unknown objects from measured electromagnetic waves in the high-frequency regime (between 300MHz and 300GHz). Due to this reason, identification of the location and shape of unknown anomaly whose values of dielectric permittivity and electric conductivity differ from the homogeneous background medium from scattering parameter data is an important and interesting research topic in microwave imaging. Many authors have proposed various remarkable algorithms for retrieving information on anomaly for example, Born iterative method for reconstructing permittivity distribution [13] and brain stroke detection [20], level-set method for shape reconstruction of unknown objects [16] and breast cancer detection [21], distorted iterated virtual experiments scheme for imaging unknown scatterers [28], conjugate gradient method for breast imaging [11], Levenberg-Marquardt technique for recovering parameter distribution [17]. We also refer to remarkable mathematical and experimental studies [1,3,12,15,26,45,46,47,48]. However, the success of iteration-based algorithms significantly depends on the priori information and initial guess, which must be close to the unknown anomaly, refer to [25,39].
For this reason, various non-iterative schemes in inverse scattering problems without a priori information about unknown anomaly have also been developed to retrieve the location and shape, such as a variational algorithm based on the inverse Fourier transform to retrieve small electromagnetic inhomogeneities [5,6], direct sampling method for localizing small electromagnetic inhomogeneities [22,23] and anomaly detection in real-world experiment [43], factorization method for crack detection [10], shape reconstruction of unknown obstacles [18] and numerical study for anomaly imaging [33], linear sampling method for imaging of unknown scatters in limited-aperture problem [8] and crack-like defects [24], MUltiple SIgnal Classification (MUSIC) algorithm for identifying small anomalies [38], fast imaging of small targets in limited-aperture measurement configuration [35] and real-world application of anomaly detection [34], orthogonality sampling method for imaging unknown targets [40], qualitative microwave imaging [2], anomaly detection in microwave imaging [37] and topological derivative strategy for imaging crack-like defects [31] and retrieving unknown scatterers in 3D [27].
The subspace migration (SM) algorithm is a well-known, non-iterative imaging technique in both inverse scattering problem and microwave imaging. It has been applied successfully to the various problems such as localization of small targets [4], identification of extended objects [9], fast imaging of curve-like cracks [29,30] and anomaly detection in microwave imaging [32,36]. Throughout several studies, it has been confirmed that SM is a fast, robust and effective technique for retrieving unknown anomaly from scattering parameter data. However, for successful application in microwave imaging, accurate values of background permittivity and conductivity must be known because the exact background wavenumber value must be applied. Generally, most researchers have used the statistical values of the background material instead of the true ones and obtained inaccurate locations and shapes of anomaly. This can be examined through the results of numerical simulations but no reliable mathematical theory explaining this phenomenon has yet been developed.
In this paper, we apply the SM to retrieve unknown anomaly from scattering parameter data without complete information about the background material; that is, inaccurate values of background permittivity and conductivity are applied. To explain the appearance of an inaccurate location and shape of an anomaly, we show that the imaging function of the SM can be written as the infinite series of Bessel function of integer order, antenna arrangement and applied inaccurate values of background permittivity and conductivity. This enables us to theoretically explain the appearance of an inaccurate location and shape of an anomaly through the SM. To confirm the theoretical results, simulation results with inaccurate values of background permittivity and conductivity are also presented.
The rest of this paper is organized as follows. In Section 2, the two-dimensional direct problem and imaging function of the SM are introduced. The structure of the imaging function with inaccurate values of background permittivity and conductivity is revealed in Section 3. In Section 4, a set of numerical simulation results with synthetic data generated by CST STUDIO SUITE is presented. Finally, a short conclusion is provided in Section 5.
Let D be a circle-like anomaly with radius α and center r⋆ that is surrounded by a circular array of antennas An with location an and |an|=R, n=1,2,⋯,N(>2). Throughout this paper, we use Ω to denote a homogeneous domain filled by matching liquid such that D⊂Ω and assume that all materials D and Ω are nonmagnetic and classified by their value of dielectric permittivity and electric conductivity at a given angular frequency of operation ω, i.e., the value of magnetic permeability is constant for every r∈Ω say, μ(r)≡μb=1.257×10−6H/m, refer to [41]. We use εb and σb to denote the background permittivity and conductivity, respectively. Analogously, let ε⋆ and σ⋆ be those of D. Then, we introduce the piecewise constant permittivity ε(r) and conductivity σ(r),
ε(r)={ε⋆forr∈D,εbforr∈Ω∖D,andσ(r)={σ⋆forr∈D,σbforr∈Ω∖D, |
respectively. With this, let kb be the background wavenumber that satisfies
k2b=ω2μb(εb−iσbω), |
and further assume that
ωεb≫σband√ε⋆εb<1+wavelength4α. | (2.1) |
With this, the time-dependent, homogeneous, linear Maxwell Equations take the form:
curl E(r,t)=−μb∂H(r,t)∂tandcurl H(r,t)=σbE(r,t)+εb∂E(r,t)∂t,r∈Ω, | (2.2) |
where E(r,t)∈R3 and H(r,t)∈R3 are the electric and magnetic fields, respectively. Here, we consider time-harmonic solutions to the (2.2) such that
E(r,t)=Re[E(r)e−iωt]andH(r,t)=Re[H(r)e−iωt],r∈Ω,t>0. |
Then, E(r)∈C3 and H(r)∈C3 satisfy
curl E(r)=iωμbH(r)andcurl H(r)=(σb−iωεb)EinΩ. | (2.3) |
For a detailed description, we refer to [7].
Let Einc(kb,r,am) be the incident electric field due to the point current density J at Am. Then, based on (2.3), it satisfies
{curl Einc(kb,r,am)=iωμbHinc(kb,am,r),curl Hinc(kb,r,am)=(σb−iωεb)Einc(kb,r,am), |
Analogously, let Etot(an,r) be the total field measured at An in the presence of D that satisfies
{curl Etot(kb,an,r)=iωμbHtot(kb,an,r),curl Htot(kb,an,r)=(σ(r)−iωε(r))Etot(kb,an,r), |
with transmission condition on ∂D.
Let S(n,m) be the S-parameter (or scattering parameter), which is defined as the ratio of the output voltage (or reflected waves) at the An antenna and the input voltage (or incident waves) at the Am (see [41] for instance). We use Sinc(n,m) and Stot(n,m) to denote the incident and total field S-parameters, respectively, due to the absence and presence of D. Correspondingly, we let Sscat(n,m)=Stot(n,m)−Sinc(n,m) be the scattered field S-parameter. Then, based on [19], Sscat(n,m) can be represented as follows:
Sscat(n,m)=−ik204ωμb∫Ω(ε(r′)−εbεb+iσ(r′)−σbωεb)Einc(kb,r′,am)⋅Etot(kb,an,r′)dr′, |
where k0 denotes the lossless background wave number that satisfies k20=ω2μbεb.
In this paper, we adopt the simulation configuration introduced in [32,43,44]. Notice that the height of microwave machine can be said to be long enough, only the z-component of the incident and total fields can be handled based on the mathematical treatment of the scattering of time-harmonic electromagnetic waves from thin infinitely long cylindrical obstacles. Correspondingly, by denoting E(z)inc and E(z)tot as the z-components of the incident and total fields, respectively, Sscat(n,m) can be written as follows:
Sscat(n,m)=−ik204ωμb∫Ω(ε(r′)−εbεb+iσ(r′)−σbωεb)E(z)inc(kb,r′,am)E(z)tot(kb,an,r′)dr′. | (2.4) |
Unfortunately, exact expression of the field E(z)tot(kb,an,r′) is unknown thus we cannot design the imaging function by using Sscat(n,m) of (2.4) directly. Since the condition (2.1) holds, we can apply the Born approximation such that (see [42] for instance)
E(z)tot(kb,an,r′)=E(z)inc(kb,an,r′)+o(α2)=−i4H(1)0(kb|a−r′|)+o(α2), |
where H(1)0 denotes the Hankel function of order zero of the first kind. Correspondingly, Sscat(n,m) of (2.4) can be written as
Sscat(n,m)=−ik204ωμb∫D(ε⋆−εbεb+iσ⋆−σbωεb)E(z)inc(kb,r′,am)E(z)tot(kb,an,r′)dr′+o(α2k20)=ik2064ωμb∫D(ε⋆−εbεb+iσ⋆−σbωεb)H(1)0(kb|am−r′|)H(1)0(kb|an−r′|)dr′+o(α2k20). | (2.5) |
Now, we introduce the imaging function. To this end, let us generate the scattering matrix K such that
K=[skip0Sscat(1,2)⋯Sscat(1,N−1)Sscat(1,N)Sscat(2,1)0⋯Sscat(2,N−1)Sscat(2,N)⋮⋮⋱⋮⋮Sscat(N,1)Sscat(N,2)⋯Sscat(N,N−1)0]. | (2.6) |
See [32] for an explanation of why the diagonal elements of K are set to zero. Since there is a single, small anomaly, the singular value decomposition (SVD) of K can be written as
K=UDV∗=N∑n=1τnUnV∗n≈τ1U1V∗1, | (2.7) |
where ∗ denotes the mark of Hermitian, τn are the singular values and Un and Vn are the left and right singular vectors of K, respectively. Then, on the basis of (2.5) and (2.7), we define a unit vector: For each r∈Ω,
W(kb,r)=F(kb,r)|F(kb,r)|,F(kb,r)=[H(1)0(kb|a1−r|),H(1)0(kb|a2−r|),…,H(1)0(kb|aN−r|)]T. | (2.8) |
With this, we introduce the following imaging function of the SM: For each r∈Ω,
F(kb,r)=|⟨W(kb,r),U1⟩⟨W(kb,r),¯V1⟩|, | (2.9) |
where ⟨U,V⟩=U∗V and ¯V1 denotes the complex conjugate of V1. Then, based on [4]
⟨W(kb,r),U1⟩≈1and⟨W(kb,r),¯V1⟩≈1whenr∈D, |
and the orthonormal property of singular vectors, the value of F(kb,r) will be close to 1 when r∈D and less than 1 at r∈Ω∖D, so the location and outline shape of D can be identified through the map of F(kb,r).
Let us emphasize that, to generate the test vector W(kb,r) of (2.8), the exact value of kb must be known, i.e., a priori information of the εb and σb must be available. However, because these values are statistical, the exact value may not be unknown. For this reason, we assume that the exact values of εb and σb are unknown and apply an alternative value ka instead of the true kb and set a unit test vector W(ka,r) from (2.8). Then, by using the imaging function F(ka,r) from (2.9), the exact location and shape of D cannot be retrieved. Fortunately, we can recognize the existence of D and the identified location is shifted in a specific direction.
In this section, we explore the structure of the imaging function F(ka,r) to explain that retrieved location of D is shifted in a specific direction and size is smaller or larger than the true anomaly. To explain this phenomenon, we explore the structure of the imaging function.
Theorem 3.1. Let θn=an/|an|=(cosθn,sinθn) and kbr′−kar=|kbr′−kar|(cosϕ,sinϕ). If an satisfies |an−r|≫{1/4|ka|,1/4|kb|} for all n=1,2,⋯,N, then F(r,ka) can be represented as follows:
F(ka,r)=N(N−1)area(D)|∫D(J0(|kbr′−kar|)+Ψ(kb,ka,r)N)2dr′−1N∫D(J0(2|kbr′−kar|)+Ψ(2kb,2ka,r)N)dr′|+o(α2k20), | (3.1) |
where Js denotes the Bessel function of order s and
Ψ(kb,ka,r)=N∑n=1∞∑s=−∞,s≠0isJs(|kbr′−kar|)eis(θn−ϕ). |
Proof. Since K≈τ1U1V∗1, we can examine that
F(kb,r)=|⟨W(ka,r),U1⟩⟨W(ka,r),¯V1⟩|=|W(ka,r)∗U1V∗1¯W(ka,r)|≈|1τ1W(ka,r)∗K¯W(ka,r)|. |
Based on the assumption |an−r|≫{1/4|ka|,1/4|kb|} for all n=1,2,⋯,N, the following asymptotic forms of the Hankel function hold (see [14,Theorem 2.5], for instance)
H(1)0(kb|r−r′|)=(1−i)e−ikb|r|√kbπ|an|eikbθn⋅r′+O(1)andH(1)0(ka|r−r′|)=(1−i)e−ika|r|√kaπ|an|eikaθn⋅r′+O(1). | (3.2) |
Then, W(ka,r) and K can be represented as
W(ka,r)=1√N[eikaθ1⋅r+O(1)eikaθ2⋅r+O(1)⋮eikaθN⋅r+O(1)] |
and
K=C[0∫Deikb(θ1+θ2)⋅r′dr′+o(α2k20)⋯∫Deikb(θ1+θN)⋅r′dr′+o(α2k20)∫Deikb(θ2+θ1)⋅r′dr′+o(α2k20)0⋯∫Deikb(θ2+θN)⋅r′dr′+o(α2k20)⋮⋮⋱⋮∫Deikb(θN+θ1)⋅r′dr′+o(α2k20)∫Deikb(θN+θ2)⋅r′dr′+o(α2k20)⋯0], |
respectively. Here the constant C is given by
C=k20e−2ikbR32kbωμbπR(ε⋆−εbεb+iσ⋆−σbωεb). |
Note that, since the following JacobiAnger expansion formula holds uniformly
eixcosθ=∞∑s=−∞isJs(x)eisθ=J0(x)+∞∑s=−∞,s≠0isJs(x)eisθ, | (3.3) |
we have for n=1,2,⋯,N
N∑n=1eiθn⋅(kbr′−kar)=N∑n=1ei|kbr′−kar|cos(θn−ϕ)=N∑n=1(J0(|kbr′−kar|)+∞∑s=−∞,s≠0isJs(|kbr′−kar|)eis(θn−ϕ))=NJ0(|kbr′−kar|)+Ψ(kb,ka,r) |
and correspondingly,
W(ka,r)∗K=C√N[∫Deikbθ1⋅r′N∑n=1(eiθn(kbr′−kar)−eiθ1(kbr′−kar))dr′+o(α2k20)∫Deikbθ2⋅r′N∑n=1(eiθn(kbr′−kar)−eiθ2(kbr′−kar))dr′+o(α2k20)⋮∫Deikbθ1⋅r′N∑n=1(eiθn(kbr′−kar)−eiθN(kbr′−kar))dr′+o(α2k20)]T=C√N[∫Deikbθ1⋅r′(NJ0(|kbr′−kar|)+Ψ(kb,ka,r)−eiθ1(kbr′−kar))dr′+o(α2k20)∫Deikbθ2⋅r′(NJ0(|kbr′−kar|)+Ψ(kb,ka,r)−eiθ2(kbr′−kar))dr′+o(α2k20)⋮∫Deikbθ1⋅r′(NJ0(|kbr′−kar|)+Ψ(kb,ka,r)−eiθN(kbr′−kar))dr′+o(α2k20)]T. |
With this, we can evaluate
W(ka,r)∗K¯W(ka,r)=CN[∫Deikbθ1⋅r′(NJ0(|kbr′−kar|)+Ψ(kb,ka,r)−eiθ1(kbr′−kar))dr′+o(α2k20)∫Deikbθ2⋅r′(NJ0(|kbr′−kar|)+Ψ(kb,ka,r)−eiθ2(kbr′−kar))dr′+o(α2k20)⋮∫DeikbθN⋅r′(NJ0(|kbr′−kar|)+Ψ(kb,ka,r)−eiθN(kbr′−kar))dr′+o(α2k20)]T[e−ikaθ1⋅r+O(1)e−ikaθ2⋅r+O(1)⋮e−ikaθN⋅r+O(1)]=CN∫DN∑n=1eiθn⋅(kbr′−kar)(NJ0(|kbr′−kar|)+Ψ(kb,ka,r)−eiθn(kbr′−kar))dr′+o(α2k20)=CN∫Deiθn⋅(kbr′−kar)(NJ0(|kbr′−kar|)+Ψ(kb,ka,r))dr′−CN∫DN∑n=1e2iθn⋅(kbr′−kar)dr′+o(α2k20)=CN∫D[(NJ0(|kbr′−kar|)+Ψ(kb,ka,r))2dr′−(NJ0(2|kbr′−kar|)+Ψ(2kb,2ka,r))]dr′+o(α2k20). |
Hence,
⟨W(ka,r),U1⟩⟨W(ka,r),¯V1⟩=CNτ1∫D(J0(|kbr′−kar|)+Ψ(kb,ka,r)N)2dr′−Cτ1∫D(J0(2|kbr′−kar|)+Ψ(2kb,2ka,r)N)dr′+o(α2k20). |
Since ⟨W(r),U1⟩⟨W(r),¯V1⟩=1, J0(|kbr′−kar|)=J0(2|kbr′−kar|)=1 and Ψ(kb,ka,r)=Ψ(2kb,2ka,r)=0 when kbr′=kar, we have
CNτ1∫Ddr′−Cτ1∫Ddr′+o(α2k20)=1impliesC=τ1(N−1)area(D)+o(α2k20). |
Therefore,
⟨W(r),U1⟩⟨W(r),¯V1⟩=N(N−1)area(D)∫D(J0(|kbr′−kar|)+Ψ(kb,ka,r)N)2dr′−1(N−1)area(D)∫D(J0(2|kbr′−kar|)+Ψ(2kb,2ka,r)N)dr′+o(α2k20). |
With this, we can obtain the structure (3.1).
From the derived structure (3.1), we can observe that since J0(|kbr′−kar|)=1 and Ψ(kb,ka,r)=0 when r=(kb/ka)r′ for r′∈D, an inaccurate location and shape of D must be retrieved through the map of F(ka,r). This is the theoretical reason why an inaccurate location and shape of the anomaly is retrieved when inaccurate values of εb and σb were applied. Further properties will be discussed in the simulation results.
Here, we present simulation results to support the result in Theorem 3.1. To this end, a circular array of N=16 antennas An is used to transmit and receive signals operated at f=1.24GHz. The location of the antennas was set to
an=0.09m(cosθn,sinθn),θn=2π(n−1)N |
and the search domain Ω was selected as a square region (−0.1m,0.1m)×(−0.1m,0.1m) with (εb,σb)=(20ε0,0.2S/m). Here, ε0=8.854×10−12F/m is the vacuum permittivity. Correspondingly, the exact value of the background wavenumber is kb=116.5273+8.4020i. For anomalies, we selected two small balls D1 and D2 with centers r1=(0.01m,0.03m) and r2=(−0.04m,−0.02m), same radii α=0.01m and material properties (ε1,σ1)=(55ε0,1.2S/m) and (ε2,σ2)=(45ε0,1.0S/m). With these settings, the measurement data Sscat(n,m) of (2.6) and the incident field data of (2.8) were generated by CST STUDIO SUITE.
Example 4.1. (Only exact value of εb is unknown) First, we consider the case where only the exact value of εb is unknown. Instead, of the application of εb, we applied alternative values εa and corresponding wavenumber
ka=ω√μb(εa−iσbω). |
Note that we already assumed that ωεb≫σb. Thus, if the condition ωεa≫σb is satisfied, the identified location becomes
r=(kbka)r′=√ωεb−iσbωεa−iσbr′≈√εbεar′for eachr′∈D1. | (4.1) |
Hence, identified anomalies will be concentrated at the origin and their retrieved sizes will be smaller than the true one when εa>εb. Otherwise, identified anomalies will be far from the origin and their retrieved sizes will be larger than the true one when εa<εb. See Figure 1 for a related illustration.
Figure 2 shows maps of F(ka,r) with various selections of εa in the presence of D1. As we discussed above, the location of the retrieved anomaly gets closer to the origin and it becomes smaller as the value of εa increases (here, εa=3εb,10εb). Otherwise, as the value of εa decreases (here, εa=0.5εb), the identified location becomes far from the origin and the size becomes larger. If εa=0.01εb i.e., the value of εa is very small, it is difficult to distinguish the D1 and artifacts.
Notice that since
Js(|kbr′−kar|)=Js(|ka||(kbka)r′−r|):=Js(|ka||r″−r|), |
due to the oscillating property of the Bessel function, several artifacts will be included in the map of F(ka,r) if |ka| is large, i.e., the value of εa is large enough compared to εb. In contrast, the map of F(ka,r) will contain no artifacts but the imaging result will be blurred if the value of εa is small enough compared to εb. This is the reason why several artifacts are included in the map of F(ka,r) when εa=10εb and why the obtained image is blurred when εa=0.1εb and εa=0.01εb. We can observe the same phenomenon in the presence of multiple anomalies D1 and D2, as shown in the Figure 3.
Example 4.2. (Only exact value of σb is unknown) Next, we consider the case where only the exact value of σb is unknown and apply an alternative one σa such that
ka=ω√μb(εb−iσaω). |
Same as the Example 4.1., if σa satisfies ωεa≫σa, the identified location becomes
r=(kbka)r′=√ωεb−iσbωεb−iσar′≈√εbεbr′=r′for eachr′∈D1. | (4.2) |
Hence, it will be possible to retrieve almost accurate shapes and locations of anomalies when σa is sufficiently small.
Figure 4 shows maps of F(ka,r) with various selections of σa in the presence of D1. In contrast to the results in Example 4.1., almost the exact location and shape of D1 were retrieved if σa<3σb, i.e., when σa was sufficiently small. The location of the retrieved anomaly gets closer to the origin and its size becomes smaller as the value of εa increases (here, εa=3εb,10εb). Unfortunately, it is very difficult to recognize D1 due to the appearance of a huge artifact with a large magnitude if σa is not small σa=10σb=2S/m.
We can observe the same phenomenon in the presence of multiple anomalies D1 and D2, as shown in the Figure 5, and conclude that it will be possible to retrieve the accurate shape and location of anomalies by choosing a very small (close to zero) value of σa when the exact value of background permittivity is known.
Example 4.3. (Identification of circle and rectangular shaped anomalies) Here, we consider the imaging of anomalies with different shapes. To this end, we applied f=1.0GHz, used N=36 antennas An, and selected D1 as a ball of Examples 1 and 2 except (ε1,σ1)=(45ε0,1.0S/m), and D2 as a square with vertices (−0.05,−0.03), (−0.03,−0.03), (−0.03,−0.01) and (−0.05,−0.01) with (ε2,σ2)=(45ε0,1.0S/m). With this configuration, the scattering parameter data were generated by using the FEKO (Feldberechnung für Körper mit beliebiger Oberfläche).
Figure 6 shows maps of F(ka,r) when εa≠εb and σa=σb. Similar to the results in Example 4.1., identified anomalies are concentrated at the origin and their retrieved sizes are smaller than the true one when εa>εb. Moreover, identified anomalies located far from the origin and their retrieved sizes are larger than the true one when εa<εb. However, opposite to the results in Examples 1 and 2, it is hard to recognize the shape of D1 and D2 due to the appearance of several artifacts in the neighborhood of anomalies.
Figure 7 shows maps of F(ka,r) when εa=εb and σa≠σb. Similar to the results in Example 4.2., we can examine that it is possible to recognize the outline shape of anomalies by choosing a very small value of σa. However, exact shape of anomalies cannot be retrieved still.
The structure of the imaging function of SM for retrieving small anomalies from scattering matrix is revealed when complete information of the background medium is not available. On the basis of its relationship with the infinite series of Bessel function of the first kind, we have theoretically confirmed why the accurate shape and location of anomalies cannot be retrieved.
The main subject of this paper is the imaging of small anomaly in two-dimensional microwave imaging. An extension to multiple, small anomalies will be carried out in forthcoming work. Moreover, the development of an effective algorithm for retrieving the exact value of background wavenumber will be an interesting research subject. Finally, we expect that the methodology presented in this paper could be applied to real-world microwave imaging with inhomogeneous background.
The author declares he has not used Artificial Intelligence (AI) tools in the creation of this article.
The author would like to acknowledge anonymous reviewers for their comments that help to increase the quality of the paper. The author is also grateful to Sangwoo Kang, Kwang-Jae Lee and Seong-Ho Son for helping in generating scattering parameter data. This work was supported by the research program of the Kookmin University.
The author declares no conflicts of interest regarding the publication of this paper.
[1] |
S. Nadeem, W. Fuzhang, F. M. Alharbi, F. Sajid, N. Abbas, A. S. El-Shafay, et al., Numerical computations for Buongiorno nano fluid model on the boundary layer flow of viscoelastic fluid towards a nonlinear stretching sheet, Alex. Eng. J., 61 (2022), 1769–1778. https://doi.org/10.1016/j.aej.2021.11.013 doi: 10.1016/j.aej.2021.11.013
![]() |
[2] |
M. Qayyum, S. Afzal, E. Ahmad, M. B. Riaz, Fractional modeling and analysis of unsteady squeezing flow of Casson nanofluid via extended He-Laplace algorithm in Liouville-Caputo sense, Alex. Eng. J., 73 (2023), 579–591. https://doi.org/10.1016/j.aej.2023.05.010 doi: 10.1016/j.aej.2023.05.010
![]() |
[3] |
A. E. Aboanber, A. A. Nahla, A. M. El-Mhlawy, O. Maher, An efficient exponential representation for solving the two-energy group point telegraph kinetics model, Ann. Nucl. Energy, 166 (2022), 108698. https://doi.org/10.1016/j.anucene.2021.108698 doi: 10.1016/j.anucene.2021.108698
![]() |
[4] |
C. Villa, A. Gerisch, M. A. J. Chaplain, A novel nonlocal partial differential equation model of endothelial progenitor cell cluster formation during the early stages of vasculogenesis, J. Theor. Biol., 534 (2022), 110963. https://doi.org/10.1016/j.jtbi.2021.110963 doi: 10.1016/j.jtbi.2021.110963
![]() |
[5] |
I. Ahmad, H. Ahmad, P. Thounthong, Y. Chu, C. Cesarano, Solution of multi-term time-fractional PDE models arising in mathematical biology and physics by local meshless method, Symmetry, 12 (2020), 1195. https://doi.org/10.1016/j.heliyon.2023.e16522 doi: 10.1016/j.heliyon.2023.e16522
![]() |
[6] |
O. D. Adeyemo, C. M. Khalique, Lie group classification of generalized variable coefficient Korteweg-de Vries equation with dual power-law nonlinearities with linear damping and dispersion in quantum field theory, Symmetry, 14 (2022), 83. https://doi.org/10.3390/sym14010083 doi: 10.3390/sym14010083
![]() |
[7] |
S. Afzal, M. Qayyum, M. B. Riaz, A. Wojciechowski, Modeling and simulation of blood flow under the influence of radioactive materials having slip with MHD and nonlinear mixed convection, Alex. Eng. J., 69 (2023), 9–24. https://doi.org/10.1016/j.aej.2023.01.013 doi: 10.1016/j.aej.2023.01.013
![]() |
[8] |
L. Guo, H. Wu, T. Zhou, Normalizing field flows: Solving forward and inverse stochastic differential equations using physics-informed flow models, J. Comput. Phys., 461 (2022), 111202. https://doi.org/10.1016/j.jcp.2022.111202 doi: 10.1016/j.jcp.2022.111202
![]() |
[9] |
S. P. Joseph, New traveling wave exact solutions to the coupled Klein-Gordon system of equations, PDE Appl. Math., 5 (2022), 100208. https://doi.org/10.1016/j.padiff.2021.100208 doi: 10.1016/j.padiff.2021.100208
![]() |
[10] |
M. Farman, A. Akgül, S. Askar, T. Botmart, A. Ahmad, H. Ahmad, Modeling and analysis of fractional order Zika model, AIMS Math., 7 (2022), 3912–3938. https://doi.org/10.3934/math.2022216 doi: 10.3934/math.2022216
![]() |
[11] | J. Liouville, Mémoire sur quelques questions de géométrie et de mécanique, et sur un nouveau genre de calcul pour résoudre ces questions, J. Éc. Polytech. Math., 1832. |
[12] | K. S. Miller, B. Ross, An Introduction to the Fractional Calculus and Fractional Differential Equations, Hoboken: Wiley, 1993. |
[13] | B. Riemann, Versuch einer allgemeinen auffassung der integration und differentiation, Gesammelte Werke, 62 (1876), 1876. |
[14] | M. Caputo, Elasticita e Dissipazione, Bologna: Zanichelli, 1969. |
[15] | H. Ahmad, N. Alam, M. Omri, New computational results for a prototype of an excitable system, Results Phys., 28 (2021). https://doi.org/10.1016/j.rinp.2021.104666 |
[16] |
M. Qayyum, E. Ahmad, S. Afzal, T. Sajid, W. Jamshed, A. Musa, et al., Fractional analysis of unsteady squeezing flow of casson fluid via homotopy perturbation method, Sci. Rep., 12 (2022), 18406. https://doi.org/10.1038/s41598-022-23239-0 doi: 10.1038/s41598-022-23239-0
![]() |
[17] |
K. S. Nisar, K. Logeswari, V. Vijayaraj, H. M. Baskonus, C. Ravichandran, Fractional order modeling the gemini virus in capsicum annuum with optimal control, Fractal Fract., 6 (2022), 61. https://doi.org/10.3390/fractalfract6020061 doi: 10.3390/fractalfract6020061
![]() |
[18] |
A. Yusuf, S. Qureshi, U. T. Mustapha, S. S. Musa, T. A. Sulaiman, Fractional modeling for improving scholastic performance of students with optimal control, Int. J. Appl. Comput. Math., 8 (2022), 37. https://doi.org/10.1007/s40819-021-01177-1 doi: 10.1007/s40819-021-01177-1
![]() |
[19] |
H. Hassani, J. A. Tenreiro Machado, Z. Avazzadeh, E. Naraghirad, S. Mehrabi, Optimal solution of the fractional-order smoking model and its public health implications, Nonlinear Dynam., 108 (2022), 2815–2831. https://doi.org/10.1007/s11071-022-07343-4 doi: 10.1007/s11071-022-07343-4
![]() |
[20] |
M. Qayyum, E. Ahmad, S. T. Saeed, H. Ahmad, S. Askar, Homotopy perturbation method-based soliton solutions of the time-fractional (2+1)-dimensional wu-zhang system describing long dispersive gravity water waves in the ocean, Front. Phys., 11 (2023), 1178154. https://doi.org/10.3389/fphy.2023.1178154 doi: 10.3389/fphy.2023.1178154
![]() |
[21] |
C. Wang, X. Zhou, X. Shi, Y. Jin, Variable fractional order sliding mode control for seismic vibration suppression of uncertain building structure, J. Vib. Eng. Tech., 10 (2021), 299–312. https://doi.org/10.1007/s42417-021-00377-9 doi: 10.1007/s42417-021-00377-9
![]() |
[22] |
I. M. Batiha, S. A. Njadat, R. M. Batyha, A. Zraiqat, A. Dababneh, S. Momani, Design fractional-order PID controllers for single-joint robot arm model, Int. J. Adv. Soft Comput. Appl., 14 (2022), 97–114. https://doi.org/10.15849/ijasca.220720.07 doi: 10.15849/ijasca.220720.07
![]() |
[23] |
M. H. Derakhshan, Existence, uniqueness, Ulam-Hyers stability and numerical simulation of solutions for variable order fractional differential equations in fluid mechanics, J. Appl. Math. Comput., 68 (2021), 403–429. https://doi.org/10.1007/s12190-021-01537-6 doi: 10.1007/s12190-021-01537-6
![]() |
[24] |
A. Cardone, D. Conte, R. D'Ambrosio, B. Paternoster, Multivalue collocation methods for ordinary and fractional differential equations, Mathematics, 10 (2022), 185. https://doi.org/10.3390/math10020185 doi: 10.3390/math10020185
![]() |
[25] |
N. A. Shah, A. Wakif, E. R. El-Zahar, T. Thumma, S. J. Yook, Heat transfers thermodynamic activity of a second-grade ternary nanofluid flow over a vertical plate with Atangana-Baleanu time-fractional integral, Alex. Eng. J., 61 (2022), 10045–10053. https://doi.org/10.1016/j.aej.2022.03.048 doi: 10.1016/j.aej.2022.03.048
![]() |
[26] |
N. P. Dong, H. V. Long, N. L. Giang, The fuzzy fractional SIQR model of computer virus propagation in wireless sensor network using Caputo Atangana–Baleanu derivatives, Fuzzy Sets Syst., 429 (2022), 28–59. https://doi.org/10.1016/j.fss.2021.04.012 doi: 10.1016/j.fss.2021.04.012
![]() |
[27] |
A. Din, F. M. Khan, Z. U. Khan, A. Yusuf, T. Munir, The mathematical study of climate change model under nonlocal fractional derivative, PDE Appl. Math., 5 (2022), 100204. https://doi.org/10.1016/j.padiff.2021.100204 doi: 10.1016/j.padiff.2021.100204
![]() |
[28] |
Y. Gurefe, Y. Pandir, T. Akturk, Analysis of exact solutions of a mathematical model by new function method, Cumhuriyet Sci. J., 43 (2022), 703–707. https://doi.org/10.17776/csj.1083033 doi: 10.17776/csj.1083033
![]() |
[29] |
M. R. Ahamed Fahim, P. R. Kundu, M. E. Islam, M. A. Akbar, M. S. Osman, Wave profile analysis of a couple of (3+1)-dimensional nonlinear evolution equations by sine-Gordon expansion approach, J. Ocean Eng. Sci., 7 (2022), 272–279. https://doi.org/10.1016/j.joes.2021.08.009 doi: 10.1016/j.joes.2021.08.009
![]() |
[30] |
B. Ghanbari, Employing Hirota's bilinear form to find novel lump waves solutions to an important nonlinear model in fluid mechanics, Results Phys., 29 (2021), 104689. https://doi.org/10.1016/j.rinp.2021.104689 doi: 10.1016/j.rinp.2021.104689
![]() |
[31] | W. Razzaq, A. Zafar, H. M. Ahmed, W. B. Rabie, Construction solitons for fractional nonlinear Schrodinger equation with β-time derivative by the new sub-equation method, J. Ocean Eng. Sci., in press, 2022. https://doi.org/10.1016/j.joes.2022.06.013 |
[32] |
G. Akram, M. Sadaf, S. Arshed, F. Sameen, Bright, dark, kink, singular and periodic soliton solutions of Lakshmanan-Porsezian-Daniel model by generalized projective riccati equations method, Optik, 241 (2021), 167051. https://doi.org/10.1016/j.ijleo.2021.167051 doi: 10.1016/j.ijleo.2021.167051
![]() |
[33] |
S. Liao, Homotopy analysis method: A new analytical technique for nonlinear problems, Commun. Nonlinear Sci. Numer. Simul., 2 (1997), 95–100. https://doi.org/10.1016/s1007-5704(97)90047-2 doi: 10.1016/s1007-5704(97)90047-2
![]() |
[34] |
S. Afzal, M. Qayyum, G. Chambashi, Heat and mass transfer with entropy optimization in hybrid nanofluid using heat source and velocity slip: a hamilton–crosser approach, Sci. Rep., 13 (2023), 12392. https://doi.org/10.1038/s41598-023-39176-5 doi: 10.1038/s41598-023-39176-5
![]() |
[35] |
T. Hayat, K. Muhammad, S. Momani, Melting heat and viscous dissipation in flow of hybrid nanomaterial: a numerical study via finite difference method, J. Therm. Anal. Calorime., 147 (2021), 6393–6401. https://doi.org/10.1007/s10973-021-10944-7 doi: 10.1007/s10973-021-10944-7
![]() |
[36] |
H. Ahmad, M. N. Khan, I. Ahmad, M. Omri, M. F. Alotaibi, A meshless method for numerical solutions of linear and nonlinear time-fractional Black-Scholes models, AIMS Math., 8 (2023), 19677–19698. https://doi.org/10.3934/math.20231003 doi: 10.3934/math.20231003
![]() |
[37] |
J. H. He, Homotopy perturbation technique, Comput. Meth. Appl. Mech. Eng., 178 (1999), 257–262. https://doi.org/10.1016/s0045-7825(99)00018-3 doi: 10.1016/s0045-7825(99)00018-3
![]() |
[38] |
R. Amin, K. Shah, H. Ahmad, A. H. Ganie, A. Abdel-Aty, T. Botmart, Haar wavelet method for solution of variable order linear fractional integro-differential equations, AIMS Math., 7 (2022), 5431–5443. https://doi.org/10.3934/math.2022301 doi: 10.3934/math.2022301
![]() |
[39] |
G. Singh, I. Singh, New laplace variational iterative technique to solve twodimensional Schr¨odinger equation, Mater. Today Proc., 62 (2022), 3995–4000. https://doi.org/10.1016/j.matpr.2022.04.585 doi: 10.1016/j.matpr.2022.04.585
![]() |
[40] |
M. Croci, G. R. de Souza, Mixed-precision explicit stabilized Runge-Kutta methods for single- and multi-scale differential equations, J. Comput. Phys., 464 (2022), 111349. https://doi.org/10.1016/j.jcp.2022.111349 doi: 10.1016/j.jcp.2022.111349
![]() |
[41] |
M. Aslam, M. Farman, H. Ahmad, T. N. Gia, A. Ahmad, S. Askar, Fractal fractional derivative on chemistry kinetics hires problem, AIMS Math., 7 (2022), 1155–1184. https://doi.org/10.3934/math.2022068 doi: 10.3934/math.2022068
![]() |
[42] |
J. H. He, M. L. Jiao, K. A. Gepreel, Y. Khan, Homotopy perturbation method for strongly nonlinear oscillators, Math. Comput. Simul., 204 (2023), 243–258. https://doi.org/10.1016/j.matcom.2022.08.005 doi: 10.1016/j.matcom.2022.08.005
![]() |
[43] |
M. Qayyum, E. Ahmad, M. B. Riaz, J. Awrejcewicz, Improved soliton solutions of generalized fifth order time-fractional KdV models: Laplace transform with homotopy perturbation algorithm, Universe, 8 (2022), 563. https://doi.org/10.3390/universe8110563 doi: 10.3390/universe8110563
![]() |
[44] | Y. Pandir, T. A¨gir, Genisletilmis deneme denklemi yöntemi ile k¨ubik lineer olmayan Schrödinger denkleminin yeni tam ccöz¨umleri, Afyon Kocatepe Uni. J. Sci. Eng., 20 (2020), 582–588. |
[45] |
K. J. Wang, G. D. Wang, Variational theory and new abundant solutions to the (1+2)-dimensional chiral nonlinear Schrödinger equation in optics, Phys. Letters A, 412 (2021), 127588. https://doi.org/10.1016/j.physleta.2021.127588 doi: 10.1016/j.physleta.2021.127588
![]() |
[46] |
M. Al-Smadi, O. A. Arqub, S. Momani, Numerical computations of coupled fractional resonant Schrödinger equations arising in quantum mechanics under conformable fractional derivative sense, Phys. Scripta, 95 (2020), 075218. https://doi.org/10.1088/1402-4896/ab96e0 doi: 10.1088/1402-4896/ab96e0
![]() |
[47] |
S. F. Tian, X. F. Wang, T. T. Zhang, W. H. Qiu, Stability analysis, solitary wave and explicit power series solutions of a (2 + 1)-dimensional nonlinear Schrödinger equation in a multicomponent plasma, Int. J. Numer. Meth. Heat Fluid Flow, 31 (2021), 1732–1748. https://doi.org/10.1108/hff-08-2020-0517 doi: 10.1108/hff-08-2020-0517
![]() |
[48] |
D. F. Li, J. L. Wang, J. W. Zhang, Unconditionally convergent l1-galerkin FEMs for nonlinear time-fractional Schrödinger equations, SIAM J. Sci. Comput., 39 (2017), A3067–A3088. https://doi.org/10.1137/16m1105700 doi: 10.1137/16m1105700
![]() |
[49] |
W. Q. Yuan, C. J. Zhang, D. F. Li, Linearized fast time-stepping schemes for time–space fractional Schrödinger equations, Phys. D Nonlinear Phenomena, 454 (2023), 133865. https://doi.org/10.1016/j.physd.2023.133865 doi: 10.1016/j.physd.2023.133865
![]() |
[50] |
K. Hosseini, E. Hincal, S. Salahshour, M. Mirzazadeh, K. Dehingia, B. J. Nath, On the dynamics of soliton waves in a generalized nonlinear Schrödinger equation, Optik, 272 (2023), 170215. https://doi.org/10.1016/j.ijleo.2022.170215 doi: 10.1016/j.ijleo.2022.170215
![]() |
[51] |
N. A. Kudryashov, Method for finding optical solitons of generalized nonlinear Schrödinger equations, Optik, 261 (2022), 169163. https://doi.org/10.1016/j.ijleo.2022.169163 doi: 10.1016/j.ijleo.2022.169163
![]() |
[52] |
W. Q. Yuan, D. F. Li, C. J. Zhang, Linearized transformed l1 galerkin FEMs with unconditional convergence for nonlinear time fractional Schrödinger equations, Numer. Math. Theory Meth. Appl., 16 (2023), 348–369. https://doi.org/10.4208/nmtma.oa-2022-0087 doi: 10.4208/nmtma.oa-2022-0087
![]() |
[53] |
T. Y. Han, Z. Li, X. Zhang, Bifurcation and new exact traveling wave solutions to time-space coupled fractional nonlinear Schrödinger equation, Phys. Letters A, 395 (2021), 127217. https://doi.org/10.1016/j.physleta.2021.127217 doi: 10.1016/j.physleta.2021.127217
![]() |
[54] |
P. F. Dai, Q. B. Wu, An efficient block Gauss-Seidel iteration method for the space fractional coupled nonlinear Schrödinger equations, Appl. Math. Letters, 117 (2021), 107116. https://doi.org/10.1016/j.aml.2021.107116 doi: 10.1016/j.aml.2021.107116
![]() |
[55] |
K. S. Nisar, S. Ahmad, A. Ullah, K. Shah, H. Alrabaiah, M. Arfan, Mathematical analysis of SIRD model of COVID-19 with Caputo fractional derivative based on real data, Results Phys., 21 (2021), 103772. https://doi.org/10.1016/j.rinp.2020.103772 doi: 10.1016/j.rinp.2020.103772
![]() |
[56] |
T. Bakkyaraj, Lie symmetry analysis of system of nonlinear fractional partial differential equations with Caputo fractional derivative, Eur. Phys. J. Plus, 135 (2020), 126. https://doi.org/10.1140/epjp/s13360-020-00170-9 doi: 10.1140/epjp/s13360-020-00170-9
![]() |
[57] |
N. H. Tuan, H. Mohammadi, S. Rezapour, A mathematical model for COVID-19 transmission by using the Caputo fractional derivative, Chaos Solitons Fract., 140 (2020), 110107. https://doi.org/10.1016/j.chaos.2020.110107 doi: 10.1016/j.chaos.2020.110107
![]() |
1. | Janghoon Jeong, Seong-Ho Son, Localization of small moving objects using only total fields in microwave imaging, 2024, 61, 22113797, 107777, 10.1016/j.rinp.2024.107777 | |
2. | Junyong Eom, Won-Kwang Park, Real-time detection of small objects in transverse electric polarization: Evaluations on synthetic and experimental datasets, 2024, 9, 2473-6988, 22665, 10.3934/math.20241104 | |
3. | Seong-Ho Son, Kwang-Jae Lee, Won-Kwang Park, Real-time tracking of moving objects from scattering matrix in real-world microwave imaging, 2024, 9, 2473-6988, 13570, 10.3934/math.2024662 | |
4. | Won-Kwang Park, On the application of subspace migration from scattering matrix with constant-valued diagonal elements in microwave imaging, 2024, 9, 2473-6988, 21356, 10.3934/math.20241037 | |
5. | Janghoon Jeong, Jang-Moon Jo, Soeng-Ho Son, 2024, Microwave Imaging Method for Object Localization without Background Measurements, 979-8-3503-7581-7, 108, 10.1109/ICITEE62483.2024.10808351 |