Poisson Equation Neumann boundaries singularity

  • #1
vector_problems
10
0
I am trying to solve the poisson equation with neumann BC's in a 2D cartesian geometry as part of a Navier-Stokes solver routine and was hoping for some help.

I am using a fast Fourier transform in the x direction and a finite difference scheme in the y. This means the poisson equation becomes

-kx^2p_{i,j}+(((p_{i,j-1})-(2*p_{i,j})+(p_{i,j+1}))/(dy^2))=RHS

with dp/dy=0 at the boundaries being enforced using ghost boundary points

this inverts easily for all wavenumbers except when kx=0 when the matrix is singular.

I was wondering if anyone had any experience dealing with this problem and any standard methods of solving such a problem?

Can you simply set p_{hat}=0 at these points and then convert back to real space?

Thanks
 

Answers and Replies

  • #2
18,931
9,194
Thanks for the post! This is an automated courtesy bump. Sorry you aren't generating responses at the moment. Do you have any further information, come to any new conclusions or is it possible to reword the post?
 
  • #3
Strum
105
4
Is p_{i,j} the Fourier transform? You should take care that the method does not produce some unwanted complex phase which might ruin the solution.
The kx=0 entry is the mean value for that "row" so you can not just set it equal to zero. I would imagine it depends a bit on your RHS how to best proceed.

PS: If you want to keep it fully spectral you could solve in the y-direction using a cosine transform.
 

Suggested for: Poisson Equation Neumann boundaries singularity

  • Last Post
Replies
3
Views
436
Replies
22
Views
550
Replies
3
Views
542
Replies
6
Views
958
Replies
1
Views
1K
Replies
8
Views
1K
Replies
6
Views
2K
Replies
1
Views
904
Replies
5
Views
1K
Top