Search NASA⌕ Search

SEARCH · Search NASA

Results for “higher order statistics”

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.

191 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

Ambulation Increases Decompression Sickness in Spacewalk Simulations

Musculoskeletal activity has the potential to both improve and compromise decompression safety. Exercise enhances inert gas elimination during oxygen breathing prior to decompression (prebreathe), but it may also promote bubble nuclei formation (nucleation), which can lead to gas phase separation and bubble growth and increase the risk of decompression sickness (DCS). The timing, pattern and intensity of musculoskeletal activity and the level of tissue supersaturation may be critical to the net effect. Understanding the relationships is important to evaluate exercise prebreathe protocols and quantify decompression risk in gravity and microgravity environments. Data gathered during NASA's Prebreathe Reduction Program (PRP) studies combined oxygen prebreathe and exercise followed by low pressure (4.3 psi; altitude equivalent of 30,300 ft [9,235 m]) microgravity simulation to produce two protocols used by astronauts preparing for extravehicular activity. Both the Phase II/CEVIS (cycle ergometer vibration isolation system) and ISLE (in-suit light exercise) trials eliminated ambulation to more closely simulate the microgravity environment. The CEVIS results (35 male, 10 female) serve as control data for this NASA/Duke study to investigate the influence of ambulation exercise on bubble formation and the subsequent risk of DCS. METHODS Four experiments will replicate the CEVIS exercise-enhanced oxygen prebreathe protocol, each with a different exception. The first of these is currently underway. Experiment 1 - Subjects complete controlled ambulation (walking in place with fixed cadence and step height) during both preflight and at 4.3 psi instead of remaining nonambulatory throughout. Experiment 2 - Subjects remain non-ambulatory during the preflight period and ambulatory at 4.3 psi. Experiment 3 - Subjects ambulate during the preflight period and remain non-ambulatory at 4.3 psi. Experiment 4 - The order of heavy and light exercise employed in the CEVIS protocol is reversed, with the light exercise occurring first (subjects remain non-ambulatory throughout). Decompression stress is assessed with non-invasive ultrasound during each of 14 epochs of a 4 hour simulated spacewalk at 4.3 psi; aural Doppler is used to monitor bubbles (Spencer grade 0-IV scale) passing through the pulmonary artery, and two-dimensional echocardiographic imaging is used to look for left ventricular gas emboli (LVGE; the presence of which is a test termination criterion). Venous blood is collected at baseline and twice following repressurization to determine if the decompression stress is correlated with microparticles (cell fragments) accumulation. The plan is to test 25-50 subjects in each experiment. Fisher Exact Tests (one-tailed) are used to compare test and control groups. Trials are suspended when the DCS or grade IV VGE observations reach 70% confidence of DCS risk >15% and grade IV VGE risk >20%. RESULTS Experiment 1 was concluded with 20 complete trials (15 male, 5 female) since the statistical outcome would not change with five additional trials. The observed DCS was significantly greater in Experiment 1 than in CEVIS trials (4/20 [20%] vs. 0/45 [0%], respectively, p=0.007), as was the frequency of peak grade IV VGE (6/21 [29%; including one additional subject that presented grade IV VGE but whose trial was ended before completion when LVGE were observed] vs. 3/45 [7%], respectively, p=0.024). Experiment 3 trials are now underway, with 11 trials completed (10 male, 1 female). Preliminary results indicate no difference in observed DCS between Experiment 3 and CEVIS trials (1/11 [9%] vs. 0/45 [0%], respectively, p=0.196), or between Experiment 3 and Experiment 1 trials (p=0.405). The frequency of peak grade IV VGE in Experiment 3 (2/11 [18%]) did not differ from CEVIS or Experiment 1 trials (p=0.251 and p=0.425, respectively). Microparticle patterns are widely variable and still under analysis. DISCUSSION The results of the Experiment 1 trials support the thesis that decompression stress is increased by ambulation exercise, given the higher incidence of DCS and grade IV VGE when compared to the non-ambulatory PRP CEVIS trials. Experiment 3 trials are incomplete, but suggest that the effect of ambulation during ground level preflight oxygen breathing alone, when subjects are undersaturated with inert gas, may not differ in risk from ambulation at both preflight and spacesuit pressures, the latter when subjects are supersaturated with inert gas. Further trials are needed to confirm the relative effects of ambulation in undersaturated vs. supersaturated states and to determine whether light exercise facilitates the removal of heavy exercise-induced nucleation (Experiment 4).

Pollock, N. W.↗
NASA NTRS2010

Selection Effects in Identifying Magnetic Clouds and the Importance of the Closest Approach Parameter

This study is motivated by the unusually low number of magnetic clouds (MCs) that are strictly identified within interplanetary coronal mass ejections (ICMEs), as observed at 1 AU; this is usually estimated to be around 30% or lower. But a looser definition of MCs may significantly increase this percentage. Another motivation is the unexpected shape of the occurrence distribution of the observers' "closest approach distances" (measured from a MC's axis, and called CA) which drops off somewhat rapidly as |CA| (in % of MC radius) approaches 100%, based on earlier studies. We suggest, for various geometrical and physical reasons, that the |CA|-distribution should be somewhere between a uniform one and the one actually observed, and therefore the 30% estimate should be higher. So we ask, When there is a failure to identify a MC within an ICME, is it occasionally due to a large |CA| passage, making MC identification more difficult, i.e., is it due to an event selection effect? In attempting to answer this question we examine WIND data to obtain an accurate distribution of the number of MCs vs. |CA| distance, whether the event is ICME-related or not, where initially a large number of cases (N=98) are considered. This gives a frequence distribution that is far from uniform, confirming earlier studies. This along with the fact that there are many ICME identification-parameters that do not depend on |CA| suggest that, indeed an MC event selection effect may explain at least part of the low ratio of (No. MCs)/(No. ICMEs). We also show that there is an acceptable geometrical and physical consistency in the relationships for both average "normalized" magnetic field intensity change and field direction change vs. |CA| within a MC, suggesting that our estimates of |CA|, B(sub 0) (magnetic field intensity on the axis), and choice of a proper "cloud coordinate" system (all needed in the analysis) are acceptably accurate. Therefore the MC fitting model (Lepping et al., 1990) is adequate, on average, for our analysis. However, this selection effect is not likely to completely answer our original question, on the unexpected ratio of MCs to ICMEs, so we must look for other factors, such as peculiarities of CME birth conditions. As a by-product of this analysis, we determine that the first order structural effects within a MC due to its interaction with the solar wind, plus the MC's usual expansion at 1 AU (i.e., the non-force free components of the MC's field) are, on average, weakly dependent on radial distance from the MC's axis; that is, in the outer reaches of a typical MC the non-force free effects show up, but even there they are rather weak. Finally, we show that it is not likely that a MC's size distribution statistically controls the occurrence distribution of the estimated |CA|s.

Lepping, R. P.↗
NASA NTRS1998

Scaling Properties of Particle Density Fields Formed in Simulated Turbulent Flows

Direct numerical simulations (DNS) of particle concentrations in fully developed 3D turbulence were carried out in order to study the nonuniform structure of the particle density field. Three steady-state turbulent fluid fields with Taylor microscale Reynolds numbers (Re(sub lambda)) of 40, 80 and 140 were generated by solving the Navier-Stokes equations with pseudospectral methods. Large scale forcing was used to drive the turbulence and maintain temporal stationarity. The response of the particles to the fluid was parameterized by the particle Stokes number St, defined as the ratio of the particle's stopping time to the mean period of eddies on the Kolmogorov scale (eta). In this paper, we consider only passive particles optimally coupled to these eddies (St approx. = 1) because of their tendency to concentrate more than particles with lesser or greater St values. The trajectories of up to 70 million particles were tracked in the equilibrated turbulent flows until the particle concentration field reached a statistically stationary state. The nonuniform structure of the concentration fields was characterized by the multifractal singularity spectrum, f(alpha), derived from measures obtained after binning particles into cells ranging from 2(eta) to 15(eta) in size. We observed strong systematic variations of f(alpha) across this scale range in all three simulations and conclude that the particle concentration field is not statistically self similar across the scale range explored. However, spectra obtained at the 2(eta), 4(eta), and 8(eta) scales of each flow case were found to be qualitatively similar. This result suggests that the local structure of the particle concentration field may be flow-Independent. The singularity spectra found for 2n-sized cells were used to predict concentration distributions in good agreement with those obtained directly from the particle data. This Singularity spectrum has a shape similar to the analogous spectrum derived for the inertial-range energy dissipation fields of experimental turbulent flows at Re(sub lambda) = 110 and 1100. Based on this agreement, and the expectation that both dissipation and particle concentration are controlled by the same cascade process, we hypothesize that singularity spectra similar to the ones found in this work provide a good characterization of the spatially averaged statistical properties of preferentially concentrated particles in higher Re(sub lambda) turbulent flows.

Hogan, Robert C.↗
NASA NTRSDate not supplied

TPSAS-NF1676L-17940-DND

The NASA/GEWEX SRB (Global Energy and Water Exchanges, Surface Radiation Budget) project produces and archives shortwave and longwave radiation budget flux estimates at the top of the atmosphere and at the Earth's surface. The latest version in the archive, Release 3.0, is available as 3-hourly, 3-hourly-monthly, daily and monthly means continuously over the period from July 1983 to December 2007 on a quasi-equal-area grid system of 44016 grid boxes. SRB Release 4 with further improvements in data quality and higher spatial resolution is being developed. The SRB shortwave/longwave fluxes at the Earth's surface from the algorithm GSW(V3.0)/GLW(V3.1) have been extensively validated against high-quality ground-based observations, in particular observed data from the Baseline Surface Radiation Network (BSRN). Comparisons with nearly 6000 site-months of both shortwave and longwave data from 52 BSRN sites show generally good agreement. In addition, the GEWEX SRB data have also been found to compare favorably with the World Radiation Data Centre (WRDC) data and the Global Energy Balance Archive (GEBA) data. In spite of the fact that the BSRN sites are scattered on all seven continents, the validation of the SRB data over the vast oceans had not been done until recently. In this paper, we present comparisons of the GEWEX-SRB data shortwave/longwave data with observations made on arrays moored in tropical oceans. Specifically, we have data from 21 buoys, or moorings, from Predictions and Research Moored Array in the Atlantic (PIRATA), 14 buoys from the Research Moored Array for African-Asian-Australian Monsoon Analysis and Prediction (RAMA) in the Indian Ocean, 20 buoys from the Tropical Atmosphere Ocean (TAO) array in the Pacific, and 3 buoys from the Woods Hole Oceanographic Institute (WHOI) (2 in the Pacific and 1 in the Atlantic). The data from these buoys span 12 years from 2000 to 2011, though not necessarily continuously. It is found that except for occasional outliers, the majority of the comparable pairs of data points agree well. Further examination of the outliers is necessary to determine whether they can be attributed to wind-driven tilt and/or aerosol accumulation on the instruments. In order to see how the ocean-based validation compares with land-based validation, we made SRB-BSRN monthly mean shortwave/longwave radiation comparisons in the tropical region only, where all the buoys are moored. The results show that in the case of shortwave radiation, the bias error from the ocean-based validation is somewhat larger than its land-based counterpart, but their RMS errors are not significantly different. In the case of longwave radiation, the ocean-based and land-based validations show nearly identical comparison statistics.

Taiping Zhang↗
NASA NTRS1998

New AFM Techniques for Investigating Molecular Growth Mechanisms of Protein Crystals

Atomic Force Microscopy (AFM) has emerged as a powerful technique for investigating protein crystal growth. Earlier AFM studies were among the first to demonstrate that these crystals grew by dislocation and 2D nucleation growth mechanisms [1]. These investigations were restricted to the micron range where only surface features, such as dislocation hillocks and 2D islands are visible. Most AFM instruments can scan at higher resolutions and have the potential to resolve individual protein molecules at nanometer ranges. Such scans are essential for determining the molecular packing arrangements on crystal faces and for probing the growth process at the molecular level. However, at this resolution the AFM tip influences the image produced, with the resulting image being a convolution of the tip shape and the surface morphology [2]. In most studies this problem is resolved by deconvoluting the image to obtain the true surface morphology. Although deconvolution routines work reasonably well for simple one- dimensional shapes, for complex surfaces this approach does not produce accurate results. In this study we devised a new approach which takes advantage of the precise molecular order of crystal surfaces, combined with the knowledge of individual molecular shapes from the crystallographic data of the protein and the AFM tip shape. This information is used to construct expected theoretical AFM images by convoluting the tip shape with the constructed crystal surface shape for a given surface packing arrangement. By comparing the images from actual AFM scans with the constructed ones for different possible surface packing arrangements, the correct packing arrangement can be conclusively determined. This approach was used in this study to determine the correct one from two possible packing arrangements on (I 10) faces of tetragonal lysozyme crystals. Another novel AFM technique was also devised to measure the dimension of individual growth units of the crystal faces. Measuring these units was not attempted before and most studies have assumed that the growth unit consisted of individual protein molecules. The linescan mode of AFM instruments allows the crystal surface to be scanned along a single line. By scanning across a growth step an image showing the motion of the step is obtained. Normally such an image shows a straight line for continuous and constant step velocity. In this study by increasing the scan rate and by decreasing the step velocity (by decreasing the supersaturation), we were able to capture images of individual growth events, shown by jump discontinuities in the step line. By suitable integration of the image the growth unit dimension in the scanned direction can be obtained. Since multiple units can be involved in the growth process it is necessary to collect a statistically relevant sample before drawing conclusions about the growth mechanism. This technique was successfully employed to obtain the dimensions of growth units for the (110) face, showing that they consisted of various aggregates corresponding to the 43 helices in the crystal structure.

Li, Huayu↗
NASA NTRS2003

A New Objective Technique for Verifying Mesoscale Numerical Weather Prediction Models

This report presents a new objective technique to verify predictions of the sea-breeze phenomenon over east-central Florida by the Regional Atmospheric Modeling System (RAMS) mesoscale numerical weather prediction (NWP) model. The Contour Error Map (CEM) technique identifies sea-breeze transition times in objectively-analyzed grids of observed and forecast wind, verifies the forecast sea-breeze transition times against the observed times, and computes the mean post-sea breeze wind direction and speed to compare the observed and forecast winds behind the sea-breeze front. The CEM technique is superior to traditional objective verification techniques and previously-used subjective verification methodologies because: It is automated, requiring little manual intervention, It accounts for both spatial and temporal scales and variations, It accurately identifies and verifies the sea-breeze transition times, and It provides verification contour maps and simple statistical parameters for easy interpretation. The CEM uses a parallel lowpass boxcar filter and a high-order bandpass filter to identify the sea-breeze transition times in the observed and model grid points. Once the transition times are identified, CEM fits a Gaussian histogram function to the actual histogram of transition time differences between the model and observations. The fitted parameters of the Gaussian function subsequently explain the timing bias and variance of the timing differences across the valid comparison domain. Once the transition times are all identified at each grid point, the CEM computes the mean wind direction and speed during the remainder of the day for all times and grid points after the sea-breeze transition time. The CEM technique performed quite well when compared to independent meteorological assessments of the sea-breeze transition times and results from a previously published subjective evaluation. The algorithm correctly identified a forecast or observed sea-breeze occurrence or absence 93% of the time during the two- month evaluation period from July and August 2000. Nearly all failures in CEM were the result of complex precipitation features (observed or forecast) that contaminated the wind field, resulting in a false identification of a sea-breeze transition. A qualitative comparison between the CEM timing errors and the subjectively determined observed and forecast transition times indicate that the algorithm performed very well overall. Most discrepancies between the CEM results and the subjective analysis were again caused by observed or forecast areas of precipitation that led to complex wind patterns. The CEM also failed on a day when the observed sea- breeze transition affected only a very small portion of the verification domain. Based on the results of CEM, the RAMS tended to predict the onset and movement of the sea-breeze transition too early and/or quickly. The domain-wide timing biases provided by CEM indicated an early bias on 30 out of 37 days when both an observed and forecast sea breeze occurred over the same portions of the analysis domain. These results are consistent with previous subjective verifications of the RAMS sea breeze predictions. A comparison of the mean post-sea breeze winds indicate that RAMS has a positive wind-speed bias for .all days, which is also consistent with the early bias in the sea-breeze transition time since the higher wind speeds resulted in a faster inland penetration of the sea breeze compared to reality.

Case, Jonathan L.↗
NASA NTRSDate not supplied

Parameters Impacting Columnated Granular Soil Pneumatic Seal Performance

The ability of a column of loose granular soil to form a pneumatic seal was investigated by varying the diameter of the soil column, the effective column height, and the level of compaction in the soil. The soil column diameter was tested at three levels using pipes with inner diameters measuring 5.08, 10.16, and 15.24 centimeters (2, 4, and 6 inches). Soil was filled in each pipe to form 15, 30, and 45 centimeter (6, 12, and 18 inch) tall soil columns. Each diameter/height configuration was also tested at three levels of soil compaction, compared by calculating the bulk density of the soil with mass and volume measurements. GRC-1a simulant was used, with approximate low/medium/high bulk densities of 1.6, 1.75, and 1.9 g/cc achieved with a combination of vibration and tamping. The order of tests for a given column diameter was randomized and repeated three times. With the top of the soil open to atmosphere at room temperature, compressed air was injected through a small diffuser at the column base with several small downward-facing holes. The number and size of these holes was scaled such that a constant total inlet orifice area to column cross sectional area ratio was maintained for each column diameter. Inlet air pressure was slowly increased via a precision regulator to preserve quasi-static equilibrium in the soil column to minimize the impact of dynamics. Air pressure was increased all the way through the static and bubbling regimes until slugging or turbulent behavior was observed in the soil to ensure that the entire static regime had been captured during data collection. A three-factor, three-level analysis of variance (ANOVA) statistical analysis was performed on the resulting data to determine the extent to which each physical parameter impacted the soil column seal performance. It was concluded that there is statistically significant evidence that column height, and the interaction between column height and diameter impact soil seal performance. In all other cases there was insufficient data to identify a statistically significant causal relationship. Additionally, plots were generated comparing experimental data to the predictive formula developed by Ogino et al. for fluidized beds in 1993. Because this model was developed for industrial spouted fluidized beds, the accuracy of its output prior to fluidization in the static seal ‘edge case’ is unknown, especially considering in this application the working gas was diffused across the column base rather than being injected through a spout. Further, the fidelity of the model had not yet been tested with lunar soil simulants. Plotting the Ogino et al. model alongside test data allows for a more intuitive sense of the impact of test parameters on soil seal performance, as well as providing a quick means to further tune this predictive model for more accurate use with static seals across granular lunar soil simulants. The model provided by Ogino et al. was further tuned using test data to determine the degree to which a spouted bed model could be applied to a slightly modified set of testing conditions: lunar soil simulants and a more diffuse, homogeneous application of pneumatic pressure. A Matlab script was created to test different values for the leading coefficient and exponents in Ogino’s formula. The script swept preset ranges, then iterated with higher resolutions across narrower ranges to converge on the optimal value for each parameter. The resulting adjusted model was compared to the original, as well as test data with noticeable improvements across the entire test domain.

Jack Stewart↗
NASA NTRS1999

The Enigmatic Longevity of Granular Materials on Mars: The Case for Geologically Episodic Dune Formation

Martian sand dunes are concentrated in vast sand seas in the circumpolar belt of the planet's northern hemisphere, but they are also pervasive over the whole planet. Their occurrence is to be expected on a super-arid planetary surface subjected to boundary layer drag from a continually active atmosphere. Whilst their occurrence is to be expected, their survival is enigmatic. But the enigma only arises if the martian system is considered similar to Earth's --where sand is moved highly frequently, more or less on a seasonal basis. Experimentally it is readily demonstrated that active sand will soon wear down to small grains and eventually diminish to below the critical sand size required to sustain dune formation. According to conventional wisdom, sand moves at higher speeds on Mars than on Earth, and if it were to move as frequently as it does on Earth, then the dune-forming sand population should have long since disappeared, given the great longevity of the martian aeolian system (Sagan coined the term "kamikaze" grains to express this disappearance). No supply of sand could keep pace with this depletion, especially in light of the fact that Mars does not have very active weathering, nor significant crustal differentiation. On Earth, plate tectonics, magmatic activity, and general crustal differentiation over geological time have produced great concentrations of quartz crystals in the continental crustal masses. Not only are these quartz grains chemically and mechanically resilient, they are about the right size for being transported by either wind or water. Add to this, the geologically recent contribution of glacial grinding, and it is easy to see why there are dune field on Earth. So what are the martian dunes composed of, and how does the material survive the eons of attrition? In addition to experimental demonstrations of sand comminution in laboratory aeolian simulations, the problem can be approached from first principles. Sagan showed that by simple considerations of material strength versus mechanical work applied to the material, comminution to sub-sand size would be inevitable. Another semi-analytical approach might be taken by considering that the archetypal aeolian sand surface texture is an irregularly pitted ("frosted") surface composed of chipping hollows approximately 10 microns in diameter, 5 microns deep. Their volume = about 250 cubic microns, or about 1/25000 of the volume of a 100 micron diameter dune grain. Because a saltating grain always strikes another grain, then two surfaces are impacted. Thus each grain undergoes two impacts for every one saltation leap, when the impact statistics are considered for a closed dune system (it can be calculated that a grain can never undergo <1 impact, and never >2 per saltation leap). Hence, if we conservatively assume that there is damage to a grain each time it bounces, but with the minimum damage of only 2 microscopic craters per impact, then approximately 12,500 impacts are required to completely eliminate the grain. Of course, it would require only a fraction of this amount to reduce the grain to below sand size. A grain will make only several tens of saltation leaps on the stoss side of a dune before becoming buried on the lee slope. The dune then has to move its full length before the grain is exhumed again for abrasion. Even with this hiatus in transport, it is easy to see that terrestrial dunes need resupplying with sand in order to survive. In recent theoretical work it has been shown that martian aeolian transport may be initiated with high-speed grains, but this converts to a lower energy dynamic transport equilibrium in which a reptation population dominates grain transport (on Earth, at least half of the flux is by reptation and creep). On Mars, therefore, average grain speeds may be lower than those on Earth, or at least comparable. This would permit greater longevity for martian sands, but it would not go far enough to solve the survival problem. It may, however, explain why martian dunes are about the same size as terrestrial dunes. If martian saltation leaps were significantly longer than on Earth (as usually assumed), then a dune's lee slope would have to be correspondingly longer in order to trap the sand; this would scale up the whole dune structure. But with shorter trajectories in a reptation population, larger dunes would be unnecessary. Additional information is contained in the original.

Marshall, J.↗
NASA NTRSDate not supplied

Examination of L-Band Brightness Temperature Forecasts in the SMAP Level-4 Soil Moisture Analysis

The NASA Soil Moisture Active Passive (SMAP) mission [4] has been providing L-band (1.4 GHz) passive microwave brightness temperature (Tb) observations since April 2015. By assimilating the Tb observations into the NASA Catchment land surface model [5] using a spatially distributed ensemble Kalman filter, the NASA Global Modeling and Assimilation Office generates the SMAP Level-4 Soil Moisture (L4_SM) product, which provides global, 3-hourly, 9-km resolution estimates of surface (0-5 cm) and root-zone (0-100 cm) soil moisture with ~2.5-day latency for use in research and applications [6]. The L4_SM product also includes estimates of soil temperature, land surface fluxes, and assimilation diagnostics such as the model forecast and observed Tb values [7]. The output from the L4_SM system is routinely monitored by the L4_SM team. Such monitoring provides valuable information; instances of unusually large Tb observation-minus-forecast (O-F) residuals can indicate events during which soil moisture conditions are poorly described in the land modeling system [7,8]. For example, repeated occurrences of very large Tb O-F values in central Australia were traced back to deficiencies in the gauge-based precipitation product used in the land modeling system through L4_SM Version 5; this discovery prompted the use, beginning in Version 6, of satellite- and gauge-based precipitation observations outside of North America [9]. Recently, a statistical analysis of the Tb O-F residuals revealed systematic errors in the “tau-omega” L-band radiative transfer model that converts the land-model simulated soil moisture and temperature into the forecast Tb prior to the L4_SM analysis. The L4_SM analysis is built on the ensemble Kalman filter and assumes unbiased forecast errors. The presence of systematic Tb forecast errors could thus adversely impact the quality of the analyzed soil moisture. In this paper, we examine the Tb O-F residuals of the latest Version 7 L4_SM data (Science Version ID Vv7030 and Vv7032) [10,11,12]. In the L4_SM Version 7 algorithm, key parameters of the L-band radiative transfer model, including soil roughness, scattering albedo, and a (seasonally varying) climatology of vegetation opacity, are obtained from the SMAP Level-2 Radiometer retrieval product (Version 5) [1]. Additionally, we also use ground measurements of surface soil moisture and soil temperature from the SMAP core validation sites as in situ reference of soil conditions [2,3]. As part of the L4_SM system calibration, the seasonally-varying bias between the model forecast Tb and the observed values is removed prior to the assimilation of the SMAP Tb observations [7]. That is, the L4_SM system is designed to only correct errors in synoptic-scale and interannual variations from the long-term mean seasonal cycle while maintaining the model’s (potentially erroneous) climatology. The Tb O-F residuals examined here are thus computed using SMAP Tb observations after they are rescaled to the mean seasonal cycle of Tb from the modeling system without SMAP data assimilation. Consequently, the long-term average of the Tb O-F residuals has a global mean of only 0.13 K and locally small values, ranging from -1 to 3 K. Despite the small time-average values of the Tb O-F residuals, the model forecast Tb was nevertheless found to exhibit undesirable systematic errors. At some locations, the time-average Tb O-F values strongly depend on surface soil moisture (SM). At the Yanco SMAP core validation site, for example, the Tb O-F residuals typically range from 5 to 15 K under dry soil moisture conditions (SM < 0.15 m3 m-3) yet are predominantly negative under wet soil moisture conditions (SM > 0.25 m3 m-3), with values ranging from 0 to -40 K. This results in soil moisture analysis increments that persistently make the soil drier under dry SM conditions and persistently make the soil wetter under wet SM conditions, suggesting an error in the dynamic range of the simulated Tb, soil moisture or soil temperature. In this paper, we describe the higher-order systematic Tb forecast errors in more detail, examine their impact on the L4_SM product quality, and explore potential avenues to further improve the L4_SM algorithm.

Rolf Reichle↗
NASA NTRS2019

Shockwave/Boundary-Layer Interaction Studies Performed in the NASA Langley 20-Inch Mach 6 Air Tunnel

This paper highlights results from a collaborative study performed by The University of Tennessee Space Institute (UTSI) and NASA Langley Research Center on the Shockwave/Boundary-Layer Interaction (SWBLI) generated by a cylindrical protuberance on a flat plate in a Mach 6 flow. The study was performed in the 20-Inch Mach 6 Air Tunnel at NASA Langley Research Center and consisted of two separate entries. In the first entry, simultaneous high-speed schlieren and high-speed pressure-sensitive paint (PSP) imaging – which was performed for the first time in the 20-Inch Mach 6 facility at NASA Langley – were performed as well as simultaneous high-speed schlieren and oil-flow imaging. In the second entry, the model configuration was modified to increase the size of the interaction region. High-speed schlieren and infrared thermography (IR) surface imaging were performed in this second entry. The goal of these tests was to characterize the SBLI in the presence of a laminar, transitional, and turbulent boundary layer using high-speed optical imaging techniques. AoA = sting angle-of-attack (°) dcylinder = cylinder diameter (mm) dtrip = cylindrical tripping element diameter (mm) Δshock = shock stand-off distance (mm) hcylinder = cylinder height (mm) htrip = cylindrical tripping element height (mm) HSS = high-speed schlieren M∞ = freestream Mach number PSP = pressure-sensitive paint Re∞ = freestream unit Reynolds number (m-1) SWBLI = shockwave/boundary-layer interaction θplate = model plate angle (°) Introduction his paper highlights two experimental entries performed in the 20-Inch Mach 6 Air Blowdown Tunnel at NASA Langley Research Center in collaboration with The University of Tennessee Space Institute (UTSI). The purpose of these entries was to characterize the dynamic shockwave/boundary-layer interaction (SWBLI) between a vertical cylinder on a flat plate and laminar, transitional (XSWBLI), and turbulent (SWTBLI) boundary layers with a freestream Mach number of 6 using non-intrusive optical diagnostics. Experiments performed by Murphree et al.1,2 were among the first to specifically characterize XSWBLI induced by a vertical cylinder on a flat plate geometry using several optical measurement techniques. Recent optical studies of XSWBLI phenomenon have been performed by UTSI at Mach 2 in their low-enthalpy blow wind tunnel3-8 and by Texas A&M University and UTSI at Mach numbers of 6 and 7 in their Adjustable Contour Expansion wind tunnel.9 The experiments described in this paper were intended to complement previous studies by expanding the freestream unit Reynolds number range, Re∞, over which the XSWBLI phenomena has been observed. Additionally these experiments, made possible under NASA’s new facility funding model under the Aeronautics Evaluation and Test Capabilities (AETC) project, promoted collaboration between university and NASA researchers. The initial entry in the 20-Inch Mach 6 Air Tunnel at NASA Langley occurred in December of 2016. Originally, testing was to occur in November of 2016 in the 31-Inch Mach 10 Air Tunnel at NASA Langley. This facility was chosen so that the XSWBLI phenomenon could be observed at much higher Mach numbers than had previously been attempted in ground test experiments. The model selected for this experiment, a 10° half-angle wedge with a sharp leading edge (described in detail in section II.B), had previously been used by Danehy et al. [10] for boundary layer transition studies using the nitric oxide planar laser-induced fluorescence (NO PLIF) flow visualization technique. In that work, it was determined that transition could be induced downstream of a single htrip = 1-mm tall, dtrip = 4-mm diameter cylindrical tripping element and that the streamwise location of the transition could be changed for a single Re∞ by changing the model angle-of-attack (AoA) (see Fig. A3 in Ref. [10] for more details). Based on the findings of that work, a decision was made to use the wedge model with the cylindrical tripping element to trip the boundary layer flow ahead of a cylindrical protuberance in order to achieve a XSWBLI. Unfortunately, the 31-Inch Mach 10 facility had been taken offline for repairs in October of 2016 and a decision was made to move the test to the 20-Inch Mach 6 facility. Since the behavior of the boundary layer with the chosen model configuration had not been studied before in that facility and the available test time was limited, the entry was considered to be exploratory and was used to collect spatially-resolved and time-resolved flow and surface visualization data that would be used to inform a second entry. Test techniques included simultaneous high-speed schlieren (HSS) captured at 160 kHz and high-speed pressure sensitive paint captured at 10 kHz as well as oil flow visualization, captured at 750 Hz. The second entry in the 20-Inch Mach 6 facility occurred in June and July of 2017. In this follow-on test, modifications to the wind tunnel model were made based on observations made during the first entry and included removing the cylindrical tripping element, increasing the size of the cylinder used to induce the SWBLI to increase the size of the interaction while simultaneously improving spatial resolution, and using a swept ramp array, similar to that described in Ref. [11], to trip the flow to turbulence. Simultaneous HSS (captured at 140 kHz, 100 kHz, and 40 kHz) and conventional IR thermography (captured at 30 Hz) imaging were performed simultaneously in this follow-on entry. This paper is intended to serve as a summary of the work performed during these two entries, to detail lessons learned from each entry, and to highlight some of the datasets acquired. Details on the experimental setup, model configuration, and techniques used are provided. Papers providing a more rigorous analysis of data acquired during the second entry, including statistical, spectral, and modal decomposition methods, can be found in Refs. [12,13]. An entry examining XSWBLI in the 31-Inch Mach 10 Blowdown Wind Tunnel facility is currently planned for mid-to-late calendar year 2019, pending the success of facility repairs. The work performed and described in this paper and the upcoming entry in the 31-Inch Mach 10 facility at NASA Langley have been made possible by NASA’s new facility funding model under the Aeronautics Evaluation and Test Capabilities (AETC) project. Wind Tunnel Facility All experiments discussed in this paper were performed in the 20-Inch Mach 6 Air Tunnel at NASA Langley Research Center. Specific details pertaining to this facility can be found in Refs. [14,15], with only a brief description of the facility provided here. For both entries, the nominal freestream unit Reynolds number was varied between 1.8×106 m-1 (0.5×106 ft-1) and 26.3×106 m-1 (8×106 ft-1). The nominal stagnation pressure was varied between 0.21 MPa and 3.33 MPa and the nominal stagnation temperature was varied between 480 K and 520 K to achieve the desired Re∞ condition. For all runs, the nominal freestream Mach number was 6. The nearly square test section is 520.7-mm (20.5-inches) wide by 508-mm (20-inches) high. Two 431.8-mm (17-inch) diameter windows made of Corning 7940, Grade 5F schlieren-quality glass serve as the side walls of the tunnel and provide optical access for the high-speed schlieren measurements. A rectangular window made of the same material as the side windows served as the top wall of the test section and provided optical access for the high-speed PSP and oil flow measurements. For the second entry, this top window was replaced with a Zinc Selenide (ZnSe) window with an anti-reflection coating capable of passing IR wavelengths between 8μm and 12μm with greater than 98% transmittance. The model was sting supported by a strut attached to a hydraulic system that allows for the model pitch angle to be adjusted between -5° to +55°. For the first entry, an initial pitch/pause sweep of the model AoA was performed to observe the resulting SWBLI. Ultimately, however, the sting pitch angle for this entry was fixed at +10.0° so that the angle of the top surface of the wedge relative to the streamwise axis of the tunnel (referred to herein as the plate angle, θplate), was θplate = 0°. For the second entry, θplate = 0° and θplate = -13.25° were initially tested with the swept ramp array (discussed in the following section) to determine which orientation produced conditions most favorable for XSWBLI to occur based on the heating signatures observed over the top surface of the model in the IR thermography images. Based on these initial tests, θplate = -13.25° was set for the remainder of the runs in the second entry. For both entries, any model changes were performed in a housing located beneath the closed test section. Prior to performing a run of the tunnel, the housing was sealed and the tunnel started. Once the appropriate freestream conditions were achieved, the model was injected into the test section using a hydraulic injection system. B. Model Geometry For all runs, a 10° half-angle (20° full-angle) wedge model with a sharp leading edge was used. The model is described in detail in Refs. [10,16]. The top surface of the sharp leading edge of the model extended 47.8 mm from its upstream-most edge to a junction with the upstream edge of a stainless steel top plate that then extended an (a) (c) (b) Fig. 1 (a) Schematic of top surface of wedge model with gas seeding insert, (b) perspective view of the model in the 20-Inch Mach 6 tunnel with centerline pressure orifices on sharp leading edge, and (c) a perspective view of the model with stainless steel (top) and SLA middle insert (bottom) during the first entry. Flow occurs from left to right.

Bathel, Brett F.↗
NASA NTRS2017

Measuring Shock Stage of ltokawa Regolith Grains by Electron Back-Scattered Diffraction and Synchrotron X-Ray Diffraction

We have been analyzing Itokawa samples in order to definitively establish the degree of shock experienced by the regolith of asteroid Itokawa, and to devise a bridge between shock determinations by standard light optical petrography, crystal structures as determined by electron and X-ray diffraction techniques. We are making measurements of olivine crystal structures and using these to elucidate critical regolith impact processes. We use electron back-scattered diffraction (EBSD) and synchrotron X-ray diffraction (SXRD). We are comparing the Itokawa samples to L and LL chondrite meteorites chosen to span the shock scale experienced by Itokawa, specifically Chainpur (LL3.4, Shock Stage 1), Semarkona (LL3.00, S2), Kilabo (LL6, S3), NWA100 (L6, S4) and Chelyabinsk (LL5, S4). In SXRD we measure the line broadening of olivine reflections as a measure of shock stage. In this presentation we concentrate on the EBSD work. We employed JSC's Supra 55 variable pressure FEG-SEM and Bruker EBSD system. We are not seeking actual strain values, but rather indirect strain-related measurements such as extent of intra-grain lattice rotation, and determining whether shock state "standards" (meteorite samples of accepted shock state, and appropriate small grain size) show strain measurements that may be statistically differentiated, using a sampling of particles (number and size range) typical of asteroid regoliths. Using our system we determined that a column pressure of 9 Pa and no C-coating on the sample was optimal. We varied camera exposure time and gain to optimize mapping performance, concluding that 320x240 pattern pixilation, frame averaging of 3, 15 kV, and low extractor voltage yielded an acceptable balance of hit rate (>90%), speed (11 fps) and map quality using an exposure time of 30 ms (gain 650). We found that there was no strong effect of step size on Grain Orientation Spread (GOS) and Grain Reference Orientation Deviation angle (GROD-a) distribution; there was some effect on grain average Kernel Average Misorientation (KAM) (reduced with smaller step size for the same grain), as expected. We monitored GOS, Maximum Orientation Spread (MOS) and GROD-a differences between whole olivine grains and sub-sampled areas, and found that there were significant differences between the whole grain dataset and subsets, as well as between subsets, likely due to sampling-related "noise". Also, in general (and logically) whole grains exhibit greater degrees of cumulative lattice rotation. Sampling size affects the apparent strain character of the grain, at least as measured by GOS, MOS and GROD-a. There were differences in the distribution frequencies of GOS and MOS between shock stages, and in plots of MOS and GOS vs. grain diameter. These results are generally consistent with those reported this year. However, it is unknown whether the differences between samples of different shock states exceeds the clustering of these values to the extent that shock stage determinations can still be made with confidence. We are investigating this by examination of meteorites with higher shock stage 4 to 5. Our research will improve our understanding of how small, primitive solar system bodies formed and evolved, and improve understanding of the processes that determine the history and future of habitability of environments on other solar system bodies. The results will directly enrich the ongoing asteroid and comet exploration missions by NASA and JAXA, and broaden our understanding of the origin and evolution of small bodies in the early solar system, and elucidate the nature of asteroid and comet regolith.

Zolensky, Michael↗