The motivation of this study is to understand the seismic anisotropy observations from various subduction regions of the world. In subduction zone backarcs both trench-parallel and trench-normal seismic anisotropy, or fast wave polarization direction of shear wave, are observed. In the mantle the general assumption is that seismic anisotropy is caused by Lattice Preferred Orientation (LPO) of olivine minerals and that the direction of anisotropy is an indicator of the direction of mantle flow. The complex pattern of seismic anisotropy observations suggests that the flow geometry in the vicinity of subduction zones differs at different subduction zones with some subduction zones having trench perpendicular flow, consistent with corner flow in the mantle wedge while other subduction zones have trench parallel flow, consistent with a mode of flow where material from the mantle wedge flows around the edges of the slab. It should be noted that the direction of LPO orientation can also be modified by the presence or absence of water, pressure, and temperature in the mantle and that it is possible that the difference in anisotropy observations reflects a difference in water content or thermal structure of back arcs. The aim of this study is to test whether the flow geometry of mantle in numerical subduction calculations can influence the direction of seismic anisotropy and if we parameters that control the pattern of flow can be identified. In this study we explicitly assume that seismic anisotropy occurs only due to plastic and dynamic re-crystallization of mantle mineral forming LPO. To approach the problem two different models are formulated. In one of the models the trench evolves self-consistently, with no prescribed artificial zones of weakness. The self-consistent model has a sticky-air layer at the top of the model domain that mimics a "free-surface." The other model has the same initial conditions but a trench-migration velocity boundary condition is imposed to the model. The mantle flow pattern for the self-consistent model is consistent with the 2D corner flow with no flow around the trench and no trench migration. However when the trench-migration velocity boundary condition is imposed, 3D flow around the mantle is observed. The stress field from these simulations are used to calculated instantaneous strain axis directions which correlate with LPO directions. The LPO orientations are measured from the models showing that the seismic-anisotropy direction is primarily trench-perpendicular for both models. Because the models have different flow patterns, the trench-perpendicular anisotropy alignment that is calculated for both the models is a bit puzzling. It could be that factors such as high temperature and non-linear rheology cause the LPO direction to align trench perpendicular in both the cases. It can also be possible that the 3D vertical flow is not strong enough to cause change in orientation of the LPO direction. From the present study it can be concluded that by looking at the LPO direction nature of mantle flow might not be predicted. This suggests that in addition to flow direction other factors such as the presence of water in mantle wedge, pressure, and high temperature due to viscous coupling modify the seismic anisotropy directions. / Master of Science

The Boosted Fast Flux Loop (BFFL) will expand the Advanced Test Reactor (ATR) at Idaho National Laboratory. Part of the BFFL is a new corrosion test cap section for testing in the ATR. The corrosion test cap section was designed with parallel channels to reduce the pressure drop and allow coolant contact with specimens. The fluid experiment conducted by Idaho State University found the pressure drop not characteristic of parallel channel flow but greater than without parallel channels. A Computation Fluid Dynamics simulation using STAR-CCM+ was conducted with the objectives of showing sufficient flow through the test cap section for a corrosion test, verifying the fluid experiment's validity, and explaining the abnormal pressure drop. The simulation used a polyhedral volume mesh and the k-e turbulent model with segregated equations. Convergence depended on a low continuity residual and an unchanging pressure drop result. The simulation showed the same pattern as the fluid experiment. The simulation provided evidence of flow through the test cap section needed for a corrosion test. The specimen holding assembly was found to be a small contributor to the pressure drop. The counterintuitive pressure drop was found to be the sum of many factors produced from the geometry of the test cap section. The inlet of the test cap section behaved as a diverging nozzle before a sudden expansion into the test cap section chamber with both creating a pressure drop. The chaotic flow inside the chamber gave rise to pressure loss from mixing. The fluid exited the chamber through a sudden contraction to a converging nozzle behaving exit, again, producing a pressure drop. By varying the flow rate in the simulation, a disturbance in the flow where the gap fluid separated into the parallel channels was found at high flow rates. At low flow rates the pressure drop anomaly was not found. The corrosion test cap section could be used in the ATR but with a higher pressure drop than desirable. The design of the corrosion test cap section created the abnormal pressure drop.

In this study a three-dimensional mathematical model is developed to simulate the flow field and mass transfer in a PEM fuel cell. In the model, the effects of the different flow channel structures in heterogeneous composite carbon fiber bipolar plates and traditional hard surface bipolar plates on the performance are studied. The results show that, the cell performance with the heterogeneous composite carbon fiber bipolar plates have better performance than that with the traditional hard surface bipolar plates, whether in the parallel flow channel structures or the serpentine flow channel structures. The reason is that, the heterogeneous composite carbon fiber ribs are porous material, so it allows the reactants and products transport uniformly even in the rib zone. This greatly improved the mass transfer and the gases distribution in the fuel cell. With the traditional bipolar plates, the reactants can only enter the reaction zone from the side of carbon cloth under ribs, so that the performance in this area under rib is relatively poor. In the simulation of the flow channel structures, we detect that, due to the single inlet serpentine flow channel have stronger convective effects that forced reactants to flow through the whole reaction zones, so it has better performance at high current density than in the singles inlet parallel flow channel. In addition, the results also show that, higher fuel stoichiometric number and operated pressure and properly humidified at anode will all improve the performance of the fuel cell.

Le transport des particules et de la chaleur dans la zone de bord des tokamaks joue un rôle déterminant à la fois sur les performances du plasma confiné et sur l’extraction de la puissance et ainsi la durée de vie des composants face au plasma. C’est dans ce contexte que s’inscrit ce travail de thèse, qui porte sur le rôle joué par les écoulements transverses au champ magnétique dans l’équilibre entre dynamique parallèle et dynamique perpendiculaire qui régit la région périphérique d’un tokamak. Ces écoulements peuvent produire des asymétries poloïdales du dépôt de chaleur et de particules sur les composants face au plasma, et plus généralement des asymétries des diverses quantités dans le plasma. Les vitesses de dérive radiale sont d’origine électrique (liées à la présence d’un champ électrique radial résultant de l’équilibre des charges) ou liées aux effets de la géométrie toroïdale induisant une inhomogénéité du champ magnétique (vitesse de gradient-courbure). Pour progresser dans la compréhension de ces phénomènes, la modélisation numérique du transport et de la turbulence en géométrie complexe est indispensable. En complément, des outils de diagnostic synthétique permettant de modéliser les processus de mesure dans les plasmas numériques sont développés pour permettre une comparaison réaliste entre modèles et expériences. La modélisation des vitesses de dérive perpendiculaire a été introduite dans le code SOLEDGE2D décrivant le transport de la densité, quantité de mouvement et énergie d’un plasma de tokamak. Nous avons d’abord étudié l’impact d’un champ électrique prescrit sur les équilibres plasma, pour comprendre les mécanismes à l’origine des asymétries du plasma et étudier l’établissement d’écoulement parallèle et d’asymétrie du dépôt de chaleur sur les composants face au plasma. Nous avons ensuite implémenté un modèle auto-consistant de résolution du potentiel électrique dans les équations fluides de SOLEDGE2D afin de comprendre l’équilibre du champ électrique et d’étudier l’effet de la configuration magnétique du tokamak et de la vitesse de gradient-courbure sur ce dernier. Dans la deuxième partie de cette thèse, un diagnostic synthétique permettant de modéliser les mesures expérimentales de rétro-diffusion Doppler a été développé et testé en vue d’être appliqué aux simulations du code fluide 3D turbulent, TOKAM3X. Ce diagnostic permet de mesurer la vitesse perpendiculaire du plasma à partir du mouvement des fluctuations de densité. Il a été utilisé ici pour comparer les asymétries de vitesse observées expérimentalement aux asymétries mesurées dans les simulations numériques. / The transport of heat and particles in the edge of tokamaks plays a key role in both the performance of the confined plasma and the extraction of power and thus the lifetime of the plasma facing components. It’s in this context that this thesis is inscribed, which focuses on the role played by the transverse magnetic field flows in the balance between parallel and perpendicular dynamic that governs the edge region of a tokamak. These flows can produce poloidal asymmetries of heat and particles deposit on plasma facing components and generally asymmetries of various amounts in plasma. The radial drift velocities are due to the presence of a radial electric field resulting from charge balance (electric drift velocity) or related to effects of the toroidal geometry inducing a magnetic field inhomogeneity (curvature drift velocity). To advance the understanding of these phenomena, numerical modeling of transport and turbulence in complex geometries is essential. In addition, synthetic diagnostic tools for modeling the measurement process in numerical plasmas are developed to enable a realistic comparison between models and experiments. Modeling of perpendicular drift velocities was introduced into the SOLEDGE2D code describing the transport of the density, momentum and energy of a tokamak plasma. We first studied the impact of a prescribed electric field on plasma equilibrium to understand the mechanisms behind plasma asymmetries and study the establishment of parallel flows and asymmetry of the heat flux on plasma facing components. Then we implemented a self-consistent model solving the electric potential in SOLEDGE2D fluid equations to understand the equilibrium of the electric field and to study the effect of the magnetic configuration of the tokamak and the curvature drift velocity on it. In the second part of this thesis, a synthetic diagnosis modeling the experimental measurements of Doppler backscattering was developed and tested in order to be applied to simulations of 3D turbulent fluid code TOKAM3X. This diagnosis measures the perpendicular velocity of the plasma from the movement of the density fluctuations. It was used to compare the perpendicular velocity asymmetries observed experimentally to asymmetries measured in numericalsimulations.