Increasing the order of convergence of iterative schemes for solving nonlinear systems

[EN] A set of multistep iterative methods with increasing order of convergence is presented, for solving systems of nonlinear equations. One of the main advantages of these schemes is to achieve high order of convergence with few Jacobian and functional evaluations, joint with the use of the same matrix of coefficients in the most of the linear systems involved in the process. Indeed, the application of the pseudocomposition technique on these proposed schemes allows us to increase their order of convergence, obtaining new high-order, efficient methods. Finally, some numerical tests are performed in order to check their practical behavior. (C) 2012 Elsevier B.V. All rights reserved.This research was supported by Ministerio de Ciencia y Tecnología MTM2011-28636-C02-02 and FONDOCYT 2011-1-B1-33 República DominicanaCordero Barbero, A.; Torregrosa Sánchez, JR.; Penkova Vassileva, M. (2013). Increasing the order of convergence of iterative schemes for solving nonlinear systems. Journal of Computational and Applied Mathematics. 252:86-94. https://doi.org/10.1016/j.cam.2012.11.024S869425

On adomian based numerical schemes for euler and navier-stokes equations, and application to aeroacoustic propagation

140 p.En esta tesis se ha desarrollado un nuevo método de integración en tiempo de tipo derivadas sucesivas (multiderivative), llamado ABS y basado en el algoritmo de Adomian. Su motivación radica en la reducción del coste de simulación para problemas en aeroacústica, muy costosos por su naturaleza transitoria y requisitos de alta precisión. El método ha sido satisfactoriamente empleado en ambas partes de un sistema híbrido, donde se distinguen la parte aerodinámica y la acústica.En la parte aerodinámica las ecuaciones de Navier-Stokes incompresibles son resueltas con unmétodo de proyección clásico. Sin embargo, la fase de predicción de velocidad ha sido modificadapara incluir el método ABS en combinación con dos métodos: una discretización espacial MAC devolúmenes finitos, y también con un método de alto orden basado en ADER. El método se ha validado respecto a los problemas (en 2D) de vórtices de Taylor-Green, y el desarrollo de vórticesde Karman en un cilindro cuadrado. La parte acústica resuelve la propagación de ondas descritaspor las ecuaciones linearizadas de Euler, empleando una discretización de Galerkin discontinua. El método se ha validado respecto a la ecuación de Burgers.El método ABS es sencillo de programar con una formulación recursiva. Los resultados demuestran que su sencillez junto con sus altas capacidades de adaptación lo convierten en un método fácilmente extensible a órdenes altos, a la vez que reduce el coste comparado con otros métodos clásicos

On Adomian Based Numerical Schemes for Euler and Navier-Stokes Equations, and Application to Aeroacoustic Propagation

In this thesis, an Adomian Based Scheme (ABS) for the compressible Navier-Stokes equations is constructed, resulting in a new multiderivative type scheme not found in the context of fluid dynamics. Moreover, this scheme is developed as a means to reduce the computational cost associated with aeroacoustic simulations, which are unsteady in nature with high-order requirements for the acoustic wave propagation. We start by constructing a set of governing equations for the hybrid computational aeroacoustics method, splitting the problem into two steps: acoustic source computation and wave propagation. The first step solves the incompressible Navier-Stokes equation using Chorin's projection method, which can be understood as a prediction-correction method. First, the velocity prediction is obtained solving the viscous Burgers' equation. Then, its divergence-free correction is performed using a pressure Poisson type projection. In the velocity prediction substep, Burgers' equation is solved using two ABS variants: a MAC type implementation, and a ``modern'' ADER method. The second step in the hybrid method, related to wave propagation, is solved combining ABS with the discontinuous Galerkin high-order approach. Described solvers are validated against several test cases: vortex shedding and Taylor-Green vortex problems for the first step, and a Gaussian wave propagation in the second case. Although ABS is a multiderivative type scheme, it is easily programmed with an elegant recursive formulation, even for the general Navier-Stokes equations. Results show that its simplicity combined with excellent adaptivity capabilities allows for a successful extension to very high-order accuracy at relatively low cost, obtaining considerable time savings in all test cases considered.Predoc Gobierno Vasc

On adomian based numerical schemes for euler and navier-stokes equations, and application to aeroacoustic propagation

140 p.En esta tesis se ha desarrollado un nuevo método de integración en tiempo de tipo derivadas sucesivas (multiderivative), llamado ABS y basado en el algoritmo de Adomian. Su motivación radica en la reducción del coste de simulación para problemas en aeroacústica, muy costosos por su naturaleza transitoria y requisitos de alta precisión. El método ha sido satisfactoriamente empleado en ambas partes de un sistema híbrido, donde se distinguen la parte aerodinámica y la acústica.En la parte aerodinámica las ecuaciones de Navier-Stokes incompresibles son resueltas con unmétodo de proyección clásico. Sin embargo, la fase de predicción de velocidad ha sido modificadapara incluir el método ABS en combinación con dos métodos: una discretización espacial MAC devolúmenes finitos, y también con un método de alto orden basado en ADER. El método se ha validado respecto a los problemas (en 2D) de vórtices de Taylor-Green, y el desarrollo de vórticesde Karman en un cilindro cuadrado. La parte acústica resuelve la propagación de ondas descritaspor las ecuaciones linearizadas de Euler, empleando una discretización de Galerkin discontinua. El método se ha validado respecto a la ecuación de Burgers.El método ABS es sencillo de programar con una formulación recursiva. Los resultados demuestran que su sencillez junto con sus altas capacidades de adaptación lo convierten en un método fácilmente extensible a órdenes altos, a la vez que reduce el coste comparado con otros métodos clásicos

A NEW ALGORITHM IN NONLINEAR ANALYSIS OF STRUCTURES USING PARTICLE SWARM OPTIMIZATION

Solving systems of nonlinear equations is a difficult problem in numerical computation. Probably the best known and most widely used algorithm to solve a system of nonlinear equations is Newton-Raphson method. A significant shortcoming of this method becomes apparent when attempting to solve problems with limit points. Once a fixed load is defined in the first step, there is no way to modify the load vector should a limit point occur within the increment. To overcome this defect, displacement control methods for passing limit points can be used. In displacement control method, the load ratio in the first step of an increment is defined so that a particular key displacement component will change only by a prescribed amount. In this paper the load ratio is obtained using particle swarm optimization (PSO) algorithm so that the complex behavior of structures can be followed, automatically. Design variable is load ratio and its unbalanced force is also considered as objective function in optimization process. Numerical results are performed under geometrical nonlinear analysis, elastic post-buckling analysis and inelastic post-buckling analysis. The efficiency and accuracy of proposed method are demonstrated by solving these examples.