Search NASA⌕ Search

SEARCH · Search NASA

Results for “Scaling and squaring method”

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 181 records · Page 10

Winter Wheat Yield Assessment from Landsat 8 and Sentinel-2 Data: Incorporating Surface Reflectance, Through Phenological Fitting, into Regression Yield Models

A combination of Landsat 8 and Sentinel-2 offers a high frequency of observations (3–5 days) at moderate spatial resolution (10–30 m), which is essential for crop yield studies. Existing methods traditionally apply vegetation indices (VIs) that incorporate surface reflectances (SRs) in two or more spectral bands into a single variable, and rarely address the incorporation of SRs into empirical regression models of crop yield. In this work, we address these issues by normalizing satellite data (both VIs and SRs) derived from NASA’s Harmonized Landsat Sentinel-2 (HLS) product, through a phenological fitting. We apply a quadratic function to fit VIs or SRs against accumulated growing degree days (AGDDs), which affects the rate of crop development. The derived phenological metrics for VIs and SRs, namely peak, area under curve (AUC), and fitting coefficients from a quadratic function, were used to build empirical regression winter wheat models at a regional scale in Ukraine for three years, 2016–2018. The best results were achieved for the model with near infrared (NIR) and red spectral bands and derived AUC, constant, linear, and quadratic coefficients of the quadratic model. The best model yielded a root mean square error (RMSE) of 0.201 t/ha (5.4%) and coefficient of determination R2 = 0.73 on cross-validation.

phenological fitting↗

Test-Anchored Vibration Response Predictions for an Acoustically Energized Curved Orthogrid Panel with Mounted Components

A rich body of vibroacoustic test data was recently generated at Marshall Space Flight Center for a curved orthogrid panel typical of launch vehicle skin structures. Several test article configurations were produced by adding component equipment of differing weights to the flight-like vehicle panel. The test data were used to anchor computational predictions of a variety of spatially distributed responses including acceleration, strain and component interface force. Transfer functions relating the responses to the input pressure field were generated from finite element based modal solutions and test-derived damping estimates. A diffuse acoustic field model was employed to describe the assumed correlation of phased input sound pressures across the energized panel. This application demonstrates the ability to quickly and accurately predict a variety of responses to acoustically energized skin panels with mounted components. Favorable comparisons between the measured and predicted responses were established. The validated models were used to examine vibration response sensitivities to relevant modeling parameters such as pressure patch density, mesh density, weight of the mounted component and model form. Convergence metrics include spectral densities and cumulative root-mean squared (RMS) functions for acceleration, velocity, displacement, strain and interface force. Minimum frequencies for response convergence were established as well as recommendations for modeling techniques, particularly in the early stages of a component design when accurate structural vibration requirements are needed relatively quickly. The results were compared with long-established guidelines for modeling accuracy of component-loaded panels. A theoretical basis for the Response/Pressure Transfer Function (RPTF) approach provides insight into trends observed in the response predictions and confirmed in the test data. The software modules developed for the RPTF method can be easily adapted for quick replacement of the diffuse acoustic field with other pressure field models; for example a turbulent boundary layer (TBL) model suitable for vehicle ascent. Wind tunnel tests have been proposed to anchor the predictions and provide new insight into modeling approaches for this type of environment. Finally, component vibration environments for design were developed from the measured and predicted responses and compared with those derived from traditional techniques such as Barrett scaling methods for unloaded and component-loaded panels.

Frady, Gregory P.↗

Test-Anchored Vibration Response Predictions for an Acoustically Energized Curved Orthogrid Panel with Mounted Components

rich body of vibroacoustic test data was recently generated at Marshall Space Flight Center for component-loaded curved orthogrid panels typical of launch vehicle skin structures. The test data were used to anchor computational predictions of a variety of spatially distributed responses including acceleration, strain and component interface force. Transfer functions relating the responses to the input pressure field were generated from finite element based modal solutions and test-derived damping estimates. A diffuse acoustic field model was applied to correlate the measured input sound pressures across the energized panel. This application quantifies the ability to quickly and accurately predict a variety of responses to acoustically energized skin panels with mounted components. Favorable comparisons between the measured and predicted responses were established. The validated models were used to examine vibration response sensitivities to relevant modeling parameters such as pressure patch density, mesh density, weight of the mounted component and model form. Convergence metrics include spectral densities and cumulative root-mean squared (RMS) functions for acceleration, velocity, displacement, strain and interface force. Minimum frequencies for response convergence were established as well as recommendations for modeling techniques, particularly in the early stages of a component design when accurate structural vibration requirements are needed relatively quickly. The results were compared with long-established guidelines for modeling accuracy of component-loaded panels. A theoretical basis for the Response/Pressure Transfer Function (RPTF) approach provides insight into trends observed in the response predictions and confirmed in the test data. The software developed for the RPTF method allows easy replacement of the diffuse acoustic field with other pressure fields such as a turbulent boundary layer (TBL) model suitable for vehicle ascent. Structural responses using a TBL model were demonstrated, and wind tunnel tests have been proposed to anchor the predictions and provide new insight into modeling approaches for this environment. Finally, design load factors were developed from the measured and predicted responses and compared with those derived from traditional techniques such as historical Mass Acceleration Curves and Barrett scaling methods for acreage and component-loaded panels.

Frady, Gregory P.↗

Validation Experiments of Incipient Turbulent Separation over an Axisymmetric Afterbody

Historically, the flow physics involved with most turbulent separated flows have presented fundamental challenges in validation between experimental and numerical approaches. As recognized by the CFD Vision 2030 study commissioned by NASA, validation of Reynolds-averaged Navier-Stokes (RANS) models and other scale-resolving methods for turbulent separated flow requires the support of advanced, high-fidelity experiments designed specifically for CFD implementation. In accordance with this effort, a new test platform, referred to as the NASA Axisymmetric Afterbody, was designed to obtain detailed information of the flow field over a smooth body, undergoing adverse pressure gradient induced separation for a fixed Reynolds number, Re = 180,000. The parametric body offers a range of flow states from fully attached,incipient separation, to fully separated flow based on variable afterbody geometries. In an initial effort to evaluate RANS turbulence model capabilities, the present configuration of the axisymmetric model features a mild adverse pressure gradient over the contoured boattail section, inducing incipient turbulent separation. Experiments include steady pressure measurements to serve as preliminary comparisons to simulation studies, which examine the effect of variable grid domains and RANS turbulence models, in an effort to understand and evaluate the critical variability between solutions for the present configuration. Results indicate potential discrepancies may be due to the effect of the square tunnel test section walls, relatively large blockage ratio, and inconsistent reference parameters. Ongoing work will focus on the experimental campaigns to obtain surface flow visualizations and high-resolution, off-body flow field measurements using Stereoscopic Particle Image Velocimetry (SPIV) and Laser Doppler Velocimetry (LDV) techniques to provide a detailed benchmark dataset to aid turbulence modelers.

axisymmetric afterbody↗

Image Correlation Pattern Optimization for Micro-Scale In-Situ Strain Measurements

The accuracy and precision of digital image correlation (DIC) is a function of three primary ingredients: image acquisition, image analysis, and the subject of the image. Development of the first two (i.e. image acquisition techniques and image correlation algorithms) has led to widespread use of DIC; however, fewer developments have been focused on the third ingredient. Typically, subjects of DIC images are mechanical specimens with either a natural surface pattern or a pattern applied to the surface. Research in the area of DIC patterns has primarily been aimed at identifying which surface patterns are best suited for DIC, by comparing patterns to each other. Because the easiest and most widespread methods of applying patterns have a high degree of randomness associated with them (e.g., airbrush, spray paint, particle decoration, etc.), less effort has been spent on exact construction of ideal patterns. With the development of patterning techniques such as microstamping and lithography, patterns can be applied to a specimen pixel by pixel from a patterned image. In these cases, especially because the patterns are reused many times, an optimal pattern is sought such that error introduced into DIC from the pattern is minimized. DIC consists of tracking the motion of an array of nodes from a reference image to a deformed image. Every pixel in the images has an associated intensity (grayscale) value, with discretization depending on the bit depth of the image. Because individual pixel matching by intensity value yields a non-unique scale-dependent problem, subsets around each node are used for identification. A correlation criteria is used to find the best match of a particular subset of a reference image within a deformed image. The reader is referred to references for enumerations of typical correlation criteria. As illustrated by Schreier and Sutton and Lu and Cary systematic errors can be introduced by representing the underlying deformation with under-matched shape functions. An important implication, as discussed by Sutton et al., is that in the presence of highly localized deformations (e.g., crack fronts), error can be reduced by minimizing the subset size. In other words, smaller subsets allow the more accurate resolution of localized deformations. Contrarily, the choice of optimal subset size has been widely studied and a general consensus is that larger subsets with more information content are less prone to random error. Thus, an optimal subset size balances the systematic error from under matched deformations with random error from measurement noise. The alternative approach pursued in the current work is to choose a small subset size and optimize the information content within (i.e., optimizing an applied DIC pattern), rather than finding an optimal subset size. In the literature, many pattern quality metrics have been proposed, e.g., sum of square intensity gradient (SSSIG), mean subset fluctuation, gray level co-occurrence, autocorrelation-based metrics, and speckle-based metrics. The majority of these metrics were developed to quantify the quality of common pseudo-random patterns after they have been applied, and were not created with the intent of pattern generation. As such, it is found that none of the metrics examined in this study are fit to be the objective function of a pattern generation optimization. In some cases, such as with speckle-based metrics, application to pixel by pixel patterns is ill-conditioned and requires somewhat arbitrary extensions. In other cases, such as with the SSSIG, it is shown that trivial solutions exist for the optimum of the metric which are ill-suited for DIC (such as a checkerboard pattern). In the current work, a multi-metric optimization method is proposed whereby quality is viewed as a combination of individual quality metrics. Specifically, SSSIG and two auto-correlation metrics are used which have generally competitive objectives. Thus, each metric could be viewed as a constraint imposed upon the others, thereby precluding the achievement of their trivial solutions. In this way, optimization produces a pattern which balances the benefits of multiple quality metrics. The resulting pattern, along with randomly generated patterns, is subjected to numerical deformations and analyzed with DIC software. The optimal pattern is shown to outperform randomly generated patterns.

Bomarito, G. F.↗

Numerical Computations of Hypersonic Boundary-Layer over Surface Irregularities

Surface irregularities such as protuberances inside a hypersonic boundary layer may lead to premature transition on the vehicle surface. Early transition in turn causes large localized surface heating that could damage the thermal protection system. Experimental measurements as well as numerical computations aimed at building a knowledge base for transition Reynolds numbers with respect to different protuberance sizes and locations have been actively pursued in recent years. This paper computationally investigates the unsteady wake development behind large isolated cylindrical roughness elements and the scaled wind-tunnel model of the trip used in a recent flight measurement during the reentry of space shuttle Discovery. An unstructured mesh, compressible flow solver based on the space-time conservation element, solution element (CESE) method is used to perform time-accurate Navier-Stokes calculations for the flow past a roughness element under several wind-tunnel conditions. For a cylindrical roughness element with a height to the boundary-layer thickness ratio from 0.8 to 2.5, the wake flow is characterized by a mushroom-shaped centerline streak and horse-shoe vortices. While time-accurate solutions converged to a steady-state for a ratio of 0.8, strong flow unsteadiness is present for a ratio of 1.3 and 2.5. Instability waves marked by distinct disturbance frequencies were found in the latter two cases. Both the centerline streak and the horse-shoe vortices become unstable downstream. The oscillatory vortices eventually reach an early breakdown stage for the largest roughness element. Spectral analyses in conjunction with the computed root mean square variations suggest that the source of the unsteadiness and instability waves in the wake region may be traced back to possible absolute instability in the front-side separation region.

Chang, Chau-Lyan↗

Application of Particle Image Velocimetry to a Study of Flow About a Multi-Element Airfoil

An experiment was performed on the flap tip vortex shed from a half span Fowler flap. This flap was mounted on a 5 foot span NACA 63(2)-215 Mod B airfoil in the 7 by 10 foot wind tunnel at NASA Ames Research Center. Several noise reduction studies were performed with this model, and the addition of the Particle Image Velocimetry (PIV) research discussed here served as a proof case of large scale PIV. The measurement plane investigated here was a cross plane region. This is cross plane relative to the freestream flow direction. The measurement plane was located at a position 18 inches downstream of the flap trailing edge. This served to prove that measurements could also be made in the more difficult cross plane direction rather than in the downstream flow direction. Lastly the PIV data was used as a practical research tool that yielded important results that could not otherwise be obtained. The flow field area measured was 40 cm by 40 cm square, and served to characterize the downstream flow characteristics of the flap tip vortex under three configurations: the baseline configuration which was the flap and the wing only; the baseline with the addition of a 3/4 span slat; and the baseline with a Flap Edge Device which was designed to reduce the noise generated at the flap. All configurations were tested at a freestream velocity of 64.84 m/s. The test resulted in average velocity fields for the three configurations tested. The velocity fields aided in verifying other testing methods on this particular experiment, and also yielded further insight into the characteristics of the flap tip vortex under the three configurations considered. The velocity data was reduced, and we were able to calculate the vorticity of the flow field. From the position of minimum vorticity the location of the center of the vortex was determined. The circulation was also calculated and aided in comparing the effects of the three configurations on the lifting characteristics of the flap.

Walker, Stephen M.↗

SkICAT: A cataloging and analysis tool for wide field imaging surveys

We describe an integrated system, SkICAT (Sky Image Cataloging and Analysis Tool), for the automated reduction and analysis of the Palomar Observatory-ST ScI Digitized Sky Survey. The Survey will consist of the complete digitization of the photographic Second Palomar Observatory Sky Survey (POSS-II) in three bands, comprising nearly three Terabytes of pixel data. SkICAT applies a combination of existing packages, including FOCAS for basic image detection and measurement and SAS for database management, as well as custom software, to the task of managing this wealth of data. One of the most novel aspects of the system is its method of object classification. Using state-of-theart machine learning classification techniques (GID3* and O-BTree), we have developed a powerful method for automatically distinguishing point sources from non-point sources and artifacts, achieving comparably accurate discrimination a full magnitude fainter than in previous Schmidt plate surveys. The learning algorithms produce decision trees for classification by examining instances of objects classified by eye on both plate and higher quality CCD data. The same techniques will be applied to perform higher-level object classification (e.g., of galaxy morphology) in the near future. Another key feature of the system is the facility to integrate the catalogs from multiple plates (and portions thereof) to construct a single catalog of uniform calibration and quality down to the faintest limits of the survey. SkICAT also provides a variety of data analysis and exploration tools for the scientific utilization of the resulting catalogs. We include initial results of applying this system to measure the counts and distribution of galaxies in two bands down to Bj is approximately 21 mag over an approximate 70 square degree multi-plate field from POSS-II. SkICAT is constructed in a modular and general fashion and should be readily adaptable to other large-scale imaging surveys.

Weir, N.↗

Ares Launch Vehicle Transonic Buffet Testing and Analysis Techniques

It is necessary to define the launch vehicle buffet loads to ensure that structural components and vehicle subsystems possess adequate strength, stress, and fatigue margins when the vehicle structural dynamic response to buffet forcing functions are considered. In order to obtain these forcing functions, the accepted method is to perform wind-tunnel testing of a rigid model instrumented with hundreds of unsteady pressure transducers designed to measure the buffet environment across the desired frequency range. The buffet wind-tunnel test program for the Ares Crew Launch Vehicle employed 3.5 percent scale rigid models of the Ares I and Ares I-X launch vehicles instrumented with 256 unsteady pressure transducers each. These models were tested at transonic conditions at the Transonic Dynamics Tunnel at NASA Langley Research Center. The ultimate deliverable of the Ares buffet test program are buffet forcing functions (BFFs) derived from integrating the measured fluctuating pressures on the rigid wind-tunnel models. These BFFs are then used as input to a multi-mode structural analysis to determine the vehicle response to buffet and the resulting buffet loads and accelerations. This paper discusses the development of the Ares I and I-X rigid buffet model test programs from the standpoint of model design, instrumentation system design, test implementation, data analysis techniques to yield final products, and presents normalized sectional buffet forcing function root-mean-squared levels.

Piatak, David J.↗

Behavior of Forged S-816 Turbine Blades in Steady-state Operation of J33-9 Turbojet Engine with Stress-rupture and Metallographic Evaluations / by F. B.garrett, C. A. Gyorgak, and J. W. Weeton

An investigation was conducted to determine the behavior of recently produced, forged S-816 turbine blades in a full-scale turbojet engine, and in particular, the scatter in performance of the alloy. The turbine blades were operated as continuously as possible at a temperature of 1500 degrees F and a centrifugal stress of 21,500 pounds per square inch. The operating lives of the turbine blades varied from 181 to 539 hours, a range of 358 hours. Stress-rupture properties of specimens cut from blade airfoils also varied considerably, as much as 1257 hours at 20,000 pounds per square inch and 1500 degrees F. Since the variability of scatter of stress-rupture data is greater than that of blade performance, the scatter is probably caused by variations in the properties of the forged blades rather than by variations caused by engine operation or installation of the blades. Metallographic examinations were made to determine possible causes of the scatter and although numerous differences in microstructures of blades were found, no consistent tendencies were observed and the findings did not permit an explanation of the scatter of blade performance. The results of the metallographic examinations and of the physical tests indirectly indicated variables in the fabricating method caused the scatter in properties.

Garrett, Floyd B↗

Performance of Several Conical Convergent-Divergent Rocket-Type Exhaust Nozzles

An investigation was conducted to obtain nozzle performance data with relatively large-scale models at pressure ratios as high as 120. Conical convergent-divergent nozzles with divergence angles alpha of 15, 25, and 29 deg. were each tested at area ratios of approximately 10, 25, and 40. Heated air (1200 F) was supplied at the nozzle inlet at pressures up to 145 pounds per square inch absolute and was exhausted into quiescent air at pressures as low as 1.2 pounds per square inch absolute. Thrust ratios for all nozzle configurations are presented over the range of pressure ratios attainable and were extrapolated when possible to design pressure ratio and beyond. Design thrust ratios decreased with increasing nozzle divergence angle according to the trend predicted by the (1 + cos alpha)/2 parameter. Decreasing the nozzle divergence angle resulted in sizable increases in thrust ratio for a given surface-area ratio (nozzle weight), particularly at low nozzle pressure ratios. Correlations of the nozzle static pressure at separation and of the average static pressure downstream of separation with various nozzle parameters permitted the calculation of thrust in the separated-flow region from unseparated static-pressure distributions. Thrust ratios calculated by this method agreed with measured values within about 1 percent.

Campbell, C. E.↗

Interferometric detection of freeze-thaw displacements of Alaskan permafrost using ERS-1 data

The possibility of making large scale (50 km) measurements of motions of the earth's surface with high resolution (10 m) and very high accuracy (1 cm) from multipass SAR interferometry was established in 1989. Other experiments have confirmed the viability and usefulness of the method. Work is underway in various groups to measure displacements from volcanic activity, seismic events, glacier motion, and in the present study, freeze-thaw cycles in Alaskan permafrost. The ground is known to move significantly in these cycles, and provided that freezing does not cause image decorrelation, it should be possible to measure both ground swelling and subsidence. The authors have obtained data from multiple passes of ERS-1 over the Toolik Lake region of northern Alaska of suitable quality for interferometry. The data are processed into images, and single interferograms are formed in the usual manner. Phase unwrapping is performed, and the multipass baselines are estimated from the images using both orbit ephemerides and scene tie points. The phases are scaled by the baseline ratio, and a double-difference interferogram (DDI) is formed. It is found that there is a residual 'saddle-shape' phase error across the image, which is postulated to be caused by a small divergence (10(exp -2) deg.) in the orbits. A simulation of a DDI from divergent orbits confirms the shape and magnitude of the error. A two-dimensional least squares fit to the error is performed, which is used to correct the DDI. The final, corrected DDI shows significant phase (altitude) changes over the period of the observation.

Werner, Charles L.↗

Mapping Global Forest Age from Forest Inventories, Biomass and Climate Data

Forest age can determine the capacity of a forest to uptake carbon from the atmosphere. However, a lack of global diagnostics that reflect the forest stage and associated disturbance regimes hampers the quantification of age-related differences in forest carbon dynamics. This study provides a new global distribution of forest age circa 2010, estimated using a machine learning approach trained with more than 40 000 plots using forest inventory, biomass and climate data. First, an evaluation against the plot-level measurements of forest age reveals that the data-driven method has a relatively good predictive capacity of classifying old-growth vs. non-old-growth (precision = 0.81 and 0.99 for old-growth and non-old-growth, respectively) forests and estimating corresponding forest age estimates (NSE = 0.6 – Nash–Sutcliffe efficiency – and RMSE = 50 years – root-mean-square error). However, there are systematic biases of overestimation in young- and underestimation in old-forest stands, respectively. Globally, we find a large variability in forest age with the old-growth forests in the tropical regions of Amazon and Congo, young forests in China, and intermediate stands in Europe. Furthermore, we find that the regions with high rates of deforestation or forest degradation (e.g. the arc of deforestation in the Amazon) are composed mainly of younger stands. Assessment of forest age in the climate space shows that the old forests are either in cold and dry regions or warm and wet regions, while young–intermediate forests span a large climatic gradient. Finally, comparing the presented forest age estimates with a series of regional products reveals differences rooted in different approaches and different in situ observations and global-scale products. Despite showing robustness in cross-validation results, additional methodological insights on further developments should as much as possible harmonize data across the different approaches. The forest age dataset presented here provides additional insights into the global distribution of forest age to better understand the global dynamics in the forest water and carbon cycles. The forest age datasets are openly available at https://doi.org/10.17871/ForestAgeBGI.2021 (Besnard et al., 2021).

Simon Besnard↗

Analysis of small-scale microwave background radiation anisotropy in the presence of foreground contamination

Many of the current round of experiments searching for anisotropies in the microwave background radiation (MBR) are confronting the problem of how to disentangle the cosmic signal from contamination due to Galactic and intergalactic foreground sources. Here we show how commonly used likelihood function techniques can be generalized to account for foreground. Specifically we set some restrictions on the spectrum of foreground contamination but allow the amplitude to vary arbitrarily. The likelihood function thus generalized gives reasonable limits on the MBR anisotropy which, in some cases, are not much less restrictive than what one would get from more detailed modeling of the foreground. Furthermore, the likelihood function is exactly the same as one would obtain by simply projecting out foreground contamination and looking at the reduced data set. We apply this generalized analysis to the recent medium-angle data sets of ACME-HEMT (Gaier et al. 1992; Schuster et al. 1993) and MAX (Meinhold et al. 1993; Gunderson et al. 1993). The resulting analysis constrains the one free parameter in the standard cold dark matter theory to be Q(sub rms-ps) = 18(sub -5 sup +8) microKelvin. This best fit value, although in striking agreement with the normalization from Cosmic Background Explorer (COBE), is not a very good fit, with an overall chi-squared/degrees of freedom = 208/168. We also argue against three commonly used methods of dealing with foreground: (1) ignoring it completely; (2) subtracting off a best-fit foreground and treating the residuals as if uncontaminated; and (3) culling data which appears to be contaminated by foreground.

Dodelson, Scott↗

Earth Radiation Imbalance from a Constellation of 66 Iridium Satellites: Climate Science Aspects

The "global warming hiatus" since the 1998 El Nino, highlighted by Meehl et al., and the resulting "missing energy" problem highlighted by Trenberth et al., has opened the door to a more fundamental view of climate change than mere surface air temperature. That new view is based on two variables which are strongly correlated: the rate of change of ocean heat content d(OHC)/dt; and Earth Radiation Imbalance (ERI) at the top of the atmosphere, whose guesstimated range is 0.4 to 0.9 Watts per square meters (this imbalance being mainly due to increasing CO2). The Argo float array is making better and better measurements of OHC. But existing satellite systems cannot measure ERI to even one significant digit. So, climate model predictions of ERI are used in place of real measurements of it, and the satellite data are tuned to the climate model predictions. Some oceanographers say "just depend on Argo for understanding the global warming hiatus and the missing energy", but we don't think this is a good idea because d(OHC)/dt and ERI have different time scales and are never perfectly correlated. We think the ERB community needs to step up to measuring ERI correctly, just as oceanographers have deployed Argo to measure OHC correctly. This talk will overview a proposed constellation of 66 Earth radiation budget instruments, hosted on Iridium satellites, that will actually be able to measure ERI to at least one significant digit, thus enabling a crucial test of climate models. This constellation will also be able to provide ERI at two-hourly time scales and 500-km spatial scales without extrapolations from uncalibrated narrowband geostationary instruments, using the highly successful methods of GRACE to obtain spatial resolution. This high time resolution would make ERI a synoptic variable like temperature, and allow studies of ERI's response to fast-evolving phenomena like dust storms and hurricanes and even brief excursions of Total Solar Irradiance. Time permitting, we will also discuss the emerging view of clear vs. cloudy and its implications for the traditional ERB approach.

Wiscombe, W.↗

Validation of High-Resolution CFD Method for Slosh Damping Extraction of Baffled Tanks

Determination of slosh damping is a very challenging task as there is no analytical solution. The damping physics involve the vorticity dissipation which requires the full solution of the nonlinear Navier-Stokes equations. As a result, previous investigations and knowledge were mainly carried out by extensive experimental studies. A Volume-Of-Fluid (VOF) based CFD program developed at NASA MSFC was applied to extract slosh damping in a baffled tank from the first principle. First, experimental data using water with subscale smooth wall tank were used as the baseline validation. CFD simulation was demonstrated to be capable of accurately predicting natural frequency and very low damping value from the smooth wall tank at different fill levels. The damping due to a ring baffle at different liquid fill levels from barrel section and into the upper dome was then investigated to understand the slosh damping physics due to the presence of a ring baffle. Based on this study, the Root-Mean-Square error of our CFD simulation in estimating slosh damping was less than 4.8%, and the maximum error was less than 8.5%. Scalability of subscale baffled tank test using water was investigated using the validated CFD tool, and it was found that unlike the smooth wall case, slosh damping with baffle is almost independent of the working fluid and it is reasonable to apply water test data to the full scale LOX tank when the damping from baffle is dominant. On the other hand, for the smooth wall, the damping value must be scaled according to the Reynolds number. Comparison of experimental data, CFD, with the classical and modified Miles equations for upper dome was made, and the limitations of these semi-empirical equations were identified.

Yang, H. Q.↗

Lidar Altimeter Measurements of Canopy Structure: Methods and Validation for Closed Canopy, Broadleaf Forests

Lidar altimeter observations of vegetated landscapes provide a time-resolved measure of laser pulse backscatter energy from canopy surfaces and the underlying ground. Airborne lidar altimeter data was acquired using the Scanning Lidar Imager of Canopies by Echo Recovery (SLICER) for a successional sequence of four, closed-canopy, deciduous forest stands in eastern Maryland. The four stands were selected so as to include a range of canopy structures of importance to forest ecosystem function, including variation in the height and roughness of the outer-most canopy surface and the vertical organization of canopy stories and gaps. The character of the SLICER backscatter signal is described and a method is developed that accounts for occlusion of the laser energy by canopy surfaces, transforming the backscatter signal to a canopy height profile (CHP) that quantitatively represents the relative vertical distribution of canopy surface area. The transformation applies an increased weighting to the backscatter amplitude as a function of closure through the canopy and assumes a horizontally random distribution of the canopy components. SLICER CHPs, averaged over areas of overlap where lidar ground tracks intersect, are shown to be highly reproducible. CHP transects across the four stands reveal spatial variations in vegetation, at the scale of the individual 10 m diameter laser footprints, within and between stands. Averaged SLICER CHPs are compared to analogous height profile results derived from ground-based sightings to plant intercepts measured on plots within the four stands. Tbe plots were located on the segments of the lidar ground tracks from which averaged SLICER CHPs were derived, and the ground observations were acquired within two weeks of the SLICER data acquisition to minimize temporal change. The differences in canopy structure between the four stands is similarly described by the SLICER and ground-based CHP results, however a Chi-square test of similarity documents differences that are statistically significant. The differences are discussed in terms of measurement properties that define the smoothness of the resulting CHPs and Lidar Altimeter Measurements of Canopy Structure - Harding et al. canopy properties that may vertically bias the CHP representations of canopy structure. The statistical differences are most likely due to the more noisy character of the ground-based CHPs, especially high in the canopy where ground-based sightings are rare resulting in an underestimate of canopy surface area and height, and to departures from the assumption of horizontal randomness which bias the CHPs toward the observer (upward for SLICER and downward for ground-based CHPs). The results demonstrate that the SLICER observations reliably provide a measure of canopy structure that reveals ecologically interesting structural variations such as those characterizing a successional sequence of closed-canopy, broadleaf forest stands.

Harding, D. J.↗

Validation Experiments for Turbulent Separation over an Axisymmetric Body of Revolution

Historically, the flow physics involved with most turbulent separated flows have presented fundamental challenges to validating numerical approaches. As recognized by the CFD Vision 2030 study commissioned by NASA, validation of Reynolds-averaged Navier-Stokes (RANS) models and other scale-resolving methods for turbulent separated flow requires data from advanced, high-fidelity experiments designed specifically for CFD implementation. In accordance with this effort, a new test platform, referred to as the NASA Axisymmetric Afterbody, was designed to obtain detailed measurements of the flow field undergoing a smooth, adverse pressure gradient induced separation for a fixed Reynolds number, Re = 180,000. The parametric body offers a range of flow states progressing from fully attached, to incipient separation, and finally to small-scale separated flow based on variable afterbody geometries. In an initial effort to evaluate RANS turbulence model capabilities, the present configurations of the axisymmetric model host a mild adverse pressure gradient over the contoured boattail section, inducing incipient turbulent separation, as well as a slightly larger adverse pressure gradient, inducing a small-scale region of turbulent separation. Experiments include steady pressure measurements and 2-D PIV to provide the preliminary dataset for simulation studies, which examine the effect of variable grid domains and RANS turbulence models. This is done in an effort to understand and evaluate the critical variability between solutions for the present model configurations. Results indicate potential discrepancies may be due to the effect of the square tunnel test section walls, relatively large blockage ratio, and slight variability in reference parameters. Ongoing work will focus on higher fidelity experimental campaigns to obtain surface flow visualizations and Stereo-Particle Image Velocimetry (SPIV) to deliver higher spatial resolution of the three-dimensional flow field to aid turbulence modelers.

Validation↗