Applying shooting method to 4nd two point value boundary problem

In summary, you can solve a 4nd differential eq. using shooting method by iterating on the two unknown parameters. You need to make good on the two boundary conditions at the far end and Newton-Raphson is used to solve the equations.
  • #1
pgioun
5
0
Hi,
I was wondering if you can apply shooting method to a 4nd differential eq.
two point value boundary problem,

specifically I want to use this method to solve Euler-Bernoulli eq.

EI y(4)(x)=f(x), y(0)=0,y'(0)=0,y(L)=0,y'(L)=0.

Normally, if you have a 2nd order two point value boundary problem,
you can solve it with shooting method by making it an initial condition problem,
''guessing'' the one of the two initial conditions, and iterating (changing your guess) until you meat the 2nd boundary condition.

But when you have a 4nd differentila eq. you should ''guess'' 2 initial conditions... so, how this would work?
Thanks
 
Physics news on Phys.org
  • #2
You have to iterate on the two unknown parameters, insteead of just one.

How to solve this type of "optimisation" problem with many unknowns is a (large) topic in numerical analysiss, but a simple method that doesn't need any mathematical analysis of the particular problem you want to solve is http://en.wikipedia.org/wiki/Nelder–Mead_method

The basic idea is similar to searching along a line. Instead of finding an interval that you includes the solution and then making it smaller, you use the two unknown parameters as coordinates to define points on a plane, find a triangular region that contains the solution, and then choose a new point on the plane for each iteration that will shrink the area of the triangle.

There are much more efficient methods (FWIW my choice of "go-to" method would be the BFGS algorithm - google for it), but whether it's worth learning about them depends on the size of the problem. It's certainly worth reducing the computing time to solve a problem from say 100 hours to 1 hour, but reducing it from 10 sec to 0.1 sec might not be worth doing, unless it takes you less than 10 seconds to write and debug the code!
 
  • #3
pgioun said:
Hi,
I was wondering if you can apply shooting method to a 4nd differential eq.
two point value boundary problem,

specifically I want to use this method to solve Euler-Bernoulli eq.

EI y(4)(x)=f(x), y(0)=0,y'(0)=0,y(L)=0,y'(L)=0.

Normally, if you have a 2nd order two point value boundary problem,
you can solve it with shooting method by making it an initial condition problem,
''guessing'' the one of the two initial conditions, and iterating (changing your guess) until you meat the 2nd boundary condition.

But when you have a 4nd differentila eq. you should ''guess'' 2 initial conditions... so, how this would work?
Thanks

You want to make good on the two boundary conditions at the far end. The error in satisfying each of these conditions is a function of the two guesses at the other end:

ε1= ε1 (g1 , g2)

ε2= ε2 (g1 , g2)

You want g1 and g2 to be such that ε1 and ε2 are zero. This can be regarded as two non-linear equations in two unknowns. You can solve these equations using Newton-Raphson. You need to integrate from one boundary to the other using the two present guesses and get the corresponding values of the two errors. Then you need to evaluate the partial derivatives of the errors with respect to the two guesses. This is done by perturbing each of the guesses at the near boundary and then integrating to evaluate numerically the changes in the errors at the other boundary. You then apply the Newton-Raphson formula, which involves solving two linear algebraic equations in two unknowns for the changes in the guesses to make the errors come out to zero. This procedure is repeated until the errors are close enough to zero.
 

FAQ: Applying shooting method to 4nd two point value boundary problem

1. How does the shooting method work for solving a 4nd two point value boundary problem?

The shooting method is a numerical technique used to solve boundary value problems by converting them into initial value problems. It involves guessing an initial value for the unknown boundary condition and using an appropriate numerical method to solve the resulting initial value problem. This process is repeated until the desired boundary condition is met.

2. What are the advantages of using the shooting method for solving boundary value problems?

The shooting method is a relatively simple and efficient way to solve boundary value problems, especially when compared to other numerical methods such as finite difference or finite element methods. It also allows for more flexibility in choosing the boundary conditions and can be easily implemented using standard numerical software packages.

3. Are there any limitations to using the shooting method for solving boundary value problems?

The shooting method may not always converge to the desired solution, especially if the initial guess for the boundary condition is not close enough. It also requires a good understanding of the underlying mathematical model and the choice of numerical method used to solve the initial value problem.

4. Can the shooting method be used for solving higher order boundary value problems?

Yes, the shooting method can be extended to solve higher order boundary value problems by converting them into systems of first order differential equations. This allows for the use of standard numerical methods such as Euler's method or Runge-Kutta methods to solve the resulting system.

5. How do you determine the convergence of the shooting method?

The convergence of the shooting method can be determined by comparing the results obtained from different initial guesses for the boundary condition. If the results are consistent and approach the desired boundary condition, then the method is considered to have converged. Additionally, the convergence can also be checked by varying the step size used in the numerical method and observing the change in the solution.

Back
Top