^{1,2}

^{1}

^{1}

^{1}

^{2}

Electromagnetic calculation plays an important role in both military and civic fields. Some methods and models proposed for calculation of electromagnetic wave propagation in a large range bring heavy burden in CPU computation and also require huge amount of memory. Using the GPU to accelerate computation and visualization can reduce the computational burden on the CPU. Based on forward ray-tracing method, a transmission particle model (TPM) for calculating electromagnetic field is presented to combine the particle method. The movement of a particle obeys the principle of the propagation of electromagnetic wave, and then the particle distribution density in space reflects the electromagnetic distribution status. The algorithm with particle transmission, movement, reflection, and diffraction is described in detail. Since the particles in TPM are completely independent, it is very suitable for the parallel computing based on GPU. Deduction verification of TPM with the electric dipole antenna as the transmission source is conducted to prove that the particle movement itself represents the variation of electromagnetic field intensity caused by diffusion. Finally, the simulation comparisons are made against the forward and backward ray-tracing methods. The simulation results verified the effectiveness of the proposed method.

Electromagnetic calculation plays an important role in both military and civic fields. It is realized that simulation of an electromagnetic environment is useful and can give much information in communication systems design and radar testing. It is of great significance to calculate the distribution status of the electromagnetic field.

Some methods and models are proposed for calculation of indoor and outdoor wave propagation, that is, methods based on the numerical solution of Maxwell’s equations such as finite-difference time-domain method (FDTD) [

To calculate electromagnetic wave in a large range, FDTD method brings heavy burden in calculation and also requires huge amount of memory [

Ray-tracing method is based on geometrical optics (GO) principle, which decides reflection or refraction through simulating the propagation route of rays [

Particle method is used to simulate the macrostate based on a huge amount of individual particle which has a process of generation, motion, and vanishing.

Based on forward ray-tracing method and particle method, a Transmission Particle Model for calculating electromagnetic distribution status in space is presented in this paper. Using the GPU to accelerate computation and visualization can reduce the computational burden on the CPU. The movement of a particle obeys the principle of the propagation of electromagnetic wave, and then the particle distribution density in space reflects the distribution status of electromagnetic distribution status. Since the particles in TPM are completely independent, it is very suitable for the parallel computing based on GPU. Updated GPU particle data is stored in display memory, so the data does not need to be transferred back to the memory when making visualization of the particle, and it can directly draw the particle movement after transmission.

In this paper, we present Transmission Particle Model and use GPU based on CUDA, realize particles transmit, updates, demise, and other operations, and achieve particle model computation and synchronization visualization.

The basic idea of TPM is that particles are generated and transmitted intermittently in the origin of electromagnetic field, and the electromagnetic distribution status in space is obtained through calculating the movement of each particle.

So far there are many methods for electromagnetic calculation, such as finite-difference time-domain method (FDTD), parabolic equation method (PE), and ray-tracing method (RT). To calculate electromagnetic wave in a large range, FDTD method brings heavy burden in calculation and also requires huge amount of memory. PE method performs well for complex terrain but has significant limitation in the calculation of propagation angle. Generally speaking, RT method is applicable to the electromagnetic calculation in a large range.

From the generation in the origin of electromagnetic field, the particle transmission density mainly depends on antenna’s angle factor

Take the origin of electromagnetic field as a vertex; partition the space with equal interval in all directions. The sampling intervals for

Then we obtain

Denote

The initial coordinates are

The movement of particle is dynamic; that is, with the elapse of time, the current properties, such as position coordinates, speed

In Figure

The distance between the particles.

Since

At time

Since

The movement speed of particle is

Correspondingly, we have

When electromagnetic particles come to the interface of two media, they reflect according to reflection theory, some of them may vanish owing to reflection, and the probability of vanishing is decided by reflection factor

The liveability after reflection is

Then the vanishing probability is

If the electric field intensity vector is parallel to incidence plane (parallel polarization), the reflection factor of uniform plane wave in [

If the electric field intensity vector is vertical to incidence plane (vertical polarization), the reflection factor of uniform plane wave in [

After reflection, some particles vanish. The transmission interval angle between survived particles and their adjacent particles has been changed to

In Geometry Optics, there are three kinds of diffractions, that is, edge diffraction, spire diffraction, and surface diffraction. Since the attenuation of spire diffraction and surface diffraction is very fast, we can ignore the two kinds of diffractions and only take edge diffraction into account.

When the particle diffraction occurs, the position of the particle needs to be modified according to the diffraction principle of electromagnetic wave. And further, the death probability of the particle is calculated based on the diffraction factor, which can be used to decide the death of the particle.

Given a distance

The liveability after diffraction is

Then the vanishing probability after diffraction is

Similar to reflection, the transmission interval angle between the survived particles after diffraction, that is,

Electric dipole is the most simple electromagnetic transmission antenna. This paper is focused on the analysis of the electromagnetic field generated by electric dipole.

In unbounded ideal media, electromagnetic wave propagates without reflection, refraction, scattering, and absorption; there is only energy loss caused by diffusion. In calculation of electromagnetic transmission in a large area, the average power density of close field

In the far-reaching electromagnetic field, the electric field vector and magnetic field vector can be expressed approximately in [

Poynting vector mean value of any point in the space is

The radiation field generated by electric dipole is not only relevant to the distance from field point to the source but also relevant to the

By (

The antenna’s angle factor of electric dipole is

According to function

Angle factor figure of electric dipole.

Consider the coordinate origin is the transmission source with sampling intervals

The distance range of these particles to the transmission source is

The distance between the transmission source and the particles transmitted at the previous time instant is

Since these particles are uniformly distributed in local area, then we have

The volume of the area is

Substituting (

Then the density of particles is

Considering (

Since

Further we have

On account of

In order to express Poynting vector mean value by the density of particles and since the antenna’s angle factor of electric dipole is

Then we can get

Substituting this into (

Consider

In recent years, due to the limit of the integrated circuit components, we encountered bottlenecks when we want to enhance the central processing unit (CPU) clock frequency on existing infrastructure. When upgrading the performance of a single processor core becomes difficult, multicore CPU was launched, and people are turning to the Graphics Processing Units (GPU). In the past GPU was primarily used for graphics rendering and display, but with the improvements of GPU performance and functional structure, GPU is not only used for graphics processing and rendering, but also starts to become a powerful assistant for CPU computing.

The current mainstream GPU has the following characteristic features: (1) supporting IEEE 32-bit floating-point arithmetic and output; (2) providing flexible programmability, dynamic branching, and loop control functions; (3) supporting multiple rendering, avoiding multiple data transfer between the CPU and GPU; and (4) supporting texture access function. It can store data to the texture easily and get index access. (5) It supports output data to the texture features, which can improve data output speed and reduce the overhead of data readback. Using the GPU to accelerate computation and rendering can reduce the computational burden on the CPU, free CPU from the heavy computing tasks, and meet the requirements of other computing.

As the inventor of the GPU, NVIDIA launched a common computing programming environment CUDA, completely shielding the underlying graphics API, which makes GPU programming almost transparent for developers and lays the foundation for wide use of GPU common computing.

The process for the calculation of electromagnetic distribution status based on CUDA with the particle transmission model can be described as follows (Figure

The process for the calculation of electromagnetic based on CUDA.

Transmit particles, initiate each of the particles, and put the transmitted particles to the particle list.

Update the properties of the particles in the list for the next time instant.

Judge if the particle is below the ground. If so, reflection or diffraction is needed. Next, the choice between reflection and diffraction is made based on whether the particle is on the edge of spire of the terrain.

If reflection is needed, modify the position of the particle according to the reflection principle of electromagnetic wave. And further, calculate the death probability of the particle based on the reflection factor and decide the death of the particle based on the vanishing probability. If diffraction is needed, the process is similar to that of particle reflection.

Judge if the particle is outside the boundary, and then eliminate the outside particles and the death particles from the list.

If the number of the listed particles is not stable (i.e., the particles within the calculation area have not reached equilibrium), then return to Steps

If the number of the particles in the list is stable, partition the space into grids and then count and output the number of particles within each grid.

Since the particle in Transmission Particle Model is completely independent, TPM is very suitable for the parallel computing based on CUDA. Updated GPU particle data is stored in display memory, so the data does not need to be transferred back to the memory when making visualization of the particle. It can directly draw the particle movement after transmission.

In order to achieve the results, we use CUDA to compute and draw directly on the display memory. As we need to transfer data between OpenGL and CUDA, we should establish buffer zones between these two groups of API for the position coordinates and speed. We can get the actual address of the device memory through the kernel function. Firstly, call cudaGraphicsMapResources() function to tell CUDA runtime mapping the OpenGL shared memory, and then call cudaGraphicsResourceGetMappedPointer() to request a pointer to the mapped resource. Thus, the OpenGL can share with CUDA buffer.

Electromagnetic particle model is based on ray-tracing method, but it is rather different from ray-tracing. In ray-tracing method, the space distances between transmission rays are equal; we calculate the initial density on different directions and need to calculate the variation of the density of electromagnetic field caused by diffusion in space. As to electromagnetic particle model, the ray is replaced by particles, the energy of each particle is a constant, the intensity of electromagnetic field is represented by the density of particles, the transmission density on different directions depends on the angle factor of the antenna, and the movement of particle accords with the attenuation of diffusion in space. Consequently, the particle movement itself represents the variation of electromagnetic field intensity caused by diffusion.

To calculate the electromagnetic field intensity at a point, forward ray-tracing algorithm adopts the method of receiving sphere to define the receiving radius and count the numbers of rays in the scope of the receiving sphere. But the definition of receiving sphere may bring in significant errors. Consider there is only direct transmission in an infinite space; one receiving sphere should receive only one ray. The receiving radius equals the distance between two rays around the point. If the receiving sphere is defined larger than the distance, there could be two rays in the receiving sphere; if the receiving sphere is defined smaller than the distance, there could be no rays in the receiving sphere. That is the reason for errors in the definition of the receiving sphere. The energy of each particle in electromagnetic particle model is less than that of a ray in ray-tracing method, thus yielding the smaller error. Of course, in this way the deficiency of significant prediction errors in forward ray-tracing method can be overcome.

Compared with backward ray-tracing, TPM can calculate the electromagnetic distribution status in the whole area at one time, but the backward ray-tracing needs to calculate each point in the area once again. In other words, TPM needs less calculation than backward ray-tracing to get electromagnetic distribution status.

To verify the effectiveness of the proposed TPM, we made a comparison between the proposed method and the ray-tracing method with a simplified terrain, and the transmission source is an electric dipole antenna in the middle of the area.

In order to compare TPM with ray-tracing method, we choose several points to calculate the propagation loss according to the ray-tracing method and the particle method, respectively. The height for calculation is 100 m. Because backward ray-tracing is commonly used, the result is calculated by backward ray-tracing as standard. The results are shown in Figure

Compare TPM with ray-tracing method.

Wireless InSite is complex electromagnetic environment predictive analysis simulation commercial software. The software is based on UTD/GTD theory and mainly uses ray-tracing method combined with the improved FDTD method to calculate the electromagnetic propagation. The calculation model includes urban canyon model, fast 3D city models, full three-dimensional model, the vertical plane model, the free-space model, urban canyon FDTD model, FDTD model with sliding windows, HATA model, COST-HATA model, and WI real-time models. By combining electricity field with the antenna pattern, it calculates the path loss, time of arrival, angle of arrival of the electromagnetic field, and the electromagnetic trend distribution.

To further verify the feasibility of the model of Transmission Particle Model, we calculated the electromagnetic distribution in the same terrain using particle transmission simulation model and Wireless InSite and displayed the electromagnetic attenuation trend extracting surfaces 5 meter far from ground. Figure

Visualization comparing TPM with Wireless InSite software (dBm).

Wireless InSite

Transmission Particle Model

The calculation time of TPM is less than the WI using RT method. The calculation time comparison is shown in Figure

Calculation time comparing TPM with Wireless InSite (RT).

When Wireless InSite software needs to calculate the electromagnetic field strength of a point, it uses the method of setting the receiver and calculates the received electromagnetic energy for the receiver. If we want to calculate the electromagnetic field distribution trend of the entire space, we need to provide intensive receiver in the space at equal intervals. Transmission Particle Model can calculate electromagnetic distribution trend in the entire space of the region. In the calculation of the electromagnetic distribution trend in the entire space, the calculation amount of Transmission Particle Model is less than that of Wireless InSite software.

Based on ray-tracing method and the particle method, this paper constructs TPM for electromagnetic field calculation. The details, including the initiation of the particles, updating the properties of the particles, the reflection and diffraction calculation, are also presented. Since the particles in TPM are completely independent, we used parallel computing based on GPU for the computation. Updated GPU particle data is stored in display memory, when making visualization of the particle, and it can directly draw the particle movement after transmission. Deduction verification of TPM is conducted based on the electromagnetic field generated by electric dipole. Finally, the simulation comparisons are made against the ray-tracing method and the commercial software WI. The simulation results verified the effectiveness of the proposed method. In summary, TPM is an effective tool for electromagnetic calculation, such as the electromagnetic distribution status in the space with complex terrain.

The authors declare that there is no conflict of interests regarding the publication of this paper.

This work is supported by Scientific Research Common Program of Beijing Municipal Commission of Education (KM201411232007), Open Research Project of the Beijing Key Laboratory of High Dynamic Navigation Technology (Grant no. HDN2014007), and the National Natural Science Foundation of China (Grant no. 61031001), Beijing Natural Science Foundation (4142017).