lundi 30 novembre 2020

Trompe à eau (version 2)

C'est la lecture du livre de Dan Simmons : ``L'abominable'' qui a été le prétexte à la réalisation de cette version de la trompe à eau. La première version est ici :

http://pstricks.blogspot.com/2016/04/ebullition-et-pression-trompe-eauwater.html

Ce roman se passe dans l’Everest, la plus haute montagne du monde. C’est une histoire très réaliste sur les conditions d’ascension à l’époque où elle se déroule (1925), les faits qui ont marqué les tentatives pour vaincre le sommet de l’Everest, les personnages qui interviennent comme Winston Churchill(alors chancelier de l’échiquier), on y rencontre même Charlie Chaplin venu projeter son film “La Ruée vers l’or” et le moins sympathique Rudolf Hess(les nazis ont un rôle important dans l'histoire). C’est un livre dont les personnages et leurs aventures restent longtemps dans la mémoire.

L’extrait suivant à servi de prétexte à la révision du schéma :
Jean-Claude Clairoux et Jake Perry(le narrateur) arrivent au camp IV :
« Le Diacre et Tejbir Norgay lèvent des yeux surpris quand nous passons la porte de la tente Whymper et secouons la neige de nos vêtements dans le petit vestibule, avant de les rejoindre. J’imagine qu’on doit offrir un drôle de spectacle avec nos capuches en duvet d’oie remontées, nos casques d’aviateur en cuir, nos cagoules qui nous cachent tout le visage, les lampes frontales encore allumées, nos lunettes couvertes de givre et les épaules pleines de neige de nos anoraks Shackleton. Les deux hommes ne s’attendaient manifestement pas à avoir de la compagnie. Ils sont penchés sur un réchaud Unna, à faire bouillir quelque chose dans une grande casserole – à la température ridiculement basse nécessaire à l’ébullition à 7 200 mètres. L’eau bout à 76 °C à cette altitude, et à 100°C au niveau de la mer. Dès qu’ils entrent en contact avec l’air froid, nos liquides « bouillis » à 76 °C refroidissent à la température du corps. »

Quelques précisions sur la trompe à eau :

D’après Henri Bouasse, page 544 de “Jets, Tubes et Canaux” (1923) Librairie Delagrave(Paris), le cône qui relie les sections R1 et R2 ne doit pas être ni trop aigu ni trop évasé. Comme valeur optima, on trouve un angle compris entre 15° et 20°. Le rapport R/r ne peut dépasser 4 ou 5.


Les fichiers .tex et .pdf sont ici :

http://manuel.luque.free.fr/Ebullition/ebullition.zip 

ou

 trompe-à-eau-ebullition

 

 



 

lundi 23 novembre 2020

Sections circulaires d'un ellipsoïde avec pst-solides3d

 La solution de la question proposée au concours pour les deux académies de Montpellier et d'Aix (année 1870) a été rédigée par M. Auguste Macé, élève de Mathématiques spéciales au lycée de Grenoble.

http://www.numdam.org/item?id=NAM_1871_2_10__17_1

L’énoncé comporte 3 parties, la première est formulée ainsi :
« I. Étant donné un ellipsoïde, on détermine les points de contact des plans tangents parallèles aux plans des sections circulaires ; on mène deux sphères tangentes chacune en deux de ces points symétriques par rapport au grand axe : trouver l’équation des surfaces de révolution du second degré circonscrites à ces deux sphères. Classification et discussion de ces surfaces. »
On a compris qu’il s’agit de déterminer, dans la première étape de cette partie, les ombilics de l’ellipsoïde d’équation : 

$ \frac{x^2}{a^2}+\frac{y^2}{b^2}+\frac{z^2}{c^2}=1\quad\text{avec } a>b>c $

 


La démonstration d'Auguste Macé commence par : ``on sait que les plans des sections cycliques sont parallèles au plan'' :

$ x=z\sqrt{\frac{\dfrac{1}{c^2}-\dfrac{1}{b^2}}{\dfrac{1}{b^2}-\dfrac{1}{a^2}}} $

Démontrons que l'équation du plan de section circulaire passant par l'origine s'écrit bien ainsi.


La section circulaire passant par l'origine est un cercle de rayon égal au demi-axe suivant $Oy$ : $b$. Pour simplifier plaçons-nous dans le plan $xOz$. La section de l'ellipsoïde dans ce plan est une ellipse de demi-grand axe $a$ et de demi-petit axe $c$.

Les 2 sections circulaires passant par $O$ sont des cercles de rayon $b$. Déterminons les coordonnées des points $A_i$.

Ellipse : $ \frac{x^2}{a^2}+\frac{z^2}{c^2}=1$

Cercle de rayon $b$ :

\[
\left\{
\begin{array}[m]{l}
x  = b\cos\theta\\[1em]
z = b\sin\theta
\end{array}
\right.
\]

On remplace $x$ et $z$ dans l'équation de l'ellipse :
\[
\frac{b^2\cos^2\theta}{a^2}+\frac{b^2\sin^2\theta}{c^2}=1 \Rightarrow \frac{\cos^2\theta}{a^2}+\frac{\sin^2\theta}{c^2}=\frac{1}{b^2}
\]
On cherche $\theta$.
\[
\frac{\cos^2\theta}{a^2}+\frac{1-\cos^2\theta}{c^2}=\frac{1}{b^2}\Rightarrow \cos^2\theta\left[\frac{1}{a^2}-\frac{1}{c^2}\right]=\frac{1}{b^2}-\frac{1}{c^2}
\]
On en déduit :
\[
\cos\theta=\pm\sqrt{\dfrac{\dfrac{1}{c^2}-\dfrac{1}{b^2}}{\dfrac{1}{c^2}-\dfrac{1}{a^2}}}\quad\text{et}\quad \sin\theta=\pm\sqrt{\dfrac{\dfrac{1}{b^2}-\dfrac{1}{a^2}}{\dfrac{1}{c^2}-\dfrac{1}{a^2}}}
\]
\[
\tan\theta=\pm\sqrt{\dfrac{\dfrac{1}{b^2}-\dfrac{1}{a^2}}{\dfrac{1}{c^2}-\dfrac{1}{b^2}}}
\]
Les équations des plans de section circulaire passant par $O$ sont :
\[
z=\pm\sqrt{\dfrac{\dfrac{1}{b^2}-\dfrac{1}{a^2}}{\dfrac{1}{c^2}-\dfrac{1}{b^2}}}x
\]
Pour les coordonnées des 4 ombilics, on se reportera à la démonstration imparable d'Auguste Macé.
 Pour $O_1$ :
\[
\left\{
\begin{array}[m]{l}
x_1 = a\sqrt{\dfrac{a^2-b^2}{a^2-c^2}}\\[1em]
y_1 =0\\
z_1 = c\sqrt{\dfrac{b^2-c^2}{a^2-c^2}}
\end{array}
\right.
\]

Voici, par exemple, le plan tangent et la normale à l'ellipsoïde en $O_1$.


Le vecteur normal unitaire aux plans de section vaut :
\[
\left\{
\begin{array}[m]{l}
n_x  =\pm \dfrac{c}{b}\sqrt{\dfrac{a^2-b^2}{a^2-c^2}}\\[1em]
n_y=0\\
n_z  = \dfrac{a}{b}\sqrt{\dfrac{b^2-c^2}{a^2-c^2}}
\end{array}
\right.
\]
Si l'on considère un point d'un des plans situé sur le rayon qui joint $O$ à un ombilic (par exemple $O_1(x_1,0,z_1)$), ses coordonnées valent :
\[
\left\{
\begin{array}[m]{l}
x_P  =kx_1\\[1em]
y_P=0\\
z_P  = \dfrac{x_Pz_1}{x_1}
\end{array}
\right.
\]
La distance de $O$ au plan de section $d=|x_Pn_x+z_pn_z|$ avec $0<d<\dfrac{ac}{b}$, $\dfrac{ac}{b}$ est la distance de $O$ aux ombilics.

Le rayon des sections circulaires vaut :
\[
r=b\sqrt{1-\left(\frac{bd}{ac}\right)^2}
\]

On peut ajuster un cylindre de même rayon sur une section circulaire de l'ellipsoïde.

Découpe de 2 sections circulaires.


Découpe de 4 sections circulaires.

On découpe une tranche de section circulaire (une rondelle) :



On enlève une rondelle et on conserve les 2 parties restantes.

Toutes les images ont été obtenues avec pst-solides3d avec quelques options remarquables mises au point par Jean-Paul Vignault (en particulier la découpe d'un objet par un plan qui a été utilisée ici).

Pour faciliter les constructions, j'ai créé un \psSolid[objet=ellipsoid], dont les les exemples d'utilisation sont dans un fichier séparé (psEllipsoid.tex, psEllipsoid.pdf). Voici quelques images de cet objet :


 Une variante intéressante, permet, avec l’option [affinagerm], de conserver la face centrale et de la rendre transparente pour que l’intérieur soit visible, ou donner une couleur différente des armatures aux faces centrales.

 



Section par un plan horizontal, on enlève la calotte supérieure :

Les fichiers sont accessibles ici :

http://manuel.luque.free.fr/ellipsoid/Ellipsoid-sections.zip

ou sur drive :

 ellipsoid-sections-circulaires

Les fichiers pdf sont les fichiers contenant les documentations. Les fichiers .tex sont éventuellement à compiler.avec la séquence habituelle : LaTeX ->DVIPS->pst2pdf




samedi 21 novembre 2020

Courbes de poursuite avec PSTricks

C’est une commande destinée à illustrer le problème de la course poursuite, elle a été écrite en collaboration avec Jürgen Gilg. Le site suivant contient une superbe étude de ce problème :

https://mathcurve.com/courbes2d/poursuite/poursuite.shtml 

qui a inspiré ce package ainsi que l’étude théorique de Jürgen Gilg.
Le couple choisi dans l’étude que propose Jürgen Gilg dans un autre document est (chien, maître), mais on peut citer d’autres couples de protagonistes comme (renard, lapin) et plus généralement (prédateur, cible).
Le chien court après son maître en le gardant en point de mire. Les deux ont des vitesses constantes, différentes ou égales et on s’intéresse à la trajectoire suivie par le chien selon le chemin suivi par le maître. Les deux partent de leurs positions respectives en même temps. 

On s’intéresse à 2 cas :

a) le maître suit un parcours rectiligne ;

b) le maître suit un parcours circulaire.

Les équations différentielles ont été résolues avec le package https://ctan.org/pkg/pst-ode

 et les animations réalisées avec https://ctan.org/pkg/animate

d'Alexander Grahn.

Tous les fichiers sont ici :

http://manuel.luque.free.fr/pst-pursuitcurve/pst-pursuitcurve-v2.zip 

http://manuel.luque.free.fr/pst-pursuitcurve/pst-pursuitcurve.zip 

ou sur drive :

pst-pursuitcurve

 Quelques animations que vous retrouverez dans les documentations :





 

 

 

 


vendredi 13 novembre 2020

Le Botafumeiro - version 2

 Voici une version de l'animation du pendule du Botafumeiro qui avait été réalisée par Jürgen Gilg :

http://pstricks.blogspot.com/2019/02/lencensoir-de-botamufeiro-juergen-gilg.html

et qu'Alexander Grahn, l'auteur des packages `animate' et `pst-ode' utilisés par Jürgen, a amélioré en utilisant des options particulières du package pst-ode.

Ainsi la durée de l'animation est découpée en intervalles de temps égaux dt et après les premiers calculs sur le premier intervalle 0 -> dt, les calculs ne sont plus effectués que sur chacun des intervalles suivants en prenant comme conditions initiales celles qui correspondent au dernier point calculé. La commande :
\pstODErestoreState{curState}
permettant de continuer l'intégration à partir du dernier point calculé.
Le calcul numérique du système d'équations différentielles sur l'intervalle suivant se poursuit avec pstOdesolve en laissant vide les conditions initiales parce qu'on continue à partir de `curState'.
Les différents points calculés sont placés à la suite de ceux calculés dans la durée précédente avec l'option append.
Dernière précision, on obtient les coordonnées du dernier point calculé(xP,yP)(à chaque étape) avec, si vous avez nommé XY la variable pour stocker les points pour pstOdesolve :
 \pstVerb{
    mark
  XY
  /yP exch def
  /xP exch def
 cleartomark
  }
Mes explications sont peut-être un peu(beaucoup) confuses et il vous sera certainement plus utile de regarder le code dans le fichier source botafumeiro-AG.tex.

D'autres exemples utilisant cette procédure seront prochainement mis en ligne sur ce blog.

Les fichiers sont ici (botafumeiro-AG.tex) : 

Botafumeiro 

et

 http://manuel.luque.free.fr/Botafumeiro/Botafumeiro.zip

Une animation Gif obtenue à partir du pdf ( \usepackage[export]{animate})  :

 



 

 

mardi 10 novembre 2020

La courbe du nageur avec PSTricks

À l'initiative de Jürgen Gilg, voici ``la courbe du nageur'' revue avec PSTricks. Cette courbe est présentée et étudiée par Rober Ferréol sur :

https://mathcurve.com/courbes2d/nageur/nageur.shtml

Cette page propose le lien :  https://www.feynmanlectures.caltech.edu/info/exercises/boat_time.html

qui contient des solutions à l’exercice de l’ouvrage “Feynman’s Tips On Physics” dont voici l’énoncé : 

« 6 - 2 

A motorboat that runs at a
constant speed V relative to the water
is operated in a straight river channel
where the water is flowing smoothly
with a constant speed R. The boat is
first sent on a round trip from its
anchor point to a point a distance d
directly upstream. It is then sent on a
round trip from its anchor point to a
point a distance d away directly across
the stream. For simplicity assume that
the boat runs the entire distance in
each case at full speed and that no
time is lost in reversing course at the
end of the outward lap. If $t_V$ is the
time the boat took to make the round
trip in line with the stream flow, $t_A$ the
time the boat took to make the round
trip across the stream, and $t_L$ the time
the boat would take to go a distance
2d on a lake.
a) What is the ratio $t_v/t_A$?
b) What is the ratio $t_A/t_L$?

»

Il propose 2 méthodes pour atteindre la rive opposée(et revenir) qu’il nomme “the crabbing method” et “the pointing method”. C'est la deuxième méthode qui est développée par Jürgen Gilg. Il reprend et détaille les calculs de Robert Ferréol dans la deuxième partie du document.

Une commande \psCrossingRiver[options](a,b) a été créée pour faciliter la représentation des trajectoires et des vecteurs vitesses.

Les fichiers sont téléchargeables ici :

http://manuel.luque.free.fr/pst-crossingriver/pst-crossingriver.zip

ou, sur Drive :

pst-crossingriver.zip 

La résolution numérique du système d'équations différentielles a été réalisée avec le package  :

https://ctan.org/pkg/pst-ode

d'Alexander Grahn.

Voici 3 animations dont les sources sont dans l'archive (F5 pour relancer l'animation) :