Hi everybody,

I'm trying to find the correct discretisation with finite differences of a neumann boundary condition at a corner node.

As an example:

dpdx = f1(x,y) and dpdy = f2(x,y)

So how to apply them in a square cavity for the node (0,0)

It's obvious that we can't apply only one of them, nor both...

If you have any idea or a reference that gives a clue, thanks in advance.

If you want to be exact, there are solutions based on perturbation methods which account for the exact corner node. For other problems, the application of BC depends on how you discretize. Regarding the cavity problem, there was a paper long ago by goldstein ( may be wrong name ? ) which discusses the issue of applying the boundary condition at the corner node.

Thanks for your answer, but could you give me some references about the use of the perturbation methods in cfd... Also, I found no paper related to Goldstein on the cavity problem...

Cheers