Simulated Method on the Release of Sediment Pollutants in the Temple Lake

[ Objective] The research aimed to study simulated method on the release of sediment pollutants in the Temple Lake. [ Method] For existing pollution status of the Temple Lake, sediments and overlying water were sampled by the integrated device. Based on the technical route, through the orthogonal test, release processes of nitrogen and phosphorus from sediment were simulated, and their release speeds were estimated under different temperatures, pH and DO in the laboratory. Relationship between environmental factor and pollutant release was explored. [ Result] The integrated device of sediment sampling and test had simple structure, cheap cost and easy operation, and sediment pollution test could be conducted without disturbing sediment after collecting sediment. Among the studied three impact factors, temperature had the maximum influence on contaminant release process, which should be considered as the important factor of future research. [ Conclusion] Compared with routine research method of sediment contaminant release from other lakes, simulated method of contaminant release from sediment in the Temple Lake had many advantaqes, which was favorable for scientific research in latter period and actual lake control.

A spatiotemporal evolution model of a short-circuit arc to a secondary arc based on the improved charge simulation method

The initial shape of the secondary arc considerably influences its subsequent shape.To establish the model for the arcing time of the secondary arc and modify the single-phase reclosing sequence,theoretical and experimental analysis of the evolution process of the short-circuit arc to the secondary arc is critical.In this study,an improved charge simulation method was used to develop the internal-space electric-field model of the short-circuit arc.The intensity of the electric field was used as an independent variable to describe the initial shape of the secondary arc.A secondary arc evolution model was developed based on this model.Moreover,the accuracy of the model was evaluated by comparison with physical experimental results.When the secondary arc current increased,the arcing time and dispersion increased.There is an overall trend of increasing arc length with increasing arcing time.Nevertheless,there is a reduction in arc length during arc ignition due to short circuits between the arc columns.Furthermore,the arcing time decreased in the range of 0°-90°as the angle between the wind direction and the x-axis increased.This work investigated the method by which short-circuit arcs evolve into secondary arcs.The results can be used to develop the secondary arc evolution model and to provide both a technical and theoretical basis for secondary arc suppression.

Multi-group ant colony algorithm based on simulated annealing method

To overcome the default of single search tendency, the ants in the colony are divided into several sub-groups. The ants in different subgroups have different trail information and expectation coefficients. The simulated annealing method is introduced to the algorithm. Through setting the temperature changing with the iterations, after each turn of tours, the solution set obtained by the ants is taken as the candidate set. The update set is obtained by adding the solutions in the candidate set to the previous update set with the probability determined by the temperature. The solutions in the candidate set are used to update the trail information. In each turn of updating, the current best solution is also used to enhance the trail information on the current best route. The trail information is reset when the algorithm is in stagnation state. The computer experiments demonstrate that the proposed algorithm has higher stability and convergence speed.

水下航行体湍流数值模拟方法研究进展综述

In this paper,we present an overview of numerical simulation methods for the flow around typical underwater vehicles at high Reynolds numbers,which highlights the dominant flow structures in different regions of interest.This overview covers the forebody,midbody,stern,wake region,and appendages and summarizes flow phenomena,including laminar-to-turbulent transition,turbulent boundary layers,flow under the influence of curvatures,wake interactions,and all associated complex vortex structures.Furthermore,the current issues and challenges of capturing these flow structures are addressed.This overview provides a deep insight into the use of numerical simulation methods,including the Reynolds-averaged Navier–Stokes(RANS)method,large eddy simulation(LES)method,and the hybrid RANS/LES method,and evaluates their applicability in capturing detailed flow features.

Assessing cutter-rock interaction during TBM tunnelling in granite:Large-scale standing rotary cutting tests and 3D DEM simulations

The widespread utilisation of tunnel boring machines(TBMs)in underground construction engineering requires a detailed investigation of the cutter-rock interaction.In this paper,we conduct a series of largescale standing rotary cutting tests on granite in conjunction with high-fidelity numerical simulations based on a particle-type discrete element method(DEM)to explore the effects of key cutting parameters on the TBM cutter performance and the distribution of cutter-rock contact stresses.The assessment results of cutter performance obtained from the cutting tests and numerical simulations reveal similar dependencies on the key cutting parameters.More specifically,the normal and rolling forces exhibit a positive correlation with penetration but are slightly influenced by the cutting radius.In contrast,the side force decreases as the cutting radius increases.Additionally,the side force shows a positive relationship with the penetration for smaller cutting radii but tends to become negative as the cutting radius increases.The cutter's relative effectiveness in rock breaking is significantly impacted by the penetration but shows little dependency on the cutting radius.Consequently,an optimal penetration is identified,leading to a low boreability index and specific energy.A combined Hertz-Weibull function is developed to fit the cutter-rock contact stress distribution obtained in DEM simulations,whereby an improved CSM(Colorado School of Mines)model is proposed by replacing the original monotonic cutting force distribution with this combined Hertz-Weibull model.The proposed model outperforms the original CSM model as demonstrated by a comparison of the estimated cutting forces with those from the tests/simulations.The findings from this work that advance our understanding of TBM cutter performance have important implications for improving the efficiency and reliability of TBM tunnelling in granite.

On the “Onion Husk” Algorithm for Approximate Solution of the Traveling Salesman Problem

The paper describes some implementation aspects of an algorithm for approximate solution of the traveling salesman problem based on the construction of convex closed contours on the initial set of points (“cities”) and their subsequent combination into a closed path (the so-called contour algorithm or “onion husk” algorithm). A number of heuristics related to the different stages of the algorithm are considered, and various variants of the algorithm based on these heuristics are analyzed. Sets of randomly generated points of different sizes (from 4 to 90 and from 500 to 10,000) were used to test the algorithms. The numerical results obtained are compared with the results of two well-known combinatorial optimization algorithms, namely the algorithm based on the branch and bound method and the simulated annealing algorithm. .

Reconstruction of 3-D digital cores using a hybrid method

A 3-D digital core describes the pore space microstructure of rocks. An X-ray micro CT scan is the most accurate and direct but costly method to obtain a 3-D digital core. In this study, we propose a hybrid method which combines sedimentation simulation and simulated annealing （SA） method to generate 3-D digital cores based on 2-D images of rocks. The method starts with the sedimentation simulation to build a 3-D digital core, which is the initial configuration for the SA method. We update the initial digital core using the SA method to match the auto-correlation function of the 2-D rock image and eventually build the final 3-D digital core. Compared with the typical SA method, the hybrid method has significantly reduced the computation time. Local porosity theory is applied to quantitatively compare the reconstructed 3-D digital cores with the X-ray micro CT 3-D images. The results indicate that the 3-D digital cores reconstructed by the hybrid method have homogeneity and geometric connectivity similar to those of the X-ray micro CT image. The formation factors and permeabilities of the reconstructed 3-D digital cores are estimated using the finite element method （FEM） and lattice Boltzmann method （LBM）, respectively. The simulated results are in good agreement with the experimental measurements. Comparison of the simulation results suggests that the digital cores reconstructed by the hybrid method more closely reflect the true transport properties than the typical SA method alone.

Effects of loading waveforms on rock damage using particle simulation method

The particle simulation method is used to study the effects of loading waveforms （i.e. square, sinusoidal and triangle waveforms） on rock damage at mesoscopic scale. Then some influencing factors on rock damage at the mesoscopic scale, such as loading frequency, stress amplitude, mean stress, confining pressure and loading sequence, are also investigated with sinusoidal waveform in detail. The related numerical results have demonstrated that： 1） the loading waveform has a certain effect on rock failure processes. The square waveform has the most damage within these waveforms, while the triangle waveform has less damage than sinusoidal waveform. In each cycle, the number of microscopic cracks increases in the loading stage, while it keeps nearly constant in the unloading stage. 2） The loading frequency, stress amplitude, mean stress, confining pressure and loading sequence have considerable effects on rock damage subjected to cyclic loading. The higher the loading frequency, stress amplitude and mean stress, the greater the damage the rock accumulated; in contrast, the lower the confining pressure, the greater the damage the rock has accumulated. 3） There is a threshold value of mean stress and stress amplitude, below which no further damage accumulated after the first few cycle loadings. 4） The high-to-low loading sequence has more damage than the low-to-high loading sequence, suggesting that the rock damage is loading-path dependent.

A grouting simulation method for quick-setting slurry in karst conduit:The sequential flow and solidification method

It is difficult to temporally and spatially track and characterize the slurry viscosity in flowing water during grouting simulation.In this study,a sequential flow and solidification(SFS)method considering the spatial-temporal evolution of slurry viscosity in flowing water in karst conduit is proposed.First,a time-dependent model for the threshold function of slurry viscosity is established.During the grouting process,the spatial-temporal evolution of slurry viscosity is revealed by tracking the diffusion behavior of the slurry injected at different times.This method is capable of describing the gradual solidification process of the slurry during grouting.Furthermore,a physical model of grouting in a karst conduit is developed.Second,the effectiveness of the SFS method in grouting simulation is verified by the experiment of grouting conduit in flowing water.The SFS method enables real-time monitoring of fluid velocity and pressure during grouting in flowing water and provides a feasible calculation method for revealing the grouting plugging mechanism in complex karst conduits at different engineering scales.In addition,it can be used to guide the design of grouting tests in flowing water,improve cost efficiency,and provide theoretical basis for optimizing grouting design and slurry selection.

Monte Carlo Analytic Hierarchy Process (MAHP) approach to selection of optimum mining method

One of the most critical and complicated steps in mine design is a selection of suitable mining method based upon geological,geotechnical,geographical,safety and economical parameters.The aim of this study is developing a Monte Carlo simulation to selection the optimum mining method by using effective and major criteria and at the same time,taking subjective judgments of decision makers into consideration.Proposed approach is based on the combination of Monte Carlo simulation with conventional Analytic Hierarchy Process(AHP).Monte Carlo simulation is used to determine the confdence level of each alternative’s score,is calculated by AHP,with the respect to the variance of decision makers’opinion.The proposed method is applied for Jajarm Bauxite Mine in Iran and eventually the most appropriate mining methods for this mine are ranked.

Risk calculation method for complex engineering system

This paper presents a rapid and simple risk calculation method for large and complex engineering systems, the simulated maximum entropy method （SMEM）, which is based on integration of the advantages of the Monte Carlo and maximum entropy methods, thus avoiding the shortcoming of the slow convergence rate of the Monte Carlo method in risk calculation. Application of SMEM in the calculation of reservoir flood discharge risk shows that this method can make full use of the known information under the same conditions and obtain the corresponding probability distribution and the risk value. It not only greatly improves the speed, compared with the Monte Carlo method, but also provides a new approach for the risk calculation in large and complex engineering systems.

APPLICATION OF THE PROBABILISTIC FRACTURE MECHANICS METHOD OF PREDICTING THE FATIGUE LIFE OF TUBULAR JOINTS

In this paper .fatigue strength of T-type tubularjoints subjected to in-planebending or out-of-plane bending load is investigated By, considering material constantsand initial crack sizes as random variables and applying Monte Carlo simulation method, we have giyen a statistical analvsis on fatigue life. Simultaneously, linearregression analyses of computed resuldts are performed and compared with that of theknown experimental data.

Elbow precision machining technology by abrasive flow based on direct Monte Carlo method

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.

Passive simulation method of turbine flow sensors based on the 6-DOF model

A passive simulation method based on the six degrees of freedom(6-DOF)model and dynamic mesh is proposed according to the working principle to study the dynamic characteristics of the turbine flow sensors.This simulation method controls the six degrees of freedom of the impeller using the user-defined functions(UDF)program so that it can only rotate under the impact of fluid.The impeller speed can be calculated in real-time,and the inlet speed can be set with time to obtain the dynamic performance of the turbine flow sensors.Based on this simulation method,three turbine flow sensors with different diameters were simulated,and the reliability of the simulation method was verified by both steady-state and unsteady-state experiments.The results show that the trend of meter factor with flow rate acquired from the simulation is close to the experimental results.The deviation between the simulation and experiment results is low,with a maximum deviation of 2.88%.In the unsteady simulation study,the impeller speed changed with the inlet velocity of the turbine flow sensor,showing good tracking performance.The passive simulation method can be used to predict the dynamic performance of the turbine flow sensor.

Ranking important nodes in complex networks by simulated annealing

In this paper, based on simulated annealing a new method to rank important nodes in complex networks is presented.First, the concept of an importance sequence（IS） to describe the relative importance of nodes in complex networks is defined. Then, a measure used to evaluate the reasonability of an IS is designed. By comparing an IS and the measure of its reasonability to a state of complex networks and the energy of the state, respectively, the method finds the ground state of complex networks by simulated annealing. In other words, the method can construct a most reasonable IS. The results of experiments on real and artificial networks show that this ranking method not only is effective but also can be applied to different kinds of complex networks.

Radiation heat transfer model for complex superalloy turbine blade in directional solidification process based on finite element method

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.

Force Control Compensation Method with Variable Load Stiffness and Damping of the Hydraulic Drive Unit Force Control System

Each joint of hydraulic drive quadruped robot is driven by the hydraulic drive unit（HDU）, and the contacting between the robot foot end and the ground is complex and variable, which increases the difficulty of force control inevitably. In the recent years, although many scholars researched some control methods such as disturbance rejection control, parameter self-adaptive control, impedance control and so on, to improve the force control performance of HDU, the robustness of the force control still needs improving. Therefore, how to simulate the complex and variable load characteristics of the environment structure and how to ensure HDU having excellent force control performance with the complex and variable load characteristics are key issues to be solved in this paper. The force control system mathematic model of HDU is established by the mechanism modeling method, and the theoretical models of a novel force control compensation method and a load characteristics simulation method under different environment structures are derived, considering the dynamic characteristics of the load stiffness and the load damping under different environment structures. Then, simulation effects of the variable load stiffness and load damping under the step and sinusoidal load force are analyzed experimentally on the HDU force control performance test platform, which provides the foundation for the force control compensation experiment research. In addition, the optimized PID control parameters are designed to make the HDU have better force control performance with suitable load stiffness and load damping, under which the force control compensation method is introduced, and the robustness of the force control system with several constant load characteristics and the variable load characteristics respectively are comparatively analyzed by experiment. The research results indicate that if the load characteristics are known, the force control compensation method presented in this paper has positive compensation effects on the load characteristics variation, i.e., this method decreases the effects of the load characteristics variation on the force control performance and enhances the force control system robustness with the constant PID parameters, thereby, the online PID parameters tuning control method which is complex needs not be adopted. All the above research provides theoretical and experimental foundation for the force control method of the quadruped robot joints with high robustness.

An element-free Galerkin method for ground penetrating radar numerical simulation

An element-free Galerkin method(EFGM) is used to solve the two-dimensional(2D) ground penetrating radar(GPR)modelling problems, due to its simple pre-processing, the absence of elements and high accuracy. Different from element-based numerical methods, this approach makes nodes free from the elemental restraint and avoids the explicit mesh discretization. First, we derived the boundary value problem for the 2D GPR simulation problems. Second, a penalty function approach and a boundary condition truncated method were used to enforce the essential and the absorbing boundary conditions, respectively. A three-layered GPR model was used to verify our element-free approach. The numerical solutions show that our solutions have an excellent agreement with solutions of a finite element method(FEM). Then, we used the EFGM to simulate one more complex model to show its capability and limitations. Simulation results show that one obvious advantage of EFGM is the absence of element mesh, which makes the method very flexible. Due to the use of MLS fitting, a key feature of EFM, is that both the dependent variable and its gradient are continuous and have high precision.

Using modified Soil Conservation Service curve number method to simulate the role of forest in flood control in the upper reach of the Tingjiang River in China

To improve flood control efficiency and increase urban resilience to flooding,the impacts of forest type change on flood control in the upper reach of the Tingjiang River(URTR) were evaluated by a modified model based on the Soil Conservation Service curve number(SCS-CN) method. Parameters of the model were selected and determined according to the comprehensive analysis of model evaluation indexes. The first simulation of forest reconstruction scenario,namely a coniferous forest covering 59.35km^2 is replaced by a broad-leaved forest showed no significant impact on the flood reduction in the URTR. The second simulation was added with 61.75km^2 bamboo forest replaced by broad-leaved forest,the reduction of flood peak discharge and flood volume could be improved significantly. Specifically,flood peak discharge of 10-year return period event was reduced to 7-year event,and the reduction rate of small flood was 21%-28%. Moreover,the flood volume was reduced by 9%-14% and 18%-35% for moderate floods and small floods,respectively. The resultssuggest that the bamboo forest reconstruction is an effective control solution for small to moderate flood in the URTR,the effect of forest conversion on flood volume is increasingly reduced as the rainfall amount increases to more extreme magnitude. Using a hydrological model with scenarios analysis is an effective simulation approach in investigating the relationship between forest type change and flood control. This method would provide reliable support for flood control and disaster mitigation in mountainous cities.

Numerical simulation of hydraulic fracturing and associated microseismicity using finite-discrete element method

Hydraulic fracturing （HF） technique has been extensively used for the exploitation of unconventional oiland gas reservoirs. HF enhances the connectivity of less permeable oil and gas-bearing rock formationsby fluid injection, which creates an interconnected fracture network and increases the hydrocarbonproduction. Meanwhile, microseismic （MS） monitoring is one of the most effective approaches to evaluatesuch stimulation process. In this paper, the combined finite-discrete element method （FDEM） isadopted to numerically simulate HF and associated MS. Several post-processing tools, includingfrequency-magnitude distribution （b-value）, fractal dimension （D-value）, and seismic events clustering,are utilized to interpret numerical results. A non-parametric clustering algorithm designed specificallyfor FDEM is used to reduce the mesh dependency and extract more realistic seismic information.Simulation results indicated that at the local scale, the HF process tends to propagate following the rockmass discontinuities; while at the reservoir scale, it tends to develop in the direction parallel to themaximum in-situ stress. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.