When can grad(u)=(b-u)N be solved, and how?

  • Thread starter Thread starter syra
  • Start date Start date
Join the discussion
Registration is free. Ask a follow-up in this thread, or start your own.
4 replies · 1K views
syra
Messages
5
Reaction score
0
b(x) is a scalar-valued function, and N(x) a vector valued-function, taking on arguments in R^n. Not actually for schoolwork, but academic nonetheless.
 
Physics news on Phys.org
If we can ignore singular behavior ([tex]\vec{N}\neq 0[/tex], [tex]b(x)-u(x) \neq 0[/tex], etc), one way to try to solve

[tex]\nabla u = (b-u)\vec{N}[/tex]

is the following. Define a new orthonormal basis

[tex]\vec{v} = \frac{1}{\sum_i N_i^2} (N_1, \ldots, N_n),[/tex]

[tex]\vec{v}^{(1)}_\perp = \frac{1}{N_1^2+N_2^2} (N_2, -N_1,0,\ldots,0).[/tex]

The other [tex]\vec{v}^{(i)}_\perp[/tex] can be constructed by Gram-Schmidt. So these are all an orthonormal set.

We then note that

[tex]\vec{v} \cdot \nabla u = (b-u)[/tex]

[tex]\vec{v}^{(i)}_\perp \cdot \nabla u = 0.[/tex]

If we can define new coordinates [tex](X,Y_i)[/tex] by inverting

[tex]\frac{\partial}{\partial X} = \vec{v} \cdot \nabla_x,~~\frac{\partial}{\partial Y_i} = \vec{v}^{(i)}_\perp \cdot \nabla_x,[/tex]

then we obtain a first-order ordinary differential equation for [tex]u(X)[/tex]:

[tex]\frac{\partial u}{\partial X} = b- u, ~~ \frac{\partial u}{\partial Y^i} =0.[/tex]

I suspect that one of the requirements for these new coordinates to exist is that [tex]N = \nabla\phi[/tex] for some scalar [tex]\phi[/tex] so that 2nd derivatives match.
 
Excuse my laziness for not using latex.

The coordinate transformation doesn't require an orthonormal set of v's, it just requires that ||v||=1/||N|| and orthogonality. The new coordinates require grad(X) = v = N/||N||^2, but we can lax this condition up by not requiring ||v||=1, in which case the final equations will have du/dX = (b-u) ||v||/||N|| instead, and there need only exist some scalar function phi such that grad(X) = phi N for some X.

In any case, the resulting differential equation is not in fact ordinary, and I still don't know how to solve it...
 
If [tex]\partial u/\partial Y^i =0[/tex], then [tex]u = u(X)[/tex], so the remaining equation is an ODE. The equations that determine the coordinate [tex]X[/tex] are still PDEs, but are likely to be simpler than the original equation for u.
 
Last edited:
Oh woops, I was being silly.