Brevet FR3155308A1 — publié le 20250516

Le procédé comprend une étape d’acquisition de données, une étape d’estimation d’une distribution thermique dans la cellule (MP) et une étape de calcul de gradient thermique à partir de la distribution thermique, l’étape d’estimation comprenant des calculs basés sur une application de l’équation de la chaleur de Fourier, avec la méthode dite « des volumes finis » et un maillage (MD) modélisant un enroulement de couches en spirale de la cellule, le maillage comprenant des mailles parallélépipédiques (M i1 , M i3 ) modélisant une partie centrale (PC) et des mailles arrondies (M i2 , M i4 ) modélisant des flancs latéraux (FG, FD) de l’enroulement. Conformément à l’invention, les mailles arrondies d’un flanc et les mailles arrondies de l’autre flanc ont des rayons de courbure différents (r max , r min ), une maille arrondie d’un flanc ayant, dans une même couche de l’enroulement, un rayon supérieur à celui d’une maille arrondie correspondante de l’autre flanc. Figure 5 PROCÉDÉ D’ESTIMATION DE LA TEMPÉRATURE D’UNE CELLULE DE STOKAGE ELECTRIQUE PRISMATIQUE DANS UNE BATTERIE ELECTRIQUE DE TYPE LITHIUM-ION La présente invention concerne de manière générale le domaine de la gestion thermique des batteries électriques de type lithium-ion incorporant des cellules de stockage électrique prismatiques. Plus particulièrement, l’invention se rapporte à un procédé d’estimation de la température dans une cellule de stockage électrique prismatique de type lithium-ion. Le procédé de l’invention est applicable notamment pour la gestion thermique des batteries électriques de traction de type lithium-ion dans des véhicules électrifiés. Dans un véhicule électrifié, tel qu’un véhicule tout électrique ou un véhicule hybride thermique-électrique, la batterie électrique de traction peut dégager une puissance calorifique conséquente lorsqu’elle est sollicitée pour la traction du véhicule. Une gestion thermique de la batterie électrique de traction est nécessaire pour maximiser sa durée de vie, sa disponibilité et ses performances. La connaissance de la température des cellules de stockage électrique de la batterie, plus précisément, de la distribution spatiale de la température dans celles-ci, permet de calculer des gradients thermiques internes dans ces cellules. Ces gradients thermiques internes sont utilisés pour la gestion thermique de la batterie, via notamment la commande de son système de refroidissement. Ils peuvent aussi être utilisés pour l’activation de modes de traction dégradés du véhicule assurant une moindre sollicitation de la batterie, pour une réduction de ses temps de recharge via une commande optimisée de la puissance électrique de charge, ainsi que pour l’anticipation des pannes et la prévention de l’emballement thermique de la batterie. Des procédés d’estimation de la température dans une cellule de stockage électrique de type lithium-ion sont connus dans l’état de la technique. Ces procédés font appel une modélisation thermique de la cellule pour obtenir une estimation de la température interne de celle-ci. Une approche classique de modélisation consiste à considérer la cellule comme un matériau homogène et isotrope. Il en résulte un modèle simple de la cellule, mais qui ne permet pas l’obtention des gradients thermiques internes dans cette dernière, en raison de l'hypothèse d'isotropie. Des modèles plus élaborés prennent en compte l'anisotropie des matériaux et autorisent l’obtention d’une distribution spatiale de la température dans la cellule, mais ces modèles restent trop simples et manquent de précision pour une estimation fiable des gradients thermiques internes. Par ailleurs, il est connu des modèles complexes basés sur l'analyse par éléments finis, dite « FEA » pour « Finite Element Analysis » en anglais, qui procurent une définition détaillée de la structure de la cellule et permettent d’atteindre une précision élevée de la distribution spatiale de la température dans la cellule. Avec de tels modèles, il est possible d’estimer valablement les gradients thermiques internes dans la cellule, mais leur complexité demande une capacité de calcul importante, ce qui limite considérablement l’utilité de ces modèles pour certaines applications en temps réel, comme notamment la gestion thermique des batteries électriques de traction dans les véhicules électrifiés. Le document CN109902372A décrit un procédé d’estimation de la température du type ci-dessus pour une cellule lithium-ion prismatique. Le procédé de CN109902372A exige une modélisation tridimensionnelle (3D) via des logiciels de conception assistée par ordinateur, comme « AutoCAD® » ou « SolidWorks® », avec une discrétisation minutieuse de l’enroulement de couches en spirale de la cellule, dit « jelly roll » en anglais. Il requiert des ressources de calcul considérables qui ne sont pas toujours disponibles, notamment dans un système de gestion de batterie de véhicule, dit « BMS » pour « Battery Management System » en anglais. Dans la littérature, il est décrit des modèles optimisés pour des cellules cylindriques, peu gourmands en temps de calcul, qui prennent en compte l'anisotropie des matériaux et sont capables d'évaluer un gradient thermique interne dans la cellule. Dans l’article « Aging effect of temperature gradients in Li-ion cells experimental and simulative investigations and the consequences on thermal battery management » , présenté lors du symposium « 26th Electric Vehicle Symposium and Exposition (EVS26), Los Angeles, California, May 6 - 9, 2012 » , les auteurs Matthias Fleckenstein et al. divulguent un procédé d’évaluation du gradient thermique interne dans cellule prismatique lithium-ion. Le procédé fait appel à la méthode des volumes finis, en utilisant l'équation de la chaleur de Fourier, et propose pour le calcul un maillage de discrétisation de la cellule prismatique en 4x5x6 mailles élémentaires. Le maillage de discrétisation décrit dans l’article susmentionné, illustré à la FIG. 1 , repose sur une représentation géométrique approximative CE d’une cellule prismatique sous la forme de trois couches en boucle superposées CS1, CS2 et CS3. L’enroulement des couches en spirale les unes sur les autres n’est pas pris en compte. L’enroulement en spirale introduit une variation du rayon de courbure dans une même couche. Le maillage de discrétisation proposé dans l’article susmentionné résulte en un manque de précision dans le calcul de la distribution spatiale de la température. De plus, le transfert de la chaleur n’est considéré ici que dans la direction radiale (cf. la flèche DR à la FIG. 1 ), entre des mailles centrales MC et des mailles externes ME. Il est souhaitable de fournir une solution ne présentant pas les inconvénients susmentionnés de la technique antérieure et adaptée pour une intégration dans un véhicule électrifié, pour estimer de manière fiable la température et les gradients thermiques dans des cellules de stockage électrique prismatique de type lithium-ion. Selon un premier aspect, l’invention concerne un procédé d’estimation de la température dans une cellule de stockage électrique prismatique de type lithium-ion comprenant une étape d’acquisition de données, une étape d’estimation d’une distribution spatiale thermique dans la cellule de stockage électrique prismatique et une étape de calcul de gradient thermique dans la cellule de stockage électrique prismatique à partir de la distribution spatiale thermique, l’étape d’estimation de distribution spatiale thermique comprenant des calculs basés sur une application de l’équation de la chaleur de Fourier, avec la méthode dite « des volumes finis » et un maillage de discrétisation prédéterminé modélisant un enroulement de couches en spirale dit « jelly roll » de la cellule de stockage électrique prismatique, le maillage de discrétisation comprenant une pluralité de mailles parallélépipédiques modélisant une partie centrale de l’enroulement et une autre pluralité des mailles arrondies modélisant des premier et deuxième flancs latéraux arrondis distants de l’enroulement. Conformément à l’invention, les mailles arrondies modélisant le premier flanc et les mailles arrondies modélisant le deuxième flanc ont des rayons de courbure différents, une maille arrondie du premier flanc ayant, dans une même couche de l’enroulement, un rayon de courbure supérieur à celui d’une maille arrondie correspondante du deuxième flanc, de façon à modéliser fidèlement l’enroulement en spirale d’une couche sur l’autre. Selon une caractéristique particulière, l’étape d’estimation de distribution spatiale thermique comprend des calculs tridimensionnels en temps réel, à partir de la distribution spatiale thermique, pour estimer et suivre une évolution d’un ou plusieurs gradients thermiques de la cellule de stockage électrique prismatique dans un espace à trois dimensions. Selon une autre caractéristique particulière, l’étape d’estimation de distribution spatiale thermique comprend une prise en compte automatique de paramètres physico-chimiques et de caractéristiques thermo-physiques des couches. L’invention concerne aussi un calculateur comprenant une mémoire stockant des instructions de programme pour la mise en œuvre du procédé tel que décrit brièvement ci-dessus. L’invention concerne aussi un ensemble d’une batterie électrique comprenant une pluralité de cellules de stockage électrique prismatique de type lithium-ion et d’un calculateur, intégré dans un véhicule, le calculateur étant un calculateur tel que décrit ci-dessus. Dans une forme de réalisation particulière pour un véhicule électrifié, la batterie électrique est une batterie électrique de traction. Dans une autre forme de réalisation particulière, le calculateur compris dans cet ensemble est un calculateur de gestion de batterie du type dit « BMS ». L’invention concerne aussi un véhicule comprenant un ensemble tel que décrit brièvement ci-dessus. D’autres avantages et caractéristiques de la présente invention apparaîtront plus clairement à la lecture de la description détaillée ci-dessous de plusieurs modes de réalisation particuliers de l’invention, en référence aux dessins annexés, dans lesquels : La FIG. 1 est une représentation tridimensionnelle d’un maillage de discrétisation de la technique antérieure d’une cellule de stockage électrique prismatique de type lithium-ion. La FIG. 2 est un bloc-diagramme montrant schématiquement de manière simplifiée une architecture électrique d’un véhicule électrique dans lequel est mis en œuvre un mode de réalisation particulier du procédé selon l’invention. La FIG. 3 est une représentation tridimensionnelle d’un maillage de discrétisation pour une cellule de stockage électrique prismatique de type lithium-ion, utilisé dans le procédé de l’invention. La FIG. 4 est un bloc-diagramme montrant des blocs fonctionnels du procédé selon l’invention. La FIG. 5 est une vue plane de face du maillage de discrétisation de la FIG. 3 montrant les rayons de courbure distincts, pris en compte par le procédé de l’invention, de mailles arrondies distantes dans une même couche de l’enroulement de couches en spirale d’une cellule de stockage électrique prismatique de type lithium-ion. La FIG. 6 est une représentation schématique d’un ensemble formé d’une cellule de stockage électrique prismatique de type lithium-ion et d’une plaque de refroidissement. La FIG. 7 montre, à titre d’exemple, une courbe de courant de charge/décharge et une courbe correspondante d’un gradient thermique estimé dans l’ensemble de la FIG. 6 . En référence aux Figs.2 à 5, il est maintenant décrit un mode de réalisation particulier du procédé selon l’invention. Dans cet exemple de réalisation, le procédé est mis en œuvre dans un véhicule électrifié sous la forme d’un véhicule tout électrique VE montré à la FIG. 2 . Le véhicule VE comprend une chaîne de traction électrique intégrant un groupe motopropulseur électrique eGMP. Le groupe motopropulseur électrique eGMP est alimenté en énergie par un pack batterie haute tension BAT_HV, à travers un convertisseur électrique réversible de type DC/AC (non représenté). Le pack batterie BAT_HV comprend ici une pluralité de cellules de stockage électrique prismatique lithium-ion MP intégrées typiquement dans des boitiers modulaires parallélépipèdiques, qui stockent l’énergie électrique, ainsi que des bus conducteurs de connexion électrique des cellules et un système de gestion thermique (non représentés). Le pack batterie BAT_HV est géré par un calculateur BMS. Le calculateur BMS est relié à un réseau de communication de données eCAN, typiquement de type « CAN », du véhicule et héberge des modules logiciels embarqués qui réalisent différentes fonctions de gestion et de mesure pour le pack batterie BAT_HV, comme par exemple une estimation d’un état de charge dit « SOC » (pour « State Of Charge » en anglais) et une estimation d’un état de santé dit « SOH » (pour « State Of Heath » en anglais). La mise en œuvre du procédé selon l’invention fait appel à un module logiciel embarqué MOD_SW qui est hébergé ici dans le calculateur BMS susmentionnée du véhicule. Comme représenté à la FIG. 2 , le module logiciel MOD_SW est implanté dans une mémoire MEM du calculateur BMS. Le calculateur BMS peut coopérer, sous la supervision du module logiciel MOD_SW, avec un ou plusieurs autres calculateurs du véhicule pour la mise en œuvre du procédé selon l’invention, en particulier ici le calculateur eVCU qui est le calculateur superviseur de la chaîne de traction électrique du véhicule. Le module logiciel MOD_SW autorise la mise en œuvre du procédé selon l’invention par l’exécution d’instructions de code de programme par un processeur (non représenté) du calculateur BMS. Comme montré schématiquement à la FIG. 2 , pour la mise en œuvre du procédé selon l’invention, le module logiciel MOD_SW exploite différentes informations INF_MD, détaillées plus bas, qui sont typiquement déjà disponibles dans le calculateur BMS, obtenues en provenance d’un autre calculateur, comme le calculateur eVCU, et/ou fournies par des capteurs de mesure du pack batterie BAT_HV. Le module logiciel MOD_SW délivre en sortie des informations de distribution spatiale tridimensionnelle de température D_T, de gradient thermique GD et d’alerte AL. A la FIG. 3 , un exemple d’un maillage de discrétisation MD utilisé par le procédé de l’invention pour une cellule de stockage électrique prismatique de type lithium-ion MP est représenté en trois dimensions (3D) dans un repère spatial orthogonal (X, Y, Z). Comme visible à la FIG. 3 , le maillage de discrétisation MD comprend une pluralité de mailles parallélépipédiques modélisant une partie centrale parallélépipédique PC de l’enroulement de couches en spirale et une pluralité des mailles arrondies modélisant des premier et deuxième flancs latéraux arrondis distants FG et FD de l’enroulement de couches en spirale. Dans l’exemple représenté, la cellule MP est formée d’un enroulement en spirale de quatre couches CC1 à CC4. Le maillage MD comprend quatre tranches T1 à T4 de même hauteur suivant l’axe Z. Dans chaque tranche T1, T2, T3 ou T4, les couches sont représentées chacune par quatre mailles adjacentes correspondantes. Avec les quatre tranches, une couche complète est donc discrétisée en 4x4 mailles. Ainsi, par exemple, comme illustrée à la FIG. 3 , dans la tranche T1, la couche CC4 est discrétisée ici en quatre mailles M 41 , M 42 , M 43 et M 44 . Les mailles M 41 et M 43 ont des formes géométriques parallélépipédiques. Les mailles M 42 et M 44 ont des formes géométriques arrondies avec des rayons de courbures respectifs distincts r max et r min , pour représenter fidèlement l’enroulement en spirale d’une couche sur l’autre. Le rayon de courbure r max de la maille M 42 est supérieur à celui r min de la maille M 44 . Le volume de la maille M 42 est donc supérieure à celui de la maille M 44 . Chacune des couches CC1 à CC4 est formée d’une pluralité de sous-couches, typiquement huit couches, représentées schématiquement dans l’agrandissement AG de la FIG. 3 . Ces sous-couches correspondent à un empilement de feuilles ayant des paramètres physico-chimiques et des caractéristiques thermo-physiques différentes, à savoir notamment, des matériaux, des épaisseurs, des densités, des capacités thermiques et des conductivités thermiques différentes. Ces sous-couches correspondent notamment à une feuille de collecteur en aluminium, une feuille de collecteur en cuivre, deux feuilles de séparateur, deux feuilles d’anode et deux feuilles de cathode. Le processus de traitement mis en œuvre dans le procédé selon l’invention est maintenant décrit ci-dessous en référence plus particulièrement à la FIG. 4 et à la FIG. 5 . Comme représenté schématiquement à la FIG. 4 , le processus de traitement comprend essentiellement quatre blocs fonctionnels B1, B2 à B4. Ces blocs fonctionnels B1 à B4 sont respectivement un bloc d’acquisition de données, un bloc de calcul de distribution spatiale thermique, un bloc d’estimation de gradient thermique et un bloc de restitution de résultats d’estimation et d’alerte. Le bloc fonctionnel B1 a pour fonction d’acquérir et stocker différentes informations nécessaires à la mise en œuvre du procédé de l’invention, comme les informations INF_MD susmentionnées. Ces informations INF_MD comprennent notamment des données de caractéristiques géométriques de l’enroulement de couches en spirale de la cellule MP, des données des caractéristiques thermo-physiques susmentionnées des sous-couches, des données de capacité et puissance électriques, des données de courant de charge/décharge et de tension de la cellule MP, des données d’état de charge « SOC » et d’état de santé « SOH », des données aux limites comme des températures mesurées et/ou estimées de surfaces de contact avec un système de refroidissement, l’air ambiant environnant, etc., et/ou d’autres données, en fonction de l’application. Des données fixes de caractéristiques géométriques, thermo-physiques et autres sont acquises et mémorisées préalablement. Les données variables (courant, tension, température du système de refroidissement, etc.) sont typiquement acquises en temps réel pendant le déroulement du processus de traitement. Le bloc de calcul de distribution spatiale thermique B2 fait appel essentiellement à quatre fonctions principales de calcul F1 à F4. Le bloc B2 délivre en sortie les informations susmentionnées de distribution spatiale tridimensionnelle de température D_T. La fonction F1 utilise l’équation de la chaleur de Fourier EQ1 et la méthode dite « des volumes finis » pour les calculs de température dans les différentes mailles. Les calculs sont effectués dans un repère tridimensionnel montré à la FIG. 5 ayant un axe axial (coordonnée l), un axe radial (coordonnée r) et un axe ortho-radial (coordonnée z). Dans l’équation EQ1, la variable T i,j représente la température dans la maille M ij , ρ ij la densité de la maille M ij et C ij sa capacité thermique massique, Q ij la chaleur totale générée par la maille M ij et V ij son volume, et λ r , λ l et λ z les conductivités thermiques dans les trois axes considérés et S r , S l et S z les surfaces de contact entre mailles dans les trois axes considérés. La fonction F2 utilise les équations EQ2 et EQ3 pour calculer respectivement la densité ρ ij de la maille M ij et sa capacité thermique massique C ij , à partir des caractéristiques thermo-physiques des sous-couches susmentionnées. Dans les équations EQ2 et EQ4, pour la sous-couche considérée, V i représente le volume de la sous-couche dans la maille M ij , ρ ij sa densité, M i sa masse et C pi sa capacité thermique. La fonction F3 est chargée du calcul des surfaces de contact entre les mailles et de leurs volumes. En référence aussi à la FIG. 5 , le calcul pour les mailles arrondies M i2 , M i4 , est détaillé et fait appel aux équations EQ4 à EQ6 d’une sous-fonction F30. Pour les mailles parallélépipédiques M i1 , M i3 , le calcul n’est pas détaillé et est réalisé par une sous-fonction F31. Les équations EQ4 donnent les rayons de courbure des mailles arrondies M i2 , M i4 , pour une couche i, i étant le numéro de la couche sous la forme d’un entier croissant depuis la couche centrale jusqu’à la couche externe de l’enroulement de couches en spirale. Dans l’équation EQ4, rmax i-1/i , rmax i/i+1 , et rmin i-1/i , rmin i/i+1 , représentent respectivement les petits et grands rayons pour les deux mailles arrondies M i2 et M i4 de la couche i, et rmax 0/i et rmin 0/i les petits rayons pour les deux mailles arrondies M 12 et M 14 de la première couche centrale i=1, sachant que d et e sont respectivement l’épaisseur de la couche i et l’épaisseur d’un espace central creux de l’enroulement. Les équations EQ5 sont utilisées pour calculer les surfaces de contact entre mailles S r , S l et S z , dans les trois axes considérés, pour les mailles arrondies M i2 et M i4 de la couche i, avec dz étant la hauteur de la maille. Les volumes V ij pour les mailles arrondies M i2 et M i4 sont données par les équations EQ6. La fonction F4 est chargée de calculer la chaleur totale Q ij produite dans la maille considérée M ij . Un modèle électrique est utilisé ici pour évaluer la production de chaleur. Le bloc B3 reçoit les informations de distribution spatiale tridimensionnelle de température D_T fournit par le bloc B2 et estime en temps réel et en trois dimensions (3D) un ou plusieurs gradients thermiques GD, selon le besoin et/ou l’application. Les informations de distribution spatiale de température D_T et de gradient thermique GD sont délivrées en entrée au bloc B4. Le bloc B4 est chargé de la restitution des informations estimées D_T et GD à un ou plusieurs systèmes destinataires et de la mise en œuvre d’une ou plusieurs stratégies préprogrammées d’alertes AL qui sont basées sur l’exploitation des informations D_T, GD. Selon l’application, les informations D_T, GD et AL seront destinées à une ou plusieurs stratégies de pilotage de système (système de refroidissement de batterie, système de gestion de batterie, groupe motopropulseur électrique, borne de recharge électrique ou autres), à un système de maintenance préventive ou autres. Des essais et simulations ont été réalisés par l’entité inventive et ont permis de vérifier la précision des estimations obtenues grâce au procédé de l’invention. En référence à la FIG. 6 et à la FIG. 7 , à titre illustratif, il a été réalisé une simulation d’estimation d’un gradient thermique GD dans une cellule MP montée sur une plaque de refroidissement RE. Le gradient thermique GD est estimé pour des conditions aux limites données et un courant de charge/décharge Ic ayant un profil spécifique. Le profil du courant Ic montré à la FIG. 7 a été déterminé de façon à étudier l'effet de celui-ci sur le gradient thermique GD pendant la charge et la décharge de la cellule MP. La cellule MP est soumise à une conduction thermique par une face inférieure en contact avec la plaque de refroidissement RE. Le gradient thermique GD considéré ici est celui entre une section inférieure de la cellule MP proche de la plaque RE et une section supérieure distante de la cellule MP. Pour cette simulation, la température de l'air ambiant est fixée à 25°C et la température de refroidissement est réglée à 30°C. La FIG. 7 montre le gradient thermique GD observé dans la cellule MP pendant une phase de décharge DCH et une phase de recharge ultérieure CH. Le gradient thermique GD atteint GD = 0,7 °C environ pendant la phase de charge CH et est plus élevé que celui observé GD = 0,2 °C pendant la phase de décharge DCH. Cela est dû principalement au fait que le courant Ic est plus important pendant la phase de recharge CH, à savoir, de 60 A environ ici alors qu’il n’est que de 20 A environ pendant la phase de décharge DCH. Cette simulation fait apparaître que la différence de température entre les sections supérieure et inférieure de la cellule MP augmente sensiblement à des taux de charge élevés, en particulier lorsque la température de la plaque de refroidissement est basse. Le procédé de l’invention en procurant une estimation précise des gradients thermiques autorise une gestion thermique efficace de la cellule MP, de façon à maintenir les performances de celle-ci et à répondre aux exigences sécuritaires. Le procédé de l’invention, contrairement aux procédés de l’état de la technique basés sur des modèles complexes, est frugal en ressources de calcul et peut être mis en œuvre dans le système de gestion de batterie d’un véhicule électrifié. Le maillage utilisé dans le procédé de l’invention offre l’avantage d’être paramétrique, ce qui le rend adaptable à différents enroulements de couches en spirale des cellules prismatiques de type lithium-ion. Il permet un enroulement des couches élémentaires les unes sur les autres, facilitant ainsi le passage de la chaleur d'une couche à une autre. Par ailleurs, le maillage du procédé de l’invention est flexible et permet de choisir le nombre de mailles selon trois axes (axe axial, axe radial et ortho-radial), ce qui apporte la possibilité d'optimiser le temps de simulation pour un suivi en temps réel de l'évolution du gradient thermique. On notera aussi que les paramètres physico-chimiques des couches élémentaires, ainsi que les paramètres du modèle électrothermique, peuvent être calculés automatiquement par le processus. L’invention ne se limite pas au mode de réalisation particulier qui a été décrit ici à titre d’exemple. L’homme du métier, selon les applications de l’invention, pourra apporter différentes modifications et variantes entrant dans le champ de protection de l’invention. Procédé d’estimation de la température dans une cellule de stockage électrique prismatique de type lithium-ion (MP) comprenant une étape d’acquisition de données (B1), une étape d’estimation d’une distribution spatiale thermique (B2) dans ladite cellule de stockage électrique prismatique (MP) et une étape de calcul de gradient thermique (B3) dans ladite cellule de stockage électrique prismatique (MP) à partir de ladite distribution spatiale thermique, ladite étape d’estimation de distribution spatiale thermique (B2) comprenant des calculs (F1) basés sur une application de l’équation de la chaleur de Fourier, avec la méthode dite « des volumes finis » et un maillage de discrétisation prédéterminé (MD) modélisant un enroulement de couches en spirale dit « jelly roll » de ladite cellule de stockage électrique prismatique (MP), ledit maillage de discrétisation (MD) comprenant une pluralité de mailles parallélépipédiques (M i1 , M i3 ) modélisant une partie centrale (PC) dudit enroulement et une autre pluralité des mailles arrondies (M i2 , M i4 ) modélisant des premier et deuxième flancs latéraux arrondis distants (FG, FD) dudit enroulement, caractérisé en ce que lesdites mailles arrondies (M i2 ) modélisant ledit premier flanc (FG) et lesdites mailles arrondies (M i4 ) modélisant ledit deuxième flanc (FD) ont des rayons de courbure différents, une maille arrondie (M i2 ) dudit premier flanc (FG) ayant, dans une même couche (i) dudit enroulement, un rayon de courbure (r max ) supérieur à celui (r min ) d’une maille arrondie correspondante (M i4 ) dudit deuxième flanc (FD), de façon à modéliser fidèlement l’enroulement en spirale d’une couche sur l’autre. Procédé selon la revendication 1, caractérisé en ce que ladite étape de calcul de gradient thermique (B3) comprend des calculs tridimensionnels en temps réel, à partir de ladite distribution spatiale thermique, pour estimer et suivre une évolution d’un ou plusieurs gradients thermiques (GD) de ladite cellule de stockage électrique prismatique (MP) dans un espace à trois dimensions. Procédé selon la revendication 1 ou 2, caractérisé en ce que ladite étape d’estimation de distribution spatiale thermique (B2) comprend une prise en compte automatique de paramètres physico-chimiques et de caractéristiques thermo-physiques (F2) desdites couches. Calculateur (BMS) caractérisé en ce qu’il comprend une mémoire (MEM) stockant des instructions de programme (MOD_SW) pour la mise en œuvre du procédé selon l'une quelconque des revendications 1 à 3. Ensemble d’une batterie électrique (BAT_HV) comprenant une pluralité de cellules de stockage électrique prismatique de type lithium-ion (MP) et d’un calculateur, intégré dans un véhicule (VE), caractérisé en ce que ledit calculateur est un calculateur (BMS) selon la revendication 4. Ensemble selon la revendication 5, intégré dans un véhicule électrifié (VE), caractérisé en ce que ladite batterie électrique est une batterie électrique de traction (BAT_HV) dudit véhicule électrifié (VE). Ensemble selon la revendication 5 ou 6, caractérisé en ce que ledit calculateur est un calculateur de gestion de batterie du type dit « BMS ». Véhicule (VE) caractérisé en ce qu’il comprend un ensemble (BAT_HV, BMS) selon l’une quelconque des revendications 5 à 7.