CFD Online Logo CFD Online URL
Home > Forums > Main CFD Forum

something wrong

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

LinkBack Thread Tools Display Modes
Old   March 3, 2001, 11:52
Default something wrong
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
Adrin Gharakhani
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
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


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
udf error srihari FLUENT 1 October 31, 2016 15:18
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
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 12:55.