A General natural boundary condition?

maistral
Messages
235
Reaction score
17
TL;DR Summary
A little clarification with regard to the natural boundary condition.
Hi, I'd like to be clarified regarding the general natural/Neumann boundary condition for a PDE.

1. The natural boundary condition is generally defined as:
242655
(1)
and can be expressed as, according to this resource:
242656
(2)

But apparently, according to https://www.researchgate.net/post/How_to_impose_natural_boundary_conditions_with_Generalized_Finite_Difference_Method_or_meshfree_collocation_method resource (posted by Dr. Fan):
242657
(3)

Which is which? Is it supposed to be positive, or negative? When should it be positive or negative?

2. If I apply the derivative boundary condition, on say, the bottom of a square plate, I can state:
242658
(4)
or for simplicity,
242659
(5)

Obviously, comparing it with (2), ny is equal to 0. Why is this so? Why is ny = 0?

Thanks!
 
Physics news on Phys.org
1. The minus sign is an obvious typo.
2. ##\vec n## is the normal vector. The normal vector's y-component on a coordinate surface of ##y## is zero.
 
  • Like
Likes maistral
Hi! Thanks for replying. Thank you very much, it did clarify a lot of things. I still have two remaining problems, however.

1. Where is the vector n exactly pointed to? Is it going towards the region, or is it going away from the region?

2. I am assuming that the complete form of the Neumann condition is used for irregular boundaries (ie. curved ones, so on and forth). As stated by the first resource, nx = cos θ, and ny = sin θ.

If I use finite element analysis, I can just take the normal vector perpendicular to the edge of the element which is definitely a no-brainer as it's a straight line.

My problem begins to appear if I begin using finite differences. How do I exactly implement this kind of boundary condition on nodes? I'm confused - nodes are, well, nodes; just a point. I'm under the impression that I can get the angle required for calculating nx and ny by taking the tangent of the irregular boundary curve at the node, then I draw the perpendicular line needed to calculate the angle. Is this even correct?
 
There is the following linear Volterra equation of the second kind $$ y(x)+\int_{0}^{x} K(x-s) y(s)\,{\rm d}s = 1 $$ with kernel $$ K(x-s) = 1 - 4 \sum_{n=1}^{\infty} \dfrac{1}{\lambda_n^2} e^{-\beta \lambda_n^2 (x-s)} $$ where $y(0)=1$, $\beta>0$ and $\lambda_n$ is the $n$-th positive root of the equation $J_0(x)=0$ (here $n$ is a natural number that numbers these positive roots in the order of increasing their values), $J_0(x)$ is the Bessel function of the first kind of zero order. I...
Are there any good visualization tutorials, written or video, that show graphically how separation of variables works? I particularly have the time-independent Schrodinger Equation in mind. There are hundreds of demonstrations out there which essentially distill to copies of one another. However I am trying to visualize in my mind how this process looks graphically - for example plotting t on one axis and x on the other for f(x,t). I have seen other good visual representations of...
Back
Top