Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “partial correlation”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Evaluation of Composite Airframe Dynamic Impact Modeling Using Hawker 4000 Fuselage Drop Test Data

Two drop tests of partial Hawker 4000 fuselage sections were conducted at the National Aeronautics and Space Administration (NASA) Langley Research Center (LaRC) to characterize the response of representative composite aerospace structure to dynamic impact loads. Test conditions were selected to induce damage into the composite structure in order to study material failure within a composite fuselage and evaluate the capability of finite element (FE) model analysis to predict that failure. The tests were simulated using FE models which were generated to isolate the effect of developmental data availability on predictive capability. FE Models of the tested fuselage sections were generated using two limited data sets. The first model configuration was reverse engineered from the test article with no information related to design or fabrication details which would be known only by the manufacturer. The second model was generated from data provided by the manufacturer but without additional material characterization test data. Models were developed using these methodologies for both fuselage sections tested. Correlation of each model to the tests conducted was evaluated in terms of damage, deformation, and cabin acceleration measurements. Correlation between the developed models and the tested fuselage sections showed that the reverse engineered model predicted the composite damage and cabin acceleration measured during test though it was limited due to lack of detail in the composite layup changes through the structure. The model developed using manufacturer specifications did not predict damage, due to limited material and component model characterization data, but it did predict acceleration on par with the reverse engineered model. Model capability and limitation sources identified were verified through correlation of a final model which was developed by combining the individual data sets. The combined model demonstrated that the addition of calibrated composite material models to accurate composite layup definitions and detailed geometry led to improved correlation of damage and acceleration response within the composite fuselage structures.

Crashworthiness↗
NASA NTRS2001

Water Vapor Tracers as Diagnostics of the Regional Hydrologic Cycle

Numerous studies suggest that local feedback of surface evaporation on precipitation, or recycling, is a significant source of water for precipitation. Quantitative results on the exact amount of recycling have been difficult to obtain in view of the inherent limitations of diagnostic recycling calculations. The current study describes a calculation of the amount of local and remote geographic sources of surface evaporation for precipitation, based on the implementation of three-dimensional constituent tracers of regional water vapor sources (termed water vapor tracers, WVT) in a general circulation model. The major limitation on the accuracy of the recycling estimates is the veracity of the numerically simulated hydrological cycle, though we note that this approach can also be implemented within the context of a data assimilation system. In the WVT approach, each tracer is associated with an evaporative source region for a prognostic three-dimensional variable that represents a partial amount of the total atmospheric water vapor. The physical processes that act on a WVT are determined in proportion to those that act on the model's prognostic water vapor. In this way, the local and remote sources of water for precipitation can be predicted within the model simulation, and can be validated against the model's prognostic water vapor. As a demonstration of the method, the regional hydrologic cycles for North America and India are evaluated for six summers (June, July and August) of model simulation. More than 50% of the precipitation in the Midwestern United States came from continental regional sources, and the local source was the largest of the regional tracers (14%). The Gulf of Mexico and Atlantic regions contributed 18% of the water for Midwestern precipitation, but further analysis suggests that the greater region of the Tropical Atlantic Ocean may also contribute significantly. In most North American continental regions, the local source of precipitation is correlated with total precipitation. There is a general positive correlation between local evaporation and local precipitation, but it can be weaker because large evaporation can occur when precipitation is inhibited. In India, the local source of precipitation is a small percentage of the precipitation owing to the dominance of the atmospheric transport of oceanic water. The southern Indian Ocean provides a key source of water for both the Indian continent and the Sahelian region.

Bosilovich, Michael G.↗
NASA NTRS1992

The structure and energetics of Cr(CO)6 and Cr(CO)5

The geometric structure of Cr(CO)6 is optimized at the modified coupled pair functional (MCPF), single and double excitation coupled-cluster (CCSD) and CCSD(T) levels of theory (including a perturbational estimate for connected triple excitations), and the force constants for the totally symmetric representation are determined. The geometry of Cr(CO)5 is partially optimized at the MCPF, CCSD, and CCSD(T) levels of theory. Comparison with experimental data shows that the CCSD(T) method gives the best results for the structures and force constants, and that remaining errors are probably due to deficiencies in the one-particle basis sets used for CO. The total binding energies of Cr(CO)6 and Cr(CO)5 are also determined at the MCPF, CCSD, and CCSD(T) levels of theory. The CCSD(T) method gives a much larger total binding energy than either the MCPF or CCSD methods. An analysis of the basis set superposition error (BSSE) at the MCPF level of treatment points out limitations in the one-particle basis used. Calculations using larger basis sets reduce the BSSE, but the total binding energy of Cr(CO)6 is still significantly smaller than the experimental value, although the first CO bond dissociation energy of Cr(CO)6 is well described. An investigation of 3s3p correlation reveals only a small effect. In the largest basis set, the total CO binding energy of Cr(CO)6 is estimated to be 140 kcal/mol at the CCSD(T) level of theory, or about 86 percent of the experimental value. The remaining discrepancy between the experimental and theoretical value is probably due to limitations in the one-particle basis, rather than limitations in the correlation treatment. In particular an additional d function and an f function on each C and O are needed to obtain quantitative results. This is underscored by the fact that even using a very large primitive set (1042 primitive functions contracted to 300 basis functions), the superposition error for the total binding energy of Cr(CO)6 is 22 kcal/mol at the MCPF level of treatment.

Barnes, Leslie A.↗
NASA NTRS1992

First observations of mesospheric dynamics with a partial reflection radar in Hawaii (22 deg N, 160 deg W)

A partial-reflection MF radar was installed at the Pacific Missile Range Facility on Kanai, Hawaii in September 1990. The wind measurements obtained with this radar system in Hawaii during the first year of operation are reviewed. The data reveal a broad spectrum of motions ranging from a mix of the equatorial mesopause semiannual oscillation (MSAO) and the annual mean cycle at middle and high latitudes in the zonal mean structure to low-frequency, tidal, and gravity wave motions at higher frequencies. The zonal mean wind structure is characterized by a downward progression of strong eastward and westward phases of the MSAO from approximately January to July. An eastward maximum of about 60/ms near 80 km during January and February is found, which descends rapidly, and a westward maximum of about 50 ms near 85 km during March and April is found, which descends much more slowly. The second MSAO cycle is greatly suppressed relative to the first due to the reversal of the correlation between this and the annual cycle at higher latitudes from July to December.

Fritts, David C.↗
NASA NTRS1994

Photosynthetic fractionation of 13C and concentrations of dissolved CO2 in the central equatorial Pacific during the last 255,000 years

Carbon isotopically based estimates of CO2 levels have been generated from a record of the photosynthetic fractionation of 13C [is equivalent to epsilon(p)] in a central equatorial Pacific sediment core that spans the last approximately 255 ka. Contents of 13C in phytoplanktonic biomass were determined by analysis of C37 alkadienones. These compounds are exclusive products of Prymnesiophyte algae which at present grow most abundantly at depths of 70-90 m in the central equatorial Pacific. A record of the isotopic composition of dissolved CO2 was constructed from isotopic analyses of the planktonic foraminifera Neogloboquadrina dutertrei, which calcifies at 70-90 m in the same region. Values of epsilon(p), derived by comparison of the organic and inorganic delta values, were transformed to yield concentrations of dissolved CO2 [is equivalent to c(e)] based on a new, site-specific calibration of the relationship between epsilon(p) and c(e). The calibration was based on reassessment of existing epsilon(p) versus c(e) data, which support a physiologically based model in which epsilon(p) is inversely related to c(e). Values of PCO2, the partial pressure of CO2 that would be in equilibrium with the estimated concentrations of dissolved CO2, were calculated using Henry's law and the temperature determined from the alkenone-unsaturation index U(K/37). Uncertainties in these values arise mainly from uncertainties about the appropriateness (particularly over time) of the site-specific relationship between epsilon(p) and 1/c(e). These are discussed in detail and it is concluded that the observed record of epsilon(p) most probably reflects significant variations in delta pCO2, the ocean-atmosphere disequilibrium, which appears to have ranged from approximately 110 microatmospheres during glacial intervals (ocean > atmosphere) to approximately 60 microatmospheres during interglacials. Fluxes of CO2 to the atmosphere would thus have been significantly larger during glacial intervals. If this were characteristic of large areas of the equatorial Pacific, then greater glacial sinks for the equatorially evaded CO2 must have existed elsewhere. Statistical analysis of air-sea pCO2 differences and other parameters revealed significant (p<0.01) inverse correlations of delta pCO2 with sea surface temperature and with the mass accumulation rate of opal. The former suggests response to the strength of upwelling, the latter may indicate either drawdown of CO2 by siliceous phytoplankton or variation of [CO2]/[Si(OH)4] ratios in upwelling waters.

NASA Discipline Exobiology↗
NASA NTRS2013

Timescale of Petrogenetic Processes Recorded in the Mount Perkins Magma System, Northern Colorado River Extension Corridor, Arizona

The Miocene Mt. Perkins Pluton is a small composite intrusive body emplaced in the shallow crust as four separate phases during the earliest stages of crustal extension. Phase 1 (oldest) consists of isotropic hornblende gabbro and a layered cumulate sequence. Phase 2 consists of quartz monzonite to quartz monzodiorite hosting mafic microgranitoid enclaves. Phase 3 is composed of quartz monzonite and is subdivided into mafic enclave-rich zones and enclave-free zones. Phase 4 consists of aphanitic dikes of mafic, intermediate and felsic compositions hosting mafic enclaves. Phases 2-4 enclaves record significant isotopic disequilibrium with surrounding granitoid host rocks, but collectively enclaves and host rocks form a cogenetic suite exhibiting systematic variations in Nd-Sr-Pb isotopes that correlate with major and trace elements. Phases 2-4 record multiple episodes of magma mingling among cogenetic hybrid magmas that formed via magma mixing and fractional crystallization at a deeper crustal. The mafic end-member was alkali basalt similar to nearby 6-4 Ma basalt with enriched OIB-like trace elements and Nd-Sr-Pb isotopes. The felsic end-member was a subalkaline crustal-derived magma. Phase 1 isotropic gabbro exhibits elemental and isotopic compositional variations at relatively constant SiO2, suggesting generation of isotropic gabbro by an open-system process involving two mafic end-members. One end-member is similar in composition to the OIB-like mafic end-member for phases 2-4; the second is similar to nearby 11-8 Ma tholeiite basalt exhibiting low epsilon (sub Nd), and depleted incompatible trace elements. Phase 1 cumulates record in situ fractional crystallization of an OIB-like mafic magma with isotopic evidence of crustal contamination by partial melts generated in adjacent Proterozoic gneiss. The Mt Perkins pluton records a complex history in a lithospheric scale magma system involving two distinct mantle-derived mafic magmas and felsic magma sourced in the crust. Mixing and fractional crystallization of these magmas at various levels in the crust generated a suite of intermediate composition magmas. U-Pb zircon SHRIMP ages of phase 1 (15.7 +/- 0.2 Ma), phase 3 (15.8 +/- 0.2 Ma) and phase 4 (15.4 +/- 0.3 Ma) document a 100-300k year timescale for petrogenetic processes recorded in the Mt Perkins magma system.

Danielson, Lisa R.↗
NASA NTRS2014

Transitional Disks Associated with Intermediate-Mass Stars: Results of the SEEDS YSO Survey

Protoplanetary disks are where planets form, grow, and migrate to produce the diversity of exoplanet systems we observe in mature systems. Disks where this process has advanced to the stage of gap opening, and in some cases central cavity formation, have been termed pre-transitional and transitional disks in the hope that they represent intermediate steps toward planetary system formation. Recent reviews have focussed on disks where the star is of solar or sub-solar mass. In contrast to the sub-millimeter where cleared central cavities predominate, at H-band some T Tauri star transitional disks resemble primordial disks in having no indication of clearing, some show a break in the radial surface brightness profile at the inner edge of the outer disk, while others have partially to fully cleared gaps or central cavities. Recently, the Meeus Group I Herbig stars, intermediate-mass PMS stars with IR spectral energy distributions often interpreted as flared disks, have been proposed to have transitional and pre-transitional disks similar to those associated with solar-mass PMS stars, based on thermal-IR imaging, and sub-millimeter interferometry. We have investigated their appearance in scattered light as part of the Strategic Exploration of Exoplanets and Disks with Subaru (SEEDS), obtaining H-band polarimetric imagery of 10 intermediate-mass stars with Meeus Group I disks. Augmented by other disks with imagery in the literature, the sample is now sufficiently large to explore how these disks are similar to and differ from T Tauri star disks. The disk morphologies seen in the Tauri disks are also found for the intermediate-mass star disks, but additional phenomena are found; a hallmark of these disks is remarkable individuality and diversity which does not simply correlate with disk mass or stellar properties, including age, including spiral arms in remnant envelopes, arms in the disk, asymmetrically and potentially variably shadowed outer disks, gaps, and one disk where only half of the disk is seen in scattered light at H. We will discuss our survey results in terms of spiral arm theory, dust trapping vortices, and systematic differences in the relative scale height of these disks compared to those around Solar-mass stars. For the disks with spiral arms we discuss the planet-hosting potential, and limits on where giant planets can be located. We also discuss the implications for imaging with extreme adaptive optics instruments. Grady is supported under NSF AST 1008440 and through the NASA Origins of Solar Systems program on NNG13PB64P. JPW is supported NSF AST 100314. 0) in marked contrast to protoplanetary disks, transitional disks exhibit wide range of structural features1) arm visibility correlated with relative scale height in disk2) asymmetric and possibly variable shadowing of outer portions some transitional disks3) confirm pre-transitional disk nature of Oph IRS 48, MWC 758, HD 169142, etc.

Protoplanetary disks↗
NASA NTRS1993

Response of the topography and gravity field on Venus to mantle upwelling beneath a chemical boundary layer

The long wavelength correlation of the gravity and topography and the large apparent depths of compensation (approximately 150-300 km) for large highland regions on Venus require significant differences between the interior structure of Earth and Venus. The morphology, geologic history, and large apparent depths of compensation for many highlands have been interpreted to indicate areas of mantle upwelling. A large apparent depth of compensation at a mantle upwelling is generally interpreted to indicate the base of the thermal boundary layer of convection. A boundary layer thickness of 150-300 km implies that the interior of Venus is presently much colder than Earth and thus tectonically less active. The recent Magellan mission has provided contradictory evidence regarding the present level of tectonic activity on Venus, prompting considerable debate. In this study, we investigate the possibility that a chemical boundary layer acts together with a thermal boundary layer to produce large apparent depths of compensation, or equivalently, large geoid-to-topography-ratios (GTR's). The crust of a planet forms through partial melting of mantle materials. Both the melt and the residuum are lower in density than unmelted (or undepleted) mantle. In the absence of vigorous plate tectonics, a thick layer of buoyant residuum, or depleted mantle, may collect beneath the lithosphere. In this scenario, the thermal lithosphere does not need to be thick and cold to match the GTR's. Cooling of the depleted layer may lead to overturn of the upper mantle and episodic resurfacing with time scales on the order of 300-500 MY, consistent with the resurfacing age of Venus.

Smrekar, Suzanne E.↗
NASA NTRS2004

Effect of high-dose irradiation on the optically stimulated luminescence of Al2O3:C

This paper examines the effect of high-dose irradiation on the optically stimulated luminescence (OSL) of Al2O3:C, principally on the shape of the OSL decay curve and on the OSL sensitivity. The effect of the degree of deep trap filling on the OSL was also studied by monitoring the sensitivity changes after doses of beta irradiation and after step-annealing of samples previously irradiated with high doses. The OSL response to dose shows a linear-supralinear-saturation behavior, with a decrease in the response for doses higher than those required for saturation. This behavior correlates with the sensitivity changes observed in the samples annealed only to 773 K, which show sensitization for doses up to 20-50 Gy and desensitization for higher doses. Data from the step-annealing study leads to the suggestion that the sensitization is caused by the filling of deep electron traps, which become thermally unstable at 1100-1200 K, whereas the desensitization is caused by the filling of deep hole traps, which become thermally unstable at 800-875 K, along with a concomitant decrease in the concentration of recombination centers (F+ -centers). Changes in the shape of the OSL decay curves are also observed at high doses, the decay becoming faster as the dose increases. These changes in the OSL decay curves are discussed in terms of multiple overlapping components, each characterized by different photoionization cross-sections. However, using numerical solutions of the rate equations for a simple model consisting of a main trap and a recombination center, it is shown that the kinetics of OSL process may also be partially responsible for the changes in the OSL curves at high doses in Al2O3:C. Finally, the implication of these results for the dosimetry of heavy charged particles is discussed. c2004 Elsevier Ltd. All rights reserved.

Heavy Ions↗
NASA NTRS2007

The Radiative Consistency of Atmospheric Infrared Sounder and Moderate Resolution Imaging Spectroradiometer Cloud Retrievals

The consistency of cloud top temperature (Tc) and effective cloud fraction (f) retrieved by the Atmospheric Infrared Sounder (AIRS)/Advanced Microwave Sounding Unit (AMSU) observation suite and the Moderate Resolution Imaging Spectroradiometer (MODIS) on the EOS-Aqua platform are investigated. Collocated AIRS and MODIS TC and f are compared via an 'effective scene brightness temperature' (Tb,e). Tb,e is calculated with partial field of view (FOV) contributions from TC and surface temperature (TS), weighted by f and 1-f, respectively. AIRS reports up to two cloud layers while MODIS reports up to one. However, MODIS reports TC, TS, and f at a higher spatial resolution than AIRS. As a result, pixel-scale comparisons of TC and f are difficult to interpret, demonstrating the need for alternatives such as Tb,e. AIRS-MODIS Tb,e differences ((Delta)Tb,e) for identical observing scenes are useful as a diagnostic for cloud quantity comparisons. The smallest values of DTb,e are for high and opaque clouds, with increasing scatter in (Delta)Tb,e for clouds of smaller opacity and lower altitude. A persistent positive bias in DTb,e is observed in warmer and low-latitude scenes, characterized by a mixture of MODIS CO2 slicing and 11-mm window retrievals. These scenes contain heterogeneous cloud cover, including mixtures of multilayered cloudiness and misplaced MODIS cloud top pressure. The spatial patterns of (Delta)Tb,e are systematic and do not correlate well with collocated AIRS-MODIS radiance differences, which are more random in nature and smaller in magnitude than (Delta)Tb,e. This suggests that the observed inconsistencies in AIRS and MODIS cloud fields are dominated by retrieval algorithm differences, instead of differences in the observed radiances. The results presented here have implications for the validation of cloudy satellite retrieval algorithms, and use of cloud products in quantitative analyses.

clouds↗
NASA NTRS2003

The Magnetic Helicity Budget of Solar Active Regions and Coronal Mass Ejections

We compute the magnetic helicity injected by transient photospheric horizontal flows in six solar active regions associated with halo coronal mass ejections (CMEs) that produced major geomagnetic storms and magnetic clouds (MCs) at 1 AU. The velocities are computed using the local correlation tracking (LCT) method. Our computations cover time intervals of 1 10-150 hr, and in four active regions the accumulated helicities due to transient flows are factors of 8-12 larger than the accumulated helicities due to differential rotation. As was first pointed out by DCmoulin and Berger, we suggest that the helicity computed with the LCT method yields not only the helicity injected from shearing motions but also the helicity coming from flux emergence. We compare the computed helicities injected into the corona with the helicities carried away by the CMEs using the MC helicity computations as proxies to the CME helicities. If we assume that the length of the MC flux tubes is I = 2 AU, then the total helicities injected into the corona are a factor of 2.94 lower than the total CME helicities. If we use the values of 1 determined by the condition for the initiation of the kink instability in the coronal flux rope or I = 0.5 AU then the total CME helicities and the total helicities injected into the corona are broadly consistent. Our study, at least partially, clears up some of the discrepancies in the helicity budget of active regions because the discrepancies appearing in our paper are much smaller than the ones reported in previous studies. However, they point out the uncertainties in the MC/CME helicity calculations and also the limitations of the LCT method, which underestimates the computed helicities.

Nindos, A.↗
NASA NTRS2017

Evaluation of the Ozone Fields in NASA's MERRA-2 Reanalysis

The assimilated ozone product from the Modern-Era Retrospective Analysis for Research and Applications, Version 2 (MERRA-2), produced at NASAs Global Modeling and Assimilation Office (GMAO) is summarized. The reanalysis begins in 1980 with the use of retrieved partial-column ozone concentrations from a series of Solar Backscatter Ultraviolet Radiometer (SBUV) instruments on NASA and NOAA spacecraft. Beginning in October 2004, retrieved ozone profiles from the Microwave Limb Sounder (MLS) and total column ozone from the Ozone Monitoring Instrument (OMI) on NASAs EOS Aura satellite are assimilated. While this change in data streams does lead to a discontinuity in the assimilated ozone fields in MERRA-2, making it not useful for studies in decadal (secular) trends in ozone, this choice was made to prioritize demonstrating the value NASAs high-quality research data in the reanalysis context. The MERRA-2 ozone is compared with independent satellite and ozonesonde data, focusing on the representation of the spatial and temporal variability of stratospheric and upper-tropospheric ozone. The comparisons show agreement within 10 (standard deviation of the difference) between MERRA-2 profiles and independent satellite data in most of the stratosphere. The agreement improves after 2004, when EOS Aura data are assimilated. The standard deviation of the differences between the lower-stratospheric and upper-tropospheric MERRA-2 ozone and ozonesondes is 11.2 and 24.5, respectively, with correlations of 0.8 and above. This is indicative of a realistic representation of the UTLS ozone variability in MERRA-2. After 2004, the upper tropospheric ozone in MERRA-2 shows a low bias compared to the sondes, but the covariance with independent observations is improved compared to earlier years. Case studies demonstrate the integrity of MERRA-2 analyses in representing important features such as tropopause folds.

Ozone↗
NASA NTRS1999

Strain Rate Dependent Deformation and Strength Modeling of a Polymer Matrix Composite Utilizing a Micromechanics Approach

Potential gas turbine applications will expose polymer matrix composites to very high strain rate loading conditions, requiring an ability to understand and predict the material behavior under extreme conditions. Specifically, analytical methods designed for these applications must have the capability of properly capturing the strain rate sensitivities and nonlinearities that are present in the material response. The Ramaswamy-Stouffer constitutive equations, originally developed to analyze the viscoplastic deformation of metals, have been modified to simulate the nonlinear deformation response of ductile, crystalline polymers. The constitutive model is characterized and correlated for two representative ductile polymers. Fiberite 977-2 and PEEK, and the computed results correlate well with experimental values. The polymer constitutive equations are implemented in a mechanics of materials based composite micromechanics model to predict the nonlinear, rate dependent deformation response of a composite ply. Uniform stress and uniform strain assumptions are applied to compute the effective stresses of a composite unit cell from the applied strains. The micromechanics equations are successfully verified for two polymer matrix composites. IM7/977-2 and AS4/PEEK. The ultimate strength of a composite ply is predicted with the Hashin failure criteria that were implemented in the composite micromechanics model. The failure stresses of the two composite material systems are accurately predicted for a variety of fiber orientations and strain rates. The composite deformation model is implemented in LS-DYNA, a commercially available transient dynamic explicit finite element code. The matrix constitutive equations are converted into an incremental form, and the model is implemented into LS-DYNA through the use of a user defined material subroutine. The deformation response of a bulk polymer and a polymer matrix composite are predicted by finite element analyses. The results compare reasonably well to experimental values, with some discrepancies. The discrepancies are at least partially caused by the method used to integrate the rate equations in the polymer constitutive model.

Goldberg, Robert K.↗
NASA NTRS2017

Weather Impact on Airport Arrival Meter Fix Throughput

Time-based flow management provides arrival aircraft schedules based on arrival airport conditions, airport capacity, required spacing, and weather conditions. In order to meet a scheduled time at which arrival aircraft can cross an airport arrival meter fix prior to entering the airport terminal airspace, air traffic controllers make regulations on air traffic. Severe weather may create an airport arrival bottleneck if one or more of airport arrival meter fixes are partially or completely blocked by the weather and the arrival demand has not been reduced accordingly. Under these conditions, aircraft are frequently being put in holding patterns until they can be rerouted. A model that predicts the weather impacted meter fix throughput may help air traffic controllers direct arrival flows into the airport more efficiently, minimizing arrival meter fix congestion. This paper presents an analysis of air traffic flows across arrival meter fixes at the Newark Liberty International Airport (EWR). Several scenarios of weather impacted EWR arrival fix flows are described. Furthermore, multiple linear regression and regression tree ensemble learning approaches for translating multiple sector Weather Impacted Traffic Indexes (WITI) to EWR arrival meter fix throughputs are examined. These weather translation models are developed and validated using the EWR arrival flight and weather data for the period of April-September in 2014. This study also compares the performance of the regression tree ensemble with traditional multiple linear regression models for estimating the weather impacted throughputs at each of the EWR arrival meter fixes. For all meter fixes investigated, the results from the regression tree ensemble weather translation models show a stronger correlation between model outputs and observed meter fix throughputs than that produced from multiple linear regression method.

Machine Learning Model.↗
NASA NTRS2002

Participation in the TOMS Science Team

Because of the nominal funding provided by this grant, some of the relevant research is partially funded by other sources. Research performed for this funding period included the following items: We have investigated errors in TOMS ozone measurements caused by the uncertainty in wavelength calibration, coupled with the ozone cross sections in the Huggins bands and their temperature dependence. Preliminary results show that 0.1 nm uncertainty in TOMS wavelength calibration at the ozone active wavelengths corresponds to approx. 1% systematic error in O3, and thus potential 1% biases among ozone trends from the various TOMS instruments. This conclusion will be revised for absolute O3 Measurements as cross sections are further investigated for inclusion in the HITRAN database at the SAO, but the potential for relative errors remains. In order to aid further comparisons among TOMS and GOME ozone measurements, we have implemented our method of direct fitting of GOME radiances (BOAS) for O3, and now obtain the best fitting precision to date for GOME O3 Columns. This will aid in future comparisons of the actual quantities measured and fitted for the two instrument types. We have made comparisons between GOME ICFA cloud fraction and cloud fraction determined from GOME data using the Ring effect in the Ca II lines. There is a strong correlation, as expected, but there are substantial systematic biases between the determinations. This study will be refined in the near future using the recently-developed GOME Cloud Retrieval Algorithm (GOMECAT). We have improved the SAO Ring effect determination to include better convolution with instrument transfer functions and inclusion of interferences by atmospheric absorbers (e.g., O3). This has been made available to the general community.

Chance, Kelly↗
NASA NTRS1990

A study of the topology of dissipating motions in direct numerical simulations of time-developing compressible and incompressible mixing layers

A preliminary investigation of the geometry of flow patterns in numerically simulated compressible and incompressible mixing layers was carried out using 3-D critical point methodology. Motions characterized by high rates of kinetic energy dissipation and/or high enstrophy were of particular interest. In the approach the partial derivatives of the velocity field are determined at every point in the flow. These are used to construct the invariants of the velocity gradient tensor and the rate-of-strain tensor (P, Q, R, and P(sub s), Q(sub s), R(sub s) respectively). For incompressible flow the first invariant is zero. For the conditions of the compressible simulation, the first invariant is found to be everywhere small, relative to the second and third invariants, and so in both cases the local topology at a point is mainly determined by the second and third invariants. The data at every grid point is used to construct scatter plots of Q versus R and Q(sub s) versus R(sub s). Most points map to a cluster near the origin in Q-R space. However, fine scale motions, that is motions which are characterized by velocity derivatives which scale with the square root of R(sub delta), tend to map to regions which lie far from the origin. Definite trends are observed for motions characterized by high enstrophy and/or high dissipation. The observed trends suggest that, for these motions, the second and third invariants of the velocity gradient and rate-of-strain tensors are strongly correlated. The second and third invariants of the rate-of-strain tensor are related by K(-Q(sub s))(exp 3/2), which is consistent with the above scaling of velocity derivatives. The quantity K appears to depend on Reynolds number with an upper limit K = 2(the square root of 3)/9 corresponding to locally axisymmetric flow. For both the compressible and incompressible mixing layer, regions corresponding to high rates of dissipation are found to be characterized by comparable magnitudes of R(sub ij)R(sub ij) and S(sub ij)S(sub ij). For the incompressible mixing layer, regions characterized by the highest values of enstrophy are found to have relatively low strain rates.

Chen, J. H.↗
NASA NTRS1994

Lick optical spectra of quasar HS 1946+7658 at 10 kilometers per second resolution Lyman-alpha forest and metal absorption systems

We present optical spectra of the most luminous known quasi stellar object (QSO) HS 1946+7658 (z(sub em) = 3.051). Our spectra have both full wavelength coverage, 3240-10570 A, and in selected regions, either high signal-to-noise ratio, SNR approximately equals 40-100, or unusually high approximately 10 km/sec resolution, and in parts of the Ly alpha forest and to the red of Ly alpha emission they are among the best published. We find 113 Ly alpha systems and six metal-line systems, three of which are new. The metal systems at z(sub abs) = 2.844 and 3.050 have complex velocity structure with four and three prominent components, respectively. We find that the system at z(sub abs) = 2.844 is a damped Ly alpha absorption (DLA) system, with a neutral hydrogen column density of log N(H I) = 20.2 +/- 0.4, and it is the cause of the Lyman limit break at lambda approximately equals 3520 A. We believe that most of the H I column density in this system is in z(sub abs) = 2.8443 component which shows the strongest low-ionization absorption lines. The metal abundance in the gas phase of the system is (M/H) approximately equals -2.6 +/- 0.3, with a best estimate of (M/H) = -2.8, with ionizaion parameter log gamma = -2.75, from a photoionization model. The ratios of the logarithmic abundances of C, O, Al, and Si are all within a factor of 2 of solar, which is important for two reasons. First, we believe that the gas abundances which we measure are close to the total abundances, because the ratio of aluminum to other elements is near cosmic, and Al is a refractory element which depletes very readily like chromium, in the interstellar medium. Second, we do not see the enhancement of O with respect to C of (O/C) approximately equals 0.5-0.9 reported in three partial Lyman limit systems by Reimers et al. (1992) and Vogel & Reimers (1993); we measure (O/C) = -0.06 for observed ions and (O/C) approximately equals 0.2 after ionization corrections, which is consistent with solar abundances. We see C II*(lambda 1335) offset by 15 km/sec with respect to C II(lambda 1334), presumably because the gas density varies from 2 to 8 cm(exp -3) with changing velocity in the DLA system. These densities imply that the damped component is 6-25 pc thick, which is reasonable for a single cloud in a cold spiral disk. They also imply that the cloud is relatively highly ionized with more C III than C II, more O III than O I, and log N(H I) = 20.72, which is 3 times the H I column. The system at z(abs) = 1.7382 is also believed to be damped with N(H I) approximately equals 10(exp 21) cm(exp -2), because we see Cr II, but its Ly alpha line will never be seen because it is below the Lyman limit of the other DLA system. We see a 2.6 sigma lack of Ly alpha forest lines well away from the QSO redshift, which may be a chance fluctuation. We also see a correlation between column density N(H I) and Doppler parameter b for 96 unsaturated Ly Alpha forest absorption lines, and although this correlation persists in the 36 Ly alpha lines which lie in regions where the SNR approximately equals 8-16, we agree with Rauch et al. (1993) that it is probably a bogus effect of low supernova remnant (SNR). The same applies to lines with very low b values: in regions where SNR less than or equal to 8 we see many Ly alpha lines which appear to have 10 less than or equal to b less than or equal to 20, but when 8 less than or equal to SNR less than or equal to 16 we see only one line with b less than or equal to 15 km/sec, and two others which we believe have b less than or equal to 20, with values of 20 and 16 km/sec. Traditional Ly alpha line samples which include all lines which have W/sigma(W) greater than or equal to 4 are not adequate to explore the distribution of the properties of individual clouds, because we need much higher (W/sigma(W)) and SNR to avoid the strong biases.

Fan, Xiao-Ming↗
NASA NTRS2010

An Experimental Study of Small-Scale Variability of Raindrop Size Distribution

An experimental study of small-scale variability of raindrop size distributions (DSDs) has been carried out at Wallops Island, Virginia. Three Joss-Waldvogel disdrometers were operated at a distance of 0.65, 1.05, and 1.70 km in a nearly straight line. The main purpose of the study was to examine the variability of DSDs and its integral parameters of liquid water content, rainfall, and reflectivity within a 2-km array: a typical size of Cartesian radar pixel. The composite DSD of rain events showed very good agreement among the disdrometers except where there were noticeable differences in midsize and large drops in a few events. For consideration of partial beam filling where the radar pixel was not completely covered by rain, a single disdrometer reported just over 10% more rainy minutes than the rainy minutes when all three disdrometers reported rainfall. Similarly two out of three disdrometers reported5%more rainy minutes than when all three were reporting rainfall. These percentages were based on a 1-min average, and were less for longer averaging periods. Considering only the minutes when all three disdrometers were reporting rainfall, just over one quarter of the observations showed an increase in the difference in rainfall with distance. This finding was based on a 15-min average and was even less for shorter averaging periods. The probability and cumulative distributions of a gamma-fitted DSD and integral rain parameters between the three disdrometers had a very good agreement and no major variability. This was mainly due to the high percentage of light stratiform rain and to the number of storms that traveled along the track of the disdrometers. At a fixed time step, however, both DSDs and integral rain parameters showed substantial variability. The standard deviation (SD) of rain rate was near 3 mm/h, while the SD of reflectivity exceeded 3 dBZ at the longest separation distance. These standard deviations were at 6-min average and were higher at shorter averaging periods. The correlations decreased with increasing separation distance. For rain rate, the correlations were higher than previous gauge-based studies. This was attributed to the differences in data processing and the difference in rainfall characteristics in different climate regions. It was also considered that the gauge sampling errors could be a factor. In this regard, gauge measurements were simulated employing existing disdrometer dataset. While a difference was noticed in cumulative distribution of rain occurrence between the simulated gauge and disdrometer observations, the correlations in simulated gauge measurements did not differ from the disdrometer measurements.

Tokay, Ali↗