home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1993 #3 / NN_1993_3.iso / spool / sci / math / numanal / 3990 < prev    next >
Encoding:
Internet Message Format  |  1993-01-28  |  1.2 KB
  1. Path: sparky!uunet!charon.amdahl.com!amdahl!rtech!sgiblab!spool.mu.edu!agate!usenet.ins.cwru.edu!magnus.acs.ohio-state.edu!bgsuvax!uoft02.utoledo.edu!lzhao
  2. From: lzhao@uoft02.utoledo.edu
  3. Newsgroups: sci.math.num-analysis
  4. Subject: Solving Navier-Stokes Eqn using Penalty Function Method
  5. Message-ID: <1993Jan27.203547.1323@uoft02.utoledo.edu>
  6. Date: 27 Jan 93 20:35:47 EST
  7. Organization: University of Toledo, Computer Services
  8. Lines: 25
  9.  
  10. Hi, 
  11. I am solving Navier-Stokes Eqn. using FEM (Penalty-function method) right 
  12. now. I have read lots of papers, but nobody mentioned how they handeled
  13. the force vector. The force vector is
  14.  
  15.     F1=line or suface integration of(tx)
  16.     F2=line or suface integration of(ty)
  17. where
  18.     tx= 2.0*mu*(du/dx)*nx+mu*((du/dy)+(dv/dx))*ny-p*nx
  19.     ty= mu*((du/dy)+(dv/dx))*nx+2.0*mu*(du/dx)*nx-p*nx
  20. and
  21.        p=-(penalty number)*(du/dx+dv/dy)
  22.  
  23. I am now using iteration method to handel this, ie first I assume a velocity
  24. profile, secondly calculate the force vector and then sove for the velocity.
  25. After this , I use the vel. got for the next iteration. But it seems that 
  26. the result does not converge.
  27.  
  28. Is there anyone can give me some help? 
  29.  
  30. Regards.
  31.  
  32. lixin
  33.  
  34. email: lzhao@uoft02.utoledo.edu    
  35.