In this paper, we propose a two-grid algorithm for solving the stream function formulation of the stationary Navies-Stokes equations. The algorithm is constructed by reducing the original system to one small, nonlinea...In this paper, we propose a two-grid algorithm for solving the stream function formulation of the stationary Navies-Stokes equations. The algorithm is constructed by reducing the original system to one small, nonlinear system on the coarse mesh space and two similar linear systems （with same stiffness matrix but different right-hand side） on the fine mesh space. The convergence analysis and error estimation of the algorithm are given for the case of conforming elements. Furthermore, the Mgorithm produces a numerical solution with the optimal asymptotic H^2-error. Finally, we give a numerical illustration to demonstrate the effectiveness of the two-grid algorithm for solving the Navier-Stokes equations.展开更多
Distributed generation (DG) is gaining in importance due to the growing demand for electrical energy and the key role it plays in reducing actual energy losses, lowering operating costs and improving voltage stability...Distributed generation (DG) is gaining in importance due to the growing demand for electrical energy and the key role it plays in reducing actual energy losses, lowering operating costs and improving voltage stability. In this paper, we propose to inject distributed power generation into a distribution system while minimizing active energy losses. This injection should be done at a grid node (which is a point where energy can be injected into or recovered from the grid) that will be considered the optimal node when total active losses in the radial distribution system are minimal. The focus is on meeting energy demand using renewable energy sources. The main criterion is the minimization of active energy losses during injection. The method used is the algorithm of bee colony (ABC) associated with Newtonian energy flow transfer equations. The method has been implemented in MATLAB for optimal node search in IEEE 14, 33 and 57 nodes networks. The active energy loss results of this hybrid algorithm were compared with the results of previous searches. This comparison shows that the proposed algorithm allows to have reduced losses with the power injected that we have found.展开更多
Based on the extraction equilibrium and mass balances in countercurrent extraction systems, a novel method was studied for dealing with the extraction equilibrium and the mass distribution in a multi-component(gamma-c...Based on the extraction equilibrium and mass balances in countercurrent extraction systems, a novel method was studied for dealing with the extraction equilibrium and the mass distribution in a multi-component(gamma-component) system. The relationships of mass distribution (x(i), y(i), i = 1, ..., lambda) between two phases were expressed by 2 lambda dimensional simultaneous equations. These simultaneous equations can be converted to a one-dimension nonlinear equation, then it was solved by Newton-Raphson algorithm within a few number of iteration. Compared with the regular calculation method for the 2 lambda dimensional simultaneous equations, Newton-Raphson algorithm can decrease the number of iteration, increase the convergence of the equations and accelerate the speed of simulation. It was verified in many multi-component systems with satisfactory results. As an example, a five-component system is demonstrated in this paper.展开更多
A new decoupled two-gird algorithm with the Newton iteration is proposed for solving the coupled Navier-Stokes/Darcy model which describes a fluid flow filtrating through porous media. Moreover the error estimate is g...A new decoupled two-gird algorithm with the Newton iteration is proposed for solving the coupled Navier-Stokes/Darcy model which describes a fluid flow filtrating through porous media. Moreover the error estimate is given, which shows that the same order of accuracy can be achieved as solving the system directly in the fine mesh when h = H2. Both theoretical analysis and numerical experiments illustrate the efficiency of the algorithm for solving the coupled problem.展开更多
Newton’s iteration is a fundamental tool for numerical solutions of systems of equations. The well-known iteration ?rapidly refines a crude initial approximation X0?to the inverse of a general nonsingular matrix. In ...Newton’s iteration is a fundamental tool for numerical solutions of systems of equations. The well-known iteration ?rapidly refines a crude initial approximation X0?to the inverse of a general nonsingular matrix. In this paper, we will extend and apply this method to n× n?structured matrices M?, in which matrix multiplication has a lower computational cost. These matrices can be represented by their short generators which allow faster computations based on the displacement operators tool. However, the length of the generators is tend to grow and the iterations do not preserve matrix structure. So, the main goal is to control the growth of the length of the short displacement generators so that we can operate with matrices of low rank and carry out the computations much faster. In order to achieve our goal, we will compress the computed approximations to the inverse to yield a superfast algorithm. We will describe two different compression techniques based on the SVD and substitution and we will analyze these approaches. Our main algorithm can be applied to more general classes of structured matrices.展开更多
The recent result of an orbit continuation algorithm has provided a rigorous method for long-term numerical integration of an orbit on the unstable manifold of a periodic solution.This algorithm is matrix-free and emp...The recent result of an orbit continuation algorithm has provided a rigorous method for long-term numerical integration of an orbit on the unstable manifold of a periodic solution.This algorithm is matrix-free and employs a combination of the Newton-Raphson method and the Krylov subspace method.Moreover,the algorithm adopts a multiple shooting method to address the problem of orbital instability due to long-term numerical integration.The algorithm is described through computing the extension of unstable manifold of a recomputed Nagata′s lowerbranch steady solution of plane Couette flow,which is an example of an exact coherent state that has recently been studied in subcritical transition to turbulence.展开更多
In this paper,we propose a Newton iterative algorithm to numerically reconstruct a locally rough surface with Dirichlet and impedance boundary conditions by near-field measurements of acoustic waves.The algorithm reli...In this paper,we propose a Newton iterative algorithm to numerically reconstruct a locally rough surface with Dirichlet and impedance boundary conditions by near-field measurements of acoustic waves.The algorithm relies on the Frechet differentiability analysis of the locally rough surface scattering problem,which is established by reducing the original model into an equivalent boundary value problem with compactly supported boundary data.With a slight modification,the algorithm can be also extended to reconstruct the local perturbation of a non-local rough surface.Finally,numerical results are presented to illustrate the effectiveness of the inversion algorithm with the multi-frequency data.展开更多
基金supported by National Foundation of Natural Science under the Grant 11071216
文摘In this paper, we propose a two-grid algorithm for solving the stream function formulation of the stationary Navies-Stokes equations. The algorithm is constructed by reducing the original system to one small, nonlinear system on the coarse mesh space and two similar linear systems （with same stiffness matrix but different right-hand side） on the fine mesh space. The convergence analysis and error estimation of the algorithm are given for the case of conforming elements. Furthermore, the Mgorithm produces a numerical solution with the optimal asymptotic H^2-error. Finally, we give a numerical illustration to demonstrate the effectiveness of the two-grid algorithm for solving the Navier-Stokes equations.
文摘Distributed generation (DG) is gaining in importance due to the growing demand for electrical energy and the key role it plays in reducing actual energy losses, lowering operating costs and improving voltage stability. In this paper, we propose to inject distributed power generation into a distribution system while minimizing active energy losses. This injection should be done at a grid node (which is a point where energy can be injected into or recovered from the grid) that will be considered the optimal node when total active losses in the radial distribution system are minimal. The focus is on meeting energy demand using renewable energy sources. The main criterion is the minimization of active energy losses during injection. The method used is the algorithm of bee colony (ABC) associated with Newtonian energy flow transfer equations. The method has been implemented in MATLAB for optimal node search in IEEE 14, 33 and 57 nodes networks. The active energy loss results of this hybrid algorithm were compared with the results of previous searches. This comparison shows that the proposed algorithm allows to have reduced losses with the power injected that we have found.
文摘Based on the extraction equilibrium and mass balances in countercurrent extraction systems, a novel method was studied for dealing with the extraction equilibrium and the mass distribution in a multi-component(gamma-component) system. The relationships of mass distribution (x(i), y(i), i = 1, ..., lambda) between two phases were expressed by 2 lambda dimensional simultaneous equations. These simultaneous equations can be converted to a one-dimension nonlinear equation, then it was solved by Newton-Raphson algorithm within a few number of iteration. Compared with the regular calculation method for the 2 lambda dimensional simultaneous equations, Newton-Raphson algorithm can decrease the number of iteration, increase the convergence of the equations and accelerate the speed of simulation. It was verified in many multi-component systems with satisfactory results. As an example, a five-component system is demonstrated in this paper.
基金supported by National Foundation of Natural Science(11471092,11326231)Zhejiang Provincial Natural Science Foundation of China(LZ13A010003)
文摘A new decoupled two-gird algorithm with the Newton iteration is proposed for solving the coupled Navier-Stokes/Darcy model which describes a fluid flow filtrating through porous media. Moreover the error estimate is given, which shows that the same order of accuracy can be achieved as solving the system directly in the fine mesh when h = H2. Both theoretical analysis and numerical experiments illustrate the efficiency of the algorithm for solving the coupled problem.
文摘Newton’s iteration is a fundamental tool for numerical solutions of systems of equations. The well-known iteration ?rapidly refines a crude initial approximation X0?to the inverse of a general nonsingular matrix. In this paper, we will extend and apply this method to n× n?structured matrices M?, in which matrix multiplication has a lower computational cost. These matrices can be represented by their short generators which allow faster computations based on the displacement operators tool. However, the length of the generators is tend to grow and the iterations do not preserve matrix structure. So, the main goal is to control the growth of the length of the short displacement generators so that we can operate with matrices of low rank and carry out the computations much faster. In order to achieve our goal, we will compress the computed approximations to the inverse to yield a superfast algorithm. We will describe two different compression techniques based on the SVD and substitution and we will analyze these approaches. Our main algorithm can be applied to more general classes of structured matrices.
文摘The recent result of an orbit continuation algorithm has provided a rigorous method for long-term numerical integration of an orbit on the unstable manifold of a periodic solution.This algorithm is matrix-free and employs a combination of the Newton-Raphson method and the Krylov subspace method.Moreover,the algorithm adopts a multiple shooting method to address the problem of orbital instability due to long-term numerical integration.The algorithm is described through computing the extension of unstable manifold of a recomputed Nagata′s lowerbranch steady solution of plane Couette flow,which is an example of an exact coherent state that has recently been studied in subcritical transition to turbulence.
文摘In this paper,we propose a Newton iterative algorithm to numerically reconstruct a locally rough surface with Dirichlet and impedance boundary conditions by near-field measurements of acoustic waves.The algorithm relies on the Frechet differentiability analysis of the locally rough surface scattering problem,which is established by reducing the original model into an equivalent boundary value problem with compactly supported boundary data.With a slight modification,the algorithm can be also extended to reconstruct the local perturbation of a non-local rough surface.Finally,numerical results are presented to illustrate the effectiveness of the inversion algorithm with the multi-frequency data.
