SEARCH · Engineering Papers
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.
Experimental work with leading edge separation bubbles is presented to clarify the issues regarding transition in separated regions. Hot-wire measurements, in the form of oscilloscope traces, turbulence intermittency and conditionally sampled velocity distributions are given. The resulting points of transition onset and spot production rates are compared to existing correlations.
A flight experiment was conducted to investigate the pressure distribution, local flow conditions, and boundary-layer transition characteristics on a large flat plate in flight at supersonic speeds up to Mach 2.0. The primary objective of the test was to characterize the local flow field in preparation for future tests of a high Reynolds number natural laminar flow test article. The tests used a F-15B testbed aircraft with a bottom centerline mounted test fixture. A second objective was to determine the boundary-layer transition characteristics on the flat plate and the effectiveness of using a simplified surface coating for future laminar flow flight tests employing infrared thermography. Boundary-layer transition was captured using an onboard infrared imaging system. The infrared imagery was captured in both analog and digital formats. Surface pressures were measured with electronically scanned pressure modules connected to 60 surface-mounted pressure orifices. The local flow field was measured with five 5-hole conical probes mounted near the leading edge of the test fixture. Flow field measurements revealed the local flow characteristics including downwash, sidewash, and local Mach number. Results also indicated that the simplified surface coating did not provide sufficient insulation from the metallic structure, which likely had a substantial effect on boundary-layer transition compared with that of an adiabatic surface. Cold wall conditions were predominant during the acceleration to maximum Mach number, and warm wall conditions were evident during the subsequent deceleration. The infrared imaging system was able to capture shock wave impingement on the surface of the flat plate in addition to indicating laminar-to-turbulent boundary-layer transition.
We present the results of a finite difference implementation of the kinetic Fokker-Planck model with an exact form of the nonlinear collisional operator, The model is time dependent and three-dimensional; one spatial dimension and two in velocity space. The spatial dimension is aligned with the local magnetic field, and the velocity space is defined by the magnitude of the velocity and the cosine of pitch angle. An important new feature of model, the concept of integration along the particle trajectories, is discussed in detail. Integration along the trajectories combined with the operator time splitting technique results in a solution scheme which accurately accounts for both the fast convection of the particles along the magnetic field lines and relatively slow collisional process. We present several tests of the model's performance and also discuss simulation results of the evolution of the plasma distribution for realistic conditions in Earth's plasmasphere under different scenarios.
We present new, high precision, high temporal resolution measurements of total column ozone (TCO) amounts derived from ground-based direct-sun irradiance measurements using our recently deployed Pandora single-grating spectrometers. Pandora's small size and portability allow deployment at multiple sites within an urban air-shed and development of a ground-based monitoring network for studying small-scale atmospheric dynamics, spatial heterogeneities in trace gas distribution, local pollution conditions, photochemical processes and interdependencies of ozone and its major precursors. Results are shown for four mid- to high-latitude sites where different Pandora instruments were used. Comparisons with a well calibrated double-grating Brewer spectrometer over a period of more than a year in Greenbelt MD showed excellent agreement and a small bias of approximately 2 DU (or, 0.6%). This was constant with slant column ozone amount over the full range of observed solar zenith angles (15-80), indicating adequate Pandora stray light correction. A small (1-2%) seasonal difference was found, consistent with sensitivity studies showing that the Pandora spectral fitting TCO retrieval has a temperature dependence of 1% per 3K, with an underestimation in temperature (e.g., during summer) resulting in an underestimation of TCO. Pandora agreed well with Aura-OMI (Ozone Measuring Instrument) satellite data, with average residuals of <1% at the different sites when the OMI view was within 50 km from the Pandora location and OMI-measured cloud fraction was <0.2. The frequent and continuous measurements by Pandora revealed significant short-term (hourly) temporal changes in TCO, not possible to capture by sun-synchronous satellites, such as OMI, alone.
Pedotransfer functions (PTFs) provide point predictions of soil hydraulic properties from more readily measured soil characteristics, yet uncertainties and biases in measurement methods, sampling distributions, and boundary conditions can limit accuracy when estimating near-saturated hydraulic conductivity (K(n)). These limitations may be particularly problematic in understudied urban landscapes that often contain altered hydraulic properties. To better treat deficiencies in PTF performance, we addressed three objectives, which were to: 1) develop PTFs to predict urban K(n), 2) assess bulk density and coarse fragments as explanatory variables; and 3) evaluate the predictive capability of these PTFs by comparing their output to measured hydraulic conductivity values from three other studies of urban soil hydraulics. We used artificial neural networks (ANN) and random forest (RF) approaches to predict urban K(n), with the training dataset including 307 tension infiltrometer tests and other measurements drawn from urban soil assessments in 11 U.S. cities. The PTFs utilized a hierarchy of inputs, starting with percentage sand, silt, clay, and then adding percentage coarse fragments and bulk density. The ANN models performed similar to the RF models, and all models exhibited similar or better predictive performance as models results collected from published articles. The inclusion of bulk density or coarse fragments did not improve accuracy over soil texture alone. Possible reasons for this result include low correlation between K(n) and bulk density and the exclusion of large voids during flow measurements with tension infiltrometers. The models have been made available as an open-source software package to encourage adoption by users working in urban systems.
There exists a positive correlation between the pH of subcellular compartments and the median isoelectric point (pI) for the associated proteomes. Proteins in the human lysosome—a highly acidic compartment in the cell—have a median pI of ∼6.5, whereas proteins in the more basic mitochondria have a median pI of ∼8.0. Proposed mechanisms reflect potential adaptations to pH. For example, enzyme active site general acid/base residue pKs are likely evolved to match environmental pH. However, such effects would be limited to a few residues on specific proteins, and might not affect the proteome at large. A protein model that considers residue burial upon folding recapitulates the correlation between proteome pI and environmental pH. This correlation can be fully described by a neutral evolution process; no functional selection is included in the model. Proteins in acidic environments incur a lower energetic penalty for burying acidic residues than basic residues, resulting in a net accumulation of acidic residues in the protein core. The inverse is true under alkaline conditions. The pI distributions of subcellular proteomes are likely not a direct result of functional adaptations to pH, but a molecular spandrel stemming from marginal stability.
This paper investigates induced oscillation (specifically subsynchronous torsional interactions) in power systems arising from photovoltaic (PV) power modulation. Time domain simulations of a modified IEEE first-benchmark system are presented and analyzed. In the IEEE benchmark system, a PV generator was added at the same bus where the synchronous generator is connected. The impact of different parameters of the PV inverters’ volt-volt-ampere reactive (volt-var) curve on SSTI is examined. The most severe impact found arises from the time response of the inverter’s reactive power change. A second-order time response with natural frequency matching the rotor speed deviation oscillation frequency and low-damping ratio quickly drives the system to an unstable region. Although it is highly unlikely that conditions necessary for distributed energy resource-induced oscillations would occur in the field, with knowledge of this vulnerability in mind, inverters' control parameters can be designed to avoid it.
Modified diglycolamide (DGA) extractants show high affinity for light lanthanides and improved separation factors compared to phosphonic acids used commercially. Previous studies with DGAs included flowsheet design and implementation into solvent extraction equipment but with low extractant concentration. Modifications on the alkyl chains of previous DGA extractants led to a competitive product for light rare earth separation with a higher extractant concentration. The scope of this project is to test a recently developed DGA extractant for application in pilot scale solvent extraction equipment with single stage testing, by narrowing down scrubbing conditions, and collecting distribution values for flowsheet planning.
This project demonstrated that solar PV can be recruited to stabilize the grid, smooth out disturbances, manage power flows, and assist circuit switching operations. It developed a radically new, layered control framework for Distributed Energy Resources (DER) to act in response to real-time, measured conditions on their local distribution circuit, rather than waiting for a price signal to indicate preferred behavior. By enabling resources to act as good citizens on the electric grid, Phasor-Based Control (PBC) facilitates arbitrarily high solar penetration levels.
This project demonstrated that solar PV can be recruited to stabilize the grid, smooth out disturbances, manage power flows, and assist circuit switching operations. It developed a radically new, layered control framework for Distributed Energy Resources (DER) to act in response to real-time, measured conditions on their local distribution circuit, rather than waiting for a price signal to indicate preferred behavior. By enabling resources to act as good citizens on the electric grid, Phasor-Based Control (PBC) facilitates arbitrarily high solar penetration levels.
Initial land cover distribution varies among Earth system models, an uncertainty in initial conditions that can substantially affect carbon and climate projections. We use the integrated Earth System Model to show that a 3.9 M km2 difference in 2005 global forest area (9–14% of total forest area) generates uncertainties in initial atmospheric CO2 concentration, terrestrial carbon, and local temperature that propagate through a future simulation following the Representative Concentration Pathway 4.5. By 2095, the initial 6 ppmv uncertainty range increases to 9 ppmv and the initial 26 PgC uncertainty range in terrestrial carbon increases to 33 PgC. The initial uncertainty range in annual average local temperature of -0.74 to 0.96 °C persists throughout the future simulation, with a seasonal maximum in Dec-Jan-Feb. These results highlight the importance of accurately characterizing historical land use and land cover to reduce overall initial condition uncertainty.
Conventional training methods for neural networks involve starting al a random location in the solution space of the network weights, navigating an error hyper surface to reach a minimum, and sometime stochastic based techniques (e.g., genetic algorithms) to avoid entrapment in a local minimum. It is further typically necessary to preprocess the data (e.g., normalization) to keep the training algorithm on course. Conversely, Bayesian based learning is an epistemological approach concerned with formally updating the plausibility of competing candidate hypotheses thereby obtaining a posterior distribution for the network weights conditioned on the available data and a prior distribution. In this paper, we developed a powerful methodology for estimating the full residual uncertainty in network weights and therefore network predictions by using a modified Jeffery's prior combined with a Metropolis Markov Chain Monte Carlo method.
Long-duration spaceflight results in sensorimotor adaptations, which cause functional deficits during gravitational transitions, such as landing on a planetary surface after long-duration microgravity exposure. Both the vestibular system and the central nervous system are affected by gravitational transitions. These systems are responsible for coordinating head and eye movements via the vestibulo-ocular reflex (VOR) and go through an adaptation period upon exposure to microgravity. Consequently, they must also re-adapt to Earth's gravitational environment upon landing. This re-adaptation causes decrements in gaze control and dynamic visual acuity, with crewmembers reporting oscillopsia and blurred vision caused by retinal slip, or the inability to keep an image focused on their retina. This is thought to drive motion sickness symptoms experienced by most crewmembers following landing. Retinal slip can be estimated by dynamic visual acuity (DVA); visual acuity while in motion. Previously, DVA has been assessed in the laboratory where subjects walked at 6.4 km/hr on a motorized treadmill. Using this method, Peters et al. (2011) found that DVA is worsened in astronauts by an average of 0.75 eye-chart lines one day after landing. However, it is believed that re-adaptation occurs quickly and that DVA might be worse immediately upon re-exposure to a gravitational environment. Since many crewmembers are unable to walk safely upon landing, it was necessary to develop a method for replicating the vertical head movements associated with walking. In addition, the use of a chair to imitate the head displacement caused by walking isolates eye-head interactions without allowing for trunk and lower-body compensation, as seen with treadmill walking (Mulavara & Bloomberg 2003). Therefore, a modality for assessing DVA in the field within a few hours of landing was developed. In this study, we validated the ability of a manually operated oscillating chair to reproduce the oscillatory frequency of walking on a treadmill. Healthy non-astronaut subjects (n=14) participated in one test session and completed three static (seated) and three dynamic (walking/oscillated) visual acuity tests. DVA was assessed using a motorized treadmill, an automated oscillating chair, and a manually operated chair, both developed in the Neuroscience Laboratory at JSC. The automated chair was motor-driven and set to oscillate vertically at 2 Hz with a vertical displacement of +/- 5 cm to simulate vertical translation while walking. The manually operated chair was oscillated vertically by a test operator to the beat of a metronome at 120 beats/min (2 Hz) and a vertical displacement of approximately +/- 5 cm. As the subject was oscillated, they were asked to discern the direction gap of Landolt-C optotypes of varying sizes and verbally reported the direction while an operator recorded their response using a gamepad. Subjects were outfitted with accelerometers (sampling rate = 128 Hz) on their head, trunk and lumbar spine. A fast Fourier transform was performed on the vertical trunk acceleration to compare the peak and spread of the distribution of oscillation frequencies for each oscillating condition. The spread of the frequency distribution for the manual chair was not significantly different from either the treadmill or the automated chair. However, all three conditions had similar non-zero standard error values, suggesting a variance in head movement frequency which may affect DVA. The average oscillation frequency of the manual chair (1.85 Hz) was significantly different (α=0.05) from that of treadmill walking (2.24 Hz), but not significantly different from that of the automated chair (1.85 Hz) and all three conditions had small standard errors (SEM = 0.04, 0.06, and 0.08 Hz for manual, treadmill, and automated respectively). This implies that both chairs oscillate at a frequency below that of treadmill walking, but are comparable to each other and reproducible across sessions. Additionally, DVA scores did not vary significantly across conditions. The smaller spread values of the oscillating chairs' frequencies indicated mitigation of variation induced by locomotor strategies, which enables better examination of the issue of VOR adaptation. Furthermore, due to the deconditioned state of crewmembers in the initial hours after landing, it is easier to transport a manual bouncing chair into the field and safer to perform a vision test while seated in a chair versus walking on a treadmill. Therefore, the manually oscillating chair has been deemed to meet and exceed the DVA testing capabilities previously obtained by treadmill walking.
This paper provides a description of procedures typically used to predict convective heating rates to hypersonic reentry vehicles using the two-layer method. These procedures were used to compute the pitch-plane heating distributions to the Apollo geometry for a wind tunnel test case and for three flight cases. Both simple engineering methods and coupled inviscid/boundary layer solutions were used to predict the heating rates. The sensitivity of the heating results in the choice of metrics, pressure distributions, boundary layer edge conditions, and wall catalycity used in the heating analysis were evaluated. Streamline metrics, pressure distributions, and boundary layer edge properties were defined from perfect gas (wind tunnel case) and chemical equilibrium and nonequilibrium (flight cases) inviscid flow-field solutions. The results of this study indicated that the use of CFD-derived metrics and pressures provided better predictions of heating when compared to wind tunnel test data. The study also showed that modeling entropy layer swallowing and ionization had little effect on the heating predictions.
This study is concerned with the analytical solutions of natural convection flows in closed-end cylindrical vessels to obtain exact solutions of temperature and velocity distribution in the laminar flow region under steady state condition. The temperature and velocity distributions, in general power series of displacements, are substituted into the three basic equations of continuity, momentum and energy. The relationships between coefficients of all powers are obtained through initial and boundary conditions and recurrence formulae. For constant wall temperature conditions the general solution of temperature and velocity may be expressed as a function of displacement, geometry ratio and Rayleigh number. Prandtl number does not enter as an independent parameter.
Flow measurements were performed on each of six swirl-can combustor module designs under simulated combustor operating conditions to find the design which exhibited a small recirculation zone, intense air mixing, and good fuel distribution in its wake. Conditions that are favorable for producing low oxides of nitrogen emissions and high combustion efficiency were investigated. The recirculation zone, the turbulence intensity and the fuel distribution pattern are obtained in the wake region of the center module of a three module array. The most promising swirl-can module design incorporates two air swirlers which discharge air in opposite directions (contraswirl), mixes the fuel and air upstream of the inner swirler, and has a flow area blockage of 64.3% for the three module array.
Necessary and sufficient condition derived for generating rational functions of n-port rectangular structure for distributed RC networks
Two interrelated aspects of the thermalization of cometary hydrogen were addressed using a three-dimensional Monte Carlo particle trajectory model: photochemical heating of the coma and the velocity distribution of cometary hydrogen leaving the inner coma that produces the extended Lyman-alpha coma. Observed cometary conditions were examined, ranging from productive comets like Kohoutek at very small heliocentric distances to intrinsically small comets like Giacobini-Zinner at 1 AU. The collisional decoupling of hydrogen yields bulk flow speeds and gas temperatures which are smaller than generally predicted in hydrodynamic models. By properly including these conditions the hydrogen velocity distribution required to produce the observed Lyman-alpha coma can be derived naturally. The two-dimensional morphology of the H coma is diagnostic of the physical conditions in the inner coma.