To design the optimum acceleration control schedule for the Adaptive Cycle Engine(ACE)in the full flight envelope,this paper establishes a direct simulation model of the ACE transient state.In this model,geometric par...To design the optimum acceleration control schedule for the Adaptive Cycle Engine(ACE)in the full flight envelope,this paper establishes a direct simulation model of the ACE transient state.In this model,geometric parameters are used to replace the component state parameters.The corresponding relationship between geometric parameters and component state parameters is determined by sensitivity analysis.The geometric variables are controlled when the geometric adjustment speed exceeds the limit,and at the same time the corresponding component state parameters are iterated.The gradient optimization algorism is used to optimize the ground acceleration process of ACE,and the control schedule in terms of operating point of compression components and corrected acceleration rate is used as the full-envelope acceleration control schedule based on the similarity principle.The acceleration control schedules of the triple-bypass mode and the double-bypass mode are designed in this paper.The acceleration processes under various flight conditions are simulated using the acceleration control schedules.Compared with the acceleration process with the linear geometric adjustment schedule,the acceleration performance of ACE is improved by the acceleration control schedule,with the impulse of the acceleration process of the triple-bypass mode being increased by 8.7%-12.3% and the impulse of the double-bypass mode acceleration process being increased by 11.8%-14.1%.展开更多
We proposed an enhanced image binarization method.The proposed solution incorporates Monte-Carlo simulation into the local thresholding method to address the essential issues with respect to complex background,spatial...We proposed an enhanced image binarization method.The proposed solution incorporates Monte-Carlo simulation into the local thresholding method to address the essential issues with respect to complex background,spatially-changed illumination,and uncertainties of block size in traditional method.The proposed method first partitions the image into square blocks that reflect local characteristics of the image.After image partitioning,each block is binarized using Otsu’s thresholding method.To minimize the influence of the block size and the boundary effect,we incorporate Monte-Carlo simulation into the binarization algorithm.Iterative calculation with varying block sizes during Monte-Carlo simulation generates a probability map,which illustrates the probability of each pixel classified as foreground.By setting a probability threshold,and separating foreground and background of the source image,the final binary image can be obtained.The described method has been tested by benchmark tests.Results demonstrate that the proposed method performs well in dealing with the complex background and illumination condition.展开更多
Direct Monte Carlo(MC) simulation is a powerful probabilistic safety assessment method for accounting dynamics of the system.But it is not efficient at simulating rare events.A biasing transition rate method based on ...Direct Monte Carlo(MC) simulation is a powerful probabilistic safety assessment method for accounting dynamics of the system.But it is not efficient at simulating rare events.A biasing transition rate method based on direct MC simulation is proposed to solve the problem in this paper.This method biases transition rates of the components by adding virtual components to them in series to increase the occurrence probability of the rare event,hence the decrease in the variance of MC estimator.Several cases are used to benchmark this method.The results show that the method is effective at modeling system failure and is more efficient at collecting evidence of rare events than the direct MC simulation.The performance is greatly improved by the biasing transition rate method.展开更多
The random distribution problem of dendrite preferred growth direction was settled by random grid method.This method was used to study the influence of forced laminar flow effect on multiple grains during solidificati...The random distribution problem of dendrite preferred growth direction was settled by random grid method.This method was used to study the influence of forced laminar flow effect on multiple grains during solidification.Taking high pure succinonitrile （SCN） undercooled melt as an example,the forced laminar flow effect on multiple grains was studied by phase-field model of single grain which coupled with flow equations at non-isothermal condition.The simulation results show that the random grid method can reasonably settle the problem of random distribution and is more effective.When the solid fraction is relatively low,melt particles flow around the downstream side of dendrite,and the flow velocity between two dendrite arms becomes high.At the stage of solidification time less than 1800Δt,every dendrite grows freely;the upstream dendrites are stronger than the downstream ones.The higher the melt flow rate,the higher the solid fraction.However,when the solid fraction is relatively high,the dendrite arm intertwins and only a little residual melt which is not encapsulated can flow;the solid fraction will gradually tend to equal to solid fraction of melt without flow.展开更多
Serious riverbank erosion,caused by scouring and soil siltation on the bank slope in the lower reaches of the Tarim River,Northwest China urgently requires a solution.Plant roots play an important role in enhancing so...Serious riverbank erosion,caused by scouring and soil siltation on the bank slope in the lower reaches of the Tarim River,Northwest China urgently requires a solution.Plant roots play an important role in enhancing soil shear strength on the slopes to maintain slope soils,but the extent of enhancement of soil shear strength by different root distribution patterns is unclear.The study used a combination of indoor experiments and numerical simulation to investigate the effects of varying plant root morphologies on the shear strength of the sandy soil in the Tarim River.The results showed that:(1)by counting the root morphology of dominant vegetation on the bank slope,we summarized the root morphology of dominant vegetation along the coast as vertical,horizontal,and claw type;(2)the shear strength of root-soil composites(RSCs)was significantly higher than that of remolded soil,and the presence of root system made the strain-softening of soil body significantly weakened so that RSCs had better mechanical properties;and(3)compared with the lateral roots,the average particle contact degree of vertical root system was higher,and the transition zone of shear strength was more prominent.Hence,vegetation with vertical root system had the best effect on soil protection and slope fixation.The results of this study have important guiding significance for prevention and control of soil erosion in the Tarim River basin,the restoration of riparian ecosystems,and the planning of water conservancy projects.展开更多
Directional roof cutting(DRC)is one of the key techniques in non-pillar coal mining with self-formed entries(NCMSE)mining method.Due to the inability to accurately measure the expansion coefficient of the goaf rock ma...Directional roof cutting(DRC)is one of the key techniques in non-pillar coal mining with self-formed entries(NCMSE)mining method.Due to the inability to accurately measure the expansion coefficient of the goaf rock mass,the implementation of this technology often encounters design challenges,leading to suboptimal results and increased costs.This paper establishes a structural analysis model of the goaf working face roof,revealing the failure mechanism of DRC,and clarifies the positive role of DRC in improving the stress of the roadway surrounding rock and reducing the subsidence of the roof through numerical simulation experiments.On this basis,the paper further analyses the roadway pressure and roof settlement under different DRC design heights,and ultimately proposes an optimized design method for the DRC height.The results indicate that the implementation of DRC can significantly optimize the stress environment of the working face roadway surrounding rock.At the same time,during the application of DRC,three scenarios may arise:insufficient,reasonable,and excessive DRC height.Insufficient height will significantly reduce the effectiveness of the technology,while excessive height has little impact on the implementation effect but will greatly increase construction costs and difficulty.Engineering verification shows that the optimized DRC design method proposed in this paper reduces the peak stress of the protective coal pillar in the roadway by 27.2%and the central subsidence of the roof by 41.8%,demonstrating excellent application results.This method provides technical support for the further promotion of NCMSE mining method.展开更多
Fluid-structure-particle interactions in three spatial dimensions happen in many environmental and engineering flows.This paper presents the parallel algorithms for the hybrid diffuse and sharp interface immersed boun...Fluid-structure-particle interactions in three spatial dimensions happen in many environmental and engineering flows.This paper presents the parallel algorithms for the hybrid diffuse and sharp interface immersed boundary(IB)method developed in our previous work.For the moving structure modeled using the sharp interface IB method,a recursive box method is developed for efficiently classifying the background grid nodes.For the particles modeled using the diffuse interface IB method,a‘master-slave’approach is adopted.For the particle-particle interaction(PPI)and particle-structure interaction(PSI),a fast algorithm for classifying the active and inactive Lagrangian points,which discretize the particle surface,is developed for the‘dry’contact approach.The results show that the proposed recursive box method can reduce the classifying time from 52seconds to 0.3 seconds.Acceptable parallel efficiency is obtained for cases with different particle concentrations.Furthermore,the lubrication model is utilized when a particle approaches a wall,enabling an accurate simulation of the rebounding phenomena in the benchmark particle-wall collision problem.At last,the capability of the proposed computational framework is demonstrated by simulating particle-laden turbulent channel flows with rough walls.展开更多
The Monte- Carlo method is used to simulate the surface fatigue crack growth rate for offshore structural steel E36-Z35, and to determine the distributions and relevance of the parameters in the Paris equation. By thi...The Monte- Carlo method is used to simulate the surface fatigue crack growth rate for offshore structural steel E36-Z35, and to determine the distributions and relevance of the parameters in the Paris equation. By this method, the time and cost of fatigue crack propagation testing can be reduced. The application of the method is demonstrated by use of four sets of fatigue crack propagation data for offshore structural steel E36-Z35. A comparison of the test data with the theoretical prediction for surface crack growth rate shows the application of the simulation method to the fatigue crack propagation tests is successful.展开更多
Atomic growth process and structure of Amorphous Electroless Coating have been studied, using Monte-Carlo simulation method. The simulation results of amorphous Ni80P20 coating show that PDFs are in accordance with pr...Atomic growth process and structure of Amorphous Electroless Coating have been studied, using Monte-Carlo simulation method. The simulation results of amorphous Ni80P20 coating show that PDFs are in accordance with practical values. The migrations of adatoms in coating's growth are different from that of solidification of amorphous materials. In some cases, the migrated adatoms in the process of growth of amorphous coating are not enough to occupy all vacancies and traps, so the amorphous coating is micro-porous. The immovable probability k and the largest migration distance of adatoms, which lie on the electroless bath components, affect the PDF, volume density and microporosity remarkably.展开更多
The new improved directional vector simulation method foranalyzing the reliability of struc- tural systems failure probabilityis researched. This paper also points out the defects of generaldirectional vector simulati...The new improved directional vector simulation method foranalyzing the reliability of struc- tural systems failure probabilityis researched. This paper also points out the defects of generaldirectional vector simulation, and gives rise to a new higheraccuracy approximate integral formula of structural systems failureprobability. A new geometric meaning of characteristic function isobtained. A new simple method of generating uniformly distributedrandom vector sample sin n-dimensional unit hyper-spherical surfaceis put forward and strictly proved. This method is easy to put intopractice. Numerical examples are given to show the applicability andeffectiveness of the suggested approach to structural systemsreliability problems.展开更多
The boundary heat flow has important significance for the microstructures of directional solidified binary alloy.Interface evolution of the directional solidified microstructure with different boundary heat flow was d...The boundary heat flow has important significance for the microstructures of directional solidified binary alloy.Interface evolution of the directional solidified microstructure with different boundary heat flow was discussed.In this study, only one interface was allowed to have heat flow, and Neumann boundary conditions were imposed at the other three interfaces.From the calculated results, it was found that different boundary heat flows will result in different microstructures.When the boundary heat flow equals to 20 W·cm-2, the growth of longitudinal side branches is accelerated and the growth of transverse side branches is restrained, and meanwhile, there is dendritic remelting in the calculation domain.When the boundary heat flow equals to 40 W·cm-2, the growths of the transverse and longitudinal side branches compete with each other, and when the boundary heat flow equals to 100-200 W·cm-2, the growth of transverse side branches dominates absolutely.The temperature field of dendritic growth was analyzed and the relation between boundary heat flow and temperature field was also investigated.展开更多
The following article has been retracted due to the investigation of complaints received against it. Mr. Mohammadali Ghorbani (corresponding author and also the last author) cheated the author’s name: Alireza Heidari...The following article has been retracted due to the investigation of complaints received against it. Mr. Mohammadali Ghorbani (corresponding author and also the last author) cheated the author’s name: Alireza Heidari. The scientific community takes a very strong view on this matter and we treat all unethical behavior such as plagiarism seriously. This paper published in Vol.3 No.3, 260-265, 2012, has been removed from this site.展开更多
For solving higher dimensional diffusion equations with an inhomogeneous diffusion coefficient,Monte Carlo（MC） techniques are considered to be more effective than other algorithms, such as finite element method or f...For solving higher dimensional diffusion equations with an inhomogeneous diffusion coefficient,Monte Carlo（MC） techniques are considered to be more effective than other algorithms, such as finite element method or finite difference method. The inhomogeneity of diffusion coefficient strongly limits the use of different numerical techniques. For better convergence, methods with higher orders have been kept forward to allow MC codes with large step size. The main focus of this work is to look for operators that can produce converging results for large step sizes. As a first step, our comparative analysis has been applied to a general stochastic problem.Subsequently, our formulization is applied to the problem of pitch angle scattering resulting from Coulomb collisions of charge particles in the toroidal devices.展开更多
The investigation was carried out on the technical problems of finishing the inner surface of elbow parts and the action mechanism of particles in elbow precision machining by abrasive flow.This work was analyzed and ...The investigation was carried out on the technical problems of finishing the inner surface of elbow parts and the action mechanism of particles in elbow precision machining by abrasive flow.This work was analyzed and researched by combining theory,numerical and experimental methods.The direct simulation Monte Carlo(DSMC)method and the finite element analysis method were combined to reveal the random collision of particles during the precision machining of abrasive flow.Under different inlet velocity,volume fraction and abrasive particle size,the dynamic pressure and turbulence flow energy of abrasive flow in elbow were analyzed,and the machining mechanism of particles on the wall and the influence of different machining parameters on the precision machining quality of abrasive flow were obtained.The test results show the order of the influence of different parameters on the quality of abrasive flow precision machining and establish the optimal process parameters.The results of the surface morphology before and after the precision machining of the inner surface of the elbow are discussed,and the surface roughness Ra value is reduced from 1.125μm to 0.295μm after the precision machining of the abrasive flow.The application of DSMC method provides special insights for the development of abrasive flow technology.展开更多
A new hybrid method, Monte-Carlo-Heat-Flux （MCHF） method, was presented to analyze the radiative heat transfer of participating medium in a three-dimensional rectangular enclosure using combined the Monte-Carlo meth...A new hybrid method, Monte-Carlo-Heat-Flux （MCHF） method, was presented to analyze the radiative heat transfer of participating medium in a three-dimensional rectangular enclosure using combined the Monte-Carlo method with the heat flux method. Its accuracy and reliability was proved by comparing the computational results with exact results from classical ＂Zone Method＂.展开更多
In this paper, dynamic simulation of a beam-like structure with a transverse open crack subjected to a random moving mass oscillator is investigated. The simultaneous effect of a crack and a random oscillator has not ...In this paper, dynamic simulation of a beam-like structure with a transverse open crack subjected to a random moving mass oscillator is investigated. The simultaneous effect of a crack and a random oscillator has not been addressed up to now. The crack in the beam at different locations and with different depths is considered as one group of damage, each as an individual imperfection. In addition, bearing immobility is considered as another type of problem in the beam. Mass, stiffness, damping and velocity of the oscillator are assumed to be random parameters. An improved perturbation technique is applied to reduce the simulation time. It was found that there is a maximum value of the variance of each uncertain parameter, in which the maximum reliability of the perturbation method can be achieved, and that this maximum value can be obtained by the Alpha-Hilber Monte-Carlo simulation method. The simulation results reveal that the mass and the velocity uncertainty cause high uncertainty in the deflection of the beam. Also, the pattern of the deflection is not affected by different random oscillator parameters, and as a result, the type of damage can be identified even with high uncertainty. Moreover, the deflection in the nodes around the mid-span of the beam provides the best information regarding the imperfections, and consequently leads to the best sensor locations in an actual experiment.展开更多
For the sake of a more accurate shell boundary and calculation of radiation heat transfer in the Directional Solidification(DS) process, a radiation heat transfer model based on the Finite Element Method(FEM)is develo...For the sake of a more accurate shell boundary and calculation of radiation heat transfer in the Directional Solidification(DS) process, a radiation heat transfer model based on the Finite Element Method(FEM)is developed in this study. Key technologies, such as distinguishing boundaries automatically, local matrix and lumped heat capacity matrix, are also stated. In order to analyze the effect of withdrawing rate on DS process,the solidification processes of a complex superalloy turbine blade in the High Rate Solidification(HRS) process with different withdrawing rates are simulated; and by comparing the simulation results, it is found that the most suitable withdrawing rate is determined to be 5.0 mm·min^(-1). Finally, the accuracy and reliability of the radiation heat transfer model are verified, because of the accordance of simulation results with practical process.展开更多
A three-dimensional finite element thermal model in direct metal laser sintering(DMLS) including the effect of powder-to-solid transition were established to predict sintering zone, which benefited the determination o...A three-dimensional finite element thermal model in direct metal laser sintering(DMLS) including the effect of powder-to-solid transition were established to predict sintering zone, which benefited the determination of suitable process parameters in DMLS. The nonlinear transient model of the metals thermal conductivity for powder-to-solid transition was developed. The model uses solid thermal properties of material in liquid-phase zone, transitional ones in sintering or sintered zone and powder ones in unsintered zones of powder bed to predict, respectively. Sintering zone boundary was estimated by maximum temperature history profile. Experiments were carried out using multi-component Cu-based metal powder. Compared experimental and predicted results, the mean error of sintering depth and width are 7.8% and 14.4%, respectively, which confirms the accuracy of the FEM prediction.展开更多
A modified direct optimization method is proposed to solve the optimal multi-revolution transfer with low-thrust between Earth-orbits. First, through parameterizing the control steering angles by costate variables, th...A modified direct optimization method is proposed to solve the optimal multi-revolution transfer with low-thrust between Earth-orbits. First, through parameterizing the control steering angles by costate variables, the search space of free parameters has been decreased. Then, in order to obtain the global optimal solution effectively and robustly, the simulated annealing and penalty function strategies were used to handle the constraints, and a GA/SQP hybrid optimization algorithm was utilized to solve the parameter optimization problem, in which, a feasible suboptimal solution obtained by GA was submitted as an initial parameter set to SQP for refinement. Comparing to the classical direct method, this novel method has fewer free parameters, needs not initial guesses, and has higher computation precision. An optimal-fuel transfer problem from LEO to GEO was taken as an example to validate the proposed approach. The results of simulation indicate that our approach is available to solve the problem of optimal muhi-revolution transfer between Earth-orbits.展开更多
Many projects have recently been carried out and proposed for observing high energy electrons since it is realized that cosmic ray electrons are very important when studying the dark matter particles and the accelerat...Many projects have recently been carried out and proposed for observing high energy electrons since it is realized that cosmic ray electrons are very important when studying the dark matter particles and the acceleration mechanism of cosmic rays. An imaging calorimeter, BETS （Balloon-borne Electron Telescope with Scintillator fiber）, has been developed for this purpose. Using pattern analysis of the shower development, the electrons can be selected from those primary cosmic ray proton events with flux heights one-tenth that of the electrons. The Monte-Carlo simulation is indispensable for the instrument design, the signal trigger and the data analysis. We present different shower simulation codes and compare the simulation results with the beam test and the flight data of BETS. We conclude that the code FLUKA2002 gives the most consistent results with the experimental data.展开更多
基金co-supported by the National Science and Technology Major Project,China(No.J2019-I-0015-0014)the National Natural Science Foundation of China(No.52372397).
文摘To design the optimum acceleration control schedule for the Adaptive Cycle Engine(ACE)in the full flight envelope,this paper establishes a direct simulation model of the ACE transient state.In this model,geometric parameters are used to replace the component state parameters.The corresponding relationship between geometric parameters and component state parameters is determined by sensitivity analysis.The geometric variables are controlled when the geometric adjustment speed exceeds the limit,and at the same time the corresponding component state parameters are iterated.The gradient optimization algorism is used to optimize the ground acceleration process of ACE,and the control schedule in terms of operating point of compression components and corrected acceleration rate is used as the full-envelope acceleration control schedule based on the similarity principle.The acceleration control schedules of the triple-bypass mode and the double-bypass mode are designed in this paper.The acceleration processes under various flight conditions are simulated using the acceleration control schedules.Compared with the acceleration process with the linear geometric adjustment schedule,the acceleration performance of ACE is improved by the acceleration control schedule,with the impulse of the acceleration process of the triple-bypass mode being increased by 8.7%-12.3% and the impulse of the double-bypass mode acceleration process being increased by 11.8%-14.1%.
基金Project(2018YFC1505401)supported by the National Key R&D Program of ChinaProject(41702310)supported by the National Natural Science Foundation of China+1 种基金Project(SKLGP2017K014)supported by the Foundation of State Key Laboratory of Geohazard Prevention and Geo-environment Protection,ChinaProject(2018JJ3644)supported by the Natural Science Foundation of Hunan Province,China
文摘We proposed an enhanced image binarization method.The proposed solution incorporates Monte-Carlo simulation into the local thresholding method to address the essential issues with respect to complex background,spatially-changed illumination,and uncertainties of block size in traditional method.The proposed method first partitions the image into square blocks that reflect local characteristics of the image.After image partitioning,each block is binarized using Otsu’s thresholding method.To minimize the influence of the block size and the boundary effect,we incorporate Monte-Carlo simulation into the binarization algorithm.Iterative calculation with varying block sizes during Monte-Carlo simulation generates a probability map,which illustrates the probability of each pixel classified as foreground.By setting a probability threshold,and separating foreground and background of the source image,the final binary image can be obtained.The described method has been tested by benchmark tests.Results demonstrate that the proposed method performs well in dealing with the complex background and illumination condition.
基金supported by the Special Projects of International Thermonuclear Experimental Reactor(2015GB116000)the Strategic Priority Research Program of Chinese Academy of Sciences(No.XDA03040000)+1 种基金the Informatizational Special Projects of Chinese Academy of Sciences(No.XXH12504-1-09)the Major/Innovative Program of Development Foundation of Hefei Center for Physical Science and Technology(No.2014FXCX004)
文摘Direct Monte Carlo(MC) simulation is a powerful probabilistic safety assessment method for accounting dynamics of the system.But it is not efficient at simulating rare events.A biasing transition rate method based on direct MC simulation is proposed to solve the problem in this paper.This method biases transition rates of the components by adding virtual components to them in series to increase the occurrence probability of the rare event,hence the decrease in the variance of MC estimator.Several cases are used to benchmark this method.The results show that the method is effective at modeling system failure and is more efficient at collecting evidence of rare events than the direct MC simulation.The performance is greatly improved by the biasing transition rate method.
基金Project(10964004) supported by the National Natural Science Foundation of ChinaProject(20070731001) supported by Research Fund for the Doctoral Program of China+1 种基金 Project(096RJZA104) supported by the Natural Science Foundation of Gansu Province,ChinaProject(SB14200801) supported by the Doctoral Fund of Lanzhou University of Technology,China
文摘The random distribution problem of dendrite preferred growth direction was settled by random grid method.This method was used to study the influence of forced laminar flow effect on multiple grains during solidification.Taking high pure succinonitrile （SCN） undercooled melt as an example,the forced laminar flow effect on multiple grains was studied by phase-field model of single grain which coupled with flow equations at non-isothermal condition.The simulation results show that the random grid method can reasonably settle the problem of random distribution and is more effective.When the solid fraction is relatively low,melt particles flow around the downstream side of dendrite,and the flow velocity between two dendrite arms becomes high.At the stage of solidification time less than 1800Δt,every dendrite grows freely;the upstream dendrites are stronger than the downstream ones.The higher the melt flow rate,the higher the solid fraction.However,when the solid fraction is relatively high,the dendrite arm intertwins and only a little residual melt which is not encapsulated can flow;the solid fraction will gradually tend to equal to solid fraction of melt without flow.
基金funded by the Key Research and Development Task of Xinjiang Uygur Autonomous Region, China (2022B03024-3)
文摘Serious riverbank erosion,caused by scouring and soil siltation on the bank slope in the lower reaches of the Tarim River,Northwest China urgently requires a solution.Plant roots play an important role in enhancing soil shear strength on the slopes to maintain slope soils,but the extent of enhancement of soil shear strength by different root distribution patterns is unclear.The study used a combination of indoor experiments and numerical simulation to investigate the effects of varying plant root morphologies on the shear strength of the sandy soil in the Tarim River.The results showed that:(1)by counting the root morphology of dominant vegetation on the bank slope,we summarized the root morphology of dominant vegetation along the coast as vertical,horizontal,and claw type;(2)the shear strength of root-soil composites(RSCs)was significantly higher than that of remolded soil,and the presence of root system made the strain-softening of soil body significantly weakened so that RSCs had better mechanical properties;and(3)compared with the lateral roots,the average particle contact degree of vertical root system was higher,and the transition zone of shear strength was more prominent.Hence,vegetation with vertical root system had the best effect on soil protection and slope fixation.The results of this study have important guiding significance for prevention and control of soil erosion in the Tarim River basin,the restoration of riparian ecosystems,and the planning of water conservancy projects.
基金funded by the National Natural Science Foundation of China(52074298)Beijing Municipal Natural Science Foundation(8232056)+1 种基金Guizhou Province science and technology plan project([2020]3008)Liulin Energy and Environment Academician Workstation(2022XDHZ12).
文摘Directional roof cutting(DRC)is one of the key techniques in non-pillar coal mining with self-formed entries(NCMSE)mining method.Due to the inability to accurately measure the expansion coefficient of the goaf rock mass,the implementation of this technology often encounters design challenges,leading to suboptimal results and increased costs.This paper establishes a structural analysis model of the goaf working face roof,revealing the failure mechanism of DRC,and clarifies the positive role of DRC in improving the stress of the roadway surrounding rock and reducing the subsidence of the roof through numerical simulation experiments.On this basis,the paper further analyses the roadway pressure and roof settlement under different DRC design heights,and ultimately proposes an optimized design method for the DRC height.The results indicate that the implementation of DRC can significantly optimize the stress environment of the working face roadway surrounding rock.At the same time,during the application of DRC,three scenarios may arise:insufficient,reasonable,and excessive DRC height.Insufficient height will significantly reduce the effectiveness of the technology,while excessive height has little impact on the implementation effect but will greatly increase construction costs and difficulty.Engineering verification shows that the optimized DRC design method proposed in this paper reduces the peak stress of the protective coal pillar in the roadway by 27.2%and the central subsidence of the roof by 41.8%,demonstrating excellent application results.This method provides technical support for the further promotion of NCMSE mining method.
基金Project supported by the National Natural Science Foundation of China(Nos.12202456 and12172360)the Basic Science Center Program for“Multiscale Problems in Nonlinear Mechanics”of the National Natural Science Foundation of China(No.11988102)the China Postdoctoral Science Foundation(No.2021M693241)。
文摘Fluid-structure-particle interactions in three spatial dimensions happen in many environmental and engineering flows.This paper presents the parallel algorithms for the hybrid diffuse and sharp interface immersed boundary(IB)method developed in our previous work.For the moving structure modeled using the sharp interface IB method,a recursive box method is developed for efficiently classifying the background grid nodes.For the particles modeled using the diffuse interface IB method,a‘master-slave’approach is adopted.For the particle-particle interaction(PPI)and particle-structure interaction(PSI),a fast algorithm for classifying the active and inactive Lagrangian points,which discretize the particle surface,is developed for the‘dry’contact approach.The results show that the proposed recursive box method can reduce the classifying time from 52seconds to 0.3 seconds.Acceptable parallel efficiency is obtained for cases with different particle concentrations.Furthermore,the lubrication model is utilized when a particle approaches a wall,enabling an accurate simulation of the rebounding phenomena in the benchmark particle-wall collision problem.At last,the capability of the proposed computational framework is demonstrated by simulating particle-laden turbulent channel flows with rough walls.
文摘The Monte- Carlo method is used to simulate the surface fatigue crack growth rate for offshore structural steel E36-Z35, and to determine the distributions and relevance of the parameters in the Paris equation. By this method, the time and cost of fatigue crack propagation testing can be reduced. The application of the method is demonstrated by use of four sets of fatigue crack propagation data for offshore structural steel E36-Z35. A comparison of the test data with the theoretical prediction for surface crack growth rate shows the application of the simulation method to the fatigue crack propagation tests is successful.
文摘Atomic growth process and structure of Amorphous Electroless Coating have been studied, using Monte-Carlo simulation method. The simulation results of amorphous Ni80P20 coating show that PDFs are in accordance with practical values. The migrations of adatoms in coating's growth are different from that of solidification of amorphous materials. In some cases, the migrated adatoms in the process of growth of amorphous coating are not enough to occupy all vacancies and traps, so the amorphous coating is micro-porous. The immovable probability k and the largest migration distance of adatoms, which lie on the electroless bath components, affect the PDF, volume density and microporosity remarkably.
文摘The new improved directional vector simulation method foranalyzing the reliability of struc- tural systems failure probabilityis researched. This paper also points out the defects of generaldirectional vector simulation, and gives rise to a new higheraccuracy approximate integral formula of structural systems failureprobability. A new geometric meaning of characteristic function isobtained. A new simple method of generating uniformly distributedrandom vector sample sin n-dimensional unit hyper-spherical surfaceis put forward and strictly proved. This method is easy to put intopractice. Numerical examples are given to show the applicability andeffectiveness of the suggested approach to structural systemsreliability problems.
文摘The boundary heat flow has important significance for the microstructures of directional solidified binary alloy.Interface evolution of the directional solidified microstructure with different boundary heat flow was discussed.In this study, only one interface was allowed to have heat flow, and Neumann boundary conditions were imposed at the other three interfaces.From the calculated results, it was found that different boundary heat flows will result in different microstructures.When the boundary heat flow equals to 20 W·cm-2, the growth of longitudinal side branches is accelerated and the growth of transverse side branches is restrained, and meanwhile, there is dendritic remelting in the calculation domain.When the boundary heat flow equals to 40 W·cm-2, the growths of the transverse and longitudinal side branches compete with each other, and when the boundary heat flow equals to 100-200 W·cm-2, the growth of transverse side branches dominates absolutely.The temperature field of dendritic growth was analyzed and the relation between boundary heat flow and temperature field was also investigated.
文摘The following article has been retracted due to the investigation of complaints received against it. Mr. Mohammadali Ghorbani (corresponding author and also the last author) cheated the author’s name: Alireza Heidari. The scientific community takes a very strong view on this matter and we treat all unethical behavior such as plagiarism seriously. This paper published in Vol.3 No.3, 260-265, 2012, has been removed from this site.
基金supported in part by the Higher Education Commission of Pakistan under PPCR programsupported by the National Magnetic Confinement Fusion Program under Grant No.2013GB104004Fundamental Research Fund for Chinese Central Universities
文摘For solving higher dimensional diffusion equations with an inhomogeneous diffusion coefficient,Monte Carlo（MC） techniques are considered to be more effective than other algorithms, such as finite element method or finite difference method. The inhomogeneity of diffusion coefficient strongly limits the use of different numerical techniques. For better convergence, methods with higher orders have been kept forward to allow MC codes with large step size. The main focus of this work is to look for operators that can produce converging results for large step sizes. As a first step, our comparative analysis has been applied to a general stochastic problem.Subsequently, our formulization is applied to the problem of pitch angle scattering resulting from Coulomb collisions of charge particles in the toroidal devices.
基金Projects(51206011,U1937201)supported by the National Natural Science Foundation of ChinaProject(20200301040RQ)supported by the Science and Technology Development Program of Jilin Province,China+1 种基金Project(JJKH20190541KJ)supported by the Education Department of Jilin Province,ChinaProject(18DY017)supported by Changchun Science and Technology Program of Changchun City,China。
文摘The investigation was carried out on the technical problems of finishing the inner surface of elbow parts and the action mechanism of particles in elbow precision machining by abrasive flow.This work was analyzed and researched by combining theory,numerical and experimental methods.The direct simulation Monte Carlo(DSMC)method and the finite element analysis method were combined to reveal the random collision of particles during the precision machining of abrasive flow.Under different inlet velocity,volume fraction and abrasive particle size,the dynamic pressure and turbulence flow energy of abrasive flow in elbow were analyzed,and the machining mechanism of particles on the wall and the influence of different machining parameters on the precision machining quality of abrasive flow were obtained.The test results show the order of the influence of different parameters on the quality of abrasive flow precision machining and establish the optimal process parameters.The results of the surface morphology before and after the precision machining of the inner surface of the elbow are discussed,and the surface roughness Ra value is reduced from 1.125μm to 0.295μm after the precision machining of the abrasive flow.The application of DSMC method provides special insights for the development of abrasive flow technology.
基金financially supported by the National Natural Science Foundation of China (No.50464004)
文摘A new hybrid method, Monte-Carlo-Heat-Flux （MCHF） method, was presented to analyze the radiative heat transfer of participating medium in a three-dimensional rectangular enclosure using combined the Monte-Carlo method with the heat flux method. Its accuracy and reliability was proved by comparing the computational results with exact results from classical ＂Zone Method＂.
文摘In this paper, dynamic simulation of a beam-like structure with a transverse open crack subjected to a random moving mass oscillator is investigated. The simultaneous effect of a crack and a random oscillator has not been addressed up to now. The crack in the beam at different locations and with different depths is considered as one group of damage, each as an individual imperfection. In addition, bearing immobility is considered as another type of problem in the beam. Mass, stiffness, damping and velocity of the oscillator are assumed to be random parameters. An improved perturbation technique is applied to reduce the simulation time. It was found that there is a maximum value of the variance of each uncertain parameter, in which the maximum reliability of the perturbation method can be achieved, and that this maximum value can be obtained by the Alpha-Hilber Monte-Carlo simulation method. The simulation results reveal that the mass and the velocity uncertainty cause high uncertainty in the deflection of the beam. Also, the pattern of the deflection is not affected by different random oscillator parameters, and as a result, the type of damage can be identified even with high uncertainty. Moreover, the deflection in the nodes around the mid-span of the beam provides the best information regarding the imperfections, and consequently leads to the best sensor locations in an actual experiment.
基金financially supported by the Program for New Century Excellent Talents in University(No.NCET-13-0229,NCET-09-0396)the National Science & Technology Key Projects of Numerical Control(No.2012ZX04010-031,2012ZX0412-011)the National High Technology Research and Development Program("863"Program)of China(No.2013031003)
文摘For the sake of a more accurate shell boundary and calculation of radiation heat transfer in the Directional Solidification(DS) process, a radiation heat transfer model based on the Finite Element Method(FEM)is developed in this study. Key technologies, such as distinguishing boundaries automatically, local matrix and lumped heat capacity matrix, are also stated. In order to analyze the effect of withdrawing rate on DS process,the solidification processes of a complex superalloy turbine blade in the High Rate Solidification(HRS) process with different withdrawing rates are simulated; and by comparing the simulation results, it is found that the most suitable withdrawing rate is determined to be 5.0 mm·min^(-1). Finally, the accuracy and reliability of the radiation heat transfer model are verified, because of the accordance of simulation results with practical process.
文摘A three-dimensional finite element thermal model in direct metal laser sintering(DMLS) including the effect of powder-to-solid transition were established to predict sintering zone, which benefited the determination of suitable process parameters in DMLS. The nonlinear transient model of the metals thermal conductivity for powder-to-solid transition was developed. The model uses solid thermal properties of material in liquid-phase zone, transitional ones in sintering or sintered zone and powder ones in unsintered zones of powder bed to predict, respectively. Sintering zone boundary was estimated by maximum temperature history profile. Experiments were carried out using multi-component Cu-based metal powder. Compared experimental and predicted results, the mean error of sintering depth and width are 7.8% and 14.4%, respectively, which confirms the accuracy of the FEM prediction.
基金Sponsored by the National Natural Science Foundation of China(Grant No.10672044)
文摘A modified direct optimization method is proposed to solve the optimal multi-revolution transfer with low-thrust between Earth-orbits. First, through parameterizing the control steering angles by costate variables, the search space of free parameters has been decreased. Then, in order to obtain the global optimal solution effectively and robustly, the simulated annealing and penalty function strategies were used to handle the constraints, and a GA/SQP hybrid optimization algorithm was utilized to solve the parameter optimization problem, in which, a feasible suboptimal solution obtained by GA was submitted as an initial parameter set to SQP for refinement. Comparing to the classical direct method, this novel method has fewer free parameters, needs not initial guesses, and has higher computation precision. An optimal-fuel transfer problem from LEO to GEO was taken as an example to validate the proposed approach. The results of simulation indicate that our approach is available to solve the problem of optimal muhi-revolution transfer between Earth-orbits.
基金Supported by the National Natural Science Foundation of China.
文摘Many projects have recently been carried out and proposed for observing high energy electrons since it is realized that cosmic ray electrons are very important when studying the dark matter particles and the acceleration mechanism of cosmic rays. An imaging calorimeter, BETS （Balloon-borne Electron Telescope with Scintillator fiber）, has been developed for this purpose. Using pattern analysis of the shower development, the electrons can be selected from those primary cosmic ray proton events with flux heights one-tenth that of the electrons. The Monte-Carlo simulation is indispensable for the instrument design, the signal trigger and the data analysis. We present different shower simulation codes and compare the simulation results with the beam test and the flight data of BETS. We conclude that the code FLUKA2002 gives the most consistent results with the experimental data.
