# LSFEM for incompressible flow

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 LinkBack Thread Tools Display Modes
 January 14, 2010, 13:14 LSFEM for incompressible flow #1 New Member   Sébastien Join Date: Jan 2010 Posts: 2 Rep Power: 0 Hi all, I try to write a least square finite element code to solve the incompressible navier stokes equations. The equations in u,w,p formulation are a system of 8 first order ODE, they can be writed as A1 * du/dx + A2 * du/dy + A3 * du/dz + A * u with u={vx,vy,vz,wx,wy,wz,p} or in short L*u = F. The least square finite element formulation give for the element matrix and the element vector Ke = integral of { transpose (L*Nj) * (L*Ni)} Fe = integral of { transpose (L*Nj) * F } It look simple but It don't works when I put it in practice. That is what I do : For each element For i=1 to number of nodes in the element I compute L*Nj I compute the matrix vector product Fe = L*Nj * F I assemble this local vector to the global vector For j=1 to number of nodes in the element I compute L*Ni I compute the matrix product betxeen the transpose of L*Nj and L*Ni I assemble this local matrix to the global matrix Do you see my mistake ? Maybe do you know a public source code in free access from the web using the LS-FEM strategy ? Best regards, Sébastien

 January 14, 2010, 18:28 #2 New Member   supeng Join Date: Jan 2010 Posts: 5 Rep Power: 8 I think I may help you about this. I do have one code to solve the incompressible N-S equations. Please contact me via Email if interested.

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post CD adapco Group Marketing CD-adapco 3 June 21, 2011 08:33 saii CFX 2 September 18, 2009 08:07 James Main CFD Forum 6 May 15, 2009 12:14 pxyz Main CFD Forum 37 July 7, 2006 08:42 Franck Main CFD Forum 3 September 4, 2003 05:57

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

 Contact Us - CFD Online - Top