How Does Changing Boundary Conditions Affect Nonlinear Shooting Methods in BVPs?

  • Thread starter Thread starter radiogaga35
  • Start date Start date
  • Tags Tags
    Nonlinear
radiogaga35
Messages
34
Reaction score
0
Given the boundary value problem (primes denote differentiation w.r.t x):
\begin{array}{l}<br /> y&#039;&#039; = f(x,y,y&#039;) \\ <br /> y(a) = \alpha \\ <br /> y(b) = \beta \\ <br /> \end{array}

the nonlinear shooting method may be implemented to solve the problem. A bisection algorithm may be used or, with a little more effort, Newton's method may be implemented (in which case one solves a fourth order, and not a second order, IVP - http://www.math.utah.edu/~pa/6620/shoot.pdf" ).

But what happens if the form of the boundary conditions changes to:
\begin{array}{l}<br /> y&#039;(a) = \alpha \\ <br /> y(b) = \beta \\ <br /> \end{array}

Is one still justified in reducing the BVP to a 2nd order IVP, but this time with initial SLOPE fixed at y&#039;(a) = \alpha and then trying different values of y(a) in order to achieve the condition y(b) = \beta? (As opposed to varying initial slope to achieve second condition).

Furthermore, if this will indeed work, then a bisection method should be easy to implement, but what about adapting Newton's method for this case? Can anyone point me to an appropriate reference that discusses this matter?

Thank you! :-)
 
Last edited by a moderator:
Physics news on Phys.org
Ok, I applied the aforementioned approach (with appropriately-adapted Newton-method implementation) to a trial problem and it worked perfectly.

In principle it seems like a sensible enough approach, but I'm not very clued up on BVP methods, so I'm not sure if there is any theoretical reason to avoid this approach? I.e. varying the initial "displacement" instead of slope. Certainly I've only ever seen the initial slope being varied in simple, single shooting methods
 
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