Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variability 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

Investigation of MMOD Impact on STS-115 Shuttle Payload Bay Door Radiator

The Orbiter radiator system consists of eight individual 4.6 m x 3.2 m panels located with four on each payload bay door. Forward panels #1 and #2 are 2.3 cm thick while the aft panels #3 and #4 have a smaller overall thickness of 1.3 cm. The honeycomb radiator panels consist of 0.028 cm thick Aluminum 2024-T81 facesheets and Al5056-H39 cores. The face-sheets are topped with 0.005 in. (0.127 mm) silver-Teflon tape. The radiators are located on the inside of the shuttle payload bay doors, which are closed during ascent and reentry, limiting damage to the on-orbit portion of the mission. Post-flight inspections at the Kennedy Space Center (KSC) following the STS-115 mission revealed a large micrometeoroid/orbital debris (MMOD) impact near the hinge line on the #4 starboard payload bay door radiator panel. The features of this impact make it the largest ever recorded on an orbiter payload bay door radiator. The general location of the damage site and the adjacent radiator panels can be seen in Figure 2. Initial measurements of the defect indicated that the hole in the facesheet was 0.108 in. (2.74 mm) in diameter. Figure 3 shows an image of the front side damage. Subsequent observations revealed exit damage on the rear facesheet. Impact damage features on the rear facesheet included a 0.03 in. diameter hole (0.76 mm), a approx.0.05 in. tall bulge (approx.1.3 mm), and a larger approx.0.2 in. tall bulge (approx.5.1 mm) that exhibited a crack over 0.27 in. (6.8 mm) long. A large approx.1 in. (25 mm) diameter region of the honeycomb core was also damaged. Refer to Figure 4 for an image of the backside damage to the panel. No damage was found on thermal blankets or payload bay door structure under the radiator panel. Figure 5 shows the front facesheet with the thermal tape removed. Ultrasound examination indicated a maximum facesheet debond extent of approximately 1 in. (25 mm) from the entry hole. X-ray examinations revealed damage to an estimated 31 honeycomb cells with an extent of 0.85 in. x 1.1 in. (21.6 x 27.9 mm). Pieces of the radiator at and surrounding the impact site were recovered during the repair procedures at KSC. They included the thermal tape, front facesheet, honeycomb core, and rear facesheet. These articles were examined at JSC using a scanning electron microscope (SEM) with an energy dispersive x-ray spectrometer (EDS). Figure 6 shows SEM images of the entry hole in the facesheet. The asymmetric height of the lip may be attributed to projectile shape and impact angle. Numerous instances of a glass-fiber organic matrix composite were observed in the facesheet tape sample. The fibers were approximately 10 micrometers in diameter and variable lengths. EDS analysis indicated a composition of Mg, Ca, Al, Si, and O. Figures 7 and 8 present images of the fiber bundles, which were believed to be circuit board material based on similarity in fiber diameter, orientation, consistency, and composition. A test program was initiated in an attempt to simulate the observed damage to the radiator facesheet and honeycomb. Twelve test shots were performed using projectiles cut from a 1.6 mm thick fiberglass circuit board substrate panel. Results from test HITF07017, shown in figures 9 and 10, correlates with the observed impact features reasonably well. The test was performed at 4.14 km/sec with an impact angle of 45 degrees using a cylindrical projectile with a diameter and length of 1.25 mm. The fiberglass circuit board material had a density of 1.65 g/cu cm, giving a projectile mass of 2.53 mg. An analysis was performed using the Bumper code to estimate the probability of impact to the shuttle from a 1.25 mm diameter particle. Table 1 shows a 1.6% chance (impact odds = 1 in 62) of a 1.25 mm or larger MMOD impact on the radiators of the vehicle during a typical ISS mission. There is a 0.4% chance (impact odds = 1 in 260) that a 1.25 mm or larger MMOD particle would impact the RCC wing leading edge and nose cap during a typical miion. Figure 11 illustrates the vulnerable areas of the wing leading edge reinforced carbon-carbon (RCC), an area of the vehicle that is very sensitive to impact damage. The highlighted red, orange, yellow, and light green areas would be expected to experience critical damage if impacted by an OD particle such as the one that hit the RH4 radiator panel on STS-115.

Hyde, J.↗

Gravity and Rossby Wave Signatures in the Tropical Troposphere and Lower Stratosphere Based on Southern Hemisphere Additional Ozonesondes (SHADOZ), 1998-2007

Prior investigations attempted to determine the relative influence of advection and convective processes on ozone and water vapor distributions in the tropical tropopause layer (TTL) through analyses of tracers, related physical parameters (e.g., outgoing long-wave radiation, precipitable water, and temperature), or with models. In this study, stable laminae in Southern Hemisphere Additional Ozonesonde Network (SHADOZ) ozone profIles from 1998 to 2007 are interpreted in terms of gravity waves (GW) or Rossby waves (RW) that are identified with vertical and quasi-horizontal displacements, respectively. Using the method of Pierce and Grant (1998) as applied by Thompson et al. (2007a, 2007b, 2010, 2011), amplitudes and frequencies in ozone laminae are compared among representative SHADOZ sites over Africa and the Pacific, Indian, and Atlantic oceans. GW signals maximize in the TTL and lower stratosphere. Depending on site and season, GW are identified in up to 90% of the soundings. GW are most prevalent over the Pacific and eastern Indian oceans, a distribution consistent with vertically propagating equatorial Kelvin waves. Ozone laminae from RW occur more often below the tropical tropopause and with lower frequency 20%). Gravity wave and Rossby wave indices (GWI, RWI) are formulated to facilitate analysis of interannual variability of wave signatures among sites. GWI is positively correlated with a standard ENSO (El Nino-Southern Oscillation) index over American Samoa (14degS, 171degW) and negatively correlated at Watukosek, Java (7.5degS, 114degE), Kuala Lumpur (3degN, 102degE), and Ascension Island (80degS, 15degW). Generally, the responses of GW and RW to ENSO are consistent with prior studies.

Thompson, Anne M.↗

Sea Ice Thickness, Freeboard, and Snow Depth products from Operation IceBridge Airborne Data

The study of sea ice using airborne remote sensing platforms provides unique capabilities to measure a wide variety of sea ice properties. These measurements are useful for a variety of topics including model evaluation and improvement, assessment of satellite retrievals, and incorporation into climate data records for analysis of interannual variability and long-term trends in sea ice properties. In this paper we describe methods for the retrieval of sea ice thickness, freeboard, and snow depth using data from a multisensor suite of instruments on NASA's Operation IceBridge airborne campaign. We assess the consistency of the results through comparison with independent data sets that demonstrate that the IceBridge products are capable of providing a reliable record of snow depth and sea ice thickness. We explore the impact of inter-campaign instrument changes and associated algorithm adaptations as well as the applicability of the adapted algorithms to the ongoing IceBridge mission. The uncertainties associated with the retrieval methods are determined and placed in the context of their impact on the retrieved sea ice thickness. Lastly, we present results for the 2009 and 2010 IceBridge campaigns, which are currently available in product form via the National Snow and Ice Data Center

sea ice thickness↗

Formaldehyde Column Density Measurements as a Suitable Pathway to Estimate Near-Surface Ozone Tendencies from Space

In support of future satellite missions that aim to address the current shortcomings in measuring air quality from space, NASA's Deriving Information on Surface Conditions from Column and Vertically Resolved Observations Relevant to Air Quality (DISCOVER-AQ) field campaign was designed to enable exploration of relationships between column measurements of trace species relevant to air quality at high spatial and temporal resolution. In the DISCOVER-AQ data set, a modest correlation (r2 = 0.45) between ozone (O3) and formaldehyde (CH2O) column densities was observed. Further analysis revealed regional variability in the O3-CH2O relationship, with Maryland having a strong relationship when data were viewed temporally and Houston having a strong relationship when data were viewed spatially. These differences in regional behavior are attributed to differences in volatile organic compound (VOC) emissions. In Maryland, biogenic VOCs were responsible for approx.28% of CH2O formation within the boundary layer column, causing CH2O to, in general, increase monotonically throughout the day. In Houston, persistent anthropogenic emissions dominated the local hydrocarbon environment, and no discernable diurnal trend in CH2O was observed. Box model simulations suggested that ambient CH2O mixing ratios have a weak diurnal trend (+/-20% throughout the day) due to photochemical effects, and that larger diurnal trends are associated with changes in hydrocarbon precursors. Finally, mathematical relationships were developed from first principles and were able to replicate the different behaviors seen in Maryland and Houston. While studies would be necessary to validate these results and determine the regional applicability of the O3-CH2O relationship, the results presented here provide compelling insight into the ability of future satellite missions to aid in monitoring near-surface air quality.

Surface Conditions from Column and Vertically Reso↗

Improving NASA GEOS Atmospheric CO2 Simulations by Calibrating CASA Surface Fluxes with an Empirical Sink

With the adoption of the Paris climate accord, efforts to monitor and understand both anthropogenic and natural carbon sources and sinks are increasing across the world. Given their low latency and global coverage, satellite observations of atmospheric carbon dioxide (CO2) are poised to make important contributions to this field. The combination of satellite data and high resolution global models can be used to monitor changes in carbon fluxes and to evaluate the consistency of nationally reported emissions estimates in support of multiple stakeholder communities. However, a consistent challenge to such work has been the high latency of surface carbon flux estimates, which are often not available for a year or more. This presentation describes the construction of surface carbon flux estimates meant to improve the near real time simulation of atmospheric CO2 with NASA's Goddard Earth Observing System (GEOS) general circulation model. The surface flux estimates begin with a collection of bottom-up fluxes which incorporate satellite measurements in their construction, e.g. vegetation indices in the Carnegie-Ames-Stanford Approach (CASA) and nighttime lights in the Open-source Data Inventory for Anthropogenic CO2 (ODIAC). From there, we take the additional step of using an empirical sink to calibrate terrestrial net biospheric exchange (NBE) to estimated values from atmospheric inversion systems. This approach removes a known, systematic bias in predicted atmospheric mixing ratios. Using these fluxes in a free running simulation, the model is able to reproduce in situ measurements with the same skill as when it uses gridded fluxes from a flux inversion system. Using these fluxes as a prior in an assimilation system, e.g. one incorporating retrievals of column CO2 from the Orbiting Carbon Observatory 2 (OCO-2), allows the analysis to capture variability in CO2 on scales that would be missed otherwise. This approach supports NASA's capability to forecast atmospheric CO2 up to two weeks in advance by leveraging a GEOS system used to produce quasi-operational weather analyses and forecasts, providing a valuable new tool to the carbon monitoring research and applications communities.

Weir, B.↗

Relating Land Cover Change to Runoff Distribution Using NASA Earth Observations in Riley County, Kansas

Riley County, Kansas, has observed increased levels of flooding, potentially due to changes in land use/land cover (LULC) and seasonal vegetation variation. This study contrasts two methods of generating runoff curve numbers (CN) from 2006-2020. (1) The traditional Soil Conservation Service CN calculation method uses a look-up table and tracked LULC to determine runoff changes. These tables allow for CN values specific to various land and crop cover types and account for various farming best practices but lack flexibility in calculations for various seasons or plant health. (2) A dynamic method employs normalized difference vegetation index (NDVI) compiled over the rainy season each year to calculate CN using seasonal vegetation. This method allows for a more precise analysis of runoff variability within and between rainy seasons because it can be updated with greater temporal detail and it captures higher spatial resolutions by using NDVI as a proxy for LULC. This study further uses inputs from the United States Geologic Survey (USGS) National Land Cover Database (NLCD), the United States Department of Agriculture (USDA) Cropland Data Layer, and Landsat imagery to create more precise LULC raster datasets including both urban cover and crop-specific land use and curve number maps of the area. Results can guide decision makers in the City of Manhattan, Riley County Department of Planning and Development, Riley County Conservation District, the Kansas Forest Service, and the Kansas Department of Health and Environment toward informed decisions on resiliency strategies to address future flooding.

Ella Griffith↗

Riley County Water Resources Project Summary - Comparing Runoff Curve Calculation Methods to Inform Local Resiliency Initiatives in Riley County, Kansas

Riley County, Kansas, has observed increased levels of flooding, potentially due to changes in land use/land cover (LULC) and seasonal vegetation variation. This study contrasts two methods of generating runoff curve numbers (CN) from 2006-2020. (1) The traditional Soil Conservation Service CN calculation method uses a look-up table and tracked LULC to determine runoff changes. These tables allow for land cover-specific CN and account for various farming techniques but lack flexibility in calculations for various seasons or plant health. (2) A dynamic method employs normalized difference vegetation index (NDVI) compiled over the rainy season each year to calculate CN using seasonal vegetation. This method allows for a more precise analysis of runoff variability within and between rainy seasons because it can be updated with greater temporal detail and captures higher spatial resolutions by using NDVI as a proxy for LULC. This study further uses inputs from the United States Geological Survey (USGS) National Land Cover Database (NLCD), the United States Department of Agriculture (USDA) Cropland Data Layer, and Landsat imagery to create more precise LULC raster datasets including both urban cover and crop-specific land use and curve number maps of the area. Results can guide decision makers in the City of Manhattan, Riley County Department of Planning and Development, Riley County Conservation District, the Kansas Forest Service, and the Kansas Department of Health and Environment toward informed decisions on resiliency strategies to address future flooding.

DEVELOP Project Summary↗

Relating Land Cover Change to Runoff Distribution Using NASA Earth Observations in Riley County, Kansas

In recent years, Riley County, Kansas experienced high levels of flooding potentially due to changes in land use and impervious surface. To analyze the increase in floods and runoff patterns over time, this study contrasted two methods of curve number (CN) calculation, which estimated runoff from rainfall events, from 2006-2020. The conventional calculation method used a look-up table and tracked land use/land cover change. This method resulted in static tables from various inputs such as soil type, surface imperviousness, slope, and crop phenology to estimate runoff versus infiltration. These tables allowed for phenology specific curve numbers and accounted for various farming techniques but required many inputs and lacked flexibility to include seasonal variations in vegetation. In contrast, the dynamic method is a more modern and simple calculation employing normalized difference vegetation index (NDVI) data from Landsat 5, Landsat 7, and Landsat 8 satellites compiled over the rainy season each year to calculate curve numbers based off seasonal vegetation. This dynamic method allowed for more precise analysis of runoff variability because NDVI is a more accurate, real-time measurement for land cover and it could be calculated with greater temporal detail. The results will allow decision makers to make more informed decisions on resiliency strategies to address future flooding.

Trista Brophy↗

Larch (Larix dahurica Turcz) growth response to climate change in the Siberian permafrost zone

Larch-dominant communities are the most extensive high-latitude forests in Eurasia and are experiencing the strongest impacts from warming temperatures. We analyzed larch (Larix dahurica Turcz) growth index (GI) response to climate change. The studied larch-dominant communities are located within the permafrost zone of Northern Siberia at the northern tree limit (ca. N 67° 38′, E 99° 07′). Methods included dendrochronology, analysis of climate variables, root zone moisture content, and satellitederived gross (GPP) and net (NPP) primary productivity. It was found that larch response to warming included a period of increased annual growth increment (GI) (from the 1970s to ca. 1995) with a follow on GI decline. Increase in GI correlated with summer air temperature, whereas an observed decrease in GI was caused by water stress (vapor pressure deficit and drought increase). Water stress impact on larch growth in permafrost was not observed before the onset of warming (ca. 1970). Water limitation was also indicated by GI dependence on soil moisture stored during the previous year. Water stress was especially pronounced for stands growing on rocky soils with low water-holding capacity. GPP of larch communities showed an increasing trend, whereas NPP stagnated. A similar pattern of GI response to climate warming has also been observed for Larix sibirica Ledeb, Pinus sibirica Du Tour, and Abies sibirica Ledeb in the forests of southern Siberia. Thus, warming in northern Siberia permafrost zone resulted in an initial increase in larch growth from the 1970s to the mid-1990s. After that time, larch growth increment has decreased. Since ca. 1990, water stress at the beginning of the vegetative period became, along with air temperature, a main factor affecting larch growth within the permafrost zone.

Larix dahurica↗

Using Boosted Decision Trees to Select High Quality Measurements in the Mu2e Experiment at Fermilab

This thesis presents the implementation and evaluation of a Boosted Decision Tree (BDT) model to improve the selection of high-quality track measurements in the Mu2e experiment at Fermilab. The Mu2e experiment is a high-energy physics experiments seeking to observe a rare theoretical physics process known as Charged Lepton Flavor Violation. A significant challenge faced by the Mu2e experiment are so-called background events, which are events whose data mimics that of the rare physics process the experiment seeks to observe. Without a mechanism to reduce background, it would be impossible to know whether Charged Lepton Flavor Violation occurred or not. To this end, high-quality track measurements must be distinguished from low-quality track measurements. A track can be conceived of as the reconstructed path of a particle that traveled through the Mu2e detector. In addition to other data, data about such tracks is stored using a C++-based framework, specific to the domain of high-energy physics, known as ROOT. A boosted decision tree model was trained using ROOT’s Toolkit For Multivariate Analysis by leveraging variables ancillary to track quality. In evaluation, the BDT achieves a ROC-AUC of 0.927 in discriminating good-quality tracks from poor-quality tracks. Such a score is indicative of both strong discrimination and strong generalization. Subsequently, it is shown that applying a BDT-based quality cut to the distribution of particle momenta significantly enhances the signal-to-background distinction for signal electrons, paving the way for improved sensitivity to Charged Lepton Flavor Violation.

Mullany, Brendan T. [Drew U.] (ORCID:0009000818888↗

Velocity structure in long period variable star atmospheres

A regression analysis of the dependence of absorption line velocities on wavelength, line strength, excitation potential, and ionization potential is presented. The method determines the region of formation of the absorption lines for a given data and wavelength region. It is concluded that the scatter which is frequently found in velocity measurements of absorption lines in long period variables is probably the result of a shock of moderate amplitude located in or near the reversing layer and that the frequently observed correlation of velocity with excitation and ionization are a result of the velocity gradients produced by this shock in the atmosphere. A simple interpretation of the signs of the coefficients of the regression analysis is presented in terms of preshock, post shock, or across the shock, together with criteria for evaluating the validity of the fit. The amplitude of the reversing layer shock is estimated from an analysis of a series of plates for four long period variable stars along with the most probable stellar velocity for these stars.

Pilachowski, C.↗

Creation of an Upper Stage Trajectory Capability Boundary to Enable Booster System Trade Space Exploration

The problem of trajectory optimization is important in all space missions. The solution of this problem enables one to specify the optimum thrust steering program which should be followed to achieve a specified mission objective, simultaneously satisfying the constraints.1 It is well known that whether or not the ascent trajectory is optimal can have a significant impact on propellant usage for a given payload, or on payload weight for the same gross vehicle weight.2 Consequently, ascent guidance commands are usually optimized in some fashion. Multi-stage vehicles add complexity to this analysis process as changes in vehicle properties in one stage propagate to the other stages through gear ratios and changes in the optimal trajectory. These effects can cause an increase in analysis time as more variables are added and convergence of the optimizer to system closure requires more analysis iterations. In this paper, an approach to simplifying this multi-stage problem through the creation of an upper stage capability boundary is presented. This work was completed as part of a larger study focused on trade space exploration for the advanced booster system that will eventually form a part of NASA s new Space Launch System.3 The approach developed leverages Design of Experiments and Surrogate Modeling4 techniques to create a predictive model of the SLS upper stage performance. The design of the SLS core stages is considered fixed for the purposes of this study, which results in trajectory parameters such as staging conditions being the only variables relevant to the upper stage. Through the creation of a surrogate model, which takes staging conditions as inputs and predicts the payload mass delivered by the SLS upper stage to a reference orbit as the response, it is possible to identify a "surface" of staging conditions which all satisfy the SLS requirement of placing 130 metric tons into low-Earth orbit (LEO).3 This identified surface represents the 130 metric ton capability boundary for the upper stage, such that if the combined first stage and boosters can achieve any one staging point on that surface, then the design is identified as feasible. With the surrogate model created, design and analysis of advanced booster concepts is streamlined, as optimization of the upper stage trajectory is no longer required in every design loop.

Walsh, Ptrick↗

Sensitivity analysis of a deployable three longeron truss beam designed for minimum member loads during deployment

The significant design variables of a deployable three longeron truss beam, designed for minimum member loads during deployment, are identified and the sensitivity of the load level in the members of the truss beam due to variations of the design variables is established. The analysis of the deployment models, developed in this paper, is carried out using a commercially available computer code, called DADS. Based on the results of the analysis, guidelines are formulated for the design of a deployable three longeron truss beam to achieve minimum loads in the members during deployment. The paper concludes with a brief discussion of the applicability of the approach, taken in this paper, to other truss configurations.

Warnaar, Dirk B.↗

An improved multiple linear regression and data analysis computer program package

NEWRAP, an improved version of a previous multiple linear regression program called RAPIER, CREDUC, and CRSPLT, allows for a complete regression analysis including cross plots of the independent and dependent variables, correlation coefficients, regression coefficients, analysis of variance tables, t-statistics and their probability levels, rejection of independent variables, plots of residuals against the independent and dependent variables, and a canonical reduction of quadratic response functions useful in optimum seeking experimentation. A major improvement over RAPIER is that all regression calculations are done in double precision arithmetic.

Sidik, S. M.↗

Analysis and testing of two-dimensional slot nozzle ejectors with variable area mixing sections

Finite difference computer techniques have been used to calculate the detailed performance of air to air two-dimensional ejectors with symmetric variable area mixing sections and coaxial coverging primary nozzles. The analysis of the primary nozzle assumed correct expansion of the flow and is suitable for subsonic and slightly supersonic velocity levels. The variation of the mixing section channel walls is assumed to be gradual so that the static pressure can be assumed uniform on planes perpendicular to the axis. A test program was run to provide two-dimensional ejector test data for verification of the computer analysis. A primary converging nozzle with a discharge geometry of 0.125 inch x 8.0 inch was supplied with 600 SCFM of air at about 35 psia and 180 F. This nozzle was combined with two mixing section geometries with throat sizes of 1.25 inch x 8.0 inch and 1.875 inch x 8.0 inch and was tested at a total of 11 operating points.

Gilbert, G. B.↗

Analysis and testing of two-dimensional slot nozzle ejectors with variable area mixing sections

Finite difference computer techniques have been used to calculate the detailed performance of air-to-air two-dimensional ejectors with symmetric variable area mixing sections and coaxial converging primary nozzles. The analysis of the primary nozzle assumed correct expansion of the flow and is suitable for subsonic and slightly supersonic velocity levels. The variation of the mixing section channel walls is assumed to be gradual so that the static pressure can be assumed uniform on planes perpendicular to the axis. A test program was run to provide two-dimensional ejector test data for verification of the computer analysis. A primary converging nozzle with a discharge geometry of 3.18 mm x 203 mm was supplied with 0.340 kg/sec of air at about 2.43 bar and 356 K. This nozzle was combined with two mixing section geometries with throat sizes of 31.8 mm x 203 mm and 47.6 mm x 203 mm and was tested at a total of 11 operating points.

Gilbert, G. B.↗

Aeroelastic analysis for helicopter rotor blades with time-variable, non-linear structural twist and multiple structural redundancy: Mathematical derivation and program user's manual

The differential equations of motion for the lateral and torsional deformations of a nonlinearly twisted rotor blade in steady flight conditions together with those additional aeroelastic features germane to composite bearingless rotors are derived. The differential equations are formulated in terms of uncoupled (zero pitch and twist) vibratory modes with exact coupling effects due to finite, time variable blade pitch and, to second order, twist. Also presented are derivations of the fully coupled inertia and aerodynamic load distributions, automatic pitch change coupling effects, structural redundancy characteristics of the composite bearingless rotor flexbeam - torque tube system in bending and torsion, and a description of the linearized equations appropriate for eigensolution analyses. Three appendixes are included presenting material appropriate to the digital computer program implementation of the analysis, program G400.

Bielawa, R. L.↗