Spelling suggestions: "subject:"numériques"" "subject:"mumériques""
291 |
Modélisation instationnaire de l'aérodynamique externe automobile / Unsteady computation of external aerodynamics flow in automotive industryDelassaux, François 20 December 2018 (has links)
La thèse a pour but de développer une méthodologie de calcul instationnaire permettant une étude qualitative et quantitative de l’aérodynamique externe d’une automobile. La première partie de l’étude est consacrée au développement de la méthodologie numérique sur les différents corps d’Ahmed à 25°, géométries simplifiées d’une automobile réelle, afin de valider les choix stratégiques de maillages et de méthodes numériques ainsi que s’assurer de la bonne résolution de l’écoulement. Les résultats numériques sont comparés aux données expérimentales obtenues au cours d’essais réalisés à la soufflerie La Ferté Vidame lors de travaux de thèse précédents. A l’issue de ce travail, la méthode hybride Delayed Detached Eddy Simulation Shear-Stress Transport (DDES SST) est sélectionnée pour la suite de l’étude au vu des meilleures performances obtenues (torseur aérodynamique, coefficient de pression, topologie d’écoulement) par rapport aux méthodes Reynolds Averaged Navier-Stokes (RANS), Scale Adaptive Simulation (SAS) et Stress-Blended Eddy Simulation (SBES). La seconde phase de l’étude consiste à adapter la méthodologie précédemment développée sur un véhicule réel, la Peugeot 308 SW. Au préalable, une base de données expérimentales conséquente a été réalisée au sein du Groupement d’Intérêt Economique Souffleries Aéroacoustiques Automobiles (GIE S2A) au cours de ces travaux. La géométrie est tout d’abord simplifiée afin de faciliter la mise en place de la méthodologie numérique : entrées d’air fermées, soubassement lissé, roues remplacées par des carénages. Les résultats obtenus sont encourageants et démontrent globalement la supériorité de la DDES par rapport aux méthodes RANS classiques. La topologie d’écoulement est mieux prédite (soubassement et sillage), même si la prédiction du coefficient de portance reste une difficulté majeure pour ce type de méthode hybride. / The main goal of this PhD is to develop an unsteady numerical method to study the external aerodynamic flow around real vehicles. The first part of the study focuses on the flow around simplified geometries, such as 25° Ahmed bodies (with sharp and rounded edges on the back of the body), in order to determine the optimal turbulence model, mesh setup and numerical parameters. Computational Fluid Dynamics (CFD) results are compared to experimental data reported in literature conducted in the La Ferté Vidame wind tunnel. Based on this study, Shear-Stress Transport Delayed Detached Eddy Simulation (SST DDES) demonstrates superiority over Reynolds Averaged Navier-Stokes (RANS), Scale Adaptive Simulation (SAS) and Stress Blended Eddy Simulation (SBES) turbulence models, regarding both drag and lift coefficients predictions, and flow topology.Secondly, the numerical procedure is adapted for a real vehicle, the Peugeot 308 SW estate car. A substantial experimental campaign was carried out in the Groupement d’Intérêt Economique Souffleries Aéroacoustiques Automobiles (GIE S2A) wind tunnel to provide data against which the numerical results are compared. Given the geometric complexity of a real vehicle, the car is simplified for this study as follows: the front air inlets are closed, the underbody is smoothed with additional panels and the wheels are replaced by fairings. DDES computations show encouraging results. A significant improvement of the flow topology is obtained with DDES compared to RANS models. However, the prediction of the lift coefficient remains a major difficulty with these hybrid methods.
|
292 |
Applications des méthodes multigrilles à l'assimilation de données en géophysique / Multigrid methods applied to data assimilation for geophysics modelsNeveu, Emilie 31 March 2011 (has links)
Depuis ces trente dernières années, les systèmes d'observation de la Terre et les modèles numériques se sont perfectionnés et complexifiés pour nous fournir toujours plus de données, réelles et numériques. Ces données, de nature très diverse, forment maintenant un ensemble conséquent d'informations précises mais hétérogènes sur les structures et la dynamique des fluides géophysiques. Dans les années 1980, des méthodes d'optimisation, capables de combiner les informations entre elles, ont permis d'estimer les paramètres des modèles numériques et d'obtenir une meilleure prévision des courants marins et atmosphériques. Ces méthodes puissantes, appelées assimilation variationnelle de données, peinent à tirer profit de la toujours plus grande complexité des informations de par le manque de puissance de calcul disponible. L'approche, que nous développons, s'intéresse à l'utilisation des méthodes multigrilles, jusque là réservées à la résolution de systèmes d'équations différentielles, pour résoudre l'assimilation haute résolution de données. Les méthodes multigrilles sont des méthodes de résolution itératives, améliorées par des corrections calculées sur des grilles de plus basses résolutions. Nous commençons par étudier dans le cas d'un modèle linéaire la robustesse de l'approche multigrille et en particulier l'effet de la correction par grille grossière. Nous dérivons ensuite les algorithmes multigrilles dans le cadre non linéaire. Les deux types d'algorithmes étudiés reposent d'une part sur la méthode de Gauss Newton multigrille et d'autre part sur une méthode sans linéarisation globale : le Full Approximation Scheme (FAS). Ceux-ci sont appliqués au problème de l'assimilation variationnelle de données dans le cadre d'une équation de Burgers 1D puis d'un modèle Shallow-water 2D. Leur comportement est analysé et comparé aux méthodes plus traditionnelles de type incrémentale ou multi-incrémentale. / For these last thirty years, earth observation and numerical models improved greatly and provide now a huge amount of accurate, yet heterogeneous, information on geophysics fluids dynamics and structures. Optimization methods from the eighties called variational data assimilation are capable of merging information from different sources. They have been used to estimate the parameters of numerical models and better forecast oceanic and atmospheric flows. Unfortunately, these powerful methods have trouble making benefit of always more complex information, suffering from the lack of available powerful calculators. The approach developed here, focuses on the use of multigrid methods, that are commonly used in the context of differential equations systems, to solve high resolution data assimilation. Multigrid methods are iterative methods improved by the use of feedback corrections evaluated on coarse resolution. First in the case of linear assimilation, we study the robustness of multigrid approach and the efficiency of the coarse grid correction step. We then apply the multigrid algorithms on a non linear 1-D Burgers equation and on a 2-D Shallow-Water model. We study two types of algorithms, the Gauss Newton Multigrid, which lays on global linearization, and the Full Approximation Scheme. Their behavior is compared to more traditional approaches as incremental and multi-incremental ones.
|
293 |
L’apport des moyens numériques aux savoir-faire artisanaux en architecture : une analyse selon les modèles de la tectonique numérique d’OxmanMnajja, Houssem Eddine 08 1900 (has links)
Les moyens numériques ont transformé les pratiques en architecture notamment en facilitant la continuité des flux d’information entre la conception et la fabrication. À travers ce mémoire, nous explorons les contributions possibles de ces moyens dans un contexte de pénurie d’artisans et de perte de savoir-faire artisanaux en architecture.
Nous nous sommes appuyés sur les modèles de la tectonique informée pour la conception à base des matériaux, documentés par R. Oxman (2014), qui regroupent un ensemble d’approches et d’outils numériques, afin d’étudier, dans un premier temps, les apports et limites de ces derniers à partir de plusieurs cas expérimentaux; c’est-à-dire expérimental. Ainsi, nous avons pu identifier certains savoir-faire artisanaux dont le déclin est possiblement compensable par les moyens numériques. Dans un deuxième temps, nous avons décelé les contributions des moyens numériques dans un projet à partir de l’analyse des exemples de la construction de la Sagrada Familia et la restauration du parlement canadien à Ottawa.
Les principaux résultats de notre recherche sont :
• Nous avons distingué que la modélisation paramétrique est plus adaptée au processus de conception tandis que la modélisation à base de relevés est plus appropriée pour la production des répliques.
• Les outils numériques de fabrication permettent d’employer plusieurs types de matériaux. Elles assurent une rapidité et une précision dans la réalisation des artéfacts. Cependant, l’intervention des artisans demeure incontournable dans la majorité des cas.
Les moyens numériques ont un grand potentiel pour compenser la perte des savoir-faire artisanaux et pour développer de nouveaux processus créatifs. Toutefois, à ce jour, ils se limitent au rôle d’assistants et devraient être perfectionnés. / Digital tools have transformed the architectural practices due to the information flow continuity between design and manufacturing. Through this these, we explore their possible contributions in the context of artisans’ shortage and craft skills’ losses in architecture.
We relied on informed tectonic models of materials based design, documented by R. Oxman (2014), which combines a variety of digital approaches and tools to study their contributions and limits based on multiple case studies. Thus, we identified some craft skills whose decline is possibly compensated by digital tools. Then, we analyzed the digital tools’ contributions in a project based on the Sagrada Familia and Canadian parliament examples.
The main results of our research are:
• Digital in architecture is based on digital models. The parametric approach is more suited to design process while models based on lasergrammetry and photogrammetry surveys are more appropriate for replica production.
• Digital manufacturing tools allow the use of several materials types. Similarly, they provide speed and accuracy. However, the artisans’ participation remains unavoidable in most cases.
Digital media have great potential to overcome the craft skills’ loss and even exploit them in new creative processes. However, to date, they are limited to the role of assistants and should be perfected.
|
294 |
Tracking and modelling small motions at grain scale in granular materials under compression by X-Ray microtomography and discrete simulations / Matériaux granulaires en compression quasi-statique : étude des petites déformations par microtomographie X et par simulation numérique discrèteKhalili, Mohamed Hassan 03 November 2016 (has links)
Le travail réalisé durant cette thèse a été motivé par l'étude des mécanismes microscopiques à l'origine du fluage dans les matériaux granulaires.%En particulier, on cherche à explorer des techniques expérimentales et numériques pour l'étude d'un tel phénomène.Dans une première partie, on cherche à mesurer les déplacements des grains dans un matériau granulaire par observations en micro-tomographie X. Une telle identification ne peut être efficacement réalisée pour des phénomènes rapides avec les méthodes classiques de corrélation d'images numériques. Une nouvelle méthode nommée emph{corrélation discrète des projections numériques} qui contourne cette difficulté est développée dans cette thèse. Cette méthode, basée sur la corrélation des projections de tomographie, permet de mesurer les déplacements avec un nombre réduit de projections (100 fois moins que les méthodes classiques), ce qui diminue énormément le temps d'acquisition nécessaire pour la mesure. La méthode, appliquée à des données expérimentales, donne une précision comparable à celles des méthodes classiques tandis que le temps d'acquisition nécessaire est réduit à quelques minutes. Une étude portant sur l'analyse des sources d'erreurs affectant la précision des résultats est également présentée.Le but de la deuxième partie est de réaliser des simulations numériques pour fournir une caractérisation de l'essai oedométrique. Différents assemblages de billes de verre légèrement poly-disperses interagissant à travers des contacts élastiques de Hertz-Mindlin et frottement de Coulomb ont été utilisés. Ces simulations ont permis d'étudier l'évolution de certains paramètres structuraux du matériau modèle, préparant ainsi le terrain pour de futures études sur le fluage. Il a été particulièrement souligné que les contacts élastiques utilisés dans ces simulations ne reproduisent pas l'irréversibilité des déformations observée dans les expériences sur des sables. Cependant, l'irréversibilité est bien visible sur le nombre de coordination et l'anisotropie. Alors que les paramètres élastiques peuvent exprimer la réponse pour des petits incréments de déformations, la compression oedometrique est belle et bien anélastique, principalement à cause de la mobilisation du frottement. Le rapport entre les contraintes horizontales et verticales (coefficient du sol au repos) n'est particulièrement constant que lorsque l'anisotropie de structure est instaurée dans l'état initial de l'assemblage. Il est par ailleurs relié à l'anisotropie interne de la structure par une formule simple. Finalement, les coefficients du tenseur élastique dépendent principalement du nombre de coordination et son anisotropie est plus liée à l'anisotropie des contacts qu'à celle des forces / The present work is motivated by the study of creep in granular materials at the microscopic scale.The first part of this thesis deals with displacement measurements by microtomography. Classical digital image correlation fails to catch time-dependent (possibly fast) phenomena such as short-term creep. A new method named emph{Discrete Digital Projection Correlation} is developed to overcome this limitation. This method requires very few projections (about 100 times less than classical methods) of the deformed state to perform the correlation and retrieve grain displacements. Therefore, the acquisition time is remarkably reduced, which allows to study time-dependent phenomena.The method is tested on experimental data. While its accuracy compares favorably to that of conventional methods, it only requires acquisition times of a few minutes. The origins of measurement errors are tracked by numerical means, on simulated grain displacements and rotations.The second part is a numerical simulation study, by the Discrete Element Method (DEM), of oedometric compression in model granular materials, carried out with a simple model material: assemblies of slightly polydisperse spherical beads interacting by Hertz-Mindlin contact elasticity and Coulomb friction. A wide variety of initial states are subject to compression, differing in density, coordination number and fabric anisotropy. Despite apparently almost reversible strains, oedometric compression proves an essentially anelastic and irreversible process,due to friction, with important internal state changes affecting coordination number and anisotropy. Elastic moduli only describe the response to very small stress increments about well equilibrated configurations. The ratio of horizontal stress to vertical stress (or coefficient of earth pressure at rest, commonly investigated in soil mechanics) only remains constant for initially anisotropic assemblies. A simple formula relates it to force and fabric anisotropy parameters, while elastic moduli are mainly sensitive to the latter. Further studies of contact network instabilities and rearrangements should pave the way to numerical investigations of creep behavior
|
295 |
Enabling pervasive applications by understanding individual and community behaviors / Nouvelles applications pervasives par la modélisation des comportements individuels et communautairesSun, Lin 12 December 2012 (has links)
Les empreintes digitales recueillies par détection systèmes offrent de nouvelles façons de percevoir les comportements d'un individu. En outre, de grandes collections d'empreintes numériques des communautés apportent de nouvelles compréhensions des comportements humains. La perception des comportements humains à partir des empreintes digitales de détection permet de construire des nouvelles applications sur les systèmes de détection. D’après les empreintes digitales recueillies avec l'accéléromètre embarqué dans les téléphones mobiles et les taxis équipés avec GPS, nous présentons ici notre travail sur la reconnaissance des comportements individuels, la capture des comportements communautaires et la démonstration des nouveaux services activés. En reconnaissant les comportements individuels, nous présentons la reconnaissance des activités physiques d'une personne avec les lectures de l'accéléromètre recueillies à partir des téléphones mobiles mis dans les poches autour de la zone pelvienne. Avec les empreintes GPS d'un taxi, nous résumons les comportements anormaux du transport des passagers pour un individu et améliorons l'efficacité de la reconnaissance de la méthode existante IBOAT. En outre, sur la base des observations dans la vie réelle, nous proposons une méthode pour détecter les événements de changement de service d’un taxi individuel. Avec des traces GPS à grande échelle et à l’aide des milliers de taxis, nous étudions les comportements anormaux pour le transport des passagers et les comportements de changement de travail et exploitons les stratégies de service de taxi. En outre, nous présentons un système intelligent de taxi comme une étude exemplaire des nouvelles applications qui s’appuie sur les comportements perçus individuelles et communautaires / The digital footprints collected from the prevailing sensing systems provide novel ways to perceive an individual's behaviors. Furthermore, large collections of digital footprints from communities bring novel understandings of human behaviors from the community perspective (community behaviors), such as investigating their characteristics and learning the hidden human intelligence. The perception of human behaviors from the sensing digital footprints enables novel applications for the sensing systems. Bases on the digital footprints collected with accelerometer-embedded mobile phones and GPS equipped taxis, in this dissertation we present our work in recognizing individual behaviors, capturing community behaviors and demonstrating the novel services enabled. With the GPS footprints of a taxi, we summarize the individual anomalous passenger delivery behaviors and improve the recognition efficiency of the existing method iBOAT by introducing an inverted index mechanism. Besides, based on the observations in real life, we propose a method to detect the work-shifting events of an individual taxi. With real-life large-scale GPS traces of thousands of taxis, we investigate the anomalous passenger delivery behaviors and work shifting behaviors from the community perspective and exploit taxi serving strategies. We find that most anomaly behaviors are intentional detours and high detour inclination won't make taxis the top players. And the spatial-temporal distribution of work shifting events in the taxi community reveals their influences. While exploiting taxi serving strategies, we propose a novel method to find the initial intentions in passenger finding. Furthermore, we present a smart taxi system as an example to demonstrate the novel applications that are enabled by the perceived individual and community behaviors
|
296 |
Architecture distribuée pour la détection d'activité dans un Espace IntelligentOudet, Jean-Philippe January 2011 (has links)
La présente étude porte sur la capacité d'améliorer la détection des Activités de la Vie Quotidienne, AVQ (ou ADL :"Activity of Daily Life") par l'utilisation de capteur [i.e. capteurs] de mouvements portés par l'occupant d'un habitat intelligent. Les données provenant de ces capteurs devraient fusionner avec les informations issues de l'appartement pour donner une information plus pertinente par le principe de synergie [21]. La solution choisie pour le dispositif porté par la personne est l'innovation principale du projet : un réseau de capteurs disposés à plusieurs endroits sur le corps, communicant sans fil entre eux et avec le contrôle de l'appartement. Les données extraites sont le mouvement relatif du corps, et plus spécifiquement des mains et du tronc, par rapport à la verticale. De par les propriétés de ces éléments - nécessairement petits, discrets - des MEMS seront utilisés pour satisfaire ces critères. Le projet repose sur la conception des dispositifs embarqués sur l'occupant dans l'optique d'en étendre les fonctionnalités à d'autres analyses tels [i.e. telles] que le son, la position dans l'environnement, le statut médical, etc. Pour prouver la faisabilité, des capteurs externes seront ajoutés pour compléter les informations de base et donc étendre la qualité des inférences sur les activités en cours. Le mouvement est une donnée facilement détectable de par sa relative simplicité de mise en oeuvre et il fournit une bonne base de travail pour étudier de façon systématique les différents points clés de l'étude : la communication, la synergie des informations, l'analyse des activités, etc.
|
297 |
Simulation numérique par éléments finis des grandes déformations des sols : application à la scarificationRenon, Nicolas 23 January 2002 (has links) (PDF)
Le déminage mécanique consiste à scarifier le sol avec une " charrue ", outil formé d'une lame en V munie de plusieurs dents, poussée par un engin : les dents déstructurent le sol et le font remonter devant la lame qui l'évacue sur le côté avec les mines qu'il contient. L'objectif de nos travaux est de mettre en œuvre la simulation numérique par éléments finis du problème fortement non linéaire issu de la modélisation de la scarification du sol. Le code d'éléments finis implicite Forge3®, dédié à la mise en forme des métaux, a été choisi comme support numérique. Il permet de prendre en compte les grandes déformations, en particulier à l'aide de son remaillage automatique. Nous avons dans un premier temps implanté dans ce code deux modèles élastoplastiques de comportement, l'un incompressible réservé aux sols fins saturés, purement cohésifs, l'autre compressible, fondé sur la notion d'état critique, pour les matériaux purement frottants ou cohésifs et frottants. Ces modèles adoucissants sont intégrés par un schéma de retour radial généralisé, au sein d'une formulation implicite du problème aux limites. Nous montrons que la matrice de raideur est non symétrique dans le cas compressible et que la symétrisation du système ne conduit pas à une approche robuste. Nous avons donc importé et testé un solveur itératif non-symétrique : Bi-CGSTAB. Nous avons validé la programmation de ces modèles sur des essais triaxiaux. Pour les comportements adoucissants, on constate des oscillations dans la relation contrainte/déformation passé le pic de contrainte. Ces difficultés numériques sont traitées par linéarisation et régularisation. Dans un deuxième temps, nous avons mis en œuvre des simulations de scarification pour différents niveaux de complexité : une dent seule, une dent + une tranche de lame, plusieurs dents, en faisant suffisamment avancer l'outil pour atteindre le régime stationnaire ; cela se révèle d'autant plus long que l'outil est large, ce qui nous amène au calcul très intensif. L'influence de paramètres géométriques comme l'inclinaison de la dent ou l'angle d'étrave du système a été mise en évidence, elle est qualitativement conforme aux observations expérimentales. L'étude de l'influence des paramètres des modèles de comportement montre la prépondérance de ceux liés à l'état critique, i.e. aux propriétés mécaniques après de grandes déformations. Enfin nous avons validé qualitativement le code en termes de modes d'écoulements et de répartition des efforts pour des outils multi-dents. La comparaison quantitative des efforts reste à affiner en revenant sur le comportement choisi et son implémentation.
|
298 |
Modélisation et méthodes d'évaluation de contrats gaziers: Approches par contrôle stochastiqueBernhart, Marie 11 March 2011 (has links) (PDF)
Le travail présenté dans cette thèse a été motivé par des problématiques posées par l'évaluation de contrats échangés sur le marché du gaz: les contrats de stockage et d'approvisionnement en gaz. Ceux-ci incorporent de l'optionalité et des contraintes, ce qui rend leur évaluation difficile dans un contexte de prix de matières premières aléatoires. L'évaluation de ces contrats mène à des problèmes de contrôle stochastique complexes: switching optimal ou contrôle impulsionnel et contrôle stochastique en grande dimension. La première partie de cette thèse est une revue relativement exhaustive de la littérature, mettant en perspective les différentes approches d'évaluation existantes. Dans une deuxième partie, nous considérons une méthode numérique de résolution de problèmes de contrôle impulsionnel basée sur leur représentation comme solution d'EDSRs à sauts contraints. Nous proposons une approximation à temps discret utilisant une pénalisation pour traiter la contrainte et donnons un taux de convergence de l'erreur introduite. Combinée avec des techniques Monte Carlo, cette méthode a été testée numériquement sur trois problèmes: gestion optimale de biomasse, évaluation d'options Swing et de contrats de stockage gaz. Dans une troisième partie, nous proposons une méthode pour l'évaluation d'options dont le payoff dépend de moyennes mobiles de prix sous-jacents. Elle utilise sur une approximation à dimension finie de la dynamique des processus de moyenne mobile, basée sur un développement en série de Laguerre tronquée. Les résultats numériques fournis incluent des exemples de contrats Swing gaziers à prix d'exercice indexés sur moyennes mobiles de prix pétroliers.
|
299 |
Conception et modélisation numérique de composants optiques en nanophotonique intégréeŚmigaj, Wojciech 22 September 2010 (has links) (PDF)
Cette thèse est consacrée à la conception et l'analyse théorique de différents composants en optique intégrée. Nous présentons un modèle de milieu effectif pour les cristaux photoniques (CPs) 2D qui rend compte des effets de surface, puis un algorithme pour la conception de réseaux antiréfléchissants grand-angle pour ces CPs. Ces réseaux permettent d'améliorer significativement la transmission à travers une lentille plate d'indice négatif. Nous proposons une nouvelle génération de circulateurs magnétooptiques compacts, fonctionnant dans un champ magnétique extérieur uniforme et constitués d'une cavité résonnante en anneaux circulaires couplée directement à des guides d'ondes standards. Nous généralisons la méthode multipolaire 2D aux matériaux gyrotropiques et la formulons sans « lattice sums » pour les structures périodiques. Enfin, nous décrivons en détail la méthode des éléments finis pour le calcul des modes propres des cavités 3D en anneaux circulaires et matériaux gyrotropiques.
|
300 |
Hydrodynamical simulations of detonations in superbursts./ Simulations hydrodynamiques de détonations dans les superbursts.Noël, Claire 19 October 2007 (has links)
In this thesis, we construct a new hydrodynamical algorithm able of handling general compressible reactive flow problems, based on a finite-volume method inspired by the original MUSCL scheme of van Leer (1979). The algorithm is of second-order in the smooth part of the flow and avoids dimensional splitting. It uses MPI to achieve parallelism, and includes an astrophysical equation of state and a nuclear reaction network. It proves to be robust to tests cases. In particular it reproduces quite well the reactive and non-reactive results obtained with two different numerical methods (Fryxell & al. 1989, Busegnies & al. 2007). Moreover the time-dependent results are in agreement with the corresponding steady state solution. This gives us confidence in applying it to an astrophysical situation which has never been studied, the propagation of a detonation in conditions relevant to superbursts. The algorithm is described in (Noel & al. 2007).
In a firt step we obtain the detonation profiles in pure carbon and in a mixture of carbon and iron. In both cases we underline the large difference between the total reaction length and the length on which some species burn. This difference leads to enormous numerical difficulties because all the length scales cannot be resolved at the same time in a single simulation. We show that the carbon detonation might be studied in a partial resolution approach like the one of Gamezo & al. (1999).
In a second step we construct a new reduced nuclear reaction network able to reproduce the energy production due to the photo-disintegrations of heavy elements, like ruthenium, which are thought to occur during superbursts in mixed H/He accreting systems. Using this new nuclear network we simulate detonations in mixture of carbon and ruthenium. An interesting feature is that, in this case, all the reaction lengths can be resolved in the same simulation. This makes the C/Ru detonations easier to study in future multi-dimensional simulations than the pure carbon ones (Noel & al. 2007b).
Finally we perform some numerical experiments which show that our algorithm is able to deal with initially inhomogeneous medium, and that the multi-dimensional simulations are attainable even if they are quite computational time consuming.
- B. Van Leer, J. Comp. Phys., 21, 101, 1979
- Fryxell, B.A., Muller, E., and Arnett, W.D., Technical report MPA 449, 1989
- Busegnies, Y., Francois, J. and Paulus, G., Shock Waves, 11, 2007
- Gamezo, V.N., Wheeler, J.C., Khokhlov, A.M., and Oran, E.S., ApJ, 512, 827, 1999
- Noël, C., Busegnies, Y., Papalexandris, M.V. & al., A&A, 470, 653, 2007
- Noël, C., Goriely, S., Busegnies, Y. & Papalexandris, M.V., submitted to A&A, 2007b
/
Un algorithme parallèle basé sur une méthode aux volumes finis inspirée du schéma MUSCL de Van Leer (1979) a été construit. Il a été développé sur base de la méthode de Lappas & al. (1999) qui permet de résoudre simultanément toutes les dimensions spatiales. Cette méthode se base sur la construction de surfaces appropriées dans l'espace-temps, le long desquelles les équations de bilan se découplent en équations plus simples à intégrer. Cet algorithme est actuellement le seul à éviter le "splitting" des dimensions spatiales. Dans les modèles conventionnels (PPM, FCT, etc.), l'intégration spatiale des équations est réalisée de manière unidimensionnelle pour chaque direction.
Un réseau de réactions nucléaires ainsi qu'une équation d'état astrophysique ont été inclus dans l'algorithme et celui-ci a ensuite été soumis à une grande variété de cas tests réactifs et non réactifs. Il a été comparé à d'autres codes généralement utilisés en astrophysique (Fryxell & al. 1989, Fryxell & al. 2000, Busegnies & al. 2007) et il reproduit correctement leurs résultats. L'algorithme est décrit dans Noël & al. (2007).
Sur base de cet algorithme, les premières simulations de détonation dans des conditions thermodynamiques représentatives des Superbursts ont été réalisées. Différentes compositions du milieu ont été envisagées (carbone pur, mélange de carbone et de fer, mélange de carbone et de cendres du processus rp). Dans la plupart des systèmes où des Superbursts ont été observés, la matière accrétée est un mélange d'hydrogène et d'hélium. Dans ce cas, des phases de combustion précédant le Superburst produisent des nucléides plus lourd que le fer (Schatz & al. 2003). Ces nucléides peuvent être photodésintégrés durant le Superburst. Pour prendre en compte ces réactions endothermiques de photodésintégration, nous avons construit un nouveau réseau réduit de réactions nucléaires qui a été incorporé dans l'algorithme hydrodynamique (Noël & al. 2007b). Ce réseau réduit reproduit globalement l'énergétique d'un réseau complet et a permis de faire la première simulation numérique de détonation dans des conditions caractéristiques de systèmes accréteurs d'un mélange hydrogène-hélium.
Finallement quelques simulations multidimensionelles préliminaires ont éte réalisées.
- Busegnies, Y., Francois, J. and Paulus, G., Shock Waves, 11, 2007
- Fryxell, B.A., Muller, E., and Arnett, W.D., Technical report MPA 449, 1989
- Fryxell, B.A., Olson, K., Ricker, P. & al., ApJS, 131, 273, 2000
- Lappas, T., Leonard, A. and Dimotakis, P.E., SIAM J. Sci. Comput., 20, 1064, 1999
- Noël, C., Busegnies, Y., Papalexandris, M.V. & al., A&A, 470, 653, 2007
- Noël, C., Goriely, S., Busegnies, Y. & Papalexandris, M.V., submitted to A&A, 2007b
- Röpke, F. K. PhD thesis, Technischen Universitat Munchen, 2003
- Schatz, H., Bildsten, L., Cumming, A. and Ouellette, M., Nuclear Physics A, 718, 247, 2003
- Van Leer, B. Comp. Phys., 21, 101, 1979
- Weinberg, N.N. and Bildsten, L., ArXiv e-prints, 0706.3062, 2007
|
Page generated in 0.0678 seconds