Please verify integral and approximation, boundary theory

Click For Summary
SUMMARY

This discussion focuses on the application of Newton's method and Taylor approximations to solve the equation $$f'''+\frac{m+1}{2}ff''+m(1-f^{'2})=0$$, which models the velocity of air over a flat plate. The user successfully derived the relationship for velocity as a function of distance from the plate, $$f'=f'(\eta)=\frac{u(y)}{u_e}$$, and is seeking verification of their integral calculation $$\int_0^\infty (1-f') dy$$. The user also outlines a MATLAB approximation method for this integral, correcting their earlier approach to ensure accurate computation of $$\int_A^B G(x) dx$$ using a summation technique.

PREREQUISITES
  • Understanding of fluid dynamics, specifically boundary layer theory.
  • Proficiency in numerical methods, particularly Newton's method and Taylor series.
  • Familiarity with MATLAB programming for numerical integration.
  • Knowledge of calculus, especially integral calculus and its applications in physics.
NEXT STEPS
  • Review the derivation of boundary layer equations in fluid dynamics.
  • Study numerical integration techniques in MATLAB, focusing on the trapezoidal and Simpson's rules.
  • Explore advanced topics in Taylor series expansions and their applications in fluid mechanics.
  • Investigate the implications of viscosity on boundary layer behavior in different fluid flow scenarios.
USEFUL FOR

Fluid dynamics researchers, mechanical engineers, MATLAB users, and anyone involved in computational fluid dynamics (CFD) who seeks to understand boundary layer behavior and numerical integration techniques.

fahraynk
Messages
185
Reaction score
5
I used Newtons method and taylor approximations to solve this equation $$f'''+\frac{m+1}{2}ff''+m(1-f^{'2})=0$$
It solves for velocity of air over a flat plate.
The velocity is a constant ##u_e## everywhere except in a boundary layer over the plate, where the velocity is a function of distance from the plate = u(y).
$$f'=f'(\eta)=\frac{u(y)}{u_e}\\\\\eta=y\sqrt{\frac{u_e}{\nu x}}$$
##\nu## is viscosity.

now that I have f', I need to calculate this integral
$$\int_0^\infty (1-f') dy $$

So my question which I would like help with is, can someone please tell me if the following work is correct :
$$\frac{dy}{d\eta}=\sqrt{\frac{\nu x}{u_e}}\\\\
\int_0^\infty (1-f') dy = \int_0^{\eta_{max}} (1-f') \frac{dy}{d\eta}d\eta = \sqrt{\frac{\nu x}{u_e}}\int_0^{\eta_{max}} (1-f') d\eta$$

##\eta_{max}## is the location where the boundary layer ends, and ##\frac{u}{u_e}=1## when ##\eta \longrightarrow \eta_{max}##, thus the integral becomes 1-1=0 at ##\eta=\eta_{max}##

Also, I am going to use this rule to approximate in MATLAB :
$$\int_A^B G(x) dx = (B-A)G(\frac{A+B}{2})\\\\
K= B-A = f'[2:N+1]-f'[1:N]\\\\
J = f'(\frac{A+B}{2}) = \frac{1}{2}(f'[1:N]+f'{2:N+1})\\\\
G(\frac{A+B}{2}) = ones(N)-J\\\\
\int_A^B G(x) dx = sum( K * (ones(N)-J) )$$
For this to work, I would have to input a specific x value.
 
Last edited:
  • Like
Likes   Reactions: berkeman
Physics news on Phys.org
I messed up K, and some other stuff, so to write over the integration :
instead of ##\int G(x)dx=(B-A)G(\frac{A+B}{2})## it should read
$$
\int_A^B G(x)dx
=\sum_{i=1}^{N}

{\frac{(B-A)}{N}G(\frac{A+B}{2})} \implies \int_0^{\infty}(1-f')dy =

sum[\frac{\eta_{max}}{N} * (ones(N)- \frac{1}{2}( f'[1:N]+f'[2:N+1] )]
$$
##\eta_{max}/N## represents B-A, the number of steps in the ##\eta## direction. Ones(N) is a vector length N of 1's. ##Ones(N) =[1,1,1...1_{N-1},1_N]##
 
Last edited:

Similar threads

  • · Replies 1 ·
Replies
1
Views
1K
  • · Replies 2 ·
Replies
2
Views
3K
  • · Replies 9 ·
Replies
9
Views
2K
Replies
4
Views
2K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 5 ·
Replies
5
Views
3K
  • · Replies 2 ·
Replies
2
Views
1K
  • · Replies 7 ·
Replies
7
Views
5K
  • · Replies 3 ·
Replies
3
Views
3K
  • · Replies 1 ·
Replies
1
Views
3K