Comment utiliser une méthode hybride deAnalyse CFDpour l'optimisation des produits et l'amélioration des performances tout en réduisant l'effort de calcul et le temps nécessaire pour obtenir les résultats
L'IAO est très efficace pour l'optimisation des produits. Dans cet article technique, nous décrivons comment une méthode hybride d'analyse de dynamique des fluides numérique a été utilisée pour améliorer la conception d'un dispositif médical afin d'améliorer sa capacité à maintenir un microclimat stable autour du patient, mais aussi pour réduire les efforts de calcul et le temps nécessaires à l'obtention de ces résultats. Un microclimat nécessite le maintien de valeurs spécifiques de température, d'humidité et de vitesse de l'air autour du patient.
Dans les appareils ouverts, cette chaleur est généralement générée et entretenue au moyen de trois jets d'air chaud humidifié positionnés l'un de chaque côté et l'autre au pied du lit.
Les paramètres de conception de cette configuration forment un problème d'optimisation où la somme des écarts relatifs de température et d'humidité doit être minimisée au niveau de la région médiane du patient, tandis que la vitesse de l'air doit être limitée à moins de 0,35 m/s et que la perte de chaleur à travers le plafond vers la pièce doit être contenue.
Dans ce Conseil en CFD Dans l'étude, la méthode hybride a été adoptée en utilisant ANSYS CFX pour calculer le champ d'écoulement et modeFRONTIER pour automatiser le processus de conception et pour conduire la chaîne de calcul vers l'identification d'une solution optimale.
Ce travail décrit une technique avancée pour améliorer les résultats d'une optimisation multi-objectifs en utilisant des technologies avancées, telles que les cartes auto-organisatrices (SOM) et les modèles de surface de réponse (RSM) de modeFRONTIER, afin d'orienter les analyses dans la direction souhaitée avec peu de calculs supplémentaires. Dans de nombreux défis d'ingénierie réels, tels que celui décrit ici, il existe un grand nombre de paramètres d'entrée couvrant un espace de grande dimension, et les évaluations de conception sont effectuées avec des exécutions de modèles numériques coûteuses en calcul. Par conséquent, il est nécessaire de réduire autant que possible les variantes de conception nécessaires à l'optimiseur pour trouver les améliorations souhaitées. Cet article propose une procédure novatrice appliquée à l'optimisation du champ d'écoulement d'un dispositif médical.
L'objectif de cette étude est d'optimiser la conception d'un appareil médical permettant de maintenir un microclimat stable autour du patient. Dans les dispositifs ouverts, tels que ceux abordés dans cet article, le microclimat peut être obtenu grâce à une disposition spéciale de jets d'air chauffés et humidifiés : la figure 1 illustre une configuration schématique, avec la surface du lit équipée de trois buses formant les jets mentionnés ci-dessus, une de chaque côté et une au pied de l'appareil. La récupération de chaleur et d'humidité s'effectue par recirculation, alimentée par aspiration côté tête de l'appareil (non illustrée ici) : cette récupération doit être maximisée afin de réduire les pertes de chaleur. Le microclimat comprend des valeurs spécifiques de température, d'humidité et de vitesse autour du patient : comme souvent dans la littérature, la position de référence est supposée être le centre du corps du patient, désigné par le point A (voir figure 1). Ces sorties sont principalement influencées par le champ d'écoulement induit par les jets d'air des buses. Par conséquent, les paramètres de conception des jets d'air, tels que la distribution de vitesse, l'angle de vitesse, la température et l'humidité, constituent un problème d'optimisation.
Le La somme des écarts relatifs de température (T) et d'humidité (h%) au point A a été définie comme un objectif à minimiser : cible = (|T_A − T_target| / T_target) + (|h%_A − h%_target| / h%_target) De plus, elle doit être inférieure à 0,025 (2,5 %). La vitesse de l'air au point A a été contrainte à être inférieure à 0,35 m/s. Le deuxième objectif est lié à l'indice de perte de chaleur : la différence entre l'enthalpie des jets d'air d'entrée et l'air quittant le domaine du modèle par le plafond vers la pièce doit être minimisée. Au total, dix variables d'entrée ont été utilisées pour paramétrer le flux d'air chaud. En particulier, un premier ensemble de variables est lié à l'amplitude de la vitesse des jets d'air de taille et de pied, un autre ensemble à leur orientation et à leur profil, tandis que les deux derniers paramètres étaient leur température et leur humidité.
La dynamique des fluides numérique (CFD) est une technique efficace pour modéliser le champ d'écoulement et ses dépendances à partir des paramètres de conception. Le domaine de calcul comprend l'espace autour du patient et un cube formant l'environnement autour de l'appareil (la pièce). À la sortie des buses, des conditions limites d'entrée (BC) avec des profils de vitesse et des distributions de température et d'humidité appropriés ont été spécifiées et liées aux valeurs des variables d'entrée. Une BC de sortie a été spécifiée au niveau de la grille d'aspiration, tandis que les surfaces solides et les limites latérales de l'appareil ont été modélisées comme des parois adiabatiques. Le plafond de la pièce a été spécifié comme un BC d'ouverture (pression = 1 atm et température = 20 °C) afin de tenir compte des pertes de chaleur et d'air dues à la flottabilité.
Un maillage tétraédrique a été généré pour le domaine de calcul, et des calculs RANS (Reynolds Averaged Navier-Stokes) ont été effectués. Dans cette étude, ANSYS CFX est utilisé pour calculer le champ d'écoulement, tandis que modeFRONTIER est utilisé pour automatiser le processus de conception et orienter la chaîne de calcul vers une solution optimale. Le processus a été mis en place dans l'environnement modeFRONTIER et est représenté par le diagramme de flux de travail de la figure 2.
Dans un premier temps, une optimisation directe avec un algorithme génétique multi-objectif (NSGA-II) a permis au concepteur d'explorer un large éventail de solutions possibles. Cependant, les résultats n'ont pas été pleinement satisfaisants : le confort (objectif cible) a été suroptimisé par rapport aux pertes thermiques. Pour des raisons de temps, seulement quelques calculs supplémentaires étaient encore disponibles ; une stratégie avancée a donc été mise en place. L'approche proposée vise à atteindre rapidement une zone souhaitée de l'ensemble de Pareto, non cartographiée par l'optimisation initiale. Une application des techniques SOM a d'abord été utilisée pour détecter les corrélations locales et, par conséquent, circonscrire la portée de la recherche dans l'espace des variables d'entrée, réduisant ainsi le problème. Après la définition de nouvelles bornes pour les variables d'entrée, un nouveau plan d'expérience (DOE) a été mis en place et évalué par des simulations CFD, puis les nouvelles données ont été interpolées par des techniques de modélisation de surface de réponse (RSM). Finalement, une exécution quasi instantanée de l'algorithme génétique a été appliquée à ces méta-modèles, et les prévisions les plus prometteuses ont été validées au moyen de nouvelles simulations CFD, livrant le résultat final au concepteur.
Outre les améliorations apportées à la conception initiale, ce type de méthodes hybrides garantit une réduction des efforts de calcul nécessaires à l'obtention de tels résultats, par rapport à une optimisation directe complète. En particulier, l'utilisation des techniques SOM permet au concepteur de sélectionner la zone la plus prometteuse pour ses objectifs de recherche. Cette technologie offre deux avantages majeurs à ce type d'applications. Premièrement, elle est encore plus rapide et pratique par rapport aux autres approches, tout en augmentant la dimension de l'espace des variables d'entrée. Deuxièmement, le processus de focalisation sur les bornes exploite non seulement les solutions calculées, mais aussi l'estimation de la densité de probabilité de l'espace des solutions.
Initialement, les bornes des 10 variables d'entrée indépendantes étaient assez larges, afin de cartographier l'espace des solutions possibles de manière exhaustive. Avec cette configuration, un algorithme multi-objectif NSGA-II a été utilisé pour piloter le code CFD dans l'exploration de nouveaux plans. La population initiale de 20 plans a été obtenue en ajoutant une séquence aléatoire de plans d'expériences (DOE) au plan original. L'algorithme génétique a été exécuté sur 15 générations, avec une probabilité de croisement de 0,9 et une probabilité de mutation de 1,0. Toutes les contraintes ont été considérées comme des fonctions de pénalité pour les objectifs assignés. La figure 3 montre l'évolution de la valeur de l'objectif cible au cours de l'optimisation. La figure 4 représente le nuage de points des configurations explorées dans les deux plans objectifs (cible et perte de chaleur). Dans cette première phase de l'optimisation, l'algorithme a pu sortir de la zone irréalisable d'où il était parti (la population initiale ne pouvait pas atteindre la contrainte assignée sur la cible).
La figure 4 met également en évidence que l'algorithme s'est concentré sur deux régions principales de solutions réalisables et a eu tendance à suroptimiser les valeurs d'humidité et de température, tandis que les pertes de chaleur n'étaient pas suffisamment améliorées. En fait, le nuage de points réalisables au bas de la figure 4 présente encore une marge appréciable par rapport à la limite assignée à l'écart par rapport à la performance idéale de climatisation (cible inférieure à 0,025). Puisque la cible et les pertes de chaleur semblent être des objectifs opposés, on peut s'attendre à trouver des configurations capables d'améliorer encore les pertes de chaleur en respectant plus étroitement la contrainte assignée à la cible. L'idée est de les rechercher de la manière la plus efficace, en exploitant les données déjà disponibles après l'étape d'optimisation avec NSGA-II, ainsi que quelques autres simulations CFD bien ciblées. À cette fin, une optimisation virtuelle appliquant une stratégie hybride SOM-RSM est proposée ici. L'avantage de cette approche par rapport à une optimisation directe plus longue (outre le nombre réduit de simulations supplémentaires nécessaires pour trouver les solutions souhaitées) est une compréhension plus approfondie des corrélations locales entre les paramètres et les objectifs d'optimisation. Ça devient encore plus clair quand on augmente le nombre de variables et de réponses impliquées.
La carte auto-organisée (SOM) est un algorithme qui projette des données de grande dimension sur une carte bidimensionnelle. Cette projection préserve la topologie des données afin que les éléments de données similaires soient mappés à des emplacements proches sur la carte.
Dans l'approche proposée, le SOM a été utilisé pour analyser les corrélations locales entre les deux objectifs (cible et indice de perte de chaleur) en exploitant les résultats de l'optimisation NSGA-II précédente. En effet, la détection de telles corrélations locales permet d'identifier une plage de paramètres d'entrée réduite incluant les conceptions les plus prometteuses pour les valeurs de sortie souhaitées. De cette manière, une campagne de recherche ultérieure pourrait être mise en place dans cet espace réduit de paramètres d'entrée.
La figure 5 présente les graphiques des composants relatifs au SOM construit. La grille de gauche rassemble l'ensemble des graphiques des composants du SOM, y compris les paramètres de conception. À droite, les graphiques des composants de la cible et de l'indice de perte de chaleur sont détaillés. Dans ces graphiques conventionnels, chaque hexagone est une « boîte » contenant un groupe de conceptions similaires, préalablement calculées par le modèle CFD ou interpolées par le réseau de neurones non supervisé du SOM. Ces conceptions sont « semblables » par rapport au vecteur complet de paramètres et de valeurs de sortie, qui les définit dans l'espace à n dimensions. Les hexagones adjacents désignent des groupes similaires. Un hexagone correspondant dans deux cartes de paramètres différentes (une d'entrée et une de sortie) représente la même « boîte », incluant donc le même ensemble de conceptions. La couleur de l'hexagone, à l'inverse, dépend de la carte elle-même : elle représente la valeur moyenne de ce paramètre au sein du groupe de conceptions sélectionné inclus dans la « boîte ». Comparer le modèle de couleur d'un ensemble correspondant d'hexagones dans deux cartes différentes signifie rechercher un ensemble de motifs avec une certaine tendance commune de paramètres (corrélation locale entre les paramètres).
Par conséquent, des corrélations locales entre plusieurs paramètres peuvent être détectées en comparant leurs cartes SOM et en recherchant les zones où les hexagones présentent des schémas de couleurs similaires. L'objectif étant de vérifier l'existence d'une corrélation locale entre la perte de chaleur et une valeur cible légèrement inférieure à 0,025, leurs cartes SOM sont comparées dans la figure 5 (à droite). Une corrélation locale des deux sorties, où la cible est légèrement inférieure à 0,025 (couleur bleu pâle), existe dans les unités SOM soulignées d'un trait plein dans la figure 5 en haut à droite, comme en témoigne le schéma de couleurs uniforme similaire. À titre d'exemple, aucune corrélation locale ne peut être évaluée pour les groupes de conception dans le carré en pointillés (figure 5 en bas à droite) : les schémas de couleurs correspondants de la cible et de la perte de chaleur ne sont pas similaires.
Après avoir sélectionné les hexagones inclus dans le cadre en trait plein à droite de la figure 5, l'espace de conception a été redéfini en limitant les paramètres de conception à la plage définie par les conceptions incluses dans ces unités SOM (représentée par la plage de couleurs des hexagones correspondants à gauche de la figure 5, ceux encadrés en vert). Chaque plage de variables d'entrée a été sensiblement réduite, dans la plupart des cas de plus de la moitié, avec une réduction moyenne de 43 %. Cela signifie que la taille du domaine de recherche est fortement réduite et concentrée uniquement sur les zones les plus prometteuses par rapport aux exigences du concepteur. Cela mène à une deuxième étape d'optimisation beaucoup plus rapide.
Une séquence de 16 configurations du DOE factoriel réduit a été mise en place afin d'ajouter uniformément peu de nouveaux points dans cet espace restreint de paramètres d'entrée. Tous les plans précédemment simulés appartenant aux unités SOM sélectionnées, et donc compris dans la plage restreinte, ont été ajoutés à la base de données de résultats du DOE factoriel réduit. Un RSM, à savoir un ensemble de fonctions de base radiales (RBF), a ensuite été entraîné avec ces données afin de modéliser chacun des objectifs du problème. Ces RSM ont ensuite été utilisées pour interpoler de nouvelles configurations, coupant ainsi le lien entre l'optimiseur et les calculs CFD coûteux en ressources.
En fait, sur de tels RSM, une nouvelle optimisation NSGA-II approfondie a été réalisée afin de rechercher le(s) optimum(s) présentant des valeurs de pertes thermiques réduites, tout en respectant la contrainte cible. Ainsi, environ 1 800 points ont été estimés en quelques secondes. Les résultats ont ensuite été filtrés selon les critères de dominance, et les modèles candidats les plus prometteurs ont été sélectionnés pour être validés par de nouvelles simulations CFD. Ce choix a été fait en tenant compte non seulement des valeurs estimées des objectifs, mais aussi de la distance de ces points par rapport aux points de consigne d'apprentissage dans l'espace des paramètres à 10 dimensions. Trois modèles ont été sélectionnés puis validés par de nouvelles simulations CFD.
L'une des conceptions validées s'est révélée très prometteuse, car elle permet une réduction supplémentaire des pertes thermiques par rapport à la meilleure conception issue de la première phase d'optimisation du NSGA-II : un gain supplémentaire de 4 %. Comme prévu, ce résultat a été obtenu grâce à une solution qui reste proche de la valeur limite de température et d'humidité du patient, tout en conservant une marge importante par rapport à la valeur de contrainte assignée de 0,025.
En d'autres mots, l'approche proposée a bien fonctionné pour atteindre la zone souhaitée de l'ensemble de Pareto, non cartographiée par l'optimisation initiale de la première phase de NSGA-II. Il convient de souligner que l'amélioration a été obtenue avec seulement 19 calculs CFD supplémentaires (DOE de 16 pour l'espace de paramètres de conception restreint, plus 3 validations), et donc en moins d'une génération supplémentaire de l'optimisation NSGA-II initiale.
La figure 6 compare les champs de température dans le plan de symétrie vertical pour la conception originale (à gauche), la conception améliorée NSGA-II (au centre) et la conception finale (à droite). La dispersion thermique réduite est due à un panache d'air chaud plus petit et à une distribution de température plus régulière (généralement plus basse) à l'intérieur du dispositif. De plus, la conception proposée est également assez originale, par rapport au résultat de la première optimisation directe, en ce qui concerne sa localisation dans l'espace des paramètres d'entrée. Cela signifie que l'approche proposée a trouvé un nouvel « îlot » optimal dans l'espace des paramètres de conception avec un « coût » de seulement 19 calculs CFD supplémentaires après les 300 calculs CFD initiaux de la première phase d'optimisation NSGA-II.
À propos de GoEngineer
GoEngineer offre des logiciels, des technologies et une expertise qui permettent aux entreprises de libérer l'innovation en matière de conception et de livrer de meilleurs produits plus rapidement. Avec plus de 40 ans d'expérience et des dizaines de milliers de clients dans les secteurs de la haute technologie, des dispositifs médicaux, de la conception de machines, de l'énergie et d'autres industries, GoEngineer fournit des solutions de conception de premier plan avec SOLIDWORKS CAD, l'impression 3D Stratasys, la numérisation 3D Creaform et Artec, CAMWorks, PLM, et plus.
Recevez notre vaste gamme de ressources techniques directement dans votre boîte de réception.
Désabonnez-vous à tout moment.