# something wrong

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

 March 3, 2001, 11:52 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.

 March 3, 2001, 15:32 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

 March 13, 2001, 10:19 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.

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post srihari FLUENT 1 October 31, 2016 15:18 pbryant OpenFOAM Bugs 1 February 14, 2011 22:17 Thomas Baumann OpenFOAM 6 December 21, 2009 11:31 Harry Qiu FLUENT 1 June 29, 2001 05:53 olivia FLUENT 1 June 23, 2001 17:06

All times are GMT -4. The time now is 12:55.