Solving Navier-Stokes Eqn using Penalty Function Method

Thread View: sci.math.num-analysis

1 messages

[T] Back to Threads [V] Tree View [G] Group Articles

1 total messages Started by lzhao Thu, 28 Jan 1993 01:35

#1124

Author: lzhao
Date: Thu, 28 Jan 1993 01:35

25 lines
1240 bytes

Hi, 
I am solving Navier-Stokes Eqn. using FEM (Penalty-function method) right 
now. I have read lots of papers, but nobody mentioned how they handeled
the force vector. The force vector is

	F1=line or suface integration of(tx)
	F2=line or suface integration of(ty)
where
	tx= 2.0*mu*(du/dx)*nx+mu*((du/dy)+(dv/dx))*ny-p*nx
	ty= mu*((du/dy)+(dv/dx))*nx+2.0*mu*(du/dx)*nx-p*nx
and
   	p=-(penalty number)*(du/dx+dv/dy)

I am now using iteration method to handel this, ie first I assume a velocity
profile, secondly calculate the force vector and then sove for the velocity.
After this , I use the vel. got for the next iteration. But it seems that 
the result does not converge.

Is there anyone can give me some help? 

Regards.

lixin

email: lzhao@uoft02.utoledo.edu

Thread Navigation

This is a paginated view of messages in the thread with full content displayed inline.

Messages are displayed in chronological order, with the original post highlighted in green.

Use pagination controls to navigate through all messages in large threads.

Back to All Threads