Return to computing page for the second course APMA0340
Return to Mathematica tutorial for the second course APMA0340
Return to the main page for the course APMA0330
Return to the main page for the course APMA0340
Return to Part III of the course APMA0330
Glossary
Heat Transfer in Boundary Layer
The algorithm for solving boundary value problems for differential equations on infinite or semi-infinite intervals include four steps. The first step is to compute high-order Taylor series expansions using an algebraic manipulation language such as Maple or Mathematica. These expansions will contain one or more unknown parameters z which will be determined by the boundary condition at infinity. The second step is to convert the Taylor expansions into diagonal Padé approximants. The boundary condition that u(x) decays to zero at infinity becomes the condition that the coefficient of the highest power of x in the numerator polynomial must be zero. The third step is to solve this equation for the free parameter z. The final step is to evaluate each of the multiple solutions of this equation for physical plausibility and convergence (as N increases). Methods for nonlinear problems are almost always iterative and need a first guess to initialize the iteration. The Padé algorithm is unusual in that it is a direct method that requires no a priori information about the solution.
Heat Transfer in Boundary Layer
We consider the heat transfer problem in steady two-dimensional boundary layer stagnation-point flow of viscous incompressible fluid toawrds a stretching/shrinking sheet in the presence of nonuniform heat flux. The governing equations of motion and the nergy equations may be written in usual notations as follows (see Krishnendu Bhattacharyya, Heat transfer in boundary layer stagnation-point flow towards a shrinking sheet with non-uniform heat flux, Chinese Physics B, 22, 2013, 074705):
The appropriate boundary conditions for the velocity components and temperature are as follows:
Upon introducing teh similarity transformations
We convert the given problem for f to a regular boundary value problem by choosing a large number L so that get
To derive the recurrence, we consider generating functions for f and for nonlinear operator:
DSolve[{u'''[x] + 1 u[0][x]*D[u[0][x], x, x] - (D[u[0][x], x])^2 == 0, u[0] == 0, u'[0] == B, u'[L] == 0}, u[x], x]
3780 B^2 L^3 x^2 + 840 L^4 x^2 - 252 L^5 x^2 - 588 B L^5 x^2 - 49 L^7 x^2 - 3360 L x^3 + 3360 B L x^3 + 3360 B^2 L^2 x^3 -
1680 L^3 x^3 + 840 B L x^4 - 840 B^2 L x^4 + 840 L^2 x^4 + 420 B L^3 x^4 + 168 x^5 - 336 B x^5 + 168 B^2 x^5 +
168 L^2 x^5 - 168 B L^2 x^5 + 42 L^4 x^5 - 56 L x^6 + 56 B L x^6 - 28 L^3 x^6 + 8 L^2 x^7)
- John P. Boyd, Padé approximant algorithm for solving nonlinear ordinary differential equation boundary value problems on an unbounded domain, Computers in Physics, 1997, Vol. 11, 299 (1997); https://doi.org/10.1063/1.168606
- Randolph Rach, Abdul‐Majid Wazwaz, Jun‐Sheng Duan, A reliable modification of the Adomian decomposition method for higher‐order nonlinear differential equations, Kybernetes, Vol. 42 Issue: 2, 2013, pp.282--308, https://doi.org/10.1108/03684921311310611
- Jun-Sheng Duan, Randolph Rach, Abdul-Majid Wazwaz, Temuer Chaolu, Zhong Wang, A new modified Adomian decomposition method and its multistage form for solving nonlinear boundary value problems with Robin boundary conditions, Applied Mathematical Modelling, 37, 2013, 8687--8708.
- Abdul-Majid Wazwaz, Randolph Rach, Jun-Sheng Duan, Adomian decomposition method for solving the Volterra integral form of the Lane–Emden equations with initial values and boundary conditions, Applied Mathematics and Computation, 219, 2013, 5004--5019.
- Holmquist, Sonia, An examination of the effectiveness of the Adomian Decomposition Method in fluid dynamic applications, University of Central Florida, Theses and Dissertations, Doctoral Dissertation (Open Access), 2007.
Return to Mathematica page
Return to the main page (APMA0330)
Return to the Part 1 (Plotting)
Return to the Part 2 (First Order ODEs)
Return to the Part 3 (Numerical Methods)
Return to the Part 4 (Second and Higher Order ODEs)
Return to the Part 5 (Series and Recurrences)
Return to the Part 6 (Laplace Transform)
Return to the Part 7 (Boundary Value Problems)