In order to solve three kinds of fuzzy programming models, i.e. fuzzy expected value model, fuzzy chance-constrained programming model, and fuzzy dependent-chance pro-gramming model, a simultaneous perturbation stocha...In order to solve three kinds of fuzzy programming models, i.e. fuzzy expected value model, fuzzy chance-constrained programming model, and fuzzy dependent-chance pro-gramming model, a simultaneous perturbation stochastic approximation algorithm is proposed by integrating neural network with fuzzy simulation. At first, fuzzy simulation is used to generate a set of input-output data. Then a neural network is trained according to the set. Finally, the trained neural network is embedded in simultaneous perturbation stochastic approximation algorithm. Simultaneous perturbation stochastic approximation algorithm is used to search the optimal solution. Two numerical examples are presented to illustrate the effectiveness of the proposed algorithm.展开更多
At present,simultaneous localization and mapping(SLAM) for an autonomous underwater vehicle(AUV)is a research hotspot.Aiming at the problem of non-linear model and non-Gaussian noise in AUV motion,an improved method o...At present,simultaneous localization and mapping(SLAM) for an autonomous underwater vehicle(AUV)is a research hotspot.Aiming at the problem of non-linear model and non-Gaussian noise in AUV motion,an improved method of variance reduction fast simultaneous localization and mapping(FastSLAM) with simulated annealing is proposed to solve the problems of particle degradation,particle depletion and particle loss in traditional FastSLAM,which lead to the reduction of AUV location estimation accuracy.The adaptive exponential fading factor is generated by the anneal function of simulated annealing algorithm to improve the effective particle number and replace resampling.By increasing the weight of small particles and decreasing the weight of large particles,the variance of particle weight can be reduced,the number of effective particles can be increased,and the accuracy of AUV location and feature location estimation can be improved to some extent by retaining more information carried by particles.The experimental results based on trial data show that the proposed simulated annealing variance reduction FastSLAM method avoids particle degradation,maintains the diversity of particles,weakened the degeneracy and improves the accuracy and stability of AUV navigation and localization system.展开更多
Upscaling of primary geological models with huge cells, especially in porous media, is the first step in fluid flow simulation. Numerical methods are often used to solve the models. The upscaling method must preserve ...Upscaling of primary geological models with huge cells, especially in porous media, is the first step in fluid flow simulation. Numerical methods are often used to solve the models. The upscaling method must preserve the important properties of the spatial distribution of the reservoir properties. An grid upscaling method based on adaptive bandwidth in kernel function is proposed according to the spatial distribution of property. This type of upscaling reduces the number of cells, while preserves the main heterogeneity features of the original fine model. The key point of the paper is upscaling two reservoir properties simultaneously. For each reservoir feature, the amount of bandwidth or optimal threshold is calculated and the results of the upscaling are obtained. Then two approaches are used to upscaling two properties simultaneously based on maximum bandwidth and minimum bandwidth. In fact, we now have a finalized upscaled model for both reservoir properties for each approach in which not only the number of their cells, but also the locations of the cells are equal. The upscaling error of the minimum bandwidth approach is less than that of the maximum bandwidth approach.展开更多
In the present paper,two-and three-dimensional numerical simulations of the flow interference between four cylinders in an in-line square arrangement at Re = 200 are performed.Assisted with the two-dimensional(2-D) nu...In the present paper,two-and three-dimensional numerical simulations of the flow interference between four cylinders in an in-line square arrangement at Re = 200 are performed.Assisted with the two-dimensional(2-D) numerical simulation,the mean and fluctuating forces,Strouhal number(St) and vortex shedding pattern in the wake for each cylinder were analyzed with the spacing ratio(L /D) ranging from 1.5 to 6.0.It was found that,four different vortex modes(viz.,flip-flopping,shielding anti-phase-synchronized,in-phasesynchronized and anti-phase-synchronized) gradually appear with the increase of the L/D ratio.The average drag coefficient of the upstream cylinders is larger than that of the downstream cylinders,while the downstream cylinders usually undergo serious fluctuating forces.When the L/D ratio ranges from 3.0 to 4.0,the dominant frequency of the drag coefficient is equal to the value of St of upstream cylinders.This indicates that a simultaneous resonance in the in-flow and cross-flow directions may occur for some single structures of a multi-body oscillating system.For the 3-D numerical simulation,the L/D and aspect ratios are kept constant as 5.0 and 10,respectively.It was found that some vortices are formed in the wake of the upstream cylinders.Besides,with the same spacing ratio,the calculated drag coefficient and lift coefficient fluctuation are slightly larger than the 2-D results,but with a phase difference.展开更多
The simulation of a high-temperature gas-cooled reactor pebble-bed module(HTR-PM) plant is discussed.This lumped parameter model has the form of a set differential algebraic equations(DAEs) that include stiff equation...The simulation of a high-temperature gas-cooled reactor pebble-bed module(HTR-PM) plant is discussed.This lumped parameter model has the form of a set differential algebraic equations(DAEs) that include stiff equations to model point neutron kinetics.The nested approach is the most common method to solve DAE,but this approach is very expensive and time-consuming due to inner iterations.This paper deals with an alternative approach in which a simultaneous solution method is used.The DAEs are discretized over a time horizon using collocation on finite elements,and Radau collocation points are applied.The resulting nonlinear algebraic equations can be solved by existing solvers.The discrete algorithm is discussed in detail;both accuracy and stability issues are considered.Finally,the simulation results are presented to validate the efficiency and accuracy of the simultaneous approach that takes much less time than the nested one.展开更多
The velocity distribution of layers from surface wave dispersion curve is a severely nonlinear program. Base on the Metropolis rule,we improved the simulated annealing algorithm to simultaneously inverse the velocitie...The velocity distribution of layers from surface wave dispersion curve is a severely nonlinear program. Base on the Metropolis rule,we improved the simulated annealing algorithm to simultaneously inverse the velocities and thicknesses using the dispersion data and identified the Moho and the bottom of lithosphere. The application to the numerical examples with 5% noise shows the velocity RMS is 1. 56% between the non-linear results and the original models when the condition of selecting method for temperature parameters and initial temperature are satisfied. Using the pure dispersions of Rayleigh wave,the nonlinear inversion has been carried out for S-wave velocities and thicknesses of the vertical profile crossing the Indian Plate,the Qinghai-Tibetan Plateau,and the Tarim Basin. It indicated that the crustal thickness is about 70 km in the Qiangtang block,while in the hinterland of the Qinghai-Tibetan Plateau the lithosphere is relatively thin(~ 130 km)from the velocity values and their offsets.展开更多
We simulate the bond and site percolation models on several three-dimensional lattices, including the diamond, body-centered cubic, and face-centered cubic lattices. As on the simple-cubic lattice [Phys. Rev. E, 2013,...We simulate the bond and site percolation models on several three-dimensional lattices, including the diamond, body-centered cubic, and face-centered cubic lattices. As on the simple-cubic lattice [Phys. Rev. E, 2013, 87(5): 052107], it is observed that in comparison with dimensionless ratios based on cluster-size distribution, certain wrapping probabilities exhibit weaker finite-size corrections and are more sensitive to the deviation from percolation threshold Pc, and thus provide a powerful means for determining Pc. We analyze the numerical data of the wrapping probabilities simultaneously such that universal parameters are shared by the aforementioned models, and thus significantly improved estimates of Pc are obtained.展开更多
基金National Natural Science Foundation of China (No.70471049)China Postdoctoral Science Foundation (No. 20060400704)
文摘In order to solve three kinds of fuzzy programming models, i.e. fuzzy expected value model, fuzzy chance-constrained programming model, and fuzzy dependent-chance pro-gramming model, a simultaneous perturbation stochastic approximation algorithm is proposed by integrating neural network with fuzzy simulation. At first, fuzzy simulation is used to generate a set of input-output data. Then a neural network is trained according to the set. Finally, the trained neural network is embedded in simultaneous perturbation stochastic approximation algorithm. Simultaneous perturbation stochastic approximation algorithm is used to search the optimal solution. Two numerical examples are presented to illustrate the effectiveness of the proposed algorithm.
基金supported by the National Science Fund of China under Grants 61603034China Postdoctoral Science Foundation under Grant 2019M653870XB+1 种基金Beijing Municipal Natural Science Foundation (3182027)Fundamental Research Funds for the Central Universities,China,FRF-GF-17-B44,and XJS191315
文摘At present,simultaneous localization and mapping(SLAM) for an autonomous underwater vehicle(AUV)is a research hotspot.Aiming at the problem of non-linear model and non-Gaussian noise in AUV motion,an improved method of variance reduction fast simultaneous localization and mapping(FastSLAM) with simulated annealing is proposed to solve the problems of particle degradation,particle depletion and particle loss in traditional FastSLAM,which lead to the reduction of AUV location estimation accuracy.The adaptive exponential fading factor is generated by the anneal function of simulated annealing algorithm to improve the effective particle number and replace resampling.By increasing the weight of small particles and decreasing the weight of large particles,the variance of particle weight can be reduced,the number of effective particles can be increased,and the accuracy of AUV location and feature location estimation can be improved to some extent by retaining more information carried by particles.The experimental results based on trial data show that the proposed simulated annealing variance reduction FastSLAM method avoids particle degradation,maintains the diversity of particles,weakened the degeneracy and improves the accuracy and stability of AUV navigation and localization system.
文摘Upscaling of primary geological models with huge cells, especially in porous media, is the first step in fluid flow simulation. Numerical methods are often used to solve the models. The upscaling method must preserve the important properties of the spatial distribution of the reservoir properties. An grid upscaling method based on adaptive bandwidth in kernel function is proposed according to the spatial distribution of property. This type of upscaling reduces the number of cells, while preserves the main heterogeneity features of the original fine model. The key point of the paper is upscaling two reservoir properties simultaneously. For each reservoir feature, the amount of bandwidth or optimal threshold is calculated and the results of the upscaling are obtained. Then two approaches are used to upscaling two properties simultaneously based on maximum bandwidth and minimum bandwidth. In fact, we now have a finalized upscaled model for both reservoir properties for each approach in which not only the number of their cells, but also the locations of the cells are equal. The upscaling error of the minimum bandwidth approach is less than that of the maximum bandwidth approach.
基金Sponsored by the National Natural Science Foundation of China(Grant No.50538020)
文摘In the present paper,two-and three-dimensional numerical simulations of the flow interference between four cylinders in an in-line square arrangement at Re = 200 are performed.Assisted with the two-dimensional(2-D) numerical simulation,the mean and fluctuating forces,Strouhal number(St) and vortex shedding pattern in the wake for each cylinder were analyzed with the spacing ratio(L /D) ranging from 1.5 to 6.0.It was found that,four different vortex modes(viz.,flip-flopping,shielding anti-phase-synchronized,in-phasesynchronized and anti-phase-synchronized) gradually appear with the increase of the L/D ratio.The average drag coefficient of the upstream cylinders is larger than that of the downstream cylinders,while the downstream cylinders usually undergo serious fluctuating forces.When the L/D ratio ranges from 3.0 to 4.0,the dominant frequency of the drag coefficient is equal to the value of St of upstream cylinders.This indicates that a simultaneous resonance in the in-flow and cross-flow directions may occur for some single structures of a multi-body oscillating system.For the 3-D numerical simulation,the L/D and aspect ratios are kept constant as 5.0 and 10,respectively.It was found that some vortices are formed in the wake of the upstream cylinders.Besides,with the same spacing ratio,the calculated drag coefficient and lift coefficient fluctuation are slightly larger than the 2-D results,but with a phase difference.
基金Project supported by the National Basic Research Program of China (No. 2009CB320603)the National Natural Science Foundation of China (Nos. 60974007 and 60934007)
文摘The simulation of a high-temperature gas-cooled reactor pebble-bed module(HTR-PM) plant is discussed.This lumped parameter model has the form of a set differential algebraic equations(DAEs) that include stiff equations to model point neutron kinetics.The nested approach is the most common method to solve DAE,but this approach is very expensive and time-consuming due to inner iterations.This paper deals with an alternative approach in which a simultaneous solution method is used.The DAEs are discretized over a time horizon using collocation on finite elements,and Radau collocation points are applied.The resulting nonlinear algebraic equations can be solved by existing solvers.The discrete algorithm is discussed in detail;both accuracy and stability issues are considered.Finally,the simulation results are presented to validate the efficiency and accuracy of the simultaneous approach that takes much less time than the nested one.
基金sponsored by the National Natural Science Foundation of China(41774069,41504047 and 41604054)
文摘The velocity distribution of layers from surface wave dispersion curve is a severely nonlinear program. Base on the Metropolis rule,we improved the simulated annealing algorithm to simultaneously inverse the velocities and thicknesses using the dispersion data and identified the Moho and the bottom of lithosphere. The application to the numerical examples with 5% noise shows the velocity RMS is 1. 56% between the non-linear results and the original models when the condition of selecting method for temperature parameters and initial temperature are satisfied. Using the pure dispersions of Rayleigh wave,the nonlinear inversion has been carried out for S-wave velocities and thicknesses of the vertical profile crossing the Indian Plate,the Qinghai-Tibetan Plateau,and the Tarim Basin. It indicated that the crustal thickness is about 70 km in the Qiangtang block,while in the hinterland of the Qinghai-Tibetan Plateau the lithosphere is relatively thin(~ 130 km)from the velocity values and their offsets.
文摘We simulate the bond and site percolation models on several three-dimensional lattices, including the diamond, body-centered cubic, and face-centered cubic lattices. As on the simple-cubic lattice [Phys. Rev. E, 2013, 87(5): 052107], it is observed that in comparison with dimensionless ratios based on cluster-size distribution, certain wrapping probabilities exhibit weaker finite-size corrections and are more sensitive to the deviation from percolation threshold Pc, and thus provide a powerful means for determining Pc. We analyze the numerical data of the wrapping probabilities simultaneously such that universal parameters are shared by the aforementioned models, and thus significantly improved estimates of Pc are obtained.