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.

The Calculation of Stress-Strain State of Anisotropic Composite Finite-Element Area with Different Boundary Conditions on the Surface

The numerical analytic research approach of stress-strain state of anisotropic composite finite element area with different boundary conditions on the surface, is represented below. The problem is solved by using a spatial model of the elasticity theory. Differential equation system in partial derivatives reduces to one-dimensional problem using spline collocation method in two coordinate directions. Boundary problem for the system of ordinary higher-order differential equation is solved by using the stable numerical technique of discrete orthogonalization.

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

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

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.

Solution of shallow-water equations using least-squares finite-element method

A least-squares finite-element method （LSFEM） for the non-conservative shallow-water equations is presented. The model is capable of handling complex topography, steady and unsteady flows, subcritical and supercritical flows, and flows with smooth and sharp gradient changes. Advantages of the model include： （1） sources terms, such as the bottom slope, surface stresses and bed frictions, can be treated easily without any special treatment; （2） upwind scheme is no needed; （3） a single approximating space can be used for all variables, and its choice of approximating space is not subject to the Ladyzhenskaya-Babuska-Brezzi （LBB） condition; and （4） the resulting system of equations is symmetric and positive-definite （SPD） which can be solved efficiently with the preconditioned conjugate gradient method. The model is verified with flow over a bump, tide induced flow, and dam-break. Computed results are compared with analytic solutions or other numerical results, and show the model is conservative and accurate. The model is then used to simulate flow past a circular cylinder. Important flow charac-teristics, such as variation of water surface around the cylinder and vortex shedding behind the cylinder are investigated. Computed results compare well with experiment data and other numerical results.

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.

Least-squares finite-element method for shallow-water equations with source terms

Numerical solution of shallow-water equations （SWE） has been a challenging task because of its nonlinear hyperbolic nature, admitting discontinuous solution, and the need to satisfy the C-property. The presence of source terms in momentum equations, such as the bottom slope and friction of bed, compounds the difficulties further. In this paper, a least-squares finite-element method for the space discretization and θ-method for the time integration is developed for the 2D non-conservative SWE including the source terms. Advantages of the method include： the source terms can be approximated easily with interpolation functions, no upwind scheme is needed, as well as the resulting system equations is symmetric and positive-definite, therefore, can be solved efficiently with the conjugate gradient method. The method is applied to steady and unsteady flows, subcritical and transcritical flow over a bump, 1D and 2D circular dam-break, wave past a circular cylinder, as well as wave past a hump. Computed results show good C-property, conservation property and compare well with exact solutions and other numerical results for flows with weak and mild gradient changes, but lead to inaccurate predictions for flows with strong gradient changes and discontinuities.

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.

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.

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.

Mesomechanics Finite-element Method for Determining the Shear Strength of Mudded Intercalation Materials

A new method regarding mesomechanics finite-element research is proposed to predict the peak shear strength of mudded intercalation materials on a mesoscopic scale. Based on geometric and mechanical parameters, along with the strain failure criteria obtained by sample＇s deformation characteristics, uniaxial compression tests on the sample were simulated through a finite-element model, which yielded values consistent with the data from the laboratory uniaxial compression tests, implying that the method is reasonable. Based on this model, a shear test was performed to calculate the peak shear strength of the mudded intercalation, consistent with values reported in the literature, thereby providing a new approach for investigating the mechanical properties of mudded intercalation 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.