This paper analyzes the dynamic characteristics of the variations of the beach volumes for three level zonesof the Yanjing Beach in the Shuidong Bay of the western Guangdong Province by using the methods of dynamic sy...This paper analyzes the dynamic characteristics of the variations of the beach volumes for three level zonesof the Yanjing Beach in the Shuidong Bay of the western Guangdong Province by using the methods of dynamic systemanalysis and the multi-dimensional spectral estimation. The results show that the variations of the beach volume arecharaCterized by the multiband oscillations with a dominant semimonth period. Upwards the low tide level, the beachtends to be stable. The estimates of the partial coherences and the partial phases indicate that the variations of thebeach volumes are mainly the results of the direct actions of the waves which are influenced by the tidal level changesand driven by the wind stress. The simulation results of the beach volume series for different beach heart zones bythreshold mixed regressive models indicate that the influence of the tide on the variations of the beach volumes is weakened and the direct actions of the wave energy and the wind stress are apparently enhanced with the increase of thebeach height.(This project was supported by the National Natural Science Foundation of China.)展开更多
Unexpected noise in reservoir stochastic simulation realization may be too high to make the realization useful, especially when there is a lack of hard data. Through discussing the uncertainties, we present two ways t...Unexpected noise in reservoir stochastic simulation realization may be too high to make the realization useful, especially when there is a lack of hard data. Through discussing the uncertainties, we present two ways to control the uncertainty ratio that is brought by the algorithm of stochastic simulation. By reasonably reducing the random value of the stochastic simulation result, the unexpected values introduced by the residual that associates with random series can be controlled. Another way when the data disperse unevenly is to control the stochastic simulation order by grouping the points that need to be simulated to make those points which can be simulated by more neighborhood hard data calculated first. Both methods do not go against the core stochastic simulation algorithm.展开更多
Optimization of long-term mine production scheduling in open pit mines deals with the management of cash flows, typically in the order of hundreds of millions of dollars. Conventional mine scheduling utilizes optimiza...Optimization of long-term mine production scheduling in open pit mines deals with the management of cash flows, typically in the order of hundreds of millions of dollars. Conventional mine scheduling utilizes optimization methods that are not capable of accounting for inherent technical uncertainties such as uncertainty in the expected ore/metal supply from the underground, acknowledged to be the most critical factor. To integrate ore/metal uncertainty into the optimization of mine production scheduling a stochastic integer programming(SIP) formulation is tested at a copper deposit. The stochastic solution maximizes the economic value of a project and minimizes deviations from production targets in the presence of ore/metal uncertainty. Unlike the conventional approach, the SIP model accounts and manages risk in ore supply, leading to a mine production schedule with a 29% higher net present value than the schedule obtained from the conventional, industry-standard optimization approach, thus contributing to improving the management and sustainable utilization of mineral resources.展开更多
The stochastic finite-fault simulation method was applied to synthesize the horizontal ground acceleration seismograms produced by the MW6.1 Ludian earthquake on August 3,2014.For this purpose,we produced first a tota...The stochastic finite-fault simulation method was applied to synthesize the horizontal ground acceleration seismograms produced by the MW6.1 Ludian earthquake on August 3,2014.For this purpose,we produced first a total of 200 kinematic source models for the Ludian event,which are characterized by the heterogeneous slip on the conjugated ruptured fault and the slip-dependent spreading of the rupture front.The results indicated that the heterogeneous slip and the spatial extent of the ruptured fault play dominant roles in the spatial distribution of ground motions in the near-fault area.The peak ground accelerations(PGAs)and 5%-damped pseudospectral accelerations(PSAs)at periods shorter than 0.5 s estimated on the resulting synthetics generally match well with the observations at stations with Joyner-Boore distances(RJB)greater than 20 km.The synthetic PGVs and PSAs at periods of 0.5 s and 0.75 s are in good agreement with predicted medians by the Yu14 model(Yu et al.,2014).However,the synthetic results are generally much lower than the predicted medians by BSSA14 model(Boore et al.,2014).Moreover,the ground motion variability caused by the randomness in the source rupture process was evaluated by these synthetics.The standard deviations of PSAs on the base-10 logarithmic scale,Sigma[log10(PSA)],are closely dependent on either the spectral period or the RJB.The Sigma[log10(PSA)]remains a constant approximately 0.55 at periods shorter than 0.1 s,and then increase continuously up to^0.13 as the period increases from 0.1 to 2.0 s.The Sigma[log10(PSA)]values at periods of 0.1‒2.0 s show the downward tendency as the RJB values increase.However,the Sigma[log10(PSA)]values at periods shorter than 0.1 s decrease as the RJB values increase up to^50 km,and then increase with the increasing RJB.Furthermore,we found that the ground-motion variability shows the significant dependence on the azimuth.展开更多
Copula functions have been widely used in stochastic simulation and prediction of streamflow.However,existing models are usually limited to single two-dimensional or three-dimensional copulas with the same bivariate b...Copula functions have been widely used in stochastic simulation and prediction of streamflow.However,existing models are usually limited to single two-dimensional or three-dimensional copulas with the same bivariate block for all months.To address this limitation,this study developed a mixed D-vine copula-based conditional quantile model that can capture temporal correlations.This model can generate streamflow by selecting different historical streamflow variables as the conditions for different months and by exploiting the conditional quantile functions of streamflows in different months with mixed D-vine copulas.The up-to-down sequential method,which couples the maximum weight approach with the Akaike information criteria and the maximum likelihood approach,was used to determine the structures of multivariate Dvine copulas.The developed model was used in a case study to synthesize the monthly streamflow at the Tangnaihai hydrological station,the inflow control station of the Longyangxia Reservoir in the Yellow River Basin.The results showed that the developed model outperformed the commonly used bivariate copula model in terms of the performance in simulating the seasonality and interannual variability of streamflow.This model provides useful information for water-related natural hazard risk assessment and integrated water resources management and utilization.展开更多
The stochastic simulation algorithm （SSA） accurately depicts spatially homogeneous wellstirred chemically reacting systems with small populations of chemical species and properly represents noise, but it is often ab...The stochastic simulation algorithm （SSA） accurately depicts spatially homogeneous wellstirred chemically reacting systems with small populations of chemical species and properly represents noise, but it is often abandoned when modeling larger systems because of its computational complexity. In this work, a twin support vector regression based stochastic simulations algorithm （TS^3A） is proposed by combining the twin support vector regression and SSA, the former is a well-known robust regression method in machine learning. Numerical results indicate that this proposed algorithm can be applied to a wide range of chemically reacting systems and obtain significant improvements on efficiency and accuracy with fewer simulating runs over the existing methods.展开更多
Stochastic finite-fault simulations are effective for simulating ground motions and are widely used in engineering to determine the impacts of ground motion and develop relevant predictive equations.In this study,the ...Stochastic finite-fault simulations are effective for simulating ground motions and are widely used in engineering to determine the impacts of ground motion and develop relevant predictive equations.In this study,the source,path,and site amplification coefficient of western Sichuan Province,China,and stochastic finite-fault simulations were used to simulate the acceleration time series,Fourier amplitude spectra,and 5%damped response spectra of 28 strong-motion stations with rupture distances within 300 km of the 2022 MS6.8 Luding earthquake.The simulation results of 14 stations at rupture distances of 45-185 km match the observation.However,the simulation results of 3 near-and 6 far-field stations at rupture distances of 12-36 km and 222-286 km,respectively,were obviously deviated from the observations.Simulation results of the near-field stations are larger than the observations at high frequencies(>6 Hz).The discrepancy likely comes from the nonlinear site effect of near-field stations,which reduced the site amplification at high frequencies.Simulation result of the far-field stations is smaller than the observation at frequencies above 1 Hz.As these stations are located close to the Longmenshan Fault Zone(LFZ),thus,we obtained a new quality factor(Q)from data of historical events and stations located around LFZ.Using the new Q value,the discrepancies of the high-frequency simulation results of the far-field stations were corrected.This result indicated that the laterally varying Q values can be used to address the impact of strong crustal lateral heterogeneity on simulation.展开更多
Presented here is an L-leap method for accelerating stochastic simulation of well-stirred chemically reacting systems, in which the number of reactions occurring in a reaction channel with the largest propensity funct...Presented here is an L-leap method for accelerating stochastic simulation of well-stirred chemically reacting systems, in which the number of reactions occurring in a reaction channel with the largest propensity function is calculated from the leap condition and the number of reactions occurring in the other reaction channels are generated by using binomial random variables during a leap. The L-leap method can better satisfy the leap condition. Numerical simulation results indicate that the L-leap method can obtain better performance than established methods.展开更多
To reveal the low growth rate of Acidithiobacillus ferrooxidans, a stochastic growth model was proposed to analyze growth curves of these bacteria in a batch culture. An algorithm was applied to simulate the bacteria ...To reveal the low growth rate of Acidithiobacillus ferrooxidans, a stochastic growth model was proposed to analyze growth curves of these bacteria in a batch culture. An algorithm was applied to simulate the bacteria population during lag and exponential phase. The results show that the model moderately fits the experimental data. Further, the mean growth constant （K） of growth curves is obtained by fitting the logarithm of the simulating population data versus the generation numbers with the different initial population number （N0） and initial mean activity of population （A0）. When No is 300 and 700 respectively, the discrepancy of K value is only 0.91%, however, A0 is 0.34 and 0.38 respectively, the discrepancy of K value is 19.53%. It suggests that the effect of A0 on the lag phase exceeds No, though both parameters could shorten the lag phase by increasing their values.展开更多
In this paper, we develop a modified accelerated stochastic simulation method for chemically reacting systems, called the ＂final all possible steps＂ （FAPS） method, which obtains the reliable statistics of all spec...In this paper, we develop a modified accelerated stochastic simulation method for chemically reacting systems, called the ＂final all possible steps＂ （FAPS） method, which obtains the reliable statistics of all species in any time during the time course with fewer simulation times. Moreover, the FAPS method can be incorporated into the leap methods, which makes the simulation of larger systems more efficient. Numerical results indicate that the proposed methods can be applied to a wide range of chemically reacting systems with a high-precision level and obtain a significant improvement on efficiency over the existing methods.展开更多
In this article, we consider the construction of a SVIR (Susceptible, Vaccinated, Infected, Recovered) stochastic compartmental model of measles. We prove that the deterministic solution is asymptotically the average ...In this article, we consider the construction of a SVIR (Susceptible, Vaccinated, Infected, Recovered) stochastic compartmental model of measles. We prove that the deterministic solution is asymptotically the average of the stochastic solution in the case of small population size. The choice of this model takes into account the random fluctuations inherent to the epidemiological characteristics of rural populations of Niger, notably a high prevalence of measles in children under 5, coupled with a very low immunization coverage.展开更多
Biochemical systems have numerous practical applications, in particular to the study of critical intracellular processes. Frequently, biochemical kinetic models depict cellular processes as systems of chemical reactio...Biochemical systems have numerous practical applications, in particular to the study of critical intracellular processes. Frequently, biochemical kinetic models depict cellular processes as systems of chemical reactions. Many biological processes in a cell are inherently stochastic, due to the existence of some low molecular amounts. These stochastic fluctuations may have a great effect on the biochemical system’s behaviour. In such cases, stochastic models are necessary to accurately describe the system’s dynamics. Biochemical systems at the cellular level may entail many species or reactions and their mathematical models may be non-linear and with multiple scales in time. In this work, we provide a numerical technique for simplifying stochastic discrete models of well-stirred biochemical systems, which ensures that the main properties of the original system are preserved. The proposed technique employs sensitivity analysis and requires solving an optimization problem. The numerical tests on several models of practical interest show that our model reduction strategy performs very well.展开更多
In this paper, we analyze the quasi-stationary distribution of the stochastic <em>SVIR</em> (Susceptible, Vaccinated, Infected, Recovered) model for the measles. The quasi-stationary distributions, as disc...In this paper, we analyze the quasi-stationary distribution of the stochastic <em>SVIR</em> (Susceptible, Vaccinated, Infected, Recovered) model for the measles. The quasi-stationary distributions, as discussed by Danoch and Seneta, have been used in biology to describe the steady state behaviour of population models which exhibit discernible stationarity before to become extinct. The stochastic <em>SVIR</em> model is a stochastic <em>SIR</em> (Susceptible, Infected, Recovered) model with vaccination and recruitment where the disease-free equilibrium is reached, regardless of the magnitude of the basic reproduction number. But the mean time until the absorption (the disease-free) can be very long. If we assume the effective reproduction number <em>R</em><em><sub>p</sub></em> < 1 or <img src="Edit_67da0b97-83f9-42ef-8a00-a13da2d59963.bmp" alt="" />, the quasi-stationary distribution can be closely approximated by geometric distribution. <em>β</em> and <em>δ</em> stands respectively, for the disease transmission coefficient and the natural rate.展开更多
The present study proposes a stochastic simulation scheme to model reactive boundaries through a position jump process which can be readily implemented into the Inhomogeneous Stochastic Simulation Algorithm by modifyi...The present study proposes a stochastic simulation scheme to model reactive boundaries through a position jump process which can be readily implemented into the Inhomogeneous Stochastic Simulation Algorithm by modifying the propensity of the diffusive jump over the reactive boundary. As compared to the literature, the present approach does not require any correction factors for the propensity. Also, the current expression relaxes the constraint on the compartment size allowing the problem to be solved with a coarser grid and therefore saves considerable computational cost. The modified algorithm is then applied to simulate three reaction-diffusion systems with reactive boundaries.展开更多
The auto-correlation function and the cross-correlation of an autonomous stochastic system with nonlinear time-delayed feedback are investigated by using the stochastic simulation method. There are prominent differenc...The auto-correlation function and the cross-correlation of an autonomous stochastic system with nonlinear time-delayed feedback are investigated by using the stochastic simulation method. There are prominent differences be- tween the roles of quadratic time-delayed feedback and cubic time-delayed feedback on the correlations of an autonomous stochastic system. Under quadratic time-delayed feedback, the nonlinear time delay fails to improve the noisy state of the autonomous stochastic system, the auto-correlation decreases monotonously to zero, and the cross-correlation increases monotonously to zero with the decay time. Under cubic time-delayed feedback, the nonlinear time delay can improve the noisy state of the autonomous stochastic system; the auto-correlation and the cross-correlation show periodical oscillation and attenuation, finally tending to zero with the decay time. Comparing the correlations of the system between with nonfinear time-delayed feedback and linear time-delayed feedback, we find that nonlinear time-delayed feedback lowers the correlation strength of the autonomous stochastic system.展开更多
Many difficult engineering problems cannot be solved by the conventional optimization techniques in practice. Direct searches that need no recourse to explicit derivatives are revived and become popular since the new ...Many difficult engineering problems cannot be solved by the conventional optimization techniques in practice. Direct searches that need no recourse to explicit derivatives are revived and become popular since the new century. In order to get a deep insight into this field, some notes on the direct searches for non-smooth optimization problems are made. The global convergence vs. local convergence and their influences on expected solutions for simulation-based stochastic optimization are pointed out. The sufficient and simple decrease criteria for step acceptance are analyzed, and why simple decrease is enough for globalization in direct searches is identified. The reason to introduce the positive spanning set and its usage in direct searches is explained. Other topics such as the generalization of direct searches to bound, linear and non-linear constraints are also briefly discussed.展开更多
Currently,energy conservation draws wide attention in industrial manufacturing systems.In recent years,many studies have aimed at saving energy consumption in the process of manufacturing and scheduling is regarded as...Currently,energy conservation draws wide attention in industrial manufacturing systems.In recent years,many studies have aimed at saving energy consumption in the process of manufacturing and scheduling is regarded as an effective approach.This paper puts forwards a multi-objective stochastic parallel machine scheduling problem with the consideration of deteriorating and learning effects.In it,the real processing time of jobs is calculated by using their processing speed and normal processing time.To describe this problem in a mathematical way,amultiobjective stochastic programming model aiming at realizing makespan and energy consumption minimization is formulated.Furthermore,we develop a multi-objective multi-verse optimization combined with a stochastic simulation method to deal with it.In this approach,the multi-verse optimization is adopted to find favorable solutions from the huge solution domain,while the stochastic simulation method is employed to assess them.By conducting comparison experiments on test problems,it can be verified that the developed approach has better performance in coping with the considered problem,compared to two classic multi-objective evolutionary algorithms.展开更多
Uncertainty in geological structural modeling, especially geological corrosion(a kind of karst cave), is a bottleneck that restricts the development and application of geological computer modeling and effect estimatio...Uncertainty in geological structural modeling, especially geological corrosion(a kind of karst cave), is a bottleneck that restricts the development and application of geological computer modeling and effect estimation. To solve this issue, a stochastic modeling method based on the random field theory is proposed in comparison with the deterministic geometric modeling method. Then the constraint random field modeling method and the random field modeling method without constrained parameters are compared and analyzed. A case study shows that the novel stochastic simulation method is an effective tool to describe the distribution characteristics of corrosion parameters and reflect the updated geological prospecting information. The influence of geological corrosion on the dam behavior can also be better analyzed by using the stochastic simulation method. At the same time, the unconfined random field ignores the sample location information and may lead to higher variability. Therefore, the constraint random field modeling method can provide a useful reference for the numerical analysis under complex geological conditions.展开更多
We report stochastic simulations of the yeast mating signal transduction pathway. The effects of intrinsic and external noise, the influence of cell-to-cell difference in the pathway capacity, and noise propagation in...We report stochastic simulations of the yeast mating signal transduction pathway. The effects of intrinsic and external noise, the influence of cell-to-cell difference in the pathway capacity, and noise propagation in the pathway have been examined. The stochastic temporal behaviour of the pathway is found to be robust to the influence of inherent fluctuations, and intrinsic noise propagates in the pathway in a uniform pattern when the yeasts are treated with pheromones of different stimulus strengths and of varied fluctuations. In agreement with recent experimental findings, extrinsic noise is found to play a more prominent role than intrinsic noise in the variability of proteins. The occurrence frequency for the reactions in the pathway are also examined and a more compact network is obtained by dropping most of the reactions of least occurrence.展开更多
The stochastic simulation method, based on the concept of control window and the numerical solution of the Langevin equation, is applied to solve the deposition problem of particles from the flowing suspensions onto a...The stochastic simulation method, based on the concept of control window and the numerical solution of the Langevin equation, is applied to solve the deposition problem of particles from the flowing suspensions onto a fiber collector. Using the Kuwabara model to characterize the flow field, the effects of Stokes number, interception parameter, packing density, particle size distribution on the collection efficioncy, and the deposition morphology of particles onto a collector are i examined. The morphology of deposit obtained in the simulated results agrees Well with experimental observations. The estimation of the initial coUection efficiency through the simulations considers that the deposited particles are in good agreement with published experimental data. In addition, the collection efficiency of particles increases in a wider particle size distribution region.展开更多
文摘This paper analyzes the dynamic characteristics of the variations of the beach volumes for three level zonesof the Yanjing Beach in the Shuidong Bay of the western Guangdong Province by using the methods of dynamic systemanalysis and the multi-dimensional spectral estimation. The results show that the variations of the beach volume arecharaCterized by the multiband oscillations with a dominant semimonth period. Upwards the low tide level, the beachtends to be stable. The estimates of the partial coherences and the partial phases indicate that the variations of thebeach volumes are mainly the results of the direct actions of the waves which are influenced by the tidal level changesand driven by the wind stress. The simulation results of the beach volume series for different beach heart zones bythreshold mixed regressive models indicate that the influence of the tide on the variations of the beach volumes is weakened and the direct actions of the wave energy and the wind stress are apparently enhanced with the increase of thebeach height.(This project was supported by the National Natural Science Foundation of China.)
文摘Unexpected noise in reservoir stochastic simulation realization may be too high to make the realization useful, especially when there is a lack of hard data. Through discussing the uncertainties, we present two ways to control the uncertainty ratio that is brought by the algorithm of stochastic simulation. By reasonably reducing the random value of the stochastic simulation result, the unexpected values introduced by the residual that associates with random series can be controlled. Another way when the data disperse unevenly is to control the stochastic simulation order by grouping the points that need to be simulated to make those points which can be simulated by more neighborhood hard data calculated first. Both methods do not go against the core stochastic simulation algorithm.
基金funded from the National Science and Engineering Research Council of Canada,Collaborative R&D Grant CRDPJ 335696 with BHP Billiton and NSERC Discovery Grant 239019 to R. Dimitrakopoulos
文摘Optimization of long-term mine production scheduling in open pit mines deals with the management of cash flows, typically in the order of hundreds of millions of dollars. Conventional mine scheduling utilizes optimization methods that are not capable of accounting for inherent technical uncertainties such as uncertainty in the expected ore/metal supply from the underground, acknowledged to be the most critical factor. To integrate ore/metal uncertainty into the optimization of mine production scheduling a stochastic integer programming(SIP) formulation is tested at a copper deposit. The stochastic solution maximizes the economic value of a project and minimizes deviations from production targets in the presence of ore/metal uncertainty. Unlike the conventional approach, the SIP model accounts and manages risk in ore supply, leading to a mine production schedule with a 29% higher net present value than the schedule obtained from the conventional, industry-standard optimization approach, thus contributing to improving the management and sustainable utilization of mineral resources.
基金supported by the Science Foundation of the Institute of Engineering Mechanics,China Earthquake Administration(No.2018B03)National Natural Science Foundation of China(No.51808514).
文摘The stochastic finite-fault simulation method was applied to synthesize the horizontal ground acceleration seismograms produced by the MW6.1 Ludian earthquake on August 3,2014.For this purpose,we produced first a total of 200 kinematic source models for the Ludian event,which are characterized by the heterogeneous slip on the conjugated ruptured fault and the slip-dependent spreading of the rupture front.The results indicated that the heterogeneous slip and the spatial extent of the ruptured fault play dominant roles in the spatial distribution of ground motions in the near-fault area.The peak ground accelerations(PGAs)and 5%-damped pseudospectral accelerations(PSAs)at periods shorter than 0.5 s estimated on the resulting synthetics generally match well with the observations at stations with Joyner-Boore distances(RJB)greater than 20 km.The synthetic PGVs and PSAs at periods of 0.5 s and 0.75 s are in good agreement with predicted medians by the Yu14 model(Yu et al.,2014).However,the synthetic results are generally much lower than the predicted medians by BSSA14 model(Boore et al.,2014).Moreover,the ground motion variability caused by the randomness in the source rupture process was evaluated by these synthetics.The standard deviations of PSAs on the base-10 logarithmic scale,Sigma[log10(PSA)],are closely dependent on either the spectral period or the RJB.The Sigma[log10(PSA)]remains a constant approximately 0.55 at periods shorter than 0.1 s,and then increase continuously up to^0.13 as the period increases from 0.1 to 2.0 s.The Sigma[log10(PSA)]values at periods of 0.1‒2.0 s show the downward tendency as the RJB values increase.However,the Sigma[log10(PSA)]values at periods shorter than 0.1 s decrease as the RJB values increase up to^50 km,and then increase with the increasing RJB.Furthermore,we found that the ground-motion variability shows the significant dependence on the azimuth.
基金supported by the National Natural Science Foundation of China(Grant No.52109010)the Postdoctoral Science Foundation of China(Grant No.2021M701047)the China National Postdoctoral Program for Innovative Talents(Grant No.BX20200113).
文摘Copula functions have been widely used in stochastic simulation and prediction of streamflow.However,existing models are usually limited to single two-dimensional or three-dimensional copulas with the same bivariate block for all months.To address this limitation,this study developed a mixed D-vine copula-based conditional quantile model that can capture temporal correlations.This model can generate streamflow by selecting different historical streamflow variables as the conditions for different months and by exploiting the conditional quantile functions of streamflows in different months with mixed D-vine copulas.The up-to-down sequential method,which couples the maximum weight approach with the Akaike information criteria and the maximum likelihood approach,was used to determine the structures of multivariate Dvine copulas.The developed model was used in a case study to synthesize the monthly streamflow at the Tangnaihai hydrological station,the inflow control station of the Longyangxia Reservoir in the Yellow River Basin.The results showed that the developed model outperformed the commonly used bivariate copula model in terms of the performance in simulating the seasonality and interannual variability of streamflow.This model provides useful information for water-related natural hazard risk assessment and integrated water resources management and utilization.
基金This work was supported by the National Natural Science Foundation of China （No.30871341）, the National High-Tech Research and Development Program of China （No.2006AA02-Z190）, the Shanghai Leading Academic Discipline Project （No.S30405）, and the Natural Science Foundation of Shanghai Normal University （No.SK200937）.
文摘The stochastic simulation algorithm （SSA） accurately depicts spatially homogeneous wellstirred chemically reacting systems with small populations of chemical species and properly represents noise, but it is often abandoned when modeling larger systems because of its computational complexity. In this work, a twin support vector regression based stochastic simulations algorithm （TS^3A） is proposed by combining the twin support vector regression and SSA, the former is a well-known robust regression method in machine learning. Numerical results indicate that this proposed algorithm can be applied to a wide range of chemically reacting systems and obtain significant improvements on efficiency and accuracy with fewer simulating runs over the existing methods.
基金supported by the Special Fund of the Institute of Geophysics,China Earthquake Administration(Nos.DQJB2B28 and DQJB22X09)the National Natural Science Foundation of China(No.52192675).
文摘Stochastic finite-fault simulations are effective for simulating ground motions and are widely used in engineering to determine the impacts of ground motion and develop relevant predictive equations.In this study,the source,path,and site amplification coefficient of western Sichuan Province,China,and stochastic finite-fault simulations were used to simulate the acceleration time series,Fourier amplitude spectra,and 5%damped response spectra of 28 strong-motion stations with rupture distances within 300 km of the 2022 MS6.8 Luding earthquake.The simulation results of 14 stations at rupture distances of 45-185 km match the observation.However,the simulation results of 3 near-and 6 far-field stations at rupture distances of 12-36 km and 222-286 km,respectively,were obviously deviated from the observations.Simulation results of the near-field stations are larger than the observations at high frequencies(>6 Hz).The discrepancy likely comes from the nonlinear site effect of near-field stations,which reduced the site amplification at high frequencies.Simulation result of the far-field stations is smaller than the observation at frequencies above 1 Hz.As these stations are located close to the Longmenshan Fault Zone(LFZ),thus,we obtained a new quality factor(Q)from data of historical events and stations located around LFZ.Using the new Q value,the discrepancies of the high-frequency simulation results of the far-field stations were corrected.This result indicated that the laterally varying Q values can be used to address the impact of strong crustal lateral heterogeneity on simulation.
基金Project supported by the National Natural Science Foundation of China (No.30571059)the National High-Tech Research and Development Program of China(No.2006AA02Z190).
文摘Presented here is an L-leap method for accelerating stochastic simulation of well-stirred chemically reacting systems, in which the number of reactions occurring in a reaction channel with the largest propensity function is calculated from the leap condition and the number of reactions occurring in the other reaction channels are generated by using binomial random variables during a leap. The L-leap method can better satisfy the leap condition. Numerical simulation results indicate that the L-leap method can obtain better performance than established methods.
基金Project(50321402) supported by the Science Fund for Creative Research Groups of China project(2004CB619204) sup-ported by the National Key Fundamental Research Development Programof China
文摘To reveal the low growth rate of Acidithiobacillus ferrooxidans, a stochastic growth model was proposed to analyze growth curves of these bacteria in a batch culture. An algorithm was applied to simulate the bacteria population during lag and exponential phase. The results show that the model moderately fits the experimental data. Further, the mean growth constant （K） of growth curves is obtained by fitting the logarithm of the simulating population data versus the generation numbers with the different initial population number （N0） and initial mean activity of population （A0）. When No is 300 and 700 respectively, the discrepancy of K value is only 0.91%, however, A0 is 0.34 and 0.38 respectively, the discrepancy of K value is 19.53%. It suggests that the effect of A0 on the lag phase exceeds No, though both parameters could shorten the lag phase by increasing their values.
基金the National Natural Science Foundation of China(No.30571059)the National High-Tech Research and Development Program of China(No.2006AA02Z190)
文摘In this paper, we develop a modified accelerated stochastic simulation method for chemically reacting systems, called the ＂final all possible steps＂ （FAPS） method, which obtains the reliable statistics of all species in any time during the time course with fewer simulation times. Moreover, the FAPS method can be incorporated into the leap methods, which makes the simulation of larger systems more efficient. Numerical results indicate that the proposed methods can be applied to a wide range of chemically reacting systems with a high-precision level and obtain a significant improvement on efficiency over the existing methods.
文摘In this article, we consider the construction of a SVIR (Susceptible, Vaccinated, Infected, Recovered) stochastic compartmental model of measles. We prove that the deterministic solution is asymptotically the average of the stochastic solution in the case of small population size. The choice of this model takes into account the random fluctuations inherent to the epidemiological characteristics of rural populations of Niger, notably a high prevalence of measles in children under 5, coupled with a very low immunization coverage.
文摘Biochemical systems have numerous practical applications, in particular to the study of critical intracellular processes. Frequently, biochemical kinetic models depict cellular processes as systems of chemical reactions. Many biological processes in a cell are inherently stochastic, due to the existence of some low molecular amounts. These stochastic fluctuations may have a great effect on the biochemical system’s behaviour. In such cases, stochastic models are necessary to accurately describe the system’s dynamics. Biochemical systems at the cellular level may entail many species or reactions and their mathematical models may be non-linear and with multiple scales in time. In this work, we provide a numerical technique for simplifying stochastic discrete models of well-stirred biochemical systems, which ensures that the main properties of the original system are preserved. The proposed technique employs sensitivity analysis and requires solving an optimization problem. The numerical tests on several models of practical interest show that our model reduction strategy performs very well.
文摘In this paper, we analyze the quasi-stationary distribution of the stochastic <em>SVIR</em> (Susceptible, Vaccinated, Infected, Recovered) model for the measles. The quasi-stationary distributions, as discussed by Danoch and Seneta, have been used in biology to describe the steady state behaviour of population models which exhibit discernible stationarity before to become extinct. The stochastic <em>SVIR</em> model is a stochastic <em>SIR</em> (Susceptible, Infected, Recovered) model with vaccination and recruitment where the disease-free equilibrium is reached, regardless of the magnitude of the basic reproduction number. But the mean time until the absorption (the disease-free) can be very long. If we assume the effective reproduction number <em>R</em><em><sub>p</sub></em> < 1 or <img src="Edit_67da0b97-83f9-42ef-8a00-a13da2d59963.bmp" alt="" />, the quasi-stationary distribution can be closely approximated by geometric distribution. <em>β</em> and <em>δ</em> stands respectively, for the disease transmission coefficient and the natural rate.
文摘The present study proposes a stochastic simulation scheme to model reactive boundaries through a position jump process which can be readily implemented into the Inhomogeneous Stochastic Simulation Algorithm by modifying the propensity of the diffusive jump over the reactive boundary. As compared to the literature, the present approach does not require any correction factors for the propensity. Also, the current expression relaxes the constraint on the compartment size allowing the problem to be solved with a coarser grid and therefore saves considerable computational cost. The modified algorithm is then applied to simulate three reaction-diffusion systems with reactive boundaries.
基金Supported by the National Natural Science Foundation of China under Grant No.11265012Yunnan Province Open Key Laboratory of Mechanics in Colleges and Universities
文摘The auto-correlation function and the cross-correlation of an autonomous stochastic system with nonlinear time-delayed feedback are investigated by using the stochastic simulation method. There are prominent differences be- tween the roles of quadratic time-delayed feedback and cubic time-delayed feedback on the correlations of an autonomous stochastic system. Under quadratic time-delayed feedback, the nonlinear time delay fails to improve the noisy state of the autonomous stochastic system, the auto-correlation decreases monotonously to zero, and the cross-correlation increases monotonously to zero with the decay time. Under cubic time-delayed feedback, the nonlinear time delay can improve the noisy state of the autonomous stochastic system; the auto-correlation and the cross-correlation show periodical oscillation and attenuation, finally tending to zero with the decay time. Comparing the correlations of the system between with nonfinear time-delayed feedback and linear time-delayed feedback, we find that nonlinear time-delayed feedback lowers the correlation strength of the autonomous stochastic system.
基金supported by the Key Foundation of Southwest University for Nationalities(09NZD001).
文摘Many difficult engineering problems cannot be solved by the conventional optimization techniques in practice. Direct searches that need no recourse to explicit derivatives are revived and become popular since the new century. In order to get a deep insight into this field, some notes on the direct searches for non-smooth optimization problems are made. The global convergence vs. local convergence and their influences on expected solutions for simulation-based stochastic optimization are pointed out. The sufficient and simple decrease criteria for step acceptance are analyzed, and why simple decrease is enough for globalization in direct searches is identified. The reason to introduce the positive spanning set and its usage in direct searches is explained. Other topics such as the generalization of direct searches to bound, linear and non-linear constraints are also briefly discussed.
文摘Currently,energy conservation draws wide attention in industrial manufacturing systems.In recent years,many studies have aimed at saving energy consumption in the process of manufacturing and scheduling is regarded as an effective approach.This paper puts forwards a multi-objective stochastic parallel machine scheduling problem with the consideration of deteriorating and learning effects.In it,the real processing time of jobs is calculated by using their processing speed and normal processing time.To describe this problem in a mathematical way,amultiobjective stochastic programming model aiming at realizing makespan and energy consumption minimization is formulated.Furthermore,we develop a multi-objective multi-verse optimization combined with a stochastic simulation method to deal with it.In this approach,the multi-verse optimization is adopted to find favorable solutions from the huge solution domain,while the stochastic simulation method is employed to assess them.By conducting comparison experiments on test problems,it can be verified that the developed approach has better performance in coping with the considered problem,compared to two classic multi-objective evolutionary algorithms.
基金Supported by Tianjin Youth Research Program of Application Foundation and Advanced Technology(No.15JCQNJC08000)the National Natural Science Foundation of China(No.51509182)+1 种基金the Foundation for Innovative Research Groups of the National Natural Science Foundation of China(No.51321065)Open Foundation from State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering,Hohai University(No.2014491211)
文摘Uncertainty in geological structural modeling, especially geological corrosion(a kind of karst cave), is a bottleneck that restricts the development and application of geological computer modeling and effect estimation. To solve this issue, a stochastic modeling method based on the random field theory is proposed in comparison with the deterministic geometric modeling method. Then the constraint random field modeling method and the random field modeling method without constrained parameters are compared and analyzed. A case study shows that the novel stochastic simulation method is an effective tool to describe the distribution characteristics of corrosion parameters and reflect the updated geological prospecting information. The influence of geological corrosion on the dam behavior can also be better analyzed by using the stochastic simulation method. At the same time, the unconfined random field ignores the sample location information and may lead to higher variability. Therefore, the constraint random field modeling method can provide a useful reference for the numerical analysis under complex geological conditions.
基金supported by the National Natural Science Foundation of China(Grant No 10774008)National Key Basic Research Program of China(Grant Nos 2007CB814800 and 2006CB910706)the support of the National Funds for Fostering Talents in Basic Science of China(Grant No J0630311)
文摘We report stochastic simulations of the yeast mating signal transduction pathway. The effects of intrinsic and external noise, the influence of cell-to-cell difference in the pathway capacity, and noise propagation in the pathway have been examined. The stochastic temporal behaviour of the pathway is found to be robust to the influence of inherent fluctuations, and intrinsic noise propagates in the pathway in a uniform pattern when the yeasts are treated with pheromones of different stimulus strengths and of varied fluctuations. In agreement with recent experimental findings, extrinsic noise is found to play a more prominent role than intrinsic noise in the variability of proteins. The occurrence frequency for the reactions in the pathway are also examined and a more compact network is obtained by dropping most of the reactions of least occurrence.
基金Shanghai Leading Academic Discipline Project,China(No.B604)
文摘The stochastic simulation method, based on the concept of control window and the numerical solution of the Langevin equation, is applied to solve the deposition problem of particles from the flowing suspensions onto a fiber collector. Using the Kuwabara model to characterize the flow field, the effects of Stokes number, interception parameter, packing density, particle size distribution on the collection efficioncy, and the deposition morphology of particles onto a collector are i examined. The morphology of deposit obtained in the simulated results agrees Well with experimental observations. The estimation of the initial coUection efficiency through the simulations considers that the deposited particles are in good agreement with published experimental data. In addition, the collection efficiency of particles increases in a wider particle size distribution region.
