Next Article in Journal
Hyperparameter Optimization of Bayesian Neural Network Using Bayesian Optimization and Intelligent Feature Engineering for Load Forecasting
Previous Article in Journal
Design and Development of a Battery State of Health Estimation Model for Efficient Battery Monitoring Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Narrow-Bandpass One-Step Leapfrog Hybrid Implicit-Explicit Algorithm with Convolutional Boundary Condition for Its Applications in Sensors

1
School of Electronic and Information Engineering, Beihang University, Beijing 100191, China
2
Beijing Key Laboratory of Microwave Sensing and Security Applications, Beihang University, Beijing 100191, China
3
Shenzhen Institute of Beihang University, Shenzhen 518000, China
4
Xi’an Institute of Space Radio Technology, Xi’an 710100, China
5
School of Electronic and Information Engineering, Nanjing University of Information Science and Technology, Nanjing 210096, China
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(12), 4445; https://doi.org/10.3390/s22124445
Submission received: 25 April 2022 / Revised: 7 June 2022 / Accepted: 9 June 2022 / Published: 12 June 2022
(This article belongs to the Section Intelligent Sensors)

Abstract

:
A large number of sensors work in the narrow bandpass circumstance. Meanwhile, some of them hold fine details merely along one and two dimensions. In order to efficiently simulate these sensors and devices, the one-step leapfrog hybrid implicit-explicit (HIE) algorithm with the complex envelope (CE) method and absorbing boundary condition is proposed in the narrow bandpass circumstance. To be more precise, absorbing boundary condition is implemented by the higher order convolutional perfectly matched layer (CPML) formulation to further enhance the absorption during the entire simulation. Numerical examples and their experiments are carried out to further illustrate the effectiveness of the proposed algorithm. The results show considerable agreement with the experiment and theory resolution. The relationship between the time step and mesh size can break the Courant–Friedrichs–Levy condition which indicates the physical size/selection mesh size. Such a condition indicates that the proposed algorithm behaviors are considerably accurate due to the rational choice in discretized mesh. It also shows decrement in simulation duration and memory consumption compared with the other algorithms. In addition, absorption performance can be improved by employing the proposed higher order CPML algorithm during the whole simulation.

1. Introduction

With the development of the electronic industry, sensors have become more important than ever before. Large numbers of sensors work in narrow bandpass circumstances or individual frequency [1,2]. Although the finite-difference time-domain (FDTD) algorithm has been extensively employed in many applications, sensor simulation by the FDTD algorithm is severely limited by the inefficient bandpass calculation. Such conditions also restrict the development of the FDTD algorithm [3]. The reason is that the FDTD algorithm is based on the lowpass-sampling-theorem which can achieve considerable performance merely in broadband simulation [4]. The time step becomes small, resulting in an unacceptable simulation duration in bandpass problems. In order to efficiently simulate bandpass problems, investigation of the alternative FDTD algorithm has become urgent. Complex envelope (CE) method which is based on the bandpass-sampling-theorem shows considerable potential in bandpass problems [5]. In the CE method, the mesh size can be chosen according to the bandwidth rather than the maximum frequency, resulting in significantl improvement, both from the aspect of efficiency and accuracy. The CE method was firstly employed into the conventional FDTD algorithm [6]. Then, it was introduced into the unconditionally stable algorithms [7,8,9,10]. Through results, it can be concluded that these algorithms can receive better behavior compared with the lowpass-sampling-theorem based algorithms in bandpass simulation.
With the development of sensors, structures have become more complex than ever before. For the fine details simulation, mesh size must be fine enough to satisfy the calculation accuracy. Due to the existence of the Courant–Friedrichs–Levy (CFL) condition, the relationship between time step and mesh size is established [11]. Thus, time step must also be fine enough to maintain the stability of the algorithm. Such conditions also result in an extremely long simulation duration. In order to alleviate such conditions, unconditionally stable algorithms are proposed including the alternating direction implicit (ADI), locally one-dimension (LOD), split-step (SS) procedures and others [12,13,14,15]. It has been testified that the unconditionally stable algorithms are merely efficient in fine details with all directions. When simulating fine details along one and two dimensions, more than six matrices should be calculated at each time step resulting in increments of memory consumption and simulation duration. Although these algorithms can overcome the CFL condition, the effectiveness of the algorithm decreases significantly. In order to alleviate such conditions, a hybrid implicit-explicit (HIE) procedure is proposed [16]. Through explicitly updating components in a single direction and implicitly updating the others, the HIE procedure can obtain considerable performance in low-dimensions [17,18]. However, the HIE procedure is still a split-field updating procedure which increases the memory and simulation duration [19]. To alleviate such conditions, the one-step leapfrog scheme is employed into the HIE procedure [20].
Because the computer cannot solve infinite computational domains, an adequate boundary condition must be employed at the boundaries [11]. To simulate infinite computational domain in finite space, perfectly matched layer (PML) is regarded as one of the most powerful and efficient absorbing boundary conditions [21]. The original PML formulation is a split-field scheme which shows limitations in absorption and resources [22]. Unsplit-field formulation is proposed to alleviate such conditions including convolutional PML (CPML), stretched coordinate PML (SC-PML) and complex-frequency-shifted PML (CFS-PML) schemes [23,24,25]. Although they hold advantages in reducing late-time reflections and absorbing low-frequency evanescent waves, reflection at the low frequency still needs to be enhanced [26]. Higher order PML formulation is proposed to not only enhance the absorption in the time domain but also absorb low-frequency propagation waves in the frequency domain [27,28]. The original higher order PML formulation must introduce six auxiliary variables during calculation which occupies large number of computational resources [29]. To alleviate such conditions, the higher order PML formulation with four auxiliary variables is proposed [30]. In addition, existing higher order formulation is mainly based on the CFS-PML scheme which shows limitation in complex media. The higher order CPML is proposed to overcome these drawbacks [31,32]. It has been testified that compared with the other schemes, it can significantly decrease the computational complexity with considerable absorption in complex structures and media.
In order to simulate complex sensors and components whose fine details merely exist along one and two dimensions, several algorithms have been proposed which still have their own limitations. The original HIE procedures with PML (NPML) and CFS-PML are proposed in two dimensions which cannot be employed into the practical employment [33,34]. Recently, the HIE procedure and its leapfrog scheme with higher order PML are proposed in three dimensions [35,36,37]. However, they can merely be employed in vacuum and magnetized plasma. The formulation for dielectric material still needs to be investigated. From another aspect, all of the existing absorbing boundary conditions for the HIE procedure and leapfrog scheme are based on the NPML and CFS-PML formulation [35,36,37,38,39]. These algorithms must be modified according to the different computational domains due to the media-dependent formulation. Computation complexity increases significantly due to the employment of higher order formulation. The reduction in computation complexity also needs to be further investigated.
In order to efficiently simulate sensors with fine details in narrow bandpass circumstances, the one-step leapfrog HIE procedure is proposed with the CE method and higher order CPML formulation. The proposed algorithm takes advantages of leapfrog HIE procedure, CE method and higher order CPML formulation, in terms of the calculating of fine details along one and two dimensions, simulating bandpass problems and enhancing absorption. Numerical examples and experiments are carried out for the illustration of effectiveness. It can be concluded that results simulated by employing the proposed algorithm show considerable agreement with the experiments and theory resolution. Compared to the other algorithms, the proposed algorithm also shows significant improvement in accuracy, efficiency and absorption.

2. Formulation

Because massive components exist in Maxwell’s equations, assuming the fine details merely locate along x-direction, E y and H z components are selected as examples for demonstration. The other components can be obtained by employing the similar approach. The Maxwell’s equations inside the dielectric material can be given as follows:
j ω ε 0 E y + σ E y = S z 1 z H x S x 1 x H z
j ω μ 0 H y = S z 1 z E x S x 1 x E z
where S η ,   η = x , y , z is the stretched coordinate variable side PML regions with higher order formulation which can be defined as:
S η = ( κ η 1 + σ η 1 α η 1 + j ω ε 0 ) ( κ η 2 + σ η 2 α η 2 + j ω ε 0 )
where σ η n ,   n = 1 , 2 and α η n are supposed to be positive real, κ η n is supposed to be real. According to the CPML formulation, by employing the partial fraction expansion to the stretched coordinate variables, results can be converted by employing the inverse Laplace transform, one obtains:
S ¯ η ( t ) = ( δ ( t ) κ η 1 + ζ η 1 ( t ) ) ( δ ( t ) κ η 2 + ζ η 2 ( t ) )
where the coefficients can be given as ζ η n = r η n e v η n t u ( t ) , r η n = σ η n / ( κ η n 2 ε 0 ) , v η n = α η n / ε 0 + σ η n / ( κ η n ε 0 ) . In addition, δ ( t ) and u ( t ) represent the unit speed function and unit impulse function. It can be observed that the original Maxwell’s equations obey the lowpass sampling-theorem which shows performance decrement with the bandpass simulation. In order to alleviate such conditions, the CE method is introduced into the Maxwell’s equations [5]. According to the CE method, the original signal can be given as the complex envelope form, one obtains:
Φ = Re { Φ ^ e j ϖ t }
where Re { } represents the real number in the complex signal, Φ and Φ ^ represent the original signal and complex envelope signal, ϖ is the center frequency of the modulated signal. According to the CE method and CPML formulation, Maxwell’s equations can be given as:
( j ϖ + t ) E ^ y + σ E ^ y = S ¯ z ( t ) * z H ^ x S ¯ x ( t ) * x H ^ z
( j ϖ + t ) H ^ z = S ¯ x ( t ) * x E ^ z S ¯ z ( t ) * z E ^ x
In order to update the equations, Equation (2) is substituted into Equation (4a,b), one obtains:
( j ϖ + t ) E ^ y + σ E ^ y = F ^ y z 12 + κ z 1 1 F ^ y z 1 + κ z 2 1 F ^ y z 2 + κ z 1 1 κ z 2 1 z H ^ x F ^ y x 12 κ x 1 1 F ^ y x 1 κ x 2 1 F ^ y x 2 + κ x 1 1 κ x 2 1 x H ^ z
( j ϖ + t ) H ^ z = G ^ z x 12 + κ x 1 1 G ^ z x 1 + κ x 2 1 G ^ z x 2 + κ x 1 1 κ x 2 1 x E ^ y G ^ z y 12 κ y 1 1 G ^ z y 1 κ y 2 1 G ^ z y 2 κ y 1 1 κ y 2 1 y E ^ x
where F and G are the auxiliary variables. By rewriting the equations into the discretized domain for updating, the results can be given according to the sub-step procedure as:
The First Time Step:
E ^ y n + 1 / 2 p 3 x δ x H ^ z n + 1 / 2 + p 3 z δ z H ^ x n + 1 / 2 = k 1 c 1 E ^ y n + p 4 z F ^ y z 12 n + p 1 z F ^ y z 1 n + p 2 z F ^ y z 2 n p 4 x F ^ y x 12 n p 1 x F ^ y x 1 n p 2 x F ^ y x 2 n
H ^ z n + 1 / 2 p 3 x δ x E ^ y n + 1 / 2 = c 1 H ^ z n + p 4 x G ^ z x 12 n + p 1 x G ^ z x 1 n + p 2 x G ^ z x 2 n p 4 y G ^ z y 12 n p 1 y G ^ z y 1 n p 2 y G ^ z y 2 n p 3 y δ y E ^ x n
The Second Time Step:
E ^ y n + 1 = k 1 c 1 E ^ y n + 1 / 2 + p 4 z F ^ y z 12 n + 1 / 2 + p 1 z F ^ y z 1 n + 1 / 2 + p 2 z F ^ y z 2 n + 1 / 2 + p 3 z δ z H ^ x n + 1 / 2 p 4 x F ^ y x 12 n + 1 / 2 p 1 x F ^ y x 1 n + 1 / 2 p 2 x F ^ y x 2 n + 1 / 2 p 3 x δ x H ^ z n + 1 / 2
H ^ z n + 1 + p 3 y δ y E ^ x n + 1 = c 1 H ^ z n + 1 / 2 + p 4 x G ^ z x 12 n + 1 / 2 + p 1 x G ^ z x 1 n + 1 / 2 + p 2 x G ^ z x 2 n + 1 / 2 p 3 x δ x E ^ y n + 1 / 2 p 4 y G ^ z y 12 n + 1 / 2 p 1 y G ^ z y 1 n + 1 / 2 p 2 y G ^ z y 2 n + 1 / 2
where the coefficients can be given as follows:
k 1 = ε 0 / ( ε 0 + σ Δ t ) , c 1 = ( 2 j ϖ Δ t ) / ( 2 + j ϖ Δ t ) , c 2 = 2 Δ t / ( 2 + j ϖ Δ t ) , c 3 = κ η 1 1 κ η 2 1 + a η 2 κ η 1 1 + a η 1 κ η 2 1 + a η 1 a η 2 , a η n = [ σ η n exp ( v η n Δ t 1 ) ] / ( σ η n κ η n + α η n κ η n 2 ) , b η n = v η n Δ t , p 1 η = c 2 k 1 c 3 κ z 2 1 , p 2 η = c 2 k 1 c 3 κ z 2 1 , p 3 η = c 2 k 1 c 3 κ z 1 1 κ z 2 1 / Δ η and p 4 η = c 2 k 1 c 3 .
According to the leapfrog HIE procedure, the first time step can be given by substituting the magnetic components into electric components, one obtains:
( 1 p 3 x 2 δ 2 x ) E ^ y n + 1 / 2 = ( k 1 c 1 p 3 z 2 δ 2 z ) E ^ y n + p 4 z F ^ y z 12 n + p 1 z F ^ y z 1 n + p 2 z F ^ y z 2 n p 4 x F ^ y x 12 n p 1 x F ^ y x 1 n p 2 x F ^ y x 2 n + p 1 x p 3 x δ x G ^ z x 1 n + p 2 x p 3 x δ x G ^ z x 2 n p 4 y p 3 x δ x G ^ z y 12 n p 1 y p 3 x δ x G ^ z y 1 n p 2 y p 3 x δ x G ^ z y 2 n p 1 z p 3 z δ z G ^ x z 1 n p 2 z p 3 z δ z G ^ x z 2 n + p 4 y p 3 z δ z G ^ x y 12 n + p 1 y p 3 z δ z G ^ x y 1 n + p 2 y p 3 z δ z G ^ x y 2 n + c 1 p 3 x δ x H ^ z n + p 4 x p 3 x δ x G ^ z x 12 n p 3 y p 3 x δ x δ y E ^ x n + 1 / 2 c 1 p 3 z δ z H ^ x n p 4 z p 3 z δ z G ^ x z 12 n p 3 y p 3 z δ z δ y E ^ z n + 1 / 2
By rewriting the second time step to the previous step of n − 1, substituting magnetic components to electric components accordingly, one obtains:
( 1 p 3 z 2 δ 2 z ) E ^ y n = ( k 1 c 1 + p 3 x 2 δ 2 x ) E ^ y n 1 / 2 + p 4 z F ^ y z 12 n 1 / 2 + p 1 z F ^ y z 1 n 1 / 2 + p 2 z F ^ y z 2 n 1 / 2 p 4 x F ^ y x 12 n 1 / 2 p 1 x F ^ y x 1 n 1 / 2 p 2 x F ^ y x 2 n 1 / 2 + p 1 x p 3 x δ x G ^ z x 1 n 1 / 2 + p 2 x p 3 x δ x G ^ z x 2 n 1 / 2 p 4 y p 3 x δ x G ^ z y 12 n 1 / 2 p 1 y p 3 x δ x G ^ z y 1 n 1 / 2 p 2 y p 3 x δ x G ^ z y 2 n 1 / 2 + p 1 x p 3 z δ z G ^ z x 1 n 1 / 2 + p 2 x p 3 z δ z G ^ z x 2 n 1 / 2 p 4 y p 3 z δ z G ^ z y 12 n 1 / 2 p 1 y p 3 z δ z G ^ z y 1 n 1 / 2 p 2 y p 3 z δ z G ^ z y 2 n 1 / 2 + p 4 x p 3 x δ x G ^ z x 12 n 1 / 2 c 1 p 3 z δ z H ^ x n p 3 y p 3 z δ z δ y E ^ x n 1 / 2 + p 4 x p 3 z δ z G ^ z x 12 n 1 / 2 + c 1 p 3 x δ x H ^ z n p 3 y p 3 z δ z δ y E ^ z n 1 / 2
According to the leapfrog scheme, by adding Equations (7) and (8) together into a single equation, after some manipulations, one obtains:
( 1 p 3 x 2 δ 2 x ) E ^ y n + 1 / 2 = ( k 1 2 c 1 2 + k 1 c 1 p 3 x 2 δ 2 x ) E ^ y n 1 / 2 + c 1 p 3 x δ x H ^ z n c 1 p 3 z δ z H ^ x n + k 1 c 1 p 4 z 2 ( F ^ y z 12 n + F ^ y z 12 n 1 ) + k 1 c 1 p 1 z 2 ( F ^ y z 1 n + F ^ y z 1 n 1 ) + k 1 c 1 p 2 z 2 ( F ^ y z 2 n + F ^ y z 2 n 1 ) k 1 c 1 p 4 x 2 ( F ^ y x 12 n + F ^ y x 12 n 1 ) k 1 c 1 p 1 x 2 ( F ^ y x 1 n + F ^ y x 1 n 1 ) k 1 c 1 p 2 x 2 ( F ^ y x 2 n + F ^ y x 2 n 1 ) k 1 c 1 p 3 y p 3 z δ z δ y ( E ^ x n + 1 / 2 + E ^ x n 1 / 2 ) k 1 c 1 p 3 y p 3 z δ z δ y ( E ^ z n + 1 / 2 + E ^ z n 1 / 2 )
The auxiliary variables can be given as, for example,
F ^ x η 1 n + 1 / 2 = a η 2 m = 0 n 1 e v η 2 m Δ t η H ^ η ˜ n m
F ^ x η 2 n + 1 / 2 = a η 1 m = 0 n 1 e v η 1 m Δ t η H ^ η ˜ n m
F ^ x η 12 n + 1 / 2 = a η 1 a η 2 m = 0 n 1 e v η 1 m Δ t m = 0 n 1 e v η 2 m Δ t η H ^ η ˜ n m
where η ˜ denotes the complement of η , for example, η = y , and η ˜ = z . In order to solve the summation term in Equation (10a–c), the recursive convolution (RC) method is employed during the calculation. According to the RC method, the results of auxiliary variables can be given as:
F ^ x η 1 n + 1 / 2 = b η 2 F ^ x η 1 n 1 / 2 + a η 2 η H ^ η ˜ n
F ^ x η 2 n + 1 / 2 = b η 1 F ^ x η 2 n 1 / 2 + a η 1 η H ^ η ˜ n
F ^ x η 1 n + 1 / 2 = b η 1 b η 2 F ^ x η 1 n 1 / 2 + a η 1 a η 2 η H ^ η ˜ n + 1 / 2
According to the convexity of Maxwell’s equations, the other auxiliary variables and components can be solved by employing a similar approach. By employing the above-mentioned equations, leapfrog HIE procedure with CPML formulation and CE method can be implemented. The entire update procedure of the described is as follows:
(1)
Explicitly update E x n + 1 / 2 , E z n + 1 / 2 , in the half-integer time step;
(2)
Explicitly update H x n , H z n in the integer time step;
(3)
Implicitly update E y n + 1 / 2 and H y n + 1 / 2 along the directions of fine details;
(4)
Explicitly update auxiliary variables inside PML regions.
The block diagram of the entire update procedure is shown in Figure 1 which corresponds to the above-mentioned procedure. In order to demonstrate the simulation duration and memory consumption of different PML algorithms, a number of operators are employed during comparison. During the algorithm implementation, a calculation of multiplication/division and addition/subtraction determine the efficiency and memory in theory. Multiplication/division operators and implicit equations occupy much more resources compared with addition/subtraction operators and explicit equations, respectively. Thus, it is important to analyze the influence of multiplication/division and addition/subtraction numbers on efficiency and resources. The compare different PML algorithms, FDTD-PML in [40], HIE-PML in [33], LHIE-CPML in [34], CE-HIE-HPML in [41], CE-LHIE-HPML and FDTD-HPML in [42] with explicit and implicit equations are employed with results that are shown in Table 1.
As can be concluded from Table 1, the explicit algorithms show the most considerable performance is due to the non-calculation of matrices. Due to the introduction of the implicit scheme, a number of multiplication/division and addition/subtraction operators show significant improvement, resulting in increments of simulation duration and consumption memory. Meanwhile, due to the employment of the higher order formulation, algorithms also show significant increments in efficiency and resources at the same time. To be more precise, the leapfrog scheme and CPML formulation show a decrement in computational complexity. Such conditions result in a memory and simulation duration decrement by CPML and leapfrog schemes compared to the CFS-PML and its original schemes, respectively. Due to the introduction of the CE method, several coefficients and operators should be calculated at each time step resulting in the increment of memory. In conclusion, the proposed algorithm can reduce the complexity of the algorithm compared with the other PML formulations based on the HIE procedure. However, in the bandpass simulation, the time step can be enlarged according to the percent bandwidth of the source. Owing to such conditions, the CE method can significantly improve the entire performance during the simulation.

3. Numerical Results and Experiments

The effectiveness of the proposed scheme can be demonstrated through employing the numerical example and experiment including a filter for the sensors system with experiments and remote sensing problems with theory resolution. The computer for the code holds the parameter of i7-11700k and DDR4 128 GB. The program is implemented by Fortran under Visual Studio. The entire computation duration can be described as the following procedure: (1) Import geometry structures with STL files; (2) Set condition for calculation including time step, mesh size and more; (3) Start Calculation; (4) Export results for TXT file.

3.1. Micro-Strip Filter for Sensors System

Within microwave components, micro-strip components can be regarded one of the most powerful and widely spread. Among some of them, an extremely thin metal patch is located on the surface of the device. Thus, by employing the unconditionally stable algorithms and conventional explicit algorithms directly, uniform extremely fine mesh must be employed within the entire computational domain which results in significant increments of memory and duration. In order to alleviate such conditions, the HIE procedure is introduced for the fine details along one and two dimensions. Here, a micro-strip filter for sensors system is employed for the illustration. The sketch picture of the structure is shown as Figure 2. Through the establishment of a rectangular coordinate system, a complex and oblique patch can be expressed and located by the specific coordinate location, as shown in Table 2. Along the vertex coordinates of the patch, the entire structure can be established in the CAD mechanical drawing. The coordinate origin is located at the middle of the left bottom edge. In order to clarify the position of each point, Figure 2a shows each point of the patch from xoy plane. The enlargement of the left side of the structure is shown in Figure 2b. The right side of the structure can be obtained through the symmetry principle.
The entire filter structure is composed of substrate and patch. The dielectric substrate with the electrical parameter of ε r = 9.4 holds the dimensions of 19.4 × 4 × 0.38 mm. The gold patches whose parameters can be founded in [43] are located on the surface of the devices. The entire model can be regarded as the center–symmetry model. The excitation port 1 and port 2 with the width and length of 0.34 mm and 0.385 mm are located at the left and bottom of the devices, respectively. The height of substrate and patch along the z-direction are 0.38 mm and 5 mm, respectively. It can be observed from Figure 2c that the entire structure holds the fine details along z-direction. By applying uniform mesh sizes directly to the entire structure, large computational domains will occur resulting in an unacceptable simulation.
The entire computational domain is shown in Figure 3. As can be observed, the entire computational domain holds the parameters of 19.8 × 4 × 3.385 mm. At the boundaries of x-directions and the top boundary of z-direction, the 10-cell-PML is employed for termination to absorb outgoing waves. The other boundaries are terminated by the perfectly electronic conductor (PEC). The source with the center frequency and bandwidth of 20 GHz and 5 GHz propagates at port 1, along the positive side along x-direction which corresponds to the percent bandwidth of % B = 25 . The observation point is located at the corner of port 2 for observing the wave reflection. In order to obtain the best absorption, both in the time domain and frequency domain, the parameters inside the PML regions are chosen to achieve such conditions. The parameters of the higher order PML algorithms are κ η 1 = 12 , α η 1 = 1.1 , m η 1 = 3 , σ η 1 max = 2.6 σ η 1 o p t , κ η 2 = 11 , α η 2 = 0.001 , m η 2 = 4 and σ η 2 max = 3.0 σ η 2 o p t , where
σ η n o p t = ( m η n + 1 ) / ( 150 π Δ )
For comparison, parameters of CFS-PML regions are κ η = 6 , α η = 1.7 , m η = 2 and σ η max = 3.2 σ η o p t . In order to maintain the accuracy of calculation, non-uniform mesh sizes are employed in the HIE algorithm with the parameters of Δ x = Δ y = 0.2 mm and Δ z = 0.0025 mm. The uniform mesh sizes of the conventional FDTD algorithm are chosen as Δ x = Δ y = Δ z = Δ = 0.0025 mm. Thus, time step in the explicit and implicit algorithms are chosen as 0.48 fs and 4.32 fs, respectively. It can be concluded that the CE method can significantly enlarge the time step in the bandpass simulation. The accuracy of the algorithm in the time domain can be evaluated by the waveform obtained by a different PML algorithm in the time domain at the observation point, shown as Figure 4.
Waveforms obtained by explicit algorithms and implicit algorithms with CFLN = 1 are almost overlapped, as shown in Figure 4a. Due to the employment of bandpass-sampling-theorem, the complex envelope method generates the envelope of waveform. Meanwhile, waveforms obtained by the CE method also overlapped with the boundaries of the waveform obtained by the explicit algorithms. Such conditions indicate that these algorithms almost hold the same accuracy. The waveforms obtained by different PML algorithms with large CFLNs are shown in Figure 4b. Compared with the waveform with CFLN = 1, waveforms show significant shifting which indicates the degeneration of calculation accuracy. Compared with the lowpass-sampling-theorem based algorithms, bandpass-sampling-theorem based algorithms show less curves shifting which indicates a better computational accuracy with the circumstance of larger CFLNs. The curve obtained by the proposed scheme is overlapped with that obtained by the CE-LHIE-CPML scheme. Such conditions indicates that these algorithms hold the same accuracy. Meanwhile, it shows the least shifting compared to the other implicit algorithms with larger CFLNs. Such conditions indicates that the proposed algorithm has the best accuracy among these implicit algorithms.
The wave reflection at the boundaries of the domain also affects the accuracy of the calculation. The wave reflection can be reflected by the relative reflection error in the time domain which can be defined as:
R d B ( t ) = 20 log 10 [ | E t ( t ) E r ( t ) | / | max { E r ( t ) } | ]
where E t ( t ) is the test solution which can be directly obtained at the observation point and E r is the reference solution. The reference solution can be calculated by the 20 times enlarged computational domain with 128-cell-PML regions. By employing such a circumstance, reflection wave at the boundaries of the domain can be ignored during the evaluation. Figure 5 shows the relative reflection error in the time domain obtained by different PML algorithms.
The relative reflection error obtained by different PML algorithms with CFLN = 1 is shown in Figure 5a. It can be concluded that the absorption decreases slightly compared with the conventional explicit algorithms due to the introduction of the matrix calculation. Thus, LHIE-PML and LHIE-CPML schemes almost hold the same performance which is inferior compared to the FDTD-PML. The performance can be significantly improved by employing the higher order formulation from both the aspect of late-time reflections and maximum relative reflection error (MRRE). In addition, by employing the CE method, absorption can be significantly improved in the bandpass circumstance. Thus, bandpass-sampling-theorem based algorithms can obtain better performance compared with the FDTD-HPML. It can be concluded that the CE-LHIE-HPML and CE-LHIE-CPML algorithms can obtain better performance compared to the CE-HIE-PML. The reason is that the LHIE algorithm can obtain better numerical dispersion compared to the original HIE algorithm, resulting in such a condition. Figure 5b shows the relative reflection error in the time domain obtained by different PML algorithms with larger CFLNs. The absorption decreases with the enlargement of CFLNs due to the increment of numerical dispersion with the enlargement of time steps. It can also be concluded that LHIE algorithm can obtain better performance compared with the HIE algorithm, which can be reflected by the LHIE-CPML, CE-LHIE-HPML and CE-LHIE-CPML algorithms compared with the HIE-PML and CE-HIE-HPML, respectively. The CE-LHIE-HPML can obtain almost the same, compared with the proposed CE-LHIE-CPML. Most importantly, performance shows its significant advantages with larger CFLNs, compared to the FDTD-PML in bandpass circumstance.
It can be concluded from Table 3 that the employment of higher order formulation and CE method results in the increment of simulation duration and computational resources due to the increment of many more coefficients, auxiliary variables, field components and matrices. Thus, compared with the PML scheme and explicit schemes, time and memory occupied by the HPML and HIE procedure increases significantly. Although the CE method slightly increases the memory consumption during the simulation, iteration steps decrease at the same time due to the enlargement of the time step. It should be noticed that CE method alleviates the time increment problem of the implicit algorithm with lower CFLNs. The proposed scheme can obtain a 13.7% time improvement compared to the FDTD-PML. With the improvement of CFLNs, simulation duration can be decreased by enlarging the time step resulting in an improvement of efficiency. With CFLN = 8, the CE method based implicit algorithm shows more than an 80% time improvement compared with the FDTD-PML. In summary, through the time domain simulation, the proposed scheme can obtain the same accuracy and absorption compared with the previous work which is based on the bandpass-sampling-theorem. Meanwhile, the memory occupation and simulation duration can be decreased significantly.
The above-mentioned discussion verifies the effectiveness of the algorithm in the time domain. In order to further testify the algorithm in the frequency domain, a filter for the sensors system is simulated and manufactured. The scattering parameter is one of the most important parameters in the evaluation. Here, return loss (S11) and transmission loss (S21) are considered in the frequency domain. The photograph of the filter and its sketch size is shown in Figure 6 which is compared with the coin with radius of 25 mm.
In order to demonstrate the performance of the filter in the frequency domain, results are measured through the following system, shown as Figure 7. Due to the fact the excitation ports hold extremely fine details, the excitation source with a probe is employed to excite the filter structure. The probe at each port is also employed to collect the reflection wave and transmission wave. At the other end of the cable, a spectrum analyzer is employed for the evaluation of performance. The return loss and transmission loss can be measured from such a system.
Figure 8 shows S11 and S21 obtained by different PML algorithms and its experiments, respectively. Through Figure 8a,c, it can be observed that these curves obtained by different PML algorithms are almost overlapped. Such conditions indicate that these algorithms almost hold the same computational accuracy during the calculation in the entire frequency band. Through Figure 8b,d, it can be observed that the curves show shifting compared to the lower CFLNs. Such conditions indicate that the accuracy of algorithms decrease with the enlargement of CFLNs. The reason is that the numerical dispersion increases with the enlargement of the time step which results in the degeneration of accuracy. The curves obtained by the proposed scheme and CE-LHIE-HPML hold the least shifting, indicating the best accuracy. It can also be concluded that compared with the lowpass-sampling-theorem based algorithms, algorithms that are based on the CE method hold better accuracy in the entire frequency band. In summary, simulation and measurement shows the same conclusion with that obtained in the time domain simulation. In addition, all of these algorithms show considerable agreement with the experiment. Such a condition indicate that the proposed algorithm can be extensively employed in practical engineering.
However, there is still error between simulation and experiment. The reason can be described from the following aspects: (1) The manufacturing process tolerance results in disagreement compared to the size in the simulation. The generation of shifting is caused by such a condition. Thus, the manufacturing process tolerance overcoming the technique becomes increasingly important; (2) The measurement of the filter system is a precision-type instrument, especially for the probe. The pollution of the devices also leads to shifting in the experiment; (3) Due to the existence of the complex and oblique patch, the conventional conformal method in the Yee’s grid becomes inaccurate. Conformal error causes occur at the oblique and complex part. Such conditions also result in the curves shifting. Thus, improvement on the conformal method for the HIE algorithm is also worth investigating.

3.2. Remote Sensing Problem with Theory Resolution

The performance of the proposed algorithm has been demonstrated by employing the above-mentioned numerical example and its corresponding experiment. Through results, its behaviors show considerable performance during the whole simulation. In order to further demonstrate the efficiency in much more complex problems, a metal sphere model is considered. The normal vector of the sphere model ranges from −180 degree to 180 degree which holds the most complexity among all of the structures. Because the metal sphere model in the sensing problem holds theory resolution, calculation accuracy can be demonstrated through the comparison. Figure 9 demonstrates the sketch picture of the remote sensing problem with the metal sphere.
The entire computational domain holds dimensions of 250 × 250 × 250 in each direction. The metal sphere model with the radius of 50 mm can be regarded as the PEC, located at the center of the domain. The modulated Gaussian pulse with the center frequency and maximum frequency of 1.5 GHz and 2.0 GHz propagates along the negative side along x-direction, whose percent bandwidth can be obtained as % B = 50 . At the boundaries of the domain, all of the boundaries are terminated by 10-cell-PML for wave reflection. The parameters of higher order PML algorithms are κ η 1 = 100 , α η 1 = 2.5 , m η 1 = 2 , σ η 1 max = 0.5 σ η 1 o p t , κ η 2 = 3 , α η 2 = 1.3 , m η 2 = 1 and σ η 2 max = 0.01 σ η 2 o p t . For comparison, parameters of CFS-PML regions are κ η = 170 , α η = 2.1 , m η = 3 and σ η max = 0.8 σ η o p t . In order to make a comparison of computational accuracy between different algorithms, uniform mesh size is employed as Δ x = Δ y = Δ z = Δ = 2.5 mm. Thus, the entire computational domain can be discretized as 100 Δ x × 100 Δ y × 100 Δ z . The time step of the lowpass-sampling-theorem based and bandpass-sampling-theorem based algorithms are 14.4 fs and 57.6 fs, respectively.
The accuracy of the algorithm can be reflected by the radar cross section (RCS) parameter in the frequency domain. Figure 10 shows the RCS obtained by different PML algorithms and CFLNs in the frequency domain and its theory resolution. As can be obtained from Figure 10a, all of the curves are overlapped. Such conditions indicate that all of the algorithms hold the same accuracy with CFLN = 1. However, the curves show shifting compared to the theory resolution. The reason is that numerical dispersion and reflection waves increase at the same time with the enlargement of CFLNs. It can be concluded that the proposed algorithm shows the least shifting. Such conditions indicate that it holds the best accuracy in the frequency domain with larger CFLNs. The effectiveness of the algorithm can also be reflected by the computational duration, consumption memory, iteration step, memory increment and time reduction obtained by different PML algorithms, as shown in Table 4.
It can be concluded from Table 4 that the bandpass-sampling-theorem can significantly improve the calculation efficiency due to the enlargement of time steps and decrement of iteration steps. Compared with the original HIE procedure, the one-step leapfrog scheme can update the components at the half-integer and integer steps, respectively. Such conditions also improve the performance of the entire algorithm. In summary, it can obtain the same conclusion as the previous example, that the proposed scheme holds the best accuracy, efficiency and resources compared to the existed algorithms.
Due to the existence of a massive absorbing boundary condition, the choice of absorbing boundary condition according to different simulation problems has raised much attention. There are several widely-spread absorbing boundary conditions including perfectly matched layer and transparent boundary conditions (TBC) [44]. It can be concluded from the simulation that higher order PML formulation can significantly enhance the absorption at low-frequency band due to the absorption of the low-frequency propagation waves. The employment of higher order PML formulation is available for problems with large amounts of low-frequency waves. For the simulation without much low-frequency waves, the unsplit CFS-PML formulation shows its advantages in late-time reflections and low-frequency evanescent waves. For the higher frequency range applications, due to the non-absorption of plasmonic waves in optical waveguide devices, TBC formulation shows its significant advantages with great potential.

4. Conclusions

In order to efficiently simulate devices for a sensors system in narrow bandpass circumstances, the one-step leapfrog HIE algorithm with the complex envelope method and higher order CPML formulation is proposed. Through the implementation, dielectric material with fine details along one and two dimensions can be solved. The numerical example and its experiment are carried out for the demonstration of effectiveness. It can be concluded from the results that the proposed algorithm shows its advantages in accuracy, efficiency, absorption and computational resources, both in the time domain and frequency domain. Meanwhile, the CFL condition can be broken in the proposed scheme which also indicates the efficiency improvement in the entire simulation. Compared with the previous work, it can maintain considerable performance with decrement of computational resources and simulation duration.
In future work, the method to reduce memory consumption and simulation duration increment caused by matrix solution in the HIE procedure can be further investigated. Meanwhile, numerical dispersion increases with the enlargement of the time step, which affects the calculation accuracy significantly. The method to decrease numerical dispersion with a larger time step is worth investigating. Finally, conformal errors still exist in the complex and oblique structures. Improvements on the conformal method for the HIE algorithm is also worth investigating.

Author Contributions

Conceptualization, H.J. and P.W.; Data curation, Y.W.; Investigation, Y.W., Y.X., H.J. and P.W.; Methodology, H.J. and P.W.; Project administration, Y.X.; Software, Y.W.; Supervision, Y.X.; Writing—review & editing, Y.X. and H.J. All authors have read and agreed to the published version of the manuscript.

Funding

His work is supported by the National Key Research and Development Program of China (2020YFB1807400), National Natural Science Foundation of China (61571022, 61971022, 61801376), Postdoctoral Fellows of “Zhuoyue” Program and National Key Laboratory Foundation of China (HTKJ2019KL504013, 61424020305).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Savelyev, D.; Kazanskiy, N. Near-Field Vortex Beams Diffraction on Surface Micro-Defects and Diffractive Axicons for Polarization State Recognition. Sensors 2021, 21, 1973. [Google Scholar] [CrossRef] [PubMed]
  2. Akib, T.B.A.; Mou, S.F.; Rahman, M.M.; Rana, M.M.; Islam, M.R.; Mehedi, I.M.; Mahmud, M.A.P.; Kouzani, A.Z. Design and Numerical Analysis of a Graphene-Coated SPR Biosensor for Rapid Detection of the Novel Coronavirus. Sensors 2021, 21, 3491. [Google Scholar] [CrossRef] [PubMed]
  3. Theparod, T.; Harnsoongnoen, S. Narrow-Band Light-Emitting Diodes (LEDs) Effects on Sunflower (Helianthus annuus) Sprouts with Remote Monitoring and Recording by Internet of Things Device. Sensors 2022, 22, 1503. [Google Scholar] [CrossRef] [PubMed]
  4. Li, Y.; Wang, N.; Lei, J.; Wang, F.; Li, C. Modeling GPR Wave Propagation in Complex Underground Structures Using Conformal ADI-FDTD Algorithm. Appl. Sci. 2022, 12, 5219. [Google Scholar] [CrossRef]
  5. Wu, P.; Xie, Y.; Jiang, H.; Natsuki, T. Complex Envelope Approximate CN-PML Algorithm with Improved Absorption. IEEE Antennas Wirel. Propag. Lett. 2020, 19, 1521–1525. [Google Scholar] [CrossRef]
  6. Wu, S.; Dong, Y.; Su, F.; Liu, L.; Chen, X. Complex Envelope Approximate Crank-Nicolson Method and Its Open Boundary Implementation for Bandpass Problem. IEEE Access 2021, 9, 45095–45103. [Google Scholar] [CrossRef]
  7. Wu, P.; Wang, X.; Xie, Y.; Jiang, H.; Natsuki, T. Bandpass Simulation of Anisotropic Magnetized Ferrite Material with Alternating Direction Implicit Scheme in Open Region. IEEE Trans. Magn. 2022, 58, 1–8. [Google Scholar] [CrossRef]
  8. Li, Z.; Tan, J.; Guo, Q.; Su, Z.; Long, Y. A Parallel CE-LOD-FDTD Model for Instrument Landing System Signal Disturbance Analyzing. IEEE Trans. Antennas Propag. 2019, 67, 2503–2512. [Google Scholar] [CrossRef]
  9. Singh, G.; Ravi, K.; Wang, Q.; Ho, S.-T. CE-ADI-FDTD method for active photonic device simulation with semiconductor/solid-state media. In Proceedings of the 2012 IEEE Asia-Pacific Conference on Antennas and Propagation, Singapore, 27–29 August 2012; pp. 11–12. [Google Scholar]
  10. Jung, K.; Teixeira, F.L. CE-ADI-FDTD analysis of photonic crystals with a degenerate band edge (DBE). In Proceedings of the 2007 IEEE Antennas and Propagation Society International Symposium, Honolulu, HI, USA, 9–15 June 2007; pp. 4449–4452. [Google Scholar]
  11. Taflove, A.; Hagness, S.C. Computational Electrodynamics: The Finite-Difference Time Domain Method, 3rd ed.; Artech House: Norwood, MA, USA, 2005. [Google Scholar]
  12. Tan, E.L. From Time-Collocated to Leapfrog Fundamental Schemes for ADI and CDI FDTD Methods. Axioms 2022, 11, 23. [Google Scholar] [CrossRef]
  13. Moradi, M.; Nayyeri, V.; Ramahi, O.M. An Unconditionally Stable Single-Field Finite-Difference Time-Domain Method for the Solution of Maxwell Equations in Three Dimensions. IEEE Trans. Antennas Propag. 2020, 68, 3859–3868. [Google Scholar] [CrossRef]
  14. Zhang, K.; Wu, Y.; Liang, W.; Wang, M.; Fan, C.; Zheng, H.; Li, E. Improved Leapfrog LOD-FDTD Method with Controlling Parameters. IEEE Microw. Wirel. Compon. Lett. 2022, 32, 269–272. [Google Scholar] [CrossRef]
  15. Wang, X.; Wang, L.; Zhuo, J.; Lu, X.; Yuan, M.; Zhou, J.; Liu, Q.H. A Hybrid CN-FDTD-SPICE Solver for Field-Circuit Analyses in Low-Frequency Wideband Problems. IEEE Trans. Compon. Packag. Manufact. Technol. 2020, 10, 1721–1728. [Google Scholar] [CrossRef]
  16. Chen, J.; Wang, J. Three-Dimensional Semi-Implicit FDTD Scheme for Calculation of shielding Effectiveness of Enclosure with Thin Slots. IEEE Trans. Electromagn. Compat. 2007, 49, 354–360. [Google Scholar] [CrossRef]
  17. Chen, J.; Wang, J. Numerical Simulation Using HIE-FDTD Method to Estimate Various Antennas with Fine Scale Structures. IEEE Trans. Antennas Propag. 2007, 55, 3603–3612. [Google Scholar] [CrossRef]
  18. Chen, J.; Wang, J. Three-Dimensional Dispersive Hybrid Implicit-Explicit Finite-Difference Time-Domain Method for Simulations of Graphene. Comput. Phys. Commun. 2016, 207, 211–216. [Google Scholar] [CrossRef]
  19. Xiao, F.; Tang, X.; Wang, L. Stability and Numerical Dispersion Analysis of a 3D Hybrid Implicit-Explicit FDTD Method. IEEE Trans. Antennas Propag. 2008, 56, 3346–3350. [Google Scholar] [CrossRef]
  20. Wang, J.; Zhou, B.; Shi, L.; Gao, C.; Chen, B. A Novel 3-D HIE-FDTD Method with One-Step Leapfrog Scheme. IEEE Trans. Micro. Theory Techn. 2014, 62, 1275–1283. [Google Scholar] [CrossRef]
  21. Yao, H.M.; Jiang, L. Machine-Learning-Based PML for the FDTD Method. IEEE Antennas Wire. Propag. Lett. 2019, 18, 192–196. [Google Scholar] [CrossRef]
  22. Berenger, J.P. Perfectly Matched Layer (PML) for Computational Electromagnetics; Morgan & Claypool: San Rafael, CA, USA, 2007. [Google Scholar]
  23. Li, J.; Wu, P. Unconditionally Stable Higher Order CNAD-PML for Left-Handed Materials. IEEE Trans. Antennas Propag. 2019, 67, 7156–7161. [Google Scholar] [CrossRef]
  24. Wu, P.Y.; Jiang, H.L.; Xie, Y.J.; Niu, L.Q. Three-Dimensional Higher Order PML Based on Alternating Direction Implicit Algorithm. IEEE Antennas Wire. Propag. Lett. 2019, 18, 2592–2596. [Google Scholar] [CrossRef]
  25. Fang, Y.; Liu, J.; Jiao, Z.; Bai, G.; Xi, X. A 3-D Stochastic FDTD Algorithm for Wave Propagation in Isotropic Cold Plasma Medium Based on Bilinear Transform. IEEE Trans. Plasma Sci. 2019, 47, 173–178. [Google Scholar] [CrossRef]
  26. Jiang, H.L.; Cui, T.J. Efficient PML Implementation for Approximate CN-FDTD Method. IEEE Antennas Wire. Propag. Lett. 2019, 18, 698–701. [Google Scholar] [CrossRef]
  27. Giannopoulos, A. Higher-Order Convolution PML (CPML) for FDTD Electromagnetic Modeling. IEEE Trans. Antennas Propag. 2020, 68, 6226–6231. [Google Scholar] [CrossRef]
  28. Chen, L.; Ozakin, M.B.; Bagci, H. A Low-Storage PML Implementation within a High-order Discontinuous Galerkin Time-Domain Method. In Proceedings of the 2020 IEEE International Symposium on Antennas and Propagation and North American Radio Science Meeting, Montreal, QC, Canada, 5–10 July 2020; pp. 1065–1066. [Google Scholar]
  29. Liu, J.; He, S.; Cheng, W. Pinpoint and Efficient DZT-based FDTD Implementations Using Optimal 2nd-order PML Truncation. In Proceedings of the 2019 Photonics & Electromagnetics Research Symposium-Fall (PIERS-Fall), Xiamen, China, 17–20 December 2019; pp. 897–902. [Google Scholar]
  30. Haolin, J.; Yongjun, X.; Peiyu, W.; Jianfeng, Z.; Liqiang, N. Unsplit-field higher-order nearly PML for arbitrary media in EM simulation. J. Syst. Eng. Electron. 2021, 32, 1–6. [Google Scholar] [CrossRef]
  31. Zhang, K.; Wang, L.; Wang, M.; Zheng, H.; Li, E. Efficient Implementation of the CPML in 3-D Hybrid Implicit-Explicit FDTD Method. In Proceedings of the 2019 IEEE International Conference on Computational Electromagnetics (ICCEM), Shanghai, China, 20–22 March 2019; pp. 1–3. [Google Scholar]
  32. Liu, Y.; Shi, L.; Chen, H.; Wang, J.; Lei, Q.; Duan, Y.; Zhang, Q.; Fu, S. One-Step Leapfrog LOD-FDTD Method with CPML for Periodic Structures. IEEE Antennas Wire. Propagat. Lett. 2019, 18, 9–13. [Google Scholar] [CrossRef]
  33. Liu, Z.; Chen, Y.; Sun, X.; Liu, Y. Implementation of CFS-PML for HIE-FDTD Method. IEEE Antennas Wirel. Propag. Lett. 2012, 11, 381–384. [Google Scholar]
  34. Van Londersele, A.; De Zutter, D.; Ginste, D.V. A Collocated 3-D HIE-FDTD Scheme with PML. IEEE Microw. Wirel. Compon. Lett. 2017, 27, 609–611. [Google Scholar] [CrossRef]
  35. Zhang, C.; Kong, Y. Nearly PML for a Novel Four-Steps Weakly Conditionally Stable HIE-FDTD Method. In Proceedings of the 2019 IEEE International Conference on Computational Electromagnetics (ICCEM), Shanghai, China, 20–22 March 2019; pp. 1–3. [Google Scholar]
  36. Wang, J.; Wang, J.; Zhou, B.; Wu, Q.; Gao, C.; Guo, F. HIE-FDTD Method with PML for 2-D Periodic Structures at Oblique Incidence. IEEE Antennas Wirel. Propag. Lett. 2016, 15, 984–987. [Google Scholar] [CrossRef]
  37. Lei, Q.; Wang, J.; Shi, L. Modified HIE-FDTD Algorithm for 2-D Periodic Structures with Fine Mesh in Either Direction. IEEE Antennas Wirel. Propag. Lett. 2017, 16, 1455–1459. [Google Scholar] [CrossRef]
  38. Dong, M.; Zhang, A.; Chen, J. Perfectly matched layer for hybrid implicit and explicit-FDTD method. In Proceedings of the 2015 IEEE MTT-S International Microwave Workshop Series on Advanced Materials and Processes for RF and THz Applications (IMWS-AMP), Suzhou, China, 1–3 July 2015; pp. 1–3. [Google Scholar]
  39. Wu, P.; Yu, H.; Xie, Y.; Jiang, H.; Natsuki, T. Efficient Enhanced Hybrid Implicit-Explicit Procedure to Gyrotropic Plasma in Open Regions with Fine Geometry Details Along Single Direction. IEEE Access 2021, 9, 77079–77089. [Google Scholar] [CrossRef]
  40. Li, J.; Yang, Q.; Niu, P.; Feng, N. Novel implementation of the CFS-PML and its validation for dispersive media. In Proceedings of the 2011 IEEE International Conference on Microwave Technology & Computational Electromagnetics, Beijing, China, 22–25 May 2011; pp. 422–425. [Google Scholar]
  41. Wu, P.; Xie, Y.; Jiang, H.; Natsuki, T. Modeling of Bandpass GPR Problem by HIE Procedure With Enhanced Absorption. IEEE Geosci. Remote Sens. Lett. 2022, 19, 1–5. [Google Scholar] [CrossRef]
  42. Liu, Y.; Yao, Y.; Li, Y.; Yu, G.; Ni, J. Parallel implementation and application of the higher-order FDTD with convolution PML. In Proceedings of the 2016 IEEE MTT-S International Conference on Numerical Electromagnetic and Multiphysics Modeling and Optimization (NEMO), Beijing, China, 27–29 July 2016; pp. 1–2. [Google Scholar]
  43. Mohsin, A.S.M.; Salim, M.B. Probing the Plasmon Coupling, Quantum Yield, and Effects of Tip Geometry of Gold Nanoparticle Using Analytical Models and FDTD Simulation. IEEE Photonics J. 2018, 10, 1–10. [Google Scholar] [CrossRef]
  44. Akbar, S.; Agarwal, N.; Gupta, S. Elliptical hybrid cladding borosilicate photonic crystal fiber design for minimum chromatic dispersion. In Proceedings of the 2014 Eleventh International Conference on Wireless and Optical Communications Networks (WOCN), Vijayawada, India, 11–13 September 2014; pp. 1–4. [Google Scholar]
Figure 1. The block diagram of the entire update procedure.
Figure 1. The block diagram of the entire update procedure.
Sensors 22 04445 g001
Figure 2. The sketch picture of the filter for sensors system: (a) top view; (b) top view enlargement; (c) front view.
Figure 2. The sketch picture of the filter for sensors system: (a) top view; (b) top view enlargement; (c) front view.
Sensors 22 04445 g002
Figure 3. The sketch picture of the entire filter computational domain.
Figure 3. The sketch picture of the entire filter computational domain.
Sensors 22 04445 g003
Figure 4. Waveform obtained by different PML algorithms at the observation point: (a) CFLN = 1; (b) CFLN = 8.
Figure 4. Waveform obtained by different PML algorithms at the observation point: (a) CFLN = 1; (b) CFLN = 8.
Sensors 22 04445 g004
Figure 5. The relative reflection error obtained by different PML algorithms in the time domain: (a) CFLN = 1; (b) CFLN = 8.
Figure 5. The relative reflection error obtained by different PML algorithms in the time domain: (a) CFLN = 1; (b) CFLN = 8.
Sensors 22 04445 g005
Figure 6. The picture of: (a) entire filter model; (b) detail size of the filter with ruler; (c) detail size with coin with the radius of 25 mm.
Figure 6. The picture of: (a) entire filter model; (b) detail size of the filter with ruler; (c) detail size with coin with the radius of 25 mm.
Sensors 22 04445 g006
Figure 7. The picture of: (a) probe operation table; (b) measured system.
Figure 7. The picture of: (a) probe operation table; (b) measured system.
Sensors 22 04445 g007
Figure 8. The scattering parameter obtained by different PML algorithms: (a) SCFLN = 1; 1with CFLN = 1; (b) S11 with CFLN = 8; (c) S21 with CFLN = 1; (d) S21 with CFLN = 8.
Figure 8. The scattering parameter obtained by different PML algorithms: (a) SCFLN = 1; 1with CFLN = 1; (b) S11 with CFLN = 8; (c) S21 with CFLN = 1; (d) S21 with CFLN = 8.
Sensors 22 04445 g008
Figure 9. The sketch picture of the metal sphere model inside the entire computational domain.
Figure 9. The sketch picture of the metal sphere model inside the entire computational domain.
Sensors 22 04445 g009
Figure 10. Radar cross section obtained by different PML algorithms: (a) CFLN = 1; (b) CFLN = 8.
Figure 10. Radar cross section obtained by different PML algorithms: (a) CFLN = 1; (b) CFLN = 8.
Sensors 22 04445 g010
Table 1. The multiplication/division and addition/subtraction operators with explicit and implicit equations in different algorithms.
Table 1. The multiplication/division and addition/subtraction operators with explicit and implicit equations in different algorithms.
PML AlgorithmAddition/SubtractionMultiplication/DivisionTotal Operators
ImplicitExplicitImplicitExplicit
FDTD-PML060042102
FDTD-HPML090078168
HIE-PML40543628158
LHIE-CPML36543228150
CE-HIE-HPML66725436228
CE-LHIE-HPML54723636198
Proposed48664036190
Table 2. The specific coordinate location of each point in Figure 2a,b on the surface of the model in (x, y, z) form (unit: mm).
Table 2. The specific coordinate location of each point in Figure 2a,b on the surface of the model in (x, y, z) form (unit: mm).
Point in
Figure 2a,b
Specific Coordinate Location
(x, y, z)
Point in
Figure 2a,b
Specific Coordinate Location
(x, y, z)
1(−9.7, −0.175, 0.005)2(−5.3, −0.175, 0.005)
3(−3.4, 0.48, 0.005)4(−3.46, 0.59, 0.005)
5(−4.4, 0.175, 0.005)6(−9.7, 0.175, 0.005)
7(−4.7, −0.25, 0.005)8(−4.64, −0.25, 0.005)
9(−3.3, 0.22, 0.005)10(−3.27, 0.18, 0.005)
11(−1.9, 0.77, 0.005)12(−2.0, 0.93, 0.005)
13(−3.1, −0.3, 0.005)14(−3, −0.45, 0.005)
15(−1.64, 0.11, 0.005)16(−1.64, 0.11, 0.005)
17(−3.1, 0.69, 0.005)18(−0.4, 0.88, 0.005)
19(−3.07, −0.29, 0.005)20(−3, −0.45, 0.005)
Table 3. The computational duration, consumption memory, iteration step, memory increment and time reduction obtained by different PML algorithms.
Table 3. The computational duration, consumption memory, iteration step, memory increment and time reduction obtained by different PML algorithms.
PML AlgorithmCFLNStepsMemory (GB)Memory
Increment (%)
Time (min)Time
Reduction (%)
FDTD-PML165,5360.5-21.9-
FDTD-HPML165,5360.98038.4−75.3
HIE-PML165,5360.98041.7−90.4
LHIE-CPML165,5360.86038.3−42.8
CE-HIE-HPML172821.418024.0−9.5
CE-LHIE-HPML172821.316022.1−9.0
CE-LHIE-CPML172821.112018.913.7
HIE-PML881920.98012.642.5
LHIE-CPML881920.86010.353.0
CE-HIE-HPML89111.41804.679.0
CE-LHIE-HPML89111.31603.584.0
CE-LHIE-CPML89111.11202.986.8
Table 4. The computational duration, consumption memory, iteration step, memory increment and time reduction obtained by different PML algorithms in the sphere model.
Table 4. The computational duration, consumption memory, iteration step, memory increment and time reduction obtained by different PML algorithms in the sphere model.
PML AlgorithmCFLNStepsMemory (GB)Memory
Increment (%)
Time (min)Time
Reduction (%)
FDTD-PML165,5360.4-4.6-
FDTD-HPML165,5360.8−10010.2−152.4
HIE-PML165,5360.8−10013.7−197.8
LHIE-CPML165,5360.7−7510.9−137.0
CE-HIE-HPML172821.1−1754.9−6.5
CE-LHIE-HPML172821.0−1504.7−2.2
CE-LHIE-CPML172820.8−1004.44.3
HIE-PML881920.8−1002.745.7
LHIE-CPML881920.7−752.350.0
CE-HIE-HPML89111.1−1750.784.8
CE-LHIE-HPML89111.0−1500.687.0
CE-LHIE-CPML89110.8−1000.491.3
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, Y.; Xie, Y.; Jiang, H.; Wu, P. Narrow-Bandpass One-Step Leapfrog Hybrid Implicit-Explicit Algorithm with Convolutional Boundary Condition for Its Applications in Sensors. Sensors 2022, 22, 4445. https://doi.org/10.3390/s22124445

AMA Style

Wang Y, Xie Y, Jiang H, Wu P. Narrow-Bandpass One-Step Leapfrog Hybrid Implicit-Explicit Algorithm with Convolutional Boundary Condition for Its Applications in Sensors. Sensors. 2022; 22(12):4445. https://doi.org/10.3390/s22124445

Chicago/Turabian Style

Wang, Yangjing, Yongjun Xie, Haolin Jiang, and Peiyu Wu. 2022. "Narrow-Bandpass One-Step Leapfrog Hybrid Implicit-Explicit Algorithm with Convolutional Boundary Condition for Its Applications in Sensors" Sensors 22, no. 12: 4445. https://doi.org/10.3390/s22124445

APA Style

Wang, Y., Xie, Y., Jiang, H., & Wu, P. (2022). Narrow-Bandpass One-Step Leapfrog Hybrid Implicit-Explicit Algorithm with Convolutional Boundary Condition for Its Applications in Sensors. Sensors, 22(12), 4445. https://doi.org/10.3390/s22124445

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop