Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “INTERCEPTION”

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,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
NASA NTRS2021

PolyFit: A C++ code for Polynomial Curve Fit with Calculation of Error Bars

In radiobiology, many dose-response results are modeled using the so-called linear-quadratic(LQ)model, which means that results are modeled as a function of dose Das 𝑅(𝐷)=𝛽0+𝛽1𝐷+𝛽2𝐷2. The coefficients 𝛽0, 𝛽1and 𝛽2are obtained from fitting a series of data points (𝑥𝑖,𝑦𝑖), which is usually done using a least-square method. The LQ and more generally the polynomial fit capability is implemented in many software that analyzes data. However, it is often convenient to do the fitting programmatically, especially when a large number of datasets should be analyzed. Furthermore, depending on the software used, some features may not be implemented. In this mini-review, I discuss the basis of polynomial fitting, including the calculation of errors on the coefficients and results, use of weighting and fixing the intercept value (the coefficient 𝛽0). A simple C++ code to perform the polynomial curve fitting is also provided. This code should be useful not only in radiobiology but in other fields of science as well.

Ianik Plante↗
NASA NTRSDate not supplied

Vapor Cooling of a Structural Skirt for a Large-Scale Hydrogen Tank

The demonstration of vapor cooling on a structural skirt was one of the main objectives of the Structural Heat Intercept, Insulation, and Vibration Evaluation Rig (SHIIVER) testing. SHIIVER consisted of a 4-meter diameter stainless steel tank with an aluminum forward skirt. The forward skirt was cooled by directing effluent vapor from the tank through two flow channels that each went 360 degrees around the interior of the skirt in a helical pattern. Flow rates, temperatures, and pressures in the system were measured allowing for the calculation of heat load removed via cooling stream and heat load reduction into the tank. Testing occurred at fill levels between 25% and 90% full using both liquid hydrogen and liquid nitrogen. Boil-off rate was varied independent of the skirt performance by adding multilayer insulation over the spray-on-foam insulation on the tank domes, while leaving the barrel insulated only with spray-on-foam. The results from the testing, which demonstrate vapor cooling reversed heat flow from the skirt into the tank and reduces total propellant heat load by as much as 19%, are analyzed and discussed.

Cryogenic Fluid Management↗
NASA NTRSDate not supplied

Results of Use of Heat Flux Sensors on Liquid Hydrogen Tanks

Heat flux sensors were used to characterize the performance of insulation on the Structural Heat Intercept, Insulation, and Vibration Evaluation Rig (SHIIVER), a large-scale test article designed to simulate upper stage cryogenic propellant tank thermal performance in simulated space environments. Usually, the insulation heat loads are derived from calculations removing all other heat sources and attributing the residual heat load to the insulation system. Testing for SHIIVER included the tank being just insulated with spray-on-foam insulation as well as covering the domes with multilayer insulation while leaving the barrel section insulated with spray-on-foam. Heat flux sensors were located at multiple locations on both domes as well as on the barrel section of the tank. Results from the SHIIVER testing using the heat flux sensors are compared to other calculated heat inputs for both liquid nitrogen and liquid hydrogen testing as a function of tank fill level. Further investigation into the transient nature of the SHIIVER testing including the heat flux sensors provided insight into heat flow patterns that may not have been otherwise seen using temperature sensors and calculated insulation heat loads. While the demonstrated uncertainties in the absolute values in the heat flux sensors are high, the values and the trends match well with other calculation methods. The results of SHIIVER allow for the use of heat flux sensors for measurements of insulation performance and dynamic system thermal response for future applications.

Cryogenic Fluid Management↗
NASA NTRSDate not supplied

Results of Use of Heat Flux Sensors on Liquid Hydrogen Tanks

ABSTRACT: Heat flux sensors were used to characterize the performance of insulation on the Structural Heat Intercept, Insulation, and Vibration Evaluation Rig (SHIIVER), a large-scale test article designed to simulate upper stage cryogenic propellant tank thermal performance in simulated space environments. Usually, the insulation heat loads are derived from calculations removing all other heat sources and attributing the residual heat load to the insulation system. Testing for SHIIVER included the tank being just insulated with spray-on-foam insulation as well as covering the domes with multilayer insulation while leaving the barrel section insulated with spray-on-foam. Heat flux sensors were located at multiple locations on both domes as well as on the barrel section of the tank. Results from the SHIIVER testing using the heat flux sensors are compared to other calculated heat inputs for both liquid nitrogen and liquid hydrogen testing as a function of tank fill level. Further investigation into the transient nature of the SHIIVER testing including the heat flux sensors provided insight into heat flow patterns that may not have been otherwise seen using temperature sensors and calculated insulation heat loads. While the demonstrated uncertainties in the absolute values in the heat flux sensors are high, the values and the trends match well with other calculation methods. The results of SHIIVER allow for the use of heat flux sensors for measurements of insulation performance and dynamic system thermal response for future applications.

Cryogenic Fluid Management↗
NASA NTRSDate not supplied

Evolution of Drop Size Distributions during Hurricane Harvey and Their Effects on Rain Intensity

Hurricane Harvey hit the Texas Gulf Coast as a major hurricane on August 25, 2017, before exiting the state as a tropical storm on August 29, 2017. Left in its wake was historic flooding over southeast Texas with some locations measuring more than 60 inches of rain. The Houston NEXRAD radar (KHGX) maintained operations for the entirety of the event and is utilized in this study to better understand the evolution of the drop size distribution (DSD) in the lower (liquid) portions of the storm and how changes in the DSD were tied to observed and measured changes in rain intensities. From the KHGX dual-polarimetric radar data several DSD parameters were retrieved, including three different rain rate retrievals, as well as the mass-weighted mean-drop diameter Dm, and the normalized intercept parameter NW of the gamma distribution. Two of the rain retrievals are popular hybrid approaches, where reflectivity, differential and specific differential phase are used. The other is an attenuation-based retrieval. In a previous study, it was found that both hybrid retrievals significantly underestimated the rainfall during the heaviest rain periods, while the attenuation-based retrieval, with proper adjustment, agreed quite well with a large network of gauges from the Harris County Flood Warning System. Using Dm and NW as proxies for drop size and concentration, respectively, we will show that the heaviest rain intensities were associated in areas of the storm where there was a large concentration of large drops present. It is envisaged that results of this study can be used to develop a near-real-time tuning of the attenuation-based retrieval, thus providing timely, and possibly life-saving warnings during inundating rainfall.

David B. Wolff↗
NASA NTRS2021

Measurements of Rainfall Rate, Drop Size Distribution, and Variability at Middle and Higher Latitudes: Application to the Combined DPR-GMI Algorithm

The Global Precipitation Measurement mission is a major U.S.–Japan joint mission to understand the physics of the Earth’s global precipitation as a key component of its weather, climate, and hydrological systems. The core satellite carries a dual-precipitation radar and an advanced microwave imager which provide measurements to retrieve the drop size distribution (DSD) and rain rates using a Combined Radar-Radiometer Algorithm (CORRA). Our objective is to validate key assumptions and parameterizations in CORRA and enable improved estimation of precipitation products, especially in the middle-to-higher latitudes in both hemispheres. The DSD parameters and statistical relationships between DSD parameters and radar measurements are a central part of the rainfall retrieval algorithm, which is complicated by regimes where DSD measurements are abysmally sparse (over the open ocean). In view of this, we have assembled optical disdrometer datasets gathered by research vessels, ground stations, and aircrafts to simulate radar observables and validate the scattering lookup tables used in CORRA. The joint use of all DSD datasets spans a large range of drop concentrations and characteristic drop diameters. The scaling normalization of DSDs defines an intercept parameter N(W), which normalizes the concentrations, and a scaling diameter D(m), which compresses or stretches the diameter coordinate axis. A major finding of this study is that a single relationship between N(W) and D(m), on average, unifies all datasets included, from stratocumulus to heavier rainfall regimes. A comparison with the N(W)–D(m) relation used as a constraint in versions 6 and 7 of CORRA highlights the scope for improvement of rainfall retrievals for small drops (D(m) < 1 mm) and large drops (D(m) > 2 mm). The normalized specific attenuation–reflectivity relationships used in the combined algorithm are also found to match well the equivalent relationships derived using DSDs from the three datasets, suggesting that the currently assumed lookup tables are not a major source of uncertainty in the combined algorithm rainfall estimates.

Viswanathan Bringi↗
NASA NTRSDate not supplied

Spectral Properties of Dust on Asteroid (101955) Bennu with the OSIRIS-REx Visible and Infrared Spectrometer (OVIRS)

The Origins, Spectral Interpretation, Resource Identification, Security-Regolith Explorer (OSIRIS-REx) spacecraft successfully collected a sample of regolith from asteroid (101955) Bennu on October 20, 2020[1]. The Touch-and-Go (TAG) spacecraft maneuver used to collect the sample mobilized material, including very fine particles (dust),from up to ~1m in depth and ~10m in diameter from the point of contact with the asteroid [1,2].TAG created a dust plume that intercepted the spacecraft and accumulated on the instruments. As a result, all instruments show degradation in optical throughput, with the degree of degradation depending on their aperture size, orientation, and position on the spacecraft [1]. In this study, we analyzed OSIRIS-REx Visible and Infrared Spectrometer (OVIRS) from before and after TAG to assess the dust on the instrument with implications for Bennu’s dust composition and scattering properties, and future visible–near-infrared (VNIR)spectral observations with the OSIRIS-REx spacecraft. OVIRS is a point spectrometer that measures reflected light at VNIR wavelengths from 0.4 to 4.3 μm with a 4 mrad circular field of view(FOV) [3]. Previous observations with OVIRS revealed widespread hydrated minerals and carbon-bearing materials on Bennu [4,5]. A ~3 μm absorption feature is observed globally in OVIRS spectra, with a band position and depth that is consistent with Mg/Fe phyllosilicates and suggests a similar degree of aqueous alteration to that experienced by some carbonaceous chondrite meteorites[4,6]. A series of weak (<5%) VNIR spectral features, potentially associated with hydrated phyllosilicates and iron oxides, are also found in OVIRS spectra (Fig. 1), though not typically observed in meteorite spectra [7]. OSIRIS-REx Thermal Emission Spectrometer (OTES) spectra are consistent with a very thin accumulation (a few to ~ten microns) of fine particles (<~65–100 microns in size)across the surface of Bennu, and particularly on the roughest boulders, an interpretation that is supported by the dust plume seen during TAG [1,8]. An analysis of TAG contamination on the OTES optics revealed a 15% decrease in throughput, and OTES observations of space with and without that contamination allow a recreation of the thermal infrared spectrum [1]. The spectrum of the contamination lacks volume scattering features and exhibits a stronger Mg-OH absorption (~16.5 μm) than average Bennu, suggesting a Mg-rich phyllosilicate composition [1,6]. Here we present the first compositional analysis of Bennu dust at VNIR wavelengths by analyzing the contaminant on OVIRS.

H H Kaplan↗
NASA NTRSDate not supplied

A Combined Al-Mg/Pb-Pb Age of the Solar System

Astrophysical models of planet formation and protoplanetary disk evolution demand precise and accurate timing of the sequence of events in the solar nebula, relative to a time t=0, usually taken to be during the short epoch of CAI (Ca-rich, Al-rich inclusion) formation. Most CAIs formed withlive26Al (mean-life τ26= 1.034 Myr [1]), with an abundance 26Al/27Al ≈ (26Al/27Al)SS= 5.23 × 10-5[2]. We adopt this as the widespread level of 26Al in the solar nebula at t=0. Assuming spatial homogeneity of 26Al, an inclusion that had less 26Al, (26Al/27Al)0, formed a time Δt26= τ26ln[(26Al/27Al)SS/ (26Al/27Al)0] after t=0.These ages are typical precise to within ±0.1 Myr. Igneous bulk meteorites and inclusions can be relatively dated by the Al-Mg chronometer, but only ifΔt26<6 Myr. The Pb-Pb system is useful as a longer relative chronometer. It yields absolute ages tPb using 207Pb/206Pb, 206Pb/204Pb, and 238U/235U ratios measured indifferent portions of a sample, assuming certain half-lives [4]. These absolute ages are uncertain to within ±9 Myr due to uncertainties in the 235U half-life[3], but times of formation ΔtPb= tCAI–tPb relative to t=0, are more precise(±0.5Myr),iftCAI can be found. Here, tCAI means the Pb-Pb age that would be measured in CAIs using the half-lives the community typically uses, if they achieved isotopic closure at t=0. Unfortunately, direct Pb-Pb dating of CAIs has not definitively determined tCAI. Based on four CAIs with canonical (26Al/27Al)0,[5,6] found tPb= 4567.30 ± 0.16 Myr. No other CAI ages with measured 238U/235U have been reported in the refereed literature, but there are hints of other CAIs with ages tPb= 4568.0 ± 0.2 Myr [7] and tPb= 4568.3 ± 0.2 Myr [8].It is unclear whether anyof these igneous type B CAIs isotopically closed at t=0 or represents tCAI. Instead of measurements, we advocate finding tCAI by minimizing the discrepancies between the Al-Mg and Pb-Pb chronometers. Assuming Δt26=ΔtPb, we find the implied t’CAI= tPb+Δt26, then define t*CAIas the weighted mean of the t’CAI. t*CAIis the best guess for the Pb-Pb age of t=0; the assumption of homogeneity is justified if the t’CAI cluster within errors around t*CAI. This statistical approach is similar to, but improves on, that of[9]. We find t*CAI= 4568.73 ± 0.16 Myr. Below we discuss our methodology and the implications of this age for CAIs, 1.4 Myr older than the reported and typically used age 4567.30±0.16 Myr. Methods: We base our estimate of t*CAIon five achondrites for which published (26Al/27Al)0and Pb-Pb ages exist: the quenched angrites D’Orbigny, Sahara 99555 (SAH 99555), and Northwest Africa (NWA) 1670; the pseudo-eucrite Asuka 881394; and the inner disk achondrite. All are “NC” (non-carbonaceous) achondrites that likely cooled quickly enough that the Al-Mg and Pb-Pb systems achieved isotopic closure simultaneously. We also considered the “CC” (carbonaceous chondrite-like) achondrites NWA 2796 and NWA 6704, butdo not include them in our fit. Al-Mg and Pb-Pb seem not to have closed simultaneously, possibly because formation in the outer disk from volatile-rich composition led to slower cooling. Of the 8 chondrules from NWA 5697 measured by [20,21], we also consider the 4 for which 238U/235U was measured: 2-C1, 5-C2, 3-C5, 11-C1.Depending on their post-formation thermal histories, the Al-Mg and Pb-Pb systems in chondrules may or may not have closed simultaneously. Table 1: (26Al/27Al)0, Pb-Pb ages of selected samplesSample(26Al/27Al)0/ 10-6RefPb-PbRefD’Orbigny3.98±0.15104563.43±0.19♮10-12SAH 995553.64±0.18104563.88±0.2712NWA16705.92±0.59104564.39±0.24*10Asuka 88139413.1±0.5613-154564.98±0.1715NWA 73253.03±0.14164563.4±2.616NWA 27963.94±0.16174562.89±0.5917NWA 67043.03±0.14184562.76±0.26192-C17.56±1.53204567.57±0.56*215-C27.04±1.51204567.54±0.52*213-C58.85±1.83204566.20±0.63*2111-C15.55±1.84204565.84±0.72*21*regression based on one subset of data points ♮weighted mean of two datasets Pb-Pb ages are proportional to the intercept of the line formed by linear regression of 207Pb/206Pb vs. 206Pb/204Pb data from various washes, leachates and residues of acid dissolution of a sample. Because contamination by terrestrial or primordial Pb is pervasive, some fractions must be excluded from regressions to ensure a fit with acceptable mean squares weighted deviation (MSWD). Usually points are excluded based on low [Pb], or low 206Pb/204Pb ratio(low radiogenic component), with single outliers identified [11,12,15,16,17]. In the starred examples (Table 1)and the case of 3 CAI Pb-Pb ages [5], up to half the points were excluded solely because did not fit a pre-determined line. This approach is vulnerable to confirmation bias and produces fits with low MSWD and too-low Pb-Pb age uncertainty. Regressing the same data points as [10], were produce the Pb-Pb age of NWA 1670 of 4564.39±0.24 Myr. But selecting other combinations of data points, other, equally valid, isochrons yield ages from 4563.77±0.21 Myr to 4564.64±0.23 Myr. Similar arguments apply to the Pb-Pb isochrons built by [21] for chondrules 2-C1 (we find 4567.33±0.44 to 4567.85±0.46 Myr), 5-C2 (4566.84±0.53 to 4567.70±0.44 Myr), 3-C5 (4565.84±0.54to 4567.04±0.54)and 11-C1 (4565.36±0.51 to 4565.74±0.45 Myr). Our adopted ages for these and NWA 1670 are listed in Table 2.Table 2. tCAI estimated from various components, using our regressions for the chondrules & NWA 1670.SampleΔt26(Myr)tPb(Myr)t’CAI(Myr)D’Orbigny5.05±0.044563.43±0.194568.48±0.19SAH 995555.14±0.054563.88±0.274569.02±0.27NWA16704.64±0.104564.21±0.634568.85±0.67Asuka 8813943.81±0.044564.98±0.174568.79±0.17NWA 73255.33±0.054563.4±2.64568.7±2.6NWA 27965.06±0.044562.89±0.594567.95±0.59NWA 67045.29±0.134562.76±0.264568.05±0.292-C12.00±0.214567.59±0.704569.59±0.725-C22.07±0.224567.23±0.914569.30±0.933-C51.84±0.214566.44±1.124568.28±1.1411-C12.32±0.344565.52±0.664567.84±0.73achondrite4568.72±0.16chondrules4568.76±0.58combined4568.73±0.16A weighted average of the five NC achondrites(or just D’Orbigny, SAH 99555 and Asuka 881394)yields t*CAI= 4568.72 ± 0.16Myr. All are consistent with this value to within 1.8σ, and MSWD=1.5. Including the 4 U-corrected chondrules, t*CAI= 4568.73± 0.16Myrwith MSWD=1.66, which is statistically significant. All chondrules and NC achondrites are consistent with this to within 1.8σ, (Figure 1).Figure 1. Al-Mg formation times after t=0 vs. Pb-Pb ages. The five NC achondrites and four chondrules are consistent with a Pb-Pb age of t=0 of 4568.7 Myr. Discussion: The data from achondrites and chondrules are consistent with a single Pb-Pb age at t=0, justifying the assumption of 26Al homogeneity. The age, 4568.7 Myr, is ≈1.4 Myr older than the commonly accepted Pb-Pb age of CAIs that formed with canonical 26Al/27Al at t=0 [3]. Others have interpreted the discrepancy to signify 26Al heterogeneity in the CAI-forming region[5,21]. We suggest instead that CAIs were exposed to transient heating events that reset the Pb-Pb system without disturbing the Al-Mg system. Notably, chondrules typically experienced transient heating at these times in the nebula [22]. If so, direct measurements of CAIs will not yield as reliable a Pb-Pb age of t=0 as statistical approaches like this and that of [9].References:[1] Auer et al. 2009. [2] Jacobsen, B et al. 2008, EPSL 272, 353-364. [3] Tissot, Fet al. 2017, GCA 213, 593-617.[4] Villa, I et al. 2016, GCA 172, 387-392.[5] Amelin, Y et al. 2010, EPSL 300, 343-350.[6] Connelly, Jet al. 2012, Science 338, 651.[7] Bouvier, Aet al. 2011, LPICo 1639, 9054. [8] Bouvier, Aand Wadhwa, M2010, Nat Geosci 3, 637-641.[9] Nyquist, Let al. 2009, GCA 73, 5115-5136.[10] Schiller et al. 2015. [11] Wadhwa & Brennecka 2012. [12] Tissot et al. 2017. [13] Nyquist et al. 2003. [14] Wadhwa et al. 2009. [15] Wimpenny et al. 2019, GCA 244, 478-501.[16] Koefoed et al. 2016, GCA 183, 31-45. [17] Bouvier, A et al. 2011, GCA 75, 5310-5323. [18] Sanborn, Met al. 2019, GCA 245, 577-596. [19] Amelin, Yet al. 2019, GCA 245, 628-642.[20] Bollard, Jet al. 2017, Sci Adv 3 ,e1700407. [21] Bollard, Jet al. 2019, GCA 260, 62-83.[22] Villeneuve, J et al. 2009, Science 325, 985

S. J. Desch↗
NASA NTRS2021

TLSLEAF: Automatic Leaf Angle Estimates From Single-Scan Terrestrial Laser Scanning

Leaf angle distribution (LAD) in forest canopies affects estimates of leaf area, light interception, and global-scale photosynthesis, but is often simplified to a single theoretical value. Here, we present TLSLeAF (Terrestrial Laser Scanning Leaf Angle Function), an automated open-source method of deriving LADs from terrestrial laser scanning. TLSLeAF produces canopy-scale leaf angle and LADs by relying on gridded laser scanning data. The approach increases processing speed, improves angle estimates, and requires minimal user input. Key features are automation, leaf–wood classification, beta parameter output, and implementation in R to increase accessibility for the ecology community. TLSLeAF precisely estimates leaf angle with minimal distance effects on angular estimates while rapidly producing LADs on a consumer-grade machine. We challenge the popular spherical LAD assumption, showing sensitivity to ecosystem type in plant area index and foliage profile estimates that translate to c. 25% and c. 11% increases in canopy net photosynthesis (c. 25%) and solar-induced chlorophyll fluorescence (c. 11%). TLSLeAF can now be applied to the vast catalog of laser scanning data already available from ecosystems around the globe. The ease of use will enable widespread adoption of the method outside of remote-sensing experts, allowing greater accessibility for addressing ecological hypotheses and large-scale ecosystem modeling efforts.

3D↗
NASA NTRS2021

Liana optical traits increase tropical forest albedo and reduce ecosystem productivity

Lianas are a key growth form in tropical forests. Their lack of self-supporting tissues and their vertical position on top of the canopy make them strong competitors of resources. A few pioneer studies have shown that liana optical traits differ on average from those of colocated trees. Those trait discrepancies were hypothesized to be responsible for the competitive advantage of lianas over trees. Yet, in the absence of reliable modelling tools, it is impossible to unravel their impact on the forest energy balance, light competition, and on the liana success in Neotropical forests. To bridge this gap, we performed a meta-analysis of the literature to gather all published liana leaf optical spectra, as well as all canopy spectra measured over different levels of liana infestation. We then used a Bayesian data assimilation framework applied to two radiative transfer models (RTMs) covering the leaf and canopy scales to derive tropical tree and liana trait distributions, which finally informed a full dynamic vegetation model. According to the RTMs inversion, lianas grew thinner, more horizontal leaves with lower pigment concentrations. Those traits made the lianas very efficient at light interception and significantly modified the forest energy balance and its carbon cycle. While forest albedo increased by 14% in the shortwave, light availability was reduced in the understorey (−30% of the PAR radiation) and soil temperature decreased by 0.5°C. Those liana-specific traits were also responsible for a significant reduction of tree (−19%) and ecosystem (−7%) gross primary productivity (GPP) while lianas benefited from them (their GPP increased by +27%). This study provides a novel mechanistic explanation to the increase in liana abundance, new evidence of the impact of lianas on forest functioning, and paves the way for the evaluation of the large-scale impacts of lianas on forest biogeochemical cycles.

Félicien Meunier↗
NASA NTRS2021

Testing the Drop-Size Distribution-Based Separation of Stratiform and Convective Rain Using Radar and Disdrometer Data from a Mid-Latitude Coastal Region

Stratiform and convective rain are associated with different microphysical processes and generally produce drop-size distributions (DSDs) with different characteristics. Previous studies using data from (a) a tropical coastal location, (b) a mid-latitude continental location with semi-arid climate, and (c) a sub-tropical continental location, found that the two rain types could be separated in the NW–Dm space, where Dm is the mass-weighted mean diameter and NW is the normalized intercept parameter. In this paper, we investigate the same separation technique using data and observations from a mid-latitude coastal region. Three-minute DSDs from disdrometer measurements are used for the NW- versus Dm-based classification and are compared with simultaneous observations from an S-band polarimetric radar 38 km away from the disdrometer site. Specifically, RHI (range-height indicator) scans over the disdrometer were used for confirmation. Results show that there was no need to modify the separation criteria from previous studies. Three-minute DSDs from the same location were used as input to scattering calculations to derive retrieval equations for NW and Dm for the S-band radar using an improved technique and applied to the RHI scans to identify convective and stratiform rain regions. Two events are shown as illustrative examples.

Merhala Thurai↗
NASA NTRSDate not supplied

Zero-Trust Architecture for Autonomous Edge Computing

We are at the apex of an aviation revolution where autonomy will play a central role in enabling complex, multi-agent systems to communicate, interact, and collaborate on a myriad of applications spanning autonomous swarms to wild-fire management. Autonomy is not an absolute but rather a spectrum ranging from a system requiring significant human intervention to one requiring little to none [1]. For example, the extreme, in the case of an autonomous aircraft, is one that operates independently in the airspace interacting with all other elements (air traffic controllers, other pilots) as if it were a human pilot. Critical to this vision is an architecture that enables autonomous agents to interact with minimal latency. Edge computing is an emerging architecture where compute and storage is pushed to the ‘edge’ of the network in order to minimize the round-trip time from agent to resource thereby mitigating the latency associated with cloud-only based approaches. Additionally, services can generate massive amounts of data (e.g., video feeds), which may require analysis in near real-time. Moving this data to the cloud for further processing may not be feasible due to latency, bandwidth, and cost. Privacy, security, and reliability can also be improved by edge computing architectures. However, this geo-distributed and dynamic* architecture complicates the establishment of unambiguous network security boundaries and can lead to vulnerabilities including man in the middle attacks, replay attacks, physical security breaches of edge nodes, signal interception, etc. This motivates the need for zero-trust architectures [2–4] which de-emphasize the notion of static network perimeters and, as the name implies, do not instill any innate trust in any particular agent. It is required that all agents must be authorized and approved in every transaction. In this paper, we present a zero-trust architecture suitable for edge-computing applications that demand significant low-latency, security, privacy, and reliability.

zero trust↗
NASA NTRS2022

Comparison of Likelihood Methods for Generalized Linear Mixed Models with Application to Quiet Supersonic Flights 2018 Data

Repeated measurement will be a feature of the survey data collected during the Quesst missionX-59 community response tests (CRT). Since each participant will report his or her categorical level of annoyance in response to multiple events, the responses from any single individual may be correlated with one another. Several models within the class of generalized linear mixed models (GLMM) are pertinent to the analysis of correlated categorical outcomes; the random intercept logistic regression model is one example. Both Bayesian and frequentist methods for fitting these models are available, with frequentist methods relying on some form of approximation (of either an integral or the integrand) that appears in the marginal likelihood function. Given several anticipated similarities of the X-59 CRT data to data collected during a past risk reduction, Quiet Supersonic Flights 2018 (QSF18), this short note is intended to create awareness. It documents an instance in which a reported population average dose-response relationship derived from QSF18 single event data was distorted by the integral approximation applied in likelihood-based methods. We review some of the available literature on the topic, compare the outputs of several different computational approaches implemented in available statistical software, and present simple corrective actions that may be useful during the Quesst mission.

dose-response model↗
NASA NTRSDate not supplied

Where Lies the Quantumness Behind Detecting Electromagnetic Waves for Frequencies from Infrared and Up?

The synthesis of Newtonian concept of corpuscles during emission with Huygens’ concept of secondary-wavelets during propagation implies that all EM radiations from quantized atoms and molecules are released as discrete amount of energies. However, they propagate out as time-finite Maxwellian light pulses. Huygens also underscored that his secondary wavelets keep propagating as independent pulses in the absence of any interacting medium, or until intercepted by an interface with a medium or a detector. Then we use the Superposition Principle and the coherence theory to derive Einstein’s photoelectric equation by summing innumerable random time-finite pulses. This process driven approach should yield the characteristic statistical variations of photoelectron current pulses, as generated by photodetectors for different kinds of light sources. Lamb & Scully originally proposed this semiclassical approach without assuming that light actually consists of time finite pulses. The quantumness remains confined within the excitation and de-excitations processes in the material particles.

Electromagnetic Waves↗
NASA NTRS2022

Assessment of Geo-Kompsat-2A Atmospheric Motion Vector Data and Its Assimilation Impact in the GEOS Atmospheric Data Assimilation System

Korea’s second geostationary meteorological satellite, Geo-Kompsat-2A (Geostationary-Korean Multi-Purpose Satellite-2A, GK2A), was successfully launched on 4 December 2018. GK2Agenerates Atmospheric Motion Vectors (AMVs) every 10 min in the full disk area. This data hasbeen disseminated via Global Telecommunication System (GTS) since 25 October 2019. This articleevaluates the quality of GK2A AMVs in the Goddard Earth Observing System (GEOS) atmosphericdata assimilation system (ADAS). The data show slow wind speed biases at 200–300 hPa and 600–800hPa in the northern and southern hemispheres. These biases are caused by observation heightassignment errors near jet streams. The Equivalent Blackbody Temperature (EBBT) method of GK2Atends to assign clouds at higher altitude, which mainly causes slow wind speed biases, especially inthe lower atmosphere. The IR/WV intercept method of GK2A assigns clouds slightly lower in theatmospheric layers below the altitude of 400 hPa, which causes positive biases. Quality control (QC)criteria to select the most suitable GK2A AMV data for assimilation are presented based on thesequality assessments. A new QC criterion utilizing height errors within the GEOS ADAS is introducedto exclude data with slow wind speed biases and large errors. GEOS forecast accuracy is slightlyimproved after assimilating GK2A AMVs along with other conventional, radiance, and satellitewinds which include AMVs made by the Himawari-8 satellite in nearly the same observational areaof GK2A. Additionally, the present work shows that GEOS forecasts can be significantly improved,especially in the tropics and southern hemisphere after assimilating GK2A data in the absence ofHimawari-8 AMVs. This study demonstrates that GK2A AMV data is a valuable data source toenhance the robustness of GEOS ADAS.

Atmospheric Motion Vectors↗
NASA NTRSDate not supplied

CFD Model Development of a Cryogenic Storage Tank Self-Pressurization in Normal Gravity and Validation against SHIIVER Experiment

Two-phase flow and heat transfer simulations with interfacial phase change of the Structural Heat Intercept, Insulation, and Vibration Evaluation Rig (SHIIVER) self-pressurization experiment were conducted using storage tank CFD model in the framework of the ANSYS Fluent CFD code. The simulations were performed for the 70% fill level case with MLI on domes and no vapor cooling. All the phase change calculations in these simulations were generated by in-house Schrage-based evaporation-condensation model. The calculations were performed using both the Volume of Fluid (VOF) and Sharp Interface multiphase 2D axisymmetric models. A number of parametric and sensitivity studies were performed to check the various aspects of the CFD model. These studies helped to understand the effects of varying several parameters on the tank pressure and temperature during self-pressurization. Turbulence modeling; turbulence damping at the interface; and using constant vs. temperature dependent fluid properties were shown to have the most profound influence on predicted tank pressures and temperatures. Current study indicates that including phase change at the interface into the computational model is crucial for accurate prediction of the tank self-pressurization process. The effect of accommodation coefficient was also studied. Tank pressure values predicted by the VOF, and Sharp Interface models are within 4% of the experimental ones.

Computational Fluid Dynamics↗
NASA NTRSDate not supplied

CFD Model Development of a Cryogenic Storage Tank Self-Pressurization in Normal Gravity and Validation against SHIIVER Experiment

Two-Phase flow and heat transfer simulations with interfacial phase change of the Structural Heat Intercept, Insulation, and Vibration Evaluation Rig (SHIIVER) self-pressurization experiment were conducted using storage tank CFD model in the framework of the ANSYS Fluent CFD code. The simulations were performed for the 70% fill level case with MLI on domes and no vapor cooling. All the phase change calculations in these simulations were generated by in-house Schrage-based evaporation-condensation model. The calculations were performed using both the VOF and Sharp Interface multiphase 2Daxisymmetric models. A number of parametric and sensitivity studies were performed to check the various aspects of the CFD model. These studies helped to understand the effects of varying several parameters on the tank pressure and temperature during self-pressurization. These parameters include: the effect of turbulence modeling and turbulence damping at the interface; the effect of constant vs. temperature dependent fluid properties; the effect of accommodation coefficient; the effect of modeling phase change at the interface and the effect of the fill level.

Computational Fluid Dynamics↗
NASA NTRS2023

Retrieval of Normalized Gamma Size Distribution Parameters Using Precipitation Imaging Package (Pip) Snowfall Observations During Ice-Pop 2018

Parameters of the normalized gamma particle size distribution (PSD) have been retrieved from the Precipitation Image Package (PIP) snowfall observations collected during the International Collaborative Experiment - PyeongChang Olympics and Paralympic (ICE-POP 2018). Two of the gamma PSD parameters, the mass weighted particle diameter (Dmass) and the normalized intercept parameter NW, have median values of 1.15-1.31 mm and 2.84-3.04 log(mm-1 m-3), respectively. This range arises from the choice of the relationship between the maximum versus equivalent diameter, Dmx-Deq, and the relationship between the Reynolds and Best numbers, Re-X. Normalization of snow water equivalent rate (SWER) and ice water content (W) by NW reduces the range in NW resulting in well fitted power law relationship, between SWER/NW and Dmass and between W/NW and Dmass. The bulk descriptors of snowfall are calculated from PIP observations and from the gamma PSD with values of the shape parameter (μ) ranging from -2 to 10. NASA’s Global Precipitation Measurement (GPM) mission, which adopted the normalized gamma PSD, assumes μ = 2 and μ = 3 in its two separate algorithms. The mean fractional bias (MFB) of the snowfall parameters changes with μ, where the functional dependence on μ depends on the specific snowfall parameter of interest. The MFB of the total concentration was underestimated by 0.23-0.34 when μ = 2 and by 0.29-0.40 when μ = 3, while the MFB of SWER had a much narrower range (-0.03 to 0.04) for the same μ values.

Snowfall↗