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.

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.

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.

Application of neural network to speed-up equilibrium calculations in compositional reservoir simulation

Compositional reservoir simulation is an important tool to model fluid flow in oil and gas reservoirs.Important investment decisions regarding oil recovery methods are based on simulation results,where hundred or even thousand of different runs are performed.In this work,a new methodology using artificial intelligence to learn the thermodynamic equilibrium is proposed.This algorithm is used to replace the classical equilibrium workflow in reservoir simulation.The new method avoids the stability test for single-phase cells in most cases and provides an accurate two-phase flash initial estimate.The classical and the new workflow are compared for a gas-oil mixing case,showing a simulation time speed-up of approximately 50%.The new method can be used in compositional reservoir simulations.

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.

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.

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.

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.

Theoretical and numerical simulation of critical gas supply of refuge chamber

The personnel in refuge chamber absorb O_2 and exhale CO_2 all the time. Supplying O_2 and removing CO_2 are the basic function of refuge chamber. After disaster occurs, the supply of the compressed air or oxygen for personnel in refuge chamber is limited. Thus, how to effectively use the compressed air and oxygen and try to improve the time of supply has a great significance. Supplying more oxygen will result in waste, while supplying less oxygen will cause its concentration to be lower, and harm life safety. This research uses the theoretical calculation and numerical simulation, finds critical gas supply for refuge chamber, and illuminates the change law of gas concentration with critical gas supply in refuge chamber,which provides theoretical guidance for effective use of compressed air and oxygen in refuge chamber.

Air flow calculation of reverse circulation drilling technique with air DTH hammer

Air DTH hammer has been successfully applied in minor-caliber solid mineral exploration,water-well drilling and other drilling areas. In order to expand the applications of the technology,the authors further studied the principle and analyzed the mechanism of reverse circulation drilling technique with air DTH hammer to get the perfect assembles of equipments by optimizing working parameters. No parameter seemed more important than the air volume because it could maintain the working performance stability. The minimum air volume is related to the parameters such as depth and pressure,which was calculated under the actual conditions. It was solved for the air injection flow tables of the air DTH Hammer working at the different pressures. According to the data tables,operators could adjust the air volume to meet the demand on this technique,which had a realistic guiding significance. So it could build up a set of systematic and complete hi-technique.

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.

A FUNCTION MINIMIZATION CALCULATION OF CYCLONE 10 MAGNETIC FIELDS

By app(?)ing function minimization calculation method, two function expressions are used to simulate the magnetic field measured for cyclotron Cyclone 10 in azimuth and radius The numerical fitting curves are consistent with magnetic field measured. In most pl(?)es, the accuracies are several thousandth, except those errors to be pointed out in paper.

Finite Element Simulation of Processes Involving Moving Heat Sources. Application to Welding and Surface Treatment

A wide range of welding and surface treatment processes involve the use of a heat source which is moving at a constant speed over the component. The numerical simulation of such processes implies a transient analysis using a very refined mesh in order to follow properly the path of the heat source. The 3D-mesh size can be very large if one consider the welds length or the heat-treated surface size in industrial components. To reduce the computational time to acceptable values, several techniques have been investigated. The first type is to use analytical methods such as Rosenthal equations. The second type of solutions consists in performing a transient analysis using adaptive meshing. But, for a large proportion of the involved processes, practical experience demonstrates the existence of quasi steady state conditions over the major part of the heat source path. Numerical algorithms have therefore been developed to directly compute the steady temperature, metallurgical phase proportion and stress distributions. This paper gives a general overview of the different numerical methods used to simulate welding and surface treatment processes with a special emphasis on the steady state calculation. The benefits and limitations of each of them are discussed and applications are presented.

Cooling Water Flow Rate Calculation for Hearth of Large Blast Furnaces

The cooling water flow rate for hearth of large blast furnaces was calculated by simulation. The results show that the cooling water flow rate shall be above 4 200m3/ h for hearth of large blast furnaces; to meet requirements of the increasing smelting intensity and to ensure the safety at the end of the first campaign,the designed maximum cooling water flow rate should be 5 900m3/ h; according to the flow distribution stability and the calculated resistance loss,hearth cooling stave pipes with the specification of 76 mm × 6 mm shall be adopted to assure the flow velocity in pipes of hearth cooling stave in the range of 1. 9- 2. 3 m / s.

Measurement and Simulation of Induced Voltage and Current on 110 KV Crossing Transmission Lines under UHV AC Transmission Lines

The induced electricity of 110 kV transmission lines which cross the UHV AC transmission lines may threaten personal safety of the maintenance staff. In this paper, field measurement of the induced voltage and induced current on a 110 kV crossing line inside Jinhua in Zhejiang province is performed. The electrostatic induced voltage on the measured line is 12.24 kV. The power frequency electromagnetic field simulation model is established, and the calculation results are consistent with the measured. Influence factors analysis shows that the electrostatic induced voltage on the 110 kV line is 12.78 kV, the electromagnetic induced voltage is 12.3 V, the induced current through ground wire is less than 1A when the UHV lines operate at full load. The induced voltage and current decrease while the crossing distance increases. Parallel lines induction is much higher than crossing lines. The electromagnetic induced voltage after ground knife-switch shut down would exceed the human safety voltage 36 V while the crossing angle is less than 30?, so the temporary ground wire must be hanged to ensure safety of the maintenance staff.