I. INTRODUCTION

Air crash has many investigations[1]. This paper studies that air crash would be caused by deference of pressures between outside and inside of passenger's cargo. Such deference of pressures breaks the window or door and thus the plane crashes.

In order to calculate pressure outside the passenger's cargo, we study three type equations: A "wind - temperature equation" ("w-t eq." in short), B "wind - pressure /density equation" ("w - P/peq." in short). C "jet - pressure /density equation" ("J - P/peq." in short).

In Section 2, we derive the "w - T eq". It states that the derivative of wind speed respect to time proportions to the derivative of temperature respect to track (space).

In Section 3, the solution of "w - T eq." is obtained by method of separating variables. And it has been checked by weather forecasting, by dimensional check, and by other model check.

In Section 4, the "w - P/peq." is derived by the combination of Boyles law and Charles' law.

In Section 5, the solution of "w - P/peq." can be obtained directly from the solution of "w - T eq." in Section 3. In Section 5.1, the air density is calculated. Where the traditional method and engineering tool bar considered the air density is the functions of temperature and highness, but no connected to wind speed. Here, the air density is connected with wind speed.

In section 6, set up the "J - P/peq." The absolute motion (the jet plane motion relative a referenced point $\mathsf{P}$ on Earth) is equal to the linking motion (wind speed motion relative $\mathsf{P}$ ) and the relative motion (jet speed motion relative wind motion).

In Section 7, the solution of "J - P/peq." is directly obtained from section 3. Two cases: 1, density of jet plane is a constant; 2, the density of the jet plane is variable.

In Section 8, as an example of calculation of pressure outside the passenger's cargo, MH370 has been used. Where it had been studied in cases of nonpowered flying[6].

In Section 10, a conclusion is made.

II. WIND-TEMPERATURE EQUATION OF A POINT IN AIR ("W-T EQ." IN SHORT)

The derivation of "w - T equation" follows [3].

According to combination of Boyles' law and Charles' law (B-C law, in short), we have:

Where $p$ = the air pressure, accuracy. $p$ is a stress tensor, with the same normal stress and zero shearing stress

${\sigma }_{xx}={\sigma }_{yy}={\sigma }_{zz}=p=$ compressive normal stress; shearing stress ${\sigma }_{xy}={\sigma }_{yz}={\sigma }_{zx}=0.$ i, j, k are unit vectors in Cartesian co-ordinates. $V=$ volume of a point $=dxdydz=A_{x}dx=A_{y}dy=A_{z}dz,$ with $dx,dy,dz\to 0$ $A_x,A_y,A_z$ are the cross section areas of the element (point), respectively. c is a constant, T is the absolute temperature. (2-2) is the matrix form of $p$

Let operator $\mathrm{\nabla }$ be an inner product on both sides of (2-1), then, we have

Where

Where $F$ = the applied force, $S$ = the displacement

By Newton's second law, we have

Where $m$ is the mass, $u$ is the wind speed, $u=\partial u/\partial t$. Substituting (2-9) into (2-8), we have

(2-10) is called the "Wind -Temperature Equation of A Point (Mass) in Air", or ("W - T equation" in short). It states that $u$ is in proportion to $\mathrm{\nabla }T$, i.e., the derivative of wind speed $u$ respect to time $t$ proportions to the derivative of temperature $T$ respect to track (space).

III. SOLUTION OF "W - TEQ" BY METHOD OF SEPARATING VARIABLES

(2-10) is a vector PDE. It is hard to solve. Changing (2-10) to scalar function (its component form) is an easy way to solve. By (2-2) -- (2-4), we have:

Where $k=c/M$. $T_x^{\prime }=\frac{\partial }{\partial x}T$, $T_y^{\prime }=\frac{\partial }{\partial y}T$, $T_z^{\prime }=\frac{\partial }{\partial z}T$.

Looking at (3-3):

Where $T=T(z,t),u_z={\dot{u}}_z(z,t),\dot{u}=\partial u/\partial t.$

Suppose that variables of $k{T}_z^{\prime }$ and $u_z$ can be separated:

Where $a^{\prime }(z)=da(z)/dz$, $b(t)$, $e(z)$ and $f(t)=df(t)/dt$ are unknown functions and have continuous derivatives. We need sufficient relations to determine these unknowns.

The dimensions of $k{T}_z^{\prime }$ and ${\dot{u}}_z$ are the same (see the following red part), therefore, let (3-4) = (3-5), and separating variables, we have

Where the function of the left term is $z$, while the function of the middle term is $t$, if they are equal to each other, they must be a constant--- $k$.

For simple, we choose two relations:

then, integrating both sides of (3-6) for $t$ and by (3-8), we have:

Where $C_t$ is a constant, determined by initial condition, i.e., $t=0$, $\ln f(0)=1$, $\ln |1|=0$, $C_t=0$. Then, (3-9) gives:

Substituting (3-7) into (3-6), then, integrating both sides of (3-6) for $z$, we have:

Where $C_z$ is an integral constant, determined by the boundary condition, i.e., $z=0,T_z=T_z(0)=1{(}^{\circ }\mathrm{C})\to a(0)=1,\to C_z=0$. Therefore,

Substituting (3-14), (3-10), into (3-4), we have

Substituting (3-11), (3-13), into (3-5), we have

Similar to (3-16), for (3-1), (3-2), we have:

Expressed by vector form:

Dimensional check: (3-19)

Dimensional check is a tool often used to check the correctness of calculation. The dimensions of each term of an equality must be the same.

The dimensional check (red part) shows that the dimensions on both sides of the equality of (3-19) are the same. The dimensions in (3-19) is correct.

a) Checking the "w - T eq." by weather forecasting (China Weather Net), and motion equation of mushroom cloud.

Every winter, the weather forecast (China Weather Net) alerts people be wear that cold wave comes with strong wind accompanied with temperature sharp dropping. These description on temperature sharp dropping company with strong wind confirms that the difference of temperature between positions in track causes strong wind. Which agrees with (2-10) very well.

Moreover, a motion equation of mushroom cloud[3] derived by different model based on modifying of Navier-Stock equation, with result well agreed to (2-10). Which again confirms that (2-10) is believable.

IV. THE "WIND - PRESSURE/DENSITY EQUATION"(W-P/ρ EQ.IN SHORT)

Dividing both sides of B - C law(combination of Boyles law and Charles law) (2-1) by M (mass), we have

Where ${\rho }_{air}=M/V$ is the density of air. Substituting (4-1) into (2-10), we have:

(4-2) is called the "W-P/ρ equation". (4-2) shows that $u$ is in proportion to $\mathrm{\nabla }(p/{\rho }_{air})$, i.e., the derivative of wind speed $u$ respect to time $t$ proportions to the derivative of $(p/{\rho }_{air})$ respect to track (space).

V. SOLUTION OF "W - P/ρ EQ." (4-2)

The solution of scalar form of (4-2) can be obtained from (3-17), (3-18), and (3-16) by changing $k=c/m$ to $k={\rho }_{air}^{-1}$, i.e.,

The solution of vector form of (4-2) is:

The solution of vector form of (4-2) can be obtained from (3-19) by changing $k=c/m$ to $k={\rho }_{air}^{-1}$, i.e.,

a) Calculation of ${\rho }_{air}$

Traditionally, methods, tools for calculation of ${\rho }_{air}={\rho }_{air}(p,z)$ as functions of $p$ and $z^{[4,5]}$, based on Boyles law and Charles law for static description, i.e., they have no connection with wind speed. However, our treatment of ${\rho }_{air}$ is different to that of traditional. It connects with wind speed by "w - p/ $\rho$ equation". (4-2), based on Boyles law, Charles' law and together with the Newton's second law. The calculation of ${\rho }_{air}$, can be found in Appendix.

VI. SET UP THE "JET - PRESSURE/DENSITY EQUATION" ("J- P/ρ EQ." IN SHORT)

A jet plane flies in a atmospheric environment. In which the absolute motion (the jet plane motion relative a referenced point P on Earth) is equal to the linking motion (wind speed motion relative P) and the relative motion (jet speed motion relative wind motion).

Let the relative motion of the jet plane speed, density of air and pressure be $vv,p{p}_v$, and ${\rho }_{air}$. Then, set up the "J- $p/\rho$ Equation" similar to (4-2), we have:

Where ${\rho }_{jp}={\rho }_{jp}(z,t)=M(t)/V_{jp}=$ Density of jet plane; $M(t)$ is the mass of the jet plane. it varies with time t (jet fuel consuming); $V_{jp}$ is the volume of the jet plane.

(6-2) shows that $vv$ is in proportion to $\mathrm{\nabla }(p_{vv}/{\rho }_{jp})$ i.e., the derivative of jet plane speed $vv$ respect to time $t$, proportions to the derivative of pressure $(p_v/{\rho }_{jp})$ respect to track (space)of the jet plane.

VII. THE SOLUTION OF "J - P/ρ EQUATION" (6-2)

1. If${\rho }_{jp}=\mathrm{const}$, then the solution of(6-2)is the same as(5-5), just by replacing$v$,$p_v$instead of$u$,$p$, respectively. That is:

The solutions of similar (5-1) - (5-4) are:

1. If${\rho }_{jp}={\rho }_{jp}(z,t)=M(t)/V_{jp}$, is a known function. Where$M(t)$depends on flying status, e.g., flying at a constant speed, then

Where for $t=0$, $c_1=M(0)=$ jet plsne mass of full fuel.; for

$t=t_1,M(t_1)=C_1-t_1{c}_2=\mathrm{jet}$ plane mass of empty fuel

The solution of (6-2) is the same as (7-1) with ${\rho }_{jp}=M(t)/V_{jp}$ inside $\mathrm{\nabla }$, i.e.,