A first application of a time domain Poincare-Krylov approach for the study of the all-important stability of periodic solutions of nonlinear power networks is reported in this work. Whilst a Newton method solves the ...
详细信息
A first application of a time domain Poincare-Krylov approach for the study of the all-important stability of periodic solutions of nonlinear power networks is reported in this work. Whilst a Newton method solves the nonlinear algebraic equations resulting from the Poincare map method, an iterative Krylov-subspace method based on a gmres algorithm is applied for solving the Newton correction equations. More importantly, a QR factorization of the Hessenberg matrix involved in the gmres least square problem is implemented using Givens rotations in order to avoid the linear growth of the computational complexity in large-scale power networks. Further, the stability of periodic solutions is determined by computing the Floquet multipliers using Ritz values and the Hessenberg matrix. Numerical tests carried out on a modified three-phase version of the IEEE 118-node system with a hydro-turbine synchronous generator and a grid-tied power converter demonstrate that speedup factors up to 8 are attainable with the Krylov-Subspace approach with respect to the standard Poincare map method. An outstanding outcome of the comparative study is that the incorporation of QR factorization to the Hessenberg least square problem outperformed the classic periodic steady-state solvers, providing further computational savings up to 50%.
Purpose - The purpose of this paper is to present a modified version of the hybrid Finite Element Method-Dirichlet Boundary Condition Iteration method for the solution of open-boundary skin effect problems. Design/met...
详细信息
Purpose - The purpose of this paper is to present a modified version of the hybrid Finite Element Method-Dirichlet Boundary Condition Iteration method for the solution of open-boundary skin effect problems. Design/methodology/approach - The modification consists of overlapping the truncation and the integration boundaries of the standard method, so that the integral equation becomes singular as in the well-known Finite Element Method-Boundary Element Method (FEM-BEM) method. The new method is called FEM-SDBCI. Assuming an unknown Dirichlet condition on the truncation boundary, the global algebraic system is constituted by the sparse FEM equations and by the dense integral equations, in which singularities arise. Analytical formulas are provided to compute these singular integrals. The global system is solved by means of a Generalized Minimal Residual iterative procedure. Findings - The proposed method leads to slightly less accurate numerical results than FEM-BEM, but the latter requires much more computing time. Practical implications - Then FEM-SDBCI appears more appropriate than FEM-BEM for applications which require a shorter computing time, for example in the stochastic optimization of electromagnetic devices. Originality/value - Note that FEM-SDBCI assumes a Dirichlet condition on the truncation boundary, whereas FEM-BEM assumes a Neumann one.
In this paper, we state in a new form the algebraic problem arising from the one-field displacement finite element method (FEM). The displacement approach, in this discrete form, can be considered as the dual approach...
详细信息
In this paper, we state in a new form the algebraic problem arising from the one-field displacement finite element method (FEM). The displacement approach, in this discrete form, can be considered as the dual approach (force or equilibrium) with subsidiary constraints. This approach dissociates the nonlinear operator to the linear ones and their sizes are linear functions of integration rule which is of interest in the case of reduced integration. This new form of the problem leads to an inexpensive improvement of FEM computations, which acts at local, elementary and global levels. We demonstrate the numerical performances of this approach which is independent of the mesh structure. Using the gmres algorithm we build, for nonsymmetric problems, a new algorithm based upon the discretized field of strain. The new algorithms proposed are more closer to the mechanical problem than the classical ones because all fields appear during the resolution process. The sizes of the different operators arising in these new forms are linear functions of integration rule, which is of great interest in the case of reduced integration. (C) 2003 Elsevier B.V. All rights reserved.
The matrix-free Newton-Krylov method that uses the gmres algorithm (an iterative algorithm for solving systems of linear algebraic equations) is used for the parametric continuation of the solitary traveling pulse sol...
详细信息
The matrix-free Newton-Krylov method that uses the gmres algorithm (an iterative algorithm for solving systems of linear algebraic equations) is used for the parametric continuation of the solitary traveling pulse solution in a three-component reaction-diffusion system. Using the results of integration on a short time interval, we replace the original system of nonlinear algebraic equations by another system that has more convenient (from the viewpoint of the spectral properties of the gmres algorithm) Jacobi matrix. The proposed parametric continuation proved to be efficient for large-scale problems, and it made it possible to thoroughly examine the dependence of localized solutions on a parameter of the model.
This paper deals with the resolution of the time-dependent velocity-vorticity-pressure formulation of Stokes problem in a two- and three-dimensional domain. The discretization is based on the spectral element method w...
详细信息
This paper deals with the resolution of the time-dependent velocity-vorticity-pressure formulation of Stokes problem in a two- and three-dimensional domain. The discretization is based on the spectral element method with respect to the space variable and Euler's implicit scheme with respect to the time variable. We implement the asymmetric linear system to solve the full spectral element discrete problem. Detailed numerical experiments lead to confirm the usefulness of the discretization.
This paper discusses the methods of imposing symmetry in the augmented system formulation (ASF) for least-squares (LS) problems. A particular emphasis is on upper Hessenberg problems, where the challenge lies in leavi...
详细信息
This paper discusses the methods of imposing symmetry in the augmented system formulation (ASF) for least-squares (LS) problems. A particular emphasis is on upper Hessenberg problems, where the challenge lies in leaving all zero-by-definition elements of the LS matrix unperturbed. Analytical solutions for optimal perturbation matrices are given, including upper Hessenberg matrices. Finally, the upper Hessenberg LS problems represented by unsymmetric ASF that indicate a normwise backward stability of the problem (which is not the case in general) are identified. It is observed that such problems normally arise from Arnoldi factorization (for example, in the generalized minimal residual (gmres) algorithm). The problem is illustrated with a number of practical (arising in the gmres algorithm) and some purpose-built' examples. Copyright (c) 2014 John Wiley & Sons, Ltd.
In this paper, the formulation of two-dimensional finite difference frequency domain (FDFD) for heterogeneous and lossy materials is presented. But the main attention is focus on the efficient formulation of the matri...
详细信息
ISBN:
(纸本)9781509028009
In this paper, the formulation of two-dimensional finite difference frequency domain (FDFD) for heterogeneous and lossy materials is presented. But the main attention is focus on the efficient formulation of the matrix solver. Contemporary, the numerical models are consisted of a large-scale sparse complex matrix, therefore the specific formulation of the solver enable to improve the properties of the proposed algorithm. Two iterative subroutines for the complex equation are presented and validated, i.e. the BiConjugate Gradient Stabilized method (BiCGStab) and the Generalized Minimum Residual method with restarts (gmres(m)). The stability and performance of these algorithms are tested using model based on typical building construction including concrete column.
In this paper, the Cauchy problem for the Helmholtz equation is investigated. The objective is to recover the missing data on some part of the boundary of a bounded domain from overspecified data on the remaining part...
详细信息
In this paper, the Cauchy problem for the Helmholtz equation is investigated. The objective is to recover the missing data on some part of the boundary of a bounded domain from overspecified data on the remaining part of the boundary. We propose a preconditioned Krylov algorithm to solve this ill-posed problem, based on the represen-tation of the solution with surface integral equations and the Steklov-Poincare operator. We give a theoretical and numerical validation of the proposed method conducted in the 3D setting. We show the fast convergence of the proposed algorithm tested on various synthetic examples. The numerical results show a high precision of the reconstruction obtained for different levels of noisy data.(c) 2022 Elsevier B.V. All rights reserved.
Purpose - This paper compares the hybrid FEM-BEM and FEM-DBCI methods for the solution of open-boundary static and quasi-static electromagnetic field problems. Design/methodology/approach - After a brief review of the...
详细信息
Purpose - This paper compares the hybrid FEM-BEM and FEM-DBCI methods for the solution of open-boundary static and quasi-static electromagnetic field problems.
Design/methodology/approach - After a brief review of the two methods (both coupling a differential equation for the interior problem with an integral equation for the exterior one), they are compared in terms of accuracy, memory and computing time requirements by means of a set of simple examples.
Findings - The comparison suggests that FEM-BEM is more accurate than FEM-DBCI but requires more computing time.
Practical implications - Then FEM-DBCI appears more appropriate for applications which require a shorter computing time, for example in the stochastic optimization of electromagnetic devices. Conversely, FEM-BEM is more appropriate in cases in which a high level of precision is required in a single computation.
Originality/value - Note that the FEM-BEM considered in this paper is a non standard one in which the nodes of the normal derivative on the truncation boundary are placed in positions different from those of the potential.
暂无评论