Machine Translated by Google

Journal des fluides et des structures 49 (2014) 574–593

Listes de contenu disponibles sur ScienceDirect

Journal des fluides et des structures

page d'accueil du journal : [Link]/locate/jfs

Calcul des coefficients de masse hydrodynamique et d'amortissement pour

un écoulement d'hélice marine cavitant à l'aide d'une méthode de panneau
M. Gaschler n , M. Abdel­Maksoud 1,2
Institut de dynamique des fluides et de théorie des navires, Université technique de Hambourg, Schwarzenbergstraße 95 C, 21073 Hambourg, Allemagne

Informations sur l'article abstrait

Historique de l'article : Le présent article traite du calcul numérique de la masse hydrodynamique et des coefficients d'amortissement
Reçu le 19 juin 2013 dans le cadre de la cavitation en nappe instationnaire sur les écoulements d'hélices marines. Dans la première
Accepté le 1er juin 2014
partie de l'article, le contexte mathématique et numérique derrière la méthode numérique est présenté. Les
Disponible en ligne le 1er juillet 2014
calculs numériques effectués dans ce travail sont basés sur une méthode de panneau d'ordre faible. Les
Mots clés: méthodes de panneau appartiennent à la classe des techniques de collocation et sont appliquées pour obtenir
Masses hydrodynamiques et amortissement une solution numérique d'un système d'équations intégrales aux limites basé sur un écoulement potentiel. Elles
Cavitation en feuille
sont adaptées à la présente application en raison de leur temps de calcul court qui les rend applicables dans le
Méthode des éléments de frontière
processus de conception d'hélices marines.
Méthode du panneau
Flux potentiel
Cadres de référence multiples De plus, deux approches différentes pour la détermination des masses hydrodynamiques et
de l'amortissement sont introduites dans ce travail. Les masses hydrodynamiques et l'amortissement
sont importants dans les études du mouvement des navires en mer et dans l'analyse des vibrations
d'un navire et de ses appendices. Les approches développées sont appliquées à un écoulement
d'hélice en mouvement de pilonnement. Les calculs sont effectués pour une hélice non rotative et
rotative dans des conditions de cavitation et de non­cavitation. Les résultats obtenus à partir des
simulations sont discutés en détail et une perspective est donnée.
& 2014 Elsevier Ltd. Tous droits réservés.

1. Introduction

Cet article est consacré à l'étude numérique de la masse hydrodynamique et des coefficients d'amortissement d'une hélice marine en tenant compte de la
cavitation instationnaire en nappe partielle. La cavitation se développe en général dans un écoulement liquide lorsque la pression dans l'écoulement tombe en
dessous de la pression de vapeur du liquide. Ce phénomène est principalement causé par des vitesses locales élevées du fluide à proximité des pales de l'hélice.
Il peut exister différents types de cavitation, par exemple la bulle, le tourbillon de pointe, la cavitation en nappe partielle et la supercavitation. Dans ce travail, seule
l'étude de la cavitation en nappe partielle est abordée. La cavitation en nappe partielle est un type de cavitation qui se produit sous la forme d'une nappe de vapeur
attachée à la surface du corps solide. C'est l'une des principales causes des fluctuations de pression élevées sur la coque d'un navire et peut entraîner des
vibrations et du bruit. De ce fait, la prédiction de la cavitation en nappe sur les pales d'hélice et l'estimation de ses effets sur les performances de l'hélice sont l'un
des aspects les plus intéressants du processus de conception d'hélice. Étant donné que les expériences sur modèle sont difficilement applicables au stade de la
conception, la simulation numérique

n
Auteur correspondant. Tél. : þ49 40 42878 6147 ; Fax : þ49 40 42878 6055.
Adresses e­mail : [Link]@[Link] (M. Gaschler), [Link]­maksoud@[Link] (M. Abdel­Maksoud).
URL : [Link] (M. Gaschler), [Link] (M. Abdel­Maksoud).
1Tél. : þ49 40 42878 6053.
2 Fax : þ49 40 42878 6055.

[Link] 0889­9746/&
2014 Elsevier Ltd. Tous droits réservés.
Machine Translated by Google

M. Gaschler, M. Abdel­Maksoud / Journal des fluides et des structures 49 (2014) 574–593 575

Nomenclature ðt1;t2;t3Þ base vectorielle du système de coordonnées local ajusté à la

surface
D diamètre de l'hélice V vitesse totale
F vecteur de force vitesse induite
g constante de gravité vitesse non perturbée de vind
V1
coefficient d'avance de l'hélice Vitesse d'entrée parallèle Vinflow
Jk coefficient de poussée de l'hélice Navire V
vitesse du navire
M vecteur de moment Champ de sillage nominal du navire Vwake
n vitesse de rotation de l'hélice x variable d'espace
n vecteur normal z1 distance à la surface de l'eau libre

p pression locale ou épaisseur de cavitation

pv pression de vapeur de l'eau M force dipolaire

p1 pression atmosphérique densité de l'eau
m coordonnées des coordonnées ajustées à la surface locale p force de la source/du puits
système nate σv nombre de cavitation
Φ potentiel total
sd:p: point de début de la cavitation en feuille
point de fermeture de la cavitation de la feuille Potentiel induit
sd:p:
S limite du domaine d'écoulement Ω Φ1 potentiel non perturbé
SB surface du corps Domaine de flux de potentiel Ω
Sud­Ouest surface de sillage Oh vitesse angulaire
SC surface de cavitation fnW n­ième dérivée de la fonction f par rapport à
partie cavitante de la surface du corps la variable W
SBC
t variable temporelle

Des outils peuvent être utilisés pour prédire la cavitation en feuille sur une hélice. Ces outils ne doivent pas seulement être fiables et stables d'une part
mais elles doivent également être rapides pour pouvoir être utilisées dans la procédure de conception. Un type de méthodes qui répondent à ces exigences sont
les méthodes de panel. Dans ce travail, l'outil de simulation pan MARE basé sur une méthode de panel tridimensionnel d'ordre faible est
utilisé (voir Bauer et Abdel­Maksoud, 2012).
Dans la section 2, les équations régissant le problème d'écoulement ainsi que les détails du modèle de cavitation en feuille sont présentés.
introduit. La section 3 présente le schéma numérique et décrit les spécificités du traitement des termes instables dans le
modèle de cavitation. Dans la section 4, deux approches différentes pour la prédiction de la masse hydrodynamique et des coefficients d'amortissement
sont introduits. Dans la dernière section, les capacités de la méthode numérique sont démontrées sur un exemple d'hélice pour
différents cas d'application. L'écoulement de l'hélice est simulé en mouvement de soulèvement soumis à un flux entrant homogène à la fois
dans des conditions de non­cavitation et de cavitation et les coefficients de masse hydrodynamique et d'amortissement sont calculés par
moyens des approches développées.

2. Description mathématique

Cette section présente les équations gouvernantes qui sont utilisées pour la mise en œuvre du schéma numérique dans le
code de simulation pan MARE. L'accent est mis ici sur la caractérisation de différents référentiels et modèles de mouvement de
les corps immergés dans l'écoulement ainsi que dans la description mathématique du phénomène de cavitation en nappe instable.

2.1. Systèmes de coordonnées et modèles de mouvement

Dans ce travail, l'écoulement au­delà d'un ou de plusieurs corps avec une vitesse de rotation et/ou de translation individuelle est
considéré. A cet effet, deux référentiels différents sont définis. Un référentiel cartésien qui est fixé dans le
l'espace et un référentiel cartésien se déplaçant avec le corps considéré. Le référentiel fixe dans l'espace est
décrit par les coordonnées X ¼ ðX; Y; ZÞ, tandis que le référentiel cartésien à corps fixe est identifié par les coordonnées
x ¼ ðx; y; zÞ et l'axe z est positif vers le haut (voir Fig. 1).
Pour le traitement du modèle de cavitation en feuille, qui sera présenté plus tard, il est plus pratique d'utiliser une surface ajustée locale.
système de coordonnées curvilignes. Ainsi, un système de coordonnées fixe supplémentaire avec des vecteurs unitaires de base non orthogonaux
t1, t2, t3 et les coordonnées locales s ¼ ðs1; s2; s3Þ sont introduites (voir Fig. 2). Un vecteur arbitraire b est converti du système de coordonnées
ajusté à la surface au système cartésien global fixé au corps en utilisant les transformations :

bglob ¼ t1 t2
bs1 ð Þ
bs2
t1 þ
t2 bs2 ð Þ t1 þt3bs3
bs1 t2 þx0 ; ð1Þ
2 2
jt1 t2j jt1 t2j

où x0 décrit l'origine cartésienne et bglob est le vecteur transformé dans les coordonnées cartésiennes fixes du corps global.
Machine Translated by Google

576 M. Gaschler, M. Abdel­Maksoud / Journal des fluides et des structures 49 (2014) 574–593

Fig. 1. Systèmes de coordonnées fixes de l'espace et du corps.

Fig. 2. Systèmes de coordonnées fixes de surface locale et de corps global.

Pour chaque corps considéré, sa vitesse de translation et/ou sa vitesse angulaire ainsi que l'axe autour duquel il tourne doivent être
défini. Pour une hélice, le système de référence fixe du corps se déplace avec une vitesse d'entrée Vinflow et tourne avec la
vitesse angulaire ω, telle que le mouvement global soit donné par

V1ðxÞ ¼ VinflowðxÞþω x; ð2Þ

où x est le vecteur de position. La vitesse d'entrée de l'hélice derrière le navire peut être uniforme ou non uniforme.
Si le flux entrant est uniforme et que la fraction de sillage est négligée, alors la vitesse du flux entrant est égale à la vitesse du navire
T
VinflowðxÞ¼ðVship; 0 ; 0Þ . Pour un afflux non uniforme, la vitesse d'afflux est la différence entre la vitesse du navire et la vitesse du courant.
T
champ de sillage nominal du navire VinflowðxÞ¼ðVship; 0; 0Þ Réveille­toi.

Dans le cas d'une vitesse de mouvement constante du corps, dans le référentiel cartésien fixé au corps, la particule fluide
le mouvement ne dépendra pas du temps, c'est­à­dire que le flux sera indépendant du temps. Étant donné que la vitesse d'entrée de l'hélice n'est pas
constante en général, un écoulement instable est considéré dans ce travail.

2.2. Équations gouvernantes

Pour la description mathématique du problème d'écoulement, un domaine de calcul Ω est défini et les équations régissant
sont mis en place dans le référentiel mobile. Dans le domaine d'écoulement Ω l'écoulement est considéré comme incompressible, non visqueux
et irrotationnel. Pour un écoulement non visqueux idéal, les équations de Navier­Stokes pour la masse et l'impulsion se simplifient en l'équation de Laplace tridimensionnelle et l'équation
de Bernoulli instable (voir Katz et Plotkin, 2001, p. 26f.) :

∂2ΦðxÞ ∂2ΦðxÞ ∂2ΦðxÞ

2Φð Þ¼ x c'est à dire c'est à dire
¼ 0; ð3Þ
∂2x ∂2 ans ∂2z

2 2 ∂ ðxÞ
pð Þ¼ x p1 j thr
Φ1ð Þj x
2
j Φð Þj x r Þρg zð Þ 1 ∂t Avec
; ð4Þ
Machine Translated by Google

M. Gaschler, M. Abdel­Maksoud / Journal des fluides et des structures 49 (2014) 574–593 577

où xAΩ est la variable spatiale, t la variable temporelle, Φ ¼ þΦ1 est le potentiel de vitesse total qui est composé du potentiel perturbé et du potentiel de courant libre
non perturbé Φ1, et p est la pression. Les constantes p1, ρ, g et z1 sont respectivement la pression atmosphérique, la densité de l'eau, la constante de gravité et la
distance à la surface de l'eau libre.
Par la troisième identité de Green, la solution de l'équation de Laplace peut être obtenue sous la forme d'une distribution des flux source et dipôle
sur la surface du corps, ajouté au flux non perturbé (voir Katz et Plotkin, 2001, p. 44f.)
∂ 1 σðxÞ
µð Þ x dSðÞx ; ð5Þ
ε ð Þ¼ x0 Z S ∂n rðx0; xÞ rðx0 ; xÞ

où S ¼ ∂Ω est la frontière du domaine d'écoulement, r ¼ Jxx0 J est la fonction de distance et σ ð∂Φ=∂n∂Φinner=∂nÞ et μ ðΦΦinnerÞ sont respectivement les forces
source et dipolaire (voir Katz et Plotkin, 2001, p. 47). Le potentiel Φinner décrit le potentiel interne à l'intérieur d'un corps solide. Comme il n'y a pas de signification physique
pour un potentiel interne, il peut être choisi arbitrairement. Nous supposons Φinner ¼ Φ1 tel que l'on obtienne σ ¼ ∂ =∂n et μ¼ .

La constante ε dans l'équation (5) dépend de l'emplacement du point fixe x0 :

8 4π si x0 se trouve à l'extérieur de la limite du corps solide ;

e¼ 2π si x0 se trouve sur la limite du corps solide ; ð6Þ
><

>: 0 si x0 se trouve à l'intérieur de la limite du corps solide :

Afin de déterminer la distribution de la source et du dipôle à la surface du corps, une approche basée sur le potentiel est choisie dans ce travail. Dans cette
approche, les points de collocation sont situés à l'intérieur du corps solide. Les équations intégrales de limite résultantes se lisent comme suit

∂ 1 σðxÞ
µð Þ x dSð Þ¼ x 0; 8x0 AΩintérieur; ð7Þ
AVECS ∂n rðx0; xÞ rðx0 ; xÞ

où x0 sont des points à l'intérieur du corps solide. L'équation intégrale de limite (7) est utilisée pour déterminer la distribution des sources et des dipôles sur la surface du
corps. Ensuite, le potentiel total Φ peut être calculé en un point arbitraire du domaine d'écoulement Ω et la distribution de vitesse dans l'écoulement peut être obtenue en
différenciant le potentiel total : ΦðxÞ ¼ VðxÞ ¼ vindðxÞþV1ðxÞ; 8xAΩ;

ð8Þ

où vind décrit la vitesse induite par un corps et V1 est la vitesse de référence non perturbée. Afin d'obtenir une solution unique de l'équation intégrale de limite (7), des
conditions aux limites doivent être spécifiées. La limite S se compose en général de la limite du corps solide SB, de la surface de sillage SW et de la limite extérieure à
l'infini S1. Les conditions aux limites sont formulées dans la sous­section suivante.

2.3. Conditions aux limites

En dehors du domaine d'écoulement Ω, noté S1, les perturbations dues aux corps immergés dans l'eau sont négligeables : limjx x0j­1 ðxÞ ¼ 0.

Français À l'intérieur du domaine d'écoulement, les conditions aux limites doivent être spécifiées le long de toutes les limites des corps. En
supposant qu'il y ait plusieurs corps dans le domaine d'écoulement, par exemple un navire et une hélice, les limites dans le domaine d'écoulement
peuvent être exprimées comme une union de plusieurs surfaces de corps : [M m ¼ 1Sm, où M est le nombre total de corps. Pour chaque corps, trois
types de limites sont distingués : la surface du corps solide SB, la surface du sillage SW et la surface de cavitation en feuille SC. Dans cet article,
seule la cavitation en feuille partielle est prise en compte, ce qui signifie qu'aucune supercavitation n'est modélisée sur la surface du sillage SW et
que la surface de la cavité SC est modélisée comme une feuille de vapeur qui est attachée à la surface du corps solide.
Tout d’abord, les conditions aux limites sur les parties non cavitantes du corps et de la surface du sillage sont introduites :

Sur les parties non cavitantes du SB, la condition limite de Neumann est appliquée :

ΦðxÞ n ¼ 0 ; 8xASB ; ð9Þ

où le vecteur n représente le vecteur normal au point x. A partir de cette condition, on obtient une relation pour la force de la source sur la partie mouillée d'un corps
porteur :

∂Φ1ðxÞ
¼ V1 n; 8xASB: σð Þ¼ x ∂n ð10Þ

Sur la surface du sillage SW, le saut de la composante de vitesse normale ainsi que la différence de force entre les côtés inférieur et supérieur de la feuille de sillage
doivent disparaître :

pþ ðxÞp ðxÞ ¼ 0 et V ðxÞ nV ðxÞ n ¼ 0 ; 8xASW : ð11Þ

c'est à dire

Des relations ci­dessus, il résulte :

σréveil
∂ þ ∂ ð12Þ
∂n
∂n ! ¼ 0;
Machine Translated by Google

578 M. Gaschler, M. Abdel­Maksoud / Journal des fluides et des structures 49 (2014) 574–593

et la force dipolaire μwake ð þ Þ doit satisfaire à l'équation :

∂μréveil
1 2 2
j Φþ j j Φj ¼ 0: ð13Þ
∂t 2

La condition sans force (11) au bord de fuite du sillage est également connue sous le nom de condition de Kutta (voir Katz et Plotkin,
(2001, p. 210 et suiv.).

Maintenant, les conditions limites utilisées pour la détermination de la feuille de cavité sont présentées. Ces conditions limites peuvent
être appliqué soit sur la surface exacte de la cavité SC , soit sur les parties cavitantes correspondantes de la surface du corps, désignées par SBC .
Cette dernière approche est appelée approche partiellement non linéaire et est utilisée dans cet article (voir Vaz, 2005).
Pour la modélisation de la feuille de cavité, deux régions sont définies : la région où la cavitation commence et la région où
la cavitation se termine. Le point de départ de la cavitation est appelé le point de détachement et est désigné par sd:p:, la fermeture de la cavité
Le point est appelé point de rattachement et est noté sr:p:. Dans le présent travail, le point de détachement est estimé par la
Critère de Villat–Brillouin (voir Kinnas et al., 2003). Le point de réattachement est calculé par une procédure itérative qui
sera présenté dans la section 3. La quantité de cavitation sur une pale d'hélice est généralement caractérisée par l'absence de dimension
nombre de cavitation σn :

p1 pv
ð14Þ
σv p
:

2
2 ðVshipÞ

où pv est la pression de vapeur de l'eau, n est la vitesse de rotation et D est le diamètre de l'hélice considérée.
La cavitation en feuillets est modélisée dans la théorie des écoulements potentiels par deux conditions limites :

La première condition est la condition limite cinématique (voir Vaz, 2005, p. 12) :
D
Th s1 ; s2 ; s3 = 0 ; 8s ¼ ð Þ s1 ; s3ASBC
s2 ; ; ð15Þ
DTF

où η est l'épaisseur de la cavité et Fðηðs1; s2;tÞ; s3Þ ¼ s3 ηðs1; s2;tÞ est une fonction de la forme de la cavité. D'après l'équation (15) , on
obtient une équation aux dérivées partielles linéaire pour l'épaisseur inconnue de la cavité η :

∂η ∂η ∂η
Vs1 ð Þ Vs2
t1 t2 þ Vs2 ð Þ t1 t2 Vs1 ¼ jt1 t2j Vs3 ∂s2 ; 8sASBC ; ð16Þ
∂s1 ∂t

où Vsi ; i ¼ 1; 2; 3, sont les composantes de vitesse de la vitesse totale V dans le système de coordonnées non orthogonales locales.
La deuxième condition qui décrit la physique de la cavitation en feuillets est la condition limite dynamique (voir Vaz, 2005,
p. 12) :

p ¼ pv; 8sASBC : ð17Þ

En utilisant l'équation (4), la condition limite dynamique peut être transformée en une formulation de type Dirichlet pour le potentiel μ
sur la partie cavitante du corps :
ffiffi

fp V1;s1 ds1; 8sASBC ; ½jt1 t2j ð18Þ

µ ¼ µ0 Z s1

où μ0 ¼ μðsd:p:;tÞ est le potentiel au point de détachement de la cavitation et la fonction f est définie par

2 ∂μ
f ¼ jV1j þ ðnDÞ 2σn þ2g s3;1 s3 þ2 V2 V2s3 : ð19Þ
∂t s2

La composante normale de la vitesse locale Vs3 sera négligée dans les calculs numériques car elle n'a pas une valeur élevée.
impact sur l'ampleur de la force dipolaire mais peut provoquer des instabilités numériques (voir Vaz, 2005, p. 16).

3. Le schéma numérique

Cette section traite de la méthode numérique utilisée pour discrétiser le problème aux limites développé dans ce dernier
section. De plus, l'algorithme de prédiction de la forme de cavitation en feuille dans un écoulement d'hélice instable est décrit.

3.1. Modèle de cavitation à feuille discrète

Pour la simulation numérique, le solveur d'écoulement potentiel interne Pan MARE est utilisé. Ce programme est basé sur une méthode de panneau
tridimensionnel où les surfaces du corps et du sillage sont discrétisées en éléments quadrilatéraux plats et les
les équations intégrales de limite sont appliquées sur un point fixe de chaque élément du panneau, appelé point de collocation (voir Fig. 3).
Les points de collocation sont désignés ci­après par xj, 8j ¼ 1 ; …;N, où N est le nombre total de panneaux de carrosserie.
sont définis comme des points centraux des panneaux qui sont légèrement décalés à l'intérieur de la carrosserie.
Machine Translated by Google

M. Gaschler, M. Abdel­Maksoud / Journal des fluides et des structures 49 (2014) 574–593 579

Fig. 3. Corps d'hélice discrétisé et géométrie du sillage.

En considérant M corps dans le domaine fluide et en désignant Nm comme le nombre de panneaux pour chaque corps M, le nombre
… ; xi ; 8i total de panneaux est N ¼ ∑M m ¼ 1 Nm. Les points centraux du corps et des panneaux de sillage sont désignés par m ¼ 1 ;
¼ 1 ; … ; Nm ; Nm þ1 ; … ; Nsillage ; m, où Nm est le nombre de panneaux du corps et Nsillage ; m le nombre de panneaux de sillage sur
le corps m.
Sur chaque élément du panneau de carrosserie, une source et un dipôle sont répartis avec une intensité constante sur un panneau.
panneaux, seuls les dipôles sont distribués. La forme discrétisée des équations intégrales aux limites (7) se lit comme suit

M Nm Nm

∑ ∑ Un μn i;j ∑ σnà Bi;j

je ð20Þ
m¼1 je ¼ 1 je ¼ 1
" # ¼ 0 pour tout xj; j ¼ 1; …;N;

où μn μðxi;tnÞ et σn
je σðxi;tnÞ sont respectivement l'intensité discrète du dipôle et de la source pour le pas de temps discret tn et
je

8 Comme
réveille­toi;m

Ai;j þ l∑¼ 1 Al;j si le panneau i est sur le côté supérieur du

bord de fuite de la bande radiale k;
Un Nk
ð21Þ
i;j réveille­toi

Ai;j ∑ l ¼ 1 Al;j si le panneau i est sur le côté inférieur du

>>>>>>>>>>>>>><

bord de fuite de la bande radiale k;

>>>>>>>>>>>>>:
Ai;j else;

où

∂ 1
∂n rðxj; xiÞ dS xi ðÞ 8i ¼ 1; …;Nm þNréveil;m;
ð22Þ
Ai;j Axi ; xj ¼Z Panneau

1
dS xi ðÞ 8i ¼ 1 ; …;Nm; ð23Þ
B;j Bxi ; xj ¼Z Panneau rðxj; xiÞ

sont les fonctions d'influence qui décrivent l'influence dipolaire ou source du panneau i sur le panneau j. Dans la définition des coefficients d'influence
An , la forme linéaire de la condition de Kutta a été utilisée

μréveil ¼ μU μL ; ð24Þ

où μU et μL sont les forces dipolaires des panneaux de carrosserie qui sont situés soit du côté aspiration, soit du côté
pression du bord de fuite.
Maintenant, pour le calcul des forces discrètes du dipôle et de la source, la différenciation suivante est effectuée :

Sur les panneaux de carrosserie non cavitants, les forces de la source sont connues à partir de la relation (10) et les forces des dipôles
sont déterminées en résolvant le système linéaire d'équations (20).
Sur les panneaux de carrosserie à cavitation, les forces dipolaires sont calculées par la forme discrète de la condition limite
dynamique (18). Les forces sources sont déterminées en résolvant le système linéaire d'équations (20).
Machine Translated by Google

580 M. Gaschler, M. Abdel­Maksoud / Journal des fluides et des structures 49 (2014) 574–593

Sur les panneaux de sillage, les forces dipolaires peuvent être calculées en appliquant la condition linéaire de Kutta (24) après la
calcul des forces dipolaires du corps au bord de fuite.

En supposant que sur le corps m il y a des panneaux cavitants Ncav;m et des panneaux non cavitants Nnoncav;m ¼ ðNm Ncav;mÞ , la courbe linéaire
le système d'équations se lit comme suit

M Noncav;m Voiture ;
Un
je;j μnoncav;n
m∑ je∑ ∑
je Bi;jσcav;n
je

¼1 ¼1 je ¼ 1
" #

M Voiture ; Noncav;m
Un
je;j μcav;n
¼ ∑ m¼1 ∑
je ¼ 1
je
∑
je ¼ 1
Bi;jσnoncav;n
je
ð25Þ
" # 8 ans ¼ 1 ; …;UN;

où les parties connues sont placées à droite et les parties inconnues à gauche du système. Une fois ce système linéaire
système d'équations est résolu, les vitesses locales peuvent être calculées et la forme de la surface du sillage peut être
ajusté. Une description détaillée de l'algorithme de résolution pour l'alignement du sillage peut être trouvée dans Katz et Plotkin (2001,
p. 211 et suiv.).

Tableau 1
Données de l'hélice P1356.

Caractéristiques Valeur

Diamètre D 0,25 (m)

Rapport de moyeu 0,180
Inclinaison θ 31,86 (vous)
Nombre de lames z 5
Section de lame NACA66
Rapport de pas P0:7=D 0,9967
Rapport de surface AE=AO 0,8
Facteur d’échelle l 31,6
Point d'opération J 0,719

Tableau 2
Données d'entrée pour les cas d'application.

Caractéristiques Notation Valeur Unité

Vitesse de rotation n 9,5 1/s

Coefficient d'avance J 0,719
Amplitude de la vitesse dans la direction z W0 0,3011 MS
Vitesse constante dans la direction x U0 JnD¼1,707625 MS
Vitesse angulaire de l'hélice Oh 2πn 59:6903 5 1/s
Nombre de tours d'hélice par oscillation k
Vitesse angulaire du mouvement de soulèvement du navire ω0 ω=k 11:9381 2,88 1/s
Pas angulaire Droite toi
Pas de temps Δt Δθ=ð360nÞ 0:000842 2,25/2,5 m

Nombre de cavitation St

Tableau 3
Masse hydrodynamique sans dimension calculée et moment d'inertie ajouté pour l'hélice non cavitante et non rotative.

Fonction Coefficients

fW f W_ f2W f WW_ de 3 nous 3

X 0,03549 0,00025 0,27996 0,00214 0,0 0,00002

ET 0,03922 0,00207 0,0 0,00003 0,03922 0,00207
AVEC 0,05176 0,02078 0,0 0,00002 0,05177 0,02078
K 0,00269 0,00002 0,02126 0,00014 0,0 0,0
M 0,05155 0,00062 0,0 0,00002 0,05155 0,00062
N 0,03254 0,00842 0,0 0,00001 0,03254 0,00842
Machine Translated by Google

M. Gaschler, M. Abdel­Maksoud / Journal des fluides et des structures 49 (2014) 574–593 581

3.2. Détermination de la forme de la cavité de la tôle

La forme de cavitation de la feuille est déterminée par une procédure itérative.

Tout d'abord, le problème de valeur limite (20) est mis en place et résolu avec uniquement des panneaux non cavitants afin d'obtenir une
distribution de pression sur la surface du corps. La distribution de pression calculée est ensuite utilisée pour faire une première estimation de la

Fig. 4. Distribution de force sans dimension d'une hélice non rotative et non cavitante. (a) Distribution de force sans dimension dans la direction x. (b) Distribution de
force sans dimension dans la direction y. (c) Distribution de force sans dimension dans la direction z.
Machine Translated by Google

582 M. Gaschler, M. Abdel­Maksoud / Journal des fluides et des structures 49 (2014) 574–593

longueur de la cavité en appliquant le critère :

pirpv 8i ¼ 1 ; ...;N : ð26Þ

Les panneaux pour lesquels ce critère est rempli sont définis comme des panneaux cavitants. Dans l'étape suivante, les forces dipolaires sur le
les panneaux cavitants sont calculés en discrétisant la condition limite dynamique (18) :
n
μcav;n
je
¼ μn0 þμ~ je : ð27Þ

La valeur de μn 0 est extrapolée à partir des trois valeurs dipolaires devant le point de détachement de la cavitation en feuillets.
n
La valeur de μ~ est calculée en discrétisant l'intégrale de l'équation (18) au moyen de la règle trapézoïdale (voir Vaz, 2005, p. 89).
je

Fig. 5. Distribution de l'impulsion sans dimension d'une hélice non rotative et non cavitante. (a) Distribution de l'impulsion sans dimension dans la direction x.
(b) Distribution d'impulsion sans dimension dans la direction y. (c) Distribution d'impulsion sans dimension dans la direction z.
Machine Translated by Google

M. Gaschler, M. Abdel­Maksoud / Journal des fluides et des structures 49 (2014) 574–593 583

Les termes instables apparaissant ∂μn=∂t et Vn s2 dans la condition limite dynamique peut entraîner des problèmes numériques puisque
leurs valeurs ne sont pas connues et doivent être estimées par une approximation. La valeur de ∂μcav;n je =∂t peut être approximé par
en utilisant l'équation (18)

∂ 1 2
1 2
μcav;n þ ðnDÞ 2σn þ jVn1je j g z1 zn1 ; ð28Þ
∂t 2 jV1;ij 2
je je

alors que la valeur de Vn je est approximée par sa valeur du pas de temps précédent Vn je
Vn1
je
.

Dans l'étape suivante, le système linéaire d'équations (25) avec des panneaux cavitants et non cavitants est mis en place et résolu
numériquement et les nouvelles vitesses et pressions sur le corps cavitant et non cavitant et sur les panneaux de sillage sont calculées.
L'épaisseur de la cavité sur les parties cavitantes du corps est calculée en résolvant l'équation aux dérivées partielles ( 16).
l'approximation numérique des dérivées spatiales ∂=∂s1 et ∂=∂s2 un schéma aux différences centrales et pour la
Pour l'approximation de la dérivée temporelle, un schéma de différences inverses du premier ordre est utilisé.
Si l'épaisseur de la cavité calculée à la fermeture de la cavité est inférieure à une tolérance donnée ou si le signe de la cavité
l'épaisseur à la fermeture de la cavité change pour deux étapes d'itération successives, l'algorithme s'arrête. Si ce n'est pas le cas, un nouveau
L'estimation de la longueur de la cavité est effectuée et la procédure d'itération est répétée jusqu'à ce que le critère d'arrêt soit satisfait pour toutes les radiales.
rayures. Le maillage de la surface du corps n'est pas re­grillé pendant la procédure de résolution.

3.3. Caractéristiques de l'hélice

En post­traitement, les caractéristiques de l'hélice peuvent être calculées. Les forces et les moments sont calculés
numériquement par
Nprop Nprop

F¼∑ dSðxÞ ; M ¼ ∑ dSðxÞ ; ð29Þ

je ¼ 1 pini Z Panneau je ¼ 1 viens xiÞ Z Panneau

où Nprop est le nombre de panneaux sur la surface de l'hélice.

Rapports entre la surface de la cavité et le volume sur une hélice avec Nprop
moteur les panneaux cavitants sont déterminés par

Moteur ∑Nprop
je ¼ 1 RPanneau dSðxÞ Moteur ∑Nprop
je ¼ 1 ηi R Panneau
dSðxÞ
Acav ¼ ; Vcav ¼ ; ð30Þ
R2 R3

où R est le rayon de l'hélice.

4. Détermination des masses hydrodynamiques et de l'amortissement d'une hélice en mouvement de soulèvement

Dans cette section, un écoulement d'hélice en mouvement de soulèvement dans des conditions de non­cavitation ainsi que de cavitation est
examiné. Lors du soulèvement, un navire est accéléré vers le haut et vers le bas le long de l'axe vertical. Afin de modéliser ce mouvement,
l'hélice est autorisée à se déplacer vers le haut et vers le bas selon un mouvement sinusoïdal prescrit du navire le long de l'axe z :

VinflowðtÞ¼ðU; V ; WÞ¼ðU0; 0 ; W0 péché ðω0tÞÞ : ð31Þ

On suppose que la vitesse angulaire ω0 est choisie de telle sorte que pendant un mouvement de haut en bas du navire,
Des rotations d'hélice kAN sont effectuées.
L'objectif est de dériver des formules pour le calcul des coefficients de force et d'impulsion hydrodynamiques pour un
hélice. Les masses hydrodynamiques et l'amortissement sont importants dans les études du mouvement des navires en mer. Hydrodynamique
Les masses sont des masses virtuelles qui sont utilisées pour décrire la force supplémentaire d'un corps accéléré dans un écoulement liquide stationnaire.
Pour un mouvement de corps instable, la traînée, qui est une force induite par la pression due à l'accélération d'un corps dans un écoulement liquide, peut
être inclus dans les calculs de flux potentiels.
En supposant que la vitesse du navire dans la direction x est constante U ¼ U0 const et que la vitesse du navire dans la direction z est a
fonction sinusoïdale dépendante du temps W ¼ W0 sin ðω0tÞ, on obtient les vecteurs de force et d'impulsion F ¼ ðX; Y; ZÞ et

Tableau 4
Masse hydrodynamique sans dimension calculée et moment d'inertie ajouté pour l'hélice non cavitante et rotative.

Fonction Coefficients

fW f W_ f2W f WW_ de 3 nous 3

X 0,03901 0,00047 0,30769 0,00372 0,0 0,0

ET 0,01407 0,00754 0,0 0,0 0,01407 0,00754
AVEC 0,16292 0,01935 0,0 0,0 0,16292 0,01935
K 0,00388 0,00005 0,03063 0,00042 0,0 0,0
M 0,03130 0,00320 0,0 0,0 0,03130 0,00320
N 0,09838 0,00817 0,0 0,0 0,09838 0,00817
Machine Translated by Google

584 M. Gaschler, M. Abdel­Maksoud / Journal des fluides et des structures 49 (2014) 574–593

M ¼ ðK; M;NÞ en fonction de W et de son accélération W_ :

X XðW ; W_Th K KdW ; W_Th

0 10 1 0 10 1
ET YðW ; W_Þ M MdW ; W_Th ð32Þ

B@
AVEC
CA ¼ B@ ZðW ; W_Þ QUE; B@
N CA ¼ B@ NordW; W_Th QUE:

À partir de cette hypothèse, deux approches différentes pour la détermination des masses ajoutées et de l'amortissement sont dérivées dans les
sous­sections suivantes.

Fig. 6. Distribution de force sans dimension d'une hélice rotative non cavitante. (a) Distribution de force sans dimension dans la direction x. (b) Distribution de force
sans dimension dans la direction y. (c) Distribution de force sans dimension dans la direction z.
Machine Translated by Google

M. Gaschler, M. Abdel­Maksoud / Journal des fluides et des structures 49 (2014) 574–593 585

4.1. Approche 1

La première approche suppose une approximation du second ordre par la série de Taylor d'une fonction f dans W et W_ on obtient , tel que

_ _
2
f W; W_ f 0 þΔWf DANS
þΔWf þ1 2 ΔW2f 2W þ2ΔWΔWf WW_
þΔW_ f 2W_ je; ð33Þ
DANS_
h

Fig. 7. Distribution de l'impulsion sans dimension d'une hélice rotative non cavitante. (a) Distribution de l'impulsion sans dimension dans la direction x.
(b) Distribution d'impulsion sans dimension dans la direction y. (c) Distribution d'impulsion sans dimension dans la direction z.
Machine Translated by Google

586 M. Gaschler, M. Abdel­Maksoud / Journal des fluides et des structures 49 (2014) 574–593

où ΔW ¼ W W0 ¼ W0ð sin ðω0tÞ1Þ et ΔW_ ¼ W_ W_ contenant une seule dérivée 0 ¼ W0ω0 cos ðω0tÞ. Seules les dérivées de vitesse et les termes
par rapport à l'accélération sont considérés, le terme f W_ W_ expansion se lit alors comme est négligé. Le résultat de Taylor

_ _

f W; W_ f 0 þΔWf DANS þΔWf DANS_

þ1 2 ΔW2f 2W þ2ΔWΔWf WW_
je

¼f 0 W0f DANS
þ1 2 W2 0f 2W Þ ω0t W0f DANS W2 0f 2W je
h h travailler péché ð h

þ cos ð þ ω0t hW0ω0f W_ W2 0ω0f WW_ travail cos ð Þ ω0t péché ðÞ ω0t W2 0ω0f WW_

2 1
þ péché ð Þ ω0t 2 W2 0f 2W ¼f 0 W0f W þ3 4 W2 0f 2W Þ ω0t W0ω0f W_ W2 0ω0f WW_ je
travail cos ð h

1 1
þ péché ð Þ ω0t W0f W2 0f 2W W2 0f 2W 2 W2 0ω0f WW_ je: ð34Þ
h Þ ω0t 4h
DANS
h h travail cos 2ð Þiþω0t
sin 2ð

Pour déterminer les coefficients f0, on effectue une analyseDANS ; …, la distribution de force est calculée par pan MARE pour plusieurs pas de temps et un
de Fourier f. L'analyse de Fourier fournit les coefficients de Fourier aN et bN d'un développement de Fourier :
1

fðtÞ ¼ a0 þ ∑ ½aN cos ðωNtÞþbN sin ðωNtÞ ; ð35Þ

N¼1

où ωN ¼ ω0N, NAN. Les termes d'accélération et de vitesse sont développés au second ordre :

fðtÞ a0 þa1 cos ðω1tÞþb1 sin ðω1tÞþa2 cos ðω2tÞþb2 sin ðω2tÞ
¼ a0 þþb1 cos ðω0tÞþa2 cos ð2ω0tÞþb2 sin ð2ω0tÞ : ð36Þ

Maintenant, on compare les coefficients des formules (34) et (36) et on obtient une relation entre la force et le moment
coefficients et les coefficients de Fourier calculés :

b14a2 a1 th2b2 4a2 2b2

f 0 ¼ a0 þb1 a2 ; f DANS
¼
; f W_
¼
;f 2W
¼
; f WW_
¼ :
ð37Þ
W0 W0o0 W20 W2 0ω0

4.2. Approche 2

Une approche mathématique alternative pour la détermination de la masse ajoutée inconnue et des coefficients d'amortissement est l'utilisation
des équations du mouvement. Les équations du mouvement d'une hélice en mouvement de pilonnement se lisent comme suit (voir Parsons et al., 1980, p. 2f.)

Max€ þDx_ ¼ Faire ; ð38Þ

T
où Ma est la matrice de masse ajoutée, D est la matrice d'amortissement et Fa ¼ ðX; Y; Z; K; M;NÞ est le vecteur qui décrit le
forces et moments hydrodynamiques supplémentaires dépendant du déplacement vertical, de la vitesse et de l'accélération de l'objet.
T
hélice. Le vecteur x décrit le déplacement ainsi que la rotation dans les directions x, y et z : x ¼ ðx; y; z; θx; θy; θzÞ . Dans notre

Tableau 5
Masse hydrodynamique sans dimension calculée et moment d'inertie ajouté pour l'hélice cavitante et rotative avec un nombre de cavitation σv¼2,5.

Fonction Coefficients

fW f W_ f2W f WW_ de 3 nous 3

X 0,04356 0,00052 0,34353 0,00409 0,0 0,0

ET 0,01346 0,00727 0,00002 0,00002 0,01346 0,00727
AVEC 0,16919 0,01941 0,00006 0,0 0,16918 0,01941
K 0,00609 0,00006 0,04802 0,00047 0,0 0,0
M 0,03201 0,00309 0,00002 0,00003 0,03201 0,00309
N 0,10201 0,00807 0,00006 0,0 0,10200 0,00807

Tableau 6
Masse hydrodynamique sans dimension calculée et moment d'inertie ajouté pour l'hélice cavitante et rotative avec un nombre de cavitation σv¼2,25.

Fonction Coefficients

fW f W_ f2W f WW_ de 3 nous 3

X 0,03700 0,00178 0,29203 0,01215 0,00002 0,00024

ET 0,01102 0,00702 0,00114 0,00009 0,01087 0,00703
AVEC 0,17094 0,02054 0,00475 0,00010 0,17034 0,02055
K 0,00619 0,00023 0,04885 0,00156 0,0 0,00004
M 0,03104 0,00318 0,00055 0,00001 0,03097 0,00318
N 0,10150 0,00863 0,00213 0,00009 0,10123 0,00864
Machine Translated by Google

M. Gaschler, M. Abdel­Maksoud / Journal des fluides et des structures 49 (2014) 574–593 587

exemple, seul le terme dépendant de la fréquence (déplacement vertical instable) est considéré, c'est­à­dire :

0 0 0 0
0 01 0 0 1 0 01 0 0 1
Avec_
WðtÞ z€ W_ ðtÞ
x_ ¼ ð39Þ
0 0 0 0
:

0 0 0 0
BBBBBBBB@
0 CCCCCCCCCA ¼ BBBBBBBB@
0 CCCCCCCCA et x€ ¼ BBBBBBBB@
0 CCCCCCCCCA ¼ BBBBBBBB@
0 CCCCCCCA

Fig. 8. Distribution de force sans dimension d'une hélice à cavitation pour deux nombres de cavitation différents : σv¼2,25 et σv¼2,5. (a) Distribution de force sans
dimension dans la direction x (b) Distribution de force sans dimension dans la direction y (c) Distribution de force sans dimension dans la direction z.
Machine Translated by Google

588 M. Gaschler, M. Abdel­Maksoud / Journal des fluides et des structures 49 (2014) 574–593

La partie dépendante de la vitesse dans l'équation (38) décrit la fonction d'amortissement, tandis que la partie dépendante de l'accélération
décrit les masses ajoutées. De cette considération, on obtient le système d'équations linéaires suivant :

m13W0ω0 cos ðω0tÞþd13W0 sin ðω0tÞ ¼ X

m63W0ω0 cos ðω0tÞþd63W0 sin ðω0tÞ ¼ N : ð40Þ

Fig. 9. Distribution de l'impulsion sans dimension d'une hélice à cavitation pour deux nombres de cavitation différents : σv¼2,25 et σv¼2,5. (a) Distribution de l'impulsion
sans dimension dans la direction x (b) Distribution de l'impulsion sans dimension dans la direction y (c) Distribution de l'impulsion sans dimension dans la direction z.
Machine Translated by Google

M. Gaschler, M. Abdel­Maksoud / Journal des fluides et des structures 49 (2014) 574–593 589

La comparaison des coefficients de l'équation (40) et de l'approximation de Fourier (36) fournit les coefficients inconnus de masse et d'amortissement :

a1
¼ mi3 ¼ b1 ; di3 ; je ¼ 1; …; 6; ð41Þ
W0o0 W0

où i désigne la ième composante du vecteur Fa. Cette approche est basée sur une approximation linéaire des distributions de
force et d'impulsion. En revanche, la première approche utilise une approximation quadratique.
Les coefficients peuvent être rendus sans dimension en les divisant par un terme approprié. Les coefficients sans dimension sont
désigné par l'exposant 0 et défini comme suit :

X0 X0 X0
1 XW
1 0 XW_ 1 0 2W 1 0 X2W 1
0 DANS
1 0 1 0 DANS_
1 1
Y0DANS EST Y0
EST_ Y02W Y2W ð42Þ
ρnD3 pD3 pD2
DANS_

B@ Z0DANS CA ¼ B@ ZW QUE; B@ Z0
DANS_
CA ¼ B@
ZW_ QUE; B@ Z02W CA ¼ B@ Z2W QUE;

et

X0
0 WW_
1 n 0 XWW_ 1
Y0
WW_ YWW_
ð43Þ
pD2
B@ Z0
WW_ CA ¼ B@
ZWW_ QUE:

Une mise à l'échelle analogue est effectuée pour les coefficients d'impulsion

K0 DANS K0 K0 2W
0 1 1 0
KW
1 0 DANS_
1 1 0 KW_
1 0 1 1 0 K2W
1
M0 MW M0
MW_ M0 2W M2W ð44Þ
rnD4 pD4 pD3
DANS DANS_

B@ N0 DANS CA ¼ B@ Nord­Ouest QUE; B@ N0

DANS_
CA ¼ B@
NO_ QUE; B@ N0 2W CA ¼ B@ N2W QUE;

Fig. 10. Formes de cavités en feuille et distributions de poussée de l'hélice P1356 à trois pas de temps différents. (a) Forme de cavité et coefficient de poussée à 1/4
d'oscillation (b) Forme de cavité et coefficient de poussée à 2/4 d'oscillation (c) Forme de cavité et coefficient de poussée à 3/4 d'oscillation.
Machine Translated by Google

590 M. Gaschler, M. Abdel­Maksoud / Journal des fluides et des structures 49 (2014) 574–593

Fig. 10. (suite)

et
K0WW_ KWW_
0 M0 1 n 0 1
WW_ MWW_ ð45Þ
N0WW_ pD3
B@ CA ¼ B@
NWW_ QUE:

5. Cas d'application : hélice P1356

Dans cette section, l'hélice P1356 est simulée en mouvement de soulèvement à l'échelle du modèle. L'hélice a été conçue pour le porte­conteneurs
KRISO (KCS) et est utilisée pour diverses études de recherche (voir KCS, 2011). Toutes les données pertinentes de l'hélice P1356 sont répertoriées dans
le tableau 1.
L'hélice est forcée de se déplacer au maximum de 101 vers le haut et vers le bas. L'amplitude de la vitesse angulaire et la
Le reste de toutes les données d'entrée pertinentes pour ce cas de test de simulation est résumé dans le tableau 2.
Les approches 1 et 2 pour la détermination des coefficients de masse hydrodynamique et de moment d'inertie, décrites
dans la section 4, sont appliquées à l'hélice P1356. Trois cas sont analysés : l'écoulement autour d'une hélice non rotative,
l'écoulement autour d'une hélice rotative et l'écoulement autour d'une hélice rotative et cavitante.

5.1. Résultats pour l'hélice non rotative et l'hélice non cavitante

Les résultats des simulations pour l'écoulement d'hélice non rotatif et non cavitant en mouvement de soulèvement sont illustrés dans le tableau
3 en termes de coefficients hydrodynamiques. Les résultats de la première approche sont répertoriés dans les quatre premières colonnes, tandis
que les résultats de la deuxième approche sont illustrés dans les deux dernières colonnes. À l'exception des coefficients des distributions de force
et de moment dans la direction x, les deux approches donnent des résultats presque identiques.
Machine Translated by Google

M. Gaschler, M. Abdel­Maksoud / Journal des fluides et des structures 49 (2014) 574–593 591

Fig. 10. (suite)

Français Les figures 4 et 5 représentent les diagrammes force–vitesse, impulsion–vitesse et force–accélération et impulsion–accélération. Ces diagrammes
montrent que la force et l'impulsion dans la direction x ne sont pas périodiques avec la fréquence d'oscillation mais avec le double de celle­ci. Les forces et les
moments dans les directions y et z sont périodiques avec la même fréquence que le mouvement de soulèvement. C'est pourquoi les coefficients f2W et f WW_
sont dominants pour la force et le moment dans la direction x et négligeables dans les directions y et z (voir tableau 3). La deuxième approche n'est pas en
mesure de calculer les coefficients f2W et f WW_ et est donc moins applicable pour le calcul des coefficients de force et de moment induits dans la direction x
en raison du mouvement dans la direction de transfert.

5.2. Résultats pour l'hélice rotative non cavitante

Les coefficients hydrodynamiques calculés pour une hélice rotative sont répertoriés dans le tableau 4 et illustrés dans les figures 6 et 7 de force et
d'impulsion . Pour le cas rotatif, la même considération que pour l'hélice non rotative peut être effectuée.
La distribution de la force et du moment dans la direction x est périodique avec le double de la fréquence d'oscillation.
Pour le cas rotatif, le coefficient d'amortissement pour la direction z est beaucoup plus élevé que pour le cas non rotatif.
en revanche, les masses ajoutées ne changent pas beaucoup.

5.3. Résultats pour l'hélice à cavitation rotative

Dans les tableaux 5 et 6, les coefficients hydrodynamiques calculés pour une hélice à cavitation avec deux cavitations différentes
les chiffres sont résumés. Les distributions de force et d'impulsion sont illustrées dans les figures 8 et 9.
Les figures 10a à c illustrent la forme et l’épaisseur de la feuille de cavité pour trois positions différentes des pales de l’hélice.
Bien que la surface et le volume de la cavité de la feuille soient plus élevés pour le nombre de cavitation le plus petit, il existe une différence négligeable entre
les distributions de force et de moment pour deux nombres de cavitation différents σ¼2,5 et σ¼2,25.
Machine Translated by Google

592 M. Gaschler, M. Abdel­Maksoud / Journal des fluides et des structures 49 (2014) 574–593

Fig. 11. Comparaison des masses hydrodynamiques et du moment d'inertie pour une hélice non rotative, rotative et cavitante. (a) Masses hydrodynamiques sans dimension. (b) Moment
d'inertie ajouté sans dimension.

La Fig. 11 montre une comparaison des coefficients de force et d'impulsion pour tous les cas simulés en termes de diagrammes à barres.
Pour l'hélice rotative, la différence entre les coefficients des cas cavitants et non cavitants est négligeable.
Les coefficients de l'hélice non rotative sont inférieurs aux résultats de l'hélice rotative.

6. Conclusion

Cet article est consacré à l'estimation numérique de la masse hydrodynamique et de l'amortissement des hélices en mouvement de
soulèvement. Deux approches différentes ont été présentées et appliquées à un exemple d'hélice. Pour cette hélice, trois cas différents ont été
étudiés : une hélice non rotative, une hélice rotative et une hélice cavitante. La cavitation en nappe instable semble avoir une influence négligeable
sur l'ampleur de la masse hydrodynamique et des coefficients d'amortissement.
A l'avenir, les résultats simulés devraient être validés par des données expérimentales. A cet effet, il est prévu de réaliser
mesures et de déterminer les masses hydrodynamiques et l'amortissement à partir des données mesurées.

Remerciements

Le premier auteur est soutenu par la Fondation allemande pour la recherche (DFG) au sein du groupe de formation à la recherche « Ports pour
« Les porte­conteneurs des générations futures ».

Références

Bauer, M., Abdel­Maksoud, M., 2012. Une méthode d'éléments de frontière basée sur le potentiel 3­D pour la modélisation et la simulation des écoulements d'hélices marines. Dans : The
Septième conférence de Vienne sur la modélisation mathématique, Vienne, Autriche, 2012.
Machine Translated by Google

M. Gaschler, M. Abdel­Maksoud / Journal des fluides et des structures 49 (2014) 574–593 593

Katz, J., Plotkin, A., 2001. Aérodynamique à basse vitesse. Cambridge University Press, Cambridge.
Kinnas, SA, Lee, HS, Young, YL, 2003. Modélisation de la cavitation en nappe instable sur les pales d'hélices marines. Revue internationale des machines tournantes 9,
263–277.
Porte­conteneurs KRISO (KCS), 2011. [Link] (consulté : novembre 2011).
Parsons, MG, Vorus, WS, Richard, E., 1980. Masse ajoutée et amortissement des hélices vibrantes. Rapport n° MA­RD­940­81037, Département de la Marine
Architecture et ingénierie marine, Collège d'ingénierie, Université du Michigan, Ann Arbor, Michigan, États­Unis.
Vaz, GNVB, 2005. Modélisation de la cavitation en feuille sur les hydroptères et les hélices marines à l'aide de méthodes d'éléments de frontière. Ingénierie aérospatiale Instituto Superior
Tecnico, Lisbonne, Portugal.