Harmonics analysis of the ITER poloidal field converter based on a piecewise method

Poloidal field（PF） converters provide controlled DC voltage and current to PF coils. The many harmonics generated by the PF converter flow into the power grid and seriously affect power systems and electric equipment. Due to the complexity of the system, the traditional integral operation in Fourier analysis is complicated and inaccurate. This paper presents a piecewise method to calculate the harmonics of the ITER PF converter. The relationship between the grid input current and the DC output current of the ITER PF converter is deduced. The grid current is decomposed into the sum of some simple functions. By calculating simple function harmonics based on the piecewise method, the harmonics of the PF converter under different operation modes are obtained.In order to examine the validity of the method, a simulation model is established based on Matlab/Simulink and a relevant experiment is implemented in the ITER PF integration test platform.Comparative results are given. The calculated results are found to be consistent with simulation and experiment. The piecewise method is proved correct and valid for calculating the system harmonics.

3D Multiphase Piecewise Constant Level Set Method Based on Graph Cut Minimization

Segmentation of three-dimensional(3D) complicated structures is of great importance for many real applications.In this work we combine graph cut minimization method with a variant of the level set idea for 3D segmentation based on the Mumford-Shah model.Compared with the traditional approach for solving the Euler-Lagrange equation we do not need to solve any partial differential equations.Instead,the minimum cut on a special designed graph need to be computed.The method is tested on data with complicated structures.It is rather stable with respect to initial value and the algorithm is nearly parameter free.Experiments show that it can solve large problems much faster than traditional approaches.

Breakage Distribution Estimation of Bauxite Based on Piecewise Linearized Breakage Rate

Laboratory tests were carried out to study the breakage kinetics of diasporic bauxite and determine its breakage distribution function. Non-first order breakage with different deceleration rates for different size intervals is found, which is most probably caused by the heterogeneity of the ore. Piecewise linearization method is proposed to describe the non-first order breakage according to its characteristics. In the method, grinding time is divided into several intervals and breakage is assumed to be first order in each interval. So, the breakage rates are calculated by taking the product of the last interval as feed and then established as a function of particle size and grinding time. Based on the predetermined breakage rate function, the breakage distribution of the ore is back-calculated from the experimental data using the population balance model (PBM). Finally, the obtained breakage parameters are validated and the simulated data are in good agreement with the experimental data. The obtained breakage distribution and the method for breakage rate description are both significant for modeling the full scale ball milling process of bauxite.

Methods of de-noising the low frequency electromagnetic data

The quality of the low frequency electromagnetic data is affected by the spike and the trend noises.Failure in removal of the spikes and the trends reduces the credibility of data explanation.Based on the analyses of the causes and characteristics of these noises,this paper presents the results of a preset statistics stacking method(PSSM)and a piecewise linear fitting method(PLFM)in de-noising the spikes and trends,respectively.The magnitudes of the spikes are either higher or lower than the normal values,which leads to distortion of the useful signal.Comparisons have been performed in removing of the spikes among the average,the statistics and the PSSM methods,and the results indicate that only the PSSM can remove the spikes successfully.On the other hand,the spectrums of the linear and nonlinear trends mainly lie in the low frequency band and can change the calculated resistivity significantly.No influence of the trends is observed when the frequency is higher than a certain threshold value.The PLSM can remove effectively both the linear and nonlinear trends with errors around 1% in the power spectrum.The proposed methods present an effective way for de-noising the spike and the trend noises in the low frequency electromagnetic data,and establish a research basis for de-noising the low frequency noises.

Two-dimensional face stability analysis in rock masses governed by the Hoek-Brown strength criterion with a new multi-horn mechanism

The face stability problem is a major concern for tunnels excavated in rock masses governed by the Hoek-Brown strength criterion.To provide an accurate prediction for the theoretical solution of the critical face pressure,this study adopts the piecewise linear method(PLM)to account for the nonlinearity of the strength envelope and proposes a new multi-horn rotational mechanism based on the Hoek-Brown strength criterion and the associative flow rule.The analytical solution of critical support pressure is derived from the energy-work balance equation in the framework of the plastic limit theorem;it is formulated as a multivariable nonlinear optimization problem relying on 2m dependent variables(m is the number of segments).Meanwhile,two classic linearized measures,the generalized tangential technique(GTT)and equivalent Mohr-Coulomb parameters method(EMM),are incorporated into the analysis for comparison.Surprisingly,the parametric study indicates a significant improvement in support pressure by up to 13%compared with the GTT,and as expected,the stability of the tunnel face is greatly influenced by the rock strength parameters.The stress distribution on the rupture surface is calculated to gain an intuitive understanding of the failure at the limit state.Although the limit analysis is incapable of calculating the true stress distribution in rock masses,a rough approximation of the stress vector on the rupture surface is permitted.In the end,sets of normalized face pressure are provided in the form of charts for a quick assessment of face stability in rock masses.

Piecewise Spectral Collocation Method for Second Order Volterra Integro-Differential Equations with Nonvanishing Delay

In this paper,the piecewise spectral-collocation method is used to solve the second-order Volterra integral differential equation with nonvanishing delay.In this collocation method,the main discontinuity point of the solution of the equation is used to divide the partitions to overcome the disturbance of the numerical error convergence caused by the main discontinuity of the solution of the equation.Derivative approximation in the sense of integral is constructed in numerical format,and the convergence of the spectral collocation method in the sense of the L¥and L2 norm is proved by the Dirichlet formula.At the same time,the error convergence also meets the effect of spectral accuracy convergence.The numerical experimental results are given at the end also verify the correctness of the theoretically proven results.

An Improved Power Flow Method to Cope with Non-smooth Constraints of Integrated Energy Systems

Integrated energy system applications can significantly improve energy efficiency.In this paper,we establish an integrated energy system containing heat,electricity and gas.The existing power flow(PF)calculation method applied to integrated energy systems(IESs)does not consider non-smooth constraints,such as the piecewise pipeline friction coefficient and generator buses reactive power limits,etc.Mixed integer nonlinear programming(MINLP)is conventionally used to deal with piecewise pipeline friction coefficients in gas network parts,but it is both complex and inefficient.Hence,we develop a piecewise linear function-based fitting method that can reduce the number of integer variables and enhanced the computational efficiency.In the electric network part,if the reactive power of the PV bus violates limits,it will be converted into a PQ bus,which is a non-differentiable and non-smooth constraint.Mixed complementarity problems are conventionally introduced to represent the PV-PQ buses type switching relationship and are addressed by the Newton-Raphson(NR)method.However,the above method is sensitive to the initial point.Here,we introduce a robust projected Levenberg-Marquardt(PLM)algorithm to cope with this issue.We demonstrate the advantages of our method and validate it both in a small-scale system and largescale network test cases.

A GENERATOR AND A SIMPLEX SOLVER FOR NETWORK PIECEWISE LINEAR PROGRAMS

This is a brief report on our recent work in network piecewise linear programming (NPLP),and it consists of two parts. In the first park, we describe a generator for NPLP problems which is derived from the classical network linear program generator NETGEN. The generator creates networks of the same topological structures as NETGEN, but each arc is associated with a convex piecewise linear cost. The purpose of this program is to provide a set of standard test problems which can be used to compare the performance of various algorithms for NPLP. In the second part,we introduce a network simplex method that directly solves a network piecewise linear program without reformulating it as a network linear program of higher dimension. Forty benchmark NPLP problems are solved by this method and a reformulation method. The computational results are in favor of the direct method and show that solving an NPLP problem is not much harder than solving a network linear program of the same dimension.

An MPI+OpenACC-Based PRM Scalar Advection Scheme in the GRAPES Model over a Cluster with Multiple CPUs and GPUs

A moisture advection scheme is an essential module of a numerical weather/climate model representing the horizontal transport of water vapor.The Piecewise Rational Method(PRM) scalar advection scheme in the Global/Regional Assimilation and Prediction System(GRAPES) solves the moisture flux advection equation based on PRM.Computation of the scalar advection involves boundary exchange,and computation of higher bandwidth requirements is complicated and time-consuming in GRAPES.Recently,Graphics Processing Units(GPUs) have been widely used to solve scientific and engineering computing problems owing to advancements in GPU hardware and related programming models such as CUDA/OpenCL and Open Accelerator(OpenACC).Herein,we present an accelerated PRM scalar advection scheme with Message Passing Interface(MPI) and OpenACC to fully exploit GPUs’ power over a cluster with multiple Central Processing Units(CPUs) and GPUs,together with optimization of various parameters such as minimizing data transfer,memory coalescing,exposing more parallelism,and overlapping computation with data transfers.Results show that about 3.5 times speedup is obtained for the entire model running at medium resolution with double precision when comparing the scheme’s elapsed time on a node with two GPUs(NVIDIA P100) and two 16-core CPUs(Intel Gold 6142).Further,results obtained from experiments of a higher resolution model with multiple GPUs show excellent scalability.

Stabilization of Nonuniform Euler-Bernoulli Beam with Locally Distributed Feedbacks

In this article, we study the stabilization problem of a nonuniform Euler-Bernoulli beam with locally distributed feedbacks. Firstly, using the semi-group theory, we establish the well-posedness of the associated closed loop system. Then by proving the uniqueness of the solution of a related ordinary differential equations, we derive the asymptotic stability of the closed loop system. Finally, by means of the piecewise frequency domain multiplier method, we prove that the corresponding closed loop system can be exponentially stabilized by only one of the two distributed feedback controls proposed in this paper.

Numerical study of air cavity characteristics of bow wave breaking of KCS ship under different speeds

Ship bow wave breaking contains complex flow mechanism,which is very important for ship performance.In this study,a practical numerical simulation scheme for bow wave breaking is proposed and the scheme is applied to the simulation of bow wave breaking of KCS ship model with Fr=0.26,0.30,0.35,0.40,analyzing the impact of speed on the bow wave breaking.The results indicate that an increase in speed leads to a significant rise in viscous pressure resistance and more pronounced bow wave breaking.Moreover,it is found that the traditional wave height function in OpenFOAM is not suitable for detailed studies of bow wave breaking.This study extracts different free surfaces through top and bottom views to further analyze the free surface overturning,droplet splashing,and cavity entrainment in bow wave breaking.Additionally,the spatial and temporal distribution of cavities at Fr=0.40 is analyzed,revealing that cavity distribution is closely related to vortex structures and exhibits a periodic pulsation characteristic of approximately 12 s.