
We discuss coupling conditions for the p-system in case of a transition from supersonic states to subsonic states. A single junction with adjacent pipes is considered where on each pipe the gas flow is governed by a general p-system. By extending the notion of demand and supply known from traffic flow analysis we obtain a constructive existence result of solutions compatible with the introduced conditions.
Citation: Martin Gugat, Michael Herty, Siegfried Müller. Coupling conditions for the transition from supersonic to subsonic fluid states[J]. Networks and Heterogeneous Media, 2017, 12(3): 371-380. doi: 10.3934/nhm.2017016
[1] | Joshua Holcomb, Nicholas Spellmon, Yingxue Zhang, Maysaa Doughan, Chunying Li, Zhe Yang . Protein crystallization: Eluding the bottleneck of X-ray crystallography. AIMS Biophysics, 2017, 4(4): 557-575. doi: 10.3934/biophy.2017.4.557 |
[2] | Domenico Lombardo . Scientific advance in biomembranes and biomimetic membranes of biophysical interest. AIMS Biophysics, 2022, 9(4): 341-345. doi: 10.3934/biophy.2022028 |
[3] | Mostean Bahreinipour, Hajar Zarei, Fariba Dashtestani, Jamal Rashidiani, Khadijeh Eskandari, Seyed Ali Moussavi Zarandi, Susan Kabudanian Ardestani, Hiroshi Watabe . Radioprotective effect of nanoceria and magnetic flower-like iron oxide microparticles on gamma radiation-induced damage in BSA protein. AIMS Biophysics, 2021, 8(2): 124-142. doi: 10.3934/biophy.2021010 |
[4] | Wei Zhang, Sheng Cao, Jessica L. Martin, Joachim D. Mueller, Louis M. Mansky . Morphology and ultrastructure of retrovirus particles. AIMS Biophysics, 2015, 2(3): 343-369. doi: 10.3934/biophy.2015.3.343 |
[5] | Jacques Fantini, Francisco J. Barrantes . How membrane lipids control the 3D structure and function of receptors. AIMS Biophysics, 2018, 5(1): 22-35. doi: 10.3934/biophy.2018.1.22 |
[6] | Missag Hagop Parseghian . What is the role of histone H1 heterogeneity? A functional model emerges from a 50 year mystery. AIMS Biophysics, 2015, 2(4): 724-772. doi: 10.3934/biophy.2015.4.724 |
[7] | Shivakumar Keerthikumar . A catalogue of human secreted proteins and its implications. AIMS Biophysics, 2016, 3(4): 563-570. doi: 10.3934/biophy.2016.4.563 |
[8] | Nicola Gaetano Gatta, Gaetano Cammarota, Vittorio Gentile . Possible roles of transglutaminases in molecular mechanisms responsible for human neurodegenerative diseases. AIMS Biophysics, 2016, 3(4): 529-545. doi: 10.3934/biophy.2016.4.529 |
[9] | Jonathan Trauth, Johannes Scheffer, Sophia Hasenjäger, Christof Taxis . Strategies to investigate protein turnover with fluorescent protein reporters in eukaryotic organisms. AIMS Biophysics, 2020, 7(2): 90-118. doi: 10.3934/biophy.2020008 |
[10] | Adam Redzej, Gabriel Waksman, Elena V Orlova . Structural studies of T4S systems by electron microscopy. AIMS Biophysics, 2015, 2(2): 184-199. doi: 10.3934/biophy.2015.2.184 |
We discuss coupling conditions for the p-system in case of a transition from supersonic states to subsonic states. A single junction with adjacent pipes is considered where on each pipe the gas flow is governed by a general p-system. By extending the notion of demand and supply known from traffic flow analysis we obtain a constructive existence result of solutions compatible with the introduced conditions.
The successful application of wireless sensor networks in monitoring, spectrum management, distributed detection, agriculture, and other fields is very important for a country's development [1,2,3]. As the underlying supporting technology of the Agricultural Internet of Things (AIoT), wireless sensor networks are widely deployed in various agricultural applications [4,5]. Under the new development direction of world agriculture, precision agriculture is gradually replacing traditional agriculture [6]. Soil moisture wireless sensor networks (SMWSNs) play a critical role in enabling precise irrigation of farmland, and SMWSNs are an essential tool for agricultural production [7,8]. SMWSNs have the characteristics of low cost, flexible structure, self-organizing network, and dynamic topology, which can accurately collect and effectively transmit soil moisture data and subsequently control the terminal to accurately adjust the amount of farmland irrigation [9,10]. By deploying SMWSNs in field conditions, farmers can promptly make more accurate agricultural production decisions based on monitoring the soil moisture information collected by sensor nodes in real time [11,12,13]. Figure 1 shows the sensor node information transmission. Task allocation is a difficult issue in wireless sensor networks. In wireless sensors, it is often necessary to perform various tasks, such as computing tasks, perception tasks, and communication tasks [14]. Task allocation reflects that in the dynamic environment of the network with limited energy, the tasks in the system are effectively allocated first, and then the specified tasks are scheduled to be executed on suitable nodes to achieve optimal performance [15]. Optimizing the task allocation of SMWSNs is of great significance for rationally allocating network resources to better complete information acquisition and precision irrigation [16]. SMWSNs are broadly applied in crop irrigation, though they also face the problems of a complex external environment and limited network resources. These problems will lead to untimely information collection and transmission, reducing the overall performance of the network, and thereby affecting agricultural production.
The problem of wireless sensor task allocation optimization has always been a hot topic discussed by researchers. Scholars have attempted numerous analyses and research on heuristic algorithms to address this NP-hard problem [17]. For example, Issac et al. [18] proposed a method using particle swarm optimization (PSO) to solve the multi-objective task allocation problem, which is more efficient than the task allocation of traditional PSO. The generation of the artificial bee colony algorithm (ABC) is inspired by the foraging activities of bees to seek out the optimal solution to the issue [19], where different activities are carried out according to their respective division of labor to achieve information sharing and communication. Hou et al. [20] designed an improved genetic algorithm (IGA) for the shortcomings of GA, which optimizes the task allocation problem by producing offspring to inherit the best characteristics of the previous generation. Shin et al. [21] designed an improved ant colony optimization algorithm (IACO) inspired by ant foraging behavior in nature. Gao et al. [22] proposed an improved firefly algorithm (IFA) based on simulated annealing and adaptive mutations to resolve task allocation problems, which have a better optimization ability and higher solution accuracy compared with FA. Although these heuristics can improve the efficiency of task allocation to some extent, they have the problem of low availability in evolutionary stagnation and local optimization. Applying it to the task problem of SMWSNs will lead to premature convergence and reduce network revenue. Chen et al. proposed an algorithm based on deep learning and developed a new activation function, the parameter-free adaptive Swish (PASwish), which can perform a more flexible nonlinear transformation of different input data; however, this algorithm is of a high complexity [23]. This paper mainly proposes DEAEBOA to obtain an improved higher network benefit and finally obtain an efficient task allocation scheme. The superiority of DEAEBOA is reflected by assigning different tasks to either different or the same number of nodes, saving sensor node resources to the greatest extent, and finally maximizing network benefits from task allocation.
In solving the problem of revenue maximization, DEAEBOA has an outstanding performance in SMWSNs, vastly heightening the overall benefit of the network. This paper has the following main contributions:
1) First, a novel DEAEBOA is proposed. Compared with the aforementioned studies, the algorithm can greatly enhance the efficiency of task allocation for SMWSNs. Meanwhile, it avoids evolutionary stagnation and significantly accelerates the convergence speed, which improves the performance of SMWSNs and effectively saves network resources; on the other hand, the method realizes the precise adjustment of farmland irrigation amount within the actual environment.
2) Second, a new adaptive operator is designed. Compared with the aforementioned studies, the algorithm can immensely improve the accuracy and performance. Moreover, a novel elite operator and differential evolution strategy are designed, which can efficaciously enhance the global search capacity and avert falling into local optimization.
3) Thirdly, a new mathematical model for task allocation optimization of SMWSNs is established. Compared with recent model studies, the algorithm can consider and optimize several fundamental factors affecting the performance of SMWSNs, including node sensing distance, number of nodes, and information acquisition advantages. In addition, the model employs a new objective function to balance the relationship between these factors and serve as an evaluation criterion for network benefits.
The general framework of this article is as follows. Section 2 covers the assignment of tasks related to SMWSNs. Section 3 details the shared task assignment model. Section 4 recommends DEAEBOA address task allocation in SMWSNs. Subsequently, the validity of DEAEBOA's results through simulation experiments is discussed in Section 5. Section 6 returns to the conclusion.
Wireless sensor networks have good stability, which is important for the time delay and security of signal transmission [24,25]. Heuristics are widely used in agricultural wireless sensor networks. In order to provide a fuller introduction to the work related to the assignment of tasks related to SMWSNs, this chapter will be discussed in two parts. First, the necessity of heuristic algorithms in the task allocation of agricultural wireless sensors is introduced. Finally, the application of heuristic algorithms in wireless sensors is probed. In the meantime, the innovative heuristic algorithm in this paper is elaborated.
In the period of agricultural modernization, the task distribution of wireless sensors in agriculture has attracted the attention of many scholars [26]. The task assignment problem is a typical multi-objective math issue, the purpose of which is to further optimize the performance of the system by designing a reasonable and effective allocation scheme to optimize the execution effect of all tasks [27]. According to the environment and resource constraints, a mathematical model of task allocation is established, and an objective function is also required to use efficient algorithms to solve the optimal task allocation scheme [28]. The methods for solving task allocation problems are mainly mathematical programming methods and heuristic algorithms. The mathematical programming method has evident advantages in solving constraint problems with a small number of constraints and variables. However, the agricultural production environment is complex and changeable, resulting in an increase in the amount of system network calculation, a low solution efficiency, and an inability to adapt to the development environment of smart agriculture [29]. Therefore, traditional mathematical methods to solve task allocation problems in SMWSNs are no longer desirable. Due to the advantages of simple and easy implementation, fast calculation speeds, and strong compatibility, the heuristic algorithm is widely used in agricultural wireless sensors. Compared with mathematical programming, heuristic algorithms can handle task allocation problems in SMWSNs well [30,31].
In recent research on task assignment in wireless sensor networks, many scholars have proposed their own solutions. Okhovvat et al. [32] proposed an analytical method based on the queuing theory to calculate the optimal assignment rate from task to actor node. To solve the task assignment problem in virtual wireless sensor networks, Raee et al. [33] designed a scalable heuristic algorithm that minimizes total energy consumption. Kori et al. [34] proposed a classified regression tree machine learning algorithm (CART) for resource allocation schemes to improve the performance and extend the service life of wireless sensor networks by adopting intelligent resource management schemes. Baniabdelghany et al. [35] designed an offline discrete particle swarm optimization (DPSO-TA) algorithm for reliable task assignment in wireless sensor networks of the Internet of Things (IoT), which effectively reduces network energy consumption. However, the algorithms mentioned above have some problems, such as large computation and slow convergence.
To more intuitively reflect the superiority of DEAEBOA, the algorithm proposed by other researchers mentioned below will be compared in three key points, as shown in Table 1.
Algorithm | Rate of convergence | Local optimization | Computational complexity |
DEAEBOA | fast | \ | low |
IACGA | fast | √ | low |
DPSO | fast | √ | low |
EIACO | fast | √ | high |
SNO | fast | \ | high |
MOA | fast | \ | high |
BOA | fast | √ | high |
The earliest heuristic algorithms proposed were GA, PSO, and ACO, in which these bionic algorithms attracted the favor of researchers once proposed. In-depth research based on these traditional heuristics led to the proposal of many new heuristics to solve task allocation problems. Li et al. [36] proposed an improved adaptive cloning genetic algorithm (IACGA) and achieved good results, though it was prone to evolutionary stagnation. Baniabdelghany et al. [35] proposed a discrete particle swarm optimization (DPSO) for distributing tasks between sensor nodes. This method saves time and cost and improves accuracy, but it is effortless to fall into local optimization. Xu et al. [37] designed an elite immune ant colony optimization (EIACO) to dispose task allocation problems. The algorithm has a better task allocation efficiency and higher convergence velocity, though the algorithm complexity is high and the performance is poor.
Sensor nodes can collaboratively perform complex tasks but are often constrained by energy and capacity. In this case, it is critical to find a task allocation scheme that significantly improves the efficiency of the network. Niccolai et al. [38] designed a new evolutionary algorithm for social network optimization (SNO) to deal with the task allocation trouble, which effectively improves the task assignment efficiency under the condition of energy and energy constraints; however, the algorithm has the matter of high computational complexity. Weikert et al. [39] proposed a multi-objective optimization algorithm (MOA) for the dynamic task allocation problem of the IoT. Compared with the existing single-objective algorithm, this method reduces the network latency and improves the network performance. However, with additional tasks and nodes, the computational complexity of the method increases too quickly. Arora et al. [40] used a butterfly optimization algorithm (BOA) to work out problems in wireless sensors, and achieved relatively good performance. However, one problem with this algorithm is that it tends to converge prematurely.
Due to aforementioned methods for solving wireless sensor task allocation, there are problems such as high computational complexity and an ease of falling into local optimization, resulting in a high time delay, waste of network resources, and low task allocation efficiency. Therefore, it is necessary to design a strategy that can avoid premature convergence and has a low computational complexity. This paper proposes a novel DEAEBOA. DEAEBOA has the characteristics of fast convergence speed and low complexity. In SMWSNs, there is a certain relationship between the collection advantages of soil moisture information collection tasks and network benefits, and DEAEBOA can obtain a task allocation scheme that maximizes network benefits. Compared with other heuristic methods, it is easier to jump out of local optimization, which can obviously elevate the task allocation efficiency of SMWSNs and save network resources. In addition, a unique task allocation model for SMWSNs is proposed, which has the preponderances of low power consumption and simple processing. DEAEBOA optimizes the SMWSNs model to improve task allocation efficiency and network performance. A simplified diagram of DEAEBOA acting on SMWSNs is shown in Figure 2.
In SMWSNs, sensor nodes are deployed in farmland areas to monitor soil moisture conditions. In agricultural production, it is more necessary to pay attention to the synergy and economic benefit between sensor nodes. Therefore, to enhance the efficiency of SMWSNs and to guarantee the stable performance of sensor networks, an excellent distribution scheme should be adopted to allocate a wide range of tasks to a limited number of sensor nodes. In addition, the symbol descriptions that appear in this section are listed in Table 2.
Notation | Statement | |
m | Number of soil moisture information collection tasks | |
n | Number of sensor nodes | |
B | Soil moisture information collection advantage matrix | |
U | Economic benefit matrix | |
D | Solution in integer form | |
Z | Solution in binary form | |
Zb | Soil moisture information collection advantage matrix specific solutions | |
F | Soil moisture information collection benefit matrix | |
benefit | Total soil moisture information collection benefit |
In this paper, an efficient task allocation model based on SMWSNs is designed. In SMWSNs, the model can be simply described as assigning m soil moisture information collection tasks to n soil moisture sensor nodes according to the economic benefit and information collection advantages. Visually, Figure 3 depicts the task assignment model.
Due to the difference in the distance between each sensor node and each farmland area, the advantages of different sensor nodes collecting soil moisture information in different areas may also be different. The predominance matrix is shown in (1):
B=[b1,1b1,2⋯b1,m−1b1,mb2,1b2,2⋯b2,m−1b2,m⋯⋯bx,y⋯⋯bn−1,1bn−1,2⋯bn−1,m−1bn−1,mbn,1bn,2⋯bn,m−1bn,m](x∈[1,n],y∈[1,m]), | (1) |
where bxy represents the soil moisture information collection advantage of the xth humidity sensor node over the yth farmland area, which is equal to the reciprocal of the distance between the yth farmland area and the xth humidity sensor node. The closer the distance, the better the soil moisture information collected. At the same time, the larger the bxy, the closer the xth humidity sensor node is to the yth area, and the greater the soil moisture information collection advantage of the xth humidity sensor node for the yth farmland area. Suppose that the advantages of farmland soil moisture information collection are evaluated and ranked, and the economic benefit is evaluated and described as (2):
U=[u1u2…uy…um−1um](y∈[1,m]), | (2) |
where uy represents the economic benefit of the jth farmland area. Then, the allocation scheme can be represented as
D=[d1d2…dk…dm−1dm](k∈[1,m],dk∈[1,n]). | (3) |
In (3), the number in D indicates which sensor node is assigned to, and D is used to store the allocation result. Section 4 describes coding schemes and particle initialization. In addition, in order to facilitate subsequent operations, the allocation scheme in binary form is given as follows:
Z=[z1,1z1,2z1,3…z1,mz2,1z2,2z2,3…z2,m………zi,j…zn,1zn,2zn,3…zn,m](zi,j∈{0,1},i∈[1,n],j∈[1,m]). | (4) |
In (4), Z represents a randomly generated binary task allocation scheme. The quantity of rows in Z corresponds to the quantity of soil moisture sensor nodes, and the quantity of columns corresponds to the quantity of soil moisture information collection tasks. zi,j=1 indicates that the ith sensor node is responsible for the jth soil moisture information acquisition task; otherwise, the opposite is true.
For example, let's take a model with 5 soil moisture information collection tasks and 3 soil moisture sensor nodes. As shown in Figure 4, the soil moisture information collection task allocation model.
B=[b11b12b13b14b15b21b22b23b24b25b31b32b33b34b35] | (5) |
U=[u1u2u3u4u5] | (6) |
In (5) and (6), the number of columns in B and U corresponds to the number of soil moisture information collection tasks, while the quantity of rows in B corresponds to the quantity of soil moisture sensors. For clarity, the allocation scheme for the example model is shown below:
D=[23123], | (7) |
Z=[001001001001001], | (8) |
where D and Z represent the allocation scheme in integer form and in binary form, respectively.
Zb=B⋅Z | (9) |
Zb=[00b1300b2100b2400b3200b35] | (10) |
In (9), Zb represents the advantage matrix corresponding to the randomly generated soil moisture information collection task allocation scheme Z. In addition, the value in Zb represents the advantage of jth soil moisture information collection task being allocated to the ith soil moisture sensor node.
F=Zb×WT | (11) |
Finally, the soil moisture information acquisition benefit matrix F is calculated in (11).
benefit=∑ni=1fi | (12) |
In (12), fi represents the value of ith in the benefit matrix F, and the total soil moisture information collection benefit of the entire system can be obtained by calculating the sum of all fi.
In SMWSNs, the task allocation scheme not only realizes the precise irrigation of farmland, but also increases the actual economic value of farmland and enormously elevates the efficiency of the system. In this way, each sensor node has a certain benefit in completing the corresponding task, and the sum of the benefit of all nodes is the total benefit of the whole system. Due to the influence of environmental factors and the complexity of practical problems in agricultural production, it is also necessary to consider the performance of sensor nodes and the load difference of the task assigned.
In order to make SMWSNs have a better network benefit, a task allocation method based on DEAEBOA is proposed. This method possesses a sturdy global search ability and can stably converge. At the same time, due to the innovative introduction of adaptive operators, differential evolution, and elite strategies, the method can effectively avoid falling into local optimum and premature convergence. Moreover, these strategies enable DEAEBOA to obtain the optimal information collection task assignment scheme and significantly improve the network benefit of SMWSNs. DEAEBOA is an evolutionary algorithm that mimics how butterflies search for food. In DEAEBOA, the optimization process of the problem is seen as a group of scented butterflies attracting each other. Each butterfly, also known as a particle, corresponds to a potential solution to a problem. Then, finding the particle with the strongest fragrance is the optimal solution to the problem. Initially, this group of particles emits different amounts of fragrance, each either moving randomly or towards the particles that emit more fragrance. Over the course of each iteration, each particle is continuously adjusted according to its optimal position and the optimal position of the population, recalculating the particle position and fragrance intensity to achieve the best results.
The process of DEAEBOA can be described as coding schemes and particle initialization, fitness calculations, particle motion, differential evolution strategies, adaptive operators, elite operators, and termination conditions.
Encoding is the first step in DEAEBOA. There are two commonly used encoding methods: one is binary encoding and the other is integer encoding. In order to optimize the network benefits of SMWSNs as much as possible, as well as to perform a wider range of searches, this article uses integer encoding. Assuming that there are x particles, n soil moisture sensors, and m tasks in the initial group, it can be described as (13). D1,D2⋯Dx are the x particles in the initial group, each of which is a vector of length m, representing a potential solution. The values in the vector represent integers between 1 and n, which is consistent with Eq (3). First, the initial array is randomly generated, rounding the initial group to limit the values to a range of 1 to n. Each time the group is updated, a rounding operation is made to limit the value to a range of 1 to n.
group=[D1D2D3⋮Dx]=[d1,1d1,2d1,3⋯d1,md2,1d2,2d2,3⋯d2,md3,1d3,2d3,3⋯d3,m⋯⋯⋯du,v⋯dx,1dx,2dx,3⋯dx,m] | (13) |
In general, a large factor influencing the convergence of DEAEBOA is the objective function. Each particle has its corresponding value and a corresponding degree of fitness. In DEAEBOA, fitness is assessed using the Eq (12). According to the above Eqs (11) and (12), it can be seen that in order to achieve the maximum benefit of the network, the benefit needs to be as large as possible, that is, Zb needs to be as large as possible.
With each particle iteration, its fitness is compared to its historical optimal solution. If the current particle fitness is high, the particle historical optimal solution is updated, which is the current particle fitness. Otherwise, leave it as is. Second, the current particle fitness is compared with the global optimal solution of this group. If the current particle fitness is high, the global optimal solution of the particle is updated and replaces the current optimal solution position; otherwise, it will not change.
In particle initialization, the initial position of each particle is different. Each particle emits odors of different intensities depending on where it is located. The odor formula is described in (14). Depending on the effect of odor intensity, each particle either moves randomly or towards particles that emit more odors. When a particle moves towards a particle with a strong odor concentration, it is called a global search. The equation is described in (15). Instead, the particles move randomly, which is a process called local search. The equation is described in (16). Conversion probability p to control the global and local search process, p∈[0,1].
f=ctIa | (14) |
xit+1=xit+(r2×g∗−xit)×fi | (15) |
xit+1=xit+(r2×xjt−xkt)×fi | (16) |
In Eq (14), f is the amount of fragrance perception and c is the sensory form. I represents the stimulus intensity, that is, the fitness function value, and a stands for fragrance, a,c∈[0,1]. ct represents the form of sensation in the t generation. In Eq (15), xit+1 represents the position of the ith individual in the t+1 generation, r∈[0,1], g∗ represents the global optimal solution, fi represents the scent perception of the ith individual. In Eq (16), xit,xjt,xkt indicates the position of the I, j, and k individuals in the t generation, respectively.
To be more in line with the laws of nature, the value of c will change accordingly during each iteration of the particle. The specific operation is as follows:
ct+1=ct+(b/(ct×Ngen)) | (17) |
where ct+1 represents the form of sensation in the t generation, b is a constant, and Ngen is the maximum number of iterations.
The differential evolution strategy consists of three parts: mutation, crossover, and selection operations. Finally, the particles with the highest adaptability are retained. The specific operation is as follows.
In each iteration, the mutation process is to select three different particles in the group and generate new particles through Eq (18):
Vi,G+1=Xr3,G+F∗(Xr1,G−Xr2,G) | (18) |
where i, r1, r2, and r3 belong to positive integers from 1 to n, i≠r1≠r2≠r3, and F is the variation probability, F∈[0,1].
We can select according to a certain probability, and the new particle produced in Eq (18) crosses the current particle; otherwise, leave it as it is:
Uj,i,G+1={vj,i,G+1ifjrand≤CRxj,i,Gotherwise | (19) |
where CR is the crossover probability factor and jrand is a random component, ensuring that at least a one-dimensional component of the new particle after the crossing is provided by the mutant particles.
Based on the fitness function value, the better of the new particles and target particles that have undergone mutation and cross-generation is selected as the next generation.
Xi,G+1={Ui,G+1iff(Ui,G+1)≤f(Xi,G)Xi,Gotherwise | (20) |
Adaptive operators have obvious advantages in dealing with complex optimization problems. In previous classical evolutionary algorithms, scholars typically set constants that remain constant throughout the lifetime of the algorithm. One advantage of this is that it can decrease the challenge of algorithm design, making the algorithm less complex; however, the algorithm tends to converge prematurely and fails to achieve the expected results. The adaptive operator designed in this paper is used to control a single parameter. We set the conversion probability in the algorithm to P. The adaptive conversion probability p controls whether particles are searched either globally or locally. If the setting is too large, the particles will not conform to the natural evolution law and the algorithm will be unstable. If the setting is too small, the evolution speed of the algorithm will be affected. Hence, we adopt an adaptive strategy for P to enhance the performance of the algorithm. For the case where P is close to zero, resulting in premature convergence of the algorithm, we design Eq (21) to quantify the conversion probability under nonlinear changes in individual fitness values:
P(i)=Pmin+Q∗(Pmax−Pmin)∗fmax−f(i)fmax−fmin | (21) |
where Pmax and Pmin are the upper limit of 1 and the lower limit of 0.7, respectively, Q is set to 0.92, fmax represents the value with the greatest fitness of all particles, fmin is the minimum value of all particles, and f (i) represents the current fitness of the ith particle.
In order to shorten the search time, this paper designs an elite operator that proposes that the best particles of the previous generation replace the worst particles of the next generation. This makes the algorithm converge faster.
During each iteration, DEAEBOA checks whether the termination conditions are met. The termination condition of the algorithm is to reach the maximum number of iterations.
In the algorithm flow shown in Figure 5, we will describe the algorithm flow in detail.
Step 1. Initialize the parameters in DEAEBOA, set the maximum quantity of iterations, cross probability, variation probability, fragrance parameter, sensory form parameter, maximum and minimum transition probability, and randomly generate the initial group. According to Eq (1), the sensor node information acquisition advantage matrix is generated. The task urgency matrix is generated by Eq (2).
Step 2. Through Eqs (11) and (12), the fitness of the old group is obtained, and the historical optimal value of the particle and the global optimal value of the group are recorded.
Step 3. Use Eq (14) to calculate the perceived amount of fragrance, and use Eqs (15) and (16) to search either globally or locally. Renew the sensory form according to Eq (17).
Step 4. To be more in line with the laws of nature, a differential evolution operation is carried out.
Step 5. Equations (11) and (12) are used to assess the fitness of a new swarm and to record the historical optimal of the particle and the global optimal of the swarm.
Step 6. Use an adaptive strategy for the conversion probability according to Eq (21).
Step 7. Execute elite strategies.
Step 8. Repeat 3 through 7 until the termination condition is met, that is, the maximum quantity of iterations is reached.
The difficulty of this analysis method lies in DEAEBOA. When the position of the particle is updated, the adaptive conversion probability is designed, which changes with the change of the number of iterations, which increases the difficulty of the algorithm, and conforms to the random law of natural change. In order to improve the convergence speed of the algorithm, the elite strategy is adopted to ensure the stability of the whole simulation experiment. Finally, the variation and selection operation of the differential evolution can avoid the local optimal situation, but increase the difficulty of the algorithm. The organic combination of the three strategies improves the performance of SMWSNs.
Suppose the population size is N, the number of iterations is MAX_GEN, and the dimension is D. Based on the operation rules of the time complexity symbol O, it can be deduced that the time complexity of DEAEBOA randomly initializing the population is O (N·D), and the time complexity of finding the local optimal is O (N·D). The time complexity of global position update change using differential evolution adaptive elite operator strategy is O (MaxIter·N·D), and the total time complexity of DEAEBOA is O (MaxIt·N·D). Therefore, the time complexity of DEAEBOA is the same as that of BOA, and the computation work will not be increased.
In this section, DEAEBOA's performance in soil moisture information collection task allocation will be tested through a simulation to achieve the precise irrigation of farmland. In simulation, DEAEBOA is compared with PSO, GA, and BOA methods in different situations, and the average of 100 experiments is used as the result of the four algorithms. The task allocation model with five different algorithms is simulated on a computer with an Intel(R) Core (TM) i7-7700H CPU, 8GB RAM, and Windows 10 operating system by using MATLAB R2022a.
During the simulation experiment, for the accuracy of the simulation, and to compare the performance of DEAEBOA with BOA, GA, PSO, and BWO, we set the population size of the five algorithms to 100, and the maximum quantity of iterations in this paper was set to 800 and 150. In DEAEBOA, we set the adaptive conversion probability range to 0.7–0.9, the mutation probability to 0.95, and the crossover probability to 0.95. In BOA, we set the conversion probability to 0.9. In GA, we set the mutation probability to 0.4. In PSO, we set the inertia weight to 0.5, and both the self-learning factor and the group-learning factor to 0.55. For DEAEBOA, BOA, GA, and PSO, Table 3 shows the main parameters of the four algorithms.
DEAEBOA | Population size | Adaptive conversion probability | Mutation probability | Crossover probability | Fragrance | Sensory form |
100 | 0.7–0.9 | 0.95 | 0.95 | 0.1 | 0.01 | |
BOA | Population size | Conversion probability | Fragrance | Sensory form | ||
100 | 0.95 | 0.1 | 0.01 | |||
PSO | Population size | Inertia weight | Self-learning factor | Group-learning factor | ||
100 | 0.5 | 0.6 | 0.6 | |||
GA | Population size | Mutation probability | ||||
100 | 0.1 | |||||
BWO | Population size | Equilibrium factor | ||||
100 | 0.5 |
To make the simulation results more scientific and researchable, in the discussion section of the experimental results, the data will be analyzed from multiple angles, displayed by simulating evolution curves, line charts, bar charts, pie charts and tables, and at the same time the charts are analyzed and discussed.
Figure 6(a)–(d) visually show the simulation results of DEAEBOA, PSO, GA, BOA, and BWO with 8, 16, 30, and 50 sensor nodes and 15, 30, 60, and 120 tasks, respectively. Overall, DEAEBOA has an improved optimized performance compared to BOA, PSO, GA, and BWO in 4 different situations. Figure 6(a)–(d) show that DEAEBOA has begun to converge by the time the algorithm runs to the 200th generation, at which point DEAEBOA's network benefit are much greater than PSO, GA, BOA, and BWO. At the same time, DEAEBOA shows more obvious advantages in the case of the increasing number of tasks and nodes, and DEAEBOA's benefit also increases more than BOA, PSO, GA, and BWO. It is clear that BOA, PSO, GA, and BWO exhibit a locally optimal state after 100 iterations. Due to the creative addition of adaptive strategies, elite operators, and differential evolution, DEAEBOA is easier to jump out of the local optimal and obtain a better benefit value. Overall, DEAEBOA can deal with the issue of information collection task allocation in SMWSNs when the conditions are the same.
The histogram in Figure 7 shows that when the sensor nodes are 8, the number of tasks is 15, 25, 40, and 50, respectively. The network benefit gap between the five algorithms can be more visualized using the histogram view, and Figure 7(a)–(d) are the result of 150 iterations. It is obvious that DEAEBOA consistently performs the best under these four conditions. When the quantity of nodes is constant, a change in the quantity of tasks will make it continue to grow, the network income becomes higher and higher, and DEAEBOA always performs the best. In Figure 7(d), the network benefit of DEAEBOA is much larger than that of BOA, PSO, GA, and BWO. For agricultural production problems, the information collection effect and economic value can be effectively improved.
Table 4 shows the network benefit for DEAEBOA, PSO, GA, BOA, and BWO with different numbers of nodes and tasks.
Number of nodes and tasks | DEAEBOA | PSO | GA | BOA | BWO |
15 tasks and 8 nodes | 6.4452 | 6.2580 | 6.0569 | 5.8687 | 5.8879 |
25 tasks and 8 nodes | 11.1075 | 10.6126 | 10.0674 | 9.9411 | 9.8579 |
40 tasks and 8 nodes | 16.9841 | 15.9036 | 15.3488 | 14.9900 | 14.9503 |
50 tasks and 8 nodes | 20.5443 | 19.136 | 18.9283 | 18.2450 | 18.1966 |
30 tasks and 16 nodes | 12.6718 | 11.9088 | 11.6171 | 11.2362 | 11.1878 |
60 tasks and 30 nodes | 23.3405 | 22.2220 | 22.0922 | 21.1183 | 21.6495 |
120 tasks and 50 nodes | 45.1898 | 42.3056 | 42.1824 | 40.7369 | 41.4927 |
Number of nodes and tasks | PSO | GA | BOA | BWO |
15 tasks and 8 nodes | 3.00% | 6.41% | 9.82% | 9.47% |
25 tasks and 8 nodes | 4.66% | 10.33% | 11.73% | 12.67% |
30 tasks and 16 nodes | 6.40% | 9.08% | 12.78% | 13.26% |
40 tasks and 8 nodes | 6.80% | 10.65% | 13.30% | 13.60% |
50 tasks and 8 nodes | 7.36% | 8.53% | 12.60% | 13.00% |
60 tasks and 30 nodes | 5.03% | 5.65% | 10.52% | 7.81% |
120 tasks and 50 nodes | 6.81% | 7.13% | 10.93% | 9.00% |
Figure 8 shows the benefit growth percentages of the five algorithms corresponding to different task numbers when the number of nodes is 8, 16, 30, and 50, respectively. Data are taken from Table 4, and Figures 6 and 7. As can be seen in Figure 8, DEAEBOA's network benefit growth percentage is greater than that of BOA, PSO, GA, and BWO for seven different tasks. It can be seen that the innovative design of adaptive operators, elite strategies, and differential evolution strategies can jump out of the local optimum, so as to achieve higher network returns and improve the task allocation efficiency of SMWSNs.
Figure 9 shows a line chart of the network benefit of the growth of the five algorithms when the quantity of nodes is 8, as the quantity of tasks augments. When the number of tasks is 25 and 40, the network benefit growth rate is highest in DEAEBOA, up to 16%.
Figure 10 displays a line chart of the network benefit growth of the five algorithms under different task numbers and node numbers. In the case shown in Figure 10, the benefit growth percentage of DEAEBOA is invariably greater than that of BOA, PSO, GA, and BWO. In summary, DEAEBOA is very available in resolving the task assignment problem.
Figure 11 shows the percentage increase in the benefit of DEAEBOA compared to BOA, PSO, GA, and BWO. We can draw according to the data generated under the number of nodes and tasks in Table 4. Obviously, in either case, DEAEBOA is always the best choice for SMWSNs.
Table 4 shows the network benefit improvement rate of DEAEBOA compared to PSO, GA, BOA, and BWO.
Figure 12 shows the average network benefit improvement rate of DEAEBOA compared to the other four algorithms when the quantity of nodes is 8 and the quantity of tasks is 15, 25, 40, and 50, respectively.
After analyzing and discussing the above simulation results, it proves the effectiveness of the designed DEAEBOA. The new DEAEBOA proposed in this paper has obvious advantages in solving the task allocation problem of SMWSNs, has the advantages of fast running speed and strong optimization ability, and provides reliable help for agricultural production.
In order to availably improve the efficiency of task allocation, for this purpose, a novel differential evolutionary adaptive elite butterfly optimization algorithm (DEAEBOA) for reliable task allocation of SMWSNs is proposed. Furthermore, a new adaptive strategy and a fresh elite operator are designed, which can enhance the global search capability and avoid falling into local optimization. Not only that, DEAEBOA's network benefit growth rate can reach up to 16% compared to the other four algorithms. Moreover, the network benefit of DEAEBOA has obvious advantages over BOA, PSO, GA, and BWO, and the network benefit rate is increased to 11.86%, 5.46%, 8.98%, and 12.18% respectively. It can be seen that the proposed DEAEBOA can effectively improve the network benefit and performance of SMWSNs; at the same time, this provides a research basis for the further development of intelligent agriculture.
Although the DEAEBOA proposed in this paper has achieved good performance through simulation experiments, proving its superiority, it still needs some improvements due to various constraints, such as the constraints of research capabilities and environmental conditions. In this paper, the sensor nodes of SMWSNs are statically and randomly distributed in the farmland monitoring area. However, with the development of precision agriculture, some application scenarios require the dynamic distribution of sensor nodes, that is, continuous and regular movement to collect and monitor farmland data. In the future, the influence of environmental factors such as temperature, noise, and obstacles will be considered, the dynamic distribution of sensor nodes will be studied, and the network will be placed in three-dimensional space so that the algorithm will be more practical and adaptable.
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
This paper was funded by the National Natural Science Foundation of China, grant number 61962053, the Shihezi University High-level Talent Research Start-up Fund Project, grant number RCZK2018C39, the project of Youth and middleaged Scientific and Technological Innovation Leading Talents Program of the Corps, grant number 2018CB006, the Corps innovative talents plan, grant number 2020CB001, the China Postdoctoral Science Foundation, grant number 220531, Funding Project for High Level Talents Research in Shihezi University, grant number RCZK2018C38, Project of Shihezi University, grant number ZZZC201915B, National Key R & D Program of China, grant number 2022ZD0115803, Science and Technology Innovation 2030–"New Generation of Artificial Intelligence" Major Project, grant number 2022ZD0115803.
The authors declare no competing interests.
[1] |
Coupling conditions for gas networks governed by the isothermal Euler equations. Netw. Heterog. Media (2006) 1: 295-314 (electronic). ![]() |
[2] |
Gas flow in pipeline networks. Netw. Heterog. Media (2006) 1: 41-56. ![]() |
[3] | A. Bressan, Hyperbolic Systems of Conservation Laws, The One-Dimensional Cauchy Problem, Oxford Lecture Series in Mathematics and its Applications, 20, Oxford University Press, Oxford, 2000. |
[4] |
Flow on networks: recent results and perspectives. European Mathematical Society-Surveys in Mathematical Sciences (2014) 1: 47-111. ![]() |
[5] |
The Cauchy problem for the Euler equations for compressible fluids. Handbook of Mathematical Fluid Dynamics (2002) 1: 421-543. ![]() |
[6] |
Traffic flow on a road network. SIAM J. Math. Anal. (2005) 36: 1862-1886 (electronic). ![]() |
[7] |
Optimal control in networks of pipes and canals. SIAM J. Control Optim. (2009) 48: 2032-2050. ![]() |
[8] |
A well posed Riemann problem for the p-system at a junction. Netw. Heterog. Media (2006) 1: 495-511. ![]() |
[9] |
On the Cauchy problem for the p-system at a junction. SIAM J. Math. Anal. (2008) 39: 1456-1471. ![]() |
[10] |
On 2×2 conservation laws at a junction. SIAM J. Math. Anal. (2008) 40: 605-622. ![]() |
[11] | Thèorie du mouvement non-permanent des eaux, avec application aux crues des rivière at à l'introduction des marèes dans leur lit.. C.R. Acad. Sci. Paris (1871) 73: 147-154. |
[12] | M. Garavello and B. Piccoli, Traffic Flow on Networks, vol. 1 of AIMS Series on Applied Mathematics, American Institute of Mathematical Sciences (AIMS), Springfield, MO, 2006, Conservation laws models. |
[13] |
E. Godlewski and P. -A. Raviart,
Numerical Approximation of Hyperbolic Systems of Conservation Laws, Applied Mathematical Sciences, 118, Springer-Verlag, New York, 1996. doi: 10.1007/978-1-4612-0713-9
![]() |
[14] |
Coupling conditions for a class of second-order models for traffic flow. SIAM J. Math. Anal. (2006) 38: 595-616. ![]() |
[15] |
Assessment of coupling conditions in water way intersections. Internat. J. Numer. Methods Fluids (2013) 71: 1438-1460. ![]() |
[16] |
A mathematical model of traffic flow on a network of unidirectional roads. SIAM J. Math. Anal. (1995) 26: 999-1017. ![]() |
[17] |
Riemann problems with a kink. SIAM J. Math. Anal. (1999) 30: 497-515 (electronic). ![]() |
[18] | S. Joana, M. Joris and T. Evangelos, Technical and Economical Characteristics of Co2 Transmission Pipeline Infrastructure, Technical report, JRC Scientic and Technical Reports, European Commission. |
[19] | Les modeles macroscopiques du traffic. Annales des Ponts. (1993) 67: 24-45. |
[20] |
Experimental assessment of scale effects affecting two-phase flow properties in hydraulic jumps. Experiments in Fluids (2008) 45: 513-521. ![]() |
[21] |
Simulation of transient flow in gas networks. Int. Journal for Numerical Methods in Fluid Dynamics (1984) 4: 13-23. ![]() |
[22] | B. Sultanian, Fluid Mechanics: An Intermediate Approach, CRC Press, 2015. |
[23] |
R. Ugarelli and V. D. Federico, Transition from supercritical to subcritical regime in free surface flow of yield stress fluids Geophys. Res. Lett. , 34 (2007), L21402. doi: 10.1029/2007GL031487
![]() |
1. | Jeffrey Lynham, Walid A. Houry, 2018, Chapter 4, 978-3-030-00736-2, 37, 10.1007/978-3-030-00737-9_4 | |
2. | Bhopal C. Mohapatra, Sameer Mirza, Aditya Bele, Channabasavaiah B. Gurumurthy, Mohsin Raza, Irfana Saleem, Matthew D. Storck, Aniruddha Sarkar, Sai Sundeep Kollala, Surendra K. Shukla, Siddesh Southekal, Kay-Uwe Wagner, Fang Qiu, Subodh M. Lele, Mansour A. Alsaleem, Emad A. Rakha, Chittibabu Guda, Pankaj K. Singh, Robert D. Cardiff, Hamid Band, Vimla Band, Ecdysoneless Overexpression Drives Mammary Tumorigenesis through Upregulation of C-MYC and Glucose Metabolism, 2022, 20, 1541-7786, 1391, 10.1158/1541-7786.MCR-22-0122 | |
3. | Sameer Mirza, Achyuth Kalluchi, Mohsin Raza, Irfana Saleem, Bhopal Mohapatra, Dhananjaya Pal, Michel M. Ouellette, Fang Qiu, Lulu Yu, Alexei Lobanov, Zhi-Ming Zheng, Ying Zhang, Mansour A. Alsaleem, Emad A. Rakha, Hamid Band, M. Jordan Rowley, Vimla Band, Ecdysoneless Protein Regulates Viral and Cellular mRNA Splicing to Promote Cervical Oncogenesis, 2022, 20, 1541-7786, 305, 10.1158/1541-7786.MCR-21-0567 | |
4. | Irfana Saleem, Sameer Mirza, Aniruddha Sarkar, Mohsin Raza, Bhopal Mohapatra, Insha Mushtaq, Jun Hyun Kim, Nitish K. Mishra, Mansour A. Alsaleem, Emad A. Rakha, Fang Qiu, Chittibabu Guda, Hamid Band, Vimla Band, The Mammalian Ecdysoneless Protein Interacts with RNA Helicase DDX39A To Regulate Nuclear mRNA Export, 2021, 41, 1098-5549, 1, 10.1128/MCB.00103-21 | |
5. | Jeffrey Lynham, Walid A. Houry, The Role of Hsp90-R2TP in Macromolecular Complex Assembly and Stabilization, 2022, 12, 2218-273X, 1045, 10.3390/biom12081045 |
Algorithm | Rate of convergence | Local optimization | Computational complexity |
DEAEBOA | fast | \ | low |
IACGA | fast | √ | low |
DPSO | fast | √ | low |
EIACO | fast | √ | high |
SNO | fast | \ | high |
MOA | fast | \ | high |
BOA | fast | √ | high |
Notation | Statement | |
m | Number of soil moisture information collection tasks | |
n | Number of sensor nodes | |
B | Soil moisture information collection advantage matrix | |
U | Economic benefit matrix | |
D | Solution in integer form | |
Z | Solution in binary form | |
Zb | Soil moisture information collection advantage matrix specific solutions | |
F | Soil moisture information collection benefit matrix | |
benefit | Total soil moisture information collection benefit |
DEAEBOA | Population size | Adaptive conversion probability | Mutation probability | Crossover probability | Fragrance | Sensory form |
100 | 0.7–0.9 | 0.95 | 0.95 | 0.1 | 0.01 | |
BOA | Population size | Conversion probability | Fragrance | Sensory form | ||
100 | 0.95 | 0.1 | 0.01 | |||
PSO | Population size | Inertia weight | Self-learning factor | Group-learning factor | ||
100 | 0.5 | 0.6 | 0.6 | |||
GA | Population size | Mutation probability | ||||
100 | 0.1 | |||||
BWO | Population size | Equilibrium factor | ||||
100 | 0.5 |
Number of nodes and tasks | DEAEBOA | PSO | GA | BOA | BWO |
15 tasks and 8 nodes | 6.4452 | 6.2580 | 6.0569 | 5.8687 | 5.8879 |
25 tasks and 8 nodes | 11.1075 | 10.6126 | 10.0674 | 9.9411 | 9.8579 |
40 tasks and 8 nodes | 16.9841 | 15.9036 | 15.3488 | 14.9900 | 14.9503 |
50 tasks and 8 nodes | 20.5443 | 19.136 | 18.9283 | 18.2450 | 18.1966 |
30 tasks and 16 nodes | 12.6718 | 11.9088 | 11.6171 | 11.2362 | 11.1878 |
60 tasks and 30 nodes | 23.3405 | 22.2220 | 22.0922 | 21.1183 | 21.6495 |
120 tasks and 50 nodes | 45.1898 | 42.3056 | 42.1824 | 40.7369 | 41.4927 |
Number of nodes and tasks | PSO | GA | BOA | BWO |
15 tasks and 8 nodes | 3.00% | 6.41% | 9.82% | 9.47% |
25 tasks and 8 nodes | 4.66% | 10.33% | 11.73% | 12.67% |
30 tasks and 16 nodes | 6.40% | 9.08% | 12.78% | 13.26% |
40 tasks and 8 nodes | 6.80% | 10.65% | 13.30% | 13.60% |
50 tasks and 8 nodes | 7.36% | 8.53% | 12.60% | 13.00% |
60 tasks and 30 nodes | 5.03% | 5.65% | 10.52% | 7.81% |
120 tasks and 50 nodes | 6.81% | 7.13% | 10.93% | 9.00% |
Algorithm | Rate of convergence | Local optimization | Computational complexity |
DEAEBOA | fast | \ | low |
IACGA | fast | √ | low |
DPSO | fast | √ | low |
EIACO | fast | √ | high |
SNO | fast | \ | high |
MOA | fast | \ | high |
BOA | fast | √ | high |
Notation | Statement | |
m | Number of soil moisture information collection tasks | |
n | Number of sensor nodes | |
B | Soil moisture information collection advantage matrix | |
U | Economic benefit matrix | |
D | Solution in integer form | |
Z | Solution in binary form | |
Zb | Soil moisture information collection advantage matrix specific solutions | |
F | Soil moisture information collection benefit matrix | |
benefit | Total soil moisture information collection benefit |
DEAEBOA | Population size | Adaptive conversion probability | Mutation probability | Crossover probability | Fragrance | Sensory form |
100 | 0.7–0.9 | 0.95 | 0.95 | 0.1 | 0.01 | |
BOA | Population size | Conversion probability | Fragrance | Sensory form | ||
100 | 0.95 | 0.1 | 0.01 | |||
PSO | Population size | Inertia weight | Self-learning factor | Group-learning factor | ||
100 | 0.5 | 0.6 | 0.6 | |||
GA | Population size | Mutation probability | ||||
100 | 0.1 | |||||
BWO | Population size | Equilibrium factor | ||||
100 | 0.5 |
Number of nodes and tasks | DEAEBOA | PSO | GA | BOA | BWO |
15 tasks and 8 nodes | 6.4452 | 6.2580 | 6.0569 | 5.8687 | 5.8879 |
25 tasks and 8 nodes | 11.1075 | 10.6126 | 10.0674 | 9.9411 | 9.8579 |
40 tasks and 8 nodes | 16.9841 | 15.9036 | 15.3488 | 14.9900 | 14.9503 |
50 tasks and 8 nodes | 20.5443 | 19.136 | 18.9283 | 18.2450 | 18.1966 |
30 tasks and 16 nodes | 12.6718 | 11.9088 | 11.6171 | 11.2362 | 11.1878 |
60 tasks and 30 nodes | 23.3405 | 22.2220 | 22.0922 | 21.1183 | 21.6495 |
120 tasks and 50 nodes | 45.1898 | 42.3056 | 42.1824 | 40.7369 | 41.4927 |
Number of nodes and tasks | PSO | GA | BOA | BWO |
15 tasks and 8 nodes | 3.00% | 6.41% | 9.82% | 9.47% |
25 tasks and 8 nodes | 4.66% | 10.33% | 11.73% | 12.67% |
30 tasks and 16 nodes | 6.40% | 9.08% | 12.78% | 13.26% |
40 tasks and 8 nodes | 6.80% | 10.65% | 13.30% | 13.60% |
50 tasks and 8 nodes | 7.36% | 8.53% | 12.60% | 13.00% |
60 tasks and 30 nodes | 5.03% | 5.65% | 10.52% | 7.81% |
120 tasks and 50 nodes | 6.81% | 7.13% | 10.93% | 9.00% |