CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Main CFD Forum

Analytical problem,finite element method (Q2Q1)

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

Reply
 
LinkBack Thread Tools Display Modes
Old   August 25, 2008, 12:42
Default Analytical problem,finite element method (Q2Q1)
  #1
Isi
Guest
 
Posts: n/a
Hello everybody,

I would be really gratefully if you could help me.(Sorry for the long email, but I try to do a good explanation of the problem)

I have to solve the incompressible Navier-Stokes equations in a square domain [0,1]x[0x1]. I calculate velocity and pressure fields with a decoupled scheme, using Q2Q1 pair of elements as spatial discretization, I am using finite element method. As temporal discretization I use Backward Euler Scheme, so itis a scheme of first order. In order to do an error analysis, I simulate an analytical case where the solution, both velocity and pressure,depends on time, something like "sin(a*pi*x)*cos(b*pi*y)*exp(-c*t)".

I use a grid of 2^4*2^4 elements, i.e. a spatial step=h=1/16, another one of 2^5*2^5 elements, i.e. h=1/32, another one of 2^6*2^6, i.e. h=1/64 and another grid of 2^7*2^7, i.e. h=1/128.

I use the pair element Q2Q1.

I use four different temporal steps: incrT=1/64,1/256,1/512 and 1/1024.

I think I have programmed everything in the right way. But the problem is that I don't obtain a good result. Well, the pressure error behaves in a good way since less refined the mesh, higher is the error. But, the velocity error doesn't follow this pattern, for example, I obtain:

Always with Q2Q1:

a)With incrT=1/64. (error with h=1/32) is the lowest. (error with h=1/128)>(error with h=1/64)>(error with h=1/32)

b)With incrT=1/256. (error with h=1/32) is the lowest. (error with h=1/128)>(error with h=1/64)>(error with h=1/32)

c)With incrT=1/512. (error with h=1/64) is the lowest (error with h=1/128)>(error with h=1/64)

d)With incrT=1/1024. (error with h=1/64) is the lowest (error with h=1/128)>(error with h=1/64)

Looking at the final solution and comparing it to the exact one, all the previous situations seem to reach a good agreement.

i) I have thought that perhaps I should try to implement a temporal discretization scheme of order two, such as Runge-Kutta 2, instead a Backward Euler, but I don't know if this would solve the problem.

ii) Besides, the most surprising thing is that the pressure (Q1 element) error has the properly behaviour.

Any advice will be welcome. Thanks in advance. Isa
  Reply With Quote

Old   August 26, 2008, 02:45
Default Re: Analytical problem,finite element method (Q2Q1
  #2
Rami
Guest
 
Posts: n/a
Isa,

I can't see how your analytic solution satisfy continuity. Maybe it is just because you have not explicitly written the expressions for u, v and p separately.

Rami
  Reply With Quote

Old   August 26, 2008, 03:23
Default Re: Analytical problem,finite element method (Q2Q1
  #3
Isi
Guest
 
Posts: n/a
Hello,

yes, you are right, I only wanted to show the general shape.

The analytical solutions are:

u=-cos(2*pi*x)sin(2*pi*y)exp(-8pi*pi*viscosity*t); v=cos(2*pi*y)sin(2*pi*x)exp(-8pi*pi*viscosity*t); p=-0.25*(cos(4*pi*x)+cos(4*pi*y))exp(-8pi*pi*viscosity*t);

Best Isa
  Reply With Quote

Old   August 26, 2008, 17:29
Default Re: Analytical problem,finite element method (Q2Q1
  #4
Vasilis
Guest
 
Posts: n/a
Hi Isi,

When you solve the Navier-Stokes eqns if you do not have a b.c. that sets your pressure, then you need to specify the pressure at a point in your domain. You could impose for example that P(0)=1.

I would not make any difference if you use the Runge-Kutta method. Actually, the Backward Euler method is a very stable method.

Finally, try to make the eqns dimensionless before you solve them. In this way, the Reynolds number will appear, and then you could find the value of the Reynolds number that causes this behavior.

  Reply With Quote

Old   August 27, 2008, 07:06
Default Re: Analytical problem,finite element method (Q2Q1
  #5
Isi
Guest
 
Posts: n/a
Hi!

I put b.c. on the velocity variable in all the walls.

Backward Euler method is stable but is only first order accurate, it's for that I had thought in implementing Runge-Kutta, because it is of second order.

The Reynolds number is 100.

Best Isi
  Reply With Quote

Old   August 27, 2008, 13:03
Default Re: Analytical problem,finite element method (Q2Q1
  #6
Vasilis
Guest
 
Posts: n/a
Hi Isi,

Is not enough to specify the velocity in all the boundaries! When you solve the Navier-Stokes eqns and you do not specify the pressure in one of you b.c you need to specify it in your domain.

Consider the simplest case of two parallel plates, where the top plate is moving with constant V and the bottom plate is motionless. The solution for the velocity profile is obvious, but does the Navier-Stokes eqns give you any information about the pressure? The answer is no, and the only thing that you learn from the solution of these eqns is that the pressure is constant, and nothing but that! Because of that you need to specify the pressure at a given point in your domain. It won't affect your solution!

Yes you are right, backward Euler is a first order method, but its order will affect the rate at which the error is changing.

A rule of thumb when you solve a multidimensional problem is to achieve convergence in the spatial domain, and then to achieve convergence in the time domain. If your solution has not converged in the spatial domain, do not solve the eqns for a different value of dt.

  Reply With Quote

Reply

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
Lattice Boltzmann method vs Finite Element Method and Finite Volume Method solemnpriest Main CFD Forum 3 August 12, 2013 11:00
Control Volume Finite Element Method gerardosrez Main CFD Forum 0 March 16, 2011 15:49
Is CFD coding in finite element method a formidable job? bearcat Main CFD Forum 0 June 18, 2010 20:53
Boundary Element Method = Boundary integral method? aiqch Main CFD Forum 0 April 3, 2010 16:29
What is C.V. based finite element method C-H Kuo Main CFD Forum 3 November 5, 1998 10:07


All times are GMT -4. The time now is 19:34.