Abstract:In this article, firstly we develop a method for a type of difference equations, applicable to solve approximately a class of first order ordinary differential equation systems. In a second step, we apply the results obtained to solve a non-linear two point boundary value problem relating a flight mechanics model. We highlight the algorithm obtained seems to be robust and of easy computational implementation.
What problem does this paper attempt to address?
This paper attempts to solve the following problems:
1. **Develop a numerical method to solve a class of first - order nonlinear ordinary differential equation systems**:
- Firstly, the paper proposes a method suitable for solving a specific type of difference equations, and then uses it to approximately solve a class of first - order nonlinear ordinary differential equation systems. Specifically, the author approximates the solution of the nonlinear problem through a series of linearization steps and iterative algorithms.
2. **Apply the proposed method to solve the two - point boundary - value problem on one side in flight mechanics**:
- The author applies the above - mentioned method to a specific flight mechanics model, which describes the climbing motion of an aircraft on a plane. This is a nonlinear two - point boundary - value problem (TPBVP), which involves the changes of the aircraft's altitude \( h \), speed \( V \), flight path angle \( \gamma \) and horizontal position \( x \).
3. **Improve computational efficiency and implementation simplicity**:
- The paper emphasizes that the proposed algorithm is not only easy to implement, but also has efficient computational performance on standard commercial personal computers. For a nonlinear problem with about 40,000 degrees of freedom, the algorithm can complete the calculation in just a few seconds.
### Specific problem description
The specific equations of the flight mechanics model are as follows:
\[
\begin{cases}
\dot{h} = V\sin(\gamma), \\
\dot{\gamma} = \frac{1}{m_f V}(T\sin(e_3)+L)-\frac{g}{V\cos(\gamma)}, \\
\dot{V} = \frac{1}{m_f}(T\cos(e_3)-D)-g\sin(\gamma), \\
\dot{x} = V\cos(\gamma),
\end{cases}
\]
The boundary conditions are:
\[
\begin{cases}
h(0) = h_0, \\
V(0) = V_0, \\
x(0) = x_0, \\
h(t_f) = h_f,
\end{cases}
\]
where:
- \( h \) is the altitude of the aircraft,
- \( V \) is the speed of the aircraft,
- \( \gamma \) is the flight path angle,
- \( x \) is the horizontal position,
- \( t_f = 40 \) seconds is the total time,
- \( m_f \) is the mass of the aircraft,
- \( S_f \) is the wing area,
- \( a \) is the angle of attack,
- \( e_3 \) is the thrust angle,
- \( R_e \) is the radius of the earth,
- \( g \) is the gravitational acceleration,
- \( \rho_0 \) is the air density at sea level,
- \( B_9 \) is the atmospheric scale height,
- \( C_L \) and \( C_D \) are the lift coefficient and drag coefficient respectively,
- \( L \) is the lift,
- \( D \) is the drag,
- \( T \) is the thrust.
By introducing variable substitution:
\[
\begin{cases}
h = u_1, \\
\gamma = u_2, \\
V = u_3, \\
x = u_4,
\end{cases}
\]
The system can be expressed as:
\[
\begin{cases}
\dot{u_1} = f_1(u), \\
\dot{u_2} = f_2(u), \\
\dot{u_3} = f_3(u), \\
\dot{u_4} = f_4(u),
\end{cases}
\]
where:
\[
\begin{cases}
f_1(u) = u_3\sin(u_2), \\
f_2(u) = \frac{1}{m_f u_3}(T(u)\sin(e_3)+L(u))-\frac{g}{u_3\cos(u_2)}, \\
f_3(