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

something wrong

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 3, 2001, 11:52
Default something wrong
  #1
JunseokKim
Guest
 
Posts: n/a
Hi, people.

I tested a poisson solver in 3D with an exact solution,

first of all, I tested laplacian with f(x,y,z) = cos(pi x) cos(pi y) cos(pi z), and I used standard 7 points finite difference laplacian, and as a source term I used exact solution -3pi^2 cos(pi x) cos(pi y) cos(pi z) and I calculate the residual,

but it is huge error. something like 29.00 even in 32 x 32 grid points, and my domain is unit cubic, and I imposed a Neumann boundary conditions.

Do you have any idea about this situation? How should I implement at boundary to reduce this big errors.
  Reply With Quote

Old   March 3, 2001, 15:32
Default Re: something wrong
  #2
Adrin Gharakhani
Guest
 
Posts: n/a
Why don't you start with the simpler Dirichlet BC, where you specify the BC exactly, say zero all around, without the approximations that Neumann brings into the table? So long as you have a source term this BC will be ok.

If you get a better convergence rate with DBC, then perhaps you need to use higher order approximations for the grads in the NBC or at least cluster the grids near the boundaries. If you still get bad results, then you have a bug!

Adrin Gharakhani
  Reply With Quote

Old   March 13, 2001, 10:19
Default Re: something wrong
  #3
yfyap
Guest
 
Posts: n/a
dear JunseokKim, i happened to solve a poisson eq. with neumann bc. please check if your poisson eq. if it satisfy the consequence of the Green theorem, failing to satisfy this theorem will cause problem in the numerical solution. if your eq. doesn't satisfy this theorem, you can modify the source term. please refer to the posting: convergent prob. in solving poisson eq(neumann BC) (118) - yf yap , Thu, 8 Feb 2001, 8:29 a.m. in this forum for more information on the matter. hope that it is helpful.

  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
Axi-symmetric case gives wrong answer pbryant OpenFOAM Bugs 1 February 14, 2011 22:17
BuoyantBoussinesqSimpleFoam and axial-symmetric results wrong mass flow Thomas Baumann OpenFOAM 6 December 21, 2009 11:31
udf error srihari FLUENT 0 February 9, 2009 10:00
Pressure contour seems wrong??? Harry Qiu FLUENT 1 June 29, 2001 05:53
What's wrong with my UDF? olivia FLUENT 1 June 23, 2001 17:06


All times are GMT -4. The time now is 21:22.