FR 3069583 A1 20190201 FR 1757152 A 20170727 DOMAINE DE L'INVENTION L'invention concerne un procédé et un dispositif de détection d'allumage d'une chambre à combustion de moteur-fusée, un procédé de démarrage de 5 moteur-fusée, et un moteur-fusée. ARRIERE-PLAN DE L'INVENTION Dans un moteur-fusée, le démarrage est une phase critique : l'allumage du moteur (c'est-à-dire la combustion des ergols) peut ne pas se produire, ce 10 qui conduit le plus souvent à l'échec de la mission du lanceur-fusée ; de manière alternative, si le moteur s'allume trop tard, alors qu'une grande quantité d'ergols s'est accumulée sous forme imbrûlée dans la chambre de combustion, leur inflammation soudaine peut susciter une explosion qui peut endommager le moteur voire le détruire. Il est donc important de pouvoir 15 déterminer de manière précise si la mise à feu d'un moteur a bien eu lieu et à quel moment elle se produit. Il existe donc un besoin pour un procédé et un dispositif permettant de déterminer, de la manière la plus fiable possible, si un moteur-fusée est allumé. 20 OBJET ET RESUME DE L'INVENTION Un premier objectif de l'invention est de satisfaire ce besoin en proposant un procédé de détection d'allumage d'une chambre à combustion de moteur-fusée qui permette de détecter de manière la plus sûre possible si 25 un moteur-fusée est allumé ou non. Cet objectif est atteint grâce à un procédé de détection d'allumage d'une chambre à combustion de moteur-fusée, le procédé comprenant les étapes suivantes, réalisées de manière itérative au moyen d'un calculateur : A) pour chacun des ergols, on acquiert des valeurs de grandeurs qui 30 permettent de calculer des valeurs moyennes PGC_HO, PGC_Hl et des écarts-types a_PGCHO, o_PGCH1 d'une pression théorique d'allumage et d'une pression théorique de non-allumage dans la chambre de combustion, qui sont les pressions dans la chambre de combustion respectivement dans l'hypothèse de non-allumage (HO) qui correspond au cas où la combustion 35 n'a pas lieu et dans l'hypothèse d'allumage (H1) qui correspond au cas où la combustion a lieu ; B) en fonction des valeurs acquises à l'étape A), on calcule les valeurs moyennes PGC_HO, PGC_H1 et les écarts-types cr_PGCHO, cr_PGCH1 de la 2 pression théorique d'allumage et de la pression théorique de non-allumage dans la chambre de combustion ; C) on détermine un statut d'allumage du moteur ; l'étape C) comprenant les étapes : Cl) on détermine des coefficients A et B permettant de déterminer les pressions limites d'allumage et de non-allumage, à partir de la probabilité de non-détection de l'allumage Pnd et de la probabilité de fausse alarme de détection d'allumage Pf, à l'aide des équations suivantes : Pnd 1 - Pnd (1) A= et B = 1 - P f P f C2) on attribue au statut d'allumage du moteur (S) une valeur allumé, non allumé ou indéterminé, indiquant respectivement que l'on considère, pour un instant considéré, que le moteur est allumé, qu'il n'est pas allumé, ou que l'on ne sait pas ; le statut d'allumage (S) étant obtenu en calculant des bornes d'allumage/non-allumage Al et B1 et une somme / à l'aide des équations suivantes à partir des coefficients A et B : ( cfpccHopccHi cf (2) A1= 21nA + 2n ln Il / )2) PGCHo N2 f V''GCi - PGCH0 ) VGC'i - pG7HiCHi (3) E = I ( - cr2 Cf 2 i=n-p (4) B1 = 2InB + 2n ln rGcHi) 0-PGCHO dans lesquelles : - PGCi est la pression mesurée dans la chambre de combustion à l'instant i correspondant à une boucle de calcul en cours ; - n est l'indice de la boucle de calcul à l'instant considéré ; - p est un entier prédéterminé ; et en affectant au statut d'allumage (S) une valeur parmi des possibilités Non-allumé, Allumé, Indéterminé, en comparant la somme / aux bornes d'allumage/non-allumage Al et B1 : - Si E B1 : on fixe la valeur du statut d'allumage (S) à la valeur « allumé ». Le procédé repose sur le test de Wald. Dans le cadre de l'invention, ce test est mis en oeuvre de manière à déterminer, avec une confiance 5 prédéterminée, la valeur du statut d'allumage S du moteur. Ce statut indique si la combustion a lieu, si la combustion n'a pas lieu, ou si l'on est dans une zone d'indétermination. Dans un mode de réalisation préférentiel, les grandeurs dont les valeurs sont acquises à l'étape A) incluent, pour chaque ergol, la pression et 10 la masse volumique de l'ergol au point d'injection dans la chambre de combustion. Cependant, d'autres valeurs peuvent éventuellement être acquises plutôt que celles-ci, à condition que ces autres valeurs permettent de calculer les pressions moyennes (PGC_HO, PGC_H1) et les écarts-types (a_PGCHO, a_PGCH1) de la pression théorique d'allumage et de non- 15 allumage dans la chambre de combustion. La détection d'allumage repose donc sur la démarche suivante : A. Génération d'hypothèses de fonctionnement On commence par évaluer deux hypothèses, correspondant aux deux 20 situations possibles : O L'hypothèse HO traduit le domaine de non-allumage de la chambre O L'hypothèse H1 traduit le domaine d'allumage. Pour chacune de ces hypothèses HO ou H1, on suppose que la 25 pression dans la chambre de combustion est une pression théorique caractérisée par une valeur moyenne et un écart type, en faisant l'hypothèse que cette pression théorique suit une loi normale. La probabilité Pk(xi) que la pression soit égale à la valeur xi, dans l'hypothèse k, est donc donnée par l'équation : 30 Pk(xi) = g_ e v211crk (xi-rnk) 2 1. 2 2 0-k dans cette formule : i est un indice correspondant à l'instant d'échantillonnage, xi est la valeur de la pression à un instant ti, m est la valeur moyenne de la pression, a est l'écart type de la pression, k est l'indice indiquant, pour les paramètres moyenne m et écart-type 4 a, quelle est l'hypothèse considérée, à savoir HO ou Hl. Ainsi : si k=0, alors mo = PGC_HO et ao = a_PGCHO ; si k=1, alors m1 = PGC_Hl et cri- a_PGCH1. Pour chacune de des hypothèses HO ou H1, la valeur moyenne 5 (PGC_HO, PGC_H1) et l'écart-type de la pression dans la chambre de combustion (a_PGCHO, o-_PGCH1) varient selon différents paramètres : o des paramètres liés aux conditions de fonctionnement, comme les débits d'alimentation (qui peuvent être estimés par exemple à travers les pressions d'injection) ; et 10 o des paramètres prédéterminés, comme le domaine d'incertitude des capteurs et le domaine d'incertitude des caractéristiques thermodynamiques B. Comparaison de ces hypothèses avec les valeurs mesurées de pression 15 dans la chambre de combustion afin de statuer sur la conformité du signal à l'une ou l'autre des situations ou bien un statut d'incertitude temporaire. On compare les valeurs de pression mesurées à l'une et l'autre des deux hypothèses en utilisant le test de Wald de la manière indiquée précédemment. Cette comparaison fait intervenir le rapport de vraisemblance 20 qui est : n P1(Xn) Fi P1(xi) Vn = - PO(Xn) P0(xi) i=n-p Dans cette équation, n est l'indice de l'itération en cours, Xn est l'évènement défini comme l'obtention de la suite spécifique des valeurs de pression xi mesurées lors des instants ti pour des itérations successives de n-p à n de l'algorithme. La fenêtre temporelle utilisée inclut les itérations n-p 25 à p et est donc glissante au cours du temps. Dans le cas présent, le rapport de vraisemblance s'écrit sous la forme de la variable définie Z plus haut par l'équation (3), et le test de Wald consiste à comparer la valeur de Z aux valeurs des paramètres Al et Bi. Si le test ne permet pas de conclure sur l'état allumé/non allumé de la 30 chambre de combustion et que l'on est donc dans l'incertitude, à l'itération suivante de l'algorithme la fenêtre d'étude est agrandie (en augmentant la valeur de p) suivant jusqu'à ce que l'on puisse statuer sur la situation dans laquelle on se trouve. Dès qu'au cours d'une itération de l'algorithme on parvient à statuer 35 sur la situation dans la chambre de combustion (allumée/non-allumée), on réinitialise la largeur de fenêtre en donnant à celle-ci une valeur po définie 5 par défaut pour l'algorithme. Le paramètre 'p' ou largeur de fenêtre correspond aux nombres de valeurs successives des paramètres de conditions de fonctionnement qui sont prises en compte pour réaliser le test de Wald. Il correspond donc à la durée 5 de la période de temps sur laquelle on se fonde pour déterminer si la situation est une situation d'allumage ou de non-allumage. L'entier 'p' est normalement choisi égal inférieur ou égal à une valeur maximale Pmax qui est telle que la pression dans la chambre de combustion avant l'itération d'indice n-Pmax n'apporte plus d'information quant à l'état 10 allumé/non allumé de la chambre à l'itération n. Dans un mode de mise en oeuvre, à l'étape B), la valeur moyenne PGC_Hl de la pression théorique d'allumage dans la chambre de combustion dans l'hypothèse où la combustion a lieu, est calculée à partir du débit total des ergols injectés dans la chambre de combustion, de la vitesse 15 caractéristique des gaz de combustion et de la section efficace du divergent de la tuyère du moteur-fusée. De préférence, cette valeur moyenne PGC_H1 de la pression théorique d'allumage est calculée grâce à la relation suivante : (Q0+QH)*C* QtotC* 20 (5) PGCHi = Seff Se f f Dans un mode de mise en oeuvre, lors de l'étape A, on acquiert de plus une pression dans la chambre de combustion ; et à l'étape B), l'écart-type a PCG_HO de la pression théorique dans la chambre de combustion 25 dans l'hypothèse (HO) de non-allumage de la chambre est donné par la formule suivante : (6) aPGC_HO = aPGCmes OÙ GPGCmes représente l'écart-type de la mesure de pression dans la chambre 30 de combustion. Cette mesure de pression peut être unique, et dès lors l'incertitude apGcrnesest l'incertitude de la mesure. Dans le cas où la mesure de la pression dans la chambre de combustion est faite de manière redondante, c'est-à-dire par exemple sous la forme de deux mesures PGC1 et PGC2, 35 l'incertitude a PGC peut être donnée par la formule : 6 1 GPGCmes = PGC1+G2PGC --i- '-' n22 V 2 Dans un mode de mise en oeuvre, lors de l'étape A, on acquiert de plus une valeur de la pression atmosphérique PA-rm à l'extérieur du moteur-fusée ; à l'étape B), en fonction des pressions et masses volumiques 5 calculées à l'étape A), on calcule des débits de chacun des ergols ; et l'écart- type a PCG_H1 de la pression théorique dans la chambre de combustion en cas d'allumage de la chambre (hypothèse H1) est calculé en réalisant les étapes suivantes : . on calcule le rapport de mélange RM du moteur-fusée et son écart-type aRm 10 par l'équation suivante : \i(-)0 (P01 - Patm) RM = \I o& (PHI - Patm) Par ailleurs, l'écart-type aRm du rapport de mélange RM du moteur-fusée vérifie : 1 CrRM = CrQ0/QH = + Rm2012H) QH2 Q° Par suite, on calcule l'écart-type aRm du rapport de mélange RM du moteur- 15 fusée par l'équation suivante : 1 (00 1 (1301 -Patm)2 n2 ) eRM = 0-2 (g) ± (PHI - Patm) ( P°I (PHI - Patm)2 -PHI \IV H . on calcule alors les écarts-types a_C* de la vitesse caractéristique C* et a Quir du débit total des ergols injectés dans la chambre de combustion par les 20 équations suivantes : ac. = f(aRm) 1 p '2 Q0 aQtot = CrQH+QO = ,2 avec pour chaque ergol 0-Qx = Ti -k GPx/ ' _i_ "'QUI Dans ces équations : 25 x représente un ergol quelconque, par exemple l'oxygène ou l'hydrogène, apoi et apHI désignent les écarts-types des pressions POA et PHI des deux 7 ergols à leurs points d'injection, kx le coefficient de perte de charge kI pour l'ergol x, aQx et apx1 , pour l'ergol x, l'écart type respectivement d'un débit et d'une pression de l'ergol x au point d'injection, 5 le rapport de mélange RM des deux ergols est égal au rapport entre le débit de comburant sur le débit de carburant à l'injection dans la chambre de combustion QO, QH et Qtot sont respectivement un débit de comburant, de carburant, et un débit cumulé des deux ergols, dans la chambre de combustion ; 10 C* est une vitesse caractéristique des gaz de combustion en sortie de la chambre de combustion ; f est une fonction prédéterminée ; . on calcule ensuite l'écart-type a PGC_H1 de la pression théorique dans la chambre de combustion, en cas d'allumage de la chambre au moyen de la 15 formule : o- PGC H1 = C*2cr2tot Qtot2a2. et par conséquent : Seff Q p \ 2 f p \ (C*(RM)) Wez) Qo aP°I) Qji aPHI) (+ i00 (P01- Patm) + JOH (PHI - Patm) ab 2 1 aPGC_Hl = Seff Dans les formules ci-dessus, le coefficient de perte de charge k1 pour un ergol I est donné par l'équation : 20 ki Q2 / p API où Q est le débit dans le circuit d'alimentation en ergol considéré, et API est la différence de pression entre le point d'injection en ergol pour l'ergol considéré et la chambre de combustion. 25 L'invention concerne également un procédé de démarrage d'un moteur-fusée, dans lequel lors du démarrage d'au moins un moteur du moteur-fusée, on détecte si l'allumage du moteur a bien eu lieu en exécutant le procédé de détection d'allumage tel que défini précédemment. Dans un mode particulier de réalisation, les différentes étapes du 30 procédé de détection d'allumage d'une chambre à combustion de moteur-fusée ou de démarrage d'un moteur-fusée selon l'invention sont déterminées par des instructions de programmes d'ordinateurs. 8 En conséquence, l'invention vise aussi un programme d'ordinateur sur un support d'informations, en particulier un support d'informations non volatiles, ce programme étant susceptible d'être mis en oeuvre par ordinateur, et comportant des instructions de code de programme pour 5 l'exécution des étapes de l'un des procédés définis précédemment, lorsque le programme est exécuté sur un ordinateur. Ce programme peut utiliser n'importe quel langage de programmation, et être sous la forme de code source, code objet, ou de code intermédiaire entre code source et code objet, tel que dans une forme partiellement compilée, ou dans n'importe 10 quelle autre forme souhaitable. L'invention vise aussi un support d'enregistrement, en particulier un support d'enregistrement non volatile, lisible par un ordinateur sur lequel est enregistré un programme d'ordinateur tel que défini précédemment. Le support d'informations ou d'enregistrement peut être n'importe quelle entité 15 ou dispositif capable de stocker le programme. Par exemple, le support peut comporter un moyen de stockage, tel qu'une ROM, par exemple un CD ROM ou une ROM de circuit microélectronique, ou encore un moyen d'enregistrement magnétique, par exemple une clé USB ou un disque dur. Alternativement, le support d'informations ou d'enregistrement peut 20 être un circuit intégré dans lequel le programme est incorporé, le circuit étant adapté pour exécuter ou pour être utilisé dans l'exécution du procédé de détection d'allumage d'une chambre à combustion de moteur-fusée ou de démarrage d'un moteur-fusée. Un deuxième objectif de l'invention est de proposer un dispositif de 25 détection d'allumage d'une chambre à combustion de moteur-fusée qui permette de détecter de manière la plus sûre possible si un moteur-fusée est allumé ou non. Cet objectif est atteint grâce à un dispositif de détection d'allumage d'une chambre à combustion de moteur-fusée, le dispositif comportant : 30 A) une unité d'acquisition configurée pour acquérir, pour chacun des ergols (H,0), des valeurs de grandeurs qui permettent de calculer des valeurs moyennes PGC_HO, PGC_H1 et des écarts-types cr_PGCHO, o_PGCH1 d'une pression théorique d'allumage et d'une pression théorique de non-allumage dans la chambre de combustion, qui sont les pressions dans la chambre de 35 combustion respectivement dans l'hypothèse de non-allumage (HO) qui correspond au cas où la combustion n'a pas lieu et dans l'hypothèse d'allumage (H1) qui correspond au cas où la combustion a lieu ; B) une unité de calcul de valeurs moyennes et d'écarts-types configurée 9 pour, en fonction des valeurs acquises par l'unité d'acquisition, calculer les valeurs moyennes PGC_HO, PGC_Hl et les écarts-types o-_PGCHO, a_PGCH1 de la pression théorique d'allumage et de la pression théorique de non-allumage dans la chambre de combustion ; 5 C) une unité de détermination de statut d'allumage du moteur (56) configurée pour déterminer un statut d'allumage du moteur (S) ; l'unité de détermination de statut d'allumage du moteur (S) étant configurée pour : Cl) déterminer des coefficients A et B permettant de déterminer les pressions limites d'allumage et de non-allumage, à partir de la probabilité de 10 non-détection de l'allumage Pnd et de la probabilité de fausse alarme de détection d'allumage Pf, à l'aide des équations suivantes : Pnd 1-Pnd (1) A = 1-Pf et B = et Pf C2) attribuer au statut d'allumage du moteur (S) une valeur allumé, non 15 allumé ou indéterminé, indiquant respectivement que l'on considère, pour un instant considéré, que le moteur est allumé, qu'il n'est pas allumé, ou que l'on ne sait pas ; le statut d'allumage (S) étant obtenu en calculant des bornes d'allumage/non-allumage Al et I31 et une somme E à l'aide des équations suivantes : (2) Ai= 21nA + 2n ln (aPGcHi) aPGcHo n Z ((PGCL - PGCH0 52 0.pGiGHCiPHG1 2 C11 (PGCi--- 2 i=n-p (3) E = )2 ) (4) B1 = 21nB + 2n In (apcc_H1) 20 où : - PGCi est la pression mesurée dans la chambre de combustion à l'instant (ti) correspondant à la boucle de calcul en cours ; - n est l'indice de la boucle de calcul à l'instant considéré ; - p est un entier prédéterminé (par exemple, p peut être égal au 25 nombre minimal d'itérations de l'algorithme pour lequel la pression dans la chambre de combustion avant l'itération n-p n'apporte plus d'information quant à l'état allumé/non allumé de la chambre à UPGC_HO 10 l'itération n-p) ; et en affectant au statut d'allumage (S) une valeur parmi des possibilités Non-allumé, Allumé, Indéterminé, en comparant la somme E aux bornes d'allumage/non-allumage Al et B1 : 5 - Si E B1 : en fixant la valeur du statut d'allumage (S) à la valeur 10 « allumé ». Dans un mode de réalisation préférentiel, les grandeurs dont l'unité d'acquisition est configurée pour acquérir les valeurs incluent, pour chaque ergol, la pression et la masse volumique de l'ergol au point d'injection dans la chambre de combustion. 15 L'unité d'acquisition, pour acquérir la masse volumique de l'ergol au point d'injection, peut ainsi par exemple être configurée pour acquérir des paramètres de l'ergol et pour calculer la masse volumique de l'ergol au point d'injection à partir des paramètres de l'ergol ainsi acquis. L'invention concerne également un moteur-fusée équipé d'un dispositif 20 de ce type. BREVE DESCRIPTION DES DESSINS L'invention sera bien comprise et ses avantages apparaîtront mieux à la lecture de la description détaillée qui suit, de modes de réalisation 25 représentés à titre d'exemples non limitatifs. La description se réfère aux dessins annexés, sur lesquels : - la figure 1 est une vue schématique d'un moteur-fusée comprenant un dispositif selon l'invention ; - la figure 2 est une vue schématique montrant l'architecture 30 matérielle d'un dispositif selon l'invention, dans un mode de réalisation ; - la figure 3 est un diagramme montrant les étapes du procédé de détection de détection d'allumage de chambre de combustion de moteur fusée selon l'invention. 35 DESCRIPTION DETAILLEE DE L'INVENTION En faisant référence aux figures 1 à 3, un exemple de mise en oeuvre du procédé selon l'invention va maintenant être décrit. Le procédé est présenté dans le cas d'un moteur-fusée consommant 11 comme ergols de l'hydrogène (H) et de l'oxygène (0) ; il peut naturellement être transposé à tout moteur-fusée, quels que soient les ergols consommés par celui-ci. La figure 1 représente un moteur-fusée 100 équipé d'un dispositif 50 5 de détection d'allumage d'une chambre à combustion de moteur-fusée conforme à l'invention. Le dispositif 50 permet de détecter le statut d'allumage de la chambre à combustion du moteur-fusée 100. Le moteur fusée 100 présenté par la figure 1 comporte une chambre 10 de combustion 10, une tuyère 20 présentant un col de tuyère 22. La chambre de combustion 10 est alimentée par un réservoir d'hydrogène 12A et un réservoir d'oxygène 12B, respectivement via des pompes 14A et 14B. Deux capteurs de pression 30A et 30B et deux capteurs de température 32A et 32B sont disposés au voisinage des points d'injection des 15 ergols dans la chambre de combustion. Les capteurs 30A et 30B mesurent les pressions des ergols au point d'injection dans la chambre de combustion 10. Un capteur 34 mesure la presssion Pc dans la chambre de combustion, et un capteur 36 mesure la pression, dite 'pression atmosphérique' PATm à l'extérieur du moteur 100. 20 Le moteur 100 comporte enfin un contrôleur 50, dont la structure et le fonctionnement seront détaillés plus loin. Le contrôleur 50 est relié aux capteurs de pression 30A, 30B et 36 et aux capteurs de température 32A, 32B de manière à acquérir les informations de pression et de température produites par ces capteurs. 25 Le procédé présenté vise à détecter l'allumage du moteur-fusée 100, c'est-à-dire le début de la combustion des ergols dans la chambre de combustion. L'information fournie par le procédé est un statut d'allumage S. Ce statut peut prendre trois valeurs : 'allumé', 'non-allumé', ou encore 30 'indéterminé' (On assimile l'allumage de la chambre de combustion à l'allumage du moteur). Le procédé est lancé pendant la procédure (ou du procédé) de démarrage du moteur-fusée 100. Dès que la chambre de combustion est alimentée avec les deux ergols 35 (hydrogène et oxygène), on lance le procédé de détection d'allumage, à un instant initial tO. On effectue alors de manière itérative les étapes A) à D) décrites précédemment ; à chaque itération, l'indice T de l'algorithme est incrémenté. 12 Le temps t est mesuré à partir de l'instant tO. A chaque boucle, l'instant correspondant à l'itération i est noté ti. A chaque boucle, on réalise donc les étapes suivantes : 5 A) Détermination de la pression (PHI,POI) et de la masse volumique (p OI, 1_-11) des ergols injectés, ' et cg_.essioriextérieure P -rm Pour chacun des ergols, on détermine la pression et la masse volumique des ergols lors de leur injection dans la chambre de combustion. De manière générale, ces grandeurs peuvent être mesurées par des capteurs ou 10 éventuellement estimées ou déterminées de manière indirecte par calcul. En l'occurrence, dans le mode de mise en oeuvre présenté : - on mesure les pressions PHI, POI des ergols à l'aide des capteurs 30A,30B ; - on mesure les températures THI, TOI des ergols à l'aide des capteurs 15 32A,32B ; - on calcule alors les masses volumiques p_0I, p_HI des ergols à partir des pressions et températures mesurées, respectivement PHI, POI et THI, TOI. Les pressions et températures des ergols sont mesurées au niveau des 20 points d'injection ou à faible distance en amont de ceux-ci sur les circuits d'alimentation en ergols de la chambre de combustion 10. Même dans le deuxième cas (mesure effectuée à faible distance en amont des injecteurs, sur les circuits d'alimentation en ergols de la chambre de combustion 10), les points de mesure sont considérés comme étant les points d'injections des 25 ergols dans la chambre de combustion. De plus, on acquiert la valeur de la pression atmosphérique PA-rm à l'extérieur du moteur-fusée par le capteur 36. B) Calcul des valeurs moyennes et des écarts-types de la pression théorique 30 dans la chambre de combustion La pression théorique désigne la pression dans l'une des hypothèses HO (non-allumage) ou Hl (allumage). B1) Valeurs moyennes de la pression théorique Les débits d'oxygène QO et d'hydrogène QH sont calculés à partir du 35 coefficient de pertes de charge k, des masses volumiques des ergols p0 et pH et des pressions POI et PHI au niveau des injecteurs d'oxygène et d'hydrogène, grâce à la relation : (1) Qxqp (PxI-Pcc)/k 13 où PGC est la pression dans la chambre de combustion, et PxI la pression au point d'injection (x désigne l'ergol : 0 pour Oxygène, H pour Hydrogène), p est la masse volumique, et k est le coefficient de pertes de charge 5 correspondant à l'injection dans la chambre, qui dépend de caractéristiques du circuit connues à l'avance. On en déduit le débit total Qtot, égal à la somme des débits d'oxygène et d'hydrogène, et le rapport de mélange RM, égal au rapport entre le débit de comburant et le débit de carburant (QO/QH). 10 La vitesse caractéristique (en m/s) des gaz, notée C*, est le coefficient issu des paramètres des gaz servant pour des calculs des débits en régime sonique. La vitesse caractéristique C* est alors déterminée en fonction du rapport de mélange RM au moyen d'une table. On utilise de plus le coefficient Seff, qui est un paramètre fixé égal à 15 la section du col de la tuyère multipliée par un coefficient d'efficacité représentant le rendement d'écoulement. Ce coefficient est inférieur et voisin de l'unité. La valeur moyenne théorique de la pression dans la chambre de combustion, dans l'hypothèse d'allumage du moteur (H1), est alors calculée 20 grâce à la relation suivante : Qtotc* (2) PGC H1 = Se f f Inversement, dans l'hypothèse de non-allumage du moteur (HO), la valeur moyenne de la pression dans la chambre de combustion est considérée, dans 25 la nécessité d'un calcul simplifié, comme étant égale à la pression de l'environnement, c'est-à-dire la pression du « vide » ou de la pression « atmosphérique » soit : (3) PGCHO = P atm OÙ Patm est la pression à l'extérieur du moteur, mesurée par le capteur de 30 pression 36. B2) Ecarts-types de la pression théorique L'écart-type a PGC_110 de la pression théorique dans la chambre de combustion en cas de non-allumage de la chambre (hypothèse HO) est donné par la formule suivante : (4) 0.PGC_HO = 0-PGCmes 35 où crPGCmes représente l'écart-type sur les mesures de pression acquises en entrée du système ; en effet, en cas de non-allumage, l'incertitude provient 14 uniquement des capteurs. La précision donnée à 3 écart-types pour un capteur est usuellement de l'ordre de 1 % de l'étendue de mesure du capteur. Par exemple, pour un capteur de pression configuré pour permettre la mesure d'une pression 5 maximale de 100 Bars, l'écart-type a associé au capteur peut être choisi égal à un tiers de 1% de l'étendue de mesure, soit 0,33 Bar. L'écart-type aPGC_Hl de la pression théorique dans la chambre de combustion en cas d'allumage de la chambre (hypothèse H1) est calculé de la façon suivante. 10 On utilise les débits (QO, QH) de chacun des ergols calculés à l'étape A). L'écart-type a PGC_Hl de la pression théorique dans la chambre de combustion en cas d'allumage de la chambre est calculé en réalisant les étapes suivantes : 15 . on calcule le rapport de mélange RM du moteur-fusée et son écart-type aRm par les équations suivantes : ,i(r)o (P01 - Patm) RM = .\/()H (PHI - Patm) 1 (00 1 cr2+ (P01- Patm)2 ) 2 OERM - -a (l E (PHI - Patm) ( P01 (PHI - Patm)2 aPHI 'k) H . on calcule alors les écarts-types a_C* de la vitesse caractéristique C*, et d cgo-r du débit total des ergols injectés dans la chambre de combustion, de la manière suivante : 20 L'écart-type a_C* est fonction de l'écart-type du rapport de mélange suivant une fonction déterminée à l'avance : ac* = f(o-Rm) L'écart-type du débit total a Qtot est calculé en fonction des écarts-types des débits respectifs des deux ergols : 25 aQtot = aQH+QO = aQ0 + aQH Pour chaque ergol, l'écart type du débit Q est déterminé en fonction de l'écart-type de la pression P par la formule : 1 p CrQx _ - iTi-k aPxI 30 15 . On calcule ensuite l'écart-type a PGC_Hl de la pression théorique dans la chambre de combustion, en cas d'allumage de la chambre au moyen de la formule : 1 OEPGC H1 = Seff (U)017.0aP°I) 4- -1_1 (-ilicYPHl) iiP\ I \2 ((p1 1 2) 4-U0 (1301 - Patm) +.j(DH (PHI - Patm)) 5 C) Détermination du statut d'allumage du moteur Le statut d'allumage du moteur (ou 'statut d'allumage') est une variable S qui indique que l'on considère, pour l'instant considéré, que le moteur est allumé, qu'il n'est pas allumé, ou que l'on ne sait pas. 10 Cl) Détermination des coefficients (A) et (B) permettant de calculer les limites de pression de non-allumage et de pression d'allumage Conformément à la théorie du test de Wald, on détermine alors des coefficients permettant le calcul des limites de pression de non-allumage (A) et de pression d'allumage (B), à partir de la probabilité de non-détection de 15 l'allumage (Pnd) et de la probabilité de fausse alarme de détection d'allumage (Pf), à l'aide des équations suivantes : A = Pnd / (1-Pf) B = (1 - Pnd)/Pf Dans le contexte de détection pour un élément en relation avec la sécurité, 20 les valeurs de la probabilité de non-détection de l'allumage (Pnd) et la probabilité de fausse alarme de détection d'allumage (PO sont des valeurs que l'on choisit en faisant un compromis entre un taux de fausse alarme (disponibilité) et de sécurité (taux de non détection). 25 C2) Détermination du statut d'allumage du moteur Le statut d'allumage S est calculé en deux étapes. On calcule d'abord les limites de pression Al et Bi, et la somme : (CfPGCH1) A1= 21nA + 2n ln aPGCHO 16 = ( 2 / vn (PGCi - PGCHO ) pGcHo VJGCi - PGCHi )2) 1 L 2 2 a 0-PGCHi i=n-p B1= 21nB + 2n ln (apcc_H1) Dans ce calcul, p est le nombre de valeurs sur lesquelles on se fonde pour déterminer si l'on est dans une situation d'allumage ou de non-allumage ; le 5 paramètre p définit donc la durée de la période ou fenêtre temporelle prise en compte dans le calcul.Pour déterminer le statut d'allumage S, on compare alors la somme z aux valeurs Al et B1 : - Si 1 131 : on fixe la valeur du statut d'allumage S à la valeur « allumé ». A chaque boucle, les termes Al, Bi, E sont recalculés sur la base des valeurs 15 actualisées de la pression et de l'écart-type des pressions théoriques respectivement dans l'hypothèse de non-allumage HO ou dans l'hypothèse d'allumage Hl. A l'issue d'une boucle, si l'on détermine que le moteur est allumé ou non-allumé, on réinitialise la valeur de p à une valeur par défaut. Si inversement, 20 pour l'itération considérée, l'algorithme conclut à une situation d'indétermination, on incrémente la valeur de p pour la boucle suivante, sans dépasser toutefois une valeur maximale. D) Contrôle de l'accumulation d'ergols imbrûlés 25 Sur la base du statut d'allumage du moteur, on calcule alors un indicateur d'allumage X. L'indicateur d'allumage est généralement égal au statut d'allumage et indique par conséquent si l'on considère que le moteur est allumé ou non ; cependant, l'indicateur d'allumage peut prendre en compte une valeur supplémentaire qui indique, dans le cas où l'on considère que le 30 moteur ne s'est pas allumé, qu'il faut interrompre la procédure d'allumage du moteur. En effet, après un certain délai, appelé délai d'allumage et noté T, si l'allumage du moteur n'a pas encore eu lieu, l'accumulation d'ergols imbrûlés CrPGC_HO 17 dans la chambre de combustion rend l'allumage dangereux pour le moteur. Par conséquent, au-delà du délai d'allumage T, il est préférable d'évacuer les ergols accumulés dans la chambre de combustion, puis de relancer un cycle complet d'allumage du moteur. 5 Aussi, lorsque le statut d'allumage S a été déterminé à l'étape C : Si à l'instant ti considéré, le délai maximal T n'est pas atteint (ti T): on affecte à l'indicateur d'allumage X une valeur signalant que l'allumage du moteur doit être interrompu et qu'une 10 procédure de nettoyage du moteur doit être réalisée. Après cette procédure de nettoyage, une nouvelle tentative (éventuelle) d'allumage du moteur pourra éventuellement être effectuée. Les valeurs sur lesquelles repose le test de Wald effectué à l'étape D, à savoir les valeurs moyennes et l'écart-type de la pression théorique suivant 15 les hypothèses HO et H1 (PGC_HO, PGC_Hl, a PCG_HOr a PCG_H1) sont recalculées à chaque itération de l'algorithme. Ainsi avantageusement, le test réalisé est actualisé à chaque itération en fonction des caractéristiques du moteur à l'instant considéré. Après cette présentation du fonctionnement du contrôleur 50 et notamment 20 de l'algorithme mis en oeuvre dans le calculateur 50, la structure matérielle de ce contrôleur va maintenant être présentée. Le contrôleur 50 comporte : A) une unité d'acquisition 52, configurée pour acquérir, pour chacun des ergols (H,0), la pression d'injection et la masse volumique de l'ergol aux points d'injection dans la chambre de combustion ; 25 B) une unité de calcul de valeurs moyennes et d'écarts-types 54 configurée pour, en fonction des pressions et masses volumiques calculées à l'étape A), calculer des valeurs moyennes PGC_HO, PGC_Hl et des écarts-types a_PCGHO, a_PCGH1 d'une pression théorique d'allumage et d'une pression théorique de non-allumage dans la chambre de combustion, qui sont les 30 pressions dans la chambre de combustion respectivement dans l'hypothèse de non-allumage (HO) qui correspond au cas où la combustion n'a pas lieu et dans l'hypothèse d'allumage (H1) qui correspond au cas où la combustion a lieu ; C) une unité de détermination de statut d'allumage du moteur 56 configurée 35 pour déterminer un statut d'allumage du moteur (S). L'unité d'acquisition 52 calcule la masse volumique de chacun des ergols à partir de paramètres acquis par des capteurs correspondants, à savoir en fonction de la température et de la pression des ergols, sur la base 18 d'abaques établies à l'avance pour les ergols contenus dans les réservoirs. L'unité de détermination de statut d'allumage du moteur (S) est configurée pour réaliser les opérations suivantes, décrites précédemment : Cl) Détermination des coefficients (A) et (B) permettant de calculer les 5 limites de pression de non-allumage et de pression d'allumage, et C2) attribuer au statut d'allumage du moteur (S) une valeur allumé, non allumé ou indéterminé, indiquant respectivement que l'on considère, pour l'instant considéré, que le moteur est allumé, qu'il n'est pas allumé, ou que l'on ne sait pas. 10 L'unité de détermination de statut d'allumage du moteur calcule le statut d'allumage S en suivant les mêmes étapes Cl et C2 que dans le procédé selon l'invention. Dans ce mode de réalisation, les modules fonctionnels 52 à 56 décrits précédemment sont des modules fonctionnels du dispositif 50 de détection 15 d'allumage d'une chambre à combustion de moteur-fusée. Les unités d'acquisition 52, de calcul de valeurs moyennes et d'écarts-types 54 et de détermination de statut d'allumage du moteur 56 sont réalisées sous forme de modules fonctionnels d'un ordinateur. Cet ordinateur qui constitue le dispositif de contrôle géométrique porte la référence 50 sur 20 les Fig.1 et 2. Le dispositif de contrôle géométrique 50 dispose de l'architecture matérielle d'un ordinateur, telle qu'illustrée schématiquement à la figure 2. Il comprend notamment un processeur 4, une mémoire vive 5, une mémoire morte 6, une mémoire flash non volatile 7, ainsi que des moyens de 25 communication 8. La mémoire morte 6 du dispositif de contrôle géométrique 50 constitue un support d'enregistrement conforme à l'invention, lisible par le processeur 4 et sur lequel est enregistré un programme d'ordinateur conforme à l'invention, comportant des instructions pour l'exécution des 30 étapes du procédé de détection d'allumage d'une chambre à combustion de moteur-fusée conforme à l'invention REVENDICATIONS1. Procédé de détection d'allumage d'une chambre à combustion (10) de moteur-fusée, le procédé comportant les étapes suivantes, réalisées de manière itérative au moyen d'un calculateur (50) : A) pour chacun des ergols (H,0), on acquiert des valeurs de grandeurs qui permettent de calculer des valeurs moyennes PGC_HO, PGC_Hl et des écarts-types a_PGCHO, a_PGCH1 d'une pression théorique d'allumage et d'une pression théorique de non-allumage dans la chambre de combustion, qui sont les pressions dans la chambre de combustion respectivement dans l'hypothèse de non-allumage (HO) qui correspond au cas où la combustion n'a pas lieu et dans l'hypothèse d'allumage (H1) qui correspond au cas où la combustion a lieu ; B) en fonction des valeurs acquises à l'étape A), on calcule les valeurs 15 moyennes PGC_HO, PGC_H1 et les écarts-types a_PGCHO, a_PGCH1 de la pression théorique d'allumage et de la pression théorique de non-allumage dans la chambre de combustion ; C) on détermine un statut d'allumage du moteur (S) ; l'étape C) comprenant les étapes : 20 C1) on détermine des coefficients A et B permettant de calculer des pressions limites d'allumage et de non-allumage, à partir d'une probabilité de non-détection de l'allumage Pnd et d'une probabilité de fausse alarme de détection d'allumage (Pf), à l'aide des équations suivantes : A = Pnd / (1-Pf) 25 B = (1 - Pnd)/Pf C2) on attribue à un statut d'allumage du moteur (S) une valeur allumé, non allumé ou indéterminé, indiquant respectivement que l'on considère, pour un instant considéré, que le moteur est allumé, qu'il n'est pas allumé, ou que l'on ne sait pas ; le statut d'allumage (S) étant obtenu : 30 . en calculant des bornes d'allumage/non-allumage Al et B1 et une somme z à l'aide des équations suivantes : Al = 21nA + 2n ln (crpccHi) 11 ((PGCi - PGCHO )2 (PGCi -pGP cGHC:11 )2) i=n-p aP2GCHo 0-2 ePGCHo E = 20 B1= 21nB + 2n ln (upGc_H1) où : - PGCi est la pression mesurée dans la chambre de combustion à 5 l'instant (ti) correspondant à une boucle de calcul en cours ; - n est l'indice de la boucle de calcul à l'instant considéré ; - p est un entier prédéterminé ; et . en affectant au statut d'allumage (S) une valeur parmi des possibilités Non-allumé, Allumé, Indéterminé, en comparant la somme E aux bornes 10 d'allumage/non-allumage Al et B1 : - Si E 131 : on fixe la valeur du statut d'allumage (S) à la valeur « allumé ». 2. Procédé de détection d'allumage selon la revendication 1, dans lequel à l'étape B), la valeur moyenne PGC_Hl de la pression théorique d'allumage dans la chambre de combustion dans l'hypothèse où la combustion a lieu est calculée à partir du débit total des ergols injectés dans la chambre de combustion (Qtot), de la vitesse caractéristique des gaz de combustion (C*) et de la section efficace (Seff) d'un divergent d'une tuyère (20) du moteur-fusée. 3. Procédé de détection d'allumage selon la revendication 2, dans lequel à l'étape B), la valeur moyenne PGC_Hl de la pression théorique d'allumage dans la chambre de combustion dans l'hypothèse où la combustion a lieu est calculée grâce à la relation suivante : Qsteo ft cf * (2) PGC_Hl = 4. Procédé de détection d'allumage selon l'une quelconque des revendications 1 à 3, dans lequel lors de l'étape A, on acquiert de plus une pression dans la chambre de 35 combustion ; et aPGC_HO 21 l'écart-type a PGC_HO de la pression théorique dans la chambre de combustion dans l'hypothèse (HO) de non-allumage de la chambre est donné par la formule suivante : GPGC_HO = GPGCmes 5 01.1 OPGCmes représente l'écart-type de la mesure de pression dans la chambre de combustion. 5. Procédé de détection d'allumage selon l'une quelconque des revendications 1 à 4, dans lequel 10 lors de l'étape A, on acquiert de plus une valeur de la pression atmosphérique PATM à l'extérieur du moteur-fusée ; à l'étape B, en fonction des pressions et masses volumiques calculées à l'étape A), on calcule des débits (Q0, QH) de chacun des ergols ; et l'écart-type a PGC_Hl de la pression théorique dans la chambre de 15 combustion en cas d'allumage de la chambre est calculé en réalisant les étapes suivantes : . on calcule un rapport de mélange RM du moteur-fusée et son écart-type aRM par les équations suivantes : .\/()0 (P01 - Patm) RM = \le (p.._ Patm) 1 Mo 1 (P01 - Patm)2 2 ) a2 j_ aRM = \/2 (i') (PHI - Patm) ( P01 ' (PHI - Patm)2 ' UUH . on calcule alors des écarts-types cr_C* de la vitesse caractéristique C* et a 20 QT0T du débit total des ergols injectés dans la chambre de combustion, par les équations suivantes : acs = RaRm) aQtot = aQH+QO ="QO ' "QH avec pour chaque ergol aQx - Pxl 25 dans lesquelles : x représente un ergol quelconque, apoi et 6pHi désignent les écarts-types des pressions POA et PHI des deux ergols à leurs points d'injection, 22 kx le coefficient de perte de charge pour l'ergol x, o-Qx et ap,a , pour l'ergol x, l'écart type respectivement d'un débit et d'une pression de l'ergol x au point d'injection, le rapport de mélange RM des deux ergols est égal au rapport entre le débit 5 de comburant sur le débit de carburant à l'injection dans la chambre de combustion, QO, QH et Qtot sont respectivement un débit de comburant, de carburant, et un débit cumulé des deux ergols, dans la chambre de combustion, C* est une vitesse caractéristique des gaz de combustion en sortie de la 10 chambre de combustion, et f est une fonction prédéterminée ; . on calcule ensuite l'écart-type a PGC_H1 de la pression théorique dans la chambre de combustion, en cas d'allumage de la chambre au moyen de la formule : 2 2 \ 1 ( f (P) i a ) ) 1 a ) + G: n PHI n H ,-(1-1 2 2 eC* - aPGC H1 = - Seff \ 15 + ( j(Z) (P01 - Patm) + (PHI - Patm)) o k H 6. Procédé de détection d'allumage selon l'une quelconque des revendications 1 à 5, dans lequel les grandeurs dont les valeurs sont acquises à l'étape A) incluent, pour chaque ergol, une pression (POI,PHI) et une masse volumique (p_OI, p_HI) de l'ergol à un point d'injection dans la 20 chambre de combustion. 7. Procédé de démarrage d'un moteur-fusée, dans lequel lors du démarrage d'au moins un moteur du moteur-fusée, on détecte si l'allumage du moteur a bien eu lieu en exécutant le procédé de détection d'allumage selon l'une 25 quelconque des revendications 1 à 6. 8. Programme d'ordinateur sur un support d'informations, ce programme étant susceptible d'être mis en oeuvre par ordinateur, et comportant des instructions de code de programme pour l'exécution des étapes d'un procédé 30 selon l'une quelconque des revendications 1 à 7, lorsque le programme est exécuté sur un ordinateur. 9. Support d'enregistrement lisible par un ordinateur sur lequel est enregistré 23 un programme d'ordinateur selon la revendication 8. 10. Dispositif de détection d'allumage d'une chambre à combustion (10) de moteur-fusée, le dispositif comportant : 5 A) une unité d'acquisition (52), configurée pour acquérir, pour chacun des ergols (H,0), des valeurs de grandeurs qui permettent de calculer des valeurs moyennes PGC_HO, PGC_Hl et des écarts-types a_PGCHO, cr_PGCH1 d'une pression théorique d'allumage et d'une pression théorique de non-allumage dans la chambre de combustion, qui sont les pressions dans la 10 chambre de combustion respectivement dans l'hypothèse de non-allumage (HO) qui correspond au cas où la combustion n'a pas lieu et dans l'hypothèse d'allumage (H1) qui correspond au cas où la combustion a lieu ; B) une unité de calcul de valeurs moyennes et d'écarts-types (54) configurée pour, en fonction des valeurs acquises par l'unité d'acquisition (52), calculer 15 les valeurs moyennes PGC_HO, PGC_Hl et les écarts-types a_PCGHO, a_PCGH1 de la pression théorique d'allumage et de la pression théorique de non-allumage dans la chambre de combustion ; C) une unité de détermination de statut d'allumage du moteur (56) configurée pour déterminer un statut d'allumage du moteur (S) ; l'unité de 20 détermination de statut d'allumage du moteur (S) étant configurée pour : Cl) déterminer des coefficients A et B permettant de calculer des pressions limites d'allumage et de non-allumage, à partir d'une probabilité de non-détection de l'allumage Pnd et d'une probabilité de fausse alarme de détection d'allumage Pf, à l'aide des équations suivantes : 25 Pnd, 1-Pnd A= - et B= = ; et 1-Pf Pf C2) attribuer à un statut d'allumage du moteur (S) une valeur allumé, non allumé ou indéterminé, indiquant respectivement que l'on considère, pour un instant considéré, que le moteur est allumé, qu'il n'est pas allumé, ou que 30 l'on ne sait pas ; le statut d'allumage (S) étant obtenu : . en calculant des bornes d'allumage/non-allumage Al et Bi et une somme E à l'aide des équations suivantes : (CrpGCH1) A1= 21nA + 2n ln CFPG CHO ) z = In ((pGct _pGP cGH0CH0 )2 (PGCt - PGCH1 )2 CI 2 2 (3-PGCHi i=n-p\ 24 B1= 21nB + 2n ln (crpcc_in) où : - PGCi est la pression mesurée dans la chambre de combustion à l'instant (ti) correspondant à la boucle de calcul en cours ; 5 - n est l'indice de la boucle de calcul à l'instant considéré ; - p est un entier prédéterminé ; . en affectant au statut d'allumage (S) une valeur parmi des possibilités Non-allumé, Allumé, Indéterminé, en comparant la somme / aux bornes d'allumage/non-allumage Al et Bi, et : 10 - Si / B1 : en fixant la valeur du statut d'allumage (S) à la valeur 15 « allumé ». 11. Dispositif de détection d'allumage selon la revendication 10, dans lequel dans lequel les grandeurs dont l'unité d'acquisition (52) est configurée pour acquérir les valeurs incluent, pour chaque ergol, une pression (POI,PHI) et 20 une masse volumique (p_OI, p_HI) de l'ergol à un point d'injection dans la chambre de combustion. 12. Moteur-fusée (100), comprenant un dispositif de détection d'allumage d'une chambre à combustion (10) de moteur-fusée selon la revendication 10 25 ou 11. CFPGC_HO