Fast time integration of PDEs combining POD and Galerkin projection based on a limited set of mesh points

Rapun Banzo, Maria Luisa and Vega De Prada, Jose Manuel and Terragni, Filippo (2015). Fast time integration of PDEs combining POD and Galerkin projection based on a limited set of mesh points. In: "Bifurcations and instabilities in fluid dynamics (BIFD 2015)", 15-17 jul 2015, Paris, Francia.

Description

Title: Fast time integration of PDEs combining POD and Galerkin projection based on a limited set of mesh points
Author/s:
  • Rapun Banzo, Maria Luisa
  • Vega De Prada, Jose Manuel
  • Terragni, Filippo
Item Type: Presentation at Congress or Conference (Lecture)
Event Title: Bifurcations and instabilities in fluid dynamics (BIFD 2015)
Event Dates: 15-17 jul 2015
Event Location: Paris, Francia
Title of Book: Proceedings of BIFD 2015
Date: 2015
Subjects:
Faculty: E.T.S. de Ingeniería Aeronáutica y del Espacio (UPM)
Department: Matemática Aplicada a la Ingeniería Aeroespacial
Creative Commons Licenses: Recognition - No derivative works - Non commercial

Full text

[img]
Preview
PDF - Requires a PDF viewer, such as GSview, Xpdf or Adobe Acrobat Reader
Download (1MB) | Preview

Abstract

Recently, an adaptive method to accelerate time dependent numerical solvers of systems of PDEs that require a high cost in computational time and memory has been proposed [3] (see also [1, 2]). The method combines on the fly such numerical solver with a proper orthogonal decomposition, from which we identify modes, a Galerkin projection (that provides a reduced system of equations), and the integration of the reduced system. The strategy is based on a truncation error estimate and a residual estimate, designed to control the truncation error and the mode truncation instability, respectively. These estimates support the selection of the appropiate time intervals in which the numerical solver is run to first construct and then update, on demand, the POD modes. Moreover, to reduce the computational effort needed at the outset to generate the initial POD subspace, information from former simulations or generic libraries (e.g. trigonometric functions or orthogonal polynomials) were also used. To improve the computational efficiency of the method presented in [3] a crucial step is to use a limited number of points (instead of the whole computational mesh used in the spatial discretization) to both perform POD and to Galerkin?project the equations. In this work we will discuss and compare several alternatives in representative examples illustrating that a suitable point selection can make the cost of the reduced order model (associated with POD, Galerkin projection and the integration of the resulting Galerkin system) negligible compared to that of the reference numerical solver.

More information

Item ID: 42547
DC Identifier: http://oa.upm.es/42547/
OAI Identifier: oai:oa.upm.es:42547
Deposited by: Memoria Investigacion
Deposited on: 28 Apr 2017 07:59
Last Modified: 28 Apr 2017 07:59
  • Logo InvestigaM (UPM)
  • Logo GEOUP4
  • Logo Open Access
  • Open Access
  • Logo Sherpa/Romeo
    Check whether the anglo-saxon journal in which you have published an article allows you to also publish it under open access.
  • Logo Dulcinea
    Check whether the spanish journal in which you have published an article allows you to also publish it under open access.
  • Logo de Recolecta
  • Logo del Observatorio I+D+i UPM
  • Logo de OpenCourseWare UPM