CONCOURS ARTS ET MÉTIERS ParisTech - ESTP - POLYTECH
Épreuve de Physique-Modélisation PC
Durée 4 h
Si, au cours de l'épreuve, un candidat repère ce qui lui semble être une erreur d'énoncé, d'une part il le signale au chef de salle, d'autre part il le signale sur sa copie et poursuit sa composition en indiquant les raisons des initiatives qu'il est amené à prendre.
L'usage de calculatrices est autorisé.
AVERTISSEMENT
Le candidat devra porter l'ensemble de ses réponses sur le cahier réponses, à l'exclusion de toute autre copie. Les résultats doivent être reportés dans les cadres prévus à cet effet.
Remarques préliminaires
Les explications des phénomènes étudiés interviennent dans la notation au même titre que les développements analytiques et les applications numériques.
Les résultats exprimés sans unité ne seront pas comptabilisés.
Plusieurs questions demandent une explication qualitative. Il est attendu des réponses claires et concises (moins de 10 lignes).
Tout résultat fourni dans l'énoncé peut être admis et utilisé par la suite, même s'il n'a pas été démontré par le(la) candidat(e).
Tout au long de l'énoncé, les paragraphes en italiques ont pour objet de fournir des informations et d'aider à la compréhension du problème mais ne contiennent pas de questions.
Remarques pour les questions de programmation
Toutes les questions d'informatique comportent une mention du type «Écrire en langage Python. . . ». Les codes doivent être écrits en langage Python.
On se limitera aux types suivant : entiers, flottants, chaines de caractères, listes, tableaux (array du module numpy) et tuples.
On se limitera aux mots clés suivants : if, elif, else, while, for, in, def, return, and, or, not, True, False, import, from, as et None.
On se limitera aux fonctions et méthodes de la bibliothèque standard suivante : print, plot, range, enumerate, len et append.
Les codes ne respectant pas les consignes précédentes ne seront pas comptabilisés.
La présentation, la lisibilité, l'orthographe, la qualité de la rédaction, la clarté et la précision des raisonnements entreront pour une part importante dans l'appréciation des copies. En particulier, les résultats non justifiés ne seront pas pris en compte. Les candidats sont invités à encadrer les résultats de leurs calculs.
Tournez la page S.V.P.
Il est interdit aux candidats de signer leur composition ou d'y mettre un signe quelconque pouvant indiquer sa provenance.
Rappel sur la notation complexe
Dans tout le sujet, on associe à une grandeur sinusoïdale la grandeur complexe avec (où ) et telle que , où Re est la partie réelle.
Données
distance Terre-Soleil :
rayon de la Terre :
rayon du Soleil :
masse du Soleil :
constante de gravitation universelle :
constante d'Avogadro :
masse d'un proton :
masse d'un électron :
charge élémentaire :
constante des gaz parfaits :
constante de Boltzmann :
vitesse de lumière dans le vide :
permittivité du vide :
perméabilité du vide :
Ce problème constitué de deux parties indépendantes s'intéresse à la production et à l'acheminement de l'énergié électrique produite par une centrale photovoltaïque. La première partie étudie l'origine de l'énergie photovoltaïque, le Soleil, ainsi que la puissance électrique fournie par une cellule photovoltaïque, La deuxième partie étudie l'acheminement de l'énergie électrique jusqu'à l'utilisateur, Les sous-parties sont pour la plupart indépendantes les unes des autres.
PREMIÈRE PARTIE Étude du Soleil et de son rayonnement
Le corps noir est un objet idéal qui absorbe toute l'énergie électromagnétique qu'il reçoit, sans en réfléchir ni en transmettre.
En l'absence d'énergie électromagnétique extérieure, un corps noir à la température d'équilibre émet un flux surfacique d'énergie électromagnétique dont la densité spectrale, c'est-è-dire le flux surfacique par unité de longueur d'onde émise, dépend uniquement de la longueur d'onde et de la température.
Le maximum de cette densité spectrale de flux surfacique est donné par la loi de Wien
avec en mètres et en kelvins. Cette dernière loi exprime le fait que pour un corps noir, le produit de la température et de la longueur d'onde du pic de la courbe est toujours égal à une constante. Cette loi très simple permet ainsi de connaître la température d'un corps assimilé à un corps noir par la seule position de son maximum.
D'après la loi de Stefan-Boltzmann, le flux surfacique d'énergie (en ) émis par le corps noir varie en fonction de la température absolue (exprimée en kelvin) selon la formule
où est la constante de Stefan-Boltzmann qui vaut environ .
On s'intéresse aux caractéristiques du rayonnement solaire (documents 1 et 2). Le candidat peut utiliser les informations des documents afin de répondre aux questions.
Dans cette partie, on assimile le Soleil à un corps noir de rayon .
A1. Estimer la valeur numérique de la température du Soleil assimilé à un corps noir.
Le raisonnement devra être explicité.
A2. Exprimer le flux surfacique d'énergie émis par le Soleil en fonction de et de . En déduire l'expression de la puissance totale rayonnée par le Soleil en fonction de et .
La distance Terre-Soleil étant très grande, les rayons solaires peuvent être considérés comme arrivant parallèlement entre eux. Ainsi, la Terre reçoit la même puissance que celle que recevrait un disque de rayon placé perpendiculairement aux rayons solaires incidents.
A3. Exprimer le flux surfacique d'énergie reçu par la Terre en fonction de et et montrer que la puissance totale reçue par la Terre, notée , s'écrit
Faire l'application numérique de .
En réalité, on mesure que la Terre reçoit au niveau du sol un flux surfacique total (intégré sur toutes les longueur d'ondes) d'énergie d'environ .
A4. Proposer une explication pour l'écart entre la valeur trouvée à la question précédente pour et celle mesurée de . On pourra s'aider des documents.
A5. A partir de , estimer l'énergie reçue en un jour par la Terre. Comparer cette valeur à la consommation journalière de l'humanité valant environ . Commenter la pertinence de développer l'énergie photovoltaïque pour assurer les besoins énergétiques de l'humanité.
B/ Estimation de la température du Soleil
On se propose de retrouver l'ordre de grandeur de la température du Soleil par un modèle thermodynamique.
Le Soleil est assimilé à une sphère de rayon , de centre et de masse . La masse volumique est supposée constante et égale à la masse volumique moyenne : c'est l'hypothèse notée . On utilisera les coordonnées sphériques (Figure 1).
B1. Montrer par des considérations d'invariances et de symétries que l'expression du champ de gravitation créé par le Soleil à une distance de son centre se met sous la forme .
B2. Exprimer la masse volumique moyenne en fonction de et . Par analogie avec l'électrostatique, utiliser le théorème de Gauss pour déterminer l'expression de lorsque en fonction de et .
On s'intéresse à un volume mésoscopique du Soleil, centré en un point situé à la distance du centre , dont la vitesse est notée dans le référentiel supposé galiléen héliocentrique et dont la pression est . On néglige la viscosité.
B3. Donner la définition d'un volume mésoscopique. Rappeler sans démonstration l'expression de la force volumique équivalente aux forces de pression. En appliquant la deuxième loi de Newton au volume mésoscopique considéré, retrouver l'équation d'Euler.
Pour simplifier l'étude, le fluide constituant le Soleil est supposé au repos, Ceci constitue l'hypothèse notée .
Figure 1 - Représentation des coordonnées sphériques. La base ( ) est orthonormée.
B4. Simplifier l'équation d'Euler précédente pour montrer que la loi de la statique des fluides à l'intérieur du Soleil s'écrit
On note Po la pression au centre du Soleil. L'expression du gradient en coordonnées sphériques est rappelée
B5. À partir de l'équation (3), déterminer l'expression de en fonction de , et . En supposant que la pression à l'extérieur du Soleil est nulle , déterminer l'expression de .
Pour déterminer la température à l'intérieur du Soleil, celui-ci est considéré comme constitué d'un gaz parfait (hypothèse ) d'hydrogène totalement ionisé, mélange équimolaire de protons et d'électrons .
B6. Justifier que la masse molaire moyenne du gaz vaut environ et effectuer l'application numérique. Quelle est de plus la relation dans ce modèle entre et la pression , la température , la constante des gaz parfaits et la masse volumique du Soleil?
La photosphère est une des couches externes de l'étoile, dont le rayon externe
correspond à la définition du rayon de l'étoile, et qui produit entre autres la lumière visible.
La lumière qui y est produite contient toutes les informations sur la température du rayonnement émis, la gravité de surface et la composition chimique de l'étoile. Pour le Soleil, la photosphère a une profondeur d'environ 400 kilomètres.
B7. Déduire des questions précédentes l'expression de la température dans le Soleil en fonction de et . Quelle est la partie la plus chaude de la photosphère? En considérant que la lumière est produite sur la couche interne de la photosphère, calculer la température à laquelle est émis le rayonnement électromagnétique du Soleil. Comparer à la valeur obtenue dans la partie A.
B8. D'après votre culture scientifique, discuter de la validité des hypothèses et . Peut-on considérer un gaz de particules chargées comme un gaz parfait? Discuter de la validité de l'hypothèse . Que dire du modèle proposé?
C/ Étude d'une centrale photovoltaïque
On s'intéresse à une centrale de taille comparable à celle de Martillac, près de Bordeaux, permettant d'alimenter en électricité une trentaine d'habitations. Elle est constituée 126 modules de 4 panneaux. Chaque panneau est formé de cellules photovoltaïques, et chaque cellule possède une taille . Cette centrale produit environ 100 kW dans de bonnes conditions d'éclairement des modules. La caractéristique d'une cellule est donnée sur la figure suivante (Figure 2).
Figure 2 - Caractéristique d'une cellule photovoltaïque en convention générateur pour plusieurs valeurs différentes de l'éclairement.
C1. Décrire précisément mais succinctement un protocole expérimental permettant de mesurer la caractéristique d'une cellule photovoltaïque telle que présentée sur la Figure 2.
On considère une cellule photovoltaïque recevant un flux . Les points expérimentaux correspondants sont enregistrés dans un environnement Python dans deux
objets de type liste : Liste_U contenant les valeurs des tensions mesurées et Liste_I les valeurs des intensités correspondantes.
C2. Écrire une fonction MaxPuissance en langage Python prenant en argument Liste_U et Liste_I et renvoyant la tension et l'intensité correspondant au maximum de la puissance fournie par la cellule photovoltaïque.
Les valeurs trouvées sont et . Les cellules photovoltaïques d'un panneau sont en série tandis que les 4 panneaux d'un module sont en parallèle.
C3. Déterminer la tension aux bornes d'un module, l'intensité traversant ce module et la puissance délivrée par ce module. Retrouver l'ordre de grandeur de 100 kW délivré par la centrale solaire de Martillac.
C4. D'après vous, en vous appuyant sur le document 4 et en argumentant, quelle est la composition des cellules utilisées dans la centrale de Martillac?
C5. Quel est l'ordre de grandeur de la puissance délivrée par la centrale lorsque, par temps partiellement nuageux, l'éclairement baisse à ? Commenter sur l'usage des centrales photovoltaïques.
Document 4 - Les différentes technologies de cellules photovoltaïques
D'après EDF
Différentes technologies entrent aujourd'hui dans la composition des installations photovoltaïques :
Le silicium cristallin
Le silicium amorphe
Le cuivre/indium/sélénium
Le cuivre/indium/gallium/sélénium
Les panneaux solaires à base de silicium cristallin sont les plus anciens. Ils se décomposent eux-mêmes en deux variantes : le monocristallin et le polycristallin. Ces deux variantes sont aujourd'hui très proches aussi bien en termes de rendement qu'en termes de coût. Le rendement d'un panneau photovoltaïque correspond à la quantité d'énergie solaire transformée par le panneau en électricité consommable, par rapport à l'énergie captée. Le rendement moyen d'un panneau cristallin du marché est de .
La souplesse mécanique du silicium amorphe lui permet d'être essentiellement utilisé dans des complexes de type «membrane solaire» ou « tôle solaire». Le rendement moyen des panneaux solaires à base de silicium amorphe est de 6 à .
Aujourd'hui, des technologies émergent à base de cuivre / indium / sélénium et de cuivre / indium / gallium / sélénium. Elles offrent de grandes perspectives en termes de coût et de rendement.
DEUXIÈME PARTIE
Transport de l'énergie électrique de la centrale au consommateur
L'énergie produite par une centrale photovoltaïque est dans certains cas directement injectée dans le réseau électrique basse tension pour pouvoir être consommée localement.
Dans toute cette partie, on s'intéresse au transport électrique monophasé de cette énergie, c'est-à-dire utilisant deux câbles électriques, principalement employé pour alimenter les zones peu denses en habitations.
Pour des raisons historiques et techniques, le réseau électrique basse tension fonctionne avec des tensions alternatives. En France, la fréquence utilisée vaut tandis que la tension efficace sur le réseau considéré est .
La transformation de la tension continue délivrée par la centrale photovoltaïque en signal électrique transportable sur le réseau est étudiée dans la partie . Dans la partie on examine le dimensionnement des câbles utilisés pour le transport de l'énergie sur le réseau. Enfin, la partie F s'intéresse au choix des caractéristiques électriques de la ligne monophasée pour l'adaptation au transport de l'énergie jusqu'aux installations domestiques. Ces trois parties sont indépendantes.
D/ Transformation en courant alternatif grâce à un onduleur
Pour pouvoir consommer l'énergie produite par la centrale, il faut transformer la tension continue supposée constante et délivrée par un module de 4 panneaux (voir la partie ) en tension alternative de fréquence et de tension efficace . Le principe de cette transformation se décompose en trois étapes, illustrées sur la Figure 3 :
tout d'abord, l'onduleur de tension autonome positionné après les panneaux photovoltaïques de la centrale transforme le signal continu en signal alternatif;
ensuite, une opération de filtrage est nécessaire pour rendre la tension de sortie de l'onduleur la plus proche possible d'un signal sinusoïdal à 50 Hz ;
la dernière étape, qui ne sera pas étudiée ici, consiste à amplifier cette tension pour que sa tension efficace soit de 230 V .
Pour réaliser la première étape, on étudie un onduleur de tension autonome à commande symétrique dans un premier temps puis à commande décalée dans un second temps.
Pour un onduleur autonome à commande symétrique, les interrupteurs représentés sur la Figure 3 s'ouvrent et se ferment en fonction du temps, noté , selon la séquence suivante, avec et un entier relatif :
éé
Figure 3 - Représentation schématique du circuit électrique, de la centrale au réseau électrique.
D1. Représenter le schéma électrique équivalent de l'onduleur lorsque et lorsque . Représenter alors l'allure de la tension en sortie de l'onduleur. Quelle est la tension efficace de ?
Cette tension n'étant pas sinusoïdale, un filtrage est nécessaire. On modélise de manière très simple le transformateur par une résistance et une inductance négligeable devant . Le circuit équivalent est représenté sur la Figure 4.
Figure 4 - Schéma du filtre équivalent.
Pour étudier l'influence de ce filtre sur chacun des harmoniques de la tension , on étudie le comportement de ce filtre lorsque la tension d'entrée est un signal sinusoïdal de fréquence .
D2. Rappeler l'expression de l'impédance d'une inductance ainsi que celle, , d'une résistance . En déduire que l'expression de la fonction de transfert est
où l'expression de est à déterminer en fonction de et .
D3. Quelle est la nature de ce filtre ? Est-il adapté pour filtrer en un signal sinusoïdal à 50 Hz ? Déterminer le gain du filtre. Exprimer la fréquence de coupure à -3 dB de ce filtre, notée , en fonction de .
On choisit pour la suite la valeur de telle que . La tension peut se décomposer en une série de Fourier selon :
avec les amplitudes des harmoniques de rang valant .
D4. En sortie du filtre, quels sont les rangs des harmoniques présents dans le signal ? L'amplitude de l'harmonique de rang trois est-elle négligeable devant celle du fondamental? Même question avec l'harmonique de rang cinq d'amplitude . Commenter.
Tout en conservant exactement le même filtre avec les mêmes composants, il est possible d'améliorer le filtrage en jouant sur la séquence d'ouverture et de fermeture des interrupteurs. Dans ce cas, l'onduleur est dit à «commande décalée».
L'ouverture et la fermeture des interrupteurs est commandée à partir de la comparaison du signal du réseau pré-existant et d'un signal de commande de forme triangulaire, d'amplitude avec et de fréquence représenté sur la Figure 5. Le fonctionnement est donné par la séquence suivante :
éééé
Il y a donc toujours deux interrupteurs ouverts et deux interrupteurs fermés.
On se propose d'étudier le principe de cet onduleur numériquement en langage python. On commence par définir les deux fonctions et servant au pilotage des interrupteurs grâce au programme suivant.
from math import sqrt, sin #racine carrée et sinus
fr = 50. #fréquence du réseau en Hertz
T = 1./fr
A = sqrt(2)*230 #amplitude en Volts
def Ur(t):
return A*sin(2*pi*fr*t)
fp = 400. #fréquence du signal de commande
Tp = 1./fp
alpha = 1.10
Ap = A*alpha
def p(t):
if t<0:
return p(-t)
n = t%Tp
if 0<=n<=Tp/2:
return Ap*(4/Tp*n-1)
else:
return Ap*(-4/Tp*(n-Tp/2)+1)
Figure 5 - Représentation des fonctions de commande et en fonction du temps pour l'onduleur à commande décalée. Ce graphique est donné à titre indicatif afin d'aider le candidat si besoin.
D5. Commenter et justifier la façon dont a été définie la fonction dans le code ci-dessus.
D6. Dans cette nouvelle séquence, on montre que si et sont fermés alors , si et sont fermés alors et sinon. On associe alors à chaque interrupteur un entier valant 0 si l'interrupteur est ouvert et 1 s'il est fermé. Élaborer une fonction Python nommée tension prenant en arguments les quatre valeurs des interrupteurs et et renvoyant la tension en sortie de l'onduleur.
D7. Élaborer une fonction Python nommée onduleur prenant comme argument un flottant représentant le temps , temps auquel sont évaluées les conditions d'ouverture de la séquence donnée précédemment à l'équation (8), et renvoyant la valeur de la tension en sortie de l'onduleur à cet instant .
Pour pouvoir tracer l'allure de la tension , on souhaite définir une liste de valeurs de temps régulièrement espacés compris entre 0 et tous les deux inclus. Pour la suite, on choisit et , où est la période du signal du réseau.
D8. Donner les instructions en langage Python pour construire les deux listes à enregistrer dans les variables respectives Liste_t et Liste_u et contenant respectivement les valeurs des temps et des tensions .
D9. Écrire alors en langage Python les commandes permettant de tracer le graphique de et de se représenter le fonctionnement de l'onduleur comme dans l'exemple page suivante de la Figure 6. On ne se souciera ni des légendes ni des axes.
D10. Quelle est la fréquence d'échantillonnage de ce signal? Justifier qualitativement pourquoi un si grand nombre de points de calcul du signal a été choisi?
Figure 6 - Représentation de sur une période et demi en sortie de l'onduleur à commande décalée étudié.
Grâce au module fftpack de la bibliothèque scipy, on réalise la transformée de Fourier de la tension . On dispose alors de deux listes de éléments contenant les fréquences et les amplitudes correspondantes. Ces deux listes, permettent de tracer le spectre de représenté sur la Figure 7.
Figure 7 - Spectre d'amplitude de la tension en fonction des fréquences.
D11. Écrire la ligne de code Python permettant d'importer le module fftpack de la bibliothèque scipy.
D12. À partir du spectre présenté sur la Figure 7, vérifier si la fréquence du fondamental, notée , est compatible avec celle déduite du signal temporel de la Figure 6, et commenter la présence des autres harmoniques.
On trouve et .
D13. Calculer les amplitudes relatives des harmoniques et après le filtre. Quelle est l'allure du signal de sortie en sortie de filtre lorsque l'onduleur est à commande décalée ? Quel avantage y-aurait-il à utiliser un onduleur à commande décalée dans ce montage ?
E/ Dimensionnement des câbles
On s'intéresse dans cette partie à la dimension des câbles employés pour transporter l'énergie électrique de la centrale au consommateur. On cherche à justifier le rayon des câbles utilisés en début du réseau basse tension.
Pour cela, on adopte le modèle de Drude : un électron libre de charge - est soumis à la force qu'exerce un champ électromagnétique et à une force de frottement visqueux, modélisant les collisions, de la forme avec la masse d'un électron, un temps de relaxation et la vitesse des électrons.
Un fil infini d'axe et de rayon est parcouru par un vecteur densité de courant . Le milieu, supposé électriquement neutre, contient électrons mobiles par unité de volume. Il est suffisamment dilué pour pouvoir négliger les interactions entre les différentes charges du milieu.
On utilise les coordonnées cylindrique de base ( ). On considère un champ électromagnétique monochromatique de la forme pour le champ électrique et pour le champ magnétique.
Le mouvement d'un électron du milieu conducteur est non relativiste et il est étudié dans le référentiel terrestre supposé galiléen.
E1. Dans le modèle présenté ci-dessus, appliquer la deuxième loi de Newton à un électron du milieu et donner l'équation différentielle vérifiée par sa vitesse . Faut-il prendre en compte le poids de l'électron? Justifier.
E2. On note la perméabilité du vide et la permittivité du vide. Écrire pour le milieu considéré les quatre équations de Maxwell.
E3. Justifier que la force magnétique subie par un électron est négligeable devant la force électrique. Simplifier alors l'équation du mouvement d'un électron.
On se place en régime permanent sinusoïdal et on note la vitesse de l'électron dans ce régime.
E4. Exprimer en fonction de et . En déduire l'expression du vecteur densité de courant en fonction des mêmes variables et de .
E5. Rappeler l'expression de la loi d'Ohm locale en fonction de et de la conductivité électrique . En déduire que la conductivité complexe en régime permanent sinusoïdal s'exprime
et donner l'expression de en fonction de et .
Un milieu conducteur tel que le câble étudié est caractérisé par un temps de relaxation de l'ordre de et une densité de porteurs de charge de l'ordre de . On rappelle que la fréquence du signal considéré est .
E6. Calculer . Par des calculs d'ordre de grandeurs, simplifier l'expression de la conductivité ainsi que l'équation de Maxwell-Ampère.
On rappelle que pour un champ vectoriel , on a la relation avec le laplacien vectoriel.
E7. Déduire des équations de Maxwell et de la loi d'Ohm locale que le vecteur densité de courant vérifie
Déterminer l'expression de et le calculer numériquement dans le cas étudié. Quel nom donne-t-on couramment à ?
La résolution de cette équation n'est pas demandée. La solution de l'équation précédente en coordonnées cylindriques conduit à un vecteur de la forme avec une certaine expression pour .
Figure 8 - À gauche: Représentation de l'épaisseur externe b. À droite : Représentation de pour les trois valeurs et .
On représente alors sur la Figure 8 l'intensité normalisée avec
l'intensité circulant dans l'épaisseur b la plus externe du conducteur et l'intensité totale.
E8. En étudiant les courbes représentatives de pour différents rayons, quelle est la zone du câble la plus sollicitée pour transporter le courant? Est-il utile de fabriquer des câbles dont le rayon vaut plusieurs fois ? Justifier.
E9. Au début du réseau basse tension, la section des câbles est . Commenter.
E10. En pratique, pour faire circuler des courants intenses, on utilise plusieurs fils de faible rayon isolés les uns des autres, plutôt qu'un seul fil de gros rayon. Justifier ce choix.
E11. Proposer une explication permettant à d'être supérieur à 1 pour certaines valeurs de , comme cela peut s'observer sur la Figure 8 dans le cas où .
On s'intéresse maintenant à une méthode numérique permettant de calculer les intégrales définissant et .
E12. Expliquer en quelques phrases et avec un schéma le principe d'intégration numérique par la méthode des rectangles.
E13. Élaborer une fonction Python nommée Rectangles prenant comme arguments une fonction , deux flottants et et un nombre de points de calcul et renvoyant la valeur approchée par la méthode des rectangles de calculée sur points.
E14. Quelle est la complexité de calcul de cet algorithme? On donnera la réponse sous la forme où est une grandeur à déterminer en fonction de .
On souhaite comparer la précision de la méthode des rectangles avec deux autres méthodes, celle des trapèzes et celle de Simpson. On note la valeur approchée de par une méthode numérique (rectangles, trapèzes ou Simpson) utilisant points de calcul. Cette comparaison est réalisée sur un exemple dont la solution exacte est simple : .
On calcule l'erreur . La Figure 9 donne l'erreur en fonction du nombre de points de calcul pour les trois méthodes considérées en échelle log-log.
E15. A partir de la Figure 9, estimer l'évolution de l'erreur en fonction de pour des valeurs de pas trop élevées. On donnera la réponse sous la forme où est une grandeur à déterminer en fonction de . Comparer les trois méthodes d'intégration numérique.
E16. Passer une certaine valeur de , on remarque que l'erreur augmente légèrement avec n. Ceci est visible sur la Figure 9 pour la méthode des trapèzes ou de Simpson. Proposer une explication à cela.
Figure 9 - Erreur des différentes méthodes en fonction du nombre de points de calcul en échelle
F/ Propagation de la tension le long d'une ligne électrique
On se propose ici d'étudier la propagation des signaux électriques dans une ligne monophasé, assimilée à deux câbles (en réalité il y en a plus, voir partie précédente) de longueur .
Dans un premier temps, un tronçon de longueur de cette ligne monophasée, que l'on considère très petite devant toutes les longueurs caractéristiques du phénomène, peut-être modélisé par une inductance et une capacité ( et étant constantes) selon le schéma de la Figure 10.
Figure 10 - Schématisation de la ligne électrique monophasée entre et à l'instant .
F1. L'ordre de grandeur de est de 50 km . L'Approximation des Régimes QuasiStationnaire (ARQS) est-elle valable en tout point du câble? Justifier. On rappelle que la fréquence du signal est .
F2. Montrer que l'application des lois de l'électrocinétique dans l'ARQS permet d'établir les deux équations suivantes
F3. En déduire que vérifie une équation de D'Alembert et exprimer la célérité des ondes associées en fonction de et .
On cherche des solutions sous la forme d'ondes planes progressives harmoniques (OPPH) de pulsation .
F4. Déterminer la relation de dispersion des ondes dans ce milieu entre la pulsation et le nombre d'onde noté . Donner la définition d'un milieu dispersif et d'un milieu absorbant. Le milieu est-il dispersif? Est-il absorbant? Les réponses devront être succinctement justifiées.
On cherche une solution sous la forme de la somme d'une OPPH incidente notée se propageant vers les positifs et d'une onde réfléchie se propageant vers les négatifs. On note respectivement et les amplitudes complexes de ces ondes. Chaque onde de tension est de plus associée à un courant et , d'amplitudes et .
On définit enfin l'impédance du milieu, qui dépend du type d'onde qui s'y propage.
F5. Donner la forme des OPPH incidentes et réfléchies en tension. En utilisant les relations de l'équation (11), exprimer l'impédance du milieu pour les ondes incidentes, notée , et pour les ondes réfléchies , en fonction de et . Quelle est l'unité des impédances et ?
La tension d'entrée du réseau basse tension, en , est où . Le circuit est fermé en sur une habitation modélisée par une simple résistance .
F6. Écrire, sans les résoudre, les deux équations correspondantes aux conditions de bord et permettant d'obtenir et .
F7. Déterminer l'expression de en fonction de et afin d'annuler l'onde réfléchie. Qualitativement, quel est l'avantage de cette situation d'un point de vue énergétique?
En réalité, il y a des pertes d'énergie dues à la nature résistive des câbles. On ajoute une résistance en série avec l'inductance et une autre résistance linéique de valeur , en parallèle avec la capacité (c'est-à-dire que gdx est l'admittance de cette résistance).
F8. Quelles sont les deux nouvelles équations différentielles couplées et du premier ordre en et en vérifiées par et ?
F9. On considère des solutions de la forme et . Montrer que est solution de l'équation différentielle
et .
Figure 11 - Schématisation de la ligne électrique monophasée avec pertes entre et à l'instant .
On note tel que et . Les solutions de l'équation précédente conduisent à
F10. Sans chercher à déterminer les expressions de et de , préciser si ce milieu est dispersif et/ou absorbant.
En choisissant judicieusement les composants et la géométrie de la ligne, il est possible de respecter la condition . Dans ce cas, on trouve et . On montre aussi que est minimal lorsque cette condition est respectée.
F11. Quels sont les avantages à choisir les paramètres des câbles tels que ?
On admet que sous la condition , la résistance permettant d'annuler l'onde réfléchie est la même que celle obtenue dans la question F7. Les ordres de grandeurs sont et .
F12. Calculer numériquement lorsque . En l'absence d'onde réfléchie, calculer la variation relative de l'amplitude de la tension aux bornes de lorsque l'on prend en compte les pertes dans la ligne par rapport au cas sans perte. Commenter.
s!u」noł słuəunnoop sәุ」de,d - ggLL LL - дSIOHJ NI
Épreuve/sous-épreuve :
L'usage de calculatrice est autorisé.
Cahier réponses
Épreuve de Physique-Modélisation
PC
Concours e3a-2017
Toutes les réponses seront portées sur ce cahier de réponses à l'exclusion de toute autre copie
NE PAS DÉGRAFER
A/ Approche descriptive du rayonnement du Soleil
A1. Estimer la valeur numérique de la température du Soleil assimilé à un corps noir. Le raisonnement devra être explicité.
A2. Exprimer le flux surfacique d'énergie émis par le Soleil en fonction de et de . En déduire l'expression de la puissance totale rayonnée par le Soleil en fonction de et .
A3. Exprimer le flux surfacique d'énergie reçu par la Terre en fonction de et et montrer que la puissance totale reçue par la Terre, notée , s'écrit
Faire l'application numérique de .
A4. Proposer une explication pour l'écart entre la valeur trouvée à la question précédente pour et celle mesurée de . On pourra s'aider des documents.
A5. À partir de , estimer l'énergie reçue en un jour par la Terre. Comparer cette valeur à la consommation journalière de l'humanité valant environ . Commenter la pertinence de développer l'énergie photovoltaïque pour assurer les besoins énergétiques de l'humanité.
B/ Estimation de la température du Soleil
B1. Montrer par des considérations d'invariances et de symétries que l'expression du champ de gravitation créé par le Soleil à une distance de son centre se met sous la forme .
B2. Exprimer la masse volumique moyenne en fonction de et . Par analogie avec l'électrostatique, utiliser le théorème de Gauss pour déterminer l'expression de lorsque en fonction de et .
B3. Donner la définition d'un volume mésoscopique. Rappeler sans démonstration l'expression de la force volumique équivalente aux forces de pression. En appliquant la deuxième loi de Newton au volume mésoscopique considéré, retrouver l'équation d'Euler.
B4. Simplifier l'équation d'Euler précédente pour montrer que la loi de la statique des fluides à l'intérieur du Soleil s'écrit
B5. À partir de l'équation (3), déterminer l'expression de en fonction de , et . En supposant que la pression à l'extérieur du Soleil est nulle , déterminer l'expression de .
B6. Justifier que la masse molaire moyenne du gaz vaut environ et effectuer l'application numérique. Quelle est de plus la relation dans ce modèle entre et la pression , la température , la constante des gaz parfaits et la masse volumique du Soleil?
B7. Déduire des questions précédentes l'expression de la température dans le Soleil en fonction de et . Quelle est la partie la plus chaude de la photosphère ? En considérant que la lumière est produite sur la couche interne de la photosphère, calculer la température à laquelle est émis le rayonnement électromagnétique du Soleil. Comparer à la valeur obtenue dans la partie A.
B8. D'après votre culture scientifique, discuter de la validité des hypothèses et . Peut-on considérer un gaz de particules chargées comme un gaz parfait? Discuter de la validité de l'hypothèse . Que dire du modèle proposé?
C/ Étude d'une centrale photovoltaïque
C1. Décrire précisément mais succinctement un protocole expérimental permettant de mesurer la caractéristique d'une cellule photovoltaïque telle que présentée sur la Figure 2.
C2. Écrire une fonction MaxPuissance en langage Python prenant en argument Liste_U et Liste_I et renvoyant la tension et l'intensité correspondant au maximum de la puissance fournie par la cellule photovoltaïque.
C3. Déterminer la tension aux bornes d'un module, l'intensité traversant ce module et la puissance délivrée par ce module. Retrouver l'ordre de grandeur de 100 kW délivré par la centrale solaire de Martillac.
C4. D'après vous, en vous appuyant sur le document 4 et en argumentant, quelle est la composition des cellules utilisées dans la centrale de Martillac ?
C5. Quel est l'ordre de grandeur de la puissance délivrée par la centrale lorsque, par temps partiellement nuageux, l'éclairement baisse à ? Commenter sur l'usage des centrales photovoltaïques.
D/ Transformation en courant alternatif grâce à un onduleur
D1. Représenter le schéma électrique équivalent de l'onduleur lorsque et lorsque . Représenter alors l'allure de la tension en sortie de l'onduleur. Quelle est la tension efficace de ?
D2. Rappeler l'expression de l'impédance d'une inductance ainsi que celle, , d'une résistance . En déduire que l'expression de la fonction de transfert est
où l'expression de est à déterminer en fonction de et .
D3. Quelle est la nature de ce filtre ? Est-il adapté pour filtrer en un signal sinusoïdal à 50 Hz ? Déterminer le gain du filtre. Exprimer la fréquence de coupure à -3 dB de ce filtre, notée , en fonction de .
D4. En sortie du filtre, quels sont les rangs des harmoniques présents dans le signal ? L'amplitude de l'harmonique de rang trois est-elle négligeable devant celle du fondamental? Même question avec l'harmonique de rang cinq d'amplitude . Commenter.
D5. Commenter et justifier la façon dont a été définie la fonction p dans le code ci-dessus.
D6. Dans cette nouvelle séquence, on montre que si et sont fermés alors , si et sont fermés alors et sinon. On associe alors à chaque interrupteur un entier valant 0 si l'interrupteur est ouvert et 1 s'il est fermé. Élaborer une fonction Python nommée tension prenant en arguments les quatre valeurs des interrupteurs et et renvoyant la tension en sortie de l'onduleur.
D7. Élaborer une fonction Python nommée onduleur prenant comme argument un flottant représentant le temps , temps auquel sont évaluées les conditions d'ouverture de la séquence donnée précédemment à l'équation (8), et renvoyant la valeur de la tension en sortie de l'onduleur à cet instant .
D8. Donner les instructions en langage Python pour construire les deux listes à enregistrer dans les variables respectives Liste_t et Liste_u et contenant respectivement les valeurs des temps et des tensions .
D9. Écrire alors en langage Python les commandes permettant de tracer le graphique de et de se représenter le fonctionnement de l'onduleur comme dans l'exemple page suivante de la Figure 6. On ne se souciera ni des légendes ni des axes.
D10. Quelle est la fréquence d'échantillonnage de ce signal? Justifier qualitativement pourquoi un si grand nombre de points de calcul du signal a été choisi ?
D11. Écrire la ligne de code Python permettant d'importer le module fftpack de la bibliothèque scipy.
D12. À partir du spectre présenté sur la Figure 7, vérifier si la fréquence du fondamental, notée , est compatible avec celle déduite du signal temporel de la Figure 6, et commenter la présence des autres harmoniques.
D13. Calculer les amplitudes relatives des harmoniques et après le filtre. Quelle est l'allure du signal de sortie en sortie de filtre lorsque l'onduleur est à commande décalée ? Quel avantage y-aurait-il à utiliser un onduleur à commande décalée dans ce montage ?
E/ Dimensionnement des câbles
E1. Dans le modèle présenté ci-dessus, appliquer la deuxième loi de Newton à un électron du milieu et donner l'équation différentielle vérifiée par sa vitesse . Faut-il prendre en compte le poids de l'électron? Justifier.
E2. On note la perméabilité du vide et la permittivité du vide. Écrire pour le milieu considéré les quatre équations de Maxwell.
E3. Justifier que la force magnétique subie par un électron est négligeable devant la force électrique. Simplifier alors l'équation du mouvement d'un électron.
E4. Exprimer en fonction de et . En déduire l'expression du vecteur densité de courant en fonction des mêmes variables et de .
E5. Rappeler l'expression de la loi d'Ohm locale en fonction de et de la conductivité électrique . En déduire que la conductivité complexe en régime permanent sinusoïdal s'exprime
et donner l'expression de en fonction de et .
E6. Calculer . Par des calculs d'ordre de grandeurs, simplifier l'expression de la conductivité ainsi que l'équation de Maxwell-Ampère.
E7. Déduire des équations de Maxwell et de la loi d'Ohm locale que le vecteur densité de courant vérifie
Déterminer l'expression de et le calculer numériquement dans le cas étudié. Quel nom donne-t-on couramment à ?
E8. En étudiant les courbes représentatives de pour différents rayons, quelle est la zone du câble la plus sollicitée pour transporter le courant? Est-il utile de fabriquer des câbles dont le rayon vaut plusieurs fois ? Justifier.
E9. Au début du réseau basse tension, la section des câbles est . Commenter.
E10. En pratique, pour faire circuler des courants intenses, on utilise plusieurs fils de faible rayon isolés les uns des autres, plutôt qu'un seul fil de gros rayon. Justifier ce choix.
E11. Proposer une explication permettant à d'être supérieur à 1 pour certaines valeurs de , comme cela peut s'observer sur la Figure 8 dans le cas où .
E12. Expliquer en quelques phrases et avec un schéma le principe d'intégration numérique par la méthode des rectangles.
E13. Élaborer une fonction Python nommée Rectangles prenant comme arguments une fonction , deux flottants et et un nombre de points de calcul et renvoyant la valeur approchée par la méthode des rectangles de calculée sur points.
E14. Quelle est la complexité de calcul de cet algorithme? On donnera la réponse sous la forme où est une grandeur à déterminer en fonction de .
E15. A partir de la Figure 9, estimer l'évolution de l'erreur en fonction de pour des valeurs de pas trop élevées. On donnera la réponse sous la forme où est une grandeur à déterminer en fonction de . Comparer les trois méthodes d'intégration numérique.
E16. Passer une certaine valeur de , on remarque que l'erreur augmente légèrement avec n. Ceci est visible sur la Figure 9 pour la méthode des trapèzes ou de Simpson. Proposer une explication à cela.
F/ Propagation de la tension le long d'une ligne électrique
F1. L'ordre de grandeur de est de 50 km . L'Approximation des Régimes QuasiStationnaire (ARQS) est-elle valable en tout point du câble? Justifier. On rappelle que la fréquence du signal est .
F2. Montrer que l'application des lois de l'électrocinétique dans l'ARQS permet d'établir les deux équations suivantes
F3. En déduire que vérifie une équation de D'Alembert et exprimer la célérité des ondes associées en fonction de et .
F4. Déterminer la relation de dispersion des ondes dans ce milieu entre la pulsation et le nombre d'onde noté . Donner la définition d'un milieu dispersif et d'un milieu absorbant. Le milieu est-il dispersif? Est-il absorbant? Les réponses devront être succinctement justifiées.
F5. Donner la forme des OPPH incidentes et réfléchies en tension. En utilisant les relations de l'équation (11), exprimer l'impédance du milieu pour les ondes incidentes, notée , et pour les ondes réfléchies , en fonction de et . Quelle est l'unité des impédances et ?
F6. Écrire, sans les résoudre, les deux équations correspondantes aux conditions de bord et permettant d'obtenir et .
F7. Déterminer l'expression de en fonction de et afin d'annuler l'onde réfléchie. Qualitativement, quel est l'avantage de cette situation d'un point de vue énergétique?
F8. Quelles sont les deux nouvelles équations différentielles couplées et du premier ordre en et en vérifiées par et ?
F9. On considère des solutions de la forme et . Montrer que est solution de l'équation différentielle
et .
F10. Sans chercher à déterminer les expressions de et de , préciser si ce milieu est dispersif et/ou absorbant.
F11. Quels sont les avantages à choisir les paramètres des câbles tels que ?
F12. Calculer numériquement lorsque . En l'absence d'onde réfléchie, calculer la variation relative de l'amplitude de la tension aux bornes de lorsque l'on prend en compte les pertes dans la ligne par rapport au cas sans perte. Commenter.
Fin de l'épreuve
E3A Modélisation PC 2017 - Version Web LaTeX | WikiPrépa | WikiPrépa