In order to check the validity of the mathematical model for analyzing the flow field in the air-agitated seed precipitation tank,a scaled down experimental apparatus was designed and the colored tracer and KCl tracer...In order to check the validity of the mathematical model for analyzing the flow field in the air-agitated seed precipitation tank,a scaled down experimental apparatus was designed and the colored tracer and KCl tracer were added in the apparatus to follow the real flow line.Virtue tracers were considered in the mathematical model and the algorithm of tracers was built.The comparison of the results between the experiment and numerical calculation shows that the time of the tracer flows out of stirring tube are 40 s in the experiment and 42 s in numerical calculated result.The transient diffusion process and the solution residence time of the numerical calculation are in good agreement with the experimental results,which indicates that the mathematical model is reliable and can be used to predict the flow field of the air-agitated seed precipitation tank.展开更多
In this paper, cold simulation experiments and numerical calculations are conducted to predict 3 D flow field aerodynamics for an oil furnace after being retrofitted due to its fuel variation. K ε model and SIMPLE ...In this paper, cold simulation experiments and numerical calculations are conducted to predict 3 D flow field aerodynamics for an oil furnace after being retrofitted due to its fuel variation. K ε model and SIMPLE program under body fit coordination (BFC) system, in which TTM non orthogonal method is used to control the irregular geometric boundary, are adopted to solve the control equations. Model tests are conducted to check the calculation results, showing that they are in agreement with each other. Three different alternatives with different side window locations are also calculated to optimize the designs. The field retrofitting results show that the combination of cold tests with numerical calculations has prosperous application in retrofitting or renewing medium and small boilers.展开更多
In this paper the computer code of AUTOUGH2 is used to carry out numerical simulation study on the Wuqing geothermal reservoir for evaluating reservoir performance upon different exploitation schemes.The simulating ca...In this paper the computer code of AUTOUGH2 is used to carry out numerical simulation study on the Wuqing geothermal reservoir for evaluating reservoir performance upon different exploitation schemes.The simulating calculation and analysis of two main cases,production without reinjection and production with reinjection,were carried out to illuminate and compare the temperature and pressure profiles with distance and time.The water level or pressure will decrease with the time of production,and the reinjection will cause the cooling effect of reinjection on the production zone.展开更多
This article describes numerical simulation of gas pipeline network operation using high-accuracy computational fluid dynamics (CFD) simulators of the modes of gas mixture transmission through long, multi-line pipelin...This article describes numerical simulation of gas pipeline network operation using high-accuracy computational fluid dynamics (CFD) simulators of the modes of gas mixture transmission through long, multi-line pipeline systems (CFD-simulator). The approach used in CFD-simulators for modeling gas mixture transmission through long, branched, multi-section pipelines is based on tailoring the full system of fluid dynamics equations to conditions of unsteady, non-isothermal processes of the gas mixture flow. Identification, in a CFD-simulator, of safe parameters for gas transmission through compressor stations amounts to finding the interior points of admissible sets described by systems of nonlinear algebraic equalities and inequalities. Such systems of equalities and inequalities comprise a formal statement of technological, design, operational and other constraints to which operation of the network equipment is subject. To illustrate the practicability of the method of numerical simulation of a gas transmission network, we compare computation results and gas flow parameters measured on-site at the gas transmission enter-prise.展开更多
The fatigue life of top tensioned risers under vortex-induced vibrations (VIVs) with consideration of the effect of internal flowing fluid on the riser is analyzed in the time domain.The long-term stress histories of ...The fatigue life of top tensioned risers under vortex-induced vibrations (VIVs) with consideration of the effect of internal flowing fluid on the riser is analyzed in the time domain.The long-term stress histories of the riser under VIVs are calculated and the mean stresses,the number of stress cycles and amplitudes are determined by the rainflow counting method.The Palmgren-Miner rule for cumulative damage theory with a specified S-N curve is used to estimate the fatigue life of the riser.The corresponding numerical programs numerical simulation of vortex-induced vibrations (NSVIV) which can be used to calculate the VIV response and fatigue life of the riser are compiled.Finally the influences of the riser's parameters such as flexural rigidity,top tension and internal flow velocity on the fatigue life of the riser are analyzed in detail and some conclusions are drawn.展开更多
This paper summarizes some of the typical riser vortex-induced vibration (VIV) problems in subsea oil and gas developments, and presents the corresponding computational fluid dynamics (CFD) time domain simula- tio...This paper summarizes some of the typical riser vortex-induced vibration (VIV) problems in subsea oil and gas developments, and presents the corresponding computational fluid dynamics (CFD) time domain simula- tion results to address these problems. First, the CFD time domain simulation approach was applied to analyze the wake field behind a stationary cylinder and a vibrating cylinder. Then a vertical riser VIV response under uniform current was studied. The VIV response time histories revealed some valuable clues that could lead to explanation of the higher harmonics. After that, a vertical riser VIV response under shear current was investigated. A 3 000 ft (1 ft=-0.304 8 m) water depth top tensioned riser was sized, and its VIV responses under uniform and shear current were studied. Then this paper continues to discuss one catenary flexible riser VIV response during normal lay. Last, the time domain simulation approach was applied to a partially submerged flexible jumper, to study the jumper VIV behavior, and dynamic motion envelopes. It was demonstrated that the time domain simulation ap- proach is able to disclose details of the flow field, vortex shedding pattern, and riser dynamic behavior, and han- dle different tvoes of risers under different Woe of currents.展开更多
As a new type of structure which has never been built, submerged floating tunnel was studied mainly by numerical simulations. To further study the seismic response of a submerged floating tunnel, the first model exper...As a new type of structure which has never been built, submerged floating tunnel was studied mainly by numerical simulations. To further study the seismic response of a submerged floating tunnel, the first model experiment of submerged floating tunnel (SFT) under the earthquake was carried out on the unique underwater shaking table in China. The experimental results show that vertical excitation induces larger response than horizontal and different inclination degrees of the tether also cause different seismic responses. Subsequently, based on the fluid-structure interaction theory, the corresponding numerical model is established. And comparing the numerical results with the experimental results, those of shaking table test. Numerical model adopted is effective for it is shown that the numerical results are basically identical with dynamic response of SFT.展开更多
This paper reports an investigation of Computational Fluid Dynamics(CFD)on the influence of injection momentum rate of premixed air and fuel on the flameless Moderate or Intense Low oxygen Dilution(MILD) combustion in...This paper reports an investigation of Computational Fluid Dynamics(CFD)on the influence of injection momentum rate of premixed air and fuel on the flameless Moderate or Intense Low oxygen Dilution(MILD) combustion in a recuperative furnace.Details of the furnace flow velocity,temperature,O2,CO2 and NOx concentrations are provided.Results obtained suggest that the flue gas recirculation plays a vital role in establishing the premixed MILD combustion.It is also revealed that there is a critical momentum rate of the fuel-air mixture below which MILD combustion does not occur.Moreover,the momentum rate appears to have less significant influence on conventional global combustion than on MILD combustion.展开更多
Testing centrifugal fan flow field by physical laboratory is difficult because the testing system is complex and the workload is heavy, and the results observed by naked-eye deviates far from the actual value. To addr...Testing centrifugal fan flow field by physical laboratory is difficult because the testing system is complex and the workload is heavy, and the results observed by naked-eye deviates far from the actual value. To address this problem, the computational fluid dynamics software FLUENT was applied to establish three-dimensional model of the centrifugal fan. The numeral model was verified by comparing simulation data to experimental data. The pressure centrifugal fan and the speed changes in distribution in centrifugal fan was simulated by computational fluid dynamics soft-ware FLUENT. The simulation results show that the gas flow velocity in the impeller increases with impeller radius increase. Static pressure gradually increases when gas from the fan access is imported through fan impeller leaving fans.展开更多
In this paper, two novel numerical computation methods are introduced which have been recently developed at Research Institute for Applied Mechanics ( R/AM ), Kyushu University, for strongly nonlinear wave-body inte...In this paper, two novel numerical computation methods are introduced which have been recently developed at Research Institute for Applied Mechanics ( R/AM ), Kyushu University, for strongly nonlinear wave-body interaction problems, such as ship motions in rough seas and resulting green-water impact on deck. The first method is the CIP-based Cartesian grid method, in which the free surface flow is treated as a multi-phase flow which is solved using a Cartesian grid. The second method is the MPS method, which is a so-called particle method and hence no grid is used. The features and calculation procedures of these numerical methods are described. One validation computation against a newly conducted experiment on a dam break problem, which is also described in this paper, is presented.展开更多
A multiphysics model for a production scale planar solid oxide fuel cell (SOFC) stack is important for the SOFC technology, but usually requires an unpractical amount of computing resource. The major cause for the h...A multiphysics model for a production scale planar solid oxide fuel cell (SOFC) stack is important for the SOFC technology, but usually requires an unpractical amount of computing resource. The major cause for the huge computing resource requirement is identified as the need to solve the cathode O2 transport and the associated electrochemistry. To overcome the technical obstacle, an analytical model for solving the O2 transport and its coupling with the electrochemistry is derived. The analytical model is used to greatly reduce the numerical mesh complexity of a multiphysics model. Numerical test shows that the analytical approximation is highly accurate and stable. A multiphysics numerical modeling tool taking advantage of the analytical solution is then developed through Fluent@. The numerical efficiency and stability of this modeling tool are further demonstrated by simulating a 30- cell stack with a production scale cell size. Detailed information about the stack performance is revealed and briefly discussed. The multiphysics modeling tool can be used to guide the stack design and select the operating parameters.展开更多
Unlike the limit equilibrium method(LEM), with which only the global safety factor of the landslide can be calculated, a local safety factor(LSF) method is proposed to evaluate the stability of different sections of a...Unlike the limit equilibrium method(LEM), with which only the global safety factor of the landslide can be calculated, a local safety factor(LSF) method is proposed to evaluate the stability of different sections of a landslide in this paper. Based on three-dimensional(3D) numerical simulation results, the local safety factor is defined as the ratio of the shear strength of the soil at an element on the slip zone to the shear stress parallel to the sliding direction at that element. The global safety factor of the landslide is defined as the weighted average of all local safety factors based on the area of the slip surface. Some example analyses show that the results computed by the LSF method agree well with those calculated by the General Limit Equilibrium(GLE) method in two-dimensional(2D) models and the distribution of the LSF in the 3D slip zone is consistent with that indicated by the observed deformation pattern of an actual landslide in China.展开更多
The fluid dynamic behavior of feeding gas (TiCI4) in an annular channel affects the combination of 02 and TiCI4 in an oxidation reactor, a key piece of equipment in titanium pigment production. The numerical procedu...The fluid dynamic behavior of feeding gas (TiCI4) in an annular channel affects the combination of 02 and TiCI4 in an oxidation reactor, a key piece of equipment in titanium pigment production. The numerical procedure was validated by a 3-dimensional gas flow in the annular channel. Applying the validated model, the flow character- istics of TiCl4 in the oxidation reactor with a tangential inlet were simulated and characterized. The flow distribu- tion with five rectifying rings of different structure was simulated and analyzed. The results showed that the rectifying ring improved the distribution uniformity of the pressure and outlet velocity. Compared to the original case without a rectifying ring, the non-uniformity of the pressure and outlet velocity could be reduced by up to 91% and 69% respectively. The rectifying ring #5, which can be installed and adjusted easily, is more effective in realizing even distribution. In addition, installation of the rectifying ring effectively reduced the circulating flow in an annular channel as well as the total energy loss.展开更多
The phenomenon of direct-contact condensation,used in steam driven jet injectors,nuclear reactor emergency core cooling systems and direct-contact heat exchangers,was investigated computationally by introducing a ther...The phenomenon of direct-contact condensation,used in steam driven jet injectors,nuclear reactor emergency core cooling systems and direct-contact heat exchangers,was investigated computationally by introducing a thermal equilibrium model for direct-contact condensation of steam in subcooled water.The condensation model presented was a two resistance model which takes care of the heat transfer process on both sides of the interface and uses a variable steam bubble diameter.The injection of supersonic steam jet in subcooled water tank was simulated using the Euler-Euler multiphase flow model of Fluent 6.3 code with the condensation model incorporated. The findings of the computational fluid dynamics(CFD) simulations were compared with the published experimental data and fairly good agreement was observed between the two,thus validating the condensation model.The results of CFD simulations for dimensionless penetration length of steam plume varies from 2.73-7.33,while the condensation heat transfer coefficient varies from 0.75-0.917 MW·(m ^2 ·K)^ -1 for water temperature in the range of 293-343 K.展开更多
The research on spatial epidemic models is a topic of considerable recent interest. In another hand, the advances in computer technology have stimulated the development of stochastic models. Metapopulation models are ...The research on spatial epidemic models is a topic of considerable recent interest. In another hand, the advances in computer technology have stimulated the development of stochastic models. Metapopulation models are spatial designs that involve movements of individuals between distinct subpopulations. The purpose of the present work has been to develop stochastic models in order to study the transmission dynamics and control of infectious diseases in metapopulations. The authors studied Susceptible-Infected-Susceptible (SIS) and Susceptible-lnfected-Recovered (SIR) epidemic schemes, using the Gillespie algorithm, Computational numerical simulations were carried in order to explore the models. The results obtained show how the dynamics of transmission and the application of control measures within each subpopulation may affect all subpopulations of the system. They also show how the distribution of control measures among subpopulations affects the efficacy of these strategies. The dynamics of the stochastic models developed in the current study follow the trends observed in the classic deterministic designs. Also, the present models exhibit fluctuating behavior. This work highlights the importance of the spatial distribution of the population in spread and control of infectious diseases. In addition, it shows how chance could play an important role in these scenarios.展开更多
基金Project(07JJ4016) supported by the Natural Science Foundation of Hunan Procvince,China
文摘In order to check the validity of the mathematical model for analyzing the flow field in the air-agitated seed precipitation tank,a scaled down experimental apparatus was designed and the colored tracer and KCl tracer were added in the apparatus to follow the real flow line.Virtue tracers were considered in the mathematical model and the algorithm of tracers was built.The comparison of the results between the experiment and numerical calculation shows that the time of the tracer flows out of stirring tube are 40 s in the experiment and 42 s in numerical calculated result.The transient diffusion process and the solution residence time of the numerical calculation are in good agreement with the experimental results,which indicates that the mathematical model is reliable and can be used to predict the flow field of the air-agitated seed precipitation tank.
文摘In this paper, cold simulation experiments and numerical calculations are conducted to predict 3 D flow field aerodynamics for an oil furnace after being retrofitted due to its fuel variation. K ε model and SIMPLE program under body fit coordination (BFC) system, in which TTM non orthogonal method is used to control the irregular geometric boundary, are adopted to solve the control equations. Model tests are conducted to check the calculation results, showing that they are in agreement with each other. Three different alternatives with different side window locations are also calculated to optimize the designs. The field retrofitting results show that the combination of cold tests with numerical calculations has prosperous application in retrofitting or renewing medium and small boilers.
文摘In this paper the computer code of AUTOUGH2 is used to carry out numerical simulation study on the Wuqing geothermal reservoir for evaluating reservoir performance upon different exploitation schemes.The simulating calculation and analysis of two main cases,production without reinjection and production with reinjection,were carried out to illuminate and compare the temperature and pressure profiles with distance and time.The water level or pressure will decrease with the time of production,and the reinjection will cause the cooling effect of reinjection on the production zone.
文摘This article describes numerical simulation of gas pipeline network operation using high-accuracy computational fluid dynamics (CFD) simulators of the modes of gas mixture transmission through long, multi-line pipeline systems (CFD-simulator). The approach used in CFD-simulators for modeling gas mixture transmission through long, branched, multi-section pipelines is based on tailoring the full system of fluid dynamics equations to conditions of unsteady, non-isothermal processes of the gas mixture flow. Identification, in a CFD-simulator, of safe parameters for gas transmission through compressor stations amounts to finding the interior points of admissible sets described by systems of nonlinear algebraic equalities and inequalities. Such systems of equalities and inequalities comprise a formal statement of technological, design, operational and other constraints to which operation of the network equipment is subject. To illustrate the practicability of the method of numerical simulation of a gas transmission network, we compare computation results and gas flow parameters measured on-site at the gas transmission enter-prise.
基金supported by the High Technology Research and Development Program of China (863 Program,Grant Nos.2006AA09Z356 and 2007AA09Z313)
文摘The fatigue life of top tensioned risers under vortex-induced vibrations (VIVs) with consideration of the effect of internal flowing fluid on the riser is analyzed in the time domain.The long-term stress histories of the riser under VIVs are calculated and the mean stresses,the number of stress cycles and amplitudes are determined by the rainflow counting method.The Palmgren-Miner rule for cumulative damage theory with a specified S-N curve is used to estimate the fatigue life of the riser.The corresponding numerical programs numerical simulation of vortex-induced vibrations (NSVIV) which can be used to calculate the VIV response and fatigue life of the riser are compiled.Finally the influences of the riser's parameters such as flexural rigidity,top tension and internal flow velocity on the fatigue life of the riser are analyzed in detail and some conclusions are drawn.
文摘This paper summarizes some of the typical riser vortex-induced vibration (VIV) problems in subsea oil and gas developments, and presents the corresponding computational fluid dynamics (CFD) time domain simula- tion results to address these problems. First, the CFD time domain simulation approach was applied to analyze the wake field behind a stationary cylinder and a vibrating cylinder. Then a vertical riser VIV response under uniform current was studied. The VIV response time histories revealed some valuable clues that could lead to explanation of the higher harmonics. After that, a vertical riser VIV response under shear current was investigated. A 3 000 ft (1 ft=-0.304 8 m) water depth top tensioned riser was sized, and its VIV responses under uniform and shear current were studied. Then this paper continues to discuss one catenary flexible riser VIV response during normal lay. Last, the time domain simulation approach was applied to a partially submerged flexible jumper, to study the jumper VIV behavior, and dynamic motion envelopes. It was demonstrated that the time domain simulation ap- proach is able to disclose details of the flow field, vortex shedding pattern, and riser dynamic behavior, and han- dle different tvoes of risers under different Woe of currents.
基金Projects(51108224,51179026) supported by the National Natural Science Foundation of China
文摘As a new type of structure which has never been built, submerged floating tunnel was studied mainly by numerical simulations. To further study the seismic response of a submerged floating tunnel, the first model experiment of submerged floating tunnel (SFT) under the earthquake was carried out on the unique underwater shaking table in China. The experimental results show that vertical excitation induces larger response than horizontal and different inclination degrees of the tether also cause different seismic responses. Subsequently, based on the fluid-structure interaction theory, the corresponding numerical model is established. And comparing the numerical results with the experimental results, those of shaking table test. Numerical model adopted is effective for it is shown that the numerical results are basically identical with dynamic response of SFT.
基金Supported by the National Natural Science Foundation of China(50936001)
文摘This paper reports an investigation of Computational Fluid Dynamics(CFD)on the influence of injection momentum rate of premixed air and fuel on the flameless Moderate or Intense Low oxygen Dilution(MILD) combustion in a recuperative furnace.Details of the furnace flow velocity,temperature,O2,CO2 and NOx concentrations are provided.Results obtained suggest that the flue gas recirculation plays a vital role in establishing the premixed MILD combustion.It is also revealed that there is a critical momentum rate of the fuel-air mixture below which MILD combustion does not occur.Moreover,the momentum rate appears to have less significant influence on conventional global combustion than on MILD combustion.
文摘Testing centrifugal fan flow field by physical laboratory is difficult because the testing system is complex and the workload is heavy, and the results observed by naked-eye deviates far from the actual value. To address this problem, the computational fluid dynamics software FLUENT was applied to establish three-dimensional model of the centrifugal fan. The numeral model was verified by comparing simulation data to experimental data. The pressure centrifugal fan and the speed changes in distribution in centrifugal fan was simulated by computational fluid dynamics soft-ware FLUENT. The simulation results show that the gas flow velocity in the impeller increases with impeller radius increase. Static pressure gradually increases when gas from the fan access is imported through fan impeller leaving fans.
文摘In this paper, two novel numerical computation methods are introduced which have been recently developed at Research Institute for Applied Mechanics ( R/AM ), Kyushu University, for strongly nonlinear wave-body interaction problems, such as ship motions in rough seas and resulting green-water impact on deck. The first method is the CIP-based Cartesian grid method, in which the free surface flow is treated as a multi-phase flow which is solved using a Cartesian grid. The second method is the MPS method, which is a so-called particle method and hence no grid is used. The features and calculation procedures of these numerical methods are described. One validation computation against a newly conducted experiment on a dam break problem, which is also described in this paper, is presented.
基金This work is supported the National Natural Science Foundation of China (No. 11374272 and No. 11574284), the National Basic Research Program of China (No.2012CB215405) and Collaborative Innovation Center of Suzhou Nano Science and Technology are gratefully acknowledged.
文摘A multiphysics model for a production scale planar solid oxide fuel cell (SOFC) stack is important for the SOFC technology, but usually requires an unpractical amount of computing resource. The major cause for the huge computing resource requirement is identified as the need to solve the cathode O2 transport and the associated electrochemistry. To overcome the technical obstacle, an analytical model for solving the O2 transport and its coupling with the electrochemistry is derived. The analytical model is used to greatly reduce the numerical mesh complexity of a multiphysics model. Numerical test shows that the analytical approximation is highly accurate and stable. A multiphysics numerical modeling tool taking advantage of the analytical solution is then developed through Fluent@. The numerical efficiency and stability of this modeling tool are further demonstrated by simulating a 30- cell stack with a production scale cell size. Detailed information about the stack performance is revealed and briefly discussed. The multiphysics modeling tool can be used to guide the stack design and select the operating parameters.
基金financially supported by the National Natural Science Foundation of China(Grant No.51178402,10902112)Department of Transportation Technology Projects(Grant No.2011318740240)the Fundamental Research Funds for the Central Universities(Grant No.2682014CX074)
文摘Unlike the limit equilibrium method(LEM), with which only the global safety factor of the landslide can be calculated, a local safety factor(LSF) method is proposed to evaluate the stability of different sections of a landslide in this paper. Based on three-dimensional(3D) numerical simulation results, the local safety factor is defined as the ratio of the shear strength of the soil at an element on the slip zone to the shear stress parallel to the sliding direction at that element. The global safety factor of the landslide is defined as the weighted average of all local safety factors based on the area of the slip surface. Some example analyses show that the results computed by the LSF method agree well with those calculated by the General Limit Equilibrium(GLE) method in two-dimensional(2D) models and the distribution of the LSF in the 3D slip zone is consistent with that indicated by the observed deformation pattern of an actual landslide in China.
基金Supported by the National Natural Science Foundation of China(21566015,51404123)the Applied Basic Research Projects of Yunnan(2013FC002,2013FC003,2015FB126)the Provincial Personnel Training Funds of Yunnan Province of China(KKSY201352109,KKPT201563013)
文摘The fluid dynamic behavior of feeding gas (TiCI4) in an annular channel affects the combination of 02 and TiCI4 in an oxidation reactor, a key piece of equipment in titanium pigment production. The numerical procedure was validated by a 3-dimensional gas flow in the annular channel. Applying the validated model, the flow character- istics of TiCl4 in the oxidation reactor with a tangential inlet were simulated and characterized. The flow distribu- tion with five rectifying rings of different structure was simulated and analyzed. The results showed that the rectifying ring improved the distribution uniformity of the pressure and outlet velocity. Compared to the original case without a rectifying ring, the non-uniformity of the pressure and outlet velocity could be reduced by up to 91% and 69% respectively. The rectifying ring #5, which can be installed and adjusted easily, is more effective in realizing even distribution. In addition, installation of the rectifying ring effectively reduced the circulating flow in an annular channel as well as the total energy loss.
文摘The phenomenon of direct-contact condensation,used in steam driven jet injectors,nuclear reactor emergency core cooling systems and direct-contact heat exchangers,was investigated computationally by introducing a thermal equilibrium model for direct-contact condensation of steam in subcooled water.The condensation model presented was a two resistance model which takes care of the heat transfer process on both sides of the interface and uses a variable steam bubble diameter.The injection of supersonic steam jet in subcooled water tank was simulated using the Euler-Euler multiphase flow model of Fluent 6.3 code with the condensation model incorporated. The findings of the computational fluid dynamics(CFD) simulations were compared with the published experimental data and fairly good agreement was observed between the two,thus validating the condensation model.The results of CFD simulations for dimensionless penetration length of steam plume varies from 2.73-7.33,while the condensation heat transfer coefficient varies from 0.75-0.917 MW·(m ^2 ·K)^ -1 for water temperature in the range of 293-343 K.
文摘The research on spatial epidemic models is a topic of considerable recent interest. In another hand, the advances in computer technology have stimulated the development of stochastic models. Metapopulation models are spatial designs that involve movements of individuals between distinct subpopulations. The purpose of the present work has been to develop stochastic models in order to study the transmission dynamics and control of infectious diseases in metapopulations. The authors studied Susceptible-Infected-Susceptible (SIS) and Susceptible-lnfected-Recovered (SIR) epidemic schemes, using the Gillespie algorithm, Computational numerical simulations were carried in order to explore the models. The results obtained show how the dynamics of transmission and the application of control measures within each subpopulation may affect all subpopulations of the system. They also show how the distribution of control measures among subpopulations affects the efficacy of these strategies. The dynamics of the stochastic models developed in the current study follow the trends observed in the classic deterministic designs. Also, the present models exhibit fluctuating behavior. This work highlights the importance of the spatial distribution of the population in spread and control of infectious diseases. In addition, it shows how chance could play an important role in these scenarios.