Random Green's Function Method for Large-Scale Electronic Structure Calculation

We report a linear-scaling random Green's function(rGF) method for large-scale electronic structure calculation. In this method, the rGF is defined on a set of random states and is efficiently calculated by projecting onto Krylov subspace. With the rGF method, the Fermi–Dirac operator can be obtained directly, avoiding the polynomial expansion to Fermi–Dirac function. To demonstrate the applicability, we implement the rGF method with the density-functional tight-binding method. It is shown that the Krylov subspace can maintain at small size for materials with different gaps at zero temperature, including H_(2)O and Si clusters. We find with a simple deflation technique that the rGF self-consistent calculation of H_(2)O clusters at T = 0 K can reach an error of~ 1 me V per H_(2)O molecule in total energy, compared to deterministic calculations. The rGF method provides an effective stochastic method for large-scale electronic structure simulation.

Utility and Application of a Versatile Analytical Method for MMF Calculation in AC Machines

A versatile analytical method(VAM) for calculating the harmonic components of the magnetomotive force(MMF) generated by diverse armature windings in AC machines has been proposed, and the versatility of this method has been established in early literature. However, its practical applications and significance in advancing the analysis of AC machines need further elaboration. This paper aims to complement VAM by augmenting its theory, offering additional insights into its conclusions, as well as demonstrating its utility in assessing armature windings and its application of calculating torque for permanent magnet synchronous machines(PMSM). This work contributes to advancing the analysis of AC machines and underscores the potential for improved design and performance optimization.

Discussion on Calculation Method of Poverty Incidence in the Exit Evaluation of Poverty-Stricken Counties and Villages

Poverty incidence is the key index that needs to be measured in the poverty exit examination and evaluation of 832 poverty-stricken counties and 128 000 poverty-stricken villages. In this paper, based on the statement of general concept and conventional calculation method of the poverty incidence, the calculation method of poverty incidence in the exit evaluation of poverty-stricken counties （also including poverty-stricken township and poverty-stricken villages） was investigated through the view of the third-party evaluation. In addition to considering the previous ＂number of planned poverty remaining population＂, the method also needed to give consideration to the exiting mistaken population, evaluation missing population. Based on the case in Yuanyang County, Honghe Hani and Yi Autonomous Prefecture, Yunnan Province, the poverty incidences of 10 exit planning villages by the end of 2017 in Yuanyang County were investigated and estimated, and suggestions were proposed to promote the precise poverty alleviation and poverty relief of the county.

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

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

Comparison study on calculation methods of bending behavior of Chinese reinforced concrete beams from 1912 to 1949

In order to study the calculation methods of bending behavior of Chinese reinforced concrete beams from 1912 to 1949, tests on the mechanical performance of 66 rebars from different modem Chinese concrete buildings, the concrete compressive strength of 12 modem Chinese concrete buildings, and the concrete cover thickness of 9 modem Chinese concrete buildings are carried out; and the actual material properties and structural conformations of modem Chinese concrete buildings are obtained. Then, the comparison on calculation methods of bending behavior including the original Chinese calculation method, the present Chinese calculation method, the present American calculation method and the present European calculation method is studied. The results show that the original Chinese calculation method of bending behavior is based on the allowable stress calculation method, and the design safety factors are 3.55 to 4. 0. In term of the calculation area of longitudinal rebars of reinforced concrete beams, without considering earthquake action, the original Chinese structural calculation method is safer than the present Chinese structural calculation method, the present European structural calculation method, and the present American structural calculation method. The results can provide support for the structural safety assessments of modem Chinese reinforced concrete buildings.

A Flexible-Segment-Model-Based Dynamics Calculation Method for Free Hanging Marine Risers in Re-Entry

In re-entry, the drilling riser hanging to the holding vessel takes on a free hanging state, waiting to be moved from the initial random position to the wellhead. For the re-entry, dynamics calculation is often done to predict the riser motion or evaluate the structural safety. A dynamics calculation method based on Flexible Segment Model （FSM） is proposed for free hanging marine risers. In FSM, a riser is discretized into a series of flexible segments. For each flexible segment, its deflection feature and external forces are analyzed independently. For the whole riser, the nonlinear governing equations are listed according to the moment equilibrium at nodes. For the solution of the nonlinear equations, a linearization iteration scheme is provided in the paper. Owing to its flexibility, each segment can match a long part of the riser body, which enables that good results can be obtained even with a small number of segments. Moreover, the linearization iteration scheme can avoid widely used Newton-Rapson iteration scheme in which the calculation stability is influenced by the initial points. The FSM-based dynamics calculation is timesaving and stable, so suitable for the shape prediction or real-time control of free hanging marine risers.

Application of Two-dimensional Viscous CE/SE Method in Calculation of Two-phase Detonation

The method of two-dimensional viscous space-time conservation element and solution element (CE/SE) can be used to calculate the gas-liquid two-phase interior flow field in pulse detonation engine (PDE). In this paper, the evolution of the detonation wave and the distribution of its physical parameters were analyzed. The numerical results show that the change of axial velocity of gas is the same as that of detonation pressure. The larger the liquid droplet radius is, the longer the time to get stable detonation wave is. The calculated results coincide with the experimented results better.

Simplified Calculation Methods for All-Vertical-Piled Wharf in Offshore Deep Water

All-vertical-piled wharf is a kind of high-piled wharf, but it is extremely different from the traditional ones in some aspects, such as the structural property, bearing characteristics, failure mechanism, and static or dynamic calculation methods. In this paper, the finite element method （FEM） and theoretical analysis method are combined to analyze the structural property, bearing behavior and failure mode of the all-vertical-piled wharf in offshore deep water, and to establish simplified calculation methods determining the horizontal static ultimate bearing capacity and the dynamic response for the all-vertical-piled wharf. Firstly, the bearing capability and failure mechanism for all-vertical-piled wharf are studied by use of FEM, and the failure criterion is put forward for all-vertical-piled wharf based on the ＇plastic hinge＇. According to the failure criterion and P-Y curve method, the simplified calculation method of the horizontal static ultimate bearing capacity for all-vertical-piled wharf is proposed, and it is verified that the simplified method is reasonable by comparison with the FEM. Secondly, the displacement dynamic magnification factor for the all-vertical-piled wharf under wave cyclic loads and ship impact loads is calculated by the FEM and the theory formula based on the single degree of freedom （SDOF） system. The results obtained by the two methods are in good agreement with each other, and the simplified calculation method of the displacement dynamic magnification factor for all-vertical-piled wharf under dynamic loads is proposed. Then the simplified calculation method determining the dynamic response for the all-vertical-piled wharf is proposed in combination with P-Y curve method. That is, the dynamic response of the structure can be obtained through the static calculation results of P-Y curve method multiplied by the displacement dynamic magnification factor. The feasibility of the simplified dynamic response method is verified by comparison with the FEM under different conditions.

Calculation of Phase Equilibria Based on the Levenberg-Marquardt Method

The Levenberg-Marquardt method, the best algorithm to obtain the least-square solution of nonlinear equations, is applied to calculate the stable phase equilibria. It can get the best combination between robustness and speed of the calculations. Its application to ternary AI-Si-Mg system is executed in detail. The calculated phase equilibria agree well with the experimental results. Furthermore, the Levenberg-Marquardt method is not sensitive to the initial values.

Resources calculation of cobalt-rich crusts with the grid subdivision and integral method

On the basis of three geological models and several orebody boundaries, a method of grid subdivision and integral has been proposed to calculate and evaluate the resources of cobalt-rich crusts on the seamounts in the central Pacific Ocean. The formulas of this method are deduced and the interface of program module is designed. The method is carried out in the software ＂Auto mapping system of submarine topography and geomorphology MBChart＂. This method and program will possibly become a potential tool to calculate the resources of seamounts and determine the target diggings for China＇ s next Five-year Plan.

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.

Analysis and management strategies for the low voltage problems of rural power grid based on the distribution network flow calculation method

For a long time, because of the lack of investment capital and enough attentions, the overall constructions of rural power grid were far behind than the urban power grid in Chongqing Jiangbei Power Company. The low voltage problems were highlighted in the rural power grid due to the characteristics of rural power grid. Using the distribution network flow calculation method, we evaluated the low voltage problems of the rural power grid which belongs to Chongqing Jiangbei Power Company. In addition, we collected the data of distribution transformers in electricity consumption peak period. Some practical management strategies were proposed by the analysis and evaluation of potential and appeared low voltage problems.

Precise integration method without inverse matrix calculation for structural dynamic equations

The precise integration method proposed for linear time-invariant homogeneous dynamic systems can provide accurate numerical results that approach an exact solution at integration points. However, difficulties arise when the algorithm is used for non-homogeneous dynamic systems due to the inverse matrix calculation required. In this paper, the structural dynamic equalibrium equations are converted into a special form, the inverse matrix calculation is replaced by the Crout decomposition method to solve the dynamic equilibrium equations, and the precise integration method without the inverse matrix calculation is obtained. The new algorithm enhances the present precise integration method by improving both the computational accuracy and efficiency. Two numerical examples are given to demonstrate the validity and efficiency of the proposed algorithm.

Parameter calculation and result storage for phase-field simulation in α-Mg dendrite growth of Mg-5-wt% Zn alloy

Parameter calculation and result storage, as two necessary steps in phase-field simulation play an important role in ensuring the accuracy of simulation results. A strategy of parameter calculation and result storage is presented for phase-field simulation in α-Mg dendrite growth of Mg-5-wt% Zn alloy under isothermal solidification. Based on the phase diagram and empirical formulas, key parameters of the phase-field model, such as equilibrium partition coefficient k, liquidus slope m, solutal diffusion coefficient in liquid Dl, and solutal diffusion coefficient in solid Ds, can be obtained.Both structured grid method and structured point method can be used to store simulation results, but using the latter method will reduce about 60% storage space and 37.5% storage time compared with the former. Finally, convergent simulation results of α-Mg dendrite growth are obtained and they are in good agreement with the experimental results about optical micrograph, which verify the accuracy of parameters and stability of storage method.

Comparative Studies on Calculation Methods of Shear Behavior of Reinforced Concrete Beams of Chinese Modern Reinforced Concrete Buildings

In order to study the calculation methods of shear behavior of reinforced concrete beams of Chinese modern reinforced concrete buildings,this paper carried out tests on the concrete compressive strength of 12 Chinese modern concrete buildings,the mechanical properties of 66 rebars from different Chinese modern concrete buildings,and the concrete cover thickness of 9 Chinese modern concrete buildings,and the actual material properties and structural configurations have been obtained. Then,the comparison on calculation methods include the Chinese original calculation method,the Chinese present calculation method,the American present calculation method and the European present calculation method is studied with case analysis method.The results show that the Chinese original calculation method of shear behavior of reinforced concrete beams is based on the allowable stress calculation method,and the design safety factors are 3. 55- 4. 00. The standard value of the compressive strength of concrete cubes is 8. 48 MPa,the standard value of the concrete tensile strength is 1.20 MPa,the standard value of the yield strength of rectangular rebars is 229. 56 MPa,and the standard value of the yield strength of round rebars is 276. 82 MPa. The average value of the concrete cover thickness of beams and columns is 35.96 mm. In term of calculation area of hoop rebars of reinforced concrete beams,without considering earthquake loads,the Chinese original structural calculation method is safer than the Chinese present structural calculation method,but is more unsafe than the American present structural calculation method and the European present structural calculation method. The results can provide the support for structural safety assessments and repair designs of Chinese modern reinforced concrete buildings.

Hybrid N-order Lagrangian Interpolation Eulerian-Lagrangian Method for Salinity Calculation

The Eulerian？Lagrangian method（ELM） has been used by many ocean models as the solution of the advection equation,but the numerical error caused by interpolation imposes restriction on its accuracy.In the present study,hybrid N-order Lagrangian interpolation ELM（Li ELM） is put forward in which the N-order Lagrangian interpolation is used at first,then the lower order Lagrangian interpolation is applied in the points where the interpolation results are abnormally higher or lower.The calculation results of a step-shaped salinity advection model are analyzed,which show that higher order（N=3？8） Li ELM can reduce the mean numerical error of salinity calculation,but the numerical oscillation error is still significant.Even number order Li ELM makes larger numerical oscillation error than its adjacent odd number order Li ELM.Hybrid N-order Li ELM can remove numerical oscillation,and it significantly reduces the mean numerical error when N is even and the current is in fixed direction,while it makes less effect on mean numerical error when N is odd or the current direction changes periodically.Hybrid odd number order Li ELM makes less mean numerical error than its adjacent even number order Li ELM when the current is in the fixed direction,while the mean numerical error decreases as N increases when the current direction changes periodically,so odd number of N may be better for application.Among various types of Hybrid N-order Li ELM,the scheme reducing N-order directly to 1st-order may be the optimal for synthetic selection of accuracy and computational efficiency.

The Method for Measuring the Loop Resistance and the Temperature Rise Calculation of GIS Conductor Pole

This paper studies the method for measuring the loop resistance of GIS conductor pole based on the super capacitor producing impulse current up to several thousand amperes. This method overcomes the limitations of conventional diagnostic method. Typical GIS conductor poles are chosen. Based on FEA and lab tests, the effect of different forms of current and contact condition, relationship between the temperature of contact and the loop resistance is researched. In full- scale testing under realistic operating conditions on the new 220 kV GIS using prototype instrumentation a very good sensitivity in an early stage was obtained.

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.

Thermal kinetic analysis of a complex process from a solid-state reaction by deconvolution procedure from a new calculation method and related thermodynamic functions of Mn_(0.90)Co_(0.05)Mg_(0.05)HPO_4?3H_2O

Three individual peaks of thermal solid-state reaction processes of the synthesized Mn0.90Co0.05Mg0.05HPO4?3H2O were observed corresponding to dehydration I,dehydration II and polycondensation processes.An alternative method for the calculation of the extent of conversion was proposed from the peak area of the individual DTG peak after applying the best fitting deconvolution function(Frazer–Suzuki function).An iterative integral isoconversional equation was used to compute the values of the apparent activation energy Eαand they were found to be 65.87,78.16 and 119.32 kJ/mol for three peaks,respectively.Each individual peak was guaranteed to be a single-step kinetic system with its unique kinetic parameters.The reaction mechanism functions were selected by the comparison between experimental and model plots.The results show that the first,second and final individual peaks were two-dimensional diffusion of spherical symmetry(D2),three-dimensional diffusion of spherical symmetry(D3)and contracting cylinder(cylindrical symmetry,R2)mechanisms.Pre-exponential factor values of 3.91×106,1.35×107 and 2.15×107 s?1 were calculated from the Eαvalues and reaction mechanisms.The corresponded standard thermodynamic functions of the transition-state(activated)complexes were determined and found to agree well with the experimental data.

A Two-Stage Method for an Approximate Calculation of Statically Indeterminate Trusses

The paper presents the principles of a method, which in two simple stages makes possible to carry out the statically calculation of values of forces acting in the fiat static indeterminate trusses. In each stage, it is considered the static determinate truss, scheme of which is obtained after remove the suitable number of members from the basic static indeterminate truss. The both intermediate statically determinate trusses are of the same clear span and they are loaded by forces of half values applied to the corresponding truss nodes. The method applies one of the typical procedures of calculation of the statically determinate trusses and then it is applied in an appropriate way the rule of superposition for obtaining the final values of forces acting in particular members of the basic truss. The values of forces calculated in this way are of a very close approximation to the force values determined in the special and complex ways being considered as the exact calculation methods. The proposed method can be useful mostly but not only for the initial structural design of such systems. The simplicity of the two-stage method justifies an assumption that it can be relatively easy and worthy to adjust to the requirements of the computer aided technology of statically calculation of the complex forms of trusses.