|
[Sponsors] |
FEM Implementation of pressure-correction scheme |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
December 3, 2006, 08:05 |
FEM Implementation of pressure-correction scheme
|
#1 |
Guest
Posts: n/a
|
Hi All,
I have been playing around with a FEM implementation of a simple pressure-correction scheme for unsteady incompressible flow. Being a newbie to CFD, I noticed some things which I do not fully understand. First of all, when using a triangular grid w/ linear shape function (both velocity and pressure), the scheme works when directly discretizing the (continuous) Laplacian by a simple Galerkin approach (using partial integration). On the other hand, when using the often-cited C^T M^-1 C matrix, the pressure field experiences pure oscillation ( there is no sign of a pressure distribution similar to what I get when directly discretizing the Laplacian). I tested this with T1P0 elements, but without any luck. I created a small MatLab implementation of the scheme, which is accessible here: http://www.geocities.com/secondresort/nsexample.zip This code uses the Galerkin discretization of the continuous Laplacian and employs equal-order interpolation for velocity and pressure. If you follow the instructions in the tunnel.m file and let the system compute the flow field (as given in the example parameters), you notice a build-up of "pressure" near the inlet at around timestep 1600 (yes, you have to wait a long time, set up some coffee inbetween). Do I take it right, that this is most probably caused by an insufficiently fine discretization near the boundary ? With best regards, Markus |
|
December 3, 2006, 15:21 |
Re: FEM Implementation of pressure-correction sche
|
#2 |
Guest
Posts: n/a
|
I found two serious bugs in the code (one should think that this is impossible in only a couple of lines) which are now corrected. I uploaded the new version, looks better now.
Regards, Markus |
|
December 4, 2006, 09:56 |
Re: FEM Implementation of pressure-correction sche
|
#3 |
Guest
Posts: n/a
|
I selected the link you provided but was unable to get the download, both yesterday and today. No data is transferred and I get the the message "There is a temporary network error preventing download of your file." Has anyone else tried to get this download?
I too use Matlab for incompressible flow, but on rectangles with a pressureless, divergence-free method. |
|
December 4, 2006, 11:00 |
Re: FEM Implementation of pressure-correction sche
|
#4 |
Guest
Posts: n/a
|
I've just downloaded the file without any problem.
|
|
January 6, 2007, 02:53 |
Segregated FEM of Steady Navier-Stokes Equation
|
#5 |
Guest
Posts: n/a
|
Hi, It would be a great if anybody can help me in this regard. I have done enough of literature survey and at present having a 20+ papers on it. I have implemented the procedure but could not get the result.
if anyone implemented kindly contact me at yenyes_kumar@yahoo.co.in Thanks in advance ! NSK |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Implementation of QUICK scheme | Romuald Skoda | Main CFD Forum | 11 | November 6, 2017 22:20 |
Divergence detected in AMG solver: pressure correction | emlejeen | FLUENT | 5 | December 15, 2016 00:47 |
changing the coefficients of pressure correction | Noel | Phoenics | 1 | April 7, 2009 09:54 |
Does star cd takes reference pressure? | monica | Siemens | 1 | April 19, 2007 12:26 |
residual in the pressure correction | George | Main CFD Forum | 2 | July 28, 2005 05:43 |