A RT-FDTD method of analyzing wireless propagation characteristics in underground mine

Efficient communication is crucial in reducing injuries and fatalities in coal mine accidents, necessitating the study of simulation methods for mine communication. When transceiver antennas are positioned close to the same side of the tunnel, the simulation results from the Ray Tracing (RT) method exhibit significant errors. Additionally, the Finite-Difference Time-Domain (FDTD) method demands substantial computational resources. In response to these challenges, we propose a RT-FDTD method, guided by the law of conservation of energy. This approach involves dividing the mine tunnel into a cuboidal region, using the RT method to calculate the electric field strength on the cuboid’s surface, and then employing this as the excitation source for the FDTD method. Subsequently, the FDTD method is used to calculate the electric field strength within the cuboid. Experimental results demonstrate that the RT-FDTD method effectively mitigates the limitations of the RT and FDTD methods, enhancing both the efficiency and accuracy of simulations in underground mine.

• We propose a RT-FDTD method to analyze the propagation characteristic analysis in underground mine, which combines RT method's high efficiency and FDTD method's precision using law of conservation of energy; and it expands application range for RT method and FDTD method in wireless propagation characteristics analysis for a mine.• The RT-FDTD method will calculate signal strength for different frequency electromagnetic waves consist- ently with results calculated by FDTD method and unit window recycling FDTD (UWR-FDTD) method as well as measured values in this paper, no matter where the transmitting and receiving antennas are located in the tunnel section.• The RT-FDTD method greatly improves the efficiency of electromagnetic wave propagation characteristics in technical mine and has high computational accuracy.
The rest of this paper is organized as follows.In "RT-FDTD model" section, we introduce the RT-FDTD model.
In "Numerical calculation and analysis" section, we conduct simulation experiments and underground measurement experiments using this model.Following this, we analyze the accuracy and efficiency of the simulation data.The paper concludes in "Conclusions" section.

RT-FDTD model RT-FDTD principle
Combining RT and FDTD method has been used to improve LED correlated color temperature uniformity 21 and set multiscale optical simulation 22 ,and the needed computational space is relatively small in these application scenarios.The RT method is one of the main analytical methods for simulating the characteristics of wireless propagation in long and straight tunnels 23 .It treats the high-frequency electromagnetic waves radiated from the transmitter as multiple rays, and expresses the electric field strength at a point in the tunnel by the sum of the direct and reflected rays arriving at that point.The RT method can accurately calculate the electromagnetic wave properties in the tunnel as long as the transmitter and receiver are not close to the same side of the tunnel wall at the same time.FDTD method calculation of wireless propagation characteristics in mine requires discretizing the whole tunnel space into multiple Yee grids.Based on the spatial and temporal relationships between the electric and magnetic fields, the electric and magnetic field strengths can be calculated for each Yee grid.The FDTD method has a high computational accuracy throughout the tunnel.In order to prevent numerical dispersion of electromagnetic waves during propagation after spatial discretization 24 , the spatial step of the Yee grid generally does not exceed 1/10 of the electromagnetic wave wavelength 25 , while the time step and spatial step also need to satisfy CFL condition 26,27 .Restricted by these two conditions, the higher the electromagnetic wave frequency, the shorter the Yee grid spatial step, and the more Yee grids need to be divided within the same tunnel.Indeed, while some methods [17][18][19] have made strides in addressing this limitation, it's important to note that there is a close relationship between computational error and spatial step size.This relationship can significantly impact the accuracy of simulations, making it a crucial factor to consider in the development and application of these methods.The direct use of the FDTD method in the tunnel to calculate the wireless propagation characteristics requires huge memory space and long simulation time.
In the RT-FDTD model, a vertical surface is set along the length of the tunnel from an area near the center of the tunnel cross-section, and this vertical surface forms a cuboid region with the tunnel roof, one side of the tunnel wall, and the tunnel floor.The vertical surface can also form a cuboid area with a plane parallel to the tunnel roof, a side of the tunnel, and a plane parallel to the tunnel floor.The cuboid region is divided into two parts: the surface of the cuboid body and the internal space of the cuboid body.The size of the grid in the internal space of the cuboid and the size of the discrete grid on the surface of the cuboid are determined by the FDTD method.The electric field strength of each discrete grid on the surface of the cuboid is calculated by RT method.The inside region of the cuboid is discretized with Yee grids, and the FDTD iterative formula is executed to calculate the electromagnetic field strength in the cuboid.According to the law of conservation of energy, the electric field strength of the entire surface of the cuboid obtained by the RT method can be used as the excitation source of the FDTD method.The RT-FDTD method can greatly reduce the memory space and simulation time required by FDTD method in the same simulation space, and ensure the calculation accuracy and calculation efficiency when the transceiver point is located near the tunnel wall at the same time.The schematic for the RT-FDTD model is shown in Fig. 1.The x direction is the tunnel width direction, the y is the tunnel height direction, and the z is the tunnel length direction.
In Fig. 1, the gray areas of the tunnel roof, floor and both sides of the tunnel are the tunnel walls.These four surfaces form a complete simulation of the tunnel using the same electromagnetic parameters.To give a clearer picture of Fig. 1, the tunnel roof is not represented in gray, only a discrete cuboid region named A FDTD is deline- ated for the area near wall L, which is discretized with a Yee grid.The calculation area for the RT method is the entire tunnel space, and the emission point T can be located anywhere within the tunnel section.The surfaces labeled S RF1 , S RF2 , S RF3 and S RF4 are the interfaces of the calculation areas of RT and FDTD.S RF1 is close to and parallel to the tunnel wall L, S RF2 is close to the emission point and parallel to the coordinate plane xoy, S RF3 is close to and parallel to the roof, and S RF4 is close to and parallel to the floor.Assuming that the receiving point R is on the face S RF1 , we calculate the electric field strength at the receiving point R using the RT method, with the ray labelled 1 indicating the direct path.The ray labelled 2 represents the first reflection path after being reflected by the floor, and the first reflection point is P 11 .The ray labeled 3 represents the second reflection path after being reflected by the wall L and R, and the second reflection points are P 21 and P 22 , respectively.The ray 4 represents the 3rd reflection path after being reflected by the wall L, R and L successively, and with the 3rd reflection points are P 31 , P 32 and P 33 respectively, where the reflection points P 31 and P 33 are located on the wall L and the reflection point P 32 is located on the wall R. The yellow grid region of the cuboid A FDTD is the FDTD method computing space, and the size of each Yee grid in this region determine the spatial resolution of the RT method to calculate the electric field strengths on the surfaces named S RF1 , S RF2 , S RF3 and S RF4 .Based on the law of conservation of energy, the electric field strengths on surfaces named S RF1 , S RF2 , S RF3 and S RF4 can be used as the excitation source for calculating the interior region of the A FDTD using the FDTD method and calculate the electric field strengths for each Yee grid by the standard FDTD iterative formulas.Let the coordinate origin be at the center of the tunnel section, the tunnel width is W, the tunnel height is H and the tunnel length is L also the coordinate range of surface S RF1 along the width (x direction), height (y direction) and length (z direction) of tunnel be satisfied where x RF1 , y RF1 , z RF1 are the coordinate values of each Yee grid in the surface S RF1 .The coordinate range of the surface S RF2 along the width (x-direction), height (y-direction) and length (z-direction) of the tunnel should be satisfied where x RF2 , y RF2 , z RF2 are the coordinate values of each Yee grid in the surface S RF2 .The coordinate range of surface S RF3 along the width (x-direction), height (y-direction) and length (z-direction) of the tunnel be satisfied where x RF3 , y RF3 , z RF3 are the coordinate values of each Yee grid in the surface S RF3 .The coordinate range of surface S RF4 along the width (x-direction), height (y-direction) and length (z-direction) of the tunnel be satisfied where x RF4 , y RF4 , z RF4 are the coordinate values of each Yee grid in the surface S RF4 .

Calculation of the RT region
According to (1), the distance between the surface S RF1 and the tunnel wall along the coordinates in the x direc- tion is at least 1/4 of the tunnel width.In this case, the distance between the transmitting point and the tunnel wall does not affect the accuracy of calculating the electric field strength on the surface S RF1 by the RT method, which avoids the shortcoming of large error when the RT method is used to calculate the transceiver antenna at the same time near the tunnel wall.Set the origin of coordinates at the center of the tunnel, and set the coordinate of the transmitting point T as (x 0 , y 0 , 0) , then the image point coordinates of the transmitting point with respect to the four tunnel walls are satisfied where p and q are integers.A positive p indicates the number of reflections of the transmitting point relative to the tunnel roof and a negative p indicates the number of reflections from the tunnel floor.A positive q indicates the number of reflections from the transmitting point relative to the tunnel wall with positive coordinates and a negative q indicates the number of reflections from the transmitting point relative to the tunnel wall with negative coordinates.Figure 2 shows a schematic diagram of solving for the partial image points of the transmitting point T on different wall surfaces.The width of the simulated tunnel is 2a and the height is 2b.
The distance between the ray radiated by the transmitting point T and the receiving point R after reflection p and q times is equal to the distance between the image point I p,q of the transmitting point and the receiving point R. Assuming that the coordinate of the receiving point R is (x, y, z), then the distance named r pq from R to the image point I p,q is satisfied According to RT theory, the Angle between incident ray and normal vector of tunnel wall, roof or floor is set as θ and θ ⊥ respectively, and the magnitude satisfies can be expressed as The incident wave electric field vector can be decomposed into horizontal and vertical polarized electric field vectors, named E t and E t⊥ (1) θ � = arccos x − x p r p,q θ ⊥ = arccos y − y q r p,q Vol.:(0123456789) where n is the normal vector of the plane where the reflection point is located, r is the incident wave direction vector.
Assuming that the same building materials are used so that the electromagnetic parameters are the same on tunnel walls.The reflection coefficients on both sides of tunnel walls and roof and floor can be obtained by using the angle between incident rays and reflective surfaces and the electromagnetic parameters of the walls, named Ŵ and Ŵ ⊥ where ε r = ε r − jσ/ωε 0 is the complex dielectric constant of tunnel wall, ε r is the relative dielectric constant of tunnel wall, ω is the angular frequency of electromagnetic wave, σ is the conductivity of tunnel wall, ε 0 is the dielectric constant of vacuum.
The RT method calculates an expression for the electric field strength at the receiving point based on the length of the direct path, the length of the reflected path for different numbers of reflections, and through Eq. ( 8).
According to the length of the direct path and the length of the reflected path with different reflection times, the expression of the electric field strength Er at the receiving point can be calculated by Eq. ( 8), that is where j is an imaginary number unit, k is the electromagnetic wave number.As shown in Fig. 1, the electric field strengths on S RF1 , S RF2 , S RF3 and S RF4 are calculated using Eq. ( 10) as the excitation source for the FDTD calcula- tion region inside the cuboid A FDTD , respectively.The spatial sampling intervals for the electric fields on each of the above surfaces obtained by the RT method are the same as the Yee grid sizes of the FDTD calculation region.

RT-FDTD calculation region
Modeling the entire tunnel using FDTD method directly requires huge memory space.Firstly, the RT method is used to calculate the electric field on a surface near the tunnel center, which is used as the excitation source of the FDTD calculation area.Then, the FDTD method is used to model the area near the tunnel wall.The RT-FDTD method avoids the direct division of Yee grid in the entire tunnel space, saves a lot of memory space and simulation time, and improves simulation efficiency.Figure 3 shows the model for calculating the internal spatial electric field strength of the cuboid A FDTD in Fig. 1 using the FDTD method.
The area surrounded by blue dashed lines in Fig. 3 is the FDTD calculation region, which mainly includes three parts, that is, Perfectly Matched Layer (PML), including tunnel space, excitation source surface.The tunnel space is filled with air medium, and the electric fields on S RF1 , S RF2 , S RF3 and S RF4 excitation sources composed ( 8) of multiple receiving points have been calculated by RT method.PML is wrapped in the outer sides of the tunnel wall L, roof and floor to absorb the electromagnetic wave transmitted to the tunnel wall and prevent it from reflecting back to the tunnel space.The excitation source surface is directly used as a boundary of the FDTD calculation region, and its electric field strength is only related to the calculation results of the RT method, so it is not necessary to set the PML in the excitation source surface.
The red lines S RF1 , S RF2 , S RF3 and S RF4 in the main view of FDTD calculation area represent the excitation sources calculated by RT method as FDTD update equation.In the top view, it can be seen that there are two surfaces, and, as FDTD excitation sources, and represents the electric field parallel to the xoy plane as FDTD excitation sources in the xoy direction.Once the x coordinate of surface A is determined, the number of grids that FDTD needs to divide in that direction can be calculated according to the space step size in the x direction of Yee grid.The number of Yee grids required in the y direction is determined by the y coordinates of surface, and the spatial step size of the Yee grids in that direction.The required number of Yee grids in the z direction is determined by the length of the simulated tunnel the difference of the length from the transmitting point to the surface and the spatial step size of the Yee grid in that direction.
The spatial step and temporal step in the FDTD iterative equations need to satisfy the CFL stability condition.
Where t is the discrete time step, x, y, z correspond to the Yee grid spatial step size in the three direc- tions of the coordinate axis respectively; µ is the permeability of electromagnetic wave propagation medium; ε is the dielectric constant of electromagnetic wave propagation medium.
The four surface excitation sources, such asS RF1 , S RF2 , S RF3 and S RF4 are directly coupled to the Yee grid inside the cuboid, and iterated along the x direction from the surface to the wall L in accordance with the temporal step.Finally, the electric field strength distribution in the whole tunnel is obtained.
The FDTD iterative equation for solving the electric field of the Yee grid inside the cuboid are where E x , E y , E z , H x , H y , H z represent the electric field strength and magnetic field strength in 3 directions such as x, y, z, respectively.The superscripts n, n + 1/2, n + 1 for the electric and magnetic field strengths denote the discrete time step of the FDTD method.The characters i

Analysis of calculation accuracy
While there are numerous electromagnetic simulation methods available, most are typically used to analyze small targets, with only a few being applied in underground mine environments.In our study, we have referred to various articles and selected the UWR-FDTD 16 , RT 14 , and RT-FDTD methods.These methods were used to calculate and compare the wireless propagation characteristics of 740 MHz electromagnetic waves.The environment for these simulations was a rectangular tunnel with dimensions of 200 m in length, 4.8 m in width, and 3.4 m in height.This was done to analyze the calculation accuracy of the RT-FDTD method.The construction materials of both sides of the tunnel wall, roof and floor are concrete, so they have the same electromagnetic parameters around the tunnel: a relative permittivity of 8 and a conductivity of 1e−2 S/m 28 .The Yee grid divided by FDTD method in the internal region the cuboid has a uniform spatial grid step in all three directions.When this spatial step is 1/10 of the electromagnetic wave wavelength, it is 0.04 m and the corresponding temporal step is 70 ps to satisfy the CFL condition according to Eq. ( 9).
The origin of coordinates is located in the center of the tunnel section, and we set their coordinate parameters according to the positions of the transmitting point and surface excitation source in Fig. 4.
Most numerical simulations or measurement experiments studying wireless propagation characteristics always set the transmitting and receiving points in the central area of the tunnel section.However, We all know that wireless communication base station antennas are installed on the walls or roof of tunnels.Therefore, when designing verification experiments, it is indeed crucial to not only consider the central area of the tunnel section but also the areas close to the tunnel walls.This approach will provide a more comprehensive understanding of the performance of the RT-FDTD method in real-world scenarios, which has significant practical implications.
Firstly, we calculated and compared the results when we placed the transmitting point in the center of a tunnel section with length, width and height of 200 m, 4.8 m and 3.4 m respectively; and we placed the receiving point close to a tunnel wall L at a distance of either 1/4 or 1/8 of the tunnel width.The coordinate of the transmitting point is (0,0,0).We set the x-coordinate value of a surface excitation source to be 0, which is half of the tunnel width.The distance from this surface excitation source S RF2 to the emission point along the z-axis is 4 m; and the y-coordinates of the surface excitation source of S RF2 and S RF3 are 0.8 and − 0.8 respectively.We used RT method to calculate the electric field strength at 4 surface excitation sources for a 740 MHz electromagnetic wave when we placed its emission point at the center of a tunnel section.Then we used these calculated electric (13) field strengths as FDTD excitation sources.We executed Eqs.(12-17) in a cuboid's internal space according to required time steps along x-axis using FDTD method to solve for electric field strength at each Yee grid in a cuboid composed of surfaces S RF1 , S RF3 , S RF4 and wall L. We used UWR-FDTD method from reference 13 to divide eight virtual windows along z-axis of a tunnel with unit window length, width and height being 25 m, 4.8 m and 3.4m respectively.We placed a transceiver antenna in a tunnel with sectional coordinates (− 1.2,0), which means its distance from wall L is o1/4 of a tunnel's width.In a same tunnel environment as UWR-FDTD's, we used RT method from reference 16 and our RT-FDTD method to calculate electric field strength at each receiving sampling point along z-axis of a tunnel.To compare with measured values, we converted calculated electric field strengths into received signal power expressed in dBm.The results are shown in Fig. 5.  www.nature.com/scientificreports/ Figure 5 shows that when we place a transmitting point at a tunnel section's center and place a receiving point at a distance of 1/4 of a tunnel's width from wall L, UWR-FDTD, RT and our RT-FDTD method calculate similar electromagnetic wave signal strength variation patterns for each longitudinal sampling in a tunnel.Since RT-FDTD method uses some RT method's calculation results, RT method and our RT-FDTD method have very similar calculation results.We keep transmitting point's position and excitation source plane's position unchanged; and place a receiving point on a tunnel section with coordinate (− 1.8,0), which means its distance from wall L is 1/8 of a tunnel's width as shown in Fig. 4. Figure 6 shows electromagnetic wave signal strength at each receiving sampling point calculated along z-axis of a tunnel.
In Fig. 6, when we place a transmitting point at a tunnel section's center and place a receiving point at a distance of 1/8 of a tunnel's width from wall L, all three calculation methods show a consistent attenuation pattern for received electromagnetic wave signal strength.
We keep transmitting point's position unchanged and place a receiving point at a distance of 0.04 m from a tunnel wall.When receiving sampling points are identical, RT method calculates a received signal power that is 12 dB higher than what UWR-FDTD method and our RT-FDTD method calculate; and UWR-FDTD method and our RT-FDTD method calculate identical received signal power variation patterns for each receiving sampling point.Figure 7 shows these results.
UWR-FDTD method's iterative process resembles traditional FDTD method's iterative process.Multipath affects signal strength calculated by this method more when it is near a tunnel wall, and large-scale fading causes more obvious amplitude fluctuations for received signal.Surface source's distance from a tunnel wall is 2.4 m.Since RT-FDTD method's excitation source mainly comes from RT method, RT-FDTD method calculates a relatively gentle amplitude fluctuation for received signal strength.RT method cannot accurately calculate wireless  www.nature.com/scientificreports/propagation characteristics in an underground mine when both transceiver antennas are near a tunnel wall simultaneously.So, we use hybrid RT method's simulation value from reference 14 and corresponding position's measured value to compare with RT-FDTD method's calculation accuracy near a tunnel wall simultaneously.We place a transmitting point near wall L with coordinate (− 2.36,0,0).Surface excitation source 's distance from wall L is 1.2 m; and its distance from transmitting point along z-axis is 4 m.Both surface excitation sources have identical y-coordinates that are 1/4 of a tunnel's height.We place each receiving point at a distance of 0.04 m from a tunnel wall along z-axis to simulate wireless propagation characteristics between two wireless communication base station antennas installed on same side of a tunnel wall.The 700MHz frequency band is a traditional broadcasting and television system frequency band in China.As a golden frequency band for conducting mobile communication services, the 700MHz frequency band has good propagation characteristics and this band was adjusted for 5G communication by the Ministry of Industry and Information Technology of China in 2020.So we conducted many measurement tests underground at Shuangma Coal Mine of Ningxia Coal Industry CO., LTD; and measured wireless propagation characteristics for electromagnetic waves with frequencies of 580 MHz and 740 MHz when we placed transmitting and receiving antennas at different positions on a tunnel section to promote 5G communication system in underground mine.
Testing tunnel's cross-section is horseshoe-shaped;and based on the research results in the reference 29 , we can equate it to a rectangular tunnel with width of 4.8 m and height of 3.4 m using area equivalence principle.Figures 8 and 9 show measurement environment and equivalent tunnel respectively.Measuring device's RF transmitting power is 30 dBm; and receiving device's sensitivity is − 92 dBm at a communication rate of 1 Mbps.www.nature.com/scientificreports/ We used a vertically polarized omnidirectional dipole antenna with characteristic impedance of 50 and gain of 1 dBi.The distance from the transmitting and receiving antennas to the floor is 1.70 m and the max distance from the transmitting antenna to the receiving antenna is 200 m along the length direction of the tunnel.We place transceiver antenna at a distance of 0.04 m from same side's tunnel wall.We compare simulation values calculated by UWR-FDTD method, mixed RT method and our RT-FDTD method with measured values.Figure 10 shows received signal strength for a 580 MHz frequency electromagnetic wave along length direction of a tunnel; and Fig. 11 shows received signal strength for a 740 MHz frequency electromagnetic wave.
Figure 10 shows that when both transceiver antenna and tunnel wall L are at a distance of 0.04 m from each other, three methods-UWR-FDTD, hybrid RT and our RT-FDTD calculate a consistent attenuation trend for 580 MHz electromagnetic wave signal strength along a tunnel's length; and measurement results agree with this trend.Figure 11 shows similar results for 740 MHz electromagnetic wave's received signal strength calculated along a tunnel's length.Analysis results show that we can use our RT-FDTD method to calculate electromagnetic wave propagation characteristics when both transceiver antennas are near a tunnel wall simultaneously; and this solves large error problem in RT method's simulation results for this scenario.

Computational efficiency
We study computational efficiency of RT-FDTD, FDTD, UWR-FDTD methods for the same simulation tunnel through theoretical analysis and simulation experiments.
An electromagnetic wave's wavelength is ; and we discretize each wavelength into N d grids with each grid's length being N d / .For a 3D space with length L, width W and height H; we divide it into N Yee grids.We need to allocate memory space for each electromagnetic field component in a grid.An electromagnetic wave's wavelength is fixed for a same frequency.To prevent numerical dispersion for an electromagnetic wave, we usually set the value of N d to be equal or greater than 10.Once we fix the value of N d , we consider it as a constant.In this case, space size to be solved determine number of Yee grids.We use RT method in our RT-FDTD method to calculate finite direct and reflection paths between transmitting point and receiving point in a tunnel's central area.We do not need to divide entire computing space into Yee grids or allocate memory space for electric and magnetic field components in these grids anymore.So theoretically our RT-FDTD method requires much less memory space than FDTD method.
We can set the distance from the surface excitation source S RF1 to the tunnel wall to be 1/4 of a tunnel's width, and we can calculate surface excitation source SRF1 by RT method.We use surface excitation source S RF2 's z-directional coordinate position as a reference electric field strength; its value should not be too large-an integer between 1 and 4 is desirable.We consider surface excitation sources S RF3 's and S RF4 's sizes comprehensively based on required calculation for transceiver antenna's position from a tunnel wall and computer memory's size.Least ideal condition is when surface excitation sources and coincide with roof and floor respectively.In this case, our RT-FDTD method saves at least 1/4 memory space compared to traditional FDTD method.We can also reduce memory space required by this method by adjusting surface excitation sources S RF3 's and S RF4 's heights.
We calculate memory space and simulation time required by FDTD, UWR-FDTD and our RT-FDTD methods on same server for five electromagnetic wave frequencies: 580 MHz, 740 MHz, 900 MHz, 2100 MHz and 2600 MHz to calculate computing efficiency of these methods.Server CPU model is Intel(R) Core (TM) i7-7700K @ 4.20 GHz with 32 GB RAM.Simulated tunnel has length of 200 m, width of 4.8 m and height of 3.4 m; and spatial steps for Yee grid along x-axis, y-axis and z-axis are all 1/10 of an electromagnetic wave's wavelength.
We set unit window length to be 1/8 of a tunnel's length when we simulate three frequency electromagnetic waves-580 MHz, 740 MHz and 900 MHz-using UWR-FDTD method; this is best division window for UWR-FDTD method and most efficient calculation.If we divide unit window by highest efficiency for UWR-FDTD method when we simulate 2100 MHz and 2600 MHz frequency electromagnetic waves, required memory space exceeds 32 GB and we cannot calculate.We set unit window length to be 1/30 of a tunnel's length in simulation experiment; and we trade simulation time for storage space to calculate 2100 MHz frequency electromagnetic wave using UWR-FDTD method.
We set distance from face excitation source to tunnel wall to be 1/4 of a tunnel's width; distance from face excitation source to emission point to be 4 m; and height of face excitation source and tunnel to be 1/4 of a tunnel's height when we use our RT-FDTD method to calculate above five frequency electromagnetic waves.
Table 1 shows required memory space and simulation time for those three methods.We estimate required memory space for FDTD method's calculation of 2100 MHz and 2600 MHz frequency electromagnetic waves from a line fit to required memory space for FDTD method's calculation of 580 MHz, 740 MHz and 900 MHz frequency electromagnetic waves in Table 1.Required memory space and simulation time for FDTD method's calculation of these three frequency electromagnetic waves are what we need to run MATLAB programs.
As can be seen from the data in Table 1, compared with the traditional FDTD method, the performance of UWR-FDTD and RT-FDTD methods has been improved in terms of memory space and simulation time.However, UWR-FDTD method still requires a huge memory space to calculate higher frequency electromagnetic waves.It cannot calculate and simulate the propagation characteristics of 2600 MHz electromagnetic wave in the tunnel.We calculate the memory size of 580 MHz, 740 MHz and 900 MHz by FDTD method, and then we do the linear fitting of these memory values.According to the fitting data, the memory space required by FDTD method to calculate 2100 MHz and 2600 MHz electromagnetic wave is 72,400 MB and 125,400 MB, respectively.Because the computing space required to divide Yee grid is greatly reduced, the memory space required by RT-FDTD method to compute 2600 MHz electromagnetic wave is increased by 1800 MB compared with that required by 580 MHz electromagnetic wave, which is about 3 times.The simulation time is 16 times longer than that of 580 MHz electromagnetic wave.Under the same conditions, the memory space required by FDTD method to calculate 2600 MHz electromagnetic wave is 27 times that of 580 MHz electromagnetic wave.Because UWR-FDTD method only divides the unit window along the tunnel axis, but does not divide the tunnel section, the memory space required by the tunnel section alone exceeds the available memory space of the server when calculating the electromagnetic wave propagation characteristics of 2600 MHz.After we adopt the strategy of increasing simulation time in exchange for memory space, the simulation time required by the UWR-FDTD ( 18) www.nature.com/scientificreports/method to calculate this frequency is 23 times that of 900 MHz.We can complete the simulation calculation of 2600 MHz electromagnetic wave within 1.6 h using RT-FDTD method.

Conclusions
1.The RT-FDTD method effectively combines the high efficiency of the RT method and the precision of the FDTD method, guided by the law of conservation of energy.This approach addresses the limitation of the RT method, which is not suitable when transmitting and receiving antennas are located on the same side of the tunnel.Moreover, it significantly reduces the demand for external resources required by the FDTD method, thereby improving the computational efficiency compared to using the FDTD method alone.This innovative method offers a promising solution for enhancing wireless communication in complex environments such as underground mine.
The RT-FDTD method combines RT method's high efficiency and FDTD method's precision using law of conservation of energy.It solved the problem that the RT method is not suitable for transmitting and receiving antennas located on the same side of the tunnel at the same time.At the same time, it greatly reduced the demand for external resources of the FDTD method and improved computational efficiency of using the FDTD method alone.2. The RT-FDTD method consistently calculates the signal strength of different frequency electromagnetic waves in line with results obtained from both the FDTD and UWR-FDTD methods, as well as measured values presented in this paper.This consistency holds true regardless of the location of the transmitting and receiving antennas.Given that communication equipment antennas are typically installed on the tunnel wall or roof, the RT-FDTD method proves particularly suitable for planning antenna positions of communication base stations in tunnels.This makes it a valuable tool for optimizing wireless communication in such complex environments.3. The propagation characteristics of a 740 MHz electromagnetic wave were calculated using the FDTD, UWR-FDTD, and RT-FDTD methods in a simulated tunnel.The RT-FDTD method proved to be significantly more efficient, requiring only 1/15 of the simulation time needed by the FDTD method and one-third of the time required by the UWR-FDTD method.For 900 MHz electromagnetic wave, the RT-FDTD method costs only 1/22 of the simulation time needed by the FDTD method and one-third of the time required by the UWR-FDTD method.Due to the high memory demands of the FDTD and UWR-FDTD methods, these methods could not be used to calculate the propagation characteristics for a 2600 MHz electromagnetic wave.However, the RT-FDTD method was able to perform the same simulation in approximately 1.6 hours.These results highlight the efficiency and practicality of the RT-FDTD method in simulating wireless propagation characteristics in underground mine.

1 IFigure 2 .
Figure 2. Different image points of the emitter T.

Figure 3 .
Figure 3.The region model calculated by RT-FDTD.
, j and k are integers, representing the grid coordinate values of Yee grid in the x, y, z directions in the cuboid.The expressions of the four coefficients in the equation are C a = 2ε−σ �t 2ε+σ �t , C bx = 2�t 2ε�x+σ �t�x , C by = 2�t 2ε�y+σ �t�y , C bz = 2�t 2ε�z+σ �t�z .The FDTD iterative equation for solving the magnetic field of the Yee grid inside the cuboid are

Figure 4 .
Figure 4. Location parameters of the transmitting point and excitation surface source.

Figure 9 .
Figure 9.The section size of actual tunnel and equivalent rectangular tunnel.

Figure 10 .
Figure 10.The signal strength of 580 MHz electromagnetic wave.

Figure 11 .
Figure 11.The signal strength of 740 MHz electromagnetic wave.

Table 1 .
Memory space and simulation time required of three methods.