Shooting method for PDE

조회 수: 22 (최근 30일)
Judah S
Judah S 2014년 11월 9일
댓글: Judah S 2014년 11월 10일
Hello all, There is a paper which I was going through (attached), where there are two equations which I would like to solve. Eq 16 and 17 based on 18 and 19 using shooting methods. Problem is the boundary conditions (infinity) and first order). I would be very thankful if you could help me solve the problem.
Thanking you.
Regards, Judah

답변 (1개)

Zoltán Csáti
Zoltán Csáti 2014년 11월 9일
During a shooting method you guess initial values from where you start solving the boundary value problem (BVP) as an initial value problem (IVP). You expect the result to be accurate when the right boundary condition (BC at infinity) is fulfilled. Since you solve an IVP, you can set the interval of integration. Choose it sufficiently large. If you use L1<L2 and there is little difference between them, then the truncated interval describes the problem of infinity well.
I solved similar boundary-layer-type equations several times. So I suggest you to try one or more of the followings:
  • Solve the BVP with bvp4c or bvp5c
  • Use Chebfun
  • Perform group analysis on them so that you can use a transformation method that can transform the BVP to an IVP (see e.g. Töpfer-transformation)
  • Use a series of transformations (see e.g. Iterative Transformation Method)
  • Code the shooting method (but be aware, this can be quite susceptible to the initial conditions)
  댓글 수: 3
Judah S
Judah S 2014년 11월 9일
I worked on the code just now and seem to approach in the right direction.
Here is my code (attached) Figure 2 seems to give similar result compared to the paper but Y axis values are not correct (due to boundary condition F'(0) = f(0) and F(infinity) = 0). I don't know how to implement these boundary conditions correctly. I would be very thankful if you could look into it.
Thank you.
Regards, Judah
Judah S
Judah S 2014년 11월 10일
Any help will be very appreciated.
Thanks,

댓글을 달려면 로그인하십시오.

카테고리

Help CenterFile Exchange에서 Boundary Conditions에 대해 자세히 알아보기

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by