Basé sur l'analyse conjointe et complémentaire de ~10 ans de données altimétriques et de profils verticaux de température et de salinité acquis par les profileurs Argo, l'objectif principal de cette thèse est d'étudier en détail les caractéristiques de surface et la structure verticale des tourbillons dans les 4 grands systèmes d'upwelling mondiaux (EBUS Pérou-Chili, Californie, Canaries et Benguela), qui partagent une dynamique à grande échelle relativement similaire. Les résultats principaux montrent que l'ensemble des tourbillons détectés sur les cartes altimétriques dans les 4 EBUS ont des propriétés physiques relativement proches, avec une forte proportion de tourbillons de faibles dimensions ayant une durée de vie relativement courte. Au contraire, la faible partie des tourbillons échantillonnés par les profileurs Argo montre des dimensions bien plus grandes pour des durées de vie plus longues. La sur-représentation de ces grandes échelles dans le jeu des tourbillons échantillonnés par les profileurs Argo est donc le biais majeur des résultats obtenus sur la structure verticale des tourbillons. L'analyse des profils moyens d'anomalies de température et de salinité acquis par les profileurs Argo dans les tourbillons révèle une forte hétérogénéité entre les 4 EBUS, mais également au sein de chacun de ces systèmes d'upwelling. Les structures verticales des tourbillons sont fortement liées à l'hydrologie et à la dynamique locale. Par exemple, les tourbillons de subsurface du PCUS occupent plutôt la partie Sud du système, alors que les tourbillons intensifiés en surface sont présents à la limite Nord de la gyre subtropicale. Dans le CALUS, la présence de Cuddies se devine à la côte, les autres sous-régions contenant des tourbillons dont les anomalies sont majoritairement intensifiées à la base de la pycnocline, avec comme particularité un changement de signe lié à la salinité des couches superficielles. Le CANUS est peuplé de tourbillons très différents en fonction de leur position par rapport à la zone frontale du Cap Vert : au Nord, des tourbillons intenses en subsurface, avec la présence de quelques Meddies ; au Sud, des tourbillons très superficiels. De même dans le BENUS, la partie au Nord de 15°S contient plutôt des tourbillons intensifiés dans les couches de surface, alors qu'au Sud du front d'Angola-Benguela, les tourbillons présentent des anomalies fortes sur une grande partie de la colonne d'eau. L'extrême Sud de cet EBUS est également le lieu de passage des Anneaux des Aiguilles. Les contributions de l'advection isopycnale et du déplacement de la colonne d'eau sur la verticale nous permettent d'affiner la description de ces différentes structures. La présence de grands types de tourbillons, à la morphologie bien distincte (cœur de surface, de subsurface, grande extension verticale, tourbillons fortement intensifiés, etc.), est confirmée par l'étude des tourbillons spécifiquement générés dans la bande côtière. L'analyse Lagrangienne de ces tourbillons nous permet également de décrire l'évolution temporelle de leur structure verticale, qui montre une homogénéité temporelle inattendue. Cette thèse présente donc plusieurs outils facilement applicables dans différentes régions océaniques pour caractériser la structure thermohaline des tourbillons et fournit pour la première fois une description des grands types de tourbillons peuplant les EBUS, soulignant la grande diversité de la mésoéchelle. / Merging ~10 years of altimetry maps and vertical profiles provided by Argo floats, we aim to study in details the eddy's surface characteristics and vertical structure in the 4 major Eastern Boundary Upwelling Systems (EBUS : Peru-Chile, California, Canaries and Benguela), sharing similar large-scale dynamics. Our main results show that the eddies detected on altimetry maps in the 4 EBUS have close physical properties, with a lot of small-scale structures (radius < 40km, amplitude < 1cm and lifetime < 30 days). In contrast, the few eddies sampled by Argo floats have larger dimensions (radius of ~90-140 km, amplitude of ~3-7 cm) and longer lifetimes (6-10 months). The major bias with the analyzed vertical structure is the over representation of these large-scale eddies. The temperature and salinity anomaly mean profiles acquired by Argo floats surfacing within eddies reveals a strong heterogeneity between each of the 4 EBUS, but also within them. The eddies' vertical structure is strongly influence by the local hydrology and dynamics. For example, the subsurface-intensified eddies of the PCUS tend to be located in the Southern part of this EBUS, while the surface-intensified eddies are preferentially located near the Northern boundary of the subtropical gyre. In the CALUS, we can identify Cuddies in some coastal sub-regions, but in this EBUS, most of the eddies are intensified at the base of the pycnocline, with a reversal of the salinity anomaly compared to the surface layers. In the CANUS, the Cape Verde frontal zone separates distinct subsurface-intensified eddies and some Meddies in the North, from the Southern part, where eddies are surface-intensified. In the same way, the Angola-Benguela Front of the BENUS separates the surface-intensified eddies in the North from strong, deep-reaching anomalies in the South. The Southern-most part of the BENUS is also a preferential pathway for the large Agulhas Rings and their associated cyclones. The respective contributions of isopycnal advection and vertical displacement improve the description of these very diverse structures. The presence of several eddy-types, with distinct thermohaline properties (surface or subsurface-intensified, deep vertical extend, intense or not, etc.) is confirmed by the study of eddies generated in the coastal area of each EBUS. Their Lagrangian analysis allows us to describe the temporal evolution of their vertical structure, which shows an unexpected temporal homogeneity. This manuscript presents different efficient tools used to analyze the surface characteristics, the thermohaline properties and the temporal evolution of mesoscale eddies in the 4 major EBUS, highlighting their diversity.

Bio-optical variability of surface waters in Northeastern Gulf of Mexico (NEGOM) was examined using satellite and in situ data. Relatively high chlorophyll-a concentration (chl>=1 mg m-3) and high colored dissolved organic mater (ag443>=0.1 m-1) were generally observed inshore, near major river mouths, and in plumes of Mississippi River water that extended offshore during the three consecutive summer seasons (1998, 1999, and 2000). River discharge dominated chlorophyll-a concentration variability inshore, particularly near major river mouths. Strong interannual variability in chlorophyll-a concentration was observed inshore from Escambia to Tampa Bay region during the winter to spring transition, which was different in 1998 compared to the winter to spring transition in 1999 and 2000. This was related to higher fresh water discharge during the 1997-1998 El Niño-Southern Oscillation event as well as strong upwelling in spring 1998. The Mississippi plume extended >500 km southeast of the Mississippi delta and up to the Florida Keys was observed for the periods extending over 14 weeks between May and September every year of the study. In general, ag443 covaried linearly and inversely with salinity inshore during spring and fall, indicating conservative mixing. The NEGOM salinity-ag443 relationship of fall 1998, i.e., Salinity=36.59-29.86*ag443 (n=8771, r2=0.86; 0.01<=ag443<=0.52, 16 <=S<=36), served as the best predictor of NEGOM salinity based on in situ ag443 observations for spring and fall seasons from all years (<3% mean percentage errors; corresponding to <1.03 psu). This may help estimate salinity from satellite ocean color data, but further testing using data from multiple years is needed to improve such relationship. While river discharge was an important source of colored dissolved organic matter (CDOM), phytoplankton blooms also contributed to CDOM formation in the NEGOM. Using a pigment index of phytoplankton taxonomic groups, the variability in biomass proportion of microphytoplankton explained up to 76% of the variability of the average of normalized phytoplankton absorption coefficients (545, 625, and 673 nm). The clorophyll-specific absorption coefficient, a*ph(440), varies by a factor of 7 (0.02-0.15 m2mg-1). Particle size and pigment composition played important roles in determining a*ph(440) variability. This must be accounted for in chlorophyll-a concentration algorithms based on aph.

In this thesis, parameterizations of non-linear interactions in quasigeostrophic (QG) flows for severely truncated models (STM) and Large Eddy Simulations (LES) are studied. Firstly, using Direct Numerical Simulations (DNS), atmospheric barotropic flows over topography are examined, and it is established that such flows exhibit multiple equilibrium states for a wide range of parameters. A STM is then constructed, consisting of the large scale zonal flow and a topographic mode. It is shown that, qualitatively, this system behaves similarly to the DNS as far as the interaction between the zonal flow and topography is concerned, and, in particular, exhibits multiple equilibrium states. By fitting the analytical form of the topographic stationary wave amplitude, obtained from the STM, to the results obtained from DNS, renormalized dissipation and rotation parameters are obtained. The usage of renormalized parameters in the STM results in better quantitative agreement with the DNS.¶ In the second type of problem, subgrid-scale parameterizations in LES are investigated with both atmospheric and oceanic parameters. This is in the context of two-level QG flows on the sphere, mostly, but not exclusively, employing a spherical harmonic triangular truncation at wavenumber 63 (T63) or higher. The methodology that is used is spectral, and is motivated by the stochastic representation of statistical closure theory, with the ‘damping’ and forcing covariance, representing backscatter, determined from the statistics of DNS. The damping and forcing covariance are formulated as 2 × 2 matrices for each wavenumber. As well as the transient subgrid tendency, the mean subgrid tendency is needed in the LES when the energy injection region is unresolved; this is also calculated from the statistics of the DNS. For comparison, a deterministic parameterization scheme consisting of 2×2 ‘damping’ parameters, which are calculated from the statistics of DNS, has been constructed. The main difference between atmospheric and oceanic flows, in this thesis, is that the atmospheric LES completely resolves the deformation scale, the energy and enstrophy injection region, and the truncation scale is spectrally distant from it, being well in the enstrophy cascade inertial range. In oceanic flows, however, the truncation scale is in the vicinity of the injection scale, at least for the parameters chosen, and is therefore not in an inertial range. A lower resolution oceanic LES at T15 is also examined, in which case the injection region is not resolved at all.¶ For atmospheric flows, it is found that, at T63, the matrix parameters are practically diagonal so that stratified atmospheric flows at these resolutions may be treated as uncoupled layers as far as subgrid-scale parameterizations are concerned. It is also found that the damping parameters are relatively independent of the (vertical) level, but the backscatter parameters are proportional to the subgrid flux in a given level. The stochastic and deterministic parameterization schemes give comparably good results relative to the DNS. For oceanic flows, it is found that the full matrix structure of the parameters must be used. Furthermore, it is found that there is a strong injection of barotropic energy from the subgrid scales, due to the unresolved, or partially resolved, baroclinic instability injection scales. It is found that the deterministic parameterization is too numerically unstable to be of use in the LES, and instead the stochastic parameterization must be used to obtain good agreement with the DNS. The subgrid tendency of the ensemble mean flow is also needed in some problems, and is found to reduce the available potential energy of the flow.

In this study, numerical simulations of a low emission lab-scale non-premixed combustor are conducted and analyzed. The objectives are to provide new insight into the physical phenomena in the SPRF (Stagnation Point Reverse Flow) combustor built in the Georgia Tech Combustion Lab, and to compare three Large Eddy Simulation (LES) combustion models (Eddy Break-Up [EBU], Steady Flamelet [SF] and Linear Eddy Model [LEM]) for non-premixed combustion. The nominal operating condition of the SPRF combustor achieves very low NOx and CO emissions by combining turbulent mixing of exhaust gases with preheated reactants and chemical kinetics. The SPRF numerical simulation focuses on capturing the complex interaction between turbulent mixing and heat release. LES simulations have been carried out for a non-reactive case in order to analyze the turbulent mixing inside the combustor. The LES results have been compared to PIV experimental data and the code has been validated. The dominating features of the operational mode of the SPRF combustor (dilution of hot products into reactants, pre-heating and pre-mixing) have been analyzed, and results from the EBU-LES, SF-LES and LEM-LES simulations have been compared. Analysis shows that the LEM-LES simulation achieves the best agreement with the observed flame structure and is the only model that captures the stabilization processes observed in the experiments. EBU-LES and SF-LES do not predict the correct flow pattern because of the inaccurate modeling of sub-grid scale mixing and turbulence-combustion interaction. Limitations of these two models for this type of combustor are discussed.

Large Eddy Simulation (LES) of bluff body stabilized premixed and partially premixed combustion close to the flammability limit is carried out in this thesis. The LES algorithm has no ad-hoc adjustable model parameters and is able to respond automatically to variations in the inflow conditions. Algorithm validation is achieved by comparison with reactive and non-reactive experimental data. In the reactive flow, two scalar closure models, Eddy Break-Up (EBULES) and Linear Eddy Mixing (LEMLES), are used and compared. Over important regions, the flame lies in the Broken Reaction Zone regime. Here, the EBU model assumptions fail. The flame thickness predicted by LEMLES is smaller and the flame is faster to respond to turbulent fluctuations, resulting in a more significant wrinkling of the flame surface. As a result, LEMLES captures better the subtle effects of the flame-turbulence interaction. Three premixed (equivalence ratio = 0.6, 0.65, and 0.75) cases are simulated. For the leaner case, the flame temperature is lower, the heat release is reduced and vorticity is stronger. As a result, the flame in this case is found to be unstable. In the rich case, the flame temperature is higher, and the spreading rate of the wake is increased due to the higher amount of heat release Partially premixed combustion is simulated for cases where the transverse profile of the inflow equivalence ratio is variable. The simulations show that for mixtures leaner in the core the vortical pattern tends towards anti-symmetry and the heat release decreases, resulting also in instability of the flame. For mixtures richer in the core, the flame displays sinusoidal flapping resulting in larger wake spreading. More accurate predictions of flame stability will require the use of detailed chemistry, raising the computational cost of the simulation. To address this issue, a novel algorithm for training Artificial Neural Networks (ANN) for prediction of the chemical source terms has been implemented and tested. Compared to earlier methods, the main advantages of the ANN method are in CPU time and disk space and memory reduction.

Liquid fragmentation phenomenon is explored from both a fundamental (fully resolved) and an engineering (modeled) perspective. The dual objectives compliment each other by providing an avenue to gain further understanding into fundamental processes of atomization as well as to use the newly acquired knowledge to address practical concerns. A compressible five-equation interface model based on a Roe-type scheme for the simulation of material boundaries between immiscible fluids with arbitrary equation of state is developed and validated. The detailed simulation model accounts for surface-tension, viscous, and body-force effects, in addition to acoustic and convective transport. The material interfaces are considered as diffused zones and a mixture model is given for this transition region. The simulation methodology combines a high-resolution discontinuity capturing method with a low-dissipation central scheme resulting in a hybrid approach for the solution of time- and space-accurate interface problems. Several multi-dimensional test cases are considered over a wide range of physical situations involving capillary, viscosity, and gravity effects with simultaneous presence of large viscosity and density ratios. The model is shown to accurately capture interface dynamics as well as to deal with dynamic appearance and disappearance of material boundaries. Simulation of atomization processes and its interaction with the flow field in practical devices is the secondary objective of this study. Three modeling requirements are identified to perform Large-Eddy Simulation (LES) of spray combustion in engineering devices. In concurrence with these requirements, LES of an experimental liquid-fueled Lean Direct Injection (LDI) combustor is performed using a subgrid mixing and combustion model. This approach has no adjustable parameters and the entire flow-path through the inlet swirl vanes is resolved. The inclusion of the atomization aspects within LES eliminates the need to specify dispersed-phase size-velocity correlations at the inflow boundary. Kelvin-Helmholtz (or aerodynamic) breakup model by Reitz is adopted for the combustor simulation. Two simulations (with and without breakup) are performed and compared with measurements of Cai et al. Time-averaged velocity prediction comparison for both gas- and liquid-phase with available data show reasonable agreement. The major impact of breakup is on the fuel evaporation in the vicinity of the injector. Further downstream, a wide range of drop sizes are recovered by the breakup simulation and produces similar spray quality as in the no-breakup case.

No description available.

In this work, the mathematical implications of merging two different turbulence modeling approaches are addressed by deriving the exact hybrid RANS/LES Navier-Stokes equations. These equations are derived by introducing an additive-filter, which linearly combines the RANS and LES operators with a blending function. The equations derived predict additional hybrid terms, which represent the interactions between RANS and LES formulations. Theoretically, the prediction of the hybrid terms demonstrates that the hybridization of the two approaches cannot be accomplished only by the turbulence model equations, as it is claimed in current hybrid RANS/LES models. The importance of the exact hybrid RANS/LES equations is demonstrated by conducting numerical calculations on a turbulent flat-plate boundary layer. Results indicate that the hybrid terms help to maintain an equilibrated model transition when the hybrid formulation switches from RANS to LES. Results also indicate, that when the hybrid terms are not included, the accuracy of the calculations strongly relies on the blending function implemented in the additive-filter. On the other hand, if the exact equations are resolved, results are only weakly affected by the characteristics of the blending function. Unfortunately, for practical applications the hybrid terms cannot be exactly computed. Consequently, a reconstruction procedure is proposed to approximate these terms. Results show, that the model proposed is able to mimic the exact hybrid terms, enhancing the accuracy of current hybrid RANS/LES approaches. In a second effort, the Two Level Simulation (TLS) approach is proposed as a near-wall model for LES. Here, TLS is first extended to compressible flows by deriving the small-scale equations required by the model. The full compressible TLS formulation and the hybrid TLS/LES approach is validated simulating the flow over a flat-plate turbulent boundary layer. Overall, results are found in reasonable agreement with experimental data and LES calculations.

A new hybrid approach to model high Reynolds number wall-bounded turbulent flows is developed based on coupling the two-level simulation (TLS) approach in the inner region with conventional large eddy simulation (LES) away from the wall. This new approach is significantly different from previous near-wall approaches for LES. In this hybrid TLS-LES approach, a very fine small-scale (SS) mesh is embedded inside the coarse LES mesh in the near-wall region. The SS equations capture fine-scale temporal and spatial variations in all three cartesian directions for all three velocity components near the wall. The TLS-LES equations are derived based on defining a new scale separation operator. The TLS-LES equations in the transition region are obtained by blending the TLS large-scale and LES equations. A new incompressible parallel flow solver is developed that accurately and reliably predicts turbulent flows using TLS-LES. The solver uses a primitive variable formulation based on an artificial compressibility approach and a dual time stepping method. The advective terms are discretized using fourth-order energy conservative finite differences. The SS equations are also integrated in parallel, which reduces the overall cost of the TLS-LES approach. The TLS-LES approach is validated and investigated for canonical channel flows, channel flow with adverse pressure gradient and asymmetric plane diffuser flow. The results suggest that the TLS-LES approach yields very reasonable predictions of most of the crucial flow features in spite of using relatively coarse grids.

By collocating 10 years (1999-2009) of remotely sensed surface turbulent heat fluxes with satellite altimetry data, we investigate the impact of ocean mesoscale eddies on the latent and sensible heat fluxes in the South Atlantic ocean. Eddies were identified using the method proposed by Chaigneau et al. (2009), which is based on closed contours of sea level anomaly. Most of the identified eddies had a radius of ~70 km and amplitude of ~5 cm. On average, in the South Atlantic, eddies play a minor role on the ocean-atmosphere heat exchange. However, in strongly energetic regions such as the Brazil-Malvinas confluence or Agulhas Current retroflection regions, eddies can account up to 20-30% of the total variance of the surface turbulent heat fluxes with averaged anomalies of ±10-20 W/m2 for both heat flux components. Cyclonic (anticyclonic) eddies, associated with negative (positive) heat fluxes anomalies tend to cool (warm) the overlying atmosphere. A composite analysis of the turbulent heat fluxes anomalies within the eddies reveals a direct relationship between the eddy amplitude and the intensity of the latent and sensible fluxes anomalies, such that large-amplitude eddies have a stronger signature in the turbulent surface heat fluxes. Heat fluxes anomalies are also much stronger near the eddy centers and decay radially to reach minimum values outside the eddies. / Uma combinação de 10 anos (1999-2009) de fluxos turbulentos de calor pela superfície, medidos a partir de satélites, e dados altimétricos de anomalia da altura da superfície do mar, foram o utilizados com objetivo de investigar o impacto de vórtices de meso-escala nos fluxos de calor sensível e latente na bacia do Atlântico Sul. Para a detecção dos vórtices foi aplicado o método proposto por Chaigneau et al. (2009), que baseia-se em contornos fechados de anomalia da altura da superfície do mar. A maior parte dos vórtices identificados possui raio de ~70 km e amplitude de ~5 cm. Em média, no Atlântico Sul, o impacto dos vórtices para as trocas de calor entre oceano e atmosfera é relativamente fraco. Entretanto, em regiões de alta variabilidade energética como na Confluência Brasil- Malvinas e na retroflecção da Corrente das Agulhas, vórtices de meso-escala podem contribuir com anomalias médias de até ±10-20 W/m2 nos fluxos turbulentos. Vórtices ciclônicos (anti-ciclônicos), associados com anomalias negativas (positivas) de fluxos de calor, tendem a esfriar (esquentar) a atmosfera adjacente. Mapas composite foram analisados para milhares de vórtices, mostrando um relação direta entre a magnitude das anomalias dos fluxos e a amplitude dos vórtices, de tal modo que vórtices de maior amplitude contribuem com maiores anomalias de calor latente e sensível. Além disso, os padrões espaciais dos composites médios revelam que as anomalias são significativamente maiores próximo ao centro dos vórtices e decaem radialmente até atingirem valores absolutos mínimos fora dos contornos dos vórtices.