[{"data":1,"prerenderedAt":175},["ShallowReactive",2],{"$f3sk9al5nxruiq":3},{"slug":4,"publishedAt":5,"category":6,"title":7,"summary":8,"image":9,"imageAlt":14,"readingMinutes":15,"serviceId":16,"keywords":17,"imageCaption":23,"sections":24,"sources":165},"cfd-vs-fea","2026-10-25","methods","CFD ou FEA : à quelle question chacun répond, et quand il faut les deux","CFD ou FEA : pourquoi les fluides se calculent en volumes finis et les solides en éléments finis, où le transfert thermique les relie, et comment une erreur CFD devient une erreur de contrainte.",{"src":10,"thumbnail":11,"width":12,"height":13},"\u002Fassets\u002Fcfd-cases\u002Fflange\u002Fimages\u002F02.webp","\u002Fbrand\u002Fstudies\u002Fflange.svg",1448,1086,"Illustration DEEB-FLOW de la déformation totale d’un modèle de bride",13,"structural",[18,19,20,21,22],"CFD ou FEA","différence entre CFD et FEA","volumes finis et éléments finis","contraintes thermiques couplage CFD FEA","interaction fluide-structure","Illustration conceptuelle issue de la bibliothèque de cas DEEB-FLOW. Elle ne constitue pas un résultat de cet article.",[25,31,46,56,66,79,88,133,146,158],{"id":26,"title":27,"paragraphs":28},"a-common-question","Une question fréquente et une meilleure réponse",[29,30],"« Faut-il faire de la CFD ou des éléments finis sur cet échangeur ? » C’est l’une des questions qu’un ingénieur calcul entend le plus souvent, et la réponse honnête est en général : « les deux, dans un ordre précis, et l’ordre compte plus que le choix ».","CFD ou FEA, en bref : la CFD (mécanique des fluides numérique) prédit comment les fluides s’écoulent et transportent la chaleur ; la FEA (analyse par éléments finis) prédit le plus souvent comment les solides se déforment, vibrent et supportent des contraintes. La CFD répond aux questions d’écoulement, de perte de charge, de mélange et de coefficients d’échange ; la FEA répond aux questions de contraintes, de déplacements, de fatigue et de fréquences propres. Là où la chaleur ou un effort passe d’un fluide à une structure, la sortie d’un modèle devient l’entrée de l’autre, et l’erreur passe avec elle. Cet article explique pourquoi les deux ont développé des mathématiques différentes, où ils se rejoignent et, chiffres à l’appui, comment une incertitude sur un résultat CFD devient une incertitude sur une contrainte.",{"id":32,"title":33,"paragraphs":34},"the-wrong-comparison","Pourquoi « CFD ou FEA » compare les mauvaises choses",[35,36,45],"L’expression met en face une application et une méthode. Comme le dit Cadence dans l’une des meilleures explications d’éditeur, la FEA est une méthode pour construire un schéma numérique, tandis que la CFD est un domaine d’application des méthodes de calcul. Les éléments finis peuvent résoudre des problèmes de fluides, et les volumes finis des problèmes de solides. Les paires honnêtes sont celles-ci :",{"table":37},{"head":38,"rows":41},[39,40],"Par application","Par méthode",[42],[43,44],"Simulation des fluides (CFD) face à la simulation des structures","La méthode des volumes finis (MVF) face à la méthode des éléments finis (MEF)","En pratique, les deux paires se superposent : la plupart des codes CFD commerciaux sont en volumes finis, et la plupart des codes de structure en éléments finis. Les raisons sont mathématiques, et elles expliquent beaucoup de choses sur la façon dont chaque type de modèle se trompe.",{"id":47,"title":48,"paragraphs":49},"finite-volumes-for-fluids","Pourquoi les fluides se résolvent surtout en volumes finis",[50,51,54,55],"Un code en volumes finis écrit la loi de conservation pour chaque maille sous forme intégrale, avec des termes d’accumulation, de convection, de diffusion et de source :",{"formula":52},[53],"d\u002Fdt ∫_V ρφ dV + ∮_A ρφ (u·n) dA = ∮_A Γ ∇φ·n dA + ∫_V S dV","Tout ce qui quitte une maille par une face entre dans sa voisine par la même face, avec le même flux calculé. Sommez sur toutes les mailles et les flux internes s’annulent exactement. Masse, quantité de mouvement et énergie sont conservées globalement à l’arrondi près, quelle que soit la qualité du maillage : c’est une propriété qu’un ingénieur peut vérifier dans le bilan de masse de chaque calcul, et qu’un modèle de fluide doit avoir. Un calcul de perte de charge qui perd 2 % de la masse n’est pas un calcul de perte de charge.","La seconde raison est la convection. Le terme convectif transporte l’information dans un seul sens, vers l’aval. Les équations ne sont pas auto-adjointes, et une méthode numérique qui traite l’amont et l’aval de façon symétrique peine. En volumes finis, le remède est naturel : interpoler les valeurs aux faces avec un biais vers l’amont, avec des limiteurs qui gardent la solution bornée. C’est de là que vient le nombre de Péclet de maille, objet de l’exemple chiffré de notre article sur l’apprentissage de la CFD : le schéma centré oscille quand le Péclet de maille dépasse 2, et le schéma amont du premier ordre cache son erreur sous forme de diffusion supplémentaire.",{"id":57,"title":58,"paragraphs":59},"finite-elements-for-solids","Pourquoi les solides se résolvent surtout en éléments finis",[60,61,64,65],"L’élasticité linéaire a une autre structure. Sa forme faible, point de départ de tout code de structure en éléments finis, vaut pour toute fonction test admissible v :",{"formula":62},[63],"∫_V σ(u) : ε(v) dV = ∫_V f·v dV + ∮_A t·v dA","L’opérateur est symétrique et défini positif. Cela a trois conséquences sur lesquelles les calculateurs de structures s’appuient tous les jours. La matrice de rigidité est symétrique, ce qui rend les solveurs efficaces. La méthode de Galerkin, où fonctions test et fonctions d’essai viennent du même espace, donne la meilleure approximation possible en norme d’énergie pour un maillage donné : c’est le lemme de Céa. Et les éléments d’ordre élevé (quadratiques, courbes, coques et poutres) s’intègrent naturellement, si bien que les concentrations de contraintes aux congés et aux trous se résolvent avec un nombre modeste d’éléments.","La mécanique des solides n’a pas de convection : elle ne rencontre jamais le problème qui a rendu les volumes finis attrayants pour les fluides. Elle a un principe variationnel : elle hérite du problème pour lequel les éléments finis sont les meilleurs.",{"id":67,"title":68,"paragraphs":69},"where-they-meet","Là où chaque méthode rencontre le problème de l’autre",[70,71,72,75,76,77,78],"Poussez l’une des méthodes sur le terrain de l’autre, et les mêmes difficultés mathématiques apparaissent sous d’autres noms.","**Les éléments finis face à la convection.** Des éléments de Galerkin standard appliqués à un écoulement dominé par la convection oscillent exactement comme le schéma centré ; sur un maillage 1D uniforme avec des éléments linéaires, les deux schémas sont identiques. Le remède classique est la formulation de Petrov-Galerkin décentrée le long des lignes de courant (SUPG) de Brooks et Hughes (1982), qui modifie les fonctions test pour ajouter de la diffusion uniquement le long des lignes de courant. Pour la convection-diffusion 1D, le paramètre de stabilisation optimal vaut :",{"formula":73},[74],"τ = (h \u002F 2|u|) · [coth(Pe_h) − 1\u002FPe_h],   Pe_h = |u| h \u002F (2Γ)","Avec lui, la solution en éléments linéaires est exacte aux nœuds. La CFD en éléments finis fonctionne bien, et plusieurs codes multiphysiques commerciaux l’utilisent. Elle demande de choisir cette stabilisation avec soin, et comme le note l’article de Cadence, la bonne forme dépend du problème.","**L’incompressibilité, des deux côtés.** Les solides quasi incompressibles (le caoutchouc, au coefficient de Poisson proche de 0,5) font verrouiller les éléments en déplacement de bas ordre : ils deviennent artificiellement raides. Le remède est une formulation mixte avec la pression comme inconnue séparée, et les espaces de pression et de déplacement doivent satisfaire une condition de compatibilité, la condition inf-sup ou LBB. L’écoulement incompressible a exactement la même contrainte. En volumes finis colocalisés, elle se manifeste par des champs de pression en damier, et le remède est l’interpolation de Rhie-Chow qu’utilise tout code CFD commercial. Les mêmes mathématiques dans un autre vocabulaire, l’un des plus jolis faits de la mécanique numérique.","**Les structures minces.** Les éléments volumiques linéaires en flexion souffrent de verrouillage en cisaillement ; les éléments de coque et de poutre existent parce qu’une plaque mince maillée en briques est à la fois coûteuse et trop raide. Un ingénieur CFD qui maille un déflecteur mince rencontre le problème inverse : comment représenter une paroi sans épaisseur dans un maillage volumique.","La leçon pour un responsable est simple. Les deux types de code peuvent être poussés sur le terrain de l’autre. Ce faisant, on quitte des valeurs par défaut éprouvées pour des choix qui demandent un spécialiste.",{"id":80,"title":81,"paragraphs":82},"heat-transfer","Le transfert thermique : là où les deux se rejoignent",[83,84],"Un modèle de structure FEA résout très bien la conduction dans le solide. Ce qu’il ne peut pas calculer, c’est le coefficient d’échange h à la surface, parce qu’il dépend de l’écoulement. Il lui faut h en entrée. Deux questions décident si une corrélation ou une valeur de manuel suffit, ou si la CFD doit fournir h.",{"list":85},[86,87],"**L’écoulement est-il assez simple pour une corrélation ?** Écoulement établi dans un tube droit, plaque plane, cylindre en écoulement transversal : des corrélations existent et leur dispersion est connue. Jets impactants, écoulement décollé derrière une marche, convection naturelle dans une enceinte, mélange de courants à des températures différentes : le h local varie d’un facteur important sur la surface, et une valeur unique est une supposition. Quand les coefficients de convection sont incertains, ou que la convection naturelle compte, utilisez la CFD.","**Le solide est-il sensible à la valeur locale ?** Le nombre de Biot, Bi = hL\u002Fk, compare la résistance à la conduction dans le solide à la résistance à la convection à sa surface. Un petit nombre de Biot signifie que le solide est presque isotherme dans son épaisseur. Cela ne veut pas dire que la contrainte est faible, comme le montre l’exemple ci-dessous.",{"id":89,"title":90,"paragraphs":91},"worked-example","Exemple chiffré : comment une erreur CFD devient une erreur de contrainte",[92,93,94,102,103,104,107,108,131,132],"Une paroi en acier au carbone de 10 mm d’épaisseur sépare un gaz chaud d’une eau de refroidissement. Valeurs typiques pour l’acier : conductivité k = 45 W\u002F(m·K), module d’Young E = 200 GPa, coefficient de dilatation α = 12 × 10⁻⁶ 1\u002FK, coefficient de Poisson ν = 0,3. Le gaz est à 500 °C avec h_g = 100 W\u002F(m²·K) ; l’eau à 50 °C avec h_w = 2 000 W\u002F(m²·K). Toutes les valeurs sont illustratives.","Les résistances thermiques du film gazeux, de la paroi et du film d’eau sont en série, par mètre carré :",{"formula":95},[96,97,98,99,100,101],"1\u002Fh_g = 1\u002F100 = 0.010000 m²·K\u002FW","t\u002Fk = 0.010\u002F45 = 0.000222 m²·K\u002FW","1\u002Fh_w = 1\u002F2000 = 0.000500 m²·K\u002FW","R = 0.010722 m²·K\u002FW","q = (500 − 50) \u002F 0.010722 ≈ 42.0 kW\u002Fm²","ΔT = q · t\u002Fk ≈ 9.3 K","Les surfaces de la paroi sont à environ 80 °C et 71 °C. Le nombre de Biot côté gaz vaut 100 × 0,010 \u002F 45 ≈ 0,02 : la paroi est presque isotherme selon tous les critères du transfert thermique.","Passons à la contrainte. Pour une plaque avec une distribution de température linéaire dans l’épaisseur, libre de se dilater dans son plan mais empêchée de fléchir, la contrainte en surface vaut :",{"formula":105},[106],"σ = E α ΔT \u002F [2 (1 − ν)] = 200×10⁹ × 12×10⁻⁶ × 9.3 \u002F 1.4 ≈ 16 MPa","Seize mégapascals pour une différence de 9 K à travers une paroi « presque isotherme ». Supposons maintenant que le h côté gaz vienne d’une corrélation, et que l’écoulement réel, avec une recirculation près d’une entrée, donne 30 % de plus ou 30 % de moins :",{"table":109},{"head":110,"rows":115},[111,112,113,114],"h côté gaz, W\u002F(m²·K)","Flux de chaleur, kW\u002Fm²","ΔT de paroi, K","Contrainte thermique, MPa",[116,121,126],[117,118,119,120],"70 (30 % de moins)","30.0","6.7","11.4",[122,123,124,125],"100 (cas de base)","42.0","9.3","16.0",[127,128,129,130],"130 (30 % de plus)","53.5","11.9","20.4","Une variation de 30 % de h déplace la contrainte de 27 à 29 %. La raison se lit dans le budget des résistances : le film gazeux porte 93 % de la résistance totale, si bien que le flux de chaleur, et avec lui l’écart de température dans la paroi et la contrainte, évoluent presque un pour un avec h_g. Quelle que soit l’incertitude que votre CFD (ou votre corrélation) porte sur le coefficient d’échange, la contrainte en hérite presque entièrement.","Passez à un acier inoxydable austénitique, de valeurs typiques k ≈ 16 W\u002F(m·K), E ≈ 193 GPa et α ≈ 16 × 10⁻⁶ 1\u002FK, et le même flux de chaleur produit un écart d’environ 25 K dans la paroi et une contrainte d’environ 56 MPa. Le côté fluide change à peine ; la contrainte fait plus que tripler. C’est pourquoi aucun des deux ingénieurs ne peut travailler seul sur une pièce chaude : l’ingénieur CFD maîtrise q, l’ingénieur structure maîtrise ce que q produit.",{"id":134,"title":135,"paragraphs":136},"coupling","Quand faut-il les coupler, et comment ?",[137,138,139,140,141],"**Le couplage unidirectionnel** suffit pour la plupart des problèmes industriels. Lancez la CFD, transférez les températures (ou les pressions) sur le maillage de structure, lancez la FEA. Cela suppose que la déformation de la structure ne modifie pas l’écoulement, ce qui est vrai quand les déplacements sont petits devant les jeux et les couches limites qui gouvernent l’écoulement. Pour une bride, une plaque tubulaire d’échangeur ou un carter de turbine sous charge permanente, c’est presque toujours le cas.","**Le couplage bidirectionnel** (interaction fluide-structure, IFS) est nécessaire quand la déformation modifie sensiblement l’écoulement : flottement, membranes souples, clapets et joints qui s’ouvrent sous pression, vaisseaux sanguins. Il est coûteux, et les schémas partitionnés, où les solveurs fluide et structure travaillent à tour de rôle, ont un piège. Quand la structure est légère devant le fluide qu’elle déplace (un panneau mince dans l’eau, une artère dans le sang), l’inertie du fluide agit comme une masse ajoutée qui déstabilise l’échange, un mécanisme analysé par Causin, Gerbeau et Nobile (2005). L’acier dans l’air le rencontre rarement. Les structures minces dans les liquides le rencontrent souvent.","**Le transfert thermique conjugué** se situe entre les deux. Un code CFD qui inclut les régions solides résout ensemble la température du fluide et celle du solide, si bien que h n’a jamais à être deviné ; le champ de température est ensuite transféré à la FEA pour les contraintes. Pour l’exemple ci-dessus, c’est la bonne démarche chaque fois que le côté gaz n’est pas un cas simple de corrélation.","Le transfert de données est l’endroit où les bons modèles se trompent. Les deux maillages ne coïncident presque jamais : le maillage CFD est fin dans la couche limite du fluide, le maillage FEA est fin aux congés. Quand vous transférez des données de l’un à l’autre :",{"list":142},[143,144,145],"**Les champs comme la température demandent un transfert cohérent :** interpolez la valeur à chaque nœud cible.","**Les charges comme les efforts et les flux de chaleur demandent un transfert conservatif :** le total transféré doit égaler le total calculé. Vérifiez-le. Une charge de pression transférée dont la résultante diffère de 5 % de l’effort CFD ajoute 5 % d’erreur à tout ce qui suit.","**En instationnaire, les instants doivent correspondre.** La contrainte thermique maximale au démarrage survient souvent tôt, quand la surface chauffe plus vite que le cœur. Transférer la température finale stabilisée la manque.",{"id":147,"title":148,"paragraphs":149},"common-mistakes","Erreurs fréquentes, telles qu’elles apparaissent à l’écran",[150],{"list":151},[152,153,154,155,156,157],"**Un modèle de structure avec un seul coefficient d’échange sur toutes les surfaces.** Un champ de température lisse, une contrainte plausible, et une erreur de la taille de la dispersion sur h.","**« La paroi est presque isotherme, la contrainte thermique est donc négligeable. »** L’exemple donne 16 MPa pour 9 K.","**Un champ de température transféré avec une couture.** L’interpolation a échoué là où les maillages divergent, souvent à un congé, c’est-à-dire là où se trouve la contrainte.","**De l’IFS bidirectionnelle par défaut.** Des jours de calcul pour une déformation de quelques microns qui ne change rien à l’écoulement.","**Un calcul IFS partitionné qui diverge pour une structure légère dans un liquide.** C’est l’instabilité de masse ajoutée. Raffiner le maillage ne la guérira pas.","**CFD et FEA livrées par des équipes différentes avec des hypothèses différentes.** Des températures de matériau différentes, des points de fonctionnement différents, et une contrainte qui n’appartient à aucun des deux modèles.",{"id":159,"title":160,"paragraphs":161},"what-i-would-do","Ce que je ferais",[162,163,164],"Pour une pièce où la chaleur passe d’un fluide à une structure, je commencerais par le budget des résistances sur papier, comme dans l’exemple. Cela prend dix minutes et dit quel côté gouverne le résultat. Si le côté fluide domine et que l’écoulement n’est pas un cas de manuel, le coefficient d’échange doit venir de la CFD, de préférence d’un modèle de transfert thermique conjugué, et le résultat de structure ne vaut que ce que vaut ce coefficient. Si le côté solide domine, une corrélation pour h peut suffire, et l’effort doit aller à la FEA.","J’utiliserais le couplage unidirectionnel sauf raison claire de penser que la déformation modifie l’écoulement, et je vérifierais chaque transfert en comparant les totaux : chaleur entrante face à chaleur transférée, effort calculé face à effort appliqué.","Et je garderais les deux modèles dans les mêmes mains, ou au moins dans la même pièce. Nos études combinent la CFD sous Ansys Fluent et l’analyse de structure sous Ansys Mechanical précisément pour cette raison, et nos deux cours, Modélisation CFD des écoulements de liquides et de gaz et Static Structural, sont conçus pour être suivis ensemble par les ingénieurs à qui l’on posera la question par laquelle cet article a commencé.",[166,169,172],{"href":167,"label":168},"https:\u002F\u002Fresources.system-analysis.cadence.com\u002Fblog\u002Fmsa2021-fea-vs-cfd-the-differences-and-applications-of-simulation-tools","Cadence: FEA vs CFD, the differences and applications of simulation tools",{"href":170,"label":171},"https:\u002F\u002Fwww.siemens.com\u002Fen-us\u002Fcontent\u002Ffea-vs-cfd\u002F","Siemens: FEA vs CFD",{"href":173,"label":174},"https:\u002F\u002Fdoi.org\u002F10.1016\u002Fj.cma.2004.12.005","Causin, Gerbeau and Nobile (2005), Computer Methods in Applied Mechanics and Engineering: added-mass effect in partitioned fluid-structure algorithms",1792989865051]