Petrov–Galerkin method

The Petrov–Galerkin method is a mathematical method used to obtain approximate solutions of partial differential equations which contain terms with odd order. In these type of problems a weak formulation with similar function space for test function and solution function is not possible. Hence the method is used in case the test function and solution function belong to different function spaces.[1]

Overview

An example of differential equation containing a term with odd order is as follows:

a(x)\dfrac{\mathrm{d}u}{\mathrm{d}x} + b(x)\dfrac{\mathrm{d}^2u}{\mathrm{d}x^2} = f(x),\quad x\in (0,L) \quad\&\quad u(0)=u_o, \left.\dfrac{\mathrm{d}u}{\mathrm{d}x}\right|_{x=L}=u_L'

If a test function v(x) is used to obtain the weak form, the final Galerkin formulation will be given as follows:

\int_0^L a(x)v(x)\dfrac{\mathrm{d}u}{\mathrm{d}x}\mathrm{d}x-\int_0^L b(x)\dfrac{\mathrm{d}v}{\mathrm{d}x}\dfrac{\mathrm{d}u}{\mathrm{d}x}\mathrm{d}x + \left[b(x)v\dfrac{\mathrm{d}u}{\mathrm{d}x}\right]_0^L = \int_0^L v(x)f(x) \, \mathrm{d}x

The term with even order (2nd term in LHS) is now symmetric, as the test function and solution function both have same order of differentiation and they both belong to H_0^1. However, there is no way the first term on LHS can be made this way. In this case the solution space H_0^1 and test function space L^2 are different and hence the usually employed Bubnov Galerkin method cannot be used.

See also

Notes

  1. J. N. Reddy: An introduction to the finite element method, 2006, Mcgraw–Hill


This article is issued from Wikipedia - version of the Tuesday, April 12, 2016. The text is available under the Creative Commons Attribution/Share Alike but additional terms may apply for the media files.