This paper presents a simple approach for improving the performance of the weighted essentially nonoscillatory(WENO) finite volume scheme on non-uniform grids. This technique relies on the reformulation of the fifthor...This paper presents a simple approach for improving the performance of the weighted essentially nonoscillatory(WENO) finite volume scheme on non-uniform grids. This technique relies on the reformulation of the fifthorder WENO-JS(WENO scheme presented by Jiang and Shu in J. Comput. Phys. 126:202–228, 1995) scheme designed on uniform grids in terms of one cell-averaged value and its left and/or right interfacial values of the dependent variable.The effect of grid non-uniformity is taken into consideration by a proper interpolation of the interfacial values. On nonuniform grids, the proposed scheme is much more accurate than the original WENO-JS scheme, which was designed for uniform grids. When the grid is uniform, the resulting scheme reduces to the original WENO-JS scheme. In the meantime,the proposed scheme is computationally much more efficient than the fifth-order WENO scheme designed specifically for the non-uniform grids. A number of numerical test cases are simulated to verify the performance of the present scheme.展开更多
The large-scale development of wind power is an important means to reduce greenhouse gas emissions, alleviate environmental pollution and improve the utilization rate of renewable energy. At the same time, large-scale...The large-scale development of wind power is an important means to reduce greenhouse gas emissions, alleviate environmental pollution and improve the utilization rate of renewable energy. At the same time, large-scale non grid connected wind power generation theory avoids the technical difficulties of wind power integration [1]. However, due to the randomness and uncontrollability of wind energy, the output power of the wind power generation system will fluctuate accordingly [2]. Therefore, the corresponding energy storage devices are arranged in the non-grid-connected wind power generation system to ensure the power quality, and it has become the key to full utilization of renewable energy. In the case of wind speed fluctuation, the DC bus control strategy of the wind turbine is proposed in this paper. It can reduce the impact on the unit converter and the power load;this ensures safe and stable operation of non-grid connected wind turbines.展开更多
A class of nonlinear parabolic equation on a polygonal domain Ω R2 is inves- tigated in this paper. We introduce a finite element method on overlapping non-matching grids for the nonlinear parabolic equation based o...A class of nonlinear parabolic equation on a polygonal domain Ω R2 is inves- tigated in this paper. We introduce a finite element method on overlapping non-matching grids for the nonlinear parabolic equation based on the partition of unity method. We give the construction and convergence analysis for the semi-discrete and the fully discrete finite element methods. Moreover, we prove that the error of the discrete variational problem has good approximation properties. Our results are valid for any spatial dimensions. A numerical example to illustrate the theoretical results is also given.展开更多
This paper presents a finite-difference(FD)method with spatially non-rectangular irregular grids to simulate the elastic wave propagation.Staggered irregular grid finite difference oper- ators with a second-order time...This paper presents a finite-difference(FD)method with spatially non-rectangular irregular grids to simulate the elastic wave propagation.Staggered irregular grid finite difference oper- ators with a second-order time and spatial accuracy are used to approximate the velocity-stress elastic wave equations.This method is very simple and the cost of computing time is not much.Complicated geometries like curved thin layers,cased borehole and nonplanar interfaces may be treated with non- rectangular irregular grids in a more flexible way.Unlike the multi-grid scheme,this method requires no interpolation between the fine and coarse grids and all grids are computed at the same spatial iteration.Compared with the rectangular irregular grid FD,the spurious diffractions from'staircase' interfaces can easily be eliminated without using finer grids.Dispersion and stability conditions of the proposed method can be established in a similar form as for the rectangular irregular grid scheme.The Higdon's absorbing boundary condition is adopted to eliminate boundary reflections.Numerical simu- lations show that this method has satisfactory stability and accuracy in simulating wave propagation near rough solid-fluid interfaces.The computation costs are less than those using a regular grid and rectangular grid FD method.展开更多
Finite difference methods have been widely employed in solving the eikonal equation so as to calculate traveltime of seismic phase. Most previous studies used regular orthogonal grid. However, much denser grid is requ...Finite difference methods have been widely employed in solving the eikonal equation so as to calculate traveltime of seismic phase. Most previous studies used regular orthogonal grid. However, much denser grid is required to sample the interfaces that are undulating in depth direction, such as the Moho and the 660 km discontinuity.Here we propose a new finite difference algorithm to solve the eikonal equation on non-orthogonal grid(irregular grid).To demonstrate its efficiency and accuracy, a test was conducted with a two-layer model. The test result suggests that the similar accuracy of a regular grid with ten times grids could achieve with our new algorithm, but the time cost is only about 0.1 times. A spherical earth model with an undulant660 km discontinuity was constructed to demonstrate the potential application of our new method. In that case, the traveltime curve fluctuation corresponds to topography. Our new algorithm is efficient in solving the first arrival times of waves associated with undulant interfaces.展开更多
In this paper a three degrees of freedom autoparametric system with limited power supply is investigated numerically. The system consists of the body, which is hung on a spring and a damper, and two pendulums connecte...In this paper a three degrees of freedom autoparametric system with limited power supply is investigated numerically. The system consists of the body, which is hung on a spring and a damper, and two pendulums connected by shape memory alloy (SMA) spring. Shape memory alloys have ability to change their material properties with temperature. A polynomial constitutive model is assumed to describe the behavior of the SMA spring. The non-ideal source of power adds one degree of freedom, so the system has four degrees of freedom. The equations of motion have been solved numerically and pseudoelastic effects associated with the martensitic phase transformation are studied. It is shown that in this type system one mode of vibrations might excite or damp another mode, and that except different kinds of periodic vibrations there may also appear chaotic vibrations. For the identification of the responses of the system's various techniques, including chaos techniques such as bifurcation diagrams and time histories, power spectral densities, Poincare maps and exponents of Lyapunov may be used.展开更多
Let g be a complex simple Lie algebra of rank ι, b the standard Borel subalgebra. An invertible map on Ь is said to preserve abelian ideals if it maps each abelian ideal to some such ideal of the same dimension. In ...Let g be a complex simple Lie algebra of rank ι, b the standard Borel subalgebra. An invertible map on Ь is said to preserve abelian ideals if it maps each abelian ideal to some such ideal of the same dimension. In this article, by using some results of Chevalley groups, the theory of root systems and root space decomposition, the author gives an explicit description on such maps of Ь.展开更多
A new scale transformation method is used in solving the Schrodinger equation. With it, the uniform grids in the discretization in conventional metho d are changed into non-uniform grids. Consequently, in some cases, ...A new scale transformation method is used in solving the Schrodinger equation. With it, the uniform grids in the discretization in conventional metho d are changed into non-uniform grids. Consequently, in some cases, the computing quantity will be greatly reduced at keeping the required accuracy. The calcul ation of the quantized inversion layer in MOS structure is used to demonstrate t he efficiency of the new method.展开更多
An expanding model of the confinement of non-ideal detonation of small charge is established on the basis of the nozzle theory.Making use of the expanding model,the analytic relationship of small charge detonation vel...An expanding model of the confinement of non-ideal detonation of small charge is established on the basis of the nozzle theory.Making use of the expanding model,the analytic relationship of small charge detonation velocity and the semi-empirical relationship of detonation pressure that both change with charge diameter and confinement condition are established.The detonation velocity and pressure of small charges are calculated and experimentally verified,and the detonation velocity deviation is less than 7% while the detonation pressure deviation is less than 9%.展开更多
In present paper, certain aspect of shock wave in non-ideal gas, when magnetic field is orthogonal to the trajectories of the gas particles and electrical conductivity is taken to be infinite, is investigated. Conside...In present paper, certain aspect of shock wave in non-ideal gas, when magnetic field is orthogonal to the trajectories of the gas particles and electrical conductivity is taken to be infinite, is investigated. Considering one-dimensional unsteady non-planer motion, basic equations, its general solution and formation of shock-wave, conservation laws and jumps conditions, variation of area of non-uniform cross section and analytical solution of strong non planer shock is obtained.展开更多
In the present work the widths of layers constituting the non-ideal superlattice are much bigger then the characteristic scales of space dispersion. In such a case the contribution of individual layers to gyrotropy ca...In the present work the widths of layers constituting the non-ideal superlattice are much bigger then the characteristic scales of space dispersion. In such a case the contribution of individual layers to gyrotropy can be regarded as independed. Thus the corresponding optical quantities can be expressed through the layers’ gyrotropic characteristics. This approach is applied to calculate the specific rotation angle of plane of polarization of light propagating through a nonideal 1D-superlattice, which varies in composition as well as in layers’ width. We carry out numerical calculation of the frequency dispersion of optical activity of a non-ideal superlattice, which includes impurity layers with point defects.展开更多
基金supported by the National Natural Science Foundation of China (Grant 11672160)the National Key Research and Development Program of China (Grant 2016YF A0401200)
文摘This paper presents a simple approach for improving the performance of the weighted essentially nonoscillatory(WENO) finite volume scheme on non-uniform grids. This technique relies on the reformulation of the fifthorder WENO-JS(WENO scheme presented by Jiang and Shu in J. Comput. Phys. 126:202–228, 1995) scheme designed on uniform grids in terms of one cell-averaged value and its left and/or right interfacial values of the dependent variable.The effect of grid non-uniformity is taken into consideration by a proper interpolation of the interfacial values. On nonuniform grids, the proposed scheme is much more accurate than the original WENO-JS scheme, which was designed for uniform grids. When the grid is uniform, the resulting scheme reduces to the original WENO-JS scheme. In the meantime,the proposed scheme is computationally much more efficient than the fifth-order WENO scheme designed specifically for the non-uniform grids. A number of numerical test cases are simulated to verify the performance of the present scheme.
文摘The large-scale development of wind power is an important means to reduce greenhouse gas emissions, alleviate environmental pollution and improve the utilization rate of renewable energy. At the same time, large-scale non grid connected wind power generation theory avoids the technical difficulties of wind power integration [1]. However, due to the randomness and uncontrollability of wind energy, the output power of the wind power generation system will fluctuate accordingly [2]. Therefore, the corresponding energy storage devices are arranged in the non-grid-connected wind power generation system to ensure the power quality, and it has become the key to full utilization of renewable energy. In the case of wind speed fluctuation, the DC bus control strategy of the wind turbine is proposed in this paper. It can reduce the impact on the unit converter and the power load;this ensures safe and stable operation of non-grid connected wind turbines.
基金Supported by the Natural Science Foundation of Hunan under Grant No. 06C713.
文摘A class of nonlinear parabolic equation on a polygonal domain Ω R2 is inves- tigated in this paper. We introduce a finite element method on overlapping non-matching grids for the nonlinear parabolic equation based on the partition of unity method. We give the construction and convergence analysis for the semi-discrete and the fully discrete finite element methods. Moreover, we prove that the error of the discrete variational problem has good approximation properties. Our results are valid for any spatial dimensions. A numerical example to illustrate the theoretical results is also given.
文摘This paper presents a finite-difference(FD)method with spatially non-rectangular irregular grids to simulate the elastic wave propagation.Staggered irregular grid finite difference oper- ators with a second-order time and spatial accuracy are used to approximate the velocity-stress elastic wave equations.This method is very simple and the cost of computing time is not much.Complicated geometries like curved thin layers,cased borehole and nonplanar interfaces may be treated with non- rectangular irregular grids in a more flexible way.Unlike the multi-grid scheme,this method requires no interpolation between the fine and coarse grids and all grids are computed at the same spatial iteration.Compared with the rectangular irregular grid FD,the spurious diffractions from'staircase' interfaces can easily be eliminated without using finer grids.Dispersion and stability conditions of the proposed method can be established in a similar form as for the rectangular irregular grid scheme.The Higdon's absorbing boundary condition is adopted to eliminate boundary reflections.Numerical simu- lations show that this method has satisfactory stability and accuracy in simulating wave propagation near rough solid-fluid interfaces.The computation costs are less than those using a regular grid and rectangular grid FD method.
基金supported in part by National Basic Research Program of China (No. 2014CB845900)Hubei Provincial Natural Science Foundation of China (No. 2014CFA005)
文摘Finite difference methods have been widely employed in solving the eikonal equation so as to calculate traveltime of seismic phase. Most previous studies used regular orthogonal grid. However, much denser grid is required to sample the interfaces that are undulating in depth direction, such as the Moho and the 660 km discontinuity.Here we propose a new finite difference algorithm to solve the eikonal equation on non-orthogonal grid(irregular grid).To demonstrate its efficiency and accuracy, a test was conducted with a two-layer model. The test result suggests that the similar accuracy of a regular grid with ten times grids could achieve with our new algorithm, but the time cost is only about 0.1 times. A spherical earth model with an undulant660 km discontinuity was constructed to demonstrate the potential application of our new method. In that case, the traveltime curve fluctuation corresponds to topography. Our new algorithm is efficient in solving the first arrival times of waves associated with undulant interfaces.
文摘In this paper a three degrees of freedom autoparametric system with limited power supply is investigated numerically. The system consists of the body, which is hung on a spring and a damper, and two pendulums connected by shape memory alloy (SMA) spring. Shape memory alloys have ability to change their material properties with temperature. A polynomial constitutive model is assumed to describe the behavior of the SMA spring. The non-ideal source of power adds one degree of freedom, so the system has four degrees of freedom. The equations of motion have been solved numerically and pseudoelastic effects associated with the martensitic phase transformation are studied. It is shown that in this type system one mode of vibrations might excite or damp another mode, and that except different kinds of periodic vibrations there may also appear chaotic vibrations. For the identification of the responses of the system's various techniques, including chaos techniques such as bifurcation diagrams and time histories, power spectral densities, Poincare maps and exponents of Lyapunov may be used.
基金Supported by the Doctor Foundation of Henan Polytechnic University(B2010-93)Supported by the National Natural Science Foundation of China(11126121)+2 种基金Supported by the Natural Science Foundation of Henan Province(112300410120)Supported by the Natural Science Research Program of Education Department of Henan Province(201lB110016)Supported by the Applied Mathematics Provincial-level Key Discipline of Henan Province of Henau Polytechuic University
文摘Let g be a complex simple Lie algebra of rank ι, b the standard Borel subalgebra. An invertible map on Ь is said to preserve abelian ideals if it maps each abelian ideal to some such ideal of the same dimension. In this article, by using some results of Chevalley groups, the theory of root systems and root space decomposition, the author gives an explicit description on such maps of Ь.
文摘A new scale transformation method is used in solving the Schrodinger equation. With it, the uniform grids in the discretization in conventional metho d are changed into non-uniform grids. Consequently, in some cases, the computing quantity will be greatly reduced at keeping the required accuracy. The calcul ation of the quantized inversion layer in MOS structure is used to demonstrate t he efficiency of the new method.
文摘An expanding model of the confinement of non-ideal detonation of small charge is established on the basis of the nozzle theory.Making use of the expanding model,the analytic relationship of small charge detonation velocity and the semi-empirical relationship of detonation pressure that both change with charge diameter and confinement condition are established.The detonation velocity and pressure of small charges are calculated and experimentally verified,and the detonation velocity deviation is less than 7% while the detonation pressure deviation is less than 9%.
文摘In present paper, certain aspect of shock wave in non-ideal gas, when magnetic field is orthogonal to the trajectories of the gas particles and electrical conductivity is taken to be infinite, is investigated. Considering one-dimensional unsteady non-planer motion, basic equations, its general solution and formation of shock-wave, conservation laws and jumps conditions, variation of area of non-uniform cross section and analytical solution of strong non planer shock is obtained.
文摘In the present work the widths of layers constituting the non-ideal superlattice are much bigger then the characteristic scales of space dispersion. In such a case the contribution of individual layers to gyrotropy can be regarded as independed. Thus the corresponding optical quantities can be expressed through the layers’ gyrotropic characteristics. This approach is applied to calculate the specific rotation angle of plane of polarization of light propagating through a nonideal 1D-superlattice, which varies in composition as well as in layers’ width. We carry out numerical calculation of the frequency dispersion of optical activity of a non-ideal superlattice, which includes impurity layers with point defects.