Many physical situations are modelled by systems of \(n\) differential equations in \(n\) unknown functions, where \(n\ge2\). The next three examples illustrate physical problems that lead to systems of differential equations. In these examples and throughout this chapter we’ll denote the independent variable by \(t\).
Example 10.1.1
: Flow
Tanks \(T_1\) and \(T_2\) contain 100 gallons and 300 gallons of salt solutions, respectively. Salt solutions are simultaneously added to both tanks from external sources, pumped from each tank to the other, and drained from both tanks (Figure 10.1.1
). A solution with \(1\) pound of salt per gallon is pumped into \(T_1\) from an external source at \(5\) gal/min, and a solution with \(2\) pounds of salt per gallon is pumped into \(T_2\) from an external source at \(4\) gal/min. The solution from \(T_1\) is pumped into \(T_2\) at 2 gal/min, and the solution from \(T_2\) is pumped into \(T_1\) at \(3\) gal/min. \(T_1\) is drained at \(6\) gal/min and \(T_2\) is drained at 3 gal/min. Let \(Q_1(t)\) and \(Q_2(t)\) be the number of pounds of salt in \(T_1\) and \(T_2\), respectively, at time \(t>0\). Derive a system of differential equations for \(Q_1\) and \(Q_2\). Assume that both mixtures are well stirred.
As in Section 4.2, let rate in and rate out denote the rates (lb/min) at which salt enters and leaves a tank; thus,
A mass \(m_1\) is suspended from a rigid support on a spring \(S_1\) and a second mass \(m_2\) is suspended from the first on a spring \(S_2\) (Figure 10.1.2
). The springs obey Hooke’s law, with spring constants \(k_1\) and \(k_2\). Internal friction causes the springs to exert damping forces proportional to the rates of change of their lengths, with damping constants \(c_1\) and \(c_2\). Let \(y_1=y_1(t)\) and \(y_2=y_2(t)\) be the displacements of the two masses from their equilibrium positions at time \(t\), measured positive upward. Derive a system of differential equations for \(y_1\) and \(y_2\), assuming that the masses of the springs are negligible and that vertical external forces \(F_1\) and \(F_2\) also act on the objects.
Solution
In equilibrium, \(S_1\) supports both \(m_1\) and \(m_2\) and \(S_2\) supports only \(m_2\). Therefore, if \(\Delta\ell_1\) and \(\Delta\ell_2\) are the elongations of the springs in equilibrium then
\[\label{eq:10.1.7} (m_1+m_2)g=k_1\Delta\ell_1\quad\mbox{ and }\quad m_2g=k_2\Delta\ell_2. \]
Figure 10.1.2
Let \(H_1\) be the Hooke’s law force acting on \(m_1\), and let \(D_1\) be the damping force on \(m_1\). Similarly, let \(H_2\) and \(D_2\) be the Hooke’s law and damping forces acting on \(m_2\). According to Newton’s second law of motion,
When the displacements are \(y_1\) and \(y_2\), the change in length of \(S_1\) is \(-y_1+\Delta\ell_1\) and the change in length of \(S_2\) is \(-y_2+y_1+\Delta\ell_2\). Both springs exert Hooke’s law forces on \(m_1\), while only \(S_2\) exerts a Hooke’s law force on \(m_2\). These forces are in directions that tend to restore the springs to their natural lengths. Therefore
\[\label{eq:10.1.9} H_1=k_1(-y_1+\Delta\ell_1)-k_2(-y_2+y_1+\Delta\ell_2)\quad\mbox{ and }\quad H_2=k_2(-y_2+y_1+\Delta\ell_2). \]
When the velocities are \(y_1'\) and \(y_2'\), \(S_1\) and \(S_2\) are changing length at the rates \(-y_1'\) and \(-y_2'+y_1'\), respectively. Both springs exert damping forces on \(m_1\), while only \(S_2\) exerts a damping force on \(m_2\). Since the force due to damping exerted by a spring is proportional to the rate of change of length of the spring and in a direction that opposes the change, it follows that
\[\label{eq:10.1.10} D_1=-c_1y_1'+c_2(y_2'-y_1')\quad\mbox{ and }\quad D_2=-c_2(y_2'-y_1'). \]
From Equations \ref{eq:10.1.8}, \ref{eq:10.1.9}, and \ref{eq:10.1.10},
Let \({\bf X}={\bf X}(t)=x(t)\,{\bf i}+y(t)\,{\bf j}+z(t)\,{\bf k}\) be the position vector at time \(t\) of an object with mass \(m\), relative to a rectangular coordinate system with origin at Earth’s center (Figure 10.1.3
).
Figure 10.1.3
According to Newton’s law of gravitation, Earth’s gravitational force \({\bf F}={\bf F}(x,y,z)\) on the object is inversely proportional to the square of the distance of the object from Earth’s center, and directed toward the center; thus,
Cancelling the common factor \(m\) and equating components on the two sides of this equation yields the system\[\label{eq:10.1.14} \begin{array}{ll} {x''}&{= - \dfrac{ gR^2x}{(x^2+y^2+z^2)^{3/2}}} \\[4pt] {y''} &{= - \dfrac{gR^2y}{(x^2+y^2+z^2)^{3/2}}} \\[4pt] {z''}&{= - \dfrac{gR^2z}{(x^2+y^2+z^2)^{3/2}}.} \end{array} \]
Rewriting Higher Order Systems as First Order Systems
is called a first order system, since the only derivatives occurring in it are first derivatives. The derivative of each of the unknowns may depend upon the independent variable and all the unknowns, but not on the derivatives of other unknowns. When we wish to emphasize the number of unknown functions in Equation \ref{eq:10.1.15} we will say that Equation \ref{eq:10.1.15} is an \(n\times n\) system.
Systems involving higher order derivatives can often be reformulated as first order systems by introducing additional unknowns. The next two examples illustrate this.
The difference in form between Equation \ref{eq:10.1.15} and Equation \ref{eq:10.1.17}, due to the way in which the unknowns are denoted in the two systems, isn’t important; Equation \ref{eq:10.1.17} is a first order system, in that each equation in Equation \ref{eq:10.1.17} expresses the first derivative of one of the unknown functions in a way that does not involve derivatives of any of the other unknowns.
Rewriting Scalar Differential Equations as Systems
In this chapter we’ll refer to differential equations involving only one unknown function as scalar differential equations. Scalar differential equations can be rewritten as systems of first order equations by the method illustrated in the next two examples.
Since systems of differential equations involving higher derivatives can be rewritten as first order systems by the method used in Examples 10.1.5
-10.1.7
, we’ll consider only first order systems.
Numerical Solution of Systems
The numerical methods that we studied in Chapter 3 can be extended to systems, and most differential equation software packages include programs to solve systems of equations. We will not go into detail on numerical methods for systems; however, for illustrative purposes we’ll describe the Runge-Kutta method for the numerical solution of the initial value problem
at equally spaced points \(t_0\), \(t_1\), …, \(t_n=b\) in an interval \([t_0,b]\). Thus,
\[t_i=t_0+ih,\quad i=0,1,\dots,n, \nonumber \]
where
\[h={b-t_0\over n}. \nonumber \]
We’ll denote the approximate values of \(y_1\) and \(y_2\) at these points by \(y_{10},y_{11},\dots,y_{1n}\) and \(y_{20},y_{21},\dots,y_{2n}\). The Runge-Kutta method computes these approximate values as follows: given \(y_{1i}\) and \(y_{2i}\), compute
for \(i=0\), …, \(n-1\). Under appropriate conditions on \(g_1\) and \(g_2\), it can be shown that the global truncation error for the Runge-Kutta method is \(O(h^4)\), as in the scalar case considered in Section 3.3.