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

How to formulate the FEM for the governing equatio

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

LinkBack Thread Tools Display Modes
Old   March 14, 2005, 15:54
Default How to formulate the FEM for the governing equatio
Posts: n/a
Hi there,

The governing equations are already in the integral form, say,

d/dt (S udv) = S (f dv).

where u is the velocity, f is the force and (S dv) is the volume integral.

I am wondering how I can derive the FEM formulation through the standard Galerkin residual weighting procedure?


  Reply With Quote

Old   March 15, 2005, 05:39
Default Re: How to formulate the FEM for the governing equ
Posts: n/a
Let's take a step back. Suppose your PDE is Lu=f, with u the unknown field, f - a known RHS and L - a differential operator. Appropriate BCs should also be prescribed, but for this simple explanation they will be ignored. The domain should be meshed with nonoverlapping elements.

You start by assuming the distribution of u in any element is spanned by the shape functions, Ni, which the spatial distribution within the element and the DOFs, Ui, usually the nodal values, i.e., u = sum (Ni*Ui) (the sum is on i).

Next, you multiply your PDE by a weight function, Wk, and integrate over the domain volume to get the weak formulation. Since the volume is meshed, the integral over the entire domain may be split to sum on all elements of volume integrals over each element. Further, for the standard Galerkin method, the weight functions are identical to the shape functions, i.e., Wk=Nk.

That's more-or-less the whole process, except for the BC treatment. There are, of course, many details involved (and God is in the details) in order to end up with a successful code. For these, a good book rather than a posting in a forum is advised.

I hope this helps,

  Reply With Quote

Old   March 15, 2005, 13:32
Default Re: How to formulate the FEM for the governing equ
Posts: n/a
Hi Rami,

Thanks kindly for your reply.

In fact, I am fully aware of the procedure you described and have been using the FEM for a number of years. However, my original question was about the governing equations that are already in the integral form before the weighting (i.e. the volume conservative form). It is no problem to span variables elementwise via shape functions. My concern is that how to do the weighting. One cannot simply move the weighting function into the integral unless using piecewise constant weighting function. However, by choosing the piecewise constant weighting function, the FEM reduces to FVM in certain sense.

My question, after rephrasing, is that can we use linear or higher-order weighting functions for the governing equations in the volume conservative form? If the answer is no, then does it mean that the FEM is not so conservative?

Thanks again.

  Reply With Quote

Old   March 16, 2005, 08:20
Default Re: How to formulate the FEM for the governing equ
Posts: n/a

i am the student of final year mechanical. sir,we are doing project on cfd,our project details shown below. project titleEVELOPMENT OF PREPROCESSOR FOR CFD ANALYSIS ON 2D DIFFUSERS(ST.CONICAL,CURVED DIFFUSERS) sir we have problem with deriving eqation for fluid flow through these diffusers .sir i will be great full to you if you kindly help me.we want to extend our project ,so that we can call our project as complete package for diffusers.we are writing programmes in matlab.
  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
I developed an FEM toolkit in Java: FuturEye nkliuyueming Main CFD Forum 7 January 29, 2016 14:28
comments on FDM, FEM, FVM, SM, SEM, DSEM, BEM kenn Main CFD Forum 2 July 18, 2004 18:28
FEM or FVM for CFD Astrid Main CFD Forum 18 December 15, 2000 01:02
How good is CFD? kai Main CFD Forum 39 April 7, 2000 12:48

All times are GMT -4. The time now is 21:05.