New Reply

how to set up Neumann boundary condition for a PDE in a coordinate-invariant form?

 
Share Thread
Nov9-12, 12:46 PM   #1
 

how to set up Neumann boundary condition for a PDE in a coordinate-invariant form?


I'm having trouble finding out how to set up Neumann (or, rather, "Robin") boundary conditions for a diffusion-type PDE. More specifically, I have a scalar function [itex]f(\boldsymbol{x}, t)[/itex] where [itex]\boldsymbol{x}[/itex] is n-dimensional vector space with some boundary region defined by [itex]A(\boldsymbol{x})=0[/itex] where A is another scalar function. Then I want to specify a mixed (Robin-type) boundary condition in the following way:
[itex]a f(\boldsymbol{x}) + b \frac{∂f}{∂\boldsymbol{n}} = c(\boldsymbol{x})[/itex],
where a and b are some scalar numbers and c is a scalar function. What I don't understand is what is the nature of the directional derivative which is multiplied by [itex]b[/itex] in my example: from Wikipedia [ http://en.wikipedia.org/wiki/Neumann_boundary_condition ] it follows that
[itex]\frac{∂f}{∂\boldsymbol{n}} = ∇f \cdot \boldsymbol{n}[/itex],
where the first term, gradient of the scalar function, is a covariant vector, and the second, the normal vector to the boundary, also seems to be not a "true" vector but a covector (covariant vector) given by ∇A [ http://en.wikipedia.org/wiki/Normal_...ensional_space ]. So their product cannot be a scalar function that I need. What is wrong about it?
PhysOrg.com science news on PhysOrg.com

>> City-life changes blackbird personalities, study shows
>> Origins of 'The Hoff' crab revealed (w/ Video)
>> Older males make better fathers: Mature male beetles work harder, care less about female infidelity
Nov12-12, 10:58 AM   #2
 
It's not clear to me what is the problem you have. But perhaps this will help...

I'm never quite sure what people mean when they refer to the gradient, but in my book grad f is a vector field and df is a covector field (differential form). The coefficients of the gradient depend on a metric whereas the differential form components are simply the partial derivatives in whatever coordinate system you are using. With that convention the normal derivative is equal to df(n).

The vector n is a proper vector and is equal to the gradient of A divided by its own length. You need a metric tensor to define n. But you should already have one since you have a diffusion equation and the divergence operator depends on a metric.
New Reply

Similar discussions for: how to set up Neumann boundary condition for a PDE in a coordinate-invariant form?
Thread Forum Replies
heat equation with Neumann boundary condition Differential Equations 0
Neumann Condition on Curved Boundary using Finite Difference Differential Equations 0
Neumann Condition on Curved Boundary Mechanical Engineering 0
Physical meaning of Neumann boundary condition General Physics 2
Boundary-value problems (Neumann condition) Differential Equations 2