← Acoustics & propulsion
Aérodynamique

Aérodynamique des petites hélices de drone aux faibles nombres de Reynolds

Revue de la littérature ouverte : Reynolds local, transition, bulles de séparation laminaire, et limites de la BEMT et de la CFD sur les petites hélices.

Original written in French — English translation coming.

Revue scientifique sur l’aérodynamique des petites hélices électriques de drones aux faibles nombres de Reynolds

Cadre physique, notations et ordres de grandeur

La littérature ouverte sur les petites hélices électriques de drones couvre surtout des diamètres d’environ 5 à 21 pouces dans les campagnes UIUC, des propulseurs encore plus petits de type micro-air-vehicle, et quelques hélices optimisées pour drones stratosphériques. Pour ces machines, le régime aérodynamique n’est pas « simplement incompressible et quasi-stationnaire » : il est dominé par des nombres de Reynolds chordaux locaux souvent compris entre 10410^4 et 10510^5, avec une forte variabilité radiale, des transitions laminaire-turbulent sensibles aux perturbations, des bulles de séparation laminaires, et des écarts parfois importants entre mesures, BEMT et CFD lorsqu’on approche des faibles charges, du haut rapport d’avance, du post-décrochage ou du très bas Reynolds. Merchant résume l’intervalle global des petites hélices testées à Wichita comme Re3×104Re \approx 3\times 10^4 à 3×1053\times 10^5 pour des diamètres de 6 à 22 pouces, tandis que Deters met en évidence, pour des micro-hélices et petites hélices, des régimes encore plus bas, jusqu’à Re1.6×104Re \approx 1.6\times 10^4 sur une 5×3.5 à 4000 rpm. Plus récemment, Grava et al. ont étudié une petite hélice de rayon 7.5 cm en vol axial/cross-flow avec un RecRe_c local variant d’environ 5000 à 15000, et Grande et al. ont documenté, pour un propulseur acoustique de référence, un ReRe de pointe allant de 4.3×1044.3\times 10^4 à 10510^5.

Les coefficients usuels sont

J=VnD,CT=Tρn2D4,CP=Pρn3D5,CQ=Qρn2D5,J=\frac{V_\infty}{nD},\qquad C_T=\frac{T}{\rho n^2 D^4},\qquad C_P=\frac{P}{\rho n^3 D^5},\qquad C_Q=\frac{Q}{\rho n^2 D^5},

avec P=2πnQP=2\pi nQ, d’où CP=2πCQC_P=2\pi C_Q, et

η=TVP=JCTCP.\eta=\frac{TV_\infty}{P}=J\,\frac{C_T}{C_P}.

UIUC emploie ces définitions standards, et précise que, pour les courbes statiques CT0C_{T0} et CP0C_{P0}, le Reynolds reporté est défini à la station r/R=0.75r/R=0.75 à partir de c75%c_{75\%} et de V75%=0.75ΩRV_{75\%}=0.75\,\Omega R. Les mêmes définitions sont reprises dans la littérature récente sur la conception BEMT d’hélices UAV.

Pour fixer les ordres de grandeur radiaux, l’APC Slow Flyer 10×7 de UIUC est un excellent cas témoin parce que sa géométrie radiale est publiée. À r/R=0.75r/R=0.75, on lit c/R=0.197c/R=0.197, soit c7525.0c_{75}\approx 25.0 mm pour D=0.254D=0.254 m; le vrillage géométrique tombe d’environ 34.934.9^\circ à r/R=0.15r/R=0.15 à 14.414.4^\circ à r/R=0.75r/R=0.75, puis 8.48.4^\circ au tip. À 3008 rpm, le fichier UIUC donne CTC_T, CPC_P et η\eta sur toute la plage de JJ; en combinant cette géométrie et la définition UIUC du Reynolds, on obtient, pour ce propulseur, un RecRe_c local qui varie typiquement d’environ 5.6×1035.6\times 10^3 près de r/R=0.15r/R=0.15 à 5.1×1045.1\times 10^4 autour de r/R=0.75r/R=0.75 en statique, puis jusqu’à 5.4×104\approx 5.4\times 10^4 à J0.8J\approx 0.8. Le Mach local reste modéré à 3008 rpm, de l’ordre de 0.02 à la racine et 0.12\approx 0.12 au tip; en revanche, pour la même hélice à 6006 rpm, le Mach de pointe monte vers 0.23 et le Re75Re_{75} dépasse 10510^5. Ces valeurs sont des calculs à partir des géométries et vitesses UIUC; elles illustrent bien que le « Reynolds de l’hélice » n’existe pas comme grandeur unique utile au sens didactique.

L’impact de la transition et du Mach local est également visible dans la pratique industrielle. APC indique que ses hélices « thin electric », « slow-fly » et multirotor combinent typiquement un profil de type Eppler E63 en emplanture et un profil de type Clark-Y vers le tip; APC précise aussi que, lorsque les effets de Mach deviennent significatifs près du tip, la société recourt à du washout de pas et/ou à une réduction de cambrure pour limiter la remontée de traînée compressible. Autrement dit, même si beaucoup de petites hélices de drones opèrent à Mtip0.25M_{tip}\lesssim 0.25, la compressibilité n’est pas totalement négligeable dès que le diamètre et la vitesse périphérique augmentent.

Le point central pour la physique des faibles Reynolds est le suivant : aux ReRe de l’ordre de 10410^410510^5, la couche limite laminaire est très propice à la séparation, à la transition dans la couche cisaillée séparée puis au recollement, ce qui forme une bulle de séparation laminaire. Cette dynamique modifie fortement les polaires sectionnelles, donc la distribution radiale de Cl/CdC_l/C_d, donc CTC_T, CPC_P, le rendement et même le bruit. La revue récente de Giacomini et al. sur les modèles de transition rappelle précisément que les phénomènes dominants à Re<5×104Re<5\times 10^4 sont la séparation précoce, les LSB et la transition retardée, et Grande et al. montrent expérimentalement que ces LSB peuvent même devenir une source acoustique identifiable sur une hélice réelle.

Résultats expérimentaux

Synthèse quantitative des bancs les plus utiles

Le socle expérimental le plus robuste reste aujourd’hui la base UIUC, qui regroupe quatre volumes de mesures en soufflerie sur des hélices de petits UAV et d’aéromodèles, avec définitions cohérentes, corrections de montage et données brutes téléchargeables. Le volume 1 couvre près de 140 hélices principalement entre 9 et 11 pouces; le volume 2 ajoute les micro-propulseurs et une analyse explicite des effets de Reynolds; les volumes 3 et 4 étendent le jeu à des hélices repliables puis à 17 APC Thin Electric fixes de 12 à 21 pouces. UIUC montre en outre une bonne concordance avec des mesures Ohio State sur APC 10×5E et 10×7E malgré des chaînes d’effort différentes.

Étude expérimentaleGéométrie et profilRégime et JJReRe / MMRésultat quantitatifIncertitude
UIUC Vol. 1, APC SF 10×72 pales, D=10D=10” (0.2540.254 m), P/D=0.7P/D=0.7, c/R0.75=0.197c/R\vert_{0.75}=0.197, vrillage notable du root au tipStatique, 2283–5987 rpmRe75Re_{75} calculé UIUC-style 3.9×104\approx 3.9\times 10^4 à 1.0×1051.0\times 10^5; MtipM_{tip} calculé 0.09\approx 0.09 à 0.23CTC_T statique croît de 0.1409 à 0.1606; CPC_P de 0.0678 à 0.0797Non reportée dans le fichier
UIUC Vol. 1, APC SF 10×7 à 3008 rpmMême géométrieJ=0.192J=0.192 à 0.911RecRe_c radial calculé 5.6×103\approx 5.6\times10^35.1×1045.1\times10^4 en statique; 8.9×103\approx 8.9\times10^35.2×1045.2\times10^4 à J=0.573J=0.573; Mtip0.12M_{tip}\approx 0.12CT=0.1257C_T=0.1257, CP=0.0681C_P=0.0681, η=0.355\eta=0.355 à J=0.192J=0.192; η\eta culmine à 0.708\approx 0.708 vers J=0.573J=0.573 avec CT=0.0607C_T=0.0607, CP=0.0492C_P=0.0492; CTC_T devient négatif à J=0.862J=0.862Non reportée dans le fichier
Deters 2014, Micro Invent 5×3.55” micro-propellerAvancement variable; 4000–8000 rpmReRe reporté 16000, 23900, 31900Les courbes CTC_T, CPC_P, η\eta montrent une augmentation nette de performance avec ReRe croissantNon précisée dans cet extrait
Deters 2014, DA4002 9×6.759”, 2 palesAvancement variable; 2000–5000 rpmReRe reporté 24100, 36200, 48400, 61800CTC_T augmente avec ReRe, et le rendement maximal se décale vers le haut quand ReRe monteNon précisée dans cet extrait
Deters 2014, comparaison géométriquement semblable 5” vs 9”5” et 9” NR640 / géométries semblablesMême ReRe approximatif5” à 6000 rpm : Re=22100Re=22100; 9” à 2000 rpm : Re=24100Re=24100; 10 ft pleine échelle : Re=1.8×106Re=1.8\times10^6À ReRe semblable, les courbes 5” et 9” sont très proches; mais la 9” et la 10 ft diffèrent fortement, avec une efficacité max 9” < 65% contre > 80% en pleine échelleNon précisée
Grande et al. 2022Propulseur 2 pales de benchmark acoustique, pales propres vs transition forcéeJ=0J=0 à 0.6Retip=4.3×104Re_{tip}=4.3\times 10^4 à 10510^5À J=0.4J=0.4 et 0.6, la transition forcée supprime le hump de bruit LSB et réduit le bruit d’environ 5 dB puis 10 dBNon reportée dans l’abstract
Grava et al. 2025Hélice 2 pales de rayon 7.5 cm, dérivée APC-96, profils NACA 4412 de r/R=0.2r/R=0.2 au tip, imprimée 3Dn=3000n=3000 rpm, cross-flow V=0V_\infty=0–15 m/s, μ=0\mu=0–0.64RecRe_c local 5000–15000; MtipM_{tip} calculé 0.07\approx 0.07La poussée et la puissance augmentent sensiblement avec la vitesse de cross-flow; les lois semi-empiriques rotor/hélicoptère deviennent applicables après recalageErreur PIV moyenne < 0.1%; stats d’ordre 2 < 5%
Comunian et al. 2026Deux petites hélices côte à côte, géométrie fixe, entre-axes minimisésJ=0J=0, 0.3, 0.45 sur SPIV; enveloppe complète mesuréeReRe des pales souvent < 10510^5Interaction latérale : baisse moyenne de rendement propulsif de 3.2%, effet plus fort aux faibles JJNon reportée dans l’abstract

Ce que montrent les mesures, au-delà des chiffres

Le résultat expérimental le plus robuste sur l’effet du Reynolds est celui de Deters : quand une petite hélice de 9 pouces est testée à des ReRe croissants, le CTC_T augmente systématiquement; et quand des géométries de 5 pouces et 9 pouces sont testées à ReRe voisins, les courbes sont très proches. En revanche, à géométrie semblable mais en passant de Re104Re\sim 10^410510^5 vers Re106Re\sim 10^6, les différences de performance deviennent majeures : sur le cas NR640, le rendement max passe de moins de 65% à un peu plus de 80% pour la pleine échelle. Cela montre que le « handicap » des petites hélices n’est pas d’abord leur diamètre nominal, mais le ReRe local auquel travaillent leurs sections.

Le cas UIUC APC 10×7 illustre bien la variation conjointe avec rpm et JJ. En statique, CTC_T et CPC_P augmentent avec rpm sur toute la plage mesurée; à rpm fixé, quand JJ augmente, CTC_T décroît presque monotoniquement, CPC_P décroît plus lentement, et η\eta présente un maximum net autour d’un JJ intermédiaire, ici voisin de 0.57 à 3008 rpm. Le passage vers CT<0C_T<0 aux très grands JJ rappelle qu’un même propulseur peut quitter son régime propulsif et entrer dans une zone de freinage/autorotation partielle si on le compare hors enveloppe.

Sur la transition, les expériences de Grande et al. sont particulièrement importantes parce qu’elles relient directement un phénomène de couche limite à une signature hélice mesurée. Entre J=0.4J=0.4 et 0.6, un hump de bruit haute fréquence est attribué au détachement tourbillonnaire de la bulle de séparation laminaire; le forçage précoce de transition supprime cette bulle et réduit le bruit de 5 à 10 dB. La leçon est double : d’une part, la LSB existe bel et bien sur des hélices réelles aux faibles Reynolds; d’autre part, « tripper » une pale n’est pas seulement une question de bruit, c’est aussi une modification de la physique de couche limite qui peut changer les coefficients aérodynamiques.

Les autres bancs ouverts confirment le message, mais avec leurs spécificités métrologiques. Le banc UBI/MAAT obtient une répétabilité jugée excellente sur trois jours différents et une incertitude typiquement inférieure à 0.3% sur CTC_T, 0.6% sur CPC_P et 1% sur η\eta, hors la zone 5–6 m/s où les fluctuations augmentent. Ce banc retrouve aussi, par comparaison, que les hélices de pas grossier ont une efficacité maximale supérieure aux hélices « fine pitch » dans son échantillon. Merchant, à Wichita, rapporte pour son système intégré une variation run-to-run généralement inférieure à 0.1%, des incertitudes typiques de l’ordre de 0.05% sur CTC_T et CQC_Q, 0.003% sur CPC_P, et un « worst case » de 0.2% sur CTC_T/CPC_P; il note aussi que la dispersion augmente lorsque le Reynolds diminue, en particulier sous 3000 rpm et à basse vitesse de soufflerie. Ces chiffres sont très utiles pédagogiquement parce qu’ils montrent que, si la physique est subtile, la métrologie peut être excellente à condition de caractériser correctement le montage.

Enfin, sur les configurations réellement « drones », les campagnes de Grava et Comunian rappellent que le vol en translation et l’interaction multi-hélices changent la nature même du problème. Pour une petite hélice isolée en cross-flow, la poussée et la puissance augmentent avec VV_\infty, et un modèle de type rotor/hélicoptère recalé peut déjà capturer les tendances. Pour deux hélices côte à côte, la pénalité moyenne de rendement est de 3.2%, plus marquée aux faibles JJ, avec un sillage plus large et une turbulence stochastique plus faible dans certaines zones. Cela signifie qu’un cours centré uniquement sur l’hélice isolée en axe ne couvre qu’une partie du problème drone réel.

Résultats numériques

Ce que le CFD prédit correctement

Le CFD apporte une vraie valeur lorsque le régime est dominé par la transition, la séparation et des géométries non triviales, mais seulement si le modèle est compatible avec le faible Reynolds. Yao et al. ont précisément construit une optimisation itérative d’hélice à faible Reynolds avec CFD transitionnel γReθ\gamma-Re_\theta, en optimisant simultanément les profils de sections et la géométrie globale. Au point de conception stratosphérique (D=2D=2 m, V=32V_\infty=32 m/s, ρ=0.0889\rho=0.0889 kg/m3^3, poussée requise 20 N), l’hélice avec Gurney flap atteint η=82.0%\eta=82.0\% à Reavg4.0×104Re_{avg}\approx 4.0\times10^4, contre 80.2% pour la version sans flap à Reavg4.2×104Re_{avg}\approx 4.2\times10^4. Le gain de rendement est donc 1.8 point au design point, et la puissance consommée est réduite d’environ 2.2% à vitesse d’avance égale. En revanche, à plus haut Reynolds (1.8×105\sim 1.8\times 10^5), cet avantage disparaît pratiquement, voire s’inverse légèrement. C’est l’une des démonstrations les plus nettes du fait qu’une « bonne » modification à Re4×104Re\sim 4\times10^4 ne reste pas forcément bonne à Re2×105Re\sim 2\times10^5.

Yao et al. fournissent aussi une leçon très utile sur la transition et la LSB. Leur CFD et leurs essais montrent que le propulseur muni d’un flap a encore environ 1% d’avantage de rendement vers J1J\approx 1 à Re4×104Re\approx 4\times10^4, alors que cet avantage n’existe plus vers Re1.6×105Re\approx 1.6\times10^5. Les auteurs relient explicitement ce comportement au fait que la LSB est importante à Re4×104Re\approx 4\times10^4 et devient beaucoup moins déterminante quand ReRe dépasse environ 1.4×1051.4\times10^5. Cela rejoint bien l’intuition issue des mesures de bruit LSB : à bas Reynolds, toute géométrie ou rugosité qui modifie la transition peut changer fortement l’aérodynamique; à Reynolds plus élevés, la même modification devient souvent surtout une source de traînée additionnelle.

Le papier de Kutty et Rajendran sur l’APC Slow Flyer 10×7 est très instructif pour les limites du CFD RANS « standard ». Avec FLUENT, un maillage tétraédrique non structuré, un modèle kωk-\omega standard et un MRF à 3008 rpm, le cas est raisonnablement prédit aux faibles rapports d’avance. Mais les écarts augmentent fortement près du grand JJ et du voisinage de la poussée nulle : sur leur comparaison à la base UIUC, l’erreur sur CTC_T atteint 46.27% à J=0.799J=0.799, l’erreur sur CPC_P 10.75% à J=0.192J=0.192, et l’efficacité est surestimée de 42.5% à J=0.799J=0.799. Leur propre conclusion est donc nuancée : cette chaîne CFD est fiable pour les faibles JJ et certains régimes de charge, mais pas pour les zones de haut JJ, de transition délicate ni de CTC_T proche de zéro.

Un autre résultat numérique important concerne la rugosité. Cooke et al., par DNS complétée par des essais sur banc de poussée, montrent qu’à très bas Rec=6500Re_c=6500 des micro-propellers rugueux peuvent donner environ +2% de poussée et -5% de puissance par rapport à une hélice lisse à Reynolds constant. C’est un résultat intéressant, mais il ne faut surtout pas le généraliser aux hélices de drones 5–12 pouces : ce cas est situé dans un régime beaucoup plus visqueux que celui de la majorité des hélices multirotors usuelles, et la rugosité y agit plus comme un déclencheur/organisateur de transition que comme un simple défaut de surface. Autrement dit, la rugosité n’a pas un « signe » universel.

Sur l’aéroélasticité, Möhren et al. montrent en simulation couplée poutre–BEM que l’effet de l’élasticité sur la poussée devient notable surtout pour les grandes hélices, les faibles charges surfaciques et les pales balayées. Dans leur espace paramétrique, l’augmentation du JJ en croisière accroît en valeur relative la perte de poussée due à la déformation; ils soulignent aussi que les essais statiques au sol sont mal adaptés pour quantifier l’impact aéroélastique en croisière. Ils recommandent, pour la conception préliminaire, de réserver les simulations couplées aux grands diamètres, aux planforms balayés, aux faibles disc loadings et aux matériaux souples. Cette conclusion est très probablement transposable aux hélices de drones à pales composites souples, mais beaucoup moins aux petites pales injectées rigides de multicoptères classiques.

Ce que le CFD rate encore souvent

La littérature récente sur la transition à faible Reynolds est assez claire : les modèles RANS sans transition explicite ratent fréquemment l’emplacement de transition, l’extension des LSB et donc les niveaux absolus de portance/traînée. Giacomini et al. rappellent explicitement que les LSB et la séparation précoce restent difficiles à prédire numériquement, surtout sous Re5×104Re\approx 5\times10^4. Les travaux de Grava sont cohérents avec cela : pour sa petite hélice de rayon 7.5 cm, les polaires XFOIL de NACA 4412 doivent être générées à Rec=5000Re_c=5000–15000 avec un modèle eNe^N, précisément parce que les polaires « standard » ne sont pas fiables d’emblée dans cette zone.

Une conclusion pédagogique forte ressort aussi du projet Chalmers 2024 : lorsque leur outil BEMT OPTOPROP utilisait des données NACA 16 issues de Reynolds 7×105\ge 7\times10^5 pour une hélice UAV travaillant en réalité autour de Remax105Re_{max}\approx 10^5, le décalage avec la CFD devenait important; l’intégration de polaires XFOIL adaptées au faible Reynolds améliorait la concordance BEMT–CFD jusqu’à 10 points de rendement. Cette étude n’est pas une validation universelle d’XFOIL, mais elle confirme qu’au faible Reynolds l’erreur dominante n’est souvent pas la formule BEM elle-même : c’est la mauvaise physique sectionnelle injectée dans le modèle.

Modèles analytiques

Équations essentielles et hypothèses

Pour l’hélice isolée en axe, les définitions sans dimension de JJ, CTC_T, CPC_P, CQC_Q et η\eta ont déjà été données plus haut. Elles supposent une vitesse de rotation nn en tr/s, un diamètre DD, une densité ρ\rho, une poussée TT, un couple QQ et une puissance PP. Elles sont robustes et communes à UIUC, aux bancs UAV modernes et aux codes de conception BEMT.

Pour le vol en cross-flow ou assimilé « forward flight », Grava rappelle la forme unidimensionnelle du disque actuateur :

T=2ρAvivi2+2Vvicosαdisk+V2.T=2\rho A v_i\sqrt{v_i^2+2V_\infty v_i\cos\alpha_{disk}+V_\infty^2}.

Dans le cas quasi transverse de son étude, on néglige cosαdisk\cos\alpha_{disk} et on obtient la forme adimensionnée

CT=2λiλi2+μ2,C_T=2\lambda_i\sqrt{\lambda_i^2+\mu^2},

λi=vi/(ΩR)\lambda_i=v_i/(\Omega R) est la vitesse induite réduite et μ=V/(ΩR)\mu=V_\infty/(\Omega R) est le rapport d’avance rotor/hélico. À grand μ\mu, on retrouve l’approximation de Glauert

CT2μλi.C_T \approx 2\mu\lambda_i.

Toujours dans ce cadre, la puissance est décomposée en puissance induite et puissance de profil :

CPi=k1CTλi,CP0=σcDM8(1+k2μ2),C_{P_i}=k_1 C_T\lambda_i,\qquad C_{P_0}=\frac{\sigma c_{D_M}}{8}\left(1+k_2\mu^2\right),

avec σ\sigma la solidité, cDMc_{D_M} un coefficient de traînée moyen de profil, et k1k_1, k2k_2 des coefficients correctifs empiriques. Ces équations sont utiles parce qu’elles distinguent proprement les pertes induites des pertes de profil; elles deviennent surtout utiles quand on cherche à expliquer pourquoi CTC_T et CPC_P ne varient pas de la même manière avec JJ et ReRe.

Validité et limites du disque actuateur

La théorie du disque actuateur reste excellente pour les bilans d’ordre de grandeur et pour relier poussée, vitesse induite et puissance minimale, mais elle repose sur des hypothèses fortes : chargement uniformisé sur le disque, absence de description de la géométrie des pales, et traitement purement moyen de l’écoulement. Elle ne peut donc pas prédire, par elle-même, ni la transition, ni les LSB, ni le stall local, ni les pertes de bout/racine, ni les effets de profil, ni les variations azimutales de charge propres au vol de drone en translation. C’est précisément pour cela que Grava doit « tuner » les lois semi-empiriques issues du rotorcraft avant de les appliquer à sa petite hélice, et pourquoi les études modernes complètent systématiquement le disque par BEMT, des polaires 2D et des corrections.

Validité et limites de la BEMT

La BEMT est le modèle analytique central de la littérature drone. Sa structure suppose qu’on découpe la pale en éléments de rayon rr, que chaque élément se comporte comme un profil 2D soumis à une vitesse résultante locale, puis que l’on boucle avec la théorie de quantité de mouvement annulaire pour retrouver les vitesses induites. Cette construction est remarquablement efficace à coût faible, ce qui explique sa présence dans les outils de conception et d’optimisation de nombreux groupes. Adkins et Liebeck l’ont replacée dans le cadre de la conception à pertes induites minimales; JBLADE l’étend par un modèle de « three-dimensional flow equilibrium »; McCrink et Gregory l’adaptent explicitement aux petits propulseurs électriques à faible Reynolds avec corrections de pertes de bout, Mach, composantes 3D et scaling Reynolds.

Ses limites, en revanche, sont très nettes aux faibles Reynolds. Premièrement, la BEMT n’est jamais meilleure que les polaires 2D qu’on lui fournit; or c’est précisément là que le faible Reynolds est le plus capricieux. Deuxièmement, l’hypothèse d’annuli quasi indépendants devient plus fragile près du moyeu, près du tip et en présence d’interactions, de cross-flow fort ou de déformations. Troisièmement, dès qu’une section entre dans un régime de décrochage profond, les fermetures standard Cl(α)C_l(\alpha), Cd(α)C_d(\alpha) deviennent peu fiables; c’est la raison pour laquelle JBLADE introduit explicitement un modèle de traînée à 9090^\circ pour mieux décrire le post-stall. Enfin, des données de profils valables à Re106Re\sim 10^6 introduites dans une BEMT pour des hélices travaillant à Re104Re\sim 10^410510^5 donnent des rendements irréalistes, comme le montre le cas Chalmers.

Corrections de pertes de bout et de moyeu

Dans la pratique, la BEMT pour hélices de drones s’appuie presque toujours sur la famille de corrections Prandtl/Goldstein. La correction de Prandtl représente l’atténuation de circulation au voisinage du tip par un facteur de perte FtipF_{tip}, classiquement appliqué au terme de quantité de mouvement annulaire; une correction analogue côté moyeu est souvent notée FhubF_{hub}, et l’on prend souvent un facteur combiné F=FtipFhubF=F_{tip}F_{hub}. Les approches de Goldstein, plus cohérentes pour le « minimum induced loss », sont préférables en théorie pour les propulseurs optimisés mais plus lourdes; Adkins–Liebeck ont justement reconstruit une conception optimale compatible avec l’analyse, et une large partie des codes modernes reste, au fond, un compromis entre Prandtl simple, Goldstein plus exact, et fermetures empiriques utiles en conception préliminaire. Pour les faibles Reynolds des drones, la hiérarchie des erreurs est cependant claire : une correction de tip légèrement imparfaite est souvent moins pénalisante qu’une polar 2D erronée par mauvaise transition.

Place des calculs constructeurs

Les calculs fabricants, et spécialement ceux d’APC, ne doivent être ni rejetés ni sacralisés. APC indique explicitement que ses fichiers de performance sont calculés, non mesurés, avec une méthode spécialisée fondée sur la théorie vortex et la géométrie réelle des pales; la société précise aussi que les résultats sont surtout utiles sur une base comparative, lorsqu’ils sont « ancrés » par des données expérimentales d’un design voisin. APC signale par ailleurs une mise à jour importante de son module en 2022 pour améliorer la précision sur une plage plus large de Reynolds. Autrement dit, les données constructeur sont très utiles en présélection, benchmarking et comparaison relative, mais elles ne remplacent pas une mesure UIUC/OSU/UBI lorsqu’on veut faire de la validation méthodique.

Recommandations empiriques

Les recommandations qui suivent ne sont pas des lois de similitude exactes; ce sont des règles de conception/study de haute confiance, issues du croisement expériences–CFD–modèles.

Premièrement, il faut concevoir au Reynolds local, pas au Reynolds « nominal » de l’hélice. La décision de choix de profil, de trip, de rugosité ou de flap doit se faire à la station radiale et au point de fonctionnement réellement visés. Les résultats de Deters, le benchmark acoustique de Grande, et l’optimisation de Yao convergent tous : ce qui marche à Re4×104Re\approx 4\times10^4 peut être neutre ou pénalisant à 1.5×1051.5\times10^52×1052\times10^5.

Deuxièmement, pour améliorer la poussée et le rendement à faible Reynolds, augmenter le Reynolds local est souvent le levier le plus robuste. En pratique, cela signifie souvent éviter les vitesses de rotation trop faibles sur les petites hélices, ne pas sous-charger excessivement la pale, et éviter les cordes trop faibles là où la vitesse tangentielle est encore modeste. Les données UIUC sur l’APC 10×7 et les campagnes Deters montrent toutes une hausse du CTC_T quand rpm et/ou ReRe montent, toutes choses égales par ailleurs.

Troisièmement, le choix du profil doit être faible-Re explicite. APC indique d’ailleurs qu’elle utilise pour ses hélices électriques et multirotors un mélange E63/Clark-Y-like, pas des profils « high-speed » génériques; à Chalmers, remplacer des données NACA 16 valables à Re>7×105Re>7\times10^5 par des polaires générées à bas Reynolds a réduit l’écart BEMT–CFD d’environ 10 points de rendement.

Quatrièmement, le vrillage reste essentiel, non seulement pour distribuer la charge mais aussi pour contrôler simultanément ReRe, incidence effective et Mach local. L’exemple APC 10×7 montre un washout marqué, et APC précise que le washout et la réduction de cambrure sont utilisés lorsque les effets de Mach deviennent significatifs près du tip. Une erreur classique en projet étudiant consiste à garder le bon diamètre et le bon pas « moyen » tout en négligeant que la torsion réelle fait une grande partie du travail aérodynamique.

Cinquièmement, la rugosité doit être traitée comme une variable de conception ciblée, pas comme un défaut toujours nuisible ni comme un trip toujours bénéfique. À Rec=6500Re_c=6500, Cooke et al. obtiennent des gains modestes de poussée et de puissance; mais l’ensemble de la littérature faible-Re sur transition/LSB, plus les observations émergentes sur l’érosion, invite à ne pas extrapoler ce bénéfice à des hélices de drones conventionnelles. Pour des drones de terrain, la question correcte n’est pas « faut-il une surface rugueuse ? », mais « où, à quelle échelle et pour quel segment de mission la rugosité modifie-t-elle favorablement la transition ? ».

Sixièmement, sur les pales souples ou composites, la simulation aéroélastique couplée n’est pas indispensable pour toutes les petites hélices, mais elle devient pertinente dès qu’on combine grand diamètre relatif, faible disc loading, matériaux souples ou planforms balayés. Möhren et al. sont très clairs sur ce point, et leurs conclusions sont cohérentes avec la littérature composite sur le bend-twist passif.

Résultats robustes, questions ouvertes et amélioration d’un cours universitaire

Résultats robustes

Le premier résultat robuste est que les petites hélices de drones opèrent sur une plage de Reynolds bien plus basse que les hélices aéronautiques classiques, souvent Rec104Re_c\sim 10^410510^5, avec une forte hétérogénéité radiale. Le deuxième est que l’augmentation de ReRe améliore généralement CTC_T et souvent le rendement maximal, toutes choses égales par ailleurs. Le troisième est que transition, LSB et séparation locale sont des mécanismes réels et mesurables, pas des raffinements théoriques. Le quatrième est que le disque actuateur et la BEMT restent indispensables, mais seulement si l’on explicite leurs hypothèses et si l’on fournit des polaires compatibles avec le régime. Le cinquième est que les calculs constructeur sont utiles pour comparer, pas pour clore une validation.

Points encore débattus

La meilleure manière de prédire la transition sur une hélice réelle très petite reste débattue. XFOIL, γReθ\gamma-Re_\theta, transition SST, URANS et LES ont chacun des domaines d’utilité, mais aucun n’est universel sur Re<2×104Re<2\times10^4. L’effet de la rugosité n’est pas universalement négatif ni positif. Le traitement du post-stall reste délicat en BEMT, d’où le recours à des extensions spécifiques comme JBLADE. Et, pour les petites hélices de drones usuelles, l’importance pratique de l’aéroélasticité reste très dépendante du matériau et de la rigidité réelle de fabrication. Enfin, la part respective des erreurs venant des polaires 2D, des pertes de bout/moyeu et du modèle d’inflow reste encore très configuration-dépendante, surtout en vol de translation ou en interaction multi-hélices.

Erreurs pédagogiques fréquentes

La plus fréquente consiste à parler d’« un Reynolds de l’hélice » sans jamais montrer sa variation radiale. La deuxième consiste à enseigner CT(J)C_T(J) et CP(J)C_P(J) comme s’ils étaient des courbes purement géométriques, alors qu’elles bougent avec rpm via le Reynolds. La troisième est d’utiliser une BEMT avec des polaires de profils mesurées à Re106Re\sim 10^6 pour analyser une hélice qui travaille à Re3×104Re\sim 3\times10^4. La quatrième est de présenter le disque actuateur comme si c’était un modèle prédictif complet plutôt qu’un modèle de bilan global. La cinquième est d’ignorer les corrections de montage, les incertitudes et la différence entre données mesurées UIUC et fichiers calculés constructeur. La sixième est d’enseigner le décrochage de section sur une hélice comme s’il se produisait partout de manière simultanée, alors qu’il est fortement localisé et couplé à la variation de ReRe, d’incidence et de vitesse hélicoïdale.

Recommandations concrètes pour améliorer un cours universitaire

Un bon cours sur les hélices de drones devrait introduire très tôt un cas fil rouge unique, par exemple l’APC 10×7 UIUC, et demander aux étudiants de reconstruire eux-mêmes c(r)c(r), β(r)\beta(r), Re(r)Re(r), M(r)M(r), CT(J)C_T(J), CP(J)C_P(J), CQ(J)C_Q(J) et η(J)\eta(J) à 3000, 4000, 5000 et 6000 rpm. Cela force à voir que le Reynolds est une distribution, pas un scalaire.

Le cours devrait ensuite séparer explicitement quatre niveaux de modèle : disque actuateur, BEMT, BEMT corrigée, CFD transitionnelle. L’objectif n’est pas d’empiler les méthodes, mais de faire identifier à chaque niveau ce qu’on gagne et ce qu’on perd : bilans énergétiques au niveau disque; distribution de charge en BEMT; traitement du tip/hub et du Reynolds en BEMT corrigée; séparation et transition en CFD. Les étudiants comprennent alors beaucoup mieux pourquoi une BEMT « bien nourrie » peut battre une CFD « mal fermée ».

Il faut aussi introduire un mini-module de transition à partir d’une section de pale typique. Par exemple : prendre un profil faible-Re, comparer une polar « naturelle » à une polar avec transition forcée, puis montrer au niveau hélice l’effet de cette différence sur la courbe η(J)\eta(J) et sur le bruit LSB. Grande et Yao donnent exactement les matériaux pédagogiques nécessaires pour cela.

Une séance de comparaison critique de sources est également essentielle : UIUC comme mesure de référence; banc UBI/MAAT pour l’incertitude et la répétabilité; APC comme calcul constructeur vortex; Kutty comme exemple de CFD utile mais limitée; Chalmers comme exemple où la qualité des polaires contrôle davantage le résultat que la sophistication du solveur. C’est cette comparaison qui évite les trois erreurs classiques des étudiants : « la CFD a toujours raison », « le constructeur a forcément raison », ou « la BEMT ne sert qu’au pré-dimensionnement ».

Questions ouvertes et limites de cette synthèse

Les jeux de données les plus riches restent encore concentrés sur certaines familles d’hélices commerciales, surtout APC. Pour plusieurs études récentes, les géométries détaillées complètes et les incertitudes point par point ne sont pas toujours ouvertes. Les effets de rugosité/érosion sur les hélices « drone » usuelles 5–12 pouces restent moins bien établis que sur les micro-propellers ou sur les profils 2D. Enfin, une partie de la littérature BEMT la plus utile demeure derrière paywall, ce qui limite l’extraction fine de certaines équations et validations. Les conclusions ci-dessus sont donc robustes sur les tendances, mais plus prudentes sur la hiérarchie exacte des erreurs selon le code ou la fermeture utilisés.

Bibliographie primaire

Les références ci-dessous sont volontairement centrées sur les sources primaires les plus utiles pour l’enseignement, la validation et la conception.

Brandt, J. B., and Selig, M. S., Propeller Performance Data at Low Reynolds Numbers, AIAA 2011-1255. DOI: 10.2514/6.2011-1255. Source stable UIUC: https://m-selig.ae.illinois.edu/pubs/BrandtSelig-2011-AIAA-2011-1255-LRN-Propellers.pdf

Deters, R. W., Ananda, G. K., and Selig, M. S., Reynolds Number Effects on the Performance of Small-Scale Propellers, AIAA 2014-2151. Source stable UIUC: https://m-selig.ae.illinois.edu/pubs/DetersAnandaSelig-2014-AIAA-2014-2151.pdf

Dantsker, O. D., Caccamo, C., Deters, R. W., and Selig, M. S., Performance Testing of APC Electric Fixed-Blade UAV Propellers, AIAA 2022-4020. Source stable UIUC: https://m-selig.ae.illinois.edu/pubs/DantskerCaccamoDetersSelig-2022-AIAA-Paper-2022-4020-APCEPropTestsVol4.pdf

UIUC Propeller Database, Vols. 1–4. URL stable: https://m-selig.ae.illinois.edu/props/propDB.html

Merchant, M. P., Propeller Performance Measurement for Low Reynolds Number Unmanned Aerial Vehicle Applications, M.S. Thesis, Wichita State University, 2005. URL stable: https://soar.wichita.edu/server/api/core/bitstreams/6ff13603-cc8e-4c29-be85-cd196639090a/content

Silvestre, M. A. R., Morgado, J., Alves, P., Santos, P., Gamboa, P., and Páscoa, J. C., Propeller Performance Measurements at Low Reynolds Numbers, 2015. URL stable: https://www.naun.org/main/NAUN/mechanics/2015/a372003-136.pdf

McCrink, M. H., and Gregory, J. W., Blade Element Momentum Modeling of Low-Reynolds Electric Propulsion Systems, Journal of Aircraft, 2017. DOI: 10.2514/1.C033622

Morgado, J., Silvestre, M. A. R., and Páscoa, J. C., Validation of New Formulations for Propeller Analysis, Journal of Propulsion and Power, 2015. DOI: 10.2514/1.B35240

Adkins, C. N., and Liebeck, R. H., Design of Optimum Propellers, Journal of Propulsion and Power, 1994. DOI: 10.2514/3.23779. Source stable PDF: https://www.humanpoweredflight.co.uk/hpfMedia/media/7/Adkins-and-Liebeck-paper.pdf

Grande, E., Ragni, D., Avallone, F., and Casalino, D., Laminar Separation Bubble Noise on a Propeller Operating at Low Reynolds Numbers, AIAA Journal, 2022. DOI: 10.2514/1.J061691

Kutty, H. A., and Rajendran, P., 3D CFD Simulation and Experimental Validation of Small APC Slow Flyer Propeller Blade, Aerospace, 2017, 4(1), 10. DOI: 10.3390/aerospace4010010

Yao, Y., Ma, D., Zhang, L., Yang, X., and Yu, Y., Aerodynamic Optimization and Analysis of Low Reynolds Number Propeller with Gurney Flap for Ultra-High-Altitude Unmanned Aerial Vehicle, Applied Sciences, 2022, 12, 3195. DOI: 10.3390/app12063195

Cooke, J. P., Campbell, M. F., Steager, E. B., Bargatin, I., Yim, M. H., and Park, G. I., Numerical and Experimental Study on the Addition of Surface Roughness to Micro-Propellers, Physics of Fluids, 2023, 35, 113607. DOI: 10.1063/5.0176690

Möhren, F., et al., On the Influence of Elasticity on Propeller Performance: A Parametric Study, CEAS Aeronautical Journal, 2023. DOI: 10.1007/s13272-023-00649-y

Grava, A., Serpieri, J., Iuso, G., Bernardos, L., and Cafiero, G., Experimental Investigation of a Small Drone Propeller Aerodynamics in Forward Flight, AIAA Journal, 2025, 63(5), 1855–1866. DOI: 10.2514/1.J064664. Preprint stable: https://iris.polito.it/retrieve/6d7414ac-105e-4112-a427-0e8ee7a5d780/GSIBC_AIAAJ_preprint.pdf

Comunian, P., Montagner, S., Serpieri, J., and Cafiero, G., Aerodynamic Interaction of Side-by-Side Propellers in Forward Flight Operating at Low-Reynolds Numbers, Aerospace Science and Technology, 2026, 171, 111577. DOI: 10.1016/j.ast.2025.111577. Preprint stable: https://iris.polito.it/retrieve/handle/11583/3006550/5efccf4e-c4ea-412d-b705-6fcef1dd6c97/1-s2.0-S1270963825016414-main.pdf

APC Propellers, Performance Data and Engineering pages. URLs stables: https://www.apcprop.com/technical-information/performance-data/ et https://www.apcprop.com/technical-information/engineering/

MaxFly

Delivery drones people accept. Built in the open, flight after flight.

Follow the R&D
made with by cmrabdu R&D IN THE OPEN · MMXXVI