1 / 28

Lesson 5

Lesson 5. Method of Weighted Residuals. Classical Solution Technique. The fundamental problem in calculus of variations is to obtain a function f(x) such that small variations in the function d f(x) will not change the original function

elan
Download Presentation

Lesson 5

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Lesson 5 Method of Weighted Residuals

  2. Classical Solution Technique The fundamental problem in calculus of variations is to obtain a function f(x) such that small variations in the function df(x) will not change the original function The variational function can be written in general form for a second-order governing equation (no first derivatives) as Where a,b, and g are prescribed values

  3. Classical solution (continued) An equation containing first-order derivatives may not have a corresponding variational function. In some cases, a pseudovariational function can be used where C=C(x)

  4. Classical solution (continued) Example: Consider a rod of length L. The equation defining heat transfer in the rod is with boundary conditions Integrating twice, one obtains Applying boundary conditions, the final result is

  5. Rayleigh-Ritz Method FEM variational approach attributed to Lord Rayleigh (1842-1919) & Walter Ritz (1878-1909) Let Assume a quadratic function with boundary conditions

  6. R-R Method (continued) Hence, Now integrate Thus

  7. R-R Method (continued) which becomes To find the value of C3 that makes J a minimum, Therefore,

  8. Variation Methods Given a function u(x), the following constraints must be met • (1) satisfies the constraints u(x1)=u1 and u(x2)=u2 • (2) is twice differentiable in x1<x<x2 • (3) minimizes the integral

  9. Variational methods (continued) Then it can be shown that u(x) is also the solution of the Euler-Lagrange equation where

  10. Variational methods (continued) For higher derivatives of u, Hence,

  11. Variational methods (continued) In 2-D, the constraints are • (1) satisfies the constraint u = u0 on G • (2) is twice differentiable in domain A(x,y) • (3) minimizes the functional and

  12. Variational methods (continued) Example: Find the functional statement for the 2-D heat diffusion equation Applying the Euler-Lagrange relation

  13. Variational methods (continued) we find that which yields the final functional form

  14. A Rayleigh-Ritz Example Begin with the equation with boundary conditions First find the variational statement (J)

  15. R-R example (continued) The variational statement is Assume a quadratic approximation with boundary conditions

  16. R-R example (continued) Thus, Now, To be a minimum, Finally

  17. The Weak Statement • Method of Weighted Residuals - one does not need a strong mathematical background to use FEM. However, one must be able to integrate. • To illustrate the MWR, let us begin with a simple example - Conduction of heat in a rod of length L with source term Q.

  18. weak statement (continued) Integrating, or This analytical solution serves as a useful benchmark for verifying the numerical approach.

  19. weak statement (continued) There are basically two ways to numerically solve this equation using the FEM: Rayleigh – Ritz Method and the Galerkin Method ( which produces a “weak” statement) Consider

  20. weak statement (continued) Since u = cifi(x,y) is an approximate function, substitution into the above equation may not satisfy the equation. We set the equation equal to an error (e) We now introduce a set of weighting functions (test functions) Wi, and construct an inner product (Wi,e) that is set to zero – this forces the error of the approximate differential equation to zero (average).

  21. weak statement (continued) Hence, Example: The inner product becomes

  22. weak statement (continued) Integrating by parts (Green-Gauss Theorem) Problem: Use Galerkin’s Method to solve

  23. weak statement (continued) The inner product is or The weak statement becomes

  24. weak statement (continued) Since we obtain This is the same as the Rayleigh-Ritz Method but no variational principle is required.

  25. Weighting Function Choices • Galerkin • Least Squares • Method of Moments

  26. Weighting Function Choices (continued) • Collocation • Sub domain

  27. Least Squares Example Let For the previous example problem

  28. Least Squares Example (continued) Thus, Some observations: (1) no integration by parts required – not a weak form (2) the Neumann boundary conditions do not appear naturally (3) in this case fi must satisfy global boundary conditions, which is difficult in most problems

More Related