A review of monitoring,calculation,and simulation methods for ground subsidence induced by coal mining

Subsidence data acquisition methods are crucial to mining subsidence research and an essential component of achieving the goal of environmentally friendly coal mining.The origin and history of the existing methods of field monitoring,calcula-tion,and simulation were introduced.It summarized and analyzed the main applications,flaws and solutions,and improve-ments of these methods.Based on this analysis,the future developing directions of subsidence data acquisition methods were prospected and suggested.The subsidence monitoring methods have evolved from conventional ground monitoring to combined methods involving ground-based,space-based,and air-based measurements.While the conventional methods are mature in technology and reliable in accuracy,emerging remote sensing technologies have obvious advantages in terms of reducing field workload and increasing data coverage.However,these remote sensing methods require further technological development to be more suitable for monitoring mining subsidence.The existing subsidence calculation methods have been applied to various geological and mining conditions,and many improvements have already been made.In the future,more attention should be paid to unifying the studies of calculation methods and mechanical principles.The simulation methods are quite dependent on the similarity of the model to the site conditions and are generally used as an auxiliary data source for subsidence studies.The cross-disciplinary studies between subsidence data acquisition methods and other technologies should be given serious consideration,as they can be expected to lead to breakthroughs in areas such as theories,devices,software,and other aspects.

Equivalent Method of Integrated Power Generation System of Wind, Photovoltaic and Energy Storage in Power Flow Calculation and Transient Simulation

针对工程实际开展风光储联合发电系统在潮流计算和机电暂态仿真中的等值方法研究，旨在为大容量风光储联合发电系统的并网仿真分析奠定基础。将潮流计算的等值分为单元机组和集电系统2部分来研究。单元机组等值采用根据不同控制模式选取不同节点类型的方法，针对集电系统等值提出基于损耗不变原则的方法。等值模型和详细模型的算例结果表明，潮流计算等值方法具有较好的精度。在机电暂态仿真动态等值中，基于实际工程计算的最严重工况分析原则，提出运行在满出力点的单机“倍乘”等值模型，为工程计算中的风光储联合发电系统动态等值提供了一种解决方案。

Simulation and Techno-Economic Performance of a Novel Charge Calculation and Melt Optimization Planning Model for Steel Making

Process algorithm, numerical model and techno-economic assessment of charge calculation and furnace bath optimization for target alloy for induction furnace-based steelmaking is presented in this study. The developed algorithm combines the make-to-order (MTO) and charge optimization planning (COP) of the steel melting shop in the production of target steel composition. Using a system-level approach, the unit operations involved in the melting process were analyzed with the purpose of initial charge calculation, prevailing alloy charge prediction and optimizing the sequence of melt chemistry modification. The model performance was established using real-time production data from a cast iron-based foundry with a 1- and 2-ton induction furnace capacity and a medium carbon-based foundry with a 10- and 15-ton induction furnace capacity. A simulation engine (CastMELT) was developed in Java IDE with a MySQL database for continuous interaction with changing process parameters to run the model for validation. The comparison between the model prediction and production results was analyzed for charge prediction, melt modification and ferroalloy optimization and possible cost savings. The model performance for elemental charge prediction and calculation purpose with respect to the charge input (at overall scrap meltdown) gave R-squared, Standard Error, Pearson correlation and Significance value of (0.934, 0.06, 0.97, 0.0003) for Carbon prediction, (0.962, 0.06, 0.98, 0.00009) for Silicon prediction, (0.999, 0.048, 0.999, 9E -11) for Manganese Prediction, and (0.997, 0.076, 0.999, 6E -7) for Chromium prediction respectively. Correlation analysis for melt modification (after charging of ferroalloy) using the model for after-alloying spark analysis compared with the target chemistry is at 99.82%. The results validate the suitability of the developed model as a functional system of induction furnace melting for combined charge calculation and melt optimization Techno-economic evaluation results showed that 0.98% - 0.25% ferroalloy saving per ton of melt is possible using the model. This brings about an annual production cost savings of 100,000 $/y in foundry A (medium carbon steel) and 20,000 $/y in foundry B (cast iron) on the use of different ferroalloy materials.

Calculation Method of Passive Residual Heat Removal Heat Exchanger and Numerical Simulation

The tube inside and outside heat transfer mechanism of Passive Residual Heat Removal Heat Exchanger (PRHR HX) was analyzed. The calculation method of this special heat exchanger under natural convection condition in In-containment Refueling Water Storage Tank (IRWST) was carried out. The single-tube coupling model three-dimensional natural circulation in the IRWST was simulated numerically using Fluent. The heat transfer and flow characteristics of the fluid in IRWST were obtained. The comparison of the results between theoretical arithmetic and numerical simulation showed that the theoretical calculation method is suitable for the heat transfer calculation of PRHR HX.

Phase equilibrium of Cd_(1-x)Zn_xS alloys studied by first-principles calculations and Monte Carlo simulations

The first-principles calculations based on density functional theory combined with cluster expansion techniques and Monte Carlo （MC） simulations were used to study the phase diagrams of both wurtzite （WZ） and zinc-blende （ZB） Cdl_xZnxS alloys. All formation energies are positive for WZ and ZB Cdl-xZnxS alloys, which means that the Cdl-xZnxS alloys are unstable and have a tendency to phase separation. For WZ and ZB Cdl_xZnxS alloys, the consolute temperatures are 655 K and 604 K, respectively, and they both have an asymmetric miscibility gap. We obtained the spatial distributions of Cd and Zn atoms in WZ and ZB Cd0.sZn0.sS alloys at different temperatures by MC simulations. We found that both WZ and ZB phases of Cdo.sZn0.sS alloy exhibit phase segregation of Cd and Zn atoms at low temperature, which is consistent with the phase diagrams.

STUDY ON MESO-MECHANICAL SIMULATION OF GROUTING FOR UNFAVORABLE GEOLOGICAL CONDITIONS

The numerical simulation program of PFC2D(Particle Flow Code in 2 Dimension)particle flow based on the flow-solid coupling principle and,on its built-in FISHTANK function library and FISH language,defines the flow equation and pressure equation of fluid domain respectively,and carries out numerical simulation calculations on the diffusion process and,on the morphology and particle displacement of slurry during the slurry injection process.By adjusting the parameters of hist,n_bond,s_bond and measure in the PFC command flow,the tracking of granular body displacement changes is achieved,and the mesoscopic mechanism such as the diffusion law of soil slurry at different depths and the change of formation porosity is revealed.The numerical calculations show that:the grouting pressure has a significant effect on the alteration and destruction of the formation structure,and the fracturing effect becomes gradually worse with increasing adhesive strength,while the porosity increases significantly with increasing grouting pressure.Based on the elastic-plastic theory of the Mohr-Colomb criterion to theoretically derive the stress field of the soil around the borehole,it is pointed out that the mechanical mechanism of annular tension and radial compression is the fundamental reason for the appearance of fracturing grouting action mode.The increase of slurry viscosity is beneficial to improve the grouting effect of fracturing-compacting grouting,while the increase of friction coefficient has little effect on the grouting effect.The comparative analysis of the laboratory tests shows that the PFC2D simulation of the grouting process is feasible.

Temperature Field in Laser Line Scanning Thermography: Analytical Calculation and Experiment

The temperature field in laser line scanning thermography is investigated comprehensively in this work,including analytical calculation and experiment.Firstly,the principle of laser line scanning thermography is analyzed.On this basis,a physical laser line scanning model is proposed.Afterwards,based on Fourier transform(FT)and segregation variablemethod(SVM),the heat conduction differential equation in laser line scanning thermography is derived in detail.The temperature field of the composite-based coatings model with defects is simulated numerically.The results show that the laser line scanning thermography can effectively detect the defects in the model.The correctness of the analytical calculation is verified by comparing the surface temperature distribution obtained by analytical calculation and numerical simulation.Additionally,an experiment is carried out and the changeable surface temperature obtained by analytical calculation is compared with the experimental results.It shows that the error of maximum temperature obtained by analytical calculation and by experiment is 8%with high accuracy,which proves the correctness of analytical calculation and enriches the theoretical foundation of laser line scanning thermography.

Simulating Calculation of Frictional Force between Mold and Billet

Assumption about the inner surface of mold is made according to the forming mechanisms of frictional force. The frictional force between billet and mold can be analysed by calculating the thermal shrinkage and the shearing strength of the meniscus.The model is used to calculate the maximum drawing speed at different conditions. The results are very satisfactory.

Computer Simulating Calculation on the Microstructure Evolutions during Hot Strip Rolling of C-Mn Steels

The program to predict the microstructure evolutions during hot strip rolling of C-M n steels has been developed in this paper, BV using this program, the microstructure changes with the processing parameters were analysed in detail. showing not only a good agreement of prediction with the measured values, but also entirely possibility to optimize hot strip rolling precess by computer simulation

Monte Carlo simulation based thinning and calculating method for noninvasive blood glucose sensing by near-infrared spectroscopy

Previous results show that the floating reference theory(FRT)is an effective tool to reduce the infuence of interference factors on noninvasive blood glucose sensing by near infrared spectros-copy(NIRS).It is the key to measure the floating reference point(FRP)precisely for the application of FRT.Monte Carlo(MC)simulation has been introduced to quantitatively in-vestigate the effects of positioning errors and light source drifts on measuring FRP.In this article,thinning and calculating method(TCM)is proposed to quantify the positioning error.Mean-while,the normalization process(NP)is developed to significantly reduce the error induced by light source drift.The results according to TCM show that 7 purm deviations in positioning can generate about 10.63%relative error in FRP.It is more noticeable that 1%fluctuation in light source intensity may lead to 12.21%relative errors.Gratifyingly,the proposed NP model can effectively reduce the error caused by light source drift.Therefore,the measurement system for FRPs must meet that the positioning error is less than 7 purm,and the light source drift is kept within 1%.Furthermore,an improvement for measurement system is proposed in order to take advantage of the NP model.

Numerical simulation of the blast wave of a multilayer composite charge

The numerical simulation of a blast wave of a multilayer composite charge is investigated.A calculation model of the near-field explosion and far-field propagation of the shock wave of a composite charge is established using the AUTODYN finite element program.Results of the near-field and far-field calculations of the shock wave respectively converge at cell sizes of 0.25-0.5 cm and 1-3 cm.The Euler--fluxcorrected transport solver is found to be suitable for the far-field calculation after mapping.A numerical simulation is conducted to study the formation,propagation,and interaction of the shock wave of the composite charge for different initiation modes.It is found that the initiation mode obviously affects the shock-wave waveform and pressure distribution of the composite charge.Additionally,it is found that the area of the overpressure distribution is greatest for internal and external simultaneous initiation,and the peak pressure of the shock wave exponentially decays,fitting the calculation formula of the peak overpressure attenuation under different initiation modes,which is obtained and verified by experiment.The difference between numerical and experimental results is less than 10%,and the peak overpressure of both internal and external initiation is 56.12% higher than that of central single-point initiation.

A NOVEL SLIGHTLY COMPRESSIBLE MODEL FOR LOW MACH NUMBER PERFECT GAS FLOW CALCULATION

By analyzing the characteristics of low Mach number perfect gas flows, a novel Slightly Compressible Model (SCM) for low Mach number perect gas flows is derived. In view of numerical calculations, this model is proved very efficient, for it is kept within thep-v frame but does not have to satisfy the time consuming divergence-free condition in order to get the incompressible Navier-Stokes equation solution. Writing the equations in the form of conservation laws, we have derived the characteristic systems which are necessary for numerical calculations. A cell-centered finite-volume method with flux difference upwind-biased schemes is used for the equation solutions and a new Exact Newton Relaxation (ENR) implicit method is developed. Various computed results are presented to validate the present model. Laminar flow solutions over a circular cylinder with wake developing and vortex shedding are presented. Results for inviscid flow over a sphere are compared in excellent agreement with the exact analytic incompressible solution. Three-dimensional viscous flow solutions over sphere and prolate spheroid are also calculated and compared well with experiments and other incompressible solutions. Finally, good convergent performances are shown for sphere viscous flows.

Determination of the anion O_3^- geometry by Franck- Condon simulations of its photoelectron spectrum

A theoretical method to calculate multidimensional Franck-Condon factors including Duschinsky effects is described and used to simulate the photoelectron spectroscopy of the anion O-3. Geometry optimization and harmonic vibrational frequency calculations have been performed on the AX～U4（）1A1 state of O3 and AX～U4（）2B1 state of O-3. Franck-Condon analyses and spectral simulation were carried out on the first photoelectron band of O-3. The theoretical spectrum obtained by employing CCSD（T）/6-311+G（2d,p） values are in excellent agreement with the observed one. In addition, the equilibrium geometry parameters, re（OO）= 0.135?5±0.000?5 nm and θe（O-O-O） =114.5±0.5°, of the AX～U4（）2B1 state of O-3, are derived by employing an iterative Franck-Condon analysis procedure in the spectral simulation.

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.

Theoretical and numerical simulation study on jet formation and penetration of different liner structures driven by electromagnetic pressure

The use of a shaped liner driven by electromagnetic force is a new means of forming jets. To study the mechanism of jet formation driven by electromagnetic force, we considered the current skin effect and the characteristics of electromagnetic loading and established a coupling model of "ElectriceMagnetic eForce" and the theoretical model of jet formation under electromagnetic force. The jet formation and penetration of conical and trumpet liners have been calculated. Then, a numerical simulation of liner collapse under electromagnetic force, jet generation, and the stretching motion were performed using an ANSYS multiphysics processor. The calculated jet velocity, jet shape, and depth of penetration were consistent with the experimental results, with a relative error of less than 10%. In addition, we calculated the jet formation of different curvature trumpet liners driven by the same loading condition and obtained the influence rule of the curvature of the liner on jet formation. Results show that the theoretical model and the ANSYS multiphysics numerical method can effectively calculate the jet formation of liners driven by electromagnetic force, and in a certain range, the greater the curvature of the liner is, the greater the jet velocity is.

Implementation of high-fidelity neutronics and thermal–hydraulic coupling calculations in HNET

To perform nuclear reactor simulations in a more realistic manner,the coupling scheme between neutronics and thermal-hydraulics was implemented in the HNET program for both steady-state and transient conditions.For simplicity,efficiency,and robustness,the matrixfree Newton/Krylov(MFNK)method was applied to the steady-state coupling calculation.In addition,the optimal perturbation size was adopted to further improve the convergence behavior of the MFNK.For the transient coupling simulation,the operator splitting method with a staggered time mesh was utilized to balance the computational cost and accuracy.Finally,VERA Problem 6 with power and boron perturbation and the NEACRP transient benchmark were simulated for analysis.The numerical results show that the MFNK method can outperform Picard iteration in terms of both efficiency and robustness for a wide range of problems.Furthermore,the reasonable agreement between the simulation results and the reference results for the NEACRP transient benchmark verifies the capability of predicting the behavior of the nuclear reactor.

Free energy calculation of single molecular interaction using Jarzynski's identity method:the case of HIV-1 protease inhibitor system

Jarzynski＇ identity （JI） method was suggested a promising tool for reconstructing free energy landscape of biomolecular interactions in numerical simulations and ex- periments. However, JI method has not yet been well tested in complex systems such as ligand-receptor molecular pairs. In this paper, we applied a huge number of steered molec- ular dynamics （SMD） simulations to dissociate the protease of human immunodeficiency type I virus （HIV-1 protease） and its inhibitors. We showed that because of intrinsic com- plexity of the ligand-receptor system, the energy barrier pre- dicted by JI method at high pulling rates is much higher than experimental results. However, with a slower pulling rate and fewer switch times of simulations, the predictions of JI method can approach to the experiments. These results sug- gested that the JI method is more appropriate for reconstruct- ing free energy landscape using the data taken from experi- ments, since the pulling rates used in experiments are often much slower than those in SMD simulations. Furthermore, we showed that a higher loading stiffness can produce higher precision of calculation of energy landscape because it yields a lower mean value and narrower bandwidth of work distri- bution in SMD simulations.

Improvement in Time Efficiency in Numerical Simulation for Solid Propellant Rocket Motors(SPRM)

The main purpose of the present work is to study the possibilities of reducing calculation time while maintaining the validity in the numerical simulation of the combustion product flow in SPRM chamber.Three ways of decreasing the calculation time-the use of numerical methods of high accuracy order,the reduction in spatial dimension of the problem,and the use of physical features of the processes in SPRM chamber while constructing a calculation model-were considered.Presented calculation data show that the use of these approaches makes it possible to reduce the time for solving the problems of SPRM simulation significantly(up to 100times).Also conclusions about the applicability of the mentioned above approaches in SPRM design were made.

Study on Surface Electric Field Simulation of High Voltage Transmission Line Assembly

For the studies in the field of high voltage power transmission, this paper has adopted the method of finite element node potential, and put forward two kinds of high pressure sensor-fixture modeling scheme for the sensor-fixture of the high voltage side, the simulation analysis shows that the sensor-fixture surface should be smooth, and should not appear the conclusion of edges and corners. While through establishing the four clamps assembly optimized model, and simulates the strain gages, fixtures and conductor surface field strength and electric field distribution in the model as a whole in turn, this paper Finally got the optimal size of fixture structure and assembly of each part reasonable location layout.