Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Direct coupling analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

Effects of Buoyancy on Laminar and Turbulent Premixed V-Flame

Turbulent combustion occurs naturally in almost all combustion systems and involves complex dynamic coupling of chemical and fluid mechanical processes. It is considered as one of the most challenging combustion research problems today. Though buoyancy has little effect on power generating systems operating under high pressures (e.g., IC engines and turbines), flames in atmospheric burners and the operation of small to medium furnaces and boilers are profoundly affected by buoyancy. Changes in burner orientation impacts on their blow-off, flash-back and extinction limits, and their range of operation, burning rate, heat transfer, and emissions. Theoretically, buoyancy is often neglected in turbulent combustion models. Yet the modeling results are routinely compared with experiments of open laboratory flames that are obviously affected by buoyancy. This inconsistency is an obstacle to reconciling experiments and theories. Consequently, a fundamental understanding of the coupling between turbulent flames and buoyancy is significant to both turbulent combustion science and applications. The overall effect of buoyancy relates to the dynamic interaction between the flame and its surrounding, i.e., the so-called elliptical problem. The overall flame shape, its flowfield, stability, and mean and local burning rates are dictated by both upstream and downstream boundary conditions. In steady propagating premixed flames, buoyancy affects the products region downstream of the flame zone. These effects are manifested upstream through the mean and fluctuating pressure fields to influence flame stretch and flame wrinkling. Intuitively, the effects buoyancy should diminish with increasing flow momentum. This is the justification for excluding buoyancy in turbulent combustion models that treats high Reynolds number flows. The objectives of our experimental research program is to elucidate flame-buoyancy coupling processes in laminar and turbulent premixed flames, and to characterize microgravity (micro g) premixed flames. The results are used to derive appropriate scaling parameters for guiding the development of theoretical models to include the effects of buoyancy. Knowledge gain from the analysis will also contribute to further understanding of the elliptical nature of premixed flames. Our current emphasis is to examine the momentum limit above which the effects of buoyancy would become insignificant. This is accomplished by comparing the flowfields and the mean properties of normal gravity flames (+g), and reversed gravity flames (-g, up-side-down flames) at different flow velocities and turbulence intensities. Microgravity (micro g) flames experiments provide the key reference data to reconcile the differences between flames in +g and -g. As flame configuration has significant impact on premixed flames characteristics we have studied axi-symmetric conical flames and plane-symmetric rod-stabilized v-flames. The two configurations produce distinct features that dictates how the flames couple with buoyancy. In a conical flame, the hot products plume completely envelopes the flame cone and shields the flame from direct interaction with the ambient air. The plume originates at the burner rim and generates a divergent flowfield. In comparison, the products region of v-flames forms between the twin flame sheets and it is convergent towards the center-plane. Interaction with ambient air is limited to the two end regions of the stabilized rod and beyond the flame sheets.

Cheng, Robert K.↗
NASA NTRS2022

Lunar Equator Regenerative Fuel Cell System Efficiency Analysis

It is expected that future NASA mission energy storage requirements can be adequately met by regenerative fuel cells (RFCs). Aerospace RFCs are likely to package proton exchange membrane electrochemical conversion devices based on hydrogen, oxygen, and water. To optimize system designs and direct further development, the most useful system metric is round-trip efficiency (RTE). Improving RTE generally increases specific energy and reduces system mass. The following analysis models the impact of several factors on RTE for RFCs scaling from 0.1 to 50 kW. The analysis demonstrates that higher electrolyzer (EZ) operational pressure most negatively impacts both RFC RTE and specific energy. Reactant storage heating and power management losses are significant for RFCs of any scale, but solenoid valve power becomes a noticeable parasitic load for smaller RFCs. Spherical gas storage vessels benefit RTE by only 2 to 4 percentage points but increase specific energy by ~50 percent. Cryogenic reactant storage is shown to depress RTE to below 15 percent in all cases. For a given RFC scale, there is an RTE benefit to specifying larger EZs and recharging at higher rates—if the RFC can be coupled with a corresponding power source. Although fuel cell (FC) operation is less impactful on RTE than EZ factors, efficiency is maximized by sizing the FCs relatively large and keeping current density low. RFC design is a balance of many interrelated factors, but this analysis provides a starting point for difficult decisions that must be made when designing a lunar equator RFC.

regenerative fuel cell↗
NASA NTRS1986

F-106 data summary and model results relative to threat criteria and protection design analysis

The NASA F-106 has acquired considerable data on the rates-of-change of electromagnetic parameters on the aircraft surface during 690 direct lightning strikes while penetrating thunderstorms at altitudes ranging from 15,000 to 40,000 feet. These in-situ measurements have provided the basis for the first statistical quantification of the lightning electromagnetic threat to aircrat appropriate for determining lightning indirect effects on aircraft. The data are presently being used in updating previous lightning criteria and standards developed over the years from ground-based measurements. The new lightning standards will, therefore, be the first which reflect actual aircraft responses measured at flight altitudes. The modeling technique developed to interpret and understand the direct strike electromagnetic data acquired on the F-106 provides a means to model the interaction of the lightning channel with the F-106. The reasonable results obtained with the model, compared to measured responses, yield confidence that the model may be credibly applied to other aircraft types and uses in the prediction of internal coupling effects in the design of lightning protection for new aircraft.

Pitts, F. L.↗
NASA NTRS1978

Infrared band intensities of saturated hydrocarbons

Kramers-Kronig analysis is applied to measured values of spectral reflectance at near-normal incidence to determine the real and the imaginary parts of the complex index of refraction for methane, ethane, propane, n-butane, n-hexane, n-heptane, and n-decane in the liquid state. The results indicate that the strengths of the characteristic bands as measured by the integral of the imaginary part are roughly constant for all the liquid alkanes except for methane. The intensity of the CH valence vibration bands in the spectra of the alkanes except methane is directly proportional to the number of CH groups per unit volume. The relations for the intensity of the bands due to CH2 and CH3 deformations are examined. Characteristic band intensities of the type established for NH4(+) and SO4(2-) groups in solutions and crystals cannot be extended to the more closely coupled CH2 and CH3 groups in alkane molecules.

Pinkley, L. W.↗
NASA NTRS1977

Design, fabrication, and evaluation of charge-coupled devices with aluminum-anodized-aluminum gates

A 4-phase, 49 1/2 bit CCD shift register was designed and fabricated using two levels of aluminum metallization with anodic Al2O3 insulation separating the layers. Test circuitry was also designed and constructed. A numerical analysis of an MOS-RC transmission line was made and results are given to characterize performance for various conductivities. The electrical design of the CCD included a low-noise dual-gate input and a balanced floating diffusion output circuit. Metallization was accomplished both by low voltage DC sputtering and thermal evaporation. The audization was according to published procedures using a buffered tartaric acid bath. Approximately 20 wafers were processed with 50 complete chips per wafer. All devices failed by shorting between the metal levels at some point. Experimental procedures eliminated temperature effects from sintering and drying, anodic oxide thickness, edge effects, photoresist stripping procedures, and metallization techniques as the primary causes of failure. It was believed from a study of SEM images that protuberances (hillocks) grow up from the first level metal through the oxide either causing a direct short or producing a weak, highly stressed insulation point which fails at low voltage. The cause of these hillocks is unknown; however, they have been observed to grow during temperature excursions to 470 C.

Gassaway, J. D.↗
NASA NTRS2007

Liquid Bismuth Propellant Flow Sensor

Quantifying the propellant mass flow rate in liquid bismuth-fed electric propulsion systems has two challenging facets. First, the flow sensors must be capable of providing a resolvable measurement at propellant mass flow rates on the order of 10 mg/see with and uncertainty of less that 5%. The second challenge has to do with the fact that the materials from which the flow sensors are fabricated must be capable of resisting any of the corrosive effects associated with the high-temperature propellant. The measurement itself is necessary in order to properly assess the performance (thrust efficiency, Isp) of thruster systems in the laboratory environment. The hotspot sensor[I] has been designed to provide the bismuth propellant mass flow rate measurement. In the hotspot sensor, a pulse of thermal energy (derived from a current pulse and associated joule heating) is applied near the inlet of the sensor. The flow is "tagged" with a thermal feature that is convected downstream by the flowing liquid metal. Downstream, a temperature measurement is performed to detect a "ripple" in the local temperature associated with the passing "hotspot" in the propellant. By measuring the time between the upstream generation and downstream detection of the thermal feature, the flow speed can be calculated using a "time of flight" analysis. In addition, the system can be calibrated by measuring the accumulated mass exiting the system as a-function of time and correlating this with the time it takes the hotspot to convect through the sensor. The primary advantage of this technique is that it doesn't depend on an absolute measurement of temperature but, instead, relies on the observation of thermal features. This makes the technique insensitive to other externally generated thermal fluctuations. In this paper, we describe experiments performed using the hotspot flow sensor aimed at quantifying the resolution of the sensor technology. Propellant is expelled onto an electronic scale to provide an independent measure of the propellant mass flow rate as a function of time. In addition, two separate detection schemes are employed. The first uses a thermocouple to directly measure temperature in the fluid. The second involves the ,use of a fiber optic coupled to a photodiode allowing for detection of an increase in light emission from the fluid as the hotspot passes. the detection location.

Polzin, Kurt A.↗
NASA NTRS1987

Analysis of a parallelized nonlinear elliptic boundary value problem solver with application to reacting flows

A parallelized finite difference code based on the Newton method for systems of nonlinear elliptic boundary value problems in two dimensions is analyzed in terms of computational complexity and parallel efficiency. An approximate cost function depending on 15 dimensionless parameters is derived for algorithms based on stripwise and boxwise decompositions of the domain and a one-to-one assignment of the strip or box subdomains to processors. The sensitivity of the cost functions to the parameters is explored in regions of parameter space corresponding to model small-order systems with inexpensive function evaluations and also a coupled system of nineteen equations with very expensive function evaluations. The algorithm was implemented on the Intel Hypercube, and some experimental results for the model problems with stripwise decompositions are presented and compared with the theory. In the context of computational combustion problems, multiprocessors of either message-passing or shared-memory type may be employed with stripwise decompositions to realize speedup of O(n), where n is mesh resolution in one direction, for reasonable n.

Keyes, David E.↗
NASA NTRS2004

Large-Scale Alfvenic Impulses on the Sun: How They Are Generated and What We Learn From Them

NASA GSFC The Sun's atmosphere hosts a wide variety of magnetosonic disturbances. These wave modes are detected, almost exclusively, by examining images of the Sun's magnetic atmosphere and looking for propagating distortions. Although none of the Sun's plasma parameters are measured directly, we derive a great deal of information from these observations. In fact, by modeling these propagating disturbances, we may be able to derive the most accurate estimates plasma parameters. From observations absorption, refraction, reflection, and coupling of numerous wave modes, we advance our knowledge of the Sun's magnetic field, temperature, density, and current. The Sun's continuous oscillation, coronal mass ejections, flares, and other dynamic phenomena can produce wave disturbances which are observable from near-Earth space. Several of these disturbances have been traced from the inner corona out into the heliosphere. From the generation of these disturbances, we are able to learn about the phenomena which create them as well as the media through which they re-propagating. The presentation will include a discussion of the generation of Alfvenic disturbances on the Sun, ways we observe these disturbances, and how recent advances in modeling and analysis have brought us closer to determining solar in situ parameters.

Thompson, Barbara↗
NASA NTRS1996

A General Interface Method for Aeroelastic Analysis of Aircraft

The aeroelastic analysis of an aircraft requires an accurate and efficient procedure to couple aerodynamics and structures. The procedure needs an interface method to bridge the gap between the aerodynamic and structural models in order to transform loads and displacements. Such an interface method is described in this report. This interface method transforms loads computed by any aerodynamic code to a structural finite element (FE) model and converts the displacements from the FE model to the aerodynamic model. The approach is based on FE technology in which virtual work is employed to transform the aerodynamic pressures into FE nodal forces. The displacements at the FE nodes are then converted back to aerodynamic grid points on the aircraft surface through the reciprocal theorem in structural engineering. The method allows both high and crude fidelities of both models and does not require an intermediate modeling. In addition, the method performs the conversion of loads and displacements directly between individual aerodynamic grid point and its corresponding structural finite element and, hence, is very efficient for large aircraft models. This report also describes the application of this aero-structure interface method to a simple wing and an MD-90 wing. The results show that the aeroelastic effect is very important. For the simple wing, both linear and nonlinear approaches are used. In the linear approach, the deformation of the structural model is considered small, and the loads from the deformed aerodynamic model are applied to the original geometry of the structure. In the nonlinear approach, the geometry of the structure and its stiffness matrix are updated in every iteration and the increments of loads from the previous iteration are applied to the new structural geometry in order to compute the displacement increments. Additional studies to apply the aero-structure interaction procedure to more complicated geometry will be conducted in the second phase of the present contract.

Tzong, T.↗
NASA NTRS1992

Stability of generally stiffened anisotropic noncircular cylinders

Continuous filament grid-stiffened structure is a stiffening concept that combines structural efficiency and damage tolerance. However, finite element design of such structures against buckling is expensive due to the complexities of the structure. An analytical model of such a structure is developed using a penalty method (artificial springs) with a first order shear deformation theory (FSDT). The buckling analysis under combined loadings is done using energy method with a penalty/Rayleigh-Ritz technique. The penalty/Rayleigh-Ritz approach is computationally less demanding when compared to the finite element solution and mesh generation. Apart from the published research works on buckling of stiffened plates and shells by finite element and finite strips, research works on buckling of stiffened plates and shells utilize three different approaches; smeared, column, and discrete approaches. The discrete approach considers the discrete effects of the stiffeners in the buckling behavior by modeling stiffeners as line of bending (EI) and torsion (GJ) stiffnesses on panel skin. Some local deformations are lost when stiffeners are modeled as (EI) and (GJ) stiffeners. This approach becomes difficult in the case of plate stiffened in more than two directions. Most of the work done using the discrete approach involved the Classical Plate Theory (CLPT) rather than the FSDT. We report on our formulation of a discrete approach coupled with a penalty formulation and FSDT.

Sobh, Nahil Atef↗
NASA NTRS2013

The Evolution of Software and Its Impact on Complex System Design in Robotic Spacecraft Embedded Systems

The growth in computer hardware performance, coupled with reduced energy requirements, has led to a rapid expansion of the resources available to software systems, driving them towards greater logical abstraction, flexibility, and complexity. This shift in focus from compacting functionality into a limited field towards developing layered, multi-state architectures in a grand field has both driven and been driven by the history of embedded processor design in the robotic spacecraft industry.The combinatorial growth of interprocess conditions is accompanied by benefits (concurrent development, situational autonomy, and evolution of goals) and drawbacks (late integration, non-deterministic interactions, and multifaceted anomalies) in achieving mission success, as illustrated by the case of the Mars Reconnaissance Orbiter. Approaches to optimizing the benefits while mitigating the drawbacks have taken the form of the formalization of requirements, modular design practices, extensive system simulation, and spacecraft data trend analysis. The growth of hardware capability and software complexity can be expected to continue, with future directions including stackable commodity subsystems, computer-generated algorithms, runtime reconfigurable processors, and greater autonomy.

software↗
NASA NTRS2016

Analysis and Characterization of Damage and Failure Utilizing a Generalized Composite Material Model Suitable for Use in Impact Problems

A material model which incorporates several key capabilities which have been identified by the aerospace community as lacking in state-of-the art composite impact models is under development. In particular, a next generation composite impact material model, jointly developed by the FAA and NASA, is being implemented into the commercial transient dynamic finite element code LS-DYNA. The material model, which incorporates plasticity, damage, and failure, utilizes experimentally based tabulated input to define the evolution of plasticity and damage and the initiation of failure as opposed to specifying discrete input parameters (such as modulus and strength). The plasticity portion of the orthotropic, three-dimensional, macroscopic composite constitutive model is based on an extension of the Tsai-Wu composite failure model into a generalized yield function with a non-associative flow rule. For the damage model, a strain equivalent formulation is utilized to allow for the uncoupling of the deformation and damage analyses. In the damage model, a semi-coupled approach is employed where the overall damage in a particular coordinate direction is assumed to be a multiplicative combination of the damage in that direction resulting from the applied loads in the various coordinate directions. Due to the fact that the plasticity and damage models are uncoupled, test procedures and methods to both characterize the damage model and to covert the material stress-strain curves from the true (damaged) stress space to the effective (undamaged) stress space have been developed. A methodology has been developed to input the experimentally determined composite failure surface in a tabulated manner. An analytical approach is then utilized to track how close the current stress state is to the failure surface.

Finite Element Method↗
NASA NTRS2012

Analysis of Organic Molecules Extracted from Mars Analogues and Influence of Their Mineralogy Using N-Methyl-N-(tert-butyldimethylsilyl)Trifluoroacetamide Derivatization Coupled with Gas Chromatography Mass Spectrometry in Preparation for the Sample Analysis at Mars Derivatization Experiment on the Mars Science Laboratory Mission

The search for complex organic molecules on Mars, including important biomolecules such as amino acids and carboxylic acids will require a chemical extraction and derivatization step to transform these organic compounds into species that are sufficiently volatile to be detected by gas chromatography mass spectrometry (GCMS). We have developed, a one-pot extraction and chemical derivatization protocol using N-methyl-N-(tert-butyldimethylsilyl)trifluoroacetamide (MTBSTFA) and dimethylformamide (DMF) for the Sample Analysis at Mars (SAM) experiment on the Mars Science Laboratory (MSL). The temperature and duration the derivatization reaction, pre-concentration of chemical derivatives, and gas chromatographic separation parameters have been optimized under SAM instrument design constraints. MTBSTFA/DMF extraction and derivatization at 300 C for several minutes of a variety of terrestrial Mars analogue materials facilitated the detection of amino acids and carboxylic acids in a surface soil sample collected from the Atacama Desert and a carbonate-rich stromatolite sample from Svalbard. However, the rapid reaction of MTBSTFA with water in several analogue materials that contained high abundances of hydrated minerals and the possible deactivation of derivatized compounds by iron oxides, as detected by XRD/XRF using the CheMin field unit Terra, proved to be highly problematic for the direct extraction of organics using MTBSTFA, The combination of pyrolysis and two different chemical derivatization methods employed by SAM should enable a wide range of organic compounds to be detected by GCMS if present on Mars,

Stalport, F.↗
NASA NTRS1993

A comparative study of multivariable robustness analysis methods as applied to integrated flight and propulsion control

Three multivariable robustness analysis methods are compared and contrasted. The focus of the analysis is on system stability and performance robustness to uncertainty in the coupling dynamics between two interacting subsystems. Of particular interest is interacting airframe and engine subsystems, and an example airframe/engine vehicle configuration is utilized in the demonstration of these approaches. The singular value (SV) and structured singular value (SSV) analysis methods are compared to a method especially well suited for analysis of robustness to uncertainties in subsystem interactions. This approach is referred to here as the interacting subsystem (IS) analysis method. This method has been used previously to analyze airframe/engine systems, emphasizing the study of stability robustness. However, performance robustness is also investigated here, and a new measure of allowable uncertainty for acceptable performance robustness is introduced. The IS methodology does not require plant uncertainty models to measure the robustness of the system, and is shown to yield valuable information regarding the effects of subsystem interactions. In contrast, the SV and SSV methods allow for the evaluation of the robustness of the system to particular models of uncertainty, and do not directly indicate how the airframe (engine) subsystem interacts with the engine (airframe) subsystem.

Schierman, John D.↗
NASA NTRS2014

Line Interference Effects Using a Refined Robert-Bonamy Formalism: the Test Case of the Isotropic Raman Spectra of Autoperturbed N2

A symmetrized version of the recently developed refined Robert-Bonamy formalism [Q. Ma, C. Boulet, and R. H. Tipping, J. Chem. Phys. 139, 034305 (2013)] is proposed. This model takes into account line coupling effects and hence allows the calculation of the off-diagonal elements of the relaxation matrix, without neglecting the rotational structure of the perturbing molecule. The formalism is applied to the isotropic Raman spectra of autoperturbed N2 for which a benchmark quantum relaxation matrix has recently been proposed. The consequences of the classical path approximation are carefully analyzed. Methods correcting for effects of inelasticity are considered. While in the right direction, these corrections appear to be too crude to provide off diagonal elements which would yield, via the sum rule, diagonal elements in good agreement with the quantum results. In order to overcome this difficulty, a re-normalization procedure is applied, which ensures that the off-diagonal elements do lead to the exact quantum diagonal elements. The agreement between the (re-normalized) semi-classical and quantum relaxation matrices is excellent, at least for the Raman spectra of N2, opening the way to the analysis of more complex molecular systems.

nitrogen↗
NASA NTRS2001

Wind-Stress Simulations and Equatorial Dynamics in an AGCM: Basic Results from a 1979-1999 Forced SST Experiment - Part 1

This is the first of a two part study examining the connection of the equatorial momentum budget in an AGCM (Atmospheric General Circulation Model), with simulated equatorial surface wind stresses over the Pacific. The AGCM used in this study forms part of a newly developed coupled forecasting system used at NASA's Seasonal- to-Interannual Prediction Project. Here we describe the model and present results from a 20-year (1979-1999) AMIP-type experiment forced with observed SSTs (Sea Surface Temperatures). Model results are compared them with available observational data sets. The climatological pattern of extra-tropical planetary waves as well as their ENSO-related variability is found to agree quite well with re-analysis estimates. The model's surface wind stress is examined in detail, and reveals a reasonable overall simulation of seasonal interannual variability, as well as seasonal mean distributions. However, an excessive annual oscillation in wind stress over the equatorial central Pacific is found. We examine the model's divergent circulation over the tropical Pacific and compare it with estimates based on re-analysis data. These comparisons are generally good, but reveal excessive upper-level convergence in the central Pacific. In Part II of this study a direct examination of individual terms in the AGCM's momentum budget is presented. We relate the results of this analysis to the model's simulation of surface wind stress.

Bacmeister, Julio T.↗
NASA NTRS2012

Simulation of Oxygen Disintegration and Mixing With Hydrogen or Helium at Supercritical Pressure

The simulation of high-pressure turbulent flows, where the pressure, p, is larger than the critical value, p(sub c), for the species under consideration, is relevant to a wide array of propulsion systems, e.g. gas turbine, diesel, and liquid rocket engines. Most turbulence models, however, have been developed for atmospheric-p turbulent flows. The difference between atmospheric-p and supercritical-p turbulence is that, in the former situation, the coupling between dynamics and thermodynamics is moderate to negligible, but for the latter it is very significant, and can dominate the flow characteristics. The reason for this stems from the mathematical form of the equation of state (EOS), which is the perfect-gas EOS in the former case, and the real-gas EOS in the latter case. For flows at supercritical pressure, p, the large eddy simulation (LES) equations consist of the differential conservation equations coupled with a real-gas EOS. The equations use transport properties that depend on the thermodynamic variables. Compared to previous LES models, the differential equations contain not only the subgrid scale (SGS) fluxes, but also new SGS terms, each denoted as a correction. These additional terms, typically assumed null for atmospheric pressure flows, stem from filtering the differential governing equations, and represent differences between a filtered term and the same term computed as a function of the filtered flow field. In particular, the energy equation contains a heat-flux correction (q-correction) that is the difference between the filtered divergence of the heat flux and the divergence of the heat flux computed as a function of the filtered flow field. In a previous study, there was only partial success in modeling the q-correction term, but in this innovation, success has been achieved by using a different modeling approach. This analysis, based on a temporal mixing layer Direct Numerical Simulation database, shows that the focus in modeling the q-correction should be on reconstructing the primitive variable gradients rather than their coefficients, and proposes the approximate deconvolution model (ADM) as an effective means of flow field reconstruction for LES heat flux calculation. Further, results for a study conducted for temporal mixing layers initially containing oxygen in the lower stream, and hydrogen or helium in the upper stream, show that, for any LES, including SGS-flux models (constant-coefficient Gradient or Scale-Similarity models, dynamic-coefficient Smagorinsky/Yoshizawa or mixed Smagorinsky/Yoshizawa/Gradient models), the inclusion of the q-correction in the LES leads to the theoretical maximum reduction of the SGS heat-flux difference. The remaining error in modeling this new subgrid term is thus irreducible.

Bellan, Josette↗
NASA NTRS2017

Influence of Mid-Latitude Cyclones on European Background Surface Ozone Investigated in Observations, MACC and MERRA-2 Reanalyses

The relationship between springtime mid-latitude cyclones and background ozone at two rural monitoring sites on the west coast of Europe -- Mace Head, Ireland and Monte Velho, Portugal -- is explored using a combination of observations and three reanalyses: 1) the European Centre for Medium-Range Weather Forecasts' (ECMWF) ERA-Interim reanalysis, 2) the Monitoring Atmospheric Composition and Climate (MACC) reanalysis and 3) NASA's Modern-Era Retrospective Analysis for Research and Applications Version-2 (MERRA-2) reanalysis. The ERA-Interim cyclone tracks are used here to establish the long-term relationship between cyclones and ozone observations (since 1988). The MACC reanalysis data set, which covers the period 2003-2012, is produced with the ECMWF integrated forecast system (IFS) model two-way coupled to a chemistry transport model (CTM). Since the MACC reanalysis uses a similar atmospheric model to ERA-Interim, MACC is used to explore the mechanisms within the case study cyclones that can influence surface ozone concentrations at Mace Head and Monte Velho. The MERRA-2 reanalysis also provides 3D distributions of ozone, although less ideal for analysis of surface ozone concentrations since MERRA-2 ozone under represents ozone variability outside the stratosphere as it does not have a detailed chemistry scheme or emission sources for the troposphere. The MERRA-2 reanalysis, which has the potential to identify more features within the cyclones as the resolution is higher than the MACC reanalysis, is used in conjunction with the MACC reanalysis to provide a measure of uncertainty to the case study analysis. We found the main source of high ozone to these two sites is from the stratosphere, which is well represented in both the MERRA-2 and the MACC reanalyses, either from direct injection into the cyclone or associated with aged airstreams from decaying downstream cyclones that can become entrained and descend toward the surface within new cyclones over the NA region.

MERRA-2↗