Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “monotonicity”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
NASA NTRS2020

Global Search and Trajectory Optimization in Interplanetary and Cislunar Space

This presentation describes active research and development in interplanetary and cislunar trajectory optimization and global search at NASA Goddard Space Flight Center. Two point and parallel direct shooting transcriptions are described, along with monotonic basin hopping and batch seed sharing. Applications to the Lucy mission are presented, as well as a variety of other interplanetary and cislunar examples.

Englander, Jacob↗
NASA NTRSDate not supplied

Optimization of the Lunar Icecube Trajectory Using Stochastic Global Search and Multi-Point Shooting

Lunar IceCube is a 6U cubesat that will launch on NASA’s Artemis 1 mission in 2021. Lunar IceCube will separate from Artemis 1 shortly after trans-lunar injection (TLI) and travel to its science orbit about the moon using its Busek Ion Thruster 3 (BIT-3) propulsion system. This paper describes a technique to rapidly design Lunar IceCube trajectories using the monotonic basin hopping (MBH) stochastic global search algorithm, along with low- and high-fidelity multi-point shooting transcriptions. This technique allows the Lunar IceCube team to rapidly adapt to changing initial conditions, spacecraft properties, and operational constraints.

optimization↗
NASA NTRSDate not supplied

Hopping with an Adaptive Hop Probability Distribution

Monotonic Basin Hopping (MBH) is a stochastic global search technique that may be used to design complex interplanetary trajectories. Previous research on MBH empirically demonstrated that a bi-polar Pareto distribution is an effective way to generate search “hops.” However, no analytical foundation exists to explain this performance. In this work we provide the beginning of that analytical foundation and also introduce a new variant of MBH that uses an adaptive probability distribution to generate the hops. The technique is demonstrated on historical interplanetary trajectory design problems.

A C Englander↗
NASA NTRS2020

Changing Physical Conditions and Lower and Upper Trophic Level Responses on the US Northeast Shelf

Sea surface temperature (SST), salinity, and chlorophyll concentration (CHL) have changed in the US Northeast Shelf ecosystem over recent decades. The changes in these parameters were distinctly marked by change points around the year 2012 resulting in a 0.83°C increase in SST, a 0.3 PSU increase in salinity, and decrease in CHL in excess of 0.4 mg m–3. Where temperature and salinity shifted in mean level around their respective change points, CHL declined in a more monotonic fashion. Modeled data suggest that the shift in CHL resulted in a greater contribution of pico- and nanophytoplankton and a decreased contribution of microphytoplankton to overall CHL. Complementary estimates of the contribution of different phytoplankton functional types suggest a diminished contribution of diatoms to the phytoplankton community. Hence, not only is there evidence of a decline in the overall primary production capacity of the ecosystem, but also evidence of a fundamental change in the size and quality of phytoplankton supporting food webs. Two ecosystem responses to the observed changes in SST, salinity, and CHL were analyzed. Both length and weight at age have declined for a number of species, and both measures of growth appear to be negatively associated with temperature and positively associated with CHL. Biomass of fish and macroinvertebrates has declined in recent years, with a decrease in pelagic species associated with a decrease in CHL, while the decline in demersal species was associated with an increase in temperature. Collectively, these ecosystem changes appear to be the result of the complex interactions of both thermal effects and changes at the base of the food web.

Kevin D. Friedland↗
NASA NTRSDate not supplied

Variational Asymptotic Homogenization of Finitely Deformed Viscoelastic-Viscoplastic Composites

The objective of this paper is to develop a constitutive model for finitely deformed viscoelastic-viscoplastic materials and a micromechanics approach to homogenizing composites consisting of such materials. The development of the constitutive model involves establishing a thermodynamic framework based on finite strain theory, developing a viscoelasticity and a viscoplasticity model based on the thermodynamic framework, developing a radial return algorithm based on a classic framework, and deriving a closed-from incremental constitutive relation in the spatial configuration. The development of the micromechanics approach involves pulling-back the above constitutive relation to the material configuration, formulating a variational statement with the resulting constitutive relation, discretizing variational statement in a finite-dimensional space, and solving the discretized variational statement using an Euler–Newton predictor–corrector method. The constitutive model is calibrated via monotonic uniaxial compression tests on a polymer, and the calibrated model is validated by comparing its predictions with the cyclic test data. It is shown capable of characterizing viscoelasticity, viscoplasticity, and complex loading paths. The micromechanics approach’s capabilities are demonstrated through homogenizing a unidirectional fiber-reinforced composite, subjected to uniaxial, biaxial, and shear loading, at different strain rates. It is demonstrated to be capable of handling rate dependence and complex loading paths. The present framework can be further improved by implementing more sophisticated viscoelasticity and viscoplasticity models in future work.

Finite element analysis↗
NASA NTRS2002

Superconducting kinetic inductance photon detectors.

We are investigating a novel superconducting detector and readout method that could lead to photon counting, energy resolving focal plane arrays. This concept is intrinsically different from STJ and TES detectors, and in principle could deliver large pixel counts, high sensitivity, and Fano-limited spectral resolution in the ptical/UV/X-ray bands. The readout uses the monotonic relation between the kinetic surface inductance L***subs***, of a superconductor and the density of quasiparticles n, which holds even at temperatures far below T***subc***. This allows a sensitive readout of the number of excess quasiparticles in the detector by monitoring the transmission phase of a resonant circuit. The most intriguing aspect of this concept is that passive frequency multiplexing could be used to read out ~ l0***super4*** detectors with a single HEMT amplifier. Single x-ray events have been observed in prototype detectors.

kinetic↗
NASA NTRS2017

Inner Front Pole Cover Erosion in the 12.5 kW HERMeS Hall Thruster Over a Range of Operating Conditions

Erosion characteristics on the cover of the inner front pole in a 12.5 kW Hall thruster were measured over a wide range of operating conditions in tests of 6 to 14 hours duration using an accelerated test method and a very sensitive radioactive tracer-based erosion diagnostic. The operating points included the nominal 300 - 600 V conditions on a constant 20.8 A throttle curve, but included additional conditions at other currents spanning the throttling envelope and measurements at varying magnetic eld strength, facility pressure, and discharge voltage oscillation amplitude. The results show that the 300 V condition produces the highest wear rates on the 20.8 A throttle curve, but that rates actually increase with decreasing current. The wear rate was insensitive to discharge voltage ripple, but increased monotonically with magnetic eld strength, particularly near the inner radius of the pole cover. The inner region was also sensitive to facility pressure, showing lower rates at a higher pressure level. Separate experiments in which the energy distributions of ions generated by the hollow cathode were measured suggest that the cathode plume may be a source of energetic ions responsible for some of the erosion trends, in addition to ions originating in the thruster plume. The Hall thruster simulation code Hall2De is able to reproduce the erosion characteristics observed at 600 V, 20.8 A, but cannot currently match the rates measured at lower voltages and currents.

Lopez Ortega, Alejandro↗
NASA NTRS2020

Arctic Snow Depth and Sea Ice Thickness From ICESat‐2 and CryoSat‐2 Freeboards: A First Examination

We present a first examination of Arctic sea ice snow depth estimates from differencing satellite lidar (ICESat‐2) and radar (CryoSat‐2) freeboards. These estimates cover the period between 14 October 2018 and the end of April 2019. Snow depth is related to freeboard differences by the refractive index/bulk density of the snow layer—the only free parameter in the approach. Area‐averaged snow depth ranges from 9 cm (on first‐year ice: 5 cm, multiyear ice: 14 cm) in late October to 19 cm (first‐year ice: 17 cm, multiyear ice: 27 cm) in April; on average, this snow is thinner over FYI. Spatial patterns and gradients of snow depth estimates compare well with reconstructions using snowfall from ERA‐Interim and ERA5, although snowfall from ERA5 is systematically higher. For all months, the results suggest that ~50% of the total freeboard is comprised of snow. Retrievals are within a few centimeters of snow depth data acquired by Operation IceBridge in April 2019. Sources of uncertainties associated with this freeboard‐differencing approach are discussed. Further, sea ice thicknesses calculated using the retrieved snow depth and a modified climatology are contrasted. Comparatively, the snow depth and calculated ice thickness using a modified climatology are higher by ~5 cm and 0.33 m, although these differences are not uniform throughout the season. Snow accumulation was slower between October and December but increased between December and January, unlike the modified climatology, which exhibited a monotonic accumulation for all months. Future opportunities for assessment and improvement of these estimates are discussed.

Arctic Ocean↗
NASA NTRS2020

Differences in Tropical High Clouds Among Reanalyses: Origins and Radiative Impacts

We examine differences among reanalysis highcloud products in the tropics, assess the impacts of these differences on radiation budgets at the top of the atmosphere and within the tropical upper troposphere and lower stratosphere (UTLS), and discuss their possible origins in the context of the reanalysis models. We focus on the ERA5 (fifthgeneration European Centre for Medium-range Weather Forecasts – ECMWF – reanalysis), ERA-Interim (ECMWF Interim Reanalysis), JRA-55 (Japanese 55-year Reanalysis), MERRA-2 (Modern-Era Retrospective Analysis for Research and Applications, Version 2), and CFSR/CFSv2 (Climate Forecast System Reanalysis/Climate Forecast System Version 2) reanalyses. As a general rule, JRA-55 produces the smallest tropical high-cloud fractions and cloud water contents among the reanalyses, while MERRA-2 produces the largest. Accordingly, long-wave cloud radiative effects are relatively weak in JRA-55 and relatively strong in MERRA-2. Only MERRA-2 and ERA5 among the reanalyses produce tropical-mean values of outgoing long-wave radiation (OLR) close to those observed, but ERA5 tends to underestimate cloud effects, while MERRA-2 tends to overestimate variability. ERA5 also produces distributions of longwave, short-wave, and total cloud radiative effects at the top of the atmosphere that are very consistent with those observed. The other reanalyses all exhibit substantial biases in at least one of these metrics, although compensation between the long-wave and short-wave effects helps to constrain biases in the total cloud radiative effect for most reanalyses. The vertical distribution of cloud water content emerges as a key difference between ERA-Interim and other reanalyses. Whereas ERA-Interim shows a monotonic decrease of cloud water content with increasing height, the other reanalyses all produce distinct anvil layers. The latter is in better agreement with observations and yields very different profiles of radiative heating in the UTLS. For example, whereas the altitude of the level of zero net radiative heating tends to be lower in convective regions than in the rest of the tropics in ERAInterim, the opposite is true for the other four reanalyses. Differences in cloud water content also help to explain systematic differences in radiative heating in the tropical lower stratosphere among the reanalyses. We discuss several ways in which aspects of the cloud and convection schemes impact the tropical environment. Discrepancies in the vertical profiles of temperature and specific humidity in convective regions are particularly noteworthy, as these variables are directly constrained by data assimilation, are widely used, and feed back to convective behaviour through their relationships with thermodynamic stability.

High Clouds↗
NASA NTRSDate not supplied

Experimental Observations of Damage States in Unnotched and Notched 3D Orthogonal Woven Coupons Loaded in Tension

This paper describes a series of unnotched and notched tensile tests on two 3D orthogonal woven architectures. Dog bone specimens were used for the unnotched tests, and open-hole specimens were used for the notched tests. For each architecture and test type, two specimens were loaded under monotonic displacement control to failure and two specimens were incrementally loaded to failure in five steps. X-ray computed tomography (CT) was used to characterize the woven architecture and defects on each specimen prior to testing. Digital image correlation (DIC) images were obtained from both sides of the specimen throughout the loading and unloading. Acoustic emission (AE) data were also continuously acquired throughout the loading. For the incremental tests, X-ray CT was used to document the damage progression after each load increment. The damage progression was analyzed based on the results from these three techniques. The results for both architectures and test specimens followed a pattern where matrix cracks formed at the surface perpendicular to the tensile load direction. New surface cracks formed continuously throughout the loading and cracks extended in depth and width. In general, fiber failure was not detected until specimen failure.

Wade C Jackson↗
NASA NTRS2021

Deep Space Observations of Terrestrial Glitter

Deep space climate observatory (DSCOVR) spacecraft drifts about the Lagrangian point ≈1.4–1.6 × 106 km from Earth, where its Earth polychromatic imaging camera (EPIC) observes the sun-lit face of the Earth every 1 to 2 hours. At any instance, there is a preferred (specular) spot on the globe, where a glint may be observed by EPIC. While monitoring reflectance at these spots (terrestrial glitter), we observe occasional intense glints originating from neither ocean surface nor cloud ice and we argue that mountain lakes high in the Andes are among the causes. We also examine time-averaged reflectance at the spots and find it exceeding that of neighbors, with the excess monotonically increasing with separation distance. This specular excess is found in all channels and is more pronounced in the latest and best-calibrated version of EPIC data, thus opening the possibility of testing geometric calibration by monitoring distant glitter.

DSCOVR↗
NASA NTRS2021

Evolution of the Horizontal Winds in Jupiter’s Great Red Spot from One Jovian Year of HST/WFC3 Maps

We measured the horizontal winds in Jupiter's Great Red Spot (GRS) using data from the WFC3/UVIS instrument on board the Hubble Space Telescope (HST). The data cover 11 epochs from 2009 to 2020. Long-term monotonic trends in size and shape previously noted from the visible cloud appearance are paralleled by changes in the high-speed ring around the vortex. The circularization of the GRS cannot be explained by changes in the horizontal wind shear of the surrounding environment. The velocity fields suggest no long-term trend in the static stability inside or outside the vortex. Instead, the changes are accompanied by a 4%–8% increase in the mean wind speeds of the high-speed ring from 2009 to 2020. Changes in the wind field coincided with the South Equatorial Belt Outbreak storms of 2016–2017, but not with 2019 "flaking" events involving detachment of red material from the main oval.

Michael H. Wong↗
NASA NTRSDate not supplied

Experimental Observations of Damage States in Unnotched and Notched 3D Orthogonal Woven Coupons Loaded in Tension

This paper describes a series of unnotched and notched tensile tests on two 3D orthogonal woven architectures loaded in the warp direction. Dogbone specimens were used for the unnotched tests, and open-hole specimens were used for the notched tests. For each architecture and test type, two specimens were loaded under monotonic displacement control to failure, and two specimens were incrementally loaded to failure in five steps. X-ray computed tomography (CT) was used to characterize the woven architecture and defects on each specimen prior to testing. Digital image correlation (DIC) images were obtained from both sides of the specimen throughout the loading and unloading. For the incremental tests, X-ray CT was used to document the damage progression after each load increment. The results for both architectures and test specimens followed a similar pattern where matrix cracks formed perpendicular to the load direction and were primarily located between the stacks of warp tows. Two-piece failure path always occurred perpendicular to the loading direction, along the through-thickness segments of the Z-tows. The combination of stress concentrations due to the specimen geometry, weave geometry, and cut tows created stress concentrations that appear to initiate failure.

Woven materials↗
NASA NTRS2021

Salvaging Data Records with Missing Data: Data Imputation using the Multivariate t Distribution

When doing multivariate data analysis, one commonobstacle is the presence of incomplete observations, i.e., observationsfor which one or more key fields are blank. Missing datais often countered by deleting entire observations that containmissing data. The negative effects of deleting entire observationsare multiple: deleting observations reduces sample size andcan also result in biased inferences even if data is missing atrandom. In addition, knowledge contained within incompleteobservations is knowledge lost when they are deleted– and theeffort spent collecting that knowledge is effort wasted. Data imputationmethods, or methods of statistically “filling-in” missingdata, can help combat small sample sizes by using the existinginformation in partially complete observations with the end goalof producing less biased and higher confidence inferences. Whena sample from a multivariate normal population is only partiallycomplete, and the missing data meets appropriate assumptions(missing at random), robust data imputation of the missing datacan be implemented with monotone data augmentation (MDA)using the multivariate t distribution.Missing data imputation is applied to data from the NASA InstrumentCost Model (NICM) using the MDA algorithm underthe assumption of having a multivariate t distribution with fixeddegrees of freedom. A sensitivity analysis to the degrees offreedom parameter is presented to demonstrate robustness ofthe multivariate t distribution when dealing with small samplesas compared to the multivariate normal distribution.

DiNicola, Michael↗
NASA NTRS2022

The Origins of Long-Term Variability in Martian Upper Atmospheric Densities

We quantify and interpret the long-term variability of dayside Martian upper thermosphere and lower exosphere densities within 180–275 km altitudes. Atmospheric CO2, N2, O, and Ar densities are from NASA Mars Atmosphere and Volatile EvolutioN (MAVEN) observations during the time period of 2015–2020 near solar minimum. These neutral measurements, together with contemporaneous solar irradiance measurements at Mars, enable disentanglement of the orbital effect (due to the annual Sun-Mars distance change with solar longitude) and the solar extreme ultraviolet (EUV) effect in atmospheric density variations. The relative importance of these two effects, which is obtained using a statistical method of Dominance Analysis, reveals the competition between the indirect effect of solar infrared (via the upward coupling from the middle atmosphere) and the direct effect of solar EUV (due to local heating). Our results show that, unlike the orbital effect which is relatively constant at low altitudes and then decreases with increasing altitude, the solar EUV effect nearly monotonically increases. These two effects are comparable at high altitudes (about 240/270/205 km for CO2/N2/O). This analysis is extended to include long-term exospheric mass density estimates near 400 km from Mars Global Surveyor and Mars Odyssey data, with a focus on representative solar cycle phases of solar minimum and maximum. It is found that near 400 km, the orbital effect is always a key driver regardless of the solar cycle phase, while the solar EUV effect plays a minor role during solar minimum and is greatly enhanced and slightly exceeds the orbital effect during solar maximum.

Xiaohua Fang↗
NASA NTRS2022

How a Realistic Magnetosphere Alters the Polarizations of Surface, Fast Magnetosonic, and Alfvén Waves

System-scale magnetohydrodynamic (MHD) waves within Earth's magnetosphere are often understood theoretically using box models. While these have been highly instructive in understanding many fundamental features of the various wave modes present, they neglect the complexities of geospace such as the inhomogeneities and curvilinear geometries present. Here, we show global MHD simulations of resonant waves impulsively excited by a solar wind pressure pulse. Although many aspects of the surface, fast magnetosonic (cavity/waveguide), and Alfvén modes present agree with the box and axially symmetric dipole models, we find some predictions for large-scale waves are significantly altered in a realistic magnetosphere. The radial ordering of fast mode turning points and Alfvén resonant locations may be reversed even with monotonic wave speeds. Additional nodes along field lines that are not present in the displacement/velocity occur in both the perpendicular and compressional components of the magnetic field. Close to the magnetopause, the perpendicular oscillations of the magnetic field have the opposite handedness to the velocity. Finally, widely used detection techniques for standing waves, both across and along the field, can fail to identify their presence. We explain how all these features arise from the MHD equations when accounting for a non-uniform background field and propose modified methods that might be applied to spacecraft observations.

ULF Waves↗
NASA NTRSDate not supplied

Opposed-flow Spreading Flames: Effect of Sub-atmospheric Pressure on Spread and Burning Rates

Flame spread over solid fuels is a canonical problem in fire science, due to its direct implications on material flammability and importance in fire development. Flames in a microgravity environment can behave very differently than on Earth, posing additional risks for spaceflight life safety. Sub-atmospheric pressures in ground-based experiments have been proposed to approximately replicate the burning behavior of solid fuels in reduced gravity conditions because of similar effects on heat and mass transfer mechanisms from the flame to the solid. In opposed-flow flame spread, the solid fuel is heated by the flame ahead of its leading edge, and this process is strongly affected by the ambient conditions. In this work, we consider flames spreading over flat thin acrylic samples exposed to a forced flow of 20 cm/s, and pressures between 30 and 100 kPa. When radiation losses and chemical kinetic effects are neglected, experimental and theoretical results suggest an independence of the spread rate on pressure. However, the spread rate shows a decreasing trend with lower pressure in the range 30-60 kPa. The mass burning rate, calculated from the samples weight measured before and after the experiments, shows a similar trend. The flame length, on the other hand, shows a non-monotonic behavior, reaching a maximum at 60 kPa. Additionally, gas emissions measured during the experiments are used to estimate the heat release rate of the spreading flames. The comparison of the heat release rate obtained from the measured emissions and the estimated mass burning rates, suggests that chemical kinetics is not driving the decrease in flame spread observed at low pressures, whereas radiation losses could be responsible in analogy to what has been observed in micro and partial gravity studies. These results could provide more information to guide future Earth-based material flammability testing for spacecraft applications.

fire↗
NASA NTRS2022

Stability of the Dual-Frequency Radar Equations and a New Method Applied to the GPM’s Dual-Frequency Precipitation Radar (DPR) Data

A new algorithm is proposed that estimates two parameters of the particle size distribution (PSD) at each range bin from the Global Precipitation Measurement’s (GPM’s) Dual-frequency Precipitation Radar data. The equation that expresses the relationship of the PSD parameters between adjacent range bins is derived. By including the attenuation effect within the bin in the discretized equation, the new algorithm alleviates the double-solution problem when attenuation within the bin is sufficiently large. The stability of the solutions to the equation depends on the value of the mean diameter D(sub m) and its gradient with respect to range in the case of liquid precipitation. If the critical diameter above which the dual-frequency ratio (DFR) of radar reflectivity factors becomes a monotonically increasing function of D(sub m), the backward processing of the equation provides a stable or moderately diverging solution, unlike the forward processing that often gives unstable solutions. To provide a set of initial conditions without using the surface reference technique (SRT) in the backward processing, an initialization method using the Hitschfeld-Bordan attenuation correction method is proposed and tested. The proposed algorithm may provide a tool for investigating the assumptions used in various algorithms.

GPM↗