CFD Online Logo CFD Online URL
Home > Forums > Main CFD Forum

Spectral methods efficiency

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Display Modes
Old   December 3, 2001, 16:45
Default Spectral methods efficiency
Radu Bondar
Posts: n/a
why spectral methods are so efficient(you need only a few points to find out the right solution) on some cases?
  Reply With Quote

Old   December 4, 2001, 14:25
Default Re: Spectral methods efficiency
Patrick Godon
Posts: n/a
Hi Radu,

Spectral method a high order methods, which means that they evaluate the numerical derivative of the variables (density, velocity, pressure, energy, whatever..) using all the points of the function on the grid. An example will help you to understand best:

a simple finite difference method use the Taylor series in order to numerically evaluate the derivatie:

du/dx = (u(j+1)-u(j-1))/(2*delta x)

or also du/dx = (u(j+1)-u(j))/delta x etc....

It uses only 2 points on the grid (j+1, j-1) to evaluate the derivative.

The spectral methods use basis functions (solutions of the Strum Liouville problem) to expand the function as a series of polynomials, then the derivative of this function is obtained by taking the derivative of the polynomial and making a new series:

ex: Fourier Series:

f(x) = sum exp(-ikx) a(k) , k=-n/2+1, -1, 0, 1, ..n/2.

df(x)/dx = sum -ik exp(-ikx) a(k)

since a(k) is itself obtained from the inverse Fourier Transform and numerically the transform is discrete, a(k) is actually a function of f(x_j) where j is the index of the grid points (x_j are the grid points). Consequently, df(x)/dx is numerically evaluated using all the points on the grid and not just 2!

Now in addition, only a few polynomials are needed to represent a smooth function, and therefore high accuracy can be obtained with only a few points when the function is smooth.

When the function is discontinue, then the spectral methods are actually not good at all, and are very oscillatory (Gibbs phenomenon), which is due to the fact than when the function is not continue, the series does not converge and you actually need an infinite number of terms to represent the function. The same will happend if you try to use a Fourier series to represent a non-periodic function (for example f(x)=x cannot be represented by a Fourier series).

  Reply With Quote


Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Spectral Methods in Fluid Dynamics Martin Bailon Main CFD Forum 3 January 31, 2011 14:39
Fourier Spectral methods Q Main CFD Forum 1 March 18, 2006 23:46
spectral methods for DNS and LES ali Main CFD Forum 1 January 16, 2006 10:47
comments on FDM, FEM, FVM, SM, SEM, DSEM, BEM kenn Main CFD Forum 2 July 18, 2004 18:28
Pressure less spectral methods Gabriel Usera Main CFD Forum 1 December 19, 2000 12:59

All times are GMT -4. The time now is 06:14.