Com a finalidade de estudar processos de fratura em concreto simples, este trabalho propõe uma escala intermediária de análise. Sendo assim, é apresentado um modelo simplificado de duas fases, o qual representa alguns elementos da microestrutra do concreto. Uma fase representa os agregados e a outra a matriz de cimento. Os agregados são considerados elástico-lineares e são representados pelos elementos finitos. Toda a não-linearidade do concreto e do processo de fratura é considerada na matriz de cimento. Tal fase é representada por interfaces coesivas, de espessura zero, entre os elementos finitos. Enquanto o modelo considera propriedades microscópicas do concreto como entrada, tem-se o comportamento macroscópico da estrutura como saída, o que é uma grande vantagem da escala proposta. Efeitos como a relação água/cimento e a densidade de agregados são introduzidos no modelo. Os resultados mostram que o modelo é capaz de captar, pelo menos qualitativamente, efeitos os quais somente são possíveis de obter com um modelo de três fases. / An intermediate scale of analysis is proposed here in order to study the fracture of plain concrete. A simplified two phases model is proposed as a concrete microstructure. One phase is the aggregate and other is the cement matrix. Aggregate is considered elastic linear and is represented by finite elements. All the non-linearity of the concrete and the fracture process is considered in the cement matrix. This phase is represented by a zero thickness cohesive interface between finite elements. While the model considers microscopic properties of the concrete as input, it delivers the macroscopic behavior of the structure as output, which is a great advantage of the proposed scale. Effects like water-cement ratio and density of aggregates are introduced in the model. Results show the model is able to capture, at least qualitatively, effects only possible to obtain by a three phases model.

Com a finalidade de estudar processos de fratura em concreto simples, este trabalho propõe uma escala intermediária de análise. Sendo assim, é apresentado um modelo simplificado de duas fases, o qual representa alguns elementos da microestrutra do concreto. Uma fase representa os agregados e a outra a matriz de cimento. Os agregados são considerados elástico-lineares e são representados pelos elementos finitos. Toda a não-linearidade do concreto e do processo de fratura é considerada na matriz de cimento. Tal fase é representada por interfaces coesivas, de espessura zero, entre os elementos finitos. Enquanto o modelo considera propriedades microscópicas do concreto como entrada, tem-se o comportamento macroscópico da estrutura como saída, o que é uma grande vantagem da escala proposta. Efeitos como a relação água/cimento e a densidade de agregados são introduzidos no modelo. Os resultados mostram que o modelo é capaz de captar, pelo menos qualitativamente, efeitos os quais somente são possíveis de obter com um modelo de três fases. / An intermediate scale of analysis is proposed here in order to study the fracture of plain concrete. A simplified two phases model is proposed as a concrete microstructure. One phase is the aggregate and other is the cement matrix. Aggregate is considered elastic linear and is represented by finite elements. All the non-linearity of the concrete and the fracture process is considered in the cement matrix. This phase is represented by a zero thickness cohesive interface between finite elements. While the model considers microscopic properties of the concrete as input, it delivers the macroscopic behavior of the structure as output, which is a great advantage of the proposed scale. Effects like water-cement ratio and density of aggregates are introduced in the model. Results show the model is able to capture, at least qualitatively, effects only possible to obtain by a three phases model.

Com a finalidade de estudar processos de fratura em concreto simples, este trabalho propõe uma escala intermediária de análise. Sendo assim, é apresentado um modelo simplificado de duas fases, o qual representa alguns elementos da microestrutra do concreto. Uma fase representa os agregados e a outra a matriz de cimento. Os agregados são considerados elástico-lineares e são representados pelos elementos finitos. Toda a não-linearidade do concreto e do processo de fratura é considerada na matriz de cimento. Tal fase é representada por interfaces coesivas, de espessura zero, entre os elementos finitos. Enquanto o modelo considera propriedades microscópicas do concreto como entrada, tem-se o comportamento macroscópico da estrutura como saída, o que é uma grande vantagem da escala proposta. Efeitos como a relação água/cimento e a densidade de agregados são introduzidos no modelo. Os resultados mostram que o modelo é capaz de captar, pelo menos qualitativamente, efeitos os quais somente são possíveis de obter com um modelo de três fases. / An intermediate scale of analysis is proposed here in order to study the fracture of plain concrete. A simplified two phases model is proposed as a concrete microstructure. One phase is the aggregate and other is the cement matrix. Aggregate is considered elastic linear and is represented by finite elements. All the non-linearity of the concrete and the fracture process is considered in the cement matrix. This phase is represented by a zero thickness cohesive interface between finite elements. While the model considers microscopic properties of the concrete as input, it delivers the macroscopic behavior of the structure as output, which is a great advantage of the proposed scale. Effects like water-cement ratio and density of aggregates are introduced in the model. Results show the model is able to capture, at least qualitatively, effects only possible to obtain by a three phases model.

The research presented in this thesis is a contribution to the modeling and understanding of the dynamic behavior of proton exchange membrane fuel cells (PEMFCs). A time-dependent, two-phase non-isothermal model of the membrane electrode assembly was developed and implemented using the finite element method. In addition to solving a phenomenological transport equation for water in the membrane, the model takes into consideration the non-equilibrium water sorption to better capture some of the dynamic characteristics of water transport in the MEA. Mass transfer using Fickian diffusion is implemented in the model. Two different models describing the electrochemical reactions in the catalyst layer including a macro-homogeneous model and an agglomerate model, are also implemented. Conservation of energy is included in the solution procedure in order to assess the impact of thermal effects on the dynamics of the transport in the MEA. For the purpose of model and concept validation, the model was first solved in a steady two-dimensional mode for a through- plane computational domain using a commercial software package, COMSOL Multiphysics version 3.2b. The impact of using a single- and two-phase modeling approaches was evaluated, and the predicted current-voltage performance characteristic are found in good agreement with the experimental data available in the literature. In addition, the developed model was benchmarked against a finite element-based in-house code for further validation and to evaluate numerical accuracy and computational performance. Transient simulations of operation under dynamic voltage sweeps are presented, and parametric studies are conducted to investigate the impact of various model, operation and transport properties on the predicted dynamic cell performance. In particular, the rate of load change, the difference in water content between the anode and cathode, and the water sorptions rate are shown to have significant impact on cell performance in unsteady operation, especially at higher current densities. Parametric studies also address the sensitivity of the model results to physical properties, highlighting the importance of accurately determining certain physical properties of the fuel cell components. Finally, the application of the model to air-breathing fuel cells provides further insight into the dynamic performance characteristic of such type of fuel cells.

L’opto-pyrotechnie (amorçage de la détonation par système optique) est l’une des innovations les plus prometteuses en termes de fiabilité, de sécurité et de performances pour les futurs lanceurs spatiaux. Le but de la thèse est d’étudier et de modéliser le premier des deux étages d’un Détonateur Opto-Pyrotechnique, constitué d’un explosif confiné dans une chambre de combustion fermée où se déroulent les premières phases d’une Transition Déflagration-Détonation. L’amorçage par laser de l’explosif puis la combustion en chambre isochore sont traités par le code EFAE, lequel est couplé au logiciel LS-DYNA qui simule la déformation et la rupture du disque de fermeture de la chambre, puis la propulsion du projectile résultant vers le second étage. En parallèle, diverses techniques expérimentales (adsorption de gaz, vélocimétrie hétérodyne, microscopie) ont mis en valeur plusieurs procédés physiques, ce qui a permis de tester le couplage entre EFAE et LS-DYNA, puis de déterminer et de hiérarchiser les paramètres affectant les critères industriels. / Opto-pyrotechnics (ignition of detonation by optical systems) is one of the most promising innovations to improve reliability, safety and performances on future space launchers. This thesis aims at studying and modeling the first stage from a two-stage opto-pyrotechnic detonator that consists of a condensed explosive confined in a closed combustion chamber, in which the beginning of a Deflagration-to-Detonation Transition occurs. The laser ignition of the explosive and its isochoric combustion are modeled by the EFAE code. This code is coupled with LS-DYNA software to deal with the deformation and the rupture of the metallic disk that closes the combustion chamber, and then with the subsequent propulsion of the projectile to the second stage. In parallel, various experimental technics (gas adsorption, photonic Doppler velocimetry, microscopy) have underlined several physical processes that allow first to test the coupling between EFAE and LS-DYNA, then to determine and classify influent parameters that affect the industrial specifications.

No description available.

Cette thèse traite de la modélisation des tsunamis, des grandes échelles de propagation aux impacts sur desstructures côtières. Un inventaire des phénomènes physiques associés est établi et des modèles adéquats sontprésentés. Une étude numérique avec le modèle de Saint-Venant est effectuée avec le développement d’uneméthode de raffinement de maillage à seuil automatique. La simplicité et les performances de l’approche sontdémontrées. Pour améliorer la précision des prévisions, un système original approchant le modèle Serre-Green-Nahgdi est investigué. Une méthode pour prendre en compte la dissipation d’énergie au déferlement estproposée. Ce modèle permet d’envisager la modélisation fine de la propagation et de l’arrivée à la côte destsunamis dispersifs et non linéaires en des temps de calcul acceptables. Les différents types d’impact sur desstructures sont modélisés grâce à un modèle diphasique compressible permettant de considérer les écoulementsà phases séparées et les milieux aérés. Pour envisager une résolution à tous les régimes, des schémas TousMach sont investigués. Un schéma Tous Mach à variation totale limitée est proposé. Grâce à cette approche,des impacts incompressibles et compressibles sont investigués avec le même modèle. Les impactsd’écoulements aérés induisent des pressions moins élevées mais sur des temps plus longs que leurshomologues en phases pures. Bien que le schéma Tous Mach proposé soit moins sujet aux oscillationsnumériques que les préconditionnements classiques de la littérature, des oscillations non physiques à basnombre de Mach sont mises en évidence sur certains cas tests. Pour finir, une méthode de couplage entremodèles de propagation et d’impact est proposée, afin de pouvoir simuler un tsunami finement avec desmodèles appropriés à chacune de ces phases. / This thesis deals with tsunami modeling, from the large propagation scales to impacts on coastal structures. Aninventory of the associated physical phenomena is given and some adequat models are presented. A numericalstudy is carried out with the Saint-Venant model with the devellopment of an automatic refinment adaptivemesh method. The resolution efficiency and simplicity is justified. To increase the accuracy forecasts, anoriginal system which approximates the Serre-Green-Naghdi model is investigated. A breaking wave methodassociated with this model is proposed. This dissipative model allows thinking about accurate dispersive nonlinear tsunamis simulations up to the coast. Several wave impacts on stuctures are investigated with a generaltwo-phase model allowing separate phases as well as aerated impact studies. The all Mach regime numericaldifficulties are investigated. A new all Mach scheme with limited total variation bound is proposed. Thanks tothis approach, incompressible and compressible impacts are investigated with the same model. Aerated impactsare shown to give smaller impact pressure but on longer time than pure fluid impacts. In spite that the all Machscheme proposed reduces the numerical oscillations of classical literature preconditionning, some unphysicaloscillations are highlighted on some test cases. At the end, a coupling method is proposed in order to accuratlymodel the propagation and the impact of a tsunami with appropriated models for each phases.

Heat transfer and solidification models were developed for use in a numerical model of a continuous caster to provide a means of predicting how the developing shell would react under variable operating conditions. Measurement data of the operating conditions leading up to a breakout occurrence were provided by an industrial collaborator and were used to define the model boundary conditions. Steady-state and transient simulations were conducted, using boundary conditions defined from time-averaged measurement data. The predicted shell profiles demonstrated good agreement with thickness measurements of a breakout shell segment – recovered from the quarter-width location. Further examination of the results with measurement data suggests pseudo-steady assumption may be inadequate for modeling shell and flow field transition period following sudden changes in casting speed. An adaptive mesh refinement procedure was established to increase refinement in areas of predicted shell growth and to remove excess refinement from regions containing only liquid. A control algorithm was developed and employed to automate the refinement procedure in a proof-of-concept simulation. The use of adaptive mesh refinement was found to decrease the total simulation time by approximately 11% from the control simulation – using a static mesh.

Dans un contexte industriel, l’utilisation de modèles diphasiques d’ordre réduit est nécessaire pour pouvoir effectuer des simulations numériques prédictives d’injection de combustible liquide dans les chambres de combustion automobiles et aéronautiques, afin de concevoir des équipements plus performants et moins polluants. Le processus d’atomisation du combustible, depuis sa sortie de l’injecteur sous un régime de phases séparées, jusqu’au brouillard de gouttelettes dispersées, est l’un des facteurs clés d’une combustion de bonne qualité. Aujourd’hui cependant, la prise en compte de toutes les échelles physiques impliquées dans ce processus nécessite une avancée majeure en termes de modélisation, de méthodes numériques et de calcul haute performance (HPC). Ces trois aspects sont abordés dans cette thèse. Premièrement, des modèles de mélange, dérivés par le principe variationnel de Hamilton et le second principe de la thermodynamique sont étudiés. Ils sont alors enrichis afin de pouvoir décrire des pulsations des interfaces au niveau de la sous-échelle. Des comparaisons avec des données expérimentales dans un contexte de milieux à bulles permettent de vérifier la cohérence physique des modèles et de valider la méthodologie. Deuxièmement, une stratégie de discrétisation est développée, basée sur une séparation d’opérateur, permettant la résolution indépendante de la partie convective des systèmes à l’aide de solveurs de Riemann approchés standards et les termes sources à l’aide d’intégrateurs d’équations différentielles ordinaires. Ces différentes méthodes répondent aux particularités des systèmes diphasiques compressibles, ainsi qu’au choix de l’utilisation de maillages adaptatifs (AMR). Pour ces derniers, une stratégie spécifique est développée : il s’agit du choix de critères de raffinement et de la projection de la solution d’une grille à une autre (plus fine ou plus grossière). Enfin, l’utilisation de l’AMR dans un cadre HPC est rendue possible grâce à la bibliothèque AMR p4est, laquelle a montré une excellente scalabilité jusqu’à plusieurs milliers de coeurs de calcul. Un code applicatif, CanoP, a été développé et permet de simuler des écoulements fluides avec des méthodes de volumes finis sur des maillages AMR. CanoP pourra être utilisé pour des futures simulations d’atomisation liquide. / In an industrial context, reduced-order two-phase models are used in predictive simulations of the liquid fuel injection in combustion chambers and help designing more efficient and less polluting devices. The combustion quality strongly depends on the atomization process, starting from the separated phase flow at the exit of the nozzle down to the cloud of fuel droplets characterized by a disperse-phase flow. Today, simulating all the physical scales involved in this process requires a major breakthrough in terms of modeling, numerical methods and high performance computing (HPC). These three aspects are addressed in this thesis. First, we are interested in mixture models, derived through Hamilton’s variational principle and the second principle of thermodynamics. We enrich these models, so that they can describe sub-scale pulsations mechanisms. Comparisons with experimental data in a context of bubbly flows enables to assess the models and the methodology. Based on a geometrical study of the interface evolution, new tracks are then proposed for further enriching the mixture models using the same methodology. Second, we propose a numerical strategy based on finite volume methods composed of an operator splitting strategy, approximate Riemann solvers for the resolution of the convective part and specific ODE solvers for the source terms. These methods have been adapted so as to handle several difficulties related to two-phase flows, like the large acoustic impedance ratio, the stiffness of the source terms and low-Mach issues. Moreover, a cell-based Adaptive Mesh Refinement (AMR) strategy is considered. This involves to develop refinement criteria, the setting of the solution values on the new grids and to adapt the standard methods for regular structured grids to non-conforming grids. Finally, the scalability of this AMR tool relies on the p4est AMR library, that shows excellent scalability on several thousands cores. A code named CanoP has been developed and enables to solve fluid dynamics equations on AMR grids. We show that CanoP can be used for future simulations of the liquid atomization.

No description available.