Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “generator parameter estimation”

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,567 records · Page 87
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Sensorimotor integration in human postural control

It is generally accepted that human bipedal upright stance is achieved by feedback mechanisms that generate an appropriate corrective torque based on body-sway motion detected primarily by visual, vestibular, and proprioceptive sensory systems. Because orientation information from the various senses is not always available (eyes closed) or accurate (compliant support surface), the postural control system must somehow adjust to maintain stance in a wide variety of environmental conditions. This is the sensorimotor integration problem that we investigated by evoking anterior-posterior (AP) body sway using pseudorandom rotation of the visual surround and/or support surface (amplitudes 0.5-8 degrees ) in both normal subjects and subjects with severe bilateral vestibular loss (VL). AP rotation of body center-of-mass (COM) was measured in response to six conditions offering different combinations of available sensory information. Stimulus-response data were analyzed using spectral analysis to compute transfer functions and coherence functions over a frequency range from 0.017 to 2.23 Hz. Stimulus-response data were quite linear for any given condition and amplitude. However, overall behavior in normal subjects was nonlinear because gain decreased and phase functions sometimes changed with increasing stimulus amplitude. "Sensory channel reweighting" could account for this nonlinear behavior with subjects showing increasing reliance on vestibular cues as stimulus amplitudes increased. VL subjects could not perform this reweighting, and their stimulus-response behavior remained quite linear. Transfer function curve fits based on a simple feedback control model provided estimates of postural stiffness, damping, and feedback time delay. There were only small changes in these parameters with increasing visual stimulus amplitude. However, stiffness increased as much as 60% with increasing support surface amplitude. To maintain postural stability and avoid resonant behavior, an increase in stiffness should be accompanied by a corresponding increase in damping. Increased damping was achieved primarily by decreasing the apparent time delay of feedback control rather than by changing the damping coefficient (i.e., corrective torque related to body-sway velocity). In normal subjects, stiffness and damping were highly correlated with body mass and moment of inertia, with stiffness always about 1/3 larger than necessary to resist the destabilizing torque due to gravity. The stiffness parameter in some VL subjects was larger compared with normal subjects, suggesting that they may use increased stiffness to help compensate for their loss. Overall results show that the simple act of standing quietly depends on a remarkably complex sensorimotor control system.

NASA Discipline Neuroscience↗
DOE OSTI2022

Conversion Kinetics and Ionic Conductivity in Na-β”-Alumina + YSZ (Naβ”AY) Sodium Solid Electrolyte via Vapor Phase Conversion Process

Sodium ion batteries have been receiving increasing attention and may see potential revival in the near future, particularly in large-scale grid energy storage coupling with wind and solar power generation, due to the abundant sodium resources, low cost, and sufficiently high energy density. Among the known sodium ion conductors, the Na-β”-alumina electrolyte remains highly attractive because of its high ionic conductivity. This study focuses on the vapor phase synthesis of a Na-β”-Alumina + YSZ (Naβ”AY) composite sodium electrolyte, which has higher mechanical strength and stability than conventional single phase β”-Alumina. The objectives are the measurement of conversion kinetics through a newly developed weight-gain based model and the determination of sodium ionic conductivity in the composite electrolyte. Starting samples contained ~70 vol% α-Alumina and ~30 vol% YSZ (3 mol% Y 2 O 3 stabilized Zirconia) with and without a thin alumina surface layer made by sintering in air at 1600 °C. The sintered samples were placed in a powder of Na-β”-alumina and heat-treated at 1250 °C for various periods. Sample dimensions and weight were measured as a function of heat treatment time. The conversion of α-Alumina in the α-Alumina + YSZ composite into Naβ”AY occurred by coupled diffusion of sodium ions through Na-β”-alumina and of oxygen ions through YSZ, effectively diffusing Na2O. From the analysis of the time dependence of sample mass and dimensions, the effective diffusion coefficient of Na 2 O through the sample, D eff , was estimated to be 1.74 x 10 -7 cm2 s -1 , and the effective interface transfer parameter, k eff , was estimated as 2.33 x 10 -6 cm s -1 . By depositing a thin alumina coating layer on top of the bulk composite, the chemical diffusion coefficient of oxygen through single phase Na-β”-alumina was estimated as 4.35 x 10 -10 cm 2 s -1 . An AC impedance measurement was performed on a fully converted Naβ”AY composite, and the conductivity of the composite electrolyte was 1.3 x 10 -1 S cm -1 at 300 °C and 1.6 x 10 -3 S cm -1 at 25 °C, indicating promising applications in solid state or molten salt batteries at low to intermediate temperatures.

36 MATERIALS SCIENCE↗
NASA NTRS1991

Demonstration of rapid and sensitive module leak certification for Space Station Freedom

A leak detection and quantification demonstration using perflurocarbon tracer (PFT) technology was successfully performed at the NASA Marshall Space Flight Center on January 25, 1991. The real-time Dual Trap Analyzer (DTA) at one-half hour after the start of the first run gave an estimated leak rate of 0.7 mL/min. This has since been refined to be 1.15 (+ or -) 0.09 mL/min. The leak rates in the next three runs were determined to be 9.8 (+ or -) 0.7, -0.4 (+ or -) 0.3, and 76 (+ or -) 6 mL/min, respectively. The theory on leak quantification in the steady-state and time-dependent modes for a single zone test facility was developed and applied to the above determinations. The laboratory PFT analysis system gave a limit-of-detection (LOD) of 0.05 fL for ocPDCH. This is the tracer of choice and is about 100-fold better than that for the DTA. Applied to leak certification, the LOD is about 0.00002 mL/s (0.000075 L/h), a 5 order-of-magnitude improvement over the original leak certification specification. Furthermore, this limit can be attained in a measurement period of 3 to 4 hours instead of days, weeks, or months. A new Leak Certification Facility is also proposed to provide for zonal (three zones) determination of leak rates. The appropriate multizone equations, their solutions, and error analysis have already been derived. A new concept of seal-integrity certification has been demonstrated for a variety of controlled leaks in the range of module leak testing. High structural integrity leaks were shown to have a linear dependence of flow on (Delta)p. The rapid determination of leak rates at different pressures is proposed and is to be determined while subjecting the module to other external force-generating parameters such as vibration, torque, solar intensity, etc.

Dietz, R. N.↗
NASA NTRS1996

Cascade Error Projection: A Learning Algorithm for Hardware Implementation

In this paper, we workout a detailed mathematical analysis for a new learning algorithm termed Cascade Error Projection (CEP) and a general learning frame work. This frame work can be used to obtain the cascade correlation learning algorithm by choosing a particular set of parameters. Furthermore, CEP learning algorithm is operated only on one layer, whereas the other set of weights can be calculated deterministically. In association with the dynamical stepsize change concept to convert the weight update from infinite space into a finite space, the relation between the current stepsize and the previous energy level is also given and the estimation procedure for optimal stepsize is used for validation of our proposed technique. The weight values of zero are used for starting the learning for every layer, and a single hidden unit is applied instead of using a pool of candidate hidden units similar to cascade correlation scheme. Therefore, simplicity in hardware implementation is also obtained. Furthermore, this analysis allows us to select from other methods (such as the conjugate gradient descent or the Newton's second order) one of which will be a good candidate for the learning technique. The choice of learning technique depends on the constraints of the problem (e.g., speed, performance, and hardware implementation); one technique may be more suitable than others. Moreover, for a discrete weight space, the theoretical analysis presents the capability of learning with limited weight quantization. Finally, 5- to 8-bit parity and chaotic time series prediction problems are investigated; the simulation results demonstrate that 4-bit or more weight quantization is sufficient for learning neural network using CEP. In addition, it is demonstrated that this technique is able to compensate for less bit weight resolution by incorporating additional hidden units. However, generation result may suffer somewhat with lower bit weight quantization.

Duong, Tuan A.↗
NASA NTRSDate not supplied

L2-Charged Particle Environment (L2-CPE)Low Energy Radiation Fluence Model

The L2 Charged Particle Environment (L2-CPE) model provides estimates of number flux and fluence for the low energy electron, proton, and alpha particle populations in the near Earth solar wind and the Earth’s distant magnetosheath and magnetotail. The model is an engineering tool for specifying radiation environments over an energy range from a few eV to a few MeV of importance to surface dose and radiation damage to thin space exposed materials and is intended for use in space system design applications. Mission fluences are obtained by simulating a spacecraft flight trajectory through time-dependent bow shock and magnetopause boundaries with dimensions and orientations driven by solar wind parameters. Monte Carlo sampling of flux environments within individual plasma regimes and/or fluence accumulated along a fight trajectory through multiple regimes is used to determine statistical variations (means and extremes) of the differential number flux and fluence environments for each of the three charged particle species. Model output is differential (in energy) number fluence and flux for the three particle species in units of particles/cm2-keV and particles/cm2-sec-keV, respectively. Users can select the energy range of interest but the current version of the code (L2-CPE Version 1.4.2d) is limited to an energy range of 1 eV to 10 MeV. Flux is integrated over angle to give flux and fluence to surfaces in the ±XGSE, ±YGSE, and/or ±ZGSE directions. Figure 1 shows the opening screen from the L2-CPE graphical user interface (GUI) with options for flux, fluence, plotting model output. The flux scene generate is not currently implemented in the GUI.

Joseph I Minow↗
DOE OSTI2019

CSP Systems Analysis - Final Project Report

The CSP Systems Analysis project was a three-year effort supporting the Concentrating Solar Power (CSP) Subprogram within the Solar Energy Technologies Office of the U.S. Department of Energy (DOE). The goal of the CSP Systems Analysis project was to provide timely and accurate CSP cost data to the DOE and predict performance and cost for emerging technologies that may impact research directions and industry investment. The project consisted of six, mostly independent tasks. Major accomplishments over the three-year period of performance were as follows. The NREL team maintained current cost estimates for CSP technologies by annual updates to the cost input parameters within NREL's System Advisor Model (SAM). SAM is a performance and cost simulation software that includes technology models for the primary CSP technologies. During the period from 2016 to 2018, the estimated cost within SAM's default CSP molten-salt power tower technology dropped from 13.0 to 11.1 cent/kWh (real dollars without incentives). Power cycle models for the supercritical carbon dioxide (sCO2) Brayton Cycle were developed and made an option within SAM's molten-salt power tower model. These models were subsequently used to evaluate the cost reduction potential of alternative salts. The partial-cooling model was found to provide the greatest cost-reduction potential due to its combination of efficiency and temperature differential across the turbine and thermal energy storage (TES) system. These studies were documented in two journal articles. NREL's optical modeling software SolTrace and SolarPILOT were made more versatile and released as open-source software to the CSP community. These tools help researchers and developers to model and optimize the optical performance of CSP collector systems. The new features were highlighted in a recorded webinar that can be accessed from the NREL SAM website. New models for solar-thermal systems, that is, applications for solar industrial process heat, were developed, validated, and installed within SAM. The analysis includes a 2018 review of linear collectors in the market that are tailored for solar process heat applications. In 2018 the NREL team also developed an annual simulation for a direct steam generation system integrated with a latent-heat TES system. This analysis was submitted for journal publication. The NREL team continued to support the exchange of information and development of CSP guidelines within the international CSP community of SolarPACES. Much of the work under this project was presented at SolarPACES, including the development of a modeling checklist to help CSP stakeholders standardize the models used for estimating CSP performance and cost. The CSP Systems Analysis project also supported DOE requests for tracking and reporting programmatic targets for the CSP Subprogram.

14 SOLAR ENERGY↗
DOE OSTI2015

A spectroscopic study of the extreme black widow PSR J1311–3430

We report on a series of spectroscopic observations of PSR J1311–3430, an extreme black-widow gamma-ray pulsar with a helium-star companion. In a previous study we estimated the neutron star mass as M{sub NS}=2.68±0.14 M{sub ⊙} (statistical error), based on limited spectroscopy and a basic (direct heating) light-curve model; however, much larger model-dependent systematics dominate the mass uncertainty. Our new spectroscopy reveals a range of complex source behavior. The variable He i companion wind emission lines can dominate broadband photometry, especially in red filters or near minimum brightness, and the wind flux should complete companion evaporation in a spin-down time. The heated companion face also undergoes dramatic flares, reaching ∼40,000 K over ∼20% of the star; this is likely powered by a magnetic field generated in the companion. The companion center-of-light radial velocity is now well measured with K{sub CoL}=615.4±5.1 km s{sup −1}. We detect non-sinusoidal velocity components due to the heated face flux distribution. Using our spectra to excise flares and wind lines, we generate substantially improved light curves for companion continuum fitting. We show that the inferred inclination and neutron star mass, however, remain sensitive to the poorly constrained heating pattern. The neutron star’s mass, M{sub NS}, is likely less than the direct heating value and could range as low as 1.8 M{sub ⊙} for extreme equatorial heating concentration. While we cannot yet pin down M{sub NS}, our data imply that an intrabinary shock reprocesses the pulsar emission and heats the companion. Improved spectra and, especially, models that include such shock heating are needed for precise parameter measurement.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2003

SNC Oxygen Fugacity Recorded in Pyroxenes and its Implications for the Oxidation State of the Martian Interior: An Experimental and Analytical Study

Knowledge of the oxidation state of a magma is critical as it is one of the parameters which controls the nature and composition of the resulting crystals. In terrestrial magmatic systems, oxygen fugacity (fo2) is known to vary by over nine orders of magnitude. With variations of this magnitude, understanding the compositional differences, phase changes, and crystallization sequence variations, caused by the magma fo2, is essential in deciphering the origin of all igneous rocks. Magmatic oxidation state is of great importance in that it reflects the degree of oxidation of the source region and can provide insight into magmatic processes, such as metasomatism, degassing, and assimilation, which may have changed them. Carmichael [1991] argues that most magmas are unlikely to have their redox states altered from those of their source region. This assumption allows for estimation of the oxidation state of planetary interiors. Conversely, it is known that the fo2 of the magma can be affected by other processes, which occur outside of the source region and therefore, the oxidation state may record those too. Processes which could overprint source region fugacities include melt dehydrogenation or other volatile loss, water or melt infiltration, or assimilation of oxidized or reduced wallrock. Understanding which of these processes is responsible for the redox state of a magma can provide crucial information regarding igneous processes and other forces active in the region. The composition of the SNC basalts and their widely varying proposed oxidation states raise some interesting questions. Do the SNC meteorites have an oxidized or reduced signature? What was the oxygen fugacity of the SNC source region at the time of melt generation? Is the fugacity calculated for the various SNC samples the fugacity of the magma source region or was it overprinted by later events? Are there different oxidation states in the Martian interior or a single one? This proposal seeks to address all of these questions.

McCanta, M. C.↗
NASA NTRS2013

Interferometric Meteor Head Echo Observations using the Southern Argentina Agile Meteor Radar (SAAMER)

A radar meteor echo is the radar scattering signature from the free-electrons in a plasma trail generated by entry of extraterrestrial particles into the atmosphere. Three categories of scattering mechanisms exist: specular, nonspecular trails, and head-echoes. Generally, there are two types of radars utilized to detect meteors. Traditional VHF meteor radars (often called all-sky1radars) primarily detect the specular reflection of meteor trails traveling perpendicular to the line of sight of the scattering trail, while High Power and Large Aperture (HPLA) radars efficiently detect meteor head-echoes and, in some cases, non-specular trails. The fact that head-echo measurements can be performed only with HPLA radars limits these studies in several ways. HPLA radars are very sensitive instruments constraining the studies to the lower masses, and these observations cannot be performed continuously because they take place at national observatories with limited allocated observing time. These drawbacks can be addressed by developing head echo observing techniques with modified all-sky meteor radars. In addition, the fact that the simultaneous detection of all different scattering mechanisms can be made with the same instrument, rather than requiring assorted different classes of radars, can help clarify observed differences between the different methodologies. In this study, we demonstrate that such concurrent observations are now possible, enabled by the enhanced design of the Southern Argentina Agile Meteor Radar (SAAMER) deployed at the Estacion Astronomica Rio Grande (EARG) in Tierra del Fuego, Argentina. The results presented here are derived from observations performed over a period of 12 days in August 2011, and include meteoroid dynamical parameter distributions, radiants and estimated masses. Overall, the SAAMER's head echo detections appear to be produced by larger particles than those which have been studied thus far using this technique.

SPECULAR↗
DOE OSTI2020

Perturbation-Based Diagnosis of False Data Injection Attack Using Distributed Energy Resources

Modern smart grid relies on various sensor measurements for its operational control. In a successful false data injection attack, the attacker manipulates the measurements from the grid sensors such that undetected errors are introduced into the estimates of the system parameters leading to catastrophic situations. This paper proposes a novel perturbation based false data injection attack detection mechanism that utilizes inverter based distributed energy resources (DERs) to create low magnitude perturbation signal in the distribution system voltage that is inconsequential to the normal grid operation. Two voltage sensitivity analysis based algorithms are designed to identify the optimal set of DERs that can create the voltage perturbation signal of desired magnitude. An analytical method of voltage sensitivity analysis is used to compute the magnitude of voltage perturbation signal at each node in a computationally efficient manner. Then, a detection mechanism is developed that checks for the presence of the perturbation sequence in each sensor measurement. A sensor measurement is deemed authentic if the voltage perturbation signal is present in the data. In case of sensor malfunction or cyber-attack, the perturbation signal will not be present in the measurement data. Performance of the proposed attack detection mechanism is validated via simulation of the IEEE 69 bus test system.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2019

Testing modified gravity at cosmological distances with LISA standard sirens

Modifications of General Relativity leave their imprint both on the cosmic expansion history through a non-trivial dark energy equation of state, and on the evolution of cosmological perturbations in the scalar and in the tensor sectors. Specifically, the modification in the tensor sector gives rise to a notion of gravitational-wave (GW) luminosity distance, different from the standard electromagnetic luminosity distance, that can be studied with standard sirens at GW detectors such as LISA or third-generation ground based experiments. Here, we discuss the predictions for modified GW propagation from some of the best studied theories of modified gravity, such as Horndeski or the more general degenerate higher order scalar-tensor (DHOST) theories, non-local infrared modifications of gravity, bigravity theories and the corresponding phenomenon of GW oscillation, as well as theories with extra or varying dimensions. We show that modified GW propagation is a completely generic phenomenon in modified gravity. We then use a simple parametrization of the effect in terms of two parameters (Ξ 0 , n ), that is shown to fit well the results from a large class of models, to study the prospects of observing modified GW propagation using supermassive black hole binaries as standard sirens with LISA. We build mock source catalogs and perform detailed Markov Chain Monte Carlo studies of the likelihood obtained from LISA standard sirens alone, as well as by combining them with CMB, BAO and SNe data to reduce the degeneracies between cosmological parameters. We find that the combination of LISA with the other cosmological datasets allows one to measure the parameter Ξ 0 that characterizes modified GW propagation to the percent level accuracy, sufficient to test several modified gravity theories. LISA standard sirens can also improve constraints on GW oscillations induced by extra field content by about three orders of magnitude relative to the current capability of ground detectors. We also update the forecasts on the accuracy on H 0 and on the dark-energy equation of state using more recent estimates for the LISA sensitivity.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1997

Ground-Support Algorithms for Simulation, Processing, and Calibration of Barnes Static Earth Sensor Measurements: Applications to Tropical Rainfall Measuring Mission Observatory

New algorithms are described covering the simulation, processing, and calibration of penetration angles of the Barnes static Earth sensor assembly (SESA) as implemented in the Goddard Space Flight Center Flight Dynamics Division ground support system for the Tropical Rainfall Measuring Mission (TRMM) Observatory. The new treatment involves a detailed analysis of the measurements by individual quadrants. It is shown that, to a good approximation, individual quadrant misalignments can be treated simply as penetration angle biases. Simple formulas suitable for real-time applications are introduced for computing quadrant-dependent effects. The simulator generates penetration angles by solving a quadratic equation with coefficients uniquely determined by the spacecraft's position and the quadrant's orientation in GeoCentric Inertial (GCI) coordinates. Measurement processing for attitude determination is based on linearized equations obtained by expanding the coefficients of the aforementioned quadratic equation as a Taylor series in both the Earth oblateness coefficient (alpha approx. 1/150) and the angle between the pointing axis and the geodetic nadir vector. A simple formula relating a measured value of the penetration angle to the deviation of the Earth-pointed axis from the geodetic nadir vector is derived. It is shown that even near the very edge of the quadrant's Field Of View (FOV), attitude errors resulting from quadratic effects are a few hundredths of a degree, which is small compared to the attitude determination accuracy requirement (0.18 degree, 3 sigma) of TRMM. Calibration of SESA measurements is complicated by a first-order filtering used in the TRMM onboard algorithm to compute penetration angles from raw voltages. A simple calibration scheme is introduced where these complications are avoided by treating penetration angles as the primary raw measurements, which are adjusted using biases and scale factors. In addition to three misalignment parameters, the calibration state vector contains only two average penetration angle biases (one per each pair of opposite quadrants) since, because of the very narrow sensor FOV (+/- 2.6 degrees), differences between biases of the penetration angles measured by opposite quadrants cannot be distinguished from roll and pitch sensor misalignments. After calibration, the estimated misalignments and average penetration angle biases are converted to the four penetration angle biases and to the yaw misalignment angle. The resultant biases and the estimated scale factors are finally used to update the coefficients necessary for onboard computations of penetration angles from measured voltages.

Natanson, G. A.↗
DOE OSTI2025

Control of machining-induced residual stress via tool geometry and process parameter modification

Distortion generated in machined, monolithic, thin-walled aerospace components due to residual stresses leads to significant material and economic waste in the manufacturing industry. Inherent residual stress (IRS) present in stock materials combines with machining-induced residual stress (MIRS) to influence the final machined part distortion. It is hypothesized that MIRS can be controlled, based on the part geometry, through deliberate cutting tool geometry and process parameter modifications to negate the effect of IRS on distortion, consequently resulting in distortion-free parts. A finite element (FE) orthogonal cutting model is developed to study how tool geometry and process parameters influence machining-induced residual stress (MIRS). Orthogonal cutting experiments are performed on Al 7075-T651 samples to measure cutting forces and MIRS. A cutting force dynamometer is used to measure forces during cutting and a novel digital image correlation (DIC) based hole drilling technique is employed to measure the near-surface residual stress (RS) in the cut samples. These data are subsequently utilized to validate the FE prediction model. Various cases of cutting simulations involving different depths of cut, tool tip radii, and rake angles are performed to study their effects on RS. Similar to prior literature, increasing the depth of cut, tool tip radius, or rake angle is found to promote the formation of near-surface tensile stresses. The competing effects of material plowing and temperature are shown to determine the type of RS at the end of the cut. Moreover, a window of variation of RS (up to ± 400 MPa) is estimated within the given range of conditions, allowing for the control of MIRS through tool and process modification.

Mathews, Ritin [ORNL] (ORCID:0000000301440828)↗
DOE OSTI2017

Speleogenesis, geometry, and topology of caves: A quantitative study of 3D karst conduits

Karst systems are hierarchically spatially organized three-dimensional (3D) networks of conduits behaving as drains for groundwater flow. Recently, geostatistical approaches proposed to generate karst networks from data and parameters stemming from analogous observed karst features. Other studies have qualitatively highlighted relationships between speleogenetic processes and cave patterns. However, few studies have been performed to quantitatively define these relationships. This paper reports a quantitative study of cave geometries and topologies that takes the underlying speleogenetic processes into account. In order to study the spatial organization of caves, a 3D numerical database was built from 26 caves, corresponding to 621 km of cumulative cave passages representative of the variety of karst network patterns. The database includes 3D speleological surveys for which the speleogenetic context is known, allowing the polygenic karst networks to be divided into 48 monogenic cave samples and classified into four cave patterns: vadose branchwork (VB), water-table cave (WTC), looping cave (LC), and angular maze (AM). Eight morphometric cave descriptors were calculated, four geometrical parameters (width-height ratio, tortuosity, curvature, and vertical index) and four topological ones (degree of node connectivity, α and γ graph indices, and ramification index) respectively. The results were validated by statistical analyses (Kruskal-Wallis test and PCA). The VB patterns are clearly distinct from AM ones and from a third group including WTC and LC. A quantitative database of cave morphology characteristics is provided, depending on their speleogenetic processes. These characteristics can be used to constrain and/or validate 3D geostatistical simulations. This study shows how important it is to relate the geometry and connectivity of cave networks to recharge and flow processes. Conversely, the approach developed here provides proxies to estimate the evolution of the vadose zone to epiphreatic and phreatic zones in limestones from the quantitative analysis of existing cave patterns.

58 GEOSCIENCES↗
NASA NTRS2016

Cislunar Near Rectilinear Halo Orbit for Human Space Exploration

In order to conduct sustained human exploration beyond Low Earth Orbit (LEO), spacecraft systems are designed to operate in a series of missions of increasing complexity. Regardless of the destination, Moon, Mars, asteroids or beyond, there is a substantial set of common objectives that must be met. Many orbit characterization studies have endeavored to evaluate the potential locations in cislunar space that are favorable for meeting common human exploration objectives in a stepwise approach. Multiple studies, by both NASA and other international space agencies, have indicated that Earth-­‐moon libration point orbits are attractive candidates for staging operations in the proving ground and beyond. In particular, the Near Rectilinear Orbit (NRO) has been demonstrated to meet multi-­‐mission and multi-­‐destination architectural constraints. However, a human mission to a selected NRO presents a variety of new challenges for mission planning. While a growing number of robotic missions have completed successful operations to various specific libration point orbits, human missions have never been conducted to orbits of this class. Human missions have unique challenges that differ significantly from robotic missions, including a lower tolerance for mission risk and additional operational constraints that are associated only with human spacecraft. In addition, neither robotic nor human missions have been operated in the NRO regime specifically, and NROs exhibit dynamical characteristics that can differ significantly as compared to other halo orbits. Finally, multi-­‐body orbits, such as libration point orbits, are identified to exist in a simplified orbit model known as the Circular Restricted Three Body Problem (CRTBP) and must then be re-­‐solved in the full ephemeris model. As a result, the behavior of multi-­‐body orbits cannot be effectively characterized within the classical two-­‐body orbit dynamics framework more familiar to the human spaceflight community. In fact, a given NRO is not identified by a set of Keplerian orbit parameters, and a valid epoch specific state vector must be first obtained from a multi-body dynamical model. In this paper, the significant performance and operational challenges of conducting human missions to the NRO are evaluated. First, a systematic process for generating full ephemeris based ballistic NROs of various families is outlined to demonstrate the relative ease in which a multi-­‐revolution orbit can be found for any epoch and for various orbit geometries. In the Earth-­‐Moon system, NROs, which are halo orbits with close passage over a lunar pole, can exist with respect to libration point 1 (L1) or libration point 2 (L2) and are either from a North or South family orbit class with respect to the ecliptic. Second, the ability to maintain the orbit over the lifetime of a habitat mission by applying a reliable station-keeping strategy is investigated. The NRO, while similar to the quasi-­‐halo orbits that the Artemis mission flew, requires an updated station keeping strategy. This is due to several dynamical differences such as the increased relative stability of the NRO compared to other halo orbits and the close passage over the lunar surface as shown in Figure 1. Multiple station-keeping strategies are being investigated to ensure a human spacecraft remains on a predictable path. As the NRO is not described in simple two-­‐body parameters, analysis must determine the best strategy for targeting a reference NRO as well as how closely a future state should be constrained. In addition, costs will be minimized by determining maneuver directionality based on an identified pattern in the optimal station-keeping solutions or an analytically derived relationship. The candidate station-keeping algorithm must be stable and robust to environmental and vehicle uncertainties as well to navigation estimation and flight control execution errors. To that end, navigation accuracies, the impact on the station-keeping execution errors as well as other vehicle uncertainties need to be assessed. Starting with Orion, current navigation accuracies are evaluated and then navigation requirements are derived assuming a desired station-keeping propellant budget. Third, the performance requirements to and from the NRO are evaluated. Important parameters for developing expected propellant costs include epoch of operation, size and type of NRO, Earth departure and return constraints, as well as abort or early-­‐return capability. Finally, rendezvous and proximity operations are vital aspects of multi-­‐mission human exploration endeavors. The ability to conduct rendezvous and the associated propellant costs are assessed as well as the impacts of various profile assumptions including the location within the NRO the rendezvous is performed. The results of these studies will influence plans for international cooperation on both nearer term proving ground missions and beyond.

Whitley, Ryan↗
NASA NTRS2018

Numerical Study of Turbulent Separation Bubbles with Varying Pressure Gradient and Reynolds Number

A family of cases each containing a small separation bubble is treated by direct numerical simulation (DNS), varying two parameters: the severity of the pressure gradients, generated by suction and blowing across the opposite boundary, and the Reynolds number. Each flow contains a well-developed entry region with essentially zero pressure gradient, and all are adjusted to have the same value for the momentum thickness, extrapolated from the entry region to the centre of the separation bubble. Combined with fully defined boundary conditions this will make comparisons with other simulations and turbulence models rigorous; we present results for a set of eight Reynolds-averaged Navier–Stokes turbulence models. Even though the largest Reynolds number is approximately 5.5 times higher than in a similar DNS study we presented in 1997, the models have difficulties matching the DNS skin friction very closely even in the zero pressure gradient, which complicates their assessment. In the rest of the domain, the separation location per se is not particularly difficult to predict, and the most definite disagreement between DNS and models is near reattachment. Curiously, the better models tend to cluster together in their predictions of pressure and skin friction even when they deviate from the DNS, although their eddy-viscosity levels are widely different in the outer region near the bubble (or they do not rely on an eddy viscosity). Stratford’s square-root law is satisfied by the velocity profiles, both at separation and reattachment. The Reynolds-number range covers a factor of two, with the Reynolds number based on the extrapolated momentum thickness equal to approximately 1500 and 3000. This allows tentative estimates of the improvements that even higher values will bring to the model comparisons. The solutions are used to assess models through pressure, skin friction and other measures; the flow fields are also used to produce effective eddy-viscosity targets for the models, thus guiding turbulence-modelling work in each region of the flow.

Coleman, G. N.↗
DOE OSTI2016

Modeling and Analysis Methods for an On-line Enrichment Monitor

The International Atomic Energy Agency (IAEA) has developed an On-Line Enrichment Monitor (OLEM) as one possible component in a new generation of safeguards measures for uranium enrichment plants. The OLEM measures 235U emissions from the UF6 gas flowing through a unit header pipe using NaI(Tl) spectrometers, and corrects for gas density changes using pressure and temperature sensors in order to determine the enrichment of the gas as a function of time. In parallel with the OLEM instrument development, a Virtual OLEM (VOLEM) software tool has been developed that is capable of producing synthetic gamma-ray, pressure, and temperature data representative of a wide range of enrichment plant operating conditions. VOLEM complements instrument development activities and allows the study of OLEM for scenarios that will be difficult or impossible to evaluate empirically. Uses of VOLEM include: investigation of hardware design options; inter-comparison of candidate gamma-ray spectral analysis and enrichment estimation algorithms; uncertainty budget analysis and performance prediction for typical and atypical operational scenarios; and testing of the OLEM data acquisition, analysis and reporting software. This paper describes the technical foundations of VOLEM and illustrates how it can be used. An overview of the nominal instrument design and deployment scenario for OLEM is provided, with emphasis on the key online-assay measurement challenge: accurately determining the portion of the total 235U signal that comes from a background that includes solid uranium deposits on the piping walls. Monte Carlo modeling tools, data analysis algorithms and uncertainty quantification methods are described. VOLEM is then used to quantitatively explore the uncertainty budgets and predicted instrument performance for a plausible range of typical plant operating parameters, and one set of candidate analysis algorithms. Additionally, a series of VOLEM case studies illustrates how an online enrichment monitor might respond in scenarios expected to be atypical, for example weak pressure transients or significant changes in background levels.

Smith, Leon E.↗
DOE OSTI2017

Modeling CMB lensing cross correlations with CLEFT

A new generation of surveys will soon map large fractions of sky to ever greater depths and their science goals can be enhanced by exploiting cross correlations between them. In this paper we study cross correlations between the lensing of the CMB and biased tracers of large-scale structure at high z. We motivate the need for more sophisticated bias models for modeling increasingly biased tracers at these redshifts and propose the use of perturbation theories, specifically Convolution Lagrangian Effective Field Theory (CLEFT). Since such signals reside at large scales and redshifts, they can be well described by perturbative approaches. We compare our model with the current approach of using scale independent bias coupled with fitting functions for non-linear matter power spectra, showing that the latter will not be sufficient for upcoming surveys. Furthermore, we illustrate our ideas by estimating σ 8 from the auto- and cross-spectra of mock surveys, finding that CLEFT returns accurate and unbiased results at high z. We report on uncertainties due to the redshift distribution of the tracers, and several avenues for future development.

79 ASTRONOMY AND ASTROPHYSICS↗