At low-Reynolds-number,the performance of airfoil is known to be greatly affected by the formation and burst of a laminar separation bubble(LSB),which requires a more precise simulation of the delicate flow structures...At low-Reynolds-number,the performance of airfoil is known to be greatly affected by the formation and burst of a laminar separation bubble(LSB),which requires a more precise simulation of the delicate flow structures.A framework based on the interior penalty discontinuous Galerkin method and large eddy simulation approach was adopted in the present study.The performances of various subgrid models,including the Smagorinsky(SM)model,the dynamic Smagorinsky(DSM)model,the wall-adapting local-eddy-viscosity(WALE)model,and the VREMAN model,have been analyzed through flow simulations of the SD7003 airfoil at a Reynolds number of 60000.It turns out that the SM model fails to predict the emergence of LSB,even modified by the Van-Driest damping function.On the contrary,the best agreement is generally achieved by the WALE model in terms of flow separation,reattachment,and transition locations,together with the aerodynamic loads.Furthermore,the influence of numerical dissipation has also been discussed through the comparison of skin friction and resolved Reynolds stresses.As numerical dissipation decreases,the prediction accuracy of the WALE model degrades.Meanwhile,nonlinear variation could be observed from the performances of the DSM model,which could be attributed to the interaction between the numerical dissipation and the subgrid model.展开更多
In this paper, on the basis of the logarithmic barrier function and KKT conditions, we propose a combined homotopy infeasible interior-point method （CHIIP） for convex nonlinear programming problems. For any convex n...In this paper, on the basis of the logarithmic barrier function and KKT conditions, we propose a combined homotopy infeasible interior-point method （CHIIP） for convex nonlinear programming problems. For any convex nonlinear programming, without strict convexity for the logarithmic barrier function, we get different solutions of the convex programming in different cases by CHIIP method.展开更多
Under the environment of electric power market, economic dispatch (ED) problem should consider network constraints, unit ramp rates, besides the basic constraints. For this problem, it is important to establish the ef...Under the environment of electric power market, economic dispatch (ED) problem should consider network constraints, unit ramp rates, besides the basic constraints. For this problem, it is important to establish the effective model and algorithm. This paper examines the decoupled conditions that affect the solution optimality to this problem. It proposes an effective model and solution method. Based on the look-ahead technique, it finds the number of time intervals to guarantee the solution optimality. Next, an efficient technique for finding the optimal solution via the interior point methods is described. Test cases, which include dispatching six units over 5 time intervals on the IEEE 30 test system with line flows and ramp constraints are presented. Results indicate that the computational effort as measured by iteration counts or execution time varies only modestly with the problem size.展开更多
This paper proposes an infeasible interior-point algorithm with full-Newton step for linear complementarity problem,which is an extension of Roos about linear optimization. The main iteration of the algorithm consists...This paper proposes an infeasible interior-point algorithm with full-Newton step for linear complementarity problem,which is an extension of Roos about linear optimization. The main iteration of the algorithm consists of a feasibility step and several centrality steps. At last,we prove that the algorithm has O(nlog n/ε) polynomial complexity,which coincides with the best known one for the infeasible interior-point algorithm at present.展开更多
Optimal adjustment algorithm for p coordinates is a generalization of the optimal pair adjustment algorithm for linear programming, which in turn is based on von Neumann’s algorithm. Its main advantages are simplicit...Optimal adjustment algorithm for p coordinates is a generalization of the optimal pair adjustment algorithm for linear programming, which in turn is based on von Neumann’s algorithm. Its main advantages are simplicity and quick progress in the early iterations. In this work, to accelerate the convergence of the interior point method, few iterations of this generalized algorithm are applied to the Mehrotra’s heuristic, which determines the starting point for the interior point method in the PCx software. Computational experiments in a set of linear programming problems have shown that this approach reduces the total number of iterations and the running time for many of them, including large-scale ones.展开更多
In this paper,we are mainly devoted to solving fixed point problems in more general nonconvex sets via an interior point homotopy method.Under suitable conditions,a constructive proof is given to prove the existence o...In this paper,we are mainly devoted to solving fixed point problems in more general nonconvex sets via an interior point homotopy method.Under suitable conditions,a constructive proof is given to prove the existence of fixed points,which can lead to an implementable globally convergent algorithm.展开更多
A new higher-order accurate space-time discontinuous Galerkin(DG)method using the interior penalty flux and discontinuous basis functions,both in space and in time,is pre-sented and fully analyzed for the second-order...A new higher-order accurate space-time discontinuous Galerkin(DG)method using the interior penalty flux and discontinuous basis functions,both in space and in time,is pre-sented and fully analyzed for the second-order scalar wave equation.Special attention is given to the definition of the numerical fluxes since they are crucial for the stability and accuracy of the space-time DG method.The theoretical analysis shows that the DG discre-tization is stable and converges in a DG-norm on general unstructured and locally refined meshes,including local refinement in time.The space-time interior penalty DG discre-tization does not have a CFL-type restriction for stability.Optimal order of accuracy is obtained in the DG-norm if the mesh size h and the time stepΔt satisfy h≅CΔt,with C a positive constant.The optimal order of accuracy of the space-time DG discretization in the DG-norm is confirmed by calculations on several model problems.These calculations also show that for pth-order tensor product basis functions the convergence rate in the L∞and L2-norms is order p+1 for polynomial orders p=1 and p=3 and order p for polynomial order p=2.展开更多
The finite-dimensional variational inequality problem (VIP) has been studied extensively in the literature because of its successful applications in many fields such as economics, transportation, regional science and ...The finite-dimensional variational inequality problem (VIP) has been studied extensively in the literature because of its successful applications in many fields such as economics, transportation, regional science and operations research. Barker and Pang[1] have given an excellent survey of theories, methods and applications of VIPs.展开更多
This work focuses on the application of the reconstruction method of differentiated backprojection （DBP）-projection onto convex sets （POCS） in the interior problem.First,we present the definition of the interior p...This work focuses on the application of the reconstruction method of differentiated backprojection （DBP）-projection onto convex sets （POCS） in the interior problem.First,we present the definition of the interior problem and real truncated Hilbert transform,and then outline the implementation steps of DBP-POCS.After that,we introduce the middle-part known condition for region of interest （ROI） accurate reconstruction and the unique condition of the interior problem,and verify the uniqueness and stability of the interior problem accurate reconstruction through numerical experiments,and then compare the results for the interior problem in reconstruction images using filtered backprojection （FBP）.In addition,the authors also design the application models of ROI reconstruction and make an initial attempt to the application of DBP-POCS method in the interior problem.展开更多
On the basis of the formulations of the logarithmic barrier function and the idea of following the path of minimizers for the logarithmic barrier family of problems the so called "centralpath" for linear pro...On the basis of the formulations of the logarithmic barrier function and the idea of following the path of minimizers for the logarithmic barrier family of problems the so called "centralpath" for linear programming, we propose a new framework of primal-dual infeasible interiorpoint method for linear programming problems. Without the strict convexity of the logarithmic barrier function, we get the following results: (a) if the homotopy parameterμcan not reach to zero,then the feasible set of these programming problems is empty; (b) if the strictly feasible set is nonempty and the solution set is bounded, then for any initial point x, we can obtain a solution of the problems by this method; (c) if the strictly feasible set is nonempty and the solution set is unbounded, then for any initial point x, we can obtain a (?)-solution; and(d) if the strictly feasible set is nonempty and the solution set is empty, then we can get the curve x(μ), which towards to the generalized solutions.展开更多
The incorrect surface current may be obtained in the vicinity of the resonant frequencies when the method of moments is used to solve either the electric or magnetic field surface integral equation. An effective met...The incorrect surface current may be obtained in the vicinity of the resonant frequencies when the method of moments is used to solve either the electric or magnetic field surface integral equation. An effective method is presented to determine the correct surface current, i.e., the correct surface current, i.e., the correct surface current is composed of the non-resonant mode current and the normalized resonant mode current multiplied by an unknown complex factor. The unknown complex factor can be obtained by employing the condition that the total field inside the conducting closed body must be zero at specified interior points. A numerical example is given for an infinitely long and perfectly conducting circular cylinder at the interior resonance, and the calculated surface currents are in good agreement with the analytical ones. The validity and accuracy of the presented method is thus verified.展开更多
Transmission line manipulations in a power system are necessary for the execution of preventative or corrective main- tenance in a network, thus ensuring the stability of the system. In this study, primal-dual interio...Transmission line manipulations in a power system are necessary for the execution of preventative or corrective main- tenance in a network, thus ensuring the stability of the system. In this study, primal-dual interior-point methods are used to minimize costs and losses in the generation and transmission of the predispatch active power flow in a hydroelectric system with previously scheduled line manipulations for preventative maintenance, over a period of twenty-four hours. The matrix structure of this problem and the modification that it imposes on the system is also broached in this study. From the computational standpoint, the effort required to solve a problem with or without line manipulations is similar, and the reasons for this are also discussed in this study. Computational results sustain our findings.展开更多
In this paper, we present a large-update primal-dual interior-point method for symmetric cone optimization（SCO） based on a new kernel function, which determines both search directions and the proximity measure betwe...In this paper, we present a large-update primal-dual interior-point method for symmetric cone optimization（SCO） based on a new kernel function, which determines both search directions and the proximity measure between the iterate and the center path. The kernel function is neither a self-regular function nor the usual logarithmic kernel function. Besides, by using Euclidean Jordan algebraic techniques, we achieve the favorable iteration complexity O（ √r（1/2）（log r）^2 log（r/ ε））, which is as good as the convex quadratic semi-definite optimization analogue.展开更多
Low-order wavefront error account for a large proportion of wave aberrations.A compensation method for low order aberration of projection lithography objective based on Interior Point Method is presented.Compensation ...Low-order wavefront error account for a large proportion of wave aberrations.A compensation method for low order aberration of projection lithography objective based on Interior Point Method is presented.Compensation model between wavefront error and degree of movable lens freedom is established.Converting over-determined system to underdetermined system,the compensation is solved by Interior Point Method(IPM).The presented method is compared with direct solve the over-determined system.Then,other algorithm GA,EA and PS is compared with IPM.Simulation and experimental results show that the presented compensation method can obtained compensation with less residuals compared with direct solve the over-determined system.Also,the presented compensation method can reduce computation time and obtain results with less residuals compare with AGA,EA and PS.Moreover,after compensation,RMS of wavefront error of the experimental lithography projection objective decrease from 56.05 nm to 17.88 nm.展开更多
In this paper, a high order compact difference scheme and a multigrid method are proposed for solving two-dimensional (2D) elliptic problems with variable coefficients and interior/boundary layers on nonuniform grids....In this paper, a high order compact difference scheme and a multigrid method are proposed for solving two-dimensional (2D) elliptic problems with variable coefficients and interior/boundary layers on nonuniform grids. Firstly, the original equation is transformed from the physical domain (with a nonuniform mesh) to the computational domain (with a uniform mesh) by using a coordinate transformation. Then, a fourth order compact difference scheme is proposed to solve the transformed elliptic equation on uniform girds. After that, a multigrid method is employed to solve the linear algebraic system arising from the difference equation. At last, the numerical experiments on some elliptic problems with interior/boundary layers are conducted to show high accuracy and high efficiency of the present method.展开更多
The generation expansion planning is one of complex mixed-integer optimization problems, which involves a large number of continuous or discrete decision variables and constraints. In this paper, an interior point wit...The generation expansion planning is one of complex mixed-integer optimization problems, which involves a large number of continuous or discrete decision variables and constraints. In this paper, an interior point with cutting plane （IP/CP） method is proposed to solve the mixed-integer optimization problem of the electrical power generation expansion planning. The IP/CP method could improve the overall efficiency of the solution and reduce the computational time. Proposed method is combined with the Bender＇s decomposition technique in order to decompose the generation expansion problem into a master investment problem and a slave operational problem. The numerical example is presented to compare with the effectiveness of the proposed algorithm.展开更多
The element-free Galerkin method is proposed to solve free vibration of rectangular plates with finite interior elastic point supports and elastically restrained edges.Based on the extended Hamilton＇s principle for t...The element-free Galerkin method is proposed to solve free vibration of rectangular plates with finite interior elastic point supports and elastically restrained edges.Based on the extended Hamilton＇s principle for the elastic dynamics system,the dimensionless equations of motion of rectangular plates with finite interior elastic point supports and the edge elastically restrained are established using the element-free Galerkin method.Through numerical calculation,curves of the natural frequency of thin plates with three edges simply supported and one edge elastically restrained,and three edges clamped and the other edge elastically restrained versus the spring constant,locations of elastic point support and the elastic stiffness of edge elastically restrained are obtained.Effects of elastic point supports and edge elastically restrained on the free vibration characteristics of the thin plates are analyzed.展开更多
It is well known that the incorrect results will be given using either the electric or magnetic field integral equation to calculate the radar cross section (RCS) of a closed body at the interior resonance. In this pa...It is well known that the incorrect results will be given using either the electric or magnetic field integral equation to calculate the radar cross section (RCS) of a closed body at the interior resonance. In this paper, an effective iterative technique is used to correct the calculated surface current density from the electric field integral equation. The radar cross section is computed for an infinite conducting circular cylinder at the interior resonance, and the obtained results are in good agreement with the analytical results. The backscattering cross section of an infinite triangular cylinder in the vicinity of a resonant frequency is also calculated. It is shown that the presence method is efficient and accurate.展开更多
基金This work was supported by the National Key R&D Program of China(Grant No.2022YFE0207000)the National Natural Science Foundation of China(Grant Nos.12372289,11972250,and 12102298)+1 种基金the China Postdoctoral Science Foundation(Grant No.2021M702443)Tianjin Natural Science Foundation(Grant No.22JCZDJC00910).
文摘At low-Reynolds-number,the performance of airfoil is known to be greatly affected by the formation and burst of a laminar separation bubble(LSB),which requires a more precise simulation of the delicate flow structures.A framework based on the interior penalty discontinuous Galerkin method and large eddy simulation approach was adopted in the present study.The performances of various subgrid models,including the Smagorinsky(SM)model,the dynamic Smagorinsky(DSM)model,the wall-adapting local-eddy-viscosity(WALE)model,and the VREMAN model,have been analyzed through flow simulations of the SD7003 airfoil at a Reynolds number of 60000.It turns out that the SM model fails to predict the emergence of LSB,even modified by the Van-Driest damping function.On the contrary,the best agreement is generally achieved by the WALE model in terms of flow separation,reattachment,and transition locations,together with the aerodynamic loads.Furthermore,the influence of numerical dissipation has also been discussed through the comparison of skin friction and resolved Reynolds stresses.As numerical dissipation decreases,the prediction accuracy of the WALE model degrades.Meanwhile,nonlinear variation could be observed from the performances of the DSM model,which could be attributed to the interaction between the numerical dissipation and the subgrid model.
文摘In this paper, on the basis of the logarithmic barrier function and KKT conditions, we propose a combined homotopy infeasible interior-point method （CHIIP） for convex nonlinear programming problems. For any convex nonlinear programming, without strict convexity for the logarithmic barrier function, we get different solutions of the convex programming in different cases by CHIIP method.
文摘Under the environment of electric power market, economic dispatch (ED) problem should consider network constraints, unit ramp rates, besides the basic constraints. For this problem, it is important to establish the effective model and algorithm. This paper examines the decoupled conditions that affect the solution optimality to this problem. It proposes an effective model and solution method. Based on the look-ahead technique, it finds the number of time intervals to guarantee the solution optimality. Next, an efficient technique for finding the optimal solution via the interior point methods is described. Test cases, which include dispatching six units over 5 time intervals on the IEEE 30 test system with line flows and ramp constraints are presented. Results indicate that the computational effort as measured by iteration counts or execution time varies only modestly with the problem size.
基金Supported by the National Natural Science Fund Finances Projects(71071119)
文摘This paper proposes an infeasible interior-point algorithm with full-Newton step for linear complementarity problem,which is an extension of Roos about linear optimization. The main iteration of the algorithm consists of a feasibility step and several centrality steps. At last,we prove that the algorithm has O(nlog n/ε) polynomial complexity,which coincides with the best known one for the infeasible interior-point algorithm at present.
文摘Optimal adjustment algorithm for p coordinates is a generalization of the optimal pair adjustment algorithm for linear programming, which in turn is based on von Neumann’s algorithm. Its main advantages are simplicity and quick progress in the early iterations. In this work, to accelerate the convergence of the interior point method, few iterations of this generalized algorithm are applied to the Mehrotra’s heuristic, which determines the starting point for the interior point method in the PCx software. Computational experiments in a set of linear programming problems have shown that this approach reduces the total number of iterations and the running time for many of them, including large-scale ones.
基金Supported by the NNSF of China(11026079)Supported by the Youth Backbone Teacher Foundation of Henan Province(173)
文摘In this paper,we are mainly devoted to solving fixed point problems in more general nonconvex sets via an interior point homotopy method.Under suitable conditions,a constructive proof is given to prove the existence of fixed points,which can lead to an implementable globally convergent algorithm.
文摘A new higher-order accurate space-time discontinuous Galerkin(DG)method using the interior penalty flux and discontinuous basis functions,both in space and in time,is pre-sented and fully analyzed for the second-order scalar wave equation.Special attention is given to the definition of the numerical fluxes since they are crucial for the stability and accuracy of the space-time DG method.The theoretical analysis shows that the DG discre-tization is stable and converges in a DG-norm on general unstructured and locally refined meshes,including local refinement in time.The space-time interior penalty DG discre-tization does not have a CFL-type restriction for stability.Optimal order of accuracy is obtained in the DG-norm if the mesh size h and the time stepΔt satisfy h≅CΔt,with C a positive constant.The optimal order of accuracy of the space-time DG discretization in the DG-norm is confirmed by calculations on several model problems.These calculations also show that for pth-order tensor product basis functions the convergence rate in the L∞and L2-norms is order p+1 for polynomial orders p=1 and p=3 and order p for polynomial order p=2.
基金The NNSF (10071031) of China and National 973 Project.
文摘The finite-dimensional variational inequality problem (VIP) has been studied extensively in the literature because of its successful applications in many fields such as economics, transportation, regional science and operations research. Barker and Pang[1] have given an excellent survey of theories, methods and applications of VIPs.
基金supported by the National Natural Science Foundation of China (Grant No.60872116)
文摘This work focuses on the application of the reconstruction method of differentiated backprojection （DBP）-projection onto convex sets （POCS） in the interior problem.First,we present the definition of the interior problem and real truncated Hilbert transform,and then outline the implementation steps of DBP-POCS.After that,we introduce the middle-part known condition for region of interest （ROI） accurate reconstruction and the unique condition of the interior problem,and verify the uniqueness and stability of the interior problem accurate reconstruction through numerical experiments,and then compare the results for the interior problem in reconstruction images using filtered backprojection （FBP）.In addition,the authors also design the application models of ROI reconstruction and make an initial attempt to the application of DBP-POCS method in the interior problem.
文摘On the basis of the formulations of the logarithmic barrier function and the idea of following the path of minimizers for the logarithmic barrier family of problems the so called "centralpath" for linear programming, we propose a new framework of primal-dual infeasible interiorpoint method for linear programming problems. Without the strict convexity of the logarithmic barrier function, we get the following results: (a) if the homotopy parameterμcan not reach to zero,then the feasible set of these programming problems is empty; (b) if the strictly feasible set is nonempty and the solution set is bounded, then for any initial point x, we can obtain a solution of the problems by this method; (c) if the strictly feasible set is nonempty and the solution set is unbounded, then for any initial point x, we can obtain a (?)-solution; and(d) if the strictly feasible set is nonempty and the solution set is empty, then we can get the curve x(μ), which towards to the generalized solutions.
文摘The incorrect surface current may be obtained in the vicinity of the resonant frequencies when the method of moments is used to solve either the electric or magnetic field surface integral equation. An effective method is presented to determine the correct surface current, i.e., the correct surface current, i.e., the correct surface current is composed of the non-resonant mode current and the normalized resonant mode current multiplied by an unknown complex factor. The unknown complex factor can be obtained by employing the condition that the total field inside the conducting closed body must be zero at specified interior points. A numerical example is given for an infinitely long and perfectly conducting circular cylinder at the interior resonance, and the calculated surface currents are in good agreement with the analytical ones. The validity and accuracy of the presented method is thus verified.
文摘Transmission line manipulations in a power system are necessary for the execution of preventative or corrective main- tenance in a network, thus ensuring the stability of the system. In this study, primal-dual interior-point methods are used to minimize costs and losses in the generation and transmission of the predispatch active power flow in a hydroelectric system with previously scheduled line manipulations for preventative maintenance, over a period of twenty-four hours. The matrix structure of this problem and the modification that it imposes on the system is also broached in this study. From the computational standpoint, the effort required to solve a problem with or without line manipulations is similar, and the reasons for this are also discussed in this study. Computational results sustain our findings.
基金Supported by the Natural Science Foundation of Hubei Province(2008CDZD47)
文摘In this paper, we present a large-update primal-dual interior-point method for symmetric cone optimization（SCO） based on a new kernel function, which determines both search directions and the proximity measure between the iterate and the center path. The kernel function is neither a self-regular function nor the usual logarithmic kernel function. Besides, by using Euclidean Jordan algebraic techniques, we achieve the favorable iteration complexity O（ √r（1/2）（log r）^2 log（r/ ε））, which is as good as the convex quadratic semi-definite optimization analogue.
文摘Low-order wavefront error account for a large proportion of wave aberrations.A compensation method for low order aberration of projection lithography objective based on Interior Point Method is presented.Compensation model between wavefront error and degree of movable lens freedom is established.Converting over-determined system to underdetermined system,the compensation is solved by Interior Point Method(IPM).The presented method is compared with direct solve the over-determined system.Then,other algorithm GA,EA and PS is compared with IPM.Simulation and experimental results show that the presented compensation method can obtained compensation with less residuals compared with direct solve the over-determined system.Also,the presented compensation method can reduce computation time and obtain results with less residuals compare with AGA,EA and PS.Moreover,after compensation,RMS of wavefront error of the experimental lithography projection objective decrease from 56.05 nm to 17.88 nm.
文摘In this paper, a high order compact difference scheme and a multigrid method are proposed for solving two-dimensional (2D) elliptic problems with variable coefficients and interior/boundary layers on nonuniform grids. Firstly, the original equation is transformed from the physical domain (with a nonuniform mesh) to the computational domain (with a uniform mesh) by using a coordinate transformation. Then, a fourth order compact difference scheme is proposed to solve the transformed elliptic equation on uniform girds. After that, a multigrid method is employed to solve the linear algebraic system arising from the difference equation. At last, the numerical experiments on some elliptic problems with interior/boundary layers are conducted to show high accuracy and high efficiency of the present method.
文摘The generation expansion planning is one of complex mixed-integer optimization problems, which involves a large number of continuous or discrete decision variables and constraints. In this paper, an interior point with cutting plane （IP/CP） method is proposed to solve the mixed-integer optimization problem of the electrical power generation expansion planning. The IP/CP method could improve the overall efficiency of the solution and reduce the computational time. Proposed method is combined with the Bender＇s decomposition technique in order to decompose the generation expansion problem into a master investment problem and a slave operational problem. The numerical example is presented to compare with the effectiveness of the proposed algorithm.
基金Project supported by the National Natural Science Foundation of China (Grant No.10872163)the Natural Science Foundation of Education Department of Shaanxi Province (Grant No.08JK394)
文摘The element-free Galerkin method is proposed to solve free vibration of rectangular plates with finite interior elastic point supports and elastically restrained edges.Based on the extended Hamilton＇s principle for the elastic dynamics system,the dimensionless equations of motion of rectangular plates with finite interior elastic point supports and the edge elastically restrained are established using the element-free Galerkin method.Through numerical calculation,curves of the natural frequency of thin plates with three edges simply supported and one edge elastically restrained,and three edges clamped and the other edge elastically restrained versus the spring constant,locations of elastic point support and the elastic stiffness of edge elastically restrained are obtained.Effects of elastic point supports and edge elastically restrained on the free vibration characteristics of the thin plates are analyzed.
基金This project was supported by the Foundation of MOE of China (No. 00179).
文摘It is well known that the incorrect results will be given using either the electric or magnetic field integral equation to calculate the radar cross section (RCS) of a closed body at the interior resonance. In this paper, an effective iterative technique is used to correct the calculated surface current density from the electric field integral equation. The radar cross section is computed for an infinite conducting circular cylinder at the interior resonance, and the obtained results are in good agreement with the analytical results. The backscattering cross section of an infinite triangular cylinder in the vicinity of a resonant frequency is also calculated. It is shown that the presence method is efficient and accurate.
