It is a significant task to predict the solar activity for space weather and solar physics. All kinds of approaches have been used to forecast solar activities, and they have been applied to many areas such as the sol...It is a significant task to predict the solar activity for space weather and solar physics. All kinds of approaches have been used to forecast solar activities, and they have been applied to many areas such as the solar dynamo of simulation and space mission planning. In this paper, we employ the long-shortterm memory(LSTM) and neural network autoregression(NNAR) deep learning methods to predict the upcoming 25 th solar cycle using the sunspot area(SSA) data during the period of May 1874 to December2020. Our results show that the 25 th solar cycle will be 55% stronger than Solar Cycle 24 with a maximum sunspot area of 3115±401 and the cycle reaching its peak in October 2022 by using the LSTM method. It also shows that deep learning algorithms perform better than the other commonly used methods and have high application value.展开更多
In this paper, a fully third-order accurate projection method for solving the incompressible Navier-Stokes equations is proposed. To construct the scheme, a continuous projection procedure is firstly presented. We the...In this paper, a fully third-order accurate projection method for solving the incompressible Navier-Stokes equations is proposed. To construct the scheme, a continuous projection procedure is firstly presented. We then derive a sufficient condition for the continuous projection equations to be temporally third-order accurate approximations of the original Navier-Stokes equations by means of the localtruncation-error-analysis technique. The continuous projection equations are discretized temporally and spatially to third-order accuracy on the staggered grids, resulting in a fully third-order discrete projection scheme. The possibility to design higher-order projection methods is thus demonstrated in the present paper. A heuristic stability analysis is performed on this projection method showing the probability of its being stable. The stability of the present scheme is further verified through numerical tests. The third-order accuracy of the present projection method is validated by several numerical test cases.展开更多
Generalized Jacobi polynomials with indexes α,β∈ R are introduced and some basic properties are established. As examples of applications,the second- and fourth-order elliptic boundary value problems with Dirichlet ...Generalized Jacobi polynomials with indexes α,β∈ R are introduced and some basic properties are established. As examples of applications,the second- and fourth-order elliptic boundary value problems with Dirichlet or Robin boundary conditions are considered,and the generalized Jacobi spectral schemes are proposed. For the diagonalization of discrete systems,the Jacobi-Sobolev orthogonal basis functions are constructed,which allow the exact solutions and the approximate solutions to be represented in the forms of infinite and truncated Jacobi series. Error estimates are obtained and numerical results are provided to illustrate the effectiveness and the spectral accuracy.展开更多
A hyperbolic Lindstedt-Poincare method is presented to determine the homoclinic solutions of a kind of nonlinear oscillators, in which critical value of the homoclinic bifurcation parameter can be determined. The gene...A hyperbolic Lindstedt-Poincare method is presented to determine the homoclinic solutions of a kind of nonlinear oscillators, in which critical value of the homoclinic bifurcation parameter can be determined. The generalized Lienard oscillator is studied in detail, and the present method＇s predictions are compared with those of Runge-Kutta method to illustrate its accuracy.展开更多
With the development of large-scale spectral surveys, fiber positioning technology has been developing rapidly. Because of the performance advantages of a four-quadrant（4Q） detector, a fiber positioning and real-tim...With the development of large-scale spectral surveys, fiber positioning technology has been developing rapidly. Because of the performance advantages of a four-quadrant（4Q） detector, a fiber positioning and real-time monitoring system based on the 4Q detector is proposed. The detection accuracy of this system is directly determined by the precision of the center of the spot. A Gaussian fitting algorithm based on the 4Q detector is studied and applied in the fiber positioning process to improve the calculated accuracy of the spot center. The relationship between the center position of the incident spot and the detector output signal is deduced. An experimental platform is built to complete the simulated experiment. Then we use the Gaussian fitting method to process experimental data, compare the fitting value with the theoretical one and calculate the corresponding error.展开更多
The three-way catalysts （TWCs） promoters Ce0.6Zr0.4- x TbxO2-y were prepared by sol-gel method. BET surface areas analysis indicated that an increase of the dopant Tb content from x = 0.05 to x = 0.15 favors an incr...The three-way catalysts （TWCs） promoters Ce0.6Zr0.4- x TbxO2-y were prepared by sol-gel method. BET surface areas analysis indicated that an increase of the dopant Tb content from x = 0.05 to x = 0.15 favors an increase of surface area from 66.8 to 80.4 m^2· g^-1 compared with the undoped sample Ce0 .6oZr0.40O2 65.1 m^2·g^- 1 after calcination at 650℃. Transmission electron microscopy （TEM） observation indicated that the doped samples have a higher thermal stability. The XRD and Raman spectra confirmed that the Ce0.6Zr0.4-xTbxO2-y cubic solid solution is formed. XPS analysis revealed that Ce and Tb mainly existed in the form of Ce^4＋ and Tb^3 ＋ , and Zr existed in the form of Zr^4＋ on the surface of the samples. The doped samples were homogenous in composition ; the introduction of Tb into the CeO2-ZrO2 promoters resuited in the formation of a solid solution, and the concentration of surface lattice oxygen was increased.展开更多
In general, triangular and quadrilateral elements are commonly applied in two-dimensional finite element methods. If they are used to compute polycrystalline materials, the cost of computation can be quite significant...In general, triangular and quadrilateral elements are commonly applied in two-dimensional finite element methods. If they are used to compute polycrystalline materials, the cost of computation can be quite significant. Polygonal elements can do well in simulation of the materials behavior and provide greater flexibility for the meshing of complex geometries. Hence, the study on the polygonal element is a very useful and necessary part in the finite element method. In this paper, an n-sided polygonal element based on quadratic spline interpolant, denoted by PS2 element, is presented using the triangular area coordinates and the B-net method. The PS2 element is conforming and can exactly model the quadratic field. It is valid for both convex and non-convex polygonal element, and insensitive to mesh distortions. In addition, no mapping or coordinate transformation is required and thus no Jacobian matrix and its inverse are evaluated. Some appropriate examples are employed to evaluate the performance of the proposed element.展开更多
Multi-train modeling and simulation plays a vital role in railway electrification during operation and planning phase. Study of peak power demand and energy consumed by each traction substation needs to be deter- mine...Multi-train modeling and simulation plays a vital role in railway electrification during operation and planning phase. Study of peak power demand and energy consumed by each traction substation needs to be deter- mined to verify that electrical energy flowing in its railway power feeding system is appropriate or not. Gauss-Seidel, conventional Newton-Raphson, and current injection methods are well-known and widely accepted as a tool for electrical power network solver in DC railway power supply study. In this paper, a simplified Newton-Raphson method has been proposed. The proposed method employs a set of current-balance equations at each electrical node instead of the conventional power-balance equation used in the conventional Newton-Raphson method. This concept can remarkably reduce execution time and computing complexity for multi-train simulation. To evaluate its use, Sukhumvit line of Bangkok transit system （BTS） of Thai- land with 21.6-km line length and 22 passenger stopping stations is set as a test system. The multi-train simulation integrated with the proposed power network solver is developed to simulate 1-h operation service of selected 5-min headway. From the obtained results, the proposed method is more efficient with approximately 18 % faster than the conventional Newton-Raphson method and just over 6 % faster than the current injection method.展开更多
This paper proposes a hybrid vertex-centered fi- nite volume/finite element method for solution of the two di- mensional （2D） incompressible Navier-Stokes equations on unstructured grids. An incremental pressure fra...This paper proposes a hybrid vertex-centered fi- nite volume/finite element method for solution of the two di- mensional （2D） incompressible Navier-Stokes equations on unstructured grids. An incremental pressure fractional step method is adopted to handle the velocity-pressure coupling. The velocity and the pressure are collocated at the node of the vertex-centered control volume which is formed by join- ing the centroid of cells sharing the common vertex. For the temporal integration of the momentum equations, an im- plicit second-order scheme is utilized to enhance the com- putational stability and eliminate the time step limit due to the diffusion term. The momentum equations are discretized by the vertex-centered finite volume method （FVM） and the pressure Poisson equation is solved by the Galerkin finite el- ement method （FEM）. The momentum interpolation is used to damp out the spurious pressure wiggles. The test case with analytical solutions demonstrates second-order accuracy of the current hybrid scheme in time and space for both veloc- ity and pressure. The classic test cases, the lid-driven cavity flow, the skew cavity flow and the backward-facing step flow, show that numerical results are in good agreement with the published benchmark solutions.展开更多
Condensation technique of degree of freedom is first proposed to improve the computational efficiency of meshfree method with Galerkin weak form for elastic dy- namic analysis. In the present method, scattered nodes w...Condensation technique of degree of freedom is first proposed to improve the computational efficiency of meshfree method with Galerkin weak form for elastic dy- namic analysis. In the present method, scattered nodes with- out connectivity are divided into several subsets by cells with arbitrary shape. Local discrete equation is established over each cell by using moving Kriging interpolation, in which the nodes that located in the cell are used for approxima- tion. Then local discrete equations can be simplified by con- densation of degree of freedom, which transfers equations of inner nodes to equations of boundary nodes based on cells. The global dynamic system equations are obtained by as- sembling all local discrete equations and are solved by using the standard implicit Newmark＇s time integration scheme. In the scheme of present method, the calculation of each cell is carried out by meshfree method, and local search is imple- mented in interpolation. Numerical examples show that the present method has high computational efficiency and good accuracy in solving elastic dynamic problems.展开更多
A multi-model integration method is proposed to develop a multi-source and heterogeneous model for short-term solar flare prediction. Different prediction models are constructed on the basis of extracted predictors fr...A multi-model integration method is proposed to develop a multi-source and heterogeneous model for short-term solar flare prediction. Different prediction models are constructed on the basis of extracted predictors from a pool of observation databases. The outputs of the base models are normal- ized first because these established models extract predictors from many data resources using different prediction methods. Then weighted integration of the base models is used to develop a multi-model integrated model （MIM）. The weight set that single models assign is optimized by a genetic algorithm. Seven base models and data from Solar and Heliospheric Observatory/Michelson Doppler Imager lon- gitudinal magnetograms are used to construct the MIM, and then its performance is evaluated by cross validation. Experimental results showed that the MIM outperforms any individual model in nearly every data group, and the richer the diversity of the base models, the better the performance of the MIM. Thus, integrating more diversified models, such as an expert system, a statistical model and a physical model, will greatly improve the performance of the MIM.展开更多
Isoparametric quadrilateral elements are widely used in the finite element method, but the accuracy of the isoparametric quadrilateral elements will drop obviously deteriorate due to mesh distortions. Spline functions...Isoparametric quadrilateral elements are widely used in the finite element method, but the accuracy of the isoparametric quadrilateral elements will drop obviously deteriorate due to mesh distortions. Spline functions have some properties of simplicity and conformality. Two 8-node quadrilateral elements have been developed using the trian- gular area coordinates and the B-net method, which can ex- actly model the quadratic field for both convex and concave quadrangles. Some appropriate examples are employed to evaluate the performance of the proposed elements. The nu- merical results show that the two spline elements can obtain solutions which are highly accurate and insensitive to mesh distortions.展开更多
According to earthquake catalog records of Fujian Seismic Network, the Tnow method and the fourstation continuous location method put forward by Jin Xing are inspected by using P-wave arrival information of the first ...According to earthquake catalog records of Fujian Seismic Network, the Tnow method and the fourstation continuous location method put forward by Jin Xing are inspected by using P-wave arrival information of the first four stations in each earthquake. It shows that the fourstation continuous location method can locate more seismic events than the Tnow method. By analyzing the results, it is concluded that the reason for this is that the Tnow method makes use of information from stations without being triggered, while some stations failed to be reflected in earthquake catalog because of discontinuous records or unclear records of seismic phases. For seismic events whose location results can be given, there is no obvious difference in location results of the two methods and positioning deviation of most seismic events is also not significant. For earthquakes outside the network, the positioning deviation may amplify as the epicentral distance enlarges, which may relate to the situation that the seismic stations are centered on one side of epicenter and the opening angle between seismic stations used for location and epicenter is small.展开更多
A precursor of BaCe0.5Zr0.4Y0. 1O3-α electrolytes was synthesized by the sol-gel method and sintered at temperature which were 150 - 250 ℃ lower than by solid state reaction. The AC impedance spectrums of electrolyt...A precursor of BaCe0.5Zr0.4Y0. 1O3-α electrolytes was synthesized by the sol-gel method and sintered at temperature which were 150 - 250 ℃ lower than by solid state reaction. The AC impedance spectrums of electrolytes were measured by AUTOLAB PGSTA30 electrochemical measuring device at different temperatures. The conductivities of the electrolytes are 1.62×10^-4 - 6.43×10^-3, 2.52×10^-5 - 3.73×10^-3S·cm^-1 in the temperature range of 350-800℃. The activity energies are 0.54 and 0.84 eV. At the same time BaCe0.9Y0.1O3-α was prepared by direct solid state reaction. The conductivity of BaCe0.9Y0.1O3-α is 1 × 10^-4- 4×10^-3 S·cm^-1 and the activation energy is 0.50 eV at the same condition. The results show that conduction of electrolyte prepared by sol-gel method is higher than the one by solid state reaction. As far as BaCe0.9Y0.1O3-α concerned, its conductivity of the Zr-substituted specimens is decreased.展开更多
We use wavelet transform to analyze the daily relative sunspot number series over solar cycles 10-23. The characteristics of some of the periods shorter than - 600-day are discussed. The results exhibit not only the v...We use wavelet transform to analyze the daily relative sunspot number series over solar cycles 10-23. The characteristics of some of the periods shorter than - 600-day are discussed. The results exhibit not only the variation of some short periods in the 14 solar cycles but also the characteristics and differences around solar peaks and valley years. The short periodic components with larger amplitude such as ～27, ～ 150 and ～360-day are obvious in some solar cycles, all of them are time-variable, also their lengths and amplitudes are variable and intermittent in time. The variable characteristics of the periods are rather different in different solar cycles.展开更多
Sunspots are the most striking and easily observed magnetic structures of the Sun,and statistical analysis of solar historical data could reveal a wealth of information on the long-term variation of solar activity cyc...Sunspots are the most striking and easily observed magnetic structures of the Sun,and statistical analysis of solar historical data could reveal a wealth of information on the long-term variation of solar activity cycle.The hand-drawn sunspot records of Yunnan Observatories,Chinese Academy of Sciences have been accumulating for more than 60 years,and nearly 16000 images have been preserved.In the future,the observation mode of recording sunspots by hand-drawing will be replaced inevitably by digital images observed either at ground or in space.To connect the hand-drawn sunspot data and the purely digital sunspot data in future,it is necessary to analyze the systematic errors of the data which are observed by the two observation modes in the period of transition.In this paper,we choose 268 round sunspots(Htype in modified Zurich sunspot classification)from the drawing of Yunnan Observatories to compare their positions and areas with the CCD observations made by Helioseismic and Magnetic Imager(HMI)on board Solar Dynamic Observatory(SDO)and Global Oscillation Network Group(GONG).We find that the latitude and longitude accuracy of hand-drawn sunspot are within-0.127 and 2.29 degree respectively,and the area accuracy is about 16.36 sunspot unit(μHem).Systematic errors apparently decrease with large sunspot.展开更多
Three-and five-minute oscillations are commonly observed in any sunspot.Because they are modulated by the internal thermal and magnetic structures of a sunspot,they could be used as an effective tool for researching s...Three-and five-minute oscillations are commonly observed in any sunspot.Because they are modulated by the internal thermal and magnetic structures of a sunspot,they could be used as an effective tool for researching sunspot seismology.In this paper,we investigate the properties of oscillations in sunspot groups with varying sizes and magnetic fields,and aim to establish the relationships between sunspot oscillations and its internal structure comparatively.We selected three groups of the unipolar sunspot with approximately axial-symmetric magnetic field and calculated their Fourier spectra based on the ultraviolet/extreme ultraviolet emission intensity variations recorded by the Solar Dynamics Observatory/Atmospheric Imaging Assembly.We found that the distribution of three-minute oscillation is defined by the joint effect of diverging magnetic field and the stratification of the sunspot atmosphere.Its distribution could be modified by any invading magnetic structures in the umbra.In contrast,the five-minute oscillations are more prominent in small spots,implying that five-minute oscillation is very closely connected with umbral dynamics.展开更多
The five-minute oscillations inside sunspots appear to be the absorption of the solar p-mode. It is a potential tool to probe a sunspot＇s sub-structure. We studied the collective properties of five-minute oscillation...The five-minute oscillations inside sunspots appear to be the absorption of the solar p-mode. It is a potential tool to probe a sunspot＇s sub-structure. We studied the collective properties of five-minute oscillations in the power and phase distribution at the sunspot＇s umbra-penumbra boundary. The azimuthal distributions of the power and phase of five-minute oscillations enclosing a sunspot＇s umbra were obtained with images taken with the Solar Dynamics Observatory/Atmospheric Imaging Assembly （SDO/AIA）. The azimuthal modes were quantified with periodogram analysis and justified with significance tests. The azimuthal nodal structures in an approximately ax- ially symmetric sunspot AR 11131 （2010 Dec 08） were investigated. Mode numbers ra = 2, 3, 4, 7, 10 were obtained in both 1700 A and 1600A bandpasses. The 1600A channel also revealed an extra mode at m = 9. In the upper atmosphere （304 A）, fewer modes were detected at m = 3, 4, 7. The azimuthal modes in the sunspot＇s low atmo- sphere could be interpreted as high-order azimuthal MHD body modes. They were detected in the power and phase of the five-minute oscillations in sunspot AR 11131 with SDO/AIA data. Fewer modes were detected in the sunspot＇s upper atmosphere.展开更多
As concepts closely related to microwave absorption properties,impedance matching and phase matching were rarely combined with material parameters to regulate properties and explore related mechanisms.In this work,red...As concepts closely related to microwave absorption properties,impedance matching and phase matching were rarely combined with material parameters to regulate properties and explore related mechanisms.In this work,reduction–diffusion method was innovatively applied to synthesize rare earth alloy Y_(2)Fe_(17).In order to regulate the electromagnetic parameters of absorbers,the Y_(2)Fe_(17)N_(3-δ)particles were coated with silica(Y_(2)Fe_(17)N_(3-δ)@SiO_(2))and absorbers with different volume fractions were prepared.The relationship between impedance matching,matching thickness,and the strongest reflection loss peak(RLmin)was presented obviously.Compared to the microwave absorption properties of Y_(2)Fe_(17)N_(3-δ)/PU absorber,Y_(2)Fe_(17)N_(3-δ)@SiO_(2)/PU absorbers are more conducive to the realization of microwave absorption material standards which are thin thickness,light weight,strong absorbing intensity,and broad bandwidth.Based on microwave frequency bands,the microwave absorption properties of the absorbers were analyzed and the related parameters were listed.As an important parameter related to perfect matching,reflection factor(√ε_(r)/μ_(r))was discussed combined with microwave amplitude attenuation.According to the origin and mathematical model of bandwidth,the formula of EAB(RL<-10 dB)was derived and simplified.The calculated bandwidths agreed well with experimental results.展开更多
基金supported by the National Natural Science Foundation of China under Grant numbers U2031202,U1731124 and U1531247the special foundation work of the Ministry of Science and Technology of the People’s Republic of China under Grant number 2014FY120300the 13th Five-year Informatization Plan of Chinese Academy of Sciences under Grant number XXH13505-04。
文摘It is a significant task to predict the solar activity for space weather and solar physics. All kinds of approaches have been used to forecast solar activities, and they have been applied to many areas such as the solar dynamo of simulation and space mission planning. In this paper, we employ the long-shortterm memory(LSTM) and neural network autoregression(NNAR) deep learning methods to predict the upcoming 25 th solar cycle using the sunspot area(SSA) data during the period of May 1874 to December2020. Our results show that the 25 th solar cycle will be 55% stronger than Solar Cycle 24 with a maximum sunspot area of 3115±401 and the cycle reaching its peak in October 2022 by using the LSTM method. It also shows that deep learning algorithms perform better than the other commonly used methods and have high application value.
基金The project supported by the China NKBRSF(2001CB409604)
文摘In this paper, a fully third-order accurate projection method for solving the incompressible Navier-Stokes equations is proposed. To construct the scheme, a continuous projection procedure is firstly presented. We then derive a sufficient condition for the continuous projection equations to be temporally third-order accurate approximations of the original Navier-Stokes equations by means of the localtruncation-error-analysis technique. The continuous projection equations are discretized temporally and spatially to third-order accuracy on the staggered grids, resulting in a fully third-order discrete projection scheme. The possibility to design higher-order projection methods is thus demonstrated in the present paper. A heuristic stability analysis is performed on this projection method showing the probability of its being stable. The stability of the present scheme is further verified through numerical tests. The third-order accuracy of the present projection method is validated by several numerical test cases.
基金the National Natural Science Foundation of China (Nos.11571238,11601332,91130014,11471312 and 91430216).
文摘Generalized Jacobi polynomials with indexes α,β∈ R are introduced and some basic properties are established. As examples of applications,the second- and fourth-order elliptic boundary value problems with Dirichlet or Robin boundary conditions are considered,and the generalized Jacobi spectral schemes are proposed. For the diagonalization of discrete systems,the Jacobi-Sobolev orthogonal basis functions are constructed,which allow the exact solutions and the approximate solutions to be represented in the forms of infinite and truncated Jacobi series. Error estimates are obtained and numerical results are provided to illustrate the effectiveness and the spectral accuracy.
基金supported by the National Natural Science Foundation of China (10672193)Sun Yat-sen University (Fu Lan Scholarship)the University of Hong Kong (CRGC grant).
文摘A hyperbolic Lindstedt-Poincare method is presented to determine the homoclinic solutions of a kind of nonlinear oscillators, in which critical value of the homoclinic bifurcation parameter can be determined. The generalized Lienard oscillator is studied in detail, and the present method＇s predictions are compared with those of Runge-Kutta method to illustrate its accuracy.
基金support by the Fundamental Research Funds for the Central Universities of China (2013/B15020271)the National Natural Science Foundation of China (1014/515029111)the National Undergraduate Training Program for Innovation and Entrepreneurship (201610294069)
文摘With the development of large-scale spectral surveys, fiber positioning technology has been developing rapidly. Because of the performance advantages of a four-quadrant（4Q） detector, a fiber positioning and real-time monitoring system based on the 4Q detector is proposed. The detection accuracy of this system is directly determined by the precision of the center of the spot. A Gaussian fitting algorithm based on the 4Q detector is studied and applied in the fiber positioning process to improve the calculated accuracy of the spot center. The relationship between the center position of the incident spot and the detector output signal is deduced. An experimental platform is built to complete the simulated experiment. Then we use the Gaussian fitting method to process experimental data, compare the fitting value with the theoretical one and calculate the corresponding error.
基金Project supported by the Natural Science Basic Research Plan in Shaanxi Province of China (2004B13)
文摘The three-way catalysts （TWCs） promoters Ce0.6Zr0.4- x TbxO2-y were prepared by sol-gel method. BET surface areas analysis indicated that an increase of the dopant Tb content from x = 0.05 to x = 0.15 favors an increase of surface area from 66.8 to 80.4 m^2· g^-1 compared with the undoped sample Ce0 .6oZr0.40O2 65.1 m^2·g^- 1 after calcination at 650℃. Transmission electron microscopy （TEM） observation indicated that the doped samples have a higher thermal stability. The XRD and Raman spectra confirmed that the Ce0.6Zr0.4-xTbxO2-y cubic solid solution is formed. XPS analysis revealed that Ce and Tb mainly existed in the form of Ce^4＋ and Tb^3 ＋ , and Zr existed in the form of Zr^4＋ on the surface of the samples. The doped samples were homogenous in composition ; the introduction of Tb into the CeO2-ZrO2 promoters resuited in the formation of a solid solution, and the concentration of surface lattice oxygen was increased.
基金supported by the National Natural Science Foundation of China (60533060, 10672032, 10726067)Science Foundation of Dalian University of Technology (SFDUT07001)
文摘In general, triangular and quadrilateral elements are commonly applied in two-dimensional finite element methods. If they are used to compute polycrystalline materials, the cost of computation can be quite significant. Polygonal elements can do well in simulation of the materials behavior and provide greater flexibility for the meshing of complex geometries. Hence, the study on the polygonal element is a very useful and necessary part in the finite element method. In this paper, an n-sided polygonal element based on quadratic spline interpolant, denoted by PS2 element, is presented using the triangular area coordinates and the B-net method. The PS2 element is conforming and can exactly model the quadratic field. It is valid for both convex and non-convex polygonal element, and insensitive to mesh distortions. In addition, no mapping or coordinate transformation is required and thus no Jacobian matrix and its inverse are evaluated. Some appropriate examples are employed to evaluate the performance of the proposed element.
文摘Multi-train modeling and simulation plays a vital role in railway electrification during operation and planning phase. Study of peak power demand and energy consumed by each traction substation needs to be deter- mined to verify that electrical energy flowing in its railway power feeding system is appropriate or not. Gauss-Seidel, conventional Newton-Raphson, and current injection methods are well-known and widely accepted as a tool for electrical power network solver in DC railway power supply study. In this paper, a simplified Newton-Raphson method has been proposed. The proposed method employs a set of current-balance equations at each electrical node instead of the conventional power-balance equation used in the conventional Newton-Raphson method. This concept can remarkably reduce execution time and computing complexity for multi-train simulation. To evaluate its use, Sukhumvit line of Bangkok transit system （BTS） of Thai- land with 21.6-km line length and 22 passenger stopping stations is set as a test system. The multi-train simulation integrated with the proposed power network solver is developed to simulate 1-h operation service of selected 5-min headway. From the obtained results, the proposed method is more efficient with approximately 18 % faster than the conventional Newton-Raphson method and just over 6 % faster than the current injection method.
基金supported by the Natural Science Foundation of China (11061021)the Program of Higher-level talents of Inner Mongolia University (SPH-IMU,Z200901004)the Scientific Research Projection of Higher Schools of Inner Mongolia(NJ10016,NJ10006)
文摘This paper proposes a hybrid vertex-centered fi- nite volume/finite element method for solution of the two di- mensional （2D） incompressible Navier-Stokes equations on unstructured grids. An incremental pressure fractional step method is adopted to handle the velocity-pressure coupling. The velocity and the pressure are collocated at the node of the vertex-centered control volume which is formed by join- ing the centroid of cells sharing the common vertex. For the temporal integration of the momentum equations, an im- plicit second-order scheme is utilized to enhance the com- putational stability and eliminate the time step limit due to the diffusion term. The momentum equations are discretized by the vertex-centered finite volume method （FVM） and the pressure Poisson equation is solved by the Galerkin finite el- ement method （FEM）. The momentum interpolation is used to damp out the spurious pressure wiggles. The test case with analytical solutions demonstrates second-order accuracy of the current hybrid scheme in time and space for both veloc- ity and pressure. The classic test cases, the lid-driven cavity flow, the skew cavity flow and the backward-facing step flow, show that numerical results are in good agreement with the published benchmark solutions.
基金supported by the National Natural Science Founda-tion of China(11272118)Open Found of State Key Laboratory of Explosion Science and Technology(KFJJ12-5M)
文摘Condensation technique of degree of freedom is first proposed to improve the computational efficiency of meshfree method with Galerkin weak form for elastic dy- namic analysis. In the present method, scattered nodes with- out connectivity are divided into several subsets by cells with arbitrary shape. Local discrete equation is established over each cell by using moving Kriging interpolation, in which the nodes that located in the cell are used for approxima- tion. Then local discrete equations can be simplified by con- densation of degree of freedom, which transfers equations of inner nodes to equations of boundary nodes based on cells. The global dynamic system equations are obtained by as- sembling all local discrete equations and are solved by using the standard implicit Newmark＇s time integration scheme. In the scheme of present method, the calculation of each cell is carried out by meshfree method, and local search is imple- mented in interpolation. Numerical examples show that the present method has high computational efficiency and good accuracy in solving elastic dynamic problems.
基金supported by the National Natural Science Foundation of China(Grant No.11078010)SOHO is a project of international cooperation between the European Space Agency(ESA) and NASA
文摘A multi-model integration method is proposed to develop a multi-source and heterogeneous model for short-term solar flare prediction. Different prediction models are constructed on the basis of extracted predictors from a pool of observation databases. The outputs of the base models are normal- ized first because these established models extract predictors from many data resources using different prediction methods. Then weighted integration of the base models is used to develop a multi-model integrated model （MIM）. The weight set that single models assign is optimized by a genetic algorithm. Seven base models and data from Solar and Heliospheric Observatory/Michelson Doppler Imager lon- gitudinal magnetograms are used to construct the MIM, and then its performance is evaluated by cross validation. Experimental results showed that the MIM outperforms any individual model in nearly every data group, and the richer the diversity of the base models, the better the performance of the MIM. Thus, integrating more diversified models, such as an expert system, a statistical model and a physical model, will greatly improve the performance of the MIM.
基金supported by the National Natural Science Foundation of China(11001037,11102037 and 11290143)the Fundamental Research Funds for the Central Universities
文摘Isoparametric quadrilateral elements are widely used in the finite element method, but the accuracy of the isoparametric quadrilateral elements will drop obviously deteriorate due to mesh distortions. Spline functions have some properties of simplicity and conformality. Two 8-node quadrilateral elements have been developed using the trian- gular area coordinates and the B-net method, which can ex- actly model the quadratic field for both convex and concave quadrangles. Some appropriate examples are employed to evaluate the performance of the proposed elements. The nu- merical results show that the two spline elements can obtain solutions which are highly accurate and insensitive to mesh distortions.
文摘According to earthquake catalog records of Fujian Seismic Network, the Tnow method and the fourstation continuous location method put forward by Jin Xing are inspected by using P-wave arrival information of the first four stations in each earthquake. It shows that the fourstation continuous location method can locate more seismic events than the Tnow method. By analyzing the results, it is concluded that the reason for this is that the Tnow method makes use of information from stations without being triggered, while some stations failed to be reflected in earthquake catalog because of discontinuous records or unclear records of seismic phases. For seismic events whose location results can be given, there is no obvious difference in location results of the two methods and positioning deviation of most seismic events is also not significant. For earthquakes outside the network, the positioning deviation may amplify as the epicentral distance enlarges, which may relate to the situation that the seismic stations are centered on one side of epicenter and the opening angle between seismic stations used for location and epicenter is small.
文摘A precursor of BaCe0.5Zr0.4Y0. 1O3-α electrolytes was synthesized by the sol-gel method and sintered at temperature which were 150 - 250 ℃ lower than by solid state reaction. The AC impedance spectrums of electrolytes were measured by AUTOLAB PGSTA30 electrochemical measuring device at different temperatures. The conductivities of the electrolytes are 1.62×10^-4 - 6.43×10^-3, 2.52×10^-5 - 3.73×10^-3S·cm^-1 in the temperature range of 350-800℃. The activity energies are 0.54 and 0.84 eV. At the same time BaCe0.9Y0.1O3-α was prepared by direct solid state reaction. The conductivity of BaCe0.9Y0.1O3-α is 1 × 10^-4- 4×10^-3 S·cm^-1 and the activation energy is 0.50 eV at the same condition. The results show that conduction of electrolyte prepared by sol-gel method is higher than the one by solid state reaction. As far as BaCe0.9Y0.1O3-α concerned, its conductivity of the Zr-substituted specimens is decreased.
基金Supported by the National Natural Science Foundation of China.
文摘We use wavelet transform to analyze the daily relative sunspot number series over solar cycles 10-23. The characteristics of some of the periods shorter than - 600-day are discussed. The results exhibit not only the variation of some short periods in the 14 solar cycles but also the characteristics and differences around solar peaks and valley years. The short periodic components with larger amplitude such as ～27, ～ 150 and ～360-day are obvious in some solar cycles, all of them are time-variable, also their lengths and amplitudes are variable and intermittent in time. The variable characteristics of the periods are rather different in different solar cycles.
基金supported by the National Natural Science Foundation of China(Grant Nos.U1731124,U1531247,11427901 and 11873089)the special foundation work of the Ministry of Science and Technology of China(Grant No.2014FY120300)+1 种基金the 13th Five-year Informatization Plan of Chinese Academy of Sciences(Grant No.XXH13505–04)the Youth Innovation Promotion Association CAS.The hand-drawing historic。
文摘Sunspots are the most striking and easily observed magnetic structures of the Sun,and statistical analysis of solar historical data could reveal a wealth of information on the long-term variation of solar activity cycle.The hand-drawn sunspot records of Yunnan Observatories,Chinese Academy of Sciences have been accumulating for more than 60 years,and nearly 16000 images have been preserved.In the future,the observation mode of recording sunspots by hand-drawing will be replaced inevitably by digital images observed either at ground or in space.To connect the hand-drawn sunspot data and the purely digital sunspot data in future,it is necessary to analyze the systematic errors of the data which are observed by the two observation modes in the period of transition.In this paper,we choose 268 round sunspots(Htype in modified Zurich sunspot classification)from the drawing of Yunnan Observatories to compare their positions and areas with the CCD observations made by Helioseismic and Magnetic Imager(HMI)on board Solar Dynamic Observatory(SDO)and Global Oscillation Network Group(GONG).We find that the latitude and longitude accuracy of hand-drawn sunspot are within-0.127 and 2.29 degree respectively,and the area accuracy is about 16.36 sunspot unit(μHem).Systematic errors apparently decrease with large sunspot.
基金supported by the National Natural Science Foundation of China(NSFC,Grant Nos.12173012 and 12111530078)the Shenzhen Technology project(GXWD20201230155427003-20200804151658001)supported by the Joint Funds of the NSFC(U1931107)。
文摘Three-and five-minute oscillations are commonly observed in any sunspot.Because they are modulated by the internal thermal and magnetic structures of a sunspot,they could be used as an effective tool for researching sunspot seismology.In this paper,we investigate the properties of oscillations in sunspot groups with varying sizes and magnetic fields,and aim to establish the relationships between sunspot oscillations and its internal structure comparatively.We selected three groups of the unipolar sunspot with approximately axial-symmetric magnetic field and calculated their Fourier spectra based on the ultraviolet/extreme ultraviolet emission intensity variations recorded by the Solar Dynamics Observatory/Atmospheric Imaging Assembly.We found that the distribution of three-minute oscillation is defined by the joint effect of diverging magnetic field and the stratification of the sunspot atmosphere.Its distribution could be modified by any invading magnetic structures in the umbra.In contrast,the five-minute oscillations are more prominent in small spots,implying that five-minute oscillation is very closely connected with umbral dynamics.
文摘The five-minute oscillations inside sunspots appear to be the absorption of the solar p-mode. It is a potential tool to probe a sunspot＇s sub-structure. We studied the collective properties of five-minute oscillations in the power and phase distribution at the sunspot＇s umbra-penumbra boundary. The azimuthal distributions of the power and phase of five-minute oscillations enclosing a sunspot＇s umbra were obtained with images taken with the Solar Dynamics Observatory/Atmospheric Imaging Assembly （SDO/AIA）. The azimuthal modes were quantified with periodogram analysis and justified with significance tests. The azimuthal nodal structures in an approximately ax- ially symmetric sunspot AR 11131 （2010 Dec 08） were investigated. Mode numbers ra = 2, 3, 4, 7, 10 were obtained in both 1700 A and 1600A bandpasses. The 1600A channel also revealed an extra mode at m = 9. In the upper atmosphere （304 A）, fewer modes were detected at m = 3, 4, 7. The azimuthal modes in the sunspot＇s low atmo- sphere could be interpreted as high-order azimuthal MHD body modes. They were detected in the power and phase of the five-minute oscillations in sunspot AR 11131 with SDO/AIA data. Fewer modes were detected in the sunspot＇s upper atmosphere.
基金Project supported by the National Key Research and Development Program of China(Grant No.2021YFB3501300)the National Natural Science Foundation of China(Grant No.51731001)the Fund from the State Key Laboratory of Baiyunobo Rare Earth Resource Researches and Comprehensive Utilization’s Key Research and Development Projects。
文摘As concepts closely related to microwave absorption properties,impedance matching and phase matching were rarely combined with material parameters to regulate properties and explore related mechanisms.In this work,reduction–diffusion method was innovatively applied to synthesize rare earth alloy Y_(2)Fe_(17).In order to regulate the electromagnetic parameters of absorbers,the Y_(2)Fe_(17)N_(3-δ)particles were coated with silica(Y_(2)Fe_(17)N_(3-δ)@SiO_(2))and absorbers with different volume fractions were prepared.The relationship between impedance matching,matching thickness,and the strongest reflection loss peak(RLmin)was presented obviously.Compared to the microwave absorption properties of Y_(2)Fe_(17)N_(3-δ)/PU absorber,Y_(2)Fe_(17)N_(3-δ)@SiO_(2)/PU absorbers are more conducive to the realization of microwave absorption material standards which are thin thickness,light weight,strong absorbing intensity,and broad bandwidth.Based on microwave frequency bands,the microwave absorption properties of the absorbers were analyzed and the related parameters were listed.As an important parameter related to perfect matching,reflection factor(√ε_(r)/μ_(r))was discussed combined with microwave amplitude attenuation.According to the origin and mathematical model of bandwidth,the formula of EAB(RL<-10 dB)was derived and simplified.The calculated bandwidths agreed well with experimental results.
