Equations différentielles : en route vers le supérieur

Un autre modèle physique - Exercice 1

45 min
70
Dans cet exercice, on aborde la résolution de l'équation différentielle non-linéaire issue de la modélisation d'une chute verticale dans un champ de pesanteur uniforme, avec des frottements aérodynamiques, c'est-à-dire proportionnel au carré de la norme du vecteur vitesse.
Question 1
Soit mm, kk et gg trois constantes réelles strictement positives.
On considère l'équation différentielle non linéaire (E)(E) suivante :
(E):mv(t)+kv2(t)=mg(E) : \,\, mv'(t) + kv^2(t) = mg
la fonction vv représente la vitesse.

Démontrer l'existence physique d'une vitesse limite vv_\ell. Préciser son expression.

Correction
Si la vitesse atteint une vitesse limite vv_\ell alors elle ne change plus de valeur et dans ce cas elle est constante. Ainsi sa dérivée vv'_\ell est nulle. Ainsi l'équation (E)(E) devient :
mv(t)+kv2=mgm v_\ell'(t) + kv_\ell^2 = mg
Soit :
m×0+kv2=mgm \times 0 + k v_\ell^2 = mg
Donc :
v2=mgkv=±mgkv_\ell^2 = \dfrac{mg}{k} \,\,\,\, \Longrightarrow \,\,\,\, v_\ell = \pm \sqrt{\dfrac{mg}{k}}
Or v>0v_\ell > 0, finalement on en déduit :
v=mgk\red{\boxed{v_\ell = \sqrt{\dfrac{mg}{k}}}}
Question 2

Montrer que l'équation différentielle (E)(E) prend la forme suivante : v(t)=G(1(v(t)v)2)v'(t) = G \left( 1 - \left( \dfrac{v(t)}{v_\ell} \right)^2 \right)GG est une constante donc vous préciserez l'expression.

Correction
On a :
mv(t)+kv2(t)=mgmv'(t) + kv^2(t) = mg
Mais v=mgkv_\ell = \sqrt{\dfrac{mg}{k}} ce qui implique que mg=kv2mg = k v{_\ell}^2. Ce qui nous donne :
mv(t)+kv2(t)=kv2mv'(t) + kv^2(t) = k v{_\ell}^2
Comme m0m \neq 0 on peut donc écrire :
v(t)+kmv2(t)=kmv2v(t)=kmv2kmv2(t)v(t)=km(v2v2(t))v(t)=kmv2(1v2(t)v2)v(t)=kmv2(1(v(t)v)2)v'(t) + \dfrac{k}{m} v^2(t) = \dfrac{k}{m} v{_\ell}^2 \,\,\,\, \Longleftrightarrow \,\,\,\, v'(t) = \dfrac{k}{m} v{_\ell}^2 - \dfrac{k}{m} v^2(t) \,\,\,\, \Longleftrightarrow \,\,\,\, v'(t) = \dfrac{k}{m} \left( v{_\ell}^2 - v^2(t) \right) \,\,\,\, \Longleftrightarrow \,\,\,\, v'(t) = \dfrac{k}{m}v{_\ell}^2 \left( 1 - \dfrac{v^2(t)}{v{_\ell}^2} \right) \,\,\,\, \Longleftrightarrow \,\,\,\, v'(t) = \dfrac{k}{m}v{_\ell}^2 \left( 1 - \left(\dfrac{v(t)}{v{_\ell}}\right)^2 \right)
On pose G=kmv2G = \dfrac{k}{m}v{_\ell}^2. Dans ce cas, on obtient bien :
v(t)=G(1(v(t)v)2)\red{\boxed{v'(t) = G \left( 1 - \left( \dfrac{v(t)}{v_\ell} \right)^2 \right)}}
Question 3

Montrer que l'on a : V(t)=K(1V2(t))V'(t) = K \left( 1 - V^2(t) \right). Dans cette équation V(t)=v(t)vV(t) = \dfrac{v(t)}{v_\ell} et KK est une autre constante réelle dont vous préciserez l'expression.

Correction
On a :
v(t)=G(1(v(t)v)2)v(t)v=Gv(1(v(t)v)2)v'(t) = G \left( 1 - \left( \dfrac{v(t)}{v_\ell} \right)^2 \right) \,\,\,\, \Longleftrightarrow \,\,\,\, \dfrac{v'(t)}{v_\ell} = \dfrac{G}{v_\ell} \left( 1 - \left( \dfrac{v(t)}{v_\ell} \right)^2 \right)
Comme vv\ell est une constante réelle, on a alors :
v(t)v=(v(t)v)\dfrac{v'(t)}{v_\ell} = \left(\dfrac{v(t)}{v_\ell}\right)'
On en déduit donc que :
(v(t)v)=Gv(1(v(t)v)2)\left(\dfrac{v(t)}{v_\ell}\right)' = \dfrac{G}{v_\ell} \left( 1 - \left( \dfrac{v(t)}{v_\ell} \right)^2 \right)
Si on pose V(t)=v(t)vV(t) = \dfrac{v(t)}{v_\ell} on obtient :
V(t)=Gv(1V2(t))V'(t) = \dfrac{G}{v_\ell} \left( 1 - V^2(t) \right)
Puis, on posant K=Gv=kmvK = \dfrac{G}{v_\ell} = \dfrac{k}{m}v{_\ell}, on trouve que :
V(t)=K(1V2(t))\red{\boxed{V'(t) = K\left( 1 - V^2(t) \right)}}
Question 4

Déterminer la valeur des deux nombres réels AA et BB qui vérifient : 11V2(t)=A1+V(t)+B1V(t)\dfrac{1}{1 - V^2(t)} = \dfrac{A}{1+V(t)} + \dfrac{B}{1-V(t)}.

Correction
On a :
11V2(t)=A1+V(t)+B1V(t)11V2(t)=A(1V(t))+B(1+V(t))(1+V(t))(1V(t))11V2(t)=AAV(t)+B+BV(t)1V2(t)\dfrac{1}{1 - V^2(t)} = \dfrac{A}{1+V(t)} + \dfrac{B}{1-V(t)} \,\,\,\, \Longleftrightarrow \,\,\,\, \dfrac{1}{1 - V^2(t)} = \dfrac{A \left(1-V(t)\right) + B\left(1+V(t) \right)}{\left(1+V(t) \right) \left(1-V(t)\right)} \,\,\,\, \Longleftrightarrow \,\,\,\, \dfrac{1}{1 - V^2(t)} = \dfrac{A-AV(t)+B+BV(t)}{1 - V^2(t)}
Ce qui nous donne par identification des numlérateurs :
11V2(t)=A+B+(BA)V(t)1V2(t)1+0V(t)1V2(t)=A+B+(BA)V(t)1V2(t)\dfrac{1}{1 - V^2(t)} = \dfrac{A+B+(B-A)V(t)}{1 - V^2(t)} \,\,\,\, \Longleftrightarrow \,\,\,\, \dfrac{1 + 0 V(t)}{1 - V^2(t)} = \dfrac{A+B+(B-A)V(t)}{1 - V^2(t)}
Ce qui implique que (BA)=0(B-A) = 0 et aussi que A+B=1A+B=1. Ce qui nous donne B=AB = A et donc 2A=12A = 1. Ainsi A=B=12A = B = \dfrac{1}{2}. D'où :
11V2(t)=121+V(t)+121V(t)\dfrac{1}{1 - V^2(t)} = \dfrac{\dfrac{1}{2}}{1+V(t)} + \dfrac{\dfrac{1}{2}}{1-V(t)}
Finalement :
11V2(t)=12(11+V(t)+11V(t))\red{\boxed{\dfrac{1}{1 - V^2(t)} = \dfrac{1}{2} \left( \dfrac{1}{1+V(t)} + \dfrac{1}{1-V(t)} \right)}}
Question 5

Démontrer que l'on a : 1V(t)1+V(t)=Qe2Kt\dfrac{1 - V(t)}{1 + V(t)} = Q \, e^{-2Kt} ou QQ est une nouvelle constante réelle donc vous préciserez l'expression.

Correction
On sait que :
V(t)=K(1V2(t))V'(t) = K\left( 1 - V^2(t) \right)
Soit encore :
V(t)1V2(t)=K\dfrac{V'(t)}{1 - V^2(t)} = K
Soit encore :
V(t)12(11+V(t)+11V(t))=K(V(t)1+V(t)+V(t)1V(t))=2KV(t)1+V(t)V(t)1V(t)=2KV'(t) \dfrac{1}{2} \left( \dfrac{1}{1+V(t)} + \dfrac{1}{1-V(t)} \right) = K \,\,\,\, \Longleftrightarrow \,\,\,\, \left( \dfrac{V'(t)}{1+V(t)} + \dfrac{V'(t)}{1-V(t)} \right) = 2K \,\,\,\, \Longleftrightarrow \,\,\,\, \dfrac{V'(t)}{1+V(t)} - \dfrac{-V'(t)}{1-V(t)} = 2K
En primitivant, on obtient :
ln(1+V(t))ln(1V(t))=2Kt+q(qR)\ln \left( 1 + V(t) \right) - \ln \left( 1 - V(t) \right) = 2Kt + q \,\,\, (q \in \mathbb{R})
En utilisant les propriétés élémentaires associées aux logarithmes, on obtient :
ln(1+V(t)1V(t))=2Kt+q(qR)ln(1+V(t)1V(t))=2Ktqln(1V(t)1+V(t))=2Ktq\ln \left( \dfrac{ 1 + V(t)}{ 1 - V(t)} \right) = 2Kt + q \,\,\, (q \in \mathbb{R}) \,\,\,\, \Longleftrightarrow \,\,\,\, -\ln \left( \dfrac{ 1 + V(t)}{ 1 - V(t)} \right) = -2Kt - q \,\,\,\, \Longleftrightarrow \,\,\,\, \ln \left( \dfrac{ 1 - V(t)}{ 1 + V(t)} \right) = -2Kt - q
En prenant l'exponentielle de chacun des deux membres, on a :
eln(1V(t)1+V(t))=e2Ktq1V(t)1+V(t)=e2Kt×eqe^{\ln \left( \frac{ 1 - V(t)}{ 1 + V(t)} \right)} = e^{-2Kt - q} \,\,\,\, \Longleftrightarrow \,\,\,\, \frac{ 1 - V(t)}{ 1 + V(t)} = e^{-2Kt} \times e^{-q}
En posant Q=eqRQ = e^{-q} \in \mathbb{R}, on obtient finalement :
1V(t)1+V(t)=Qe2Kt\red{\boxed{\dfrac{1 - V(t)}{1 + V(t)} = Q \, e^{-2Kt}}}
Question 6

On note par v0=v(t=0)v_0 = v(t=0). Montrer que : v(t)=v(1vv0v+v0e2kmvt1+vv0v+v0e2kmvt)v(t) = v_\ell \, \left( \dfrac{1 - \dfrac{v_\ell - v_0}{v_\ell + v_0} e^{-2\frac{k}{m} v_\ell t}}{1 + \dfrac{v_\ell - v_0}{v_\ell + v_0} e^{-2\frac{k}{m} v_\ell t}} \right)

Correction
On a :
1V(t)1+V(t)=Qe2Kt\dfrac{1 - V(t)}{1 + V(t)} = Q \, e^{-2Kt}
Soit :
1V(t)=(1+V(t))Qe2Kt1V(t)=Qe2Kt+V(t)Qe2Kt1Qe2Kt=V(t)+V(t)Qe2Kt1Qe2Kt=V(t)(1+Qe2Kt)1 - V(t) = \left( 1 + V(t) \right) Q \, e^{-2Kt} \,\,\,\, \Longleftrightarrow \,\,\,\, 1 - V(t) = Q \, e^{-2Kt} + V(t) Q \, e^{-2Kt} \,\,\,\, \Longleftrightarrow \,\,\,\, 1 - Q \, e^{-2Kt} = V(t) + V(t) Q \, e^{-2Kt} \,\,\,\, \Longleftrightarrow \,\,\,\, 1 - Q \, e^{-2Kt} = V(t) \left( 1 + Q \, e^{-2Kt} \right)
Ce qui nous donne :
V(t)=1Qe2Kt1+Qe2Ktv(t)v=1Qe2Kt1+Qe2KtV(t) = \dfrac{1 - Q \, e^{-2Kt}}{ 1 + Q \, e^{-2Kt}} \,\,\,\, \Longleftrightarrow \,\,\,\, \dfrac{v(t)}{v_\ell} = \dfrac{1 - Q \, e^{-2Kt}}{ 1 + Q \, e^{-2Kt}}
Utilisons maintenant la condition initiale v0=v(t=0)v_0 = v(t=0). On a alors :
v(t=0)v=1Qe2K×01+Qe2K×0v0v=1Qe01+Qe0v0=1Q1+Qv0(1+Q)=1Qv0v+v0vQ=1QQ+v0vQ=1v0v\dfrac{v(t=0)}{v_\ell} = \dfrac{1 - Q \, e^{-2K\times 0}}{ 1 + Q \, e^{-2K\times 0}} \,\,\,\, \Longleftrightarrow \,\,\,\, \dfrac{v_0}{v_\ell} = \dfrac{1 - Q \, e^{0}}{ 1 + Q \, e^{0}} \,\,\,\, \Longleftrightarrow \,\,\,\, v_0 = \dfrac{1 - Q}{ 1 + Q} \,\,\,\, \Longleftrightarrow \,\,\,\, v_0 (1+Q) = 1 - Q \,\,\,\, \Longleftrightarrow \,\,\,\, \dfrac{v_0}{v_\ell} + \dfrac{v_0}{v_\ell} Q = 1 - Q \,\,\,\, \Longleftrightarrow \,\,\,\, Q + \dfrac{v_0}{v_\ell} Q = 1 - \dfrac{v_0}{v_\ell}
Soit encore :
Q(1+v0v)=1v0vQ(v+v0v)=vv0vQ \left( 1 + \dfrac{v_0}{v_\ell} \right) = 1 - \dfrac{v_0}{v_\ell} \,\,\,\, \Longleftrightarrow \,\,\,\, Q \left( \dfrac{v_\ell + v_0}{v_\ell} \right) = \dfrac{v_\ell - v_0}{v_\ell}
Comme v0v_\ell \neq 0, on peut donc simplifier et obtenir :
Q(v+v0)=vv0Q (v_\ell + v_0) = v_\ell - v_0
Donc Q=vv0v+v0Q = \dfrac{v_\ell - v_0}{v_\ell + v_0}. On a alors :
v(t)v=1vv0v+v0e2Kt1+vv0v+v0e2Ktv(t)=v(1vv0v+v0e2Kt1+vv0v+v0e2Kt)\dfrac{v(t)}{v_\ell} = \dfrac{1 - \dfrac{v_\ell - v_0}{v_\ell + v_0} \, e^{-2Kt}}{ 1 + \dfrac{v_\ell - v_0}{v_\ell + v_0} \, e^{-2Kt}} \,\,\,\, \Longleftrightarrow \,\,\,\, v(t) = v_\ell \left( \dfrac{1 - \dfrac{v_\ell - v_0}{v_\ell + v_0} \, e^{-2Kt}}{ 1 + \dfrac{v_\ell - v_0}{v_\ell + v_0} \, e^{-2Kt}} \right)
Or, on sait que K=kmvK = \dfrac{k}{m}v{_\ell}. Ce qui nous donne donc :
v(t)=v(1vv0v+v0e2kmvt1+vv0v+v0e2kmvt)\red{\boxed{v(t) = v_\ell \, \left( \dfrac{1 - \dfrac{v_\ell - v_0}{v_\ell + v_0} e^{-2\frac{k}{m} v_\ell t}}{1 + \dfrac{v_\ell - v_0}{v_\ell + v_0} e^{-2\frac{k}{m} v_\ell t}} \right)}}