Change search
CiteExportLink to record
Permanent link

Direct link
Cite
Citation style
  • apa
  • harvard1
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • Other style
More styles
Language
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Other locale
More languages
Output format
  • html
  • text
  • asciidoc
  • rtf
Aspects of the recursive projection method applied to flow calculations
KTH, School of Computer Science and Communication (CSC), Numerical Analysis and Computer Science, NADA.
2005 (English)Doctoral thesis, comprehensive summary (Other scientific)
Abstract [en]

In this thesis, we have investigated the Recursive Projection Method, RPM, as an accelerator for computations of both steady and unsteady flows, and as a stabilizer in a bifurcation analysis.

The criterion of basis extraction is discussed. It can be interpreted as a tolerance for the accuracy of the eigenspace spanned by the identified basis, alternatively it can be viewed as a criterion when the approximative Krylov sequence becomes numerically rank deficient.

Steady state calculations were performed on two different turbulent test-cases; a 2D supersonic nozzle flow with the Spalart-Allmaras 1-equation model and a 2D sub-sonic airfoil simulation using the κ - ε model. RPM accelerated the test-cases with a factor between 2 and 5.

In multi-scale problems, it is often of interest to model the macro-scale behavior, still retaining the essential features of the full systems. The ``coarse time stepper'' is a heuristic approach for circumventing the analytical derivation of models. The system studied here is a linear lattice of non-linear reaction sites coupled by diffusion. After reformulation of the time-evolution equation as a fixed-point scheme, RPM coupled with arc-length continuation is used to calculate the bifurcation diagrams of the effective (but analytically unavailable) equation.

Within the frame-work of dual time-stepping, a common approach in unsteady CFD-simulation, RPM is used to accelerate the convergence. Two test-cases were investigated; the von Karman vortex-street behind a cylinder at Re=100, and the periodic shock oscillation of a symmetric airfoil at M ∞ = 0.76 with a Reynolds number Re=11 x 106.

It was believed that once a basis had been identified, it could be retained for several steps. The simulations usually showed that the basis could only be retained for one step.

The need for updating the basis motivates the use of Krylov methods. The most common method is the (Block-) Arnoldi algorithm. As the iteration proceeds, Krylov methods become increasingly expensive and restart is required. Two different restart algorithm were tested. The first is that of Lehoucq and Maschhoff, which uses a shifted QR iteration, the second is a block extension of the single-vector Arnoldi method due to Stewart. A flexible hybrid algorithm is derived combining the best features of the two.

Place, publisher, year, edition, pages
Stockholm: KTH , 2005. , ix, 26 p.
Series
Trita-NA, ISSN 0348-2952 ; 0444
Keyword [en]
Datorteknik, applied mechanics, computer science, aerospace
Keyword [sv]
Datorteknik
National Category
Computer Engineering
Identifiers
URN: urn:nbn:se:kth:diva-101ISBN: 91-7283-940-6 (print)OAI: oai:DiVA.org:kth-101DiVA: diva2:6781
Public defence
2005-01-20, Sal L1, Drottning Kristinas väg 30, Stockholm, 10:15
Opponent
Supervisors
Note

QC 20101015

Available from: 2005-02-02 Created: 2005-02-02 Last updated: 2012-09-21Bibliographically approved
List of papers
1. The Recursive Projection Method Applied to Steady-state CFD Calculations
Open this publication in new window or tab >>The Recursive Projection Method Applied to Steady-state CFD Calculations
2004 (English)Report (Other academic)
Place, publisher, year, edition, pages
Stockholm: KTH, 2004. 27 p.
Series
Trita-NA, ISSN 0348-2952 ; 04445
National Category
Computer Engineering
Identifiers
urn:nbn:se:kth:diva-4904 (URN)
Note
QC 20101015Available from: 2005-02-02 Created: 2005-02-02 Last updated: 2010-10-15Bibliographically approved
2. Equation-free, effective computation for discrete systems: a time stepper based approach
Open this publication in new window or tab >>Equation-free, effective computation for discrete systems: a time stepper based approach
Show others...
2005 (English)In: International Journal of Bifurcation and Chaos in Applied Sciences and Engineering, ISSN 0218-1274, Vol. 15, no 3, 975-996 p.Article in journal (Refereed) Published
Abstract [en]

We propose a computer-assisted approach to studying the effective continuum behavior of spatially discrete evolution equations. The advantage of the approach is that the "coarse model" (the continuum, effective equation) need not be explicitly constructed. The method only uses a time-integration code for the discrete problem and judicious choices of initial data and integration times; our bifurcation computations are based on the so-called Recursive Projection Method (RPM) with arc-length continuation [Shroff & Keller, 1993]. The technique is used to monitor features of the genuinely discrete problem such as the pinning of coherent structures and its results are compared to quasi-continuum approaches such as the ones based on Pade approximations.

Keyword
equation-free methods, homogenization, discrete problems, bifurcation, pinning condition
National Category
Computer Engineering
Identifiers
urn:nbn:se:kth:diva-4905 (URN)10.1142/S0218127405012399 (DOI)000229230400016 ()2-s2.0-21144445075 (Scopus ID)
Note
Workshop on Modeling and Computations in Dynamical Systems Budapest, HUNGARY, OCT, 2003 QC 20101015Available from: 2005-02-02 Created: 2005-02-02 Last updated: 2017-12-06Bibliographically approved
3. Recursive Projection Method for efficient unsteady CFD simulations
Open this publication in new window or tab >>Recursive Projection Method for efficient unsteady CFD simulations
2004 (English)In: / [ed] P. Neittaanmäki, T. Rossi, K. Majava, and O. Pironneau, 2004, 1-20 p.Conference paper, Published paper (Refereed)
Abstract [en]

The Recursive Projection Method (RPM) has been implemented into the unstructured grid CFD code EDGE to accelerate the inner-loop convergence of dual time stepping. The method tries to identify the slowly converging subspace and applies Newton iterations in this subspace together with a fixed point scheme in the complement. The method has been employed to compute the steady and unsteady viscous flow around a circular cylinder for a Reynolds number of 100. When converging to machine accuracy, RPM accelerated the convergence of the steady-state solution by a factor of 2.5. The time-accurate simulations were accelerated by a factor of about two.

Keyword
recursive Projection Method, unsteady aerodynamics, CFD, convergence, dual time stepping, vortex street
National Category
Computer Engineering
Identifiers
urn:nbn:se:kth:diva-4925 (URN)2-s2.0-84893505605 (Scopus ID)
Conference
European Congress on Computational Methods in Applied Sciences and Engineering, ECCOMAS 2004; Jyvaskyla; Finland; 24-28 July 2004
Note

QC 20141210

Available from: 2005-02-11 Created: 2005-02-11 Last updated: 2014-12-10Bibliographically approved
4. Evaluation of the Recursive Projection Method for Efficient Unsteady Turbulent CFD Simulations
Open this publication in new window or tab >>Evaluation of the Recursive Projection Method for Efficient Unsteady Turbulent CFD Simulations
2004 (English)In: 24th INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES, 2004, 1-13 p.Conference paper, Published paper (Refereed)
Abstract [en]

The Recursive Projection Method (RPM) hasbeen implemented into an unstructured CFD code to improve the efficiency of dual time steppingfor unsteady turbulent CFD simulations.RPM is a combined implicit-explicit method that enhances convergence. It can easily be implementedinto existing codes and the solver’s existing acceleration techniques can be used withoutchange. The method has been evaluated by computing the periodic self-induced shock oscillations over an 18% thick biconvex airfoil at0◦ angle of attack, a Mach number of 0.76 anda Reynolds number of 11 million. On average,RPM accelerated the convergence of the innerloop of dual time stepping to a predefined convergencecriterion by a factor of about 2.5.

Keyword
Recursive Projection Method, CFD, dual time stepping, convergence, unsteady flow, buffet
National Category
Computer Engineering Vehicle Engineering
Identifiers
urn:nbn:se:kth:diva-4907 (URN)
Conference
24th INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES
Note
QC 20101015Available from: 2005-02-02 Created: 2005-02-02 Last updated: 2010-10-19Bibliographically approved
5. New Implementations of the Implicitly Restarted Block Arnoldi Method
Open this publication in new window or tab >>New Implementations of the Implicitly Restarted Block Arnoldi Method
2004 (English)Report (Other academic)
Place, publisher, year, edition, pages
Stockholm: KTH, 2004. 39 p.
Series
Trita-NA, ISSN 0348-2952 ; 0446
National Category
Computer Engineering
Identifiers
urn:nbn:se:kth:diva-4908 (URN)
Note
QC 20101015Available from: 2005-02-02 Created: 2005-02-02 Last updated: 2010-10-15Bibliographically approved

Open Access in DiVA

fulltext(393 kB)1608 downloads
File information
File name FULLTEXT01.pdfFile size 393 kBChecksum MD5
c48a74d9212531c50dc3fbfe5dba802c4ccfa74c87760838c730fbb6253dc431b12fb803
Type fulltextMimetype application/pdf

Search in DiVA

By author/editor
Möller, Joakim
By organisation
Numerical Analysis and Computer Science, NADA
Computer Engineering

Search outside of DiVA

GoogleGoogle Scholar
Total: 1608 downloads
The number of downloads is the sum of all downloads of full texts. It may include eg previous versions that are now no longer available

isbn
urn-nbn

Altmetric score

isbn
urn-nbn
Total: 850 hits
CiteExportLink to record
Permanent link

Direct link
Cite
Citation style
  • apa
  • harvard1
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • Other style
More styles
Language
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Other locale
More languages
Output format
  • html
  • text
  • asciidoc
  • rtf