Prediction and Analysis of Nuclear Explosion Radioactive Pollutant Diffusion Model

: This study presents a model for the dispersion of radioactive smoke clouds from a nuclear weapon explosion. A model based on a modiﬁed Settlement model is chosen to simulate the dispersion of radioactive contaminants from a nuclear explosion in the atmosphere. The arrival time and dose rate of radioactive fallout at various distances in the downwind direction are given for different equivalents of the surface explosion and typical meteorological conditions. Thus, the prediction of the dispersion of radioactive contaminants from a nuclear explosion can be achieved under the conditions of known nuclear explosion equivalence and local meteorological parameters. This provides a theoretical basis for the estimation of the affected environment and the input of rescue forces after the explosion.


Introduction
Nuclear weapons are the most powerful type ever manufactured and possessed by mankind. The world's current stockpile of nuclear weapons is sufficient to destroy the entire planet [1,2]. On 22 January 2021, the United Nations Treaty on the Prohibition of Nuclear Weapons (TPNW) was enforced, making the development, possession, and use of nuclear weapons illegal [3]. Despite relentless global efforts to prevent the proliferation of nuclear weapons [4,5], as nuclear technology continues to advance, the trend is becoming pronounced toward nuclear confrontation and proliferation at the regional level, and the activation of nuclear terrorism is exposing national security to new threats. If a nuclear war were to break out, it would cause global climate change [6] and unrecoverable disturbances to the marine environment [7], thus affecting agriculture, fisheries, and animal husbandry, and ultimately causing a globalized food problem [8,9] with devastating consequences for the world. In a future war, if nuclear weapons are used in a military conflict, no human beings would be spared worldwide.
In addition to the hot fireballs and shockwave damage that results directly from the detonation of a nuclear weapon, it can also be a serious hazard, i.e., the spread of the radioactive smoke cloud from a nuclear explosion. The dispersal of radioactive contaminants from nuclear explosions refers to the transport of radioactive 137 Cs, 131 I, 90 Sr, and other particles [10][11][12] over time and space, contaminating the ground, water, air, and various objects [13]. This process can last for months, years, or even longer, starting at the moment of the explosion [10,12,14]. Before the first nuclear tests were conducted, scientists had already predicted that such tests would cause radioactive contamination [15]. Many nuclear weapons tests throughout history have resulted in severe radioactive contamination of indigenous peoples in the vicinity of the test sites, and the alarming rates of thyroid and other cancers have increased the importance of research into radioactive fallout [16][17][18][19].
The radioactive smoke cloud produced by a nuclear explosion carries most of the remaining radioactive products after nuclear fission or nuclear fusion, and its drifting under the action of the atmosphere causes serious ground-level radioactive contamination.
This phenomenon is crucial to study the movement and dispersion law of the radioactive smoke cloud of a nuclear explosion in the atmospheric environment. The atmospheric transport processes of radionuclides are affected by weather and meteorological conditions, such that general radiological dispersion models need to consider the driving factors of the meteorological field [20,21]. Leelössy et al. [22] provided a comprehensive review of models for radionuclide dispersion with respect to physical mechanisms, numerical methods, and practical applications. The mainstream atmospheric dispersion algorithms for radionuclides are often based on Gaussian plume models [23][24][25]. Liu et al. [26] improved the Gaussian plume model using ground reflection coefficients and corrected heights to predict the pollutant leak points and concentrations using an inverse calculation method.
For radioactive dispersion in nuclear explosions, the Gaussian plume model is not applicable due to the discontinuity of the source phase, and new computational models have been developed [27]. The earliest WSEG models [28] were based on the "smearing" of smoke clouds, which is equivalent to a cloud of smoke dragging across the ground, causing ground subsidence. Because of its simplicity and ease of application, the WSEG model has been widely used in damage assessment studies for several years. Subsequently, its prediction accuracy has been improved [29,30]. The DNAF-1 model, based on WSEG, is suitable for subsidence simulations of small-yield nuclear explosions and has better accuracy than WSEG [31].
The Defense Land Fallout Interpretative Code (DELFIC) model [32,33] is the standard type for subsidence forecasting in the USA. It describes the entire particle formation and subsidence process in terms of physical principles and takes into account the various factors that influence particle transport. However, this model is very complex because of several factors. Therefore, simplified models have been developed: the SEER model [34], the SEER II model [35], the Hotspot model [36], and the SIMFIC model [37]. The characteristics of these models are that the radioactive smoke cloud of a nuclear explosion is a rising bubble and the height of the smoke cloud is approximately proportional to the square root of the rise time, which is crucial to establishing the vertical trajectory of the particles, while its horizontal movement depends entirely on the wind field.
In order to predict the concentration distribution of radioactive contaminants in the atmosphere from a nuclear explosion in real time, we improved the Settlement model based on MATLAB, using time segmentation processing and coordinate conversion methods such that it can calculate the dispersion of the radioactive contaminants from a nuclear explosion and the change in wind direction. Thus, the application scenario is in line with the real situation of the dispersion of radioactive contaminants from a nuclear explosion. The model can predict changes in the radiation area on the ground from the nuclear explosion and meteorological parameters, and alert personnel downwind, directing them to evacuate the danger zone or take appropriate protection and concealment, thereby reducing the damage caused by the radioactive fallout from the nuclear explosion.

Mathematical Principles of Diffusion Models
In order to obtain accurate and rapid simulation calculation results, different mathematical models should be established for different cases of pollutant dispersion in the atmosphere [22]. After a nuclear explosion occurs, a large number of radionuclides are released. The common diffusion models for the diffusion of hazardous, radioactive substances are the Lagrange model, Gaussian model, and Settlement model. The Lagrange model is a set of first-order stochastic modes describing the motion of many particles. It is a widely used and mature model for radionuclide dispersion simulation, especially in the Fukushima nuclear accident in Japan [38][39][40][41]. The model focuses on the simulation of radionuclide dispersion trajectories. The Gaussian model is based on a non-random variation of the plume with wind speed in the x-direction, a Gaussian distribution in the y-direction, and vertical z-direction at a steady state [42][43][44][45]. The Settlement model is theoretically simple and fast and can provide rapid simulation results in nuclear emergency situations, but it is based on the simplification of the ideal situation, such that the accuracy of the calculation results is slightly inferior compared to other models [44].

Stable Cloud Geometry Data
After detonation, the body of the nuclear bomb releases enormous energy due to atomic fission or fusion reactions, the temperature in the nuclear reaction area rises to tens of thousands of degrees Celsius, and the pressure rises to tens of billions of atmospheres. Under such high temperatures and pressure, the projectile instantly vaporizes into a plasma gas cloud [46], heating and pressurizing the surrounding cold air. In addition, the high temperature and pressure of the gas mass with rapid outward expansion will form a huge fireball that exists for a short period. When the nuclear explosion produced by the fireball is extinguished, it rises and expands rapidly, forms a "mushroom" smoke cloud (Figure 1), and an intense vortex motion is generated inside the smoke cloud [47][48][49]. The strong suction and coiling effect trap the ground dust and other materials in the smoke cloud [50,51]. When the smoke cloud reaches its maximum height, the temperature and pressure inside the cloud are similar to those of the surrounding atmosphere, the vortex motion has disappeared, and the upward motion of the cloud has stopped, thus entering the turbulent separation phase [52], which is called "Stable Cloud".
Pollutants 2022, 2, FOR PEER REVIEW 3 variation of the plume with wind speed in the x-direction, a Gaussian distribution in the y-direction, and vertical z-direction at a steady state [42][43][44][45]. The Settlement model is theoretically simple and fast and can provide rapid simulation results in nuclear emergency situations, but it is based on the simplification of the ideal situation, such that the accuracy of the calculation results is slightly inferior compared to other models [44].

Stable Cloud Geometry Data
After detonation, the body of the nuclear bomb releases enormous energy due to atomic fission or fusion reactions, the temperature in the nuclear reaction area rises to tens of thousands of degrees Celsius, and the pressure rises to tens of billions of atmospheres. Under such high temperatures and pressure, the projectile instantly vaporizes into a plasma gas cloud [46], heating and pressurizing the surrounding cold air. In addition, the high temperature and pressure of the gas mass with rapid outward expansion will form a huge fireball that exists for a short period. When the nuclear explosion produced by the fireball is extinguished, it rises and expands rapidly, forms a "mushroom" smoke cloud (Figure 1), and an intense vortex motion is generated inside the smoke cloud [47][48][49]. The strong suction and coiling effect trap the ground dust and other materials in the smoke cloud [50,51]. When the smoke cloud reaches its maximum height, the temperature and pressure inside the cloud are similar to those of the surrounding atmosphere, the vortex motion has disappeared, and the upward motion of the cloud has stopped, thus entering the turbulent separation phase [52], which is called "Stable Cloud". The height of the cloud cap and the height of the top in the geometry of the stable radioactive smoke cloud can be expressed as follows: Among them, represents the nuclear explosion energy equivalent (kt), the unit of and is meters, and the parameters are as follows: The radius fitting formula of the stable cloud cap is: The height of the cloud cap H B and the height of the top H τ in the geometry of the stable radioactive smoke cloud can be expressed as follows: Among them, W represents the nuclear explosion energy equivalent (kt), the unit of H B and H τ is meters, and the parameters are as follows: The radius fitting formula of the stable cloud cap is: 2 (2) Pollutants 2023, 3 46 where W denotes the nuclear explosion energy equivalent (kt).

Scale Distribution and Activity Distribution of Radioactive Particles
The size distribution of radioactive particles produced during a nuclear explosion is very large, ranging from sub-micron to millimeter scale, and their settling velocity varies considerably, necessitating scale gradation calculations. Typically, the particle scale distribution is log-normal. As different radioactive particles are produced at various time points, the activity varies at different particle scales, according to Bridgman's method, i.e., the particle activity scale distribution is a weighted average of the activity surface distribution and the bulk distribution. The log-normal distribution function has many advantages and is chosen as the scale distribution function in the model.
where: D is the particle diameter; N(D)dD indicates the number of particles scaled between D and D + dD as a percentage of the total; D is the geometric mean particle diameter; S is the geometric standard deviation. The different sizes of radioactive particles in nuclear smoke have varied activity levels, which could be attributed to the different processes and ways in which radioactive particles are formed in the smoke cloud. In the early stages of the nuclear explosion, the high temperatures generated by the bomb cause the bomb body and its surrounding material to vaporize completely and then condense into radioactive particles. This particulate radioactivity is uniformly distributed in a bulk distribution, which has a low particle scale, in addition to a large amount of environmental material, such as soil, which is trapped in it as the nuclear smoke cloud rises [46]. The melted and vaporized radioactive material condenses as it cools and adheres to the surface of these environmental materials, resulting in a surface distribution of these particles, which has a large particle scale.
Due to the different conditions of the explosion, the impact of particle activity distribution is related to factors such as the equivalent burst height and burst area soil composition. Typically, the higher the burst formation of particles, the smaller the geometric mean diameter, the greater the possibility of a radioactive body distribution, and the greater the possibility of forming a radioactive surface distribution. Therefore, due to the air explosion, most of the radioactive material is in the air for a prolonged period, while in the ground explosion, most of the radioactive material quickly settles to the ground.
In this study, Bridgman's methods are used [30], i.e., the activity-scaled distribution of a particle is a weighted average of the activity surface and body distributions, such that the number-scaled distribution of a particle is log-normal. Thus, it can be seen that the surface and body distributions of the particle activity are log-normal, except that the corresponding geometric mean diameters become lnD + 2ln 2 S and lnD + 3ln 2 S , respectively, and the geometric standard deviation is S. Next, we have: where A(D)dD represents the activity in the interval from diameter D to D + dD as a percentage of the total activity of all particles, and As(D) and Av(D) are the surface and bulk distributions of particle activity, respectively: C 1 = 0.32, C 2 = 0.68. In order to simplify the calculation, the particle scale is divided into several intervals, and the percentage of the total activity in any particle size interval can be obtained by integrating the particle size: activity expression for the particle size in the interval [D i , D i+1 ], and the share of the total activity of the radioactive particles F(i) satisfies: Pollutants 2023, 3 47 In this study, equal mass grading was used, i.e., the mass of particles in each subclass is identical. The percentage of total radioactivity and gravitational settling velocity calculated for the median particle radius of the 20 particle size intervals of the ground burst are listed in Table 1. Table 1. Percentage of total radioactivity and gravitational settling velocity in the radius of the median particle of the ground burst. The horizontal transport and dispersion of radioactive particles is based on the Monte Carlo model of smoke clouds; its dispersion is based on Gifford's theory of relative tropospheric dispersion, which is characterized by the stochastic nature of atmospheric motion and the transient nature of explosive smoke cloud generation. The motion of the smoke center could be computed as described previously in the following equation:

Median
Among these, (x n+1 , y n+1 ) is the central position of the puff at time t + ∆t, (x n , y n ) is the position at time and t, and (u n+1 , v n+1 ) is the wind speed at time t + ∆t, calculated by the following formula: Among them, (u n , v n ) is the wind speed at time t, R(∆t) is the correlation coefficient, and (σ u , σ v ) is the turbulent velocity variance. Currently, there is no regular model result of the turbulent velocity variance in the troposphere, and the assumption that the degree of turbulence is not too large can be used, i.e., σ u,v ∼ u x,y /10, where u x,y is the average wind speed of the layer where the smoke cloud is located. The correlation coefficient is in exponential form, i.e., R(∆t) = e −∆t/t L , where t L is the Lagrangian time scale of tropospheric turbulence.
Typically, atmospheric nuclear explosions and radioactive fallout occur in the troposphere, where the tropospheric atmosphere inevitably contributes to nuclear smoke clouds with the passage of time. However, the laws of tropospheric atmospheric dispersion have not been well theorized, especially for contrail dispersion. Gifford et al. [53] proposed a formula for the variance of the horizontal diffusion of puffs under uniform turbulent conditions via random force theory: where σ 0 is the variance of the initial time of the smoke cloud, which is assumed to be R/4.3, K is the vortex diffusion coefficient, and t L is the Lagrangian time scale of large-scale atmospheric motion. c is related to the characteristics of the original source, obtained by comparing the results of the similarity theory using the following formula by Gifford et al. [53]: Since the diffusion coefficient K and the Lagrange time scale are difficult to determine, statistical methods were applied, and K could be roughly estimated by the following formula:

Sedimentation Rate Function
The settling rate fitting nuclear explosion radioactive smoke cloud dispersion calculation model is based on the settling rate function, considering lateral, vertical, and turbulent dispersion effects, calculating the accumulation of radioactive settling in time, and then fitting the ground radiation dose by equivalent dose rate. The sedimentation rate function selected for this model is as follows: where, g(t) is the unit activity sedimentation rate function (s −1 ), t is the time (s), and α is a dimensionless parameter: T is the time scale (s) that can pass through the maximum time point of g(t): t max = 30W 0.41556 ; 10 −3 ≤ W ≤ 1kt t max = 30W 0.65407 ; 1 < W ≤ 10 3 t max = 893.616W 0.16273 ; 10 3 < W ≤ 10 4 t max = 2497.18W 0.05115 ; 10 4 < W ≤ 10 5

Far-Field Correction
The above settling rate function fits the simulated data in the near field or early stages of settling of the radioactive smoke cloud; however, there is a large gap in the far field, such that a far-field correction has to be added.

Vertical Diffusion
The settling velocity of the particles of nominal particle size at different heights can be expressed as: (14) where f 0 = 1.6538 (m/s) and ς = 2.9 × 10 −5 m −1 . Then, the time for the particle to reach the ground is: To make it easier to give the longitudinal diffusion parameters of a nuclear explosion, first define a diffusion parameter related to the explosion yield σ W : Then, the longitudinal diffusion parameters at different times and downwind positions are:

Turbulent Diffusion Correction
Considering the correction for longitudinal diffusion by atmospheric turbulence, the result is:

Lateral Diffusion Parameters
The lateral diffusion is mainly due to the influence of wind shear and its diffusion parameters can be expressed as follows: where S y denotes the approximate vertical wind shear (m/s). Thus, the equation is as follows:

Sedimentation Rate Diffusion Factor
Based on the parameters provided above, the settling rate diffusion factor at different arrival times and locations can be expressed as follows: Pollutants 2023, 3

50
V is the equivalent wind speed (m/s). For the upwind direction, a correction is required, multiplied by the wind correction factor F up above.
If W ≤ 10, b = 0.0176 and c 1 = 570; else if W > 10, b = 0.08045w −0.66 and c 1 = −8179.82 + 3800lnW. Settling is a cumulative effect in time and this model does not consider resuspension; hence, the settling diffusion factor xq m −2 is expressed as:

Programming
The program was written using MATLAB language to simulate the dispersion of radioactive contaminants from a nuclear explosion. Figure 2 shows the flowchart of the program; the input file contains the initial module of a nuclear explosion and the meteorological monitoring module. The initial module of the nuclear explosion includes the explosion equivalent, the coordinates of the burst center, and the exit velocity. The meteorological monitoring module includes the meteorological parameters, such as wind speed, wind direction, and temperature, and transmits them to the simulation calculation module. The simulation calculation module processes the source data, meteorological data, the calculation area size, and other parameters, and substitutes them into the mathematical model for the concentration distribution calculation. The final results of the theoretical numerical prediction of the dose concentration distribution of the radioactive smoke cloud deposition in the atmosphere of the nuclear explosion are obtained.

Sedimentation Rate Diffusion Factor Validation
The nuclear explosion radioactive dispersion prediction model DNAF-1 [31] is a com-

Sedimentation Rate Diffusion Factor Validation
The nuclear explosion radioactive dispersion prediction model DNAF-1 [31] is a computational tool for rapid assessment of the effects of radioactive fallout from a nuclear explosion on the surface based on experimental and simulation data patterns summarized in the USA in the early 1980s [32,33,37]. In order to achieve a rapid and effective evaluation of the nuclear explosion radioactive contaminant dispersion model, compare its accuracy, and verify its practical application. The major computational parameters are referred to for comparison. In the rapid deposition model of the nuclear explosion radioactive smoke cloud, the model for calculating the settling rate dispersion factor is the core of the calculation that determines the whole model algorithm; hence, a direct comparison of its consistency can illustrate the accuracy of the model. Figure 3 shows the distribution of the settling rate diffusion factor with time for the different models at 1 KT explosion equivalent. Compared to the case without adding the far-field approximation correction, each model is calibrated to a specific degree; the comparison reveals that the settling rate diffusion factor of our model is similar to that of the DELFIC model and of the DNAF-1 model with a certain degree of confidence.

Real Historical Data Verification
In this study, the historical nuclear explosion experimental data shown in Table 2 verified the accuracy of the model. This model was then compared to the DELFIC, DNAF-1, and WSEG-10 models and then with the real observed data. The contour points generated by the model at the same dose were plotted by the visualization software R, as shown in Figure 4. Their heat maps were supplemented with blank values on the grid points using linear interpolation. The data were selected from the three historical nuclear explosions under code names Sugar, Koon, and Zuni, with explosive yields of 1.2, 150, and 3380 kilotons, respectively, to verify the prediction accuracy of the model for large, medium, and small yields. The inter-model validation and comparison with real recorded observations can be obtained in medium-and small-yield nuclear explosion tests. Our model showed good prediction accuracy, indicating considerable credibility and practical value. However, in large yields, similar to the other three models, our model showed large deviations from

Real Historical Data Verification
In this study, the historical nuclear explosion experimental data shown in Table 2 verified the accuracy of the model. This model was then compared to the DELFIC, DNAF-1, and WSEG-10 models and then with the real observed data. The contour points generated by the model at the same dose were plotted by the visualization software R, as shown in Figure 4. Their heat maps were supplemented with blank values on the grid points using linear interpolation. The data were selected from the three historical nuclear explosions under code names Sugar, Koon, and Zuni, with explosive yields of 1.2, 150, and 3380 kilotons, respectively, to verify the prediction accuracy of the model for large, medium, and small yields.

General Validation
Suppose that a nuclear explosion equivalent to 100 kilotons occurs at a certain location and the meteorological conditions of the explosion are an average wind speed of 6 m/s, an average wind direction of 90°, and an average wind shear of 0. The nuclear explosion and meteorological parameters are input into our model, and the predicted contamination dose of radioactive fallout from 1 to 6 h after the explosion is presented by the visualization method described in Figure 5. The radioactive smoke cloud deposition from a nuclear explosion has the following characteristics: • The maximum radioactive dose and the area of radioactive contamination increase with time in 1-6 h; • A hot line of radioactive dose extends continuously along the downwind direction; • The general trend of dose distribution along the direction of the hot line is a gradual decrease with increasing distance; • The predicted map of radioactive deposition dose varies for each nuclear explosion; • Radioactive deposition is strongly influenced by the equivalent dose, meteorological conditions, topographic environment, and blast mode. The inter-model validation and comparison with real recorded observations can be obtained in medium-and small-yield nuclear explosion tests. Our model showed good prediction accuracy, indicating considerable credibility and practical value. However, in large yields, similar to the other three models, our model showed large deviations from the real recorded observations, which might be due to the fact that large-yield nuclear explosions have a high altitude of the radioactive smoke cloud, which spreads over a wide area and is severely affected by the complex meteorological conditions; therefore, all models are less effective in predicting the cloud. Owing to the fact that the Zuni experiment was a nuclear test conducted on Bikini Island, the actual dose detection records were affected by the constraints of the marine environment and the number of ships, mobility, and observation errors, and their observations were not fully accurate isochronous distribution values of the actual dose dispersion; hence, the credibility of the observation data was limited and could only be used as a reference.

General Validation
Suppose that a nuclear explosion equivalent to 100 kilotons occurs at a certain location and the meteorological conditions of the explosion are an average wind speed of 6 m/s, an average wind direction of 90 • , and an average wind shear of 0. The nuclear explosion and meteorological parameters are input into our model, and the predicted contamination dose of radioactive fallout from 1 to 6 h after the explosion is presented by the visualization method described in Figure 5. The radioactive smoke cloud deposition from a nuclear explosion has the following characteristics:

•
The maximum radioactive dose and the area of radioactive contamination increase with time in 1-6 h; • A hot line of radioactive dose extends continuously along the downwind direction; • The general trend of dose distribution along the direction of the hot line is a gradual decrease with increasing distance; • The predicted map of radioactive deposition dose varies for each nuclear explosion; • Radioactive deposition is strongly influenced by the equivalent dose, meteorological conditions, topographic environment, and blast mode.
Pollutants 2022, 2, FOR PEER REVIEW 12 Figure 5. Predicted dose dispersion of radioactive hazardous substances from 1 to 6 h after a nuclear explosion.
Therefore, through the above historical data model comparison verification and general verification, it is reasonable to judge that our model as the prediction of nuclear explosion radioactive smoke cloud deposition in terms of time distribution or accuracy has certain reference value and can be used as one of the theoretical bases for nuclear explosion hazard prediction judgment analysis.

Conclusions
Nuclear explosion radioactive fallout is the nuclear reaction after the fission or fusion production in the condensation of particles formed in the explosion of the fireball of buoyancy to rise to high altitude [46,49]. Then in the atmosphere under the action of diffusion, settlement to the ground, the formation of a vast area of radioactive contamination. The effect of nuclear explosion radioactive fallout in the nuclear weapons destruction effect is the largest and the duration is the longest [11,12]. Numerical computational models are the most effective tools for simulating radioactive fallout from nuclear explosions, but model accuracy is still an important factor limiting the accuracy of predictions [22]. However, the prediction accuracy of nuclear explosion radioactive smoke cloud deposition is Therefore, through the above historical data model comparison verification and general verification, it is reasonable to judge that our model as the prediction of nuclear explosion radioactive smoke cloud deposition in terms of time distribution or accuracy has certain reference value and can be used as one of the theoretical bases for nuclear explosion hazard prediction judgment analysis.

Conclusions
Nuclear explosion radioactive fallout is the nuclear reaction after the fission or fusion production in the condensation of particles formed in the explosion of the fireball of buoyancy to rise to high altitude [46,49]. Then in the atmosphere under the action of diffusion, settlement to the ground, the formation of a vast area of radioactive contamination. The effect of nuclear explosion radioactive fallout in the nuclear weapons destruction effect is the largest and the duration is the longest [11,12]. Numerical computational models are the most effective tools for simulating radioactive fallout from nuclear explosions, but model accuracy is still an important factor limiting the accuracy of predictions [22]. However, the prediction accuracy of nuclear explosion radioactive smoke cloud deposition is not only related to the diffusion model, but also closely related to the accuracy of the input nuclear explosion equivalent, coupling with the atmospheric diffusion model, etc. In this paper, a hybrid model calculation program for nuclear explosion radioactive smoke cloud deposition with variable wind direction was developed using MATLAB based on the deposition model, which can simulate and calculate the dose distribution of ground-level radioactive deposition with time under the scenario of changing wind conditions. Based on this, the effects of wind speed and atmospheric stability on the dispersion of radioactive contaminants are simulated and analyzed. And through historical data validation, model comparison, and general validation it can be concluded that our model can predict the occurrence range and intensity of radioactive deposition based on the explosion equivalent and the meteorological data at the moment of explosion, and it can be applied to develop the corresponding preventive protection and rescue measures with a certain degree of confidence.
The limitation of this paper is that the subsidence model itself is a simplified model for predicting the dispersion of radioactive smoke clouds, which is fast in calculation but sacrifices some prediction accuracy. In future work, we can try to combine the Gaussian smoke cluster model, the Lagrangian model, or the Eulerian hybrid model for comprehensive testing and comparison, which will hopefully further improve the prediction accuracy of radioactive smoke cloud dispersion of nuclear explosion.