Linear scan computed tomography （LCT） is of great benefit to online industrial scanning and security inspection due to its characteristics of straight-line source trajectory and high scanning speed. However, in prac...Linear scan computed tomography （LCT） is of great benefit to online industrial scanning and security inspection due to its characteristics of straight-line source trajectory and high scanning speed. However, in practical applications of LCT, there are challenges to image reconstruction due to limited-angle and insufficient data. In this paper, a new reconstruction algorithm based on total-variation （TV） minimization is developed to reconstruct images from limited-angle and insufficient data in LCT. The main idea of our approach is to reformulate a TV problem as a linear equality constrained problem where the objective function is separable, and then minimize its augmented Lagrangian function by using alternating direction method （ADM） to solve subproblems. The proposed method is robust and efficient in the task of reconstruction by showing the convergence of ADM. The numerical simulations and real data reconstructions show that the proposed reconstruction method brings reasonable performance and outperforms some previous ones when applied to an LCT imaging problem.展开更多
The task of dividing corrupted-data into their respective subspaces can be well illustrated,both theoretically and numerically,by recovering low-rank and sparse-column components of a given matrix.Generally,it can be ...The task of dividing corrupted-data into their respective subspaces can be well illustrated,both theoretically and numerically,by recovering low-rank and sparse-column components of a given matrix.Generally,it can be characterized as a matrix and a 2,1-norm involved convex minimization problem.However,solving the resulting problem is full of challenges due to the non-smoothness of the objective function.One of the earliest solvers is an 3-block alternating direction method of multipliers(ADMM)which updates each variable in a Gauss-Seidel manner.In this paper,we present three variants of ADMM for the 3-block separable minimization problem.More preciously,whenever one variable is derived,the resulting problems can be regarded as a convex minimization with 2 blocks,and can be solved immediately using the standard ADMM.If the inner iteration loops only once,the iterative scheme reduces to the ADMM with updates in a Gauss-Seidel manner.If the solution from the inner iteration is assumed to be exact,the convergence can be deduced easily in the literature.The performance comparisons with a couple of recently designed solvers illustrate that the proposed methods are effective and competitive.展开更多
The principle of direct method used in optimal control problem is introduced. Details of applying this method to flight trajectory generation are presented including calculation of velocity and controls histories. And...The principle of direct method used in optimal control problem is introduced. Details of applying this method to flight trajectory generation are presented including calculation of velocity and controls histories. And capabilities of flight and propulsion systems are considered also. Combined with digital terrain map technique, the direct method is applied to the three dimensional trajectory optimization for low altitude penetration, and simplex algorithm is used to solve the parameters in optimization. For the small number of parameters, the trajectory can be optimized in real time on board.展开更多
Since the connection of small-scale wind farms to distribution networks,power grid voltage stability has been reduced with increasing wind penetration in recent years,owing to the variable reactive power consumption o...Since the connection of small-scale wind farms to distribution networks,power grid voltage stability has been reduced with increasing wind penetration in recent years,owing to the variable reactive power consumption of wind generators.In this study,a two-stage reactive power optimization method based on the alternating direction method of multipliers(ADMM)algorithm is proposed for achieving optimal reactive power dispatch in wind farm-integrated distribution systems.Unlike existing optimal reactive power control methods,the proposed method enables distributed reactive power flow optimization with a two-stage optimization structure.Furthermore,under the partition concept,the consensus protocol is not needed to solve the optimization problems.In this method,the influence of the wake effect of each wind turbine is also considered in the control design.Simulation results for a mid-voltage distribution system based on MATLAB verified the effectiveness of the proposed method.展开更多
This paper applies the singular integral operators, singular quadrature operators and discretization matrices associated with singular integral equations with Cauchy kernels, which are established in [1], to give a un...This paper applies the singular integral operators, singular quadrature operators and discretization matrices associated with singular integral equations with Cauchy kernels, which are established in [1], to give a unified framework for various collocation methods of numerical solutions of singular integral equations with Cauchy kernels. Under the framework, the coincidence of the direct quadrature method and the indirect quadrature method is very simple and obvious.展开更多
This paper presents a design of optimal controllers with respect to a meaningful cost function to force an underactuated omni-directional intelligent navigator （ODIN） under unknown constant environmental loads to tr...This paper presents a design of optimal controllers with respect to a meaningful cost function to force an underactuated omni-directional intelligent navigator （ODIN） under unknown constant environmental loads to track a reference trajectory in two-dimensional space. Motivated by the vehicle＇s steering practice, the yaw angle regarded as a virtual control plus the surge thrust force are used to force the position of the vehicle to globally track its reference trajectory. The control design is based on several recent results developed for inverse optimal control and stability analysis of nonlinear systems, a new design of bounded disturbance observers, and backstepping and Lyapunov＇s direct methods. Both state- and output-feedback control designs are addressed. Simulations are included to illustrate the effectiveness of the proposed results.展开更多
This paper determines the exact error order on optimization of adaptive direct methods of approximate solution of the class of Fredholm integral equations of the second kind with kernel belonging to the anisotropic So...This paper determines the exact error order on optimization of adaptive direct methods of approximate solution of the class of Fredholm integral equations of the second kind with kernel belonging to the anisotropic Sobolev classes, and also gives an optimal algorithm.展开更多
In this study, we use the direct discontinuous Galerkin method to solve the generalized Burgers-Fisher equation. The method is based on the direct weak formulation of the Burgers-Fisher equation. The two adjacent cell...In this study, we use the direct discontinuous Galerkin method to solve the generalized Burgers-Fisher equation. The method is based on the direct weak formulation of the Burgers-Fisher equation. The two adjacent cells are jointed by a numerical flux that includes the convection numerical flux and the diffusion numerical flux. We solve the ordinary differential equations arising in the direct Galerkin method by using the strong stability preserving Runge^Kutta method. Numerical results are compared with the exact solution and the other results to show the accuracy and reliability of the method.展开更多
A direct method to find the first integral for two-dimensional autonomous system in polar coordinates is suggested. It is shown that if the equation of motion expressed by differential 1-forms for a given autonomous H...A direct method to find the first integral for two-dimensional autonomous system in polar coordinates is suggested. It is shown that if the equation of motion expressed by differential 1-forms for a given autonomous Hamiltonian system is multiplied by a set of multiplicative functions, then the general expression of the first integral can be obtained, An example is given to illustrate the application of the results.展开更多
This paper studies the generalized Kawahara equation in terms of the approximate homotopy symmetry method and the approximate homotopy direct method. Using both methods it obtains the similarity reduction solutions an...This paper studies the generalized Kawahara equation in terms of the approximate homotopy symmetry method and the approximate homotopy direct method. Using both methods it obtains the similarity reduction solutions and similarity reduction equations of different orders, showing that the approximate homotopy direct method yields more general approximate similarity reductions than the approximate homotopy symmetry method. The homotopy series solutions to the generalized Kawahara equation are consequently derived.展开更多
Silicoaluminophosphates SAPO-11 molecular sieves with small particle size and hierarchical pores were synthesized using the directing agent method.The effect of crystallization time on the particle structure,morpholog...Silicoaluminophosphates SAPO-11 molecular sieves with small particle size and hierarchical pores were synthesized using the directing agent method.The effect of crystallization time on the particle structure,morphology,pore structure properties,and acid properties of SAPO-11 molecular sieves were investigated.Unlike the SAPO-11 molecular sieves synthesized with the conventional method,the results of XRD,SEM,BET and NH3-TPD analyses showed that the SAPO-11 molecular sieves synthesized by the directing agent method in a shorter crystallization time exhibited fine and uniform morphology.By increasing the crystallization time,the particle size of SAPO-11 molecular sieve was significantly reduced,and the mesoporous structure(intercrystalline pores)was formed.Furthermore,the external specific surface area and the total specific surface area reached 81.7 m^2/g and 192.0 m^2/g,respectively,which effectively reduced the pore mass transfer resistance and significantly increased the number of acid sites.The results of n-dodecane hydroisomerization revealed that the Pt/SAPO-11 prepared with the novel method exhibited higher catalytic activity and better hydroisomerization selectivity than that synthesized by the conventional hydrothermal method.Thus,the small particle molecular sieve showed a promising industrial application prospect to be used as catalyst support.展开更多
Gas holdup is one of the key parameters in flotation process. Gas holdup as measured by a differential pressure method was investigated and the relative errors compared to the average gas holdup from the volume expans...Gas holdup is one of the key parameters in flotation process. Gas holdup as measured by a differential pressure method was investigated and the relative errors compared to the average gas holdup from the volume expansion method. The errors were used to establish optimum measurement positions. The results show that the measurement position should be in the middle of the column and in the region half way from the center to the wall (the half-radius). The gas holdup along the axial direction is lower at the bottom and higher at the top of the floatation column. The gas holdup along the radial direction is lower near the wall and higher near the center of the flotation column. The average gas holdup measure- ment can be replaced by regional gas holdup values.展开更多
One of the essential points of the direct-method single-wavelength anomalous diffraction （SAD） phasing for proteins is to express the bimodal SAD phase distribution by the sum of two Gaussian functions peaked respec...One of the essential points of the direct-method single-wavelength anomalous diffraction （SAD） phasing for proteins is to express the bimodal SAD phase distribution by the sum of two Gaussian functions peaked respectively at φh″＋｜△φh｜ and φh″-｜△φh｜. The probability for △φh being positive （P＋） can be derived based on the Cochran distribution in direct methods. Hence the SAD phase ambiguity can be resolved by multiplying the Gaussian function peaked at φh″＋｜△φh｜ with P＋ and multiplying the Gaussian function peaked at φh″-｜△φh｜ with P_ （=1- P＋）. The direct-method SAD h phasing has been proved powerful in breaking SAD phase ambiguities, in particular when anomalous-scattering signals are weak. However, the approximation of bimodal phase distributions by the sum of two Gaussian functions introduces considerable errors. In this paper we show that a much better approximation can be achieved by replacing the two Gaussian functions with two von Mises distributions. Test results showed that this leads to significant improvement on the efficiency of direct-method SAD-phasing.展开更多
In this article, by introducing characteristic singular integral operator and associate singular integral equations （SIEs）, the authors discuss the direct method of solution for a class of singular integral equation...In this article, by introducing characteristic singular integral operator and associate singular integral equations （SIEs）, the authors discuss the direct method of solution for a class of singular integral equations with certain analytic inputs. They obtain both the conditions of solvability and the solutions in closed form. It is noteworthy that the method is different from the classical one that is due to Lu.展开更多
This paper proposes a hybrid decoupled power flow method for balanced power distribution systems with distributed generation sources. The method formulates the power flow equations in active power and reactive power d...This paper proposes a hybrid decoupled power flow method for balanced power distribution systems with distributed generation sources. The method formulates the power flow equations in active power and reactive power decoupled form with polar coordinates. Second-order terms are included in the active power mismatch iteration, and constant Jacobian and Hessian matrices are used. A hybrid direct and indirect solution technique is used to achieve efficiency and robustness of the algorithm. Active power correction is solved by means of a sparse lower triangular and upper triangular （LU） decomposition algorithm with partial pivoting, and the reactive power correction is solved by means of restarted generalized minimal residual algorithm with an incomplete LU pre-conditioner. Typical distribution generation models and distribution load models are included. The impact of zero-impedance branches is explicitly modeled through reconfiguring of the adjacent branches with impedances. Numerical examples on a sample distribution system with widespread photovoltaic installations are given to demonstrate the effectiveness of the proposed method.展开更多
So far, Lou＇s direct perturbation method has been applied successfully to solve the nonlinear SchrSdinger equation（NLSE） hierarchy, such as the NLSE, the coupled NLSE, the critical NLSE, and the derivative NLSE. Bu...So far, Lou＇s direct perturbation method has been applied successfully to solve the nonlinear SchrSdinger equation（NLSE） hierarchy, such as the NLSE, the coupled NLSE, the critical NLSE, and the derivative NLSE. But to our knowledge, this method for other types of perturbed nonlinear evolution equations has still been lacking. In this paper, Lou＇s direct perturbation method is applied to the study of perturbed complex Burgers equation. By this method, we calculate not only the zero-order adiabatic solution, but also the first order modification.展开更多
The general interpolation mentioned in this a rticle provides an effective way for reducing the amount of calculation of direc t optimal exploration. It has been testified by real case calculations that the interpolat...The general interpolation mentioned in this a rticle provides an effective way for reducing the amount of calculation of direc t optimal exploration. It has been testified by real case calculations that the interpolation is not only reliable but also can save the amount of calculation by nearly 36%. Large amount of calculation and lacking strict theoretical bas is has been th e two disadvantage of direct method by new. If this defect is not overcome, they will not only s eriously affect the application of this method, but also hinder its further rese arch. Based on sufficient calculation practice, this article has made a primary discussion about the theory and method of reducing the amount of calculation, an d has achieved some satisfactory results.展开更多
The approximate direct reduction method is applied to the perturbed mKdV equation with weak fourth order dispersion and weak dissipation. The similarity reduction solutions of different orders conform to formal cohere...The approximate direct reduction method is applied to the perturbed mKdV equation with weak fourth order dispersion and weak dissipation. The similarity reduction solutions of different orders conform to formal coherence, accounting for infinite series reduction solutions to the original equation and general formulas of similarity reduction equations. Painleve Ⅱ type equations, hyperbolic secant and Jacobi elliptic function solutions are obtained for zeroorder similarity reduction equations. Higher order similarity reduction equations are linear variable coefficient ordinary differential equations.展开更多
In this work, we present the direct discontinuous Galerkin （DDG） method for the one-dimensional coupled nonlinear Schr5dinger （CNLS） equation. We prove that the new discontinuous Galerkin method preserves the disc...In this work, we present the direct discontinuous Galerkin （DDG） method for the one-dimensional coupled nonlinear Schr5dinger （CNLS） equation. We prove that the new discontinuous Galerkin method preserves the discrete mass conservations corresponding to the properties of the CNLS system. The ordinary differential equations obtained by the DDG space discretization is solved via a third-order stabilized Runge Kutta method. Numerical experiments show that the new DDG scheme gives stable and less diffusive results and has excellent long-time numerical behaviors for the CNLS equations.展开更多
Crack propagation path in two-directionally graded composites was investigated by the finite element method.A graded extended finite element method(XFEM)was employed to calculate displacement and stress fields in crac...Crack propagation path in two-directionally graded composites was investigated by the finite element method.A graded extended finite element method(XFEM)was employed to calculate displacement and stress fields in cracked graded structures.And a post-processing subroutine of interaction energy integral was implemented to extract the mixed-mode stress intensity factors(SIFs).The maximum hoop stress(MHS)criterion was adopted to predict crack growth direction based on the assumption of local homogenization of asymptotic crack-tip fields in graded materials.Effects of material nonhomogeneous parameters on crack propagation paths were also discussed in detail.It is shown that the present method can provide relatively accurate predictions of crack paths in two-directionally graded composites.Crack propagates in the decreasing direction of effective Young′s modulus.The shape and steepness of property gradient perpendicular to the crack surface have great influences on crack paths.Through redesigning material property reasonably,crack growth in graded material can be changed to improve mechanical behaviours of cracked structures.展开更多
基金the National High Technology Research and Development Program of China(Grant No.2012AA011603)
文摘Linear scan computed tomography （LCT） is of great benefit to online industrial scanning and security inspection due to its characteristics of straight-line source trajectory and high scanning speed. However, in practical applications of LCT, there are challenges to image reconstruction due to limited-angle and insufficient data. In this paper, a new reconstruction algorithm based on total-variation （TV） minimization is developed to reconstruct images from limited-angle and insufficient data in LCT. The main idea of our approach is to reformulate a TV problem as a linear equality constrained problem where the objective function is separable, and then minimize its augmented Lagrangian function by using alternating direction method （ADM） to solve subproblems. The proposed method is robust and efficient in the task of reconstruction by showing the convergence of ADM. The numerical simulations and real data reconstructions show that the proposed reconstruction method brings reasonable performance and outperforms some previous ones when applied to an LCT imaging problem.
基金Supported by the National Natural Science Foundation of China(Grant No.11971149,11871381)Natural Science Foundation of Henan Province for Youth(Grant No.202300410146)。
文摘The task of dividing corrupted-data into their respective subspaces can be well illustrated,both theoretically and numerically,by recovering low-rank and sparse-column components of a given matrix.Generally,it can be characterized as a matrix and a 2,1-norm involved convex minimization problem.However,solving the resulting problem is full of challenges due to the non-smoothness of the objective function.One of the earliest solvers is an 3-block alternating direction method of multipliers(ADMM)which updates each variable in a Gauss-Seidel manner.In this paper,we present three variants of ADMM for the 3-block separable minimization problem.More preciously,whenever one variable is derived,the resulting problems can be regarded as a convex minimization with 2 blocks,and can be solved immediately using the standard ADMM.If the inner iteration loops only once,the iterative scheme reduces to the ADMM with updates in a Gauss-Seidel manner.If the solution from the inner iteration is assumed to be exact,the convergence can be deduced easily in the literature.The performance comparisons with a couple of recently designed solvers illustrate that the proposed methods are effective and competitive.
文摘The principle of direct method used in optimal control problem is introduced. Details of applying this method to flight trajectory generation are presented including calculation of velocity and controls histories. And capabilities of flight and propulsion systems are considered also. Combined with digital terrain map technique, the direct method is applied to the three dimensional trajectory optimization for low altitude penetration, and simplex algorithm is used to solve the parameters in optimization. For the small number of parameters, the trajectory can be optimized in real time on board.
基金support of The National Key Research and Development Program of China(Basic Research Class)(No.2017YFB0903000)the National Natural Science Foundation of China(No.U1909201)。
文摘Since the connection of small-scale wind farms to distribution networks,power grid voltage stability has been reduced with increasing wind penetration in recent years,owing to the variable reactive power consumption of wind generators.In this study,a two-stage reactive power optimization method based on the alternating direction method of multipliers(ADMM)algorithm is proposed for achieving optimal reactive power dispatch in wind farm-integrated distribution systems.Unlike existing optimal reactive power control methods,the proposed method enables distributed reactive power flow optimization with a two-stage optimization structure.Furthermore,under the partition concept,the consensus protocol is not needed to solve the optimization problems.In this method,the influence of the wake effect of each wind turbine is also considered in the control design.Simulation results for a mid-voltage distribution system based on MATLAB verified the effectiveness of the proposed method.
基金NNSF of China, SF of SEC of China and SF of Wuhan University.
文摘This paper applies the singular integral operators, singular quadrature operators and discretization matrices associated with singular integral equations with Cauchy kernels, which are established in [1], to give a unified framework for various collocation methods of numerical solutions of singular integral equations with Cauchy kernels. Under the framework, the coincidence of the direct quadrature method and the indirect quadrature method is very simple and obvious.
基金Supported in Part by the Australian Research Council under Grant DP0988424
文摘This paper presents a design of optimal controllers with respect to a meaningful cost function to force an underactuated omni-directional intelligent navigator （ODIN） under unknown constant environmental loads to track a reference trajectory in two-dimensional space. Motivated by the vehicle＇s steering practice, the yaw angle regarded as a virtual control plus the surge thrust force are used to force the position of the vehicle to globally track its reference trajectory. The control design is based on several recent results developed for inverse optimal control and stability analysis of nonlinear systems, a new design of bounded disturbance observers, and backstepping and Lyapunov＇s direct methods. Both state- and output-feedback control designs are addressed. Simulations are included to illustrate the effectiveness of the proposed results.
基金Project supported by the Natural Science Foundation of China(10371009)Research Fund for the Doctoral Program Higher Education
文摘This paper determines the exact error order on optimization of adaptive direct methods of approximate solution of the class of Fredholm integral equations of the second kind with kernel belonging to the anisotropic Sobolev classes, and also gives an optimal algorithm.
基金Project supported by the National Natural Science Foundation of China (Grant Nos. 61105130 and 61175124)
文摘In this study, we use the direct discontinuous Galerkin method to solve the generalized Burgers-Fisher equation. The method is based on the direct weak formulation of the Burgers-Fisher equation. The two adjacent cells are jointed by a numerical flux that includes the convection numerical flux and the diffusion numerical flux. We solve the ordinary differential equations arising in the direct Galerkin method by using the strong stability preserving Runge^Kutta method. Numerical results are compared with the exact solution and the other results to show the accuracy and reliability of the method.
文摘A direct method to find the first integral for two-dimensional autonomous system in polar coordinates is suggested. It is shown that if the equation of motion expressed by differential 1-forms for a given autonomous Hamiltonian system is multiplied by a set of multiplicative functions, then the general expression of the first integral can be obtained, An example is given to illustrate the application of the results.
基金Project supported by the National Natural Science Foundations of China(Grant Nos.10735030,10475055,10675065 and 90503006)the National Basic Research Program of China(Grant No.2007CB814800)
文摘This paper studies the generalized Kawahara equation in terms of the approximate homotopy symmetry method and the approximate homotopy direct method. Using both methods it obtains the similarity reduction solutions and similarity reduction equations of different orders, showing that the approximate homotopy direct method yields more general approximate similarity reductions than the approximate homotopy symmetry method. The homotopy series solutions to the generalized Kawahara equation are consequently derived.
基金We thank the National Natural Science Fund of China(2016-Z0030)Natural Science Foundation of Liaoning Province(L2017 LQN008,L2019014).
文摘Silicoaluminophosphates SAPO-11 molecular sieves with small particle size and hierarchical pores were synthesized using the directing agent method.The effect of crystallization time on the particle structure,morphology,pore structure properties,and acid properties of SAPO-11 molecular sieves were investigated.Unlike the SAPO-11 molecular sieves synthesized with the conventional method,the results of XRD,SEM,BET and NH3-TPD analyses showed that the SAPO-11 molecular sieves synthesized by the directing agent method in a shorter crystallization time exhibited fine and uniform morphology.By increasing the crystallization time,the particle size of SAPO-11 molecular sieve was significantly reduced,and the mesoporous structure(intercrystalline pores)was formed.Furthermore,the external specific surface area and the total specific surface area reached 81.7 m^2/g and 192.0 m^2/g,respectively,which effectively reduced the pore mass transfer resistance and significantly increased the number of acid sites.The results of n-dodecane hydroisomerization revealed that the Pt/SAPO-11 prepared with the novel method exhibited higher catalytic activity and better hydroisomerization selectivity than that synthesized by the conventional hydrothermal method.Thus,the small particle molecular sieve showed a promising industrial application prospect to be used as catalyst support.
基金supports for this work provided by the NationalKey Technology R&D Program in the 11th Five-Year Plan of China(No. 2008BAB31B03)
文摘Gas holdup is one of the key parameters in flotation process. Gas holdup as measured by a differential pressure method was investigated and the relative errors compared to the average gas holdup from the volume expansion method. The errors were used to establish optimum measurement positions. The results show that the measurement position should be in the middle of the column and in the region half way from the center to the wall (the half-radius). The gas holdup along the axial direction is lower at the bottom and higher at the top of the floatation column. The gas holdup along the radial direction is lower near the wall and higher near the center of the flotation column. The average gas holdup measure- ment can be replaced by regional gas holdup values.
基金Project supported by the Innovation Foundation of the Chinese Academy of Sciences and by the National Basic Research Program of China(Grant No.2002CB713801)
文摘One of the essential points of the direct-method single-wavelength anomalous diffraction （SAD） phasing for proteins is to express the bimodal SAD phase distribution by the sum of two Gaussian functions peaked respectively at φh″＋｜△φh｜ and φh″-｜△φh｜. The probability for △φh being positive （P＋） can be derived based on the Cochran distribution in direct methods. Hence the SAD phase ambiguity can be resolved by multiplying the Gaussian function peaked at φh″＋｜△φh｜ with P＋ and multiplying the Gaussian function peaked at φh″-｜△φh｜ with P_ （=1- P＋）. The direct-method SAD h phasing has been proved powerful in breaking SAD phase ambiguities, in particular when anomalous-scattering signals are weak. However, the approximation of bimodal phase distributions by the sum of two Gaussian functions introduces considerable errors. In this paper we show that a much better approximation can be achieved by replacing the two Gaussian functions with two von Mises distributions. Test results showed that this leads to significant improvement on the efficiency of direct-method SAD-phasing.
基金Project was supported by RFDP of Higher Education and NNSF of China, SF of Wuhan University
文摘In this article, by introducing characteristic singular integral operator and associate singular integral equations （SIEs）, the authors discuss the direct method of solution for a class of singular integral equations with certain analytic inputs. They obtain both the conditions of solvability and the solutions in closed form. It is noteworthy that the method is different from the classical one that is due to Lu.
文摘This paper proposes a hybrid decoupled power flow method for balanced power distribution systems with distributed generation sources. The method formulates the power flow equations in active power and reactive power decoupled form with polar coordinates. Second-order terms are included in the active power mismatch iteration, and constant Jacobian and Hessian matrices are used. A hybrid direct and indirect solution technique is used to achieve efficiency and robustness of the algorithm. Active power correction is solved by means of a sparse lower triangular and upper triangular （LU） decomposition algorithm with partial pivoting, and the reactive power correction is solved by means of restarted generalized minimal residual algorithm with an incomplete LU pre-conditioner. Typical distribution generation models and distribution load models are included. The impact of zero-impedance branches is explicitly modeled through reconfiguring of the adjacent branches with impedances. Numerical examples on a sample distribution system with widespread photovoltaic installations are given to demonstrate the effectiveness of the proposed method.
基金Project supported by the National Natural Science Foundation of China (Grant Nos 10575087 and 10875106)
文摘So far, Lou＇s direct perturbation method has been applied successfully to solve the nonlinear SchrSdinger equation（NLSE） hierarchy, such as the NLSE, the coupled NLSE, the critical NLSE, and the derivative NLSE. But to our knowledge, this method for other types of perturbed nonlinear evolution equations has still been lacking. In this paper, Lou＇s direct perturbation method is applied to the study of perturbed complex Burgers equation. By this method, we calculate not only the zero-order adiabatic solution, but also the first order modification.
文摘The general interpolation mentioned in this a rticle provides an effective way for reducing the amount of calculation of direc t optimal exploration. It has been testified by real case calculations that the interpolation is not only reliable but also can save the amount of calculation by nearly 36%. Large amount of calculation and lacking strict theoretical bas is has been th e two disadvantage of direct method by new. If this defect is not overcome, they will not only s eriously affect the application of this method, but also hinder its further rese arch. Based on sufficient calculation practice, this article has made a primary discussion about the theory and method of reducing the amount of calculation, an d has achieved some satisfactory results.
基金supported by the National Natural Science Foundations of China (Grant Nos 10735030,10475055,10675065 and 90503006)National Basic Research Program of China (Grant No 2007CB814800)+2 种基金PCSIRT (Grant No IRT0734)the Research Fund of Postdoctoral of China (Grant No 20070410727)Specialized Research Fund for the Doctoral Program of Higher Education of China (Grant No 20070248120)
文摘The approximate direct reduction method is applied to the perturbed mKdV equation with weak fourth order dispersion and weak dissipation. The similarity reduction solutions of different orders conform to formal coherence, accounting for infinite series reduction solutions to the original equation and general formulas of similarity reduction equations. Painleve Ⅱ type equations, hyperbolic secant and Jacobi elliptic function solutions are obtained for zeroorder similarity reduction equations. Higher order similarity reduction equations are linear variable coefficient ordinary differential equations.
基金Project supported by the National Natural Science Foundation of China （Grant No 11171038）.
文摘In this work, we present the direct discontinuous Galerkin （DDG） method for the one-dimensional coupled nonlinear Schr5dinger （CNLS） equation. We prove that the new discontinuous Galerkin method preserves the discrete mass conservations corresponding to the properties of the CNLS system. The ordinary differential equations obtained by the DDG space discretization is solved via a third-order stabilized Runge Kutta method. Numerical experiments show that the new DDG scheme gives stable and less diffusive results and has excellent long-time numerical behaviors for the CNLS equations.
基金supported by the Fundamental Research Funds for the Central Universities(No.NS2016003)
文摘Crack propagation path in two-directionally graded composites was investigated by the finite element method.A graded extended finite element method(XFEM)was employed to calculate displacement and stress fields in cracked graded structures.And a post-processing subroutine of interaction energy integral was implemented to extract the mixed-mode stress intensity factors(SIFs).The maximum hoop stress(MHS)criterion was adopted to predict crack growth direction based on the assumption of local homogenization of asymptotic crack-tip fields in graded materials.Effects of material nonhomogeneous parameters on crack propagation paths were also discussed in detail.It is shown that the present method can provide relatively accurate predictions of crack paths in two-directionally graded composites.Crack propagates in the decreasing direction of effective Young′s modulus.The shape and steepness of property gradient perpendicular to the crack surface have great influences on crack paths.Through redesigning material property reasonably,crack growth in graded material can be changed to improve mechanical behaviours of cracked structures.
