Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Average Accuracy”

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 1,783 records · Page 99
All sourcesNASA NTRSDOE OSTI
NASA NTRS2011

Analysis of Oxygen Spectral Lines in the 1.27 Micron Band for the ASCENDS Mission

The ASCENDS mission requires simultaneous laser remote sensing of CO2 and O2 in order to convert CO2 column number densities to average column CO2 mixing ratios (XCO2). As such, the CO2 column number density and the O2 column number density will be utilized to derive the average XCO2 column. NASA Langley Research Center, working with its partners, is developing O2 lidar technology in the 1.26-1.27- m band for surface pressure measurements. The O2 model optical depth calculation is very sensitive to knowledge of the transmitted wavelengths and to the choice of Voigt input parameters. Modeling using the HITRAN database is being carried out to establish the evolution of candidate O2 absorption lines as a function of atmospheric parameters such as altitude, temperature, and pressure. Preliminary results indicate limitations of the Voigt profile and the need to utilize more advanced models which take into account line mixing, line narrowing, and speed dependence. In this paper, we evaluate alternative lineshape models to establish the optimum lineshapes which better account for the variability of individual O2 absorption lines at various atmospheric conditions. The combination of our modeling efforts with accurate laboratory measurements is anticipated to aid in achieving the desired CO2 mixing ratio measurement accuracy requirement of for the ASCENDS mission.

Pliutau, Denis↗
NASA NTRS2016

Considering Combined or Separated Roughness and Vegetation Effects in Soil Moisture Retrievals

For more than six years, the Soil Moisture and Ocean Salinity (SMOS) mission has provided multi angular and full-polarization brightness temperature (TB) measurements at L-band. Geophysical products such as soil moisture (SM) and vegetation optical depth at nadir (tau(sub nad)) are retrieved by an operational algorithm using TB observations at different angles of incidence and polarizations. However, the quality of the retrievals depends on several surface effects, such as vegetation, soil roughness and texture, etc. In the microwave forward emission model used in the retrievals (L-band Microwave Emission Model, L-MEB),soil roughness is modeled with a semi-empirical equation using four main parameters (Q(sub r), H(sub r), N(sub rp), with p = H or V polarizations). At present, these parameters are calibrated with data provided by airborne studies and in situ measurements made at a local scale that is not necessarily representative of the large SMOS footprints (43 km on average) at global scale. In this study, we evaluate the impact of the calibrated values of N(sub rp) and H(sub r) on the SM and tau(sub nad) retrievals based on SMOS TB measurements (SMOS Level 3 product) over the Soil Climate Analysis Network (SCAN) network located in North America over five years (2011-2015). In this study, Qr was set equal to zero and we assumed that N(sub rH)= N(sub rV). The retrievals were performed by varying N(sub rp) from −1 to 2 by steps of 1 and H(sub r) from 0 to 0.6 by steps of 0.1. At satellite scale, the results show that combining vegetation and roughness effects in a single parameter provides the best results in terms of soil moisture retrievals, as evaluated against the in situ SM data. Even though our retrieval approach was very simplified, as we did not account for pixel heterogeneity, the accuracy we obtained in the SM retrievals was almost systematically better than those of the Level 3 product. Improved results were also obtained in terms of optical depth retrievals. These new results may have key consequences in terms of calibration of roughness effects within the algorithms of the SMOS (ESA) and the SMAP (NASA) space missions.

full-polarization↗
NASA NTRS1993

Validation of Special Sensor Microwave Imager monthly-mean wind speed from July 1987 to December 1989

Since July 1987, the Remote Sensing Systems (Santa Rosa, California) routinely computes from first principles the wind speed from Special Sensor Microwave Imager (SSM/I) measurements of the intensity of microwave radiation emitted at the ocean surface. The accuracy of monthly-mean SSM/I wind speeds is determined by comparisons with moored-buoy wind measurements, which were recorded by four different institutions at 44 sites in the Atlantic and Pacific Oceans during July 1987 to December 1989. All results for 1988 were virtually identical with 1989. The range of monthly mean moored-buoy wind speeds was 2-10 m/s. The rms difference of 697 monthly-mean matchups of the composite 1988 and 1989 data set was 1.2 m/s. The rms differences were smaller in the equatorial zone and higher in middle latitudes, in accord with the monthly standard deviation. At middle latitudes, the time series of rms differences displayed an annual cycle. In the equatorial zone, the agreement between SSM/I and in situ data was better in regions with a lesser amount of clouds, and vice versa. For SSM/I monthly standard deviations of 1-2, 2-3, and 3-4 m/s, the average absolute values of the monthly-mean difference between SSM/I and moored-buoy wind speeds were 0.6, 0.9, and 1.4 m/s, respectively.

Halpern, David↗
NASA NTRS2000

Prospective Evaluation of Thoracic Ultrasound in the Detection of Pneumothorax

Introduction: Pneumothorax (PTX) occurs commonly in trauma patients and is confirmed by examination and radiography. Thoracic ultrasound (VIS) has been suggested as an alternative method for rapidly diagnosing PTX when X-ray is unavailable as in rural, military, or space flight settings; however, its accuracy and specificity are not known. Methods: We evaluated the accuracy of thoracic U/S detection of PTX compared to radiography in stable, emergency patients with a high suspicion of PTX at a Level-l trauma center over a 6-month period. Following University and NASA Institutional Review Board approval, informed consent was obtained from patients with penetrating or blunt chest trauma, or with a history consistent with PTX. Whenever possible, the presence or absence of the " lung sliding" sign or the "comet tail" artifact were determined by U/S in both hemithoraces by residents instructed in thoracic U/S before standard radiologic verification of PTX. Results were recorded on data sheets for comparison to standard radiography. Results: Thoracic VIS had a 94% sensitivity; two PTX could not be reliably diagnosed due to subcutaneous air; the true negative rate was 100%. In one patient, the VIS exam was positive while X ray did not confirm PTX; a follow-up film 1 hour later demonstrated a small PTX. The average time for bilateral thoracic VIS examination was 2 to 3 minutes. Conclusions: Thoracic ultrasound reliably diagnoses pneumothorax. Presence of the "lung sliding" sign conclusively excludes pneumothorax. Expansion of the FAST examination to include the thorax should be investigated.

Schwarz, K. W.↗
NASA NTRSDate not supplied

Transonic Lift and Drag Predictions using Wall Modelled Large Eddy Simulations

Wall-modelled Large Eddy Simulations (WMLES) of the NASA Common Research Model (CRM) at transonic conditions and various angles of attacks leading up to and including shock-induced flow separation are performed using the LAVA computational framework. The simulations are shown to accurately predict the lift curve slope and the onset of separation characterized by the break in the pitching moment. Furthermore, careful assessment of skin-friction drag at cruise condition and its subsequent decrease with increasing angles of attack are shown to be in agreement with viscous sublayer resolving Reynolds Averaged Navier Stokes(RANS) simulations. The small differences between WMLES and RANS appear to be of the same order as differences seen between two RANS models at the cruise-point. Some sensitivity is observed to the coefficient used in the subgrid scale model, although this can be reconciled by noting the low chord- and shock-incidence Reynolds numbers (Rex≈106) seen in the outboard regions of the wing, along with the uncertainties associated with tripping and the numerical transition that occurs near the leading edge. Among the primary drawbacks seen in the WMLES predictions is the subdued increase in wave-drag with increasing angle-of-attack and the shock intensity when compared with experimental data. Preliminary assessment suggests that additional span- and stream-wise grid resolution is likely needed in mid- and out-board portions of the wing to better resolve the shock-induced separated flow and to further investigate prediction accuracy of unsteady temporal characteristics of the problem.

ARMD↗
NASA NTRS1992

Phase noise in pulsed Doppler lidar and limitations on achievable single-shot velocity accuracy

The smaller sampling volumes afforded by Doppler lidars compared to radars allows for spatial resolutions at and below some sheer and turbulence wind structure scale sizes. This has brought new emphasis on achieving the optimum product of wind velocity and range resolutions. Several recent studies have considered the effects of amplitude noise, reduction algorithms, and possible hardware related signal artifacts on obtainable velocity accuracy. We discuss here the limitation on this accuracy resulting from the incoherent nature and finite temporal extent of backscatter from aerosols. For a lidar return from a hard (or slab) target, the phase of the intermediate frequency (IF) signal is random and the total return energy fluctuates from shot to shot due to speckle; however, the offset from the transmitted frequency is determinable with an accuracy subject only to instrumental effects and the signal to noise ratio (SNR), the noise being determined by the LO power in the shot noise limited regime. This is not the case for a return from a media extending over a range on the order of or greater than the spatial extent of the transmitted pulse, such as from atmospheric aerosols. In this case, the phase of the IF signal will exhibit a temporal random walk like behavior. It will be uncorrelated over times greater than the pulse duration as the transmitted pulse samples non-overlapping volumes of scattering centers. Frequency analysis of the IF signal in a window similar to the transmitted pulse envelope will therefore show shot-to-shot frequency deviations on the order of the inverse pulse duration reflecting the random phase rate variations. Like speckle, these deviations arise from the incoherent nature of the scattering process and diminish if the IF signal is averaged over times greater than a single range resolution cell (here the pulse duration). Apart from limiting the high SNR performance of a Doppler lidar, this shot-to-shot variance in velocity estimates has a practical impact on lidar design parameters. In high SNR operation, for example, a lidar's efficiency in obtaining mean wind measurements is determined by its repetition rate and not pulse energy or average power. In addition, this variance puts a practical limit on the shot-to-shot hard target performance required of a lidar.

Mcnicholl, P.↗
NASA NTRS1997

Snow Parameters Derived from Microwave Measurements During the BOREAS Winter Field Campaign

Passive microwave data have been used to infer the snow-covered area and snow water equivalent (SWE) over forested areas, but the accuracy of these retrieved snow parameters cannot be easily validated for heterogeneous vegetated regions. The Boreal Ecosystem-Atmosphere Study Winter Field Campaign provided the opportunity to study the effect of boreal forests on snow parameter retrieval in detail. Microwave radiometers (18, 37, and 92 GHz) were flown on board the Canadian National Aeronautical Establishment's Twin Otter. Flight lines covered both the southern study area near Prince Albert and the northern study area near Thompson, Canada. During the 1994 winter campaign, extensive ground-based snow cover information, including depth, density, and grain size, was collected along most of the flight lines, jointly by U.S. and Canadian investigators. Satellite data collected by the special sensor microwave imager are also used for comparison. Preliminary results reconfirmed the relationship between microwave brightness temperature and SWE. However, the effect of forest cover observed by the aircraft sensors is different from that of the satellite observations. This is probably due to the difference in footprint averaging. There were also several flight lines flown over Candle Lake and Waskesiu Lake to assess lake ice signatures. Preliminary results show the thickness of the lake ice may be inferred from the airborne microwave observations. The microwave signature relationship between lake ice and snow matches the results from radiative transfer calculations.

Chang, A. T. C.↗
NASA NTRS2010

Improving Forecast Skill by Assimilation of AIRS Temperature Soundings

AIRS was launched on EOS Aqua on May 4, 2002, together with AMSU-A and HSB, to form a next generation polar orbiting infrared and microwave atmospheric sounding system. The primary products of AIRS/AMSU-A are twice daily global fields of atmospheric temperature-humidity profiles, ozone profiles, sea/land surface skin temperature, and cloud related parameters including OLR. The AIRS Version 5 retrieval algorithm, is now being used operationally at the Goddard DISC in the routine generation of geophysical parameters derived from AIRS/AMSU data. A major innovation in Version 5 is the ability to generate case-by-case level-by-level error estimates delta T(p) for retrieved quantities and the use of these error estimates for Quality Control. We conducted a number of data assimilation experiments using the NASA GEOS-5 Data Assimilation System as a step toward finding an optimum balance of spatial coverage and sounding accuracy with regard to improving forecast skill. The model was run at a horizontal resolution of 0.5 deg. latitude X 0.67 deg longitude with 72 vertical levels. These experiments were run during four different seasons, each using a different year. The AIRS temperature profiles were presented to the GEOS-5 analysis as rawinsonde profiles, and the profile error estimates delta (p) were used as the uncertainty for each measurement in the data assimilation process. We compared forecasts analyses generated from the analyses done by assimilation of AIRS temperature profiles with three different sets of thresholds; Standard, Medium, and Tight. Assimilation of Quality Controlled AIRS temperature profiles significantly improve 5-7 day forecast skill compared to that obtained without the benefit of AIRS data in all of the cases studied. In addition, assimilation of Quality Controlled AIRS temperature soundings performs better than assimilation of AIRS observed radiances. Based on the experiments shown, Tight Quality Control of AIRS temperature profile performs best on the average from the perspective of improving Global 7 day forecast skill.

Susskind, Joel↗
NASA NTRS2004

Multisensor Arrays for Greater Reliability and Accuracy

Arrays of multiple, nominally identical sensors with sensor-output-processing electronic hardware and software are being developed in order to obtain accuracy, reliability, and lifetime greater than those of single sensors. The conceptual basis of this development lies in the statistical behavior of multiple sensors and a multisensor-array (MSA) algorithm that exploits that behavior. In addition, advances in microelectromechanical systems (MEMS) and integrated circuits are exploited. A typical sensor unit according to this concept includes multiple MEMS sensors and sensor-readout circuitry fabricated together on a single chip and packaged compactly with a microprocessor that performs several functions, including execution of the MSA algorithm. In the MSA algorithm, the readings from all the sensors in an array at a given instant of time are compared and the reliability of each sensor is quantified. This comparison of readings and quantification of reliabilities involves the calculation of the ratio between every sensor reading and every other sensor reading, plus calculation of the sum of all such ratios. Then one output reading for the given instant of time is computed as a weighted average of the readings of all the sensors. In this computation, the weight for each sensor is the aforementioned value used to quantify its reliability. In an optional variant of the MSA algorithm that can be implemented easily, a running sum of the reliability value for each sensor at previous time steps as well as at the present time step is used as the weight of the sensor in calculating the weighted average at the present time step. In this variant, the weight of a sensor that continually fails gradually decreases, so that eventually, its influence over the output reading becomes minimal: In effect, the sensor system "learns" which sensors to trust and which not to trust. The MSA algorithm incorporates a criterion for deciding whether there remain enough sensor readings that approximate each other sufficiently closely to constitute a majority for the purpose of quantifying reliability. This criterion is, simply, that if there do not exist at least three sensors having weights greater than a prescribed minimum acceptable value, then the array as a whole is deemed to have failed.

Immer, Christopher↗
DOE OSTI2023

Measurement of the muon anomalous precession frequency at the Muon g − 2 Experiment at Fermilab

The anomalous magnetic moment of the muon, $a_\mu = \frac{g-2}{2}$ is the fractional deviation of the muon $g$-factor from the value of 2. It emerges as the cumulative effect of the virtual particles participating in the muon interaction with a magnetic field via quantum loop corrections. Its value encodes all the possible interactions between the virtual particles and, for this reason, represents an important test of the Standard Model (SM). In particular, any deviation from the SM theoretical evaluation could be due to new physics contributions. The new Muon $g-2$ (E989) Experiment at Fermilab is currently operating to repeat and improve the previous E821 experiment at Brookhaven National Laboratory (BNL), aiming to reduce the experimental error by a factor of 4 to the final accuracy of 140 parts per billion (ppb). On April 7th, 2021, the E989 collaboration published the first result based on the first year of data taking (Run-1), measuring $a_\mu = 0.001~165~920~40(54)$ with a precision of 460 ppb. The measured value is consistent with the BNL measurement and strengthens the long-standing tension with the data-driven SM prediction to a combined discrepancy of 4.2$\sigma$. On the theory side, however, new efforts involving lattice-QCD techniques are starting to question the current consensus on the theoretical prediction, demanding new improvements on both the experimental and theoretical sides. The E989 collaboration is now finalizing the analysis of Run-2 and Run-3 data and a new publication is expected in the first half of 2023 with a combined statistical uncertainty of 200 ppb. The anomalous magnetic moment $a_\mu$ is measured as the ratio between the muon spin anomalous precession frequency, $\omega_a$, and the average magnetic field experienced by the muons as they circulate in the storage ring. This thesis presents a precession frequency analysis of the Run-1 data and an evaluation of the related systematic uncertainties. A new positron reconstruction developed for the analysis of the subsequent data-taking periods, aiming to reduce some of the major systematic uncertainties of the $\omega_a$ measurement, is presented. The author's involvement in the production of the Run-2/5 data and in the precise calibration of the detectors is discussed. Finally, the complete Run-1 $a_\mu$ result is presented.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

A passive model for the evolution of subgrid-scale instabilities in turbulent flow regimes

In this work, a zero-dimensional model for the development of turbulent mixing layers driven by Rayleigh–Taylor (RT) or Richtmyer–Meshkov (RM) instabilities is presented, for application in Reynolds-Averaged Navier–Stokes (RANS) simulations. The model approximates the behavior of the turbulent mixing layer at a passively advected tracer particle that sits on the material interface of interest in the simulation, and is used to provide initial conditions for the RANS model once the mixing layer width has developed to a size that may be resolved on the computational mesh. The model relies on approximations to the Besnard–Harlow–Rauenzahn (BHR) RANS closure equations and is thus consistent the turbulent development of the mixing layer, unlike previous models of this type which have focused on the early-time laminar development of RT or RM instabilities. These instabilities will often grow through an initially laminar state before developing into turbulence, and criteria for switching from an existing laminar model to the proposed turbulent model are considered. The proposed model reproduces the mixing layer development of RT and RM instabilities seen in mixed width resolving BHR RANS simulations to a reasonable degree of accuracy, and relaxes the need for the RANS simulation to resolve the mixed width as soon as the mixing layer becomes turbulent, which may be impractical in cases where the mixing layer grows through a wide range of scales.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Proton isovector helicity PDF at NNLO and the twist-3 moment $\tilde{d}$ 2 from lattice QCD at physical quark masses

We present a lattice quantum chromodynamics calculation of the 𝑥-dependent isovector quark helicity parton distribution function (PDF) of the proton in the large momentum effective theory (LaMET) framework. Through operator product expansion (OPE) we also extract the $\tilde{d}$ 2 moment of the twist-3 PDF 𝑔 𝑇 ⁡(𝑥) for the first time in the $\overline{MS}$ scheme, which is proportional to the average color Lorentz force experienced by the quark in the proton. This calculation is performed on a lattice of spacing 𝑎 =0.076 fm at physical quark masses. The quasi-PDF matrix elements are measured in proton states boosted to momenta 𝑃 𝑧 ={0,0.25,1.02,1.53} GeV. We first extract the lowest few helicity PDF moments from the renormalization-group (RG) invariant ratios of the matrix elements with OPE. Combined with the matrix elements relevant for 𝑔 𝑇 ⁡(𝑥), we obtain $\tilde{d}$$^{u-d}_2$⁡(2 GeV) =0.0024⁢(46) at next-to-leading order in $\overline{MS}$. Then, the helicity quasi-PDF matrix elements are renormalized in the hybrid scheme with linear renormalon resummation and Fourier transformed to the 𝑥-space after an asymptotic extrapolation. The quasi-PDF is perturbatively matched to the $\overline{MS}$ PDF with RG and threshold resummations at next-to-leading power and next-to-next-to-leading logarithmic accuracies. After resummations, we determine the PDF in the region 𝑥 ∈[0.25,0.75]. The end-point regions are then parametrized, combined with the LaMET prediction at moderate 𝑥, and fitted to the short-distance matrix elements in coordinate space.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Influence of intermolecular interactions on the infrared complex indices of refraction for binary liquid mixtures

This paper investigates the accuracy of deriving the composite optical constants of binary mixtures from only the complex indices of refraction of the neat materials. These optical constants enable the reflectance spectra of the binary mixtures to be modeled for multiple scenarios (e.g., different substrates, thicknesses, volume ratios), which is important for contact and standoff chemical detection. Using volume fractions, each mixture's complex index of refraction was approximated via three different mixing rules. To explore the impact of intermolecular interactions, these predictions are tested by experimental measurements for two representative sets of binary mixtures: (1) tributyl phosphate combined with n-dodecane, a non-polar medium, to represent mixtures which primarily interact via dispersion forces and (2) tributyl phosphate and 1-butanol to represent mixtures with polar functional groups that can also interact via dipole–dipole interactions, including hydrogen bonding. The residuals and the root-mean-square error between the experimental and calculated index values are computed and demonstrate that for miscible liquids in which the average geometry of the cross-interactions can be considered isotropic (e.g., dispersion), the refractive indices of the mixtures can be modeled using composite n and k values derived from volume fractions of the neat liquids. Conversely, in spectral regions where the geometry of the cross-interactions is more restricted and anisotropic (e.g., hydrogen bonding), the calculated n and k values vary from the measured values. The impact of these interactions on the reflectance spectra are then compared by modeling a thin film of the binary mixtures on an aluminum substrate using both the measured and the mathematically computed indices of refraction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2005

Development of a Coherent Lidar for Aiding Precision Soft Landing on Planetary Bodies

Coherent lidar can play a critical role in future planetary exploration missions by providing key guidance, navigation, and control (GNC) data necessary for navigating planetary landers to the pre-selected site and achieving autonomous safe soft-landing. Although the landing accuracy has steadily improved over time to approximately 35 km for the recent Mars Exploration Rovers due to better approach navigation, a drastically different guidance, navigation and control concept is required to meet future mission requirements. For example, future rovers will require better than 6 km landing accuracy for Mars and better than 1 km for the Moon plus maneuvering capability to avoid hazardous terrain features. For this purpose, an all-fiber coherent lidar is being developed to address the call for advancement of entry, descent, and landing technologies. This lidar will be capable of providing precision range to the ground and approach velocity data, and in the case of landing on Mars, it will also measure the atmospheric wind and density. The lidar obtains high resolution range information from a frequency modulated-continuous wave (FM-CW) laser beam whose instantaneous frequency varies linearly with time, and the ground vector velocity is directly extracted from the Doppler frequency shift. Utilizing the high concentration of aerosols in the Mars atmosphere (approx. two order of magnitude higher than the Earth), the lidar can measure wind velocity with a few watts of optical power. Operating in 1.57 micron wavelength regime, the lidar can use the differential absorption (DIAL) technique to measure the average CO2 concentration along the laser beam using, that is directly proportional to the Martian atmospheric density. Employing fiber optics components allows for the lidar multi-functional operation while facilitating a highly efficient, compact and reliable design suitable for integration into a spacecraft with limited mass, size, and power resources.

Amzajerdian, Farzin↗
NASA NTRS1986

A variational assimilation method for satellite and conventional data: Development of basic model for diagnosis of cyclone systems

A three-dimensional diagnostic model for the assimilation of satellite and conventional meteorological data is developed with the variational method of undetermined multipliers. Gridded fields of data from different type, quality, location, and measurement source are weighted according to measurement accuracy and merged using least squares criteria so that the two nonlinear horizontal momentum equations, the hydrostatic equation, and an integrated continuity equation are satisfied. The model is used to compare multivariate variational objective analyses with and without satellite data with initial analyses and the observations through criteria that were determined by the dynamical constraints, the observations, and pattern recognition. It is also shown that the diagnoses of local tendencies of the horizontal velocity components are in good comparison with the observed patterns and tendencies calculated with unadjusted data. In addition, it is found that the day-night difference in TOVS biases are statistically different (95% confidence) at most levels. Also developed is a hybrid nonlinear sigma vertical coordinate that eliminates hydrostatic truncation error in the middle and upper troposphere and reduces truncation error in the lower troposphere. Finally, it is found that the technique used to grid the initial data causes boundary effects to intrude into the interior of the analysis a distance equal to the average separation between observations.

Achtemeier, G. L.↗
NASA NTRS2013

Influence of Computational Drop Representation in LES of a Droplet-Laden Mixing Layer

Multiphase turbulent flows are encountered in many practical applications including turbine engines or natural phenomena involving particle dispersion. Numerical computations of multiphase turbulent flows are important because they provide a cheaper alternative to performing experiments during an engine design process or because they can provide predictions of pollutant dispersion, etc. Two-phase flows contain millions and sometimes billions of particles. For flows with volumetrically dilute particle loading, the most accurate method of numerically simulating the flow is based on direct numerical simulation (DNS) of the governing equations in which all scales of the flow including the small scales that are responsible for the overwhelming amount of dissipation are resolved. DNS, however, requires high computational cost and cannot be used in engineering design applications where iterations among several design conditions are necessary. Because of high computational cost, numerical simulations of such flows cannot track all these drops. The objective of this work is to quantify the influence of the number of computational drops and grid spacing on the accuracy of predicted flow statistics, and to possibly identify the minimum number, or, if not possible, the optimal number of computational drops that provide minimal error in flow prediction. For this purpose, several Large Eddy Simulation (LES) of a mixing layer with evaporating drops have been performed by using coarse, medium, and fine grid spacings and computational drops, rather than physical drops. To define computational drops, an integer NR is introduced that represents the ratio of the number of existing physical drops to the desired number of computational drops; for example, if NR=8, this means that a computational drop represents 8 physical drops in the flow field. The desired number of computational drops is determined by the available computational resources; the larger NR is, the less computationally intensive is the simulation. A set of first order and second order flow statistics, and of drop statistics are extracted from LES predictions and are compared to results obtained by filtering a DNS database. First order statistics such as Favre averaged stream-wise velocity, Favre averaged vapor mass fraction, and the drop stream-wise velocity, are predicted accurately independent of the number of computational drops and grid spacing. Second order flow statistics depend both on the number of computational drops and on grid spacing. The scalar variance and turbulent vapor flux are predicted accurately by the fine mesh LES only when NR is less than 32, and by the coarse mesh LES reasonably accurately for all NR values. This is attributed to the fact that when the grid spacing is coarsened, the number of drops in a computational cell must not be significantly lower than that in the DNS.

Bellan, Josette↗
NASA NTRS2022

Assessment of ICESat-2’s Horizontal Accuracy Using Precisely Surveyed Terrains in McMurdo Dry Valleys, Antarctica

This article presents an assessment of the horizontal accuracy and precision of the laser altimetry observations collected by NASA’s Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2) mission. We selected the terrain-matching method to determine the position of laser altimeter profiles within a precisely known surface, represented by a digital elevation model (DEM). We took this classical approach a step further, approximated the DEM by planar surfaces, and calculated the optimal position of the laser profile by minimizing the square sum of the elevation differences between reference DEMs and ICESat-2 profiles. We found the highly accurate DEMs of the McMurdo Dry Valleys (DV), Antarctica, ideal for this research because of their stable landscape and rugged topography. We computed the 3-D shift parameters of 379 different laser altimeter profiles along two reference ground tracks collected within the first two years of the mission. Analyzing these results revealed a total geolocation error (mean +1 σ ) of 4.93 m for version 3 and 4.66 m for version 4 data. These numbers are the averages of the six beams, expressed as mean +1 σ and lie well within the mission requirement of 6.5 m.

calibration↗