Search NASA⌕ Search

SEARCH · Search NASA

Results for “approximate control variates”

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 289 records · Page 16

Flow Effects on the Flammability Diagrams of Solid Fuels

A research program is currently underway with the final objective of developing a fundamental understanding of the controlling mechanisms underlying the flammability diagrams of solid combustible materials and their derived fire properties. Given that there is a high possibility of an accidental fire occurring in a space-based facility, understanding the fire properties of materials that will be used in such facilities is of critical importance. With this purpose, the flammability diagrams of the materials, as those produced by the Lateral Ignition and Flame Spread Test (LIFT) apparatus and by a new forced flow device, the Forced Flow Ignition and Flame Spread Test (FIST) apparatus, will be obtained. The specific objective of the program is to apply the new flammability apparatus, which will more accurately reflect the potential ambient conditions of space-based environments, to the characterization of the materials for space applications. This paper presents a parametric study of oxidizer flow effects on the ignition curve of the flammability diagrams of PMMA. The dependence of the ignition delay time on the external radiant flux and either the sample width (LIFT) or the flow velocity (FIST) has been studied. Although preliminary, the results indicate that natural and forced convection flow changes, affect the characteristics of the ignition curves of the flammability diagrams. The major effect on the ignition time appears to be due to convective transfer variations at the fuel surface. At high radiant fluxes or high flow velocities, however, it appears that gas phase processes become increasingly important, affecting the overall ignition delay time. A numerical analysis of the solid fuel heating and pyrolysis has also been developed. The theoretical predictions approximate the experiments well for conditions in which the gas phase induction time is negligible.

Cordova, J. L.↗

Implementation of Active Thermal Control (ATC) for the Soil Moisture Active and Passive (SMAP) Radiometer

NASA's Earth Observing Soil Moisture Active and Passive (SMAP) Mission is scheduled to launch in November 2014 into a 685 kilometer near-polar, sun-synchronous orbit. SMAP will provide comprehensive global mapping measurements of soil moisture and freeze/thaw state in order to enhance understanding of the processes that link the water, energy, and carbon cycles. The primary objectives of SMAP are to improve worldwide weather and flood forecasting, enhance climate prediction, and refine drought and agriculture monitoring during its three year mission. The SMAP instrument architecture incorporates an L-band radar and an L-band radiometer which share a common feed horn and parabolic mesh reflector. The instrument rotates about the nadir axis at approximately 15 revolutions per minute, thereby providing a conically scanning wide swath antenna beam that is capable of achieving global coverage within three days. In order to make the necessary precise surface emission measurements from space, the electronics and hardware associated with the radiometer must meet tight short-term (instantaneous and orbital) and long-term (monthly and mission) thermal stabilities. Maintaining these tight thermal stabilities is quite challenging because the sensitive electronics are located on a fast spinning platform that can either be in full sunlight or total eclipse, thus exposing them to a highly transient environment. A passive design approach was first adopted early in the design cycle as a low-cost solution. With careful thermal design efforts to cocoon and protect all sensitive components, all stability requirements were met passively. Active thermal control (ATC) was later added after the instrument Preliminary Design Review (PDR) to mitigate the threat of undetected gain glitches, not for thermal-stability reasons. Gain glitches are common problems with radiometers during missions, and one simple way to avoid gain glitches is to use the in-flight set point programmability that ATC affords to operate the radiometer component away from the problematic temperature zone. A simple ThermXL model (10 nodes) was developed to exercise quick trade studies among various proposed control algorithms: Modified P control vs. PI control. The ThermXL results were then compared with the detailed Thermal Desktop (TD) model for corroboration. Once done, the simple ThermXL model was used to evaluate parameter effects such as temperature digitization, heater size and gain margin, time step, and voltage variation of power supply on the ATC performance. A Modified P control algorithm was implemented into the instrument flight electronics based on the ThermXL results. The thermal short-term stability margin decreased by 10 percent with ATC and a wide temperature error band (plus or minus 0.1 degrees Centigrade) compared to the original passive thermal design. However, a tighter temperature error band (plus or minus 0.1 degrees Centigrade) increased the thermal short-term stability margin by a factor of three over the passive thermal design. The current ATC design provides robust thermal control, tighter stability, and greater in-flight flexibility even though its implementation was prompted by non-thermal performance concerns.

Mikhaylov, Rebecca↗

Average motion, structure and orientation of the distant magnetotail determined from remote sensing of the edge of the plasma sheet boundary layer with E greater than 35 keV ions

We study gradients of the energetic ion intesity observed at the edge of the plasma sheet boundary layer (PSBL) by the energetic ion anisotropy spectrometer (EPAS) on International Sun Earth Explorer 3 (ISEE 3). In particular, we have determined the velocity of the boundary relative to the spacecraft in the direction perpendicular to the tail axis and the angle which the boundary normal makes to the spacecraft spin axis for 1160 PSBL encounters at X(sub GSM) greater than -240 R(sub E). By asuming that, on average, the edge of the PSBL is parallel to the cross-tail current sheet, we are then able to determine a number of properties of the structure, orientation and motion of the deep geomagnetic tail. We conclude the following: (1) Most crossings of the edge of the PSBL are caused by transverse motuion of the entire tail induced by solar wind direction variations, although some are caused by reconfiguration of the tail due to geomagnetic activity. (2) The typical velocity of the PSBL (and hence of the tail) in the direction perpendicular to the tail axis is 50-85 km/s. (3) The average twist of the tail is near zero, with the edge of the PSBL (and by inference the cross-tail current sheet) lying parallel to the ecliptic plane (however, large twists are found in individual events and the distribution of twists is broad, with one standard deviation of approximately 50 deg. (4) The width of the distribution decreases with downtail distance. (5) The variation of the distributions with cross-tail position reveals that this decreas in width is most likely due to the edge of the PSBL being concave, or significantly flared at the tail flanks, in the near-Earth region. (6) During days on which the Interplanetary Magnetic Field (IMF) has 'away' sector structure, the north lobe of the trail is twisted on average towards dawn by 7.0 +/-2.4 deg. (7) During days on which the IMF has 'toward' sector structure, the north lobe is tilted towards duskby 3.8 +/- 2.3. (8) A subset of events for which IMP 8 solar wind data are available show that, for southward IMF BH(sub z) the tail has a mean twist of -12.3 +/- 5.0 deg for IMF B(sub Y) greater than 0 and 5.5 +/- 3.8 deg for IMF B(sub Y) less than 0 (positive twist angles correspond to a tilt of the northern lobe towards dusk). (9) For northward IMF B(sub z) the tail has a twist of -23.9 +/- 5.0 deg for IMF B(sub Y) greater than 0 and 13.4 +/- 6.0 deg for IMF B(sub y) less than 0. Hence the tail appears more twisted on average for the IMF B(sub Z) northward case. (10) The distribution of tail twist is wider for lower levels of geomagnetic activity, indicating that the tail is able to twist more at lower levels of activity. (11) The data set reveals no evident effect of the earth's dipole wobble; tail orientation appears to be controlled by the solar wind and IMF, such that the GSE coordinate system may be appropriate for the study of field and plasma structures in the distant tail region.

Owen, C. J.↗

Exploiting the Temperature Dependence of Magnetic Susceptibility to Control Convection in Fundamental Studies of Solidification Phenomena

It is well known that convection is a dominant mass transport mechanism when materials are solidified on Earth's surface. This convection is caused by gradients in density (and therefore gravitational force) that are brought about by gradients in temperature, composition or both. Diffusion of solute is therefore dwarfed by convection and the study of fundamental parameters, such as dendrite tip shape and growth velocity in the absence of convection is nearly impossible. Significant experimental work has therefore been carried out in orbiting laboratories with the intent of minimizing convection by minimizing gravity. One of the best known experiments of this kind is the Isothermal Dendritic Growth Experiment (IDGE), supported by NASA. Naturally such experiments are costly and one objective of the present investigation is to develop an experimental method whereby convection can be- halted, in solidification and other experiments, on the surface. A second objective is to use the method to minimize convection resulting from the residual accelerations suffered by experiments in microgravity. The method to be used to minimize convection relies on the dependence of the magnetic susceptibility of a fluid on temperature or composition (whichever is driving convection). All materials experience a force when placed in a magnetic field gradient. The direction and magnitude of that force depend on the magnetic susceptibility of the material. Consequently the force will vary if the susceptibility varies with temperature or composition. With a magnetic field gradient in the right direction (typically upward) and of the right magnitude, this variation in the magnetic force can be made to exactly cancel the variation in the gravitational force. Expressed another way, normal buoyancy is exactly countered by a "magnetic buoyancy". To demonstrate the principle, a solution of MnC12 in water has been used. First the variation of the susceptibility of this paramagnetic solution with temperature and concentration was measured. Then a "cell", containing this solution and 50mm long by 15mm high by 155mm wide, was placed in a superconducting magnet at Marshall Space Flight Center. The magnetic field was measured at various positions within the bore of the magnet using a Hall effect probe. In this way, a position was found where the magnetic field gradient was predominantly upward; the magnitude of the gradient could then be adjusted by adjusting the current of the magnet. The ends of the cell consisted of machined copper blocks maintained at controlled temperatures by circulating water from constant temperature baths. The walls of the cell were of rectangular section glass tubing so that the cell contents could be seen. Velocities arising from thermal gradients within the cell were measured by particle image velocimetry (PIV). Particles used for this purpose were silver-coated hollow glass spheres of micrometers diameter and nearly the same density as the solution. A central vertical plane of the cell was illuminated by a laser beam passing through a cylindrical lens. Digital images of the particles were captured on a CCD camera and fed to a computer so that frame-to-frame movements of particles traveling with the fluid were captured. These images were employed to compute velocity maps using commercial PIV software. In a typical experiment the cold end of the cell was maintained at 10C and the warm end at 30 C. With no current in the magnet, i.e.- with natural -convection allowed to occur, the fluid was observed to circulate with an average speed of approximately 0.3 millimeters per second. It was visually apparent that this circulation was diminished as the current was increased. At currents of approximately 20A the flow was halted, to within the precision of the PIV measurements. At yet higher currents the convection was reversed with the hotter solution sinking and the cooler solution rising. At 40A this reversed convection had speeds averaging 0.43 millimeters per second. The measurements of susceptibility and density allow an estimate of the field gradient necessary to halt convection in the experiment. That estimate was 7.8T (squared) per meter and the convection was observed to halt in the magnet at a current giving 7.21T (squared) per meter from the magnetic field measurements. Calculations of the flow have been carried out using the computational fluid dynamics software FLUENT and show good agreement with the measurements.

Seybert, C. D.↗

Counterdiabatic Driving with Performance Guarantees

Counterdiabatic (CD) driving has the potential to speed up adiabatic quantum state preparation by suppressing unwanted excitations. However, existing approaches either require intractable classical computations or are based on approximations that do not have performance guarantees. We propose and analyze a nonvariational, system-agnostic CD expansion method and analytically show that it converges exponentially quickly in the expansion order. In finite systems, the required resources scale inversely with the spectral gap, which we argue is asymptotically optimal. To extend our method to the thermodynamic limit and suppress errors stemming from high-frequency transitions, we leverage finite-time adiabatic protocols. In particular, we show that a time determined by the quantum speed limit is sufficient to prepare the desired ground state, without the need to optimize the adiabatic trajectory. Numerical tests of our method on the quantum Ising chain show that our method can outperform state-of-the-art variational CD approaches.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Lead-Lag Control for Helicopter Vibration and Noise Reduction

As a helicopter transitions from hover to forward flight, the main rotor blades experience an asymmetry in flow field around the azimuth, with the blade section tangential velocities increasing on the advancing side and decreasing on the retreating side. To compensate for the reduced dynamic pressure on the retreating side, the blade pitch angles over this part of the rotor disk are increased. Eventually, a high enough forward speed is attained to produce compressibility effects on the advancing side of the rotor disk and stall on the retreating side. The onset of these two phenomena drastically increases the rotor vibratory loads and power requirements, thereby effectively establishing a limit on the maximum achievable forward speed. The alleviation of compressibility and stall (and the associated decrease in vibratory loads and power) would potentially result in an increased maximum forward speed. In the past, several methods have been examined and implemented to reduce the vibratory hub loads. Some of these methods are aimed specifically at alleviating vibration at very high flight speeds and increasing the maximum flight speed, while others focus on vibration reduction within the conventional flight envelope. Among the later are several types passive as well as active schemes. Passive schemes include a variety of vibration absorbers such as mechanical springs, pendulums, and bifilar absorbers. These mechanism are easy to design and maintain, but incur significant weight and drag penalties. Among the popular active control schemes in consideration are Higher Harmonic Control (HHC) and Individual Blade Control (IBC). HHC uses a conventional swash plate to generate a multi-cyclic pitch input to the blade. This requires actuators capable of sufficiently high power and bandwidth, increasing the cost and weight of the aircraft. IBC places actuators in the rotating reference frame, requiring the use of slip rings capable of transferring enough power to the actuators. Both schemes cause an increase in pitch link loads. Trailing Edge Flap (TEF) deployment can also used to generate unsteady aerodynamic forces and moments that counter the original vibratory loads, and thereby reduce rotor vibrations. While the vibrations absorbers, HHC, IBC, and TEF concepts discussed above attempt to reduce the vibratory loads, they do not specifically address the phenomena causing the vibrations at high advance ratios. One passive method that attempts to directly alleviate compressibility and stall, instead of reducing the ensuing vibrations, is the use of advanced tip designs. Taper, sweep, anhedral, and the manipulation of other geometric properties of the blade tips can reduce the severity of stall and compressibility effects , as well as reduce rotor power. A completely different approach to solve these problems is the tiltrotor configuration. As the forward velocity of the aircraft increases, the rotors, in this case, are tilted forward until they are perpendicular to the flow and act as propellers. This eliminates the edgewise flow encountered by conventional rotors and circumvents all the problems associated with flow asymmetry. However, the success involves a tremendous increase in cost and complexity of the aircraft. Another possible approach that has been proposed for the alleviation of vibratory loads at high forward flight speeds involves the use of controlled lead-lag motions to reduce the asymmetry in flow. A correctly phased 1/rev controlled lag motion could be introduced such that it produces a backward velocity on the advancing side and a forward velocity on the retreating side, to delay compressibility effects and stall to a higher advance ratio. Using a large enough lead-lag amplitude, the tip velocities could be reduced to levels encountered in hover. This concept was examined by two groups in the 1950's and early 1960's. In the United States, the Research Labs Division of United Aircraft developed a large lead-lag motion rotor, meant to achieve lag motion amplitudes up to 45 degrees. In order to reduce the required actuation force, the blade hinges were moved to 40% of the blade radius to increase the rotating lag frequency to approximately 1/rev. The blade hinges were redesigned to produce a flap-lag coupling so the large flapwise aerodynamic loads could be exploited to actuate the blades in the lag direction. A wind tunnel test of this rotor concept revealed actuation and blade motion scheduling problems. The project was eventually discontinued due to these problems and high blade stresses. Around the same time, at Boelkow in Germany, a similar lead-lag rotor program was conducted under the leadership of Hans Derschmidt. Here, too, the blade hinges were moved outboard to 34% radius to reduce the actuation loads. The main difference between this and the United Aircraft program was the use of a mechanical actuation scheme with maximum lead-lag motions of 400. This program was also discontinued for unclear reasons. The present study is directed toward conducting a comprehensive analytical examination to evaluate the effectiveness of controlled lead-lag motions in reducing vibratory hub loads and increasing maximum flight speed. Since both previous studies on this subject were purely experimental, only a limited data set and physical understanding of the problem was obtained. With the currently available analytical models and computational resources, the present effort is geared toward developing an in-depth physical understanding of the precise underlying mechanisms by which vibration reduction may be achieved. Additionally, in recognition of the fact that large amplitude lead-lag motions would - (i) be difficult to implement, and (ii) produce very large blade stresses; the present study examines the potential of only moderate-to-small lead-lag motions for reduction of vibratory hub loads. Using such an approach, the emphasis is not on eliminating the periodic variations in tangential velocity at the blade tip, but at best reducing these variations slightly so that compressibility and stall are delayed to slightly higher advance ratios. This study was conducted in two steps. In the first step, a hingeless helicopter rotor was modeled using rigid blades undergoing flap-lag-torsion rotations about spring restrained hinges and bearings. This model was then modified by separating the lead-lag degree of freedom into two components, a free and a prescribed motion. Using this model, a parametric study of the effect of phase and amplitude of a prescribed lead-lag motion on hub vibration was conducted. The data gathered was analyzed to obtain an understanding of the basic physics of the problem and show the capability of this method to reduce vibration and expand the flight envelope. In the second half of the study, the similar analysis was conducted using an elastic blade model to confirm the effects predicted by the simpler model.

Gandhi, Farhan↗

AeroCom Phase III Multi-Model Evaluation of the Aerosol Life Cycle and Optical Properties Using Ground and Space-Based Remote Sensing as Well as Surface In Situ Observations

Within the framework of the AeroCom (Aerosol Comparisons between Observations and Models) initiative, the state-of-the-art modelling of aerosol optical properties is assessed from 14 global models participating in the phase III control experiment (AP3). The models are similar to CMIP6/AerChemMIP Earth System Models (ESMs) and provide a robust multi-model ensemble. Inter-model spread of aerosol species lifetimes and emissions appears to be similar to that of mass extinction coefficients (MECs), suggesting that aerosol optical depth (AOD) uncertainties are associated with a broad spectrum of parameterised aerosol processes. Total AOD is approximately the same as in AeroCom phase I (AP1) simulations. However, we find a 50 % decrease in the optical depth (OD) of black carbon (BC), attributable to a combination of decreased emissions and lifetimes. Relative contributions from sea salt (SS) and dust (DU) have shifted from being approximately equal in AP1 to SS contributing about 2∕3 of the natural AOD in AP3. This shift is linked with a decrease in DU mass burden, a lower DU MEC, and a slight decrease in DU lifetime, suggesting coarser DU particle sizes in AP3 compared to AP1. Relative to observations, the AP3 ensemble median and most of the participating models underestimate all aerosol optical properties investigated, that is, total AOD as well as fine and coarse AOD (AODf, AODc), Ångström exponent (AE), dry surface scattering (SCdry), and absorption (ACdry) coefficients. Compared to AERONET, the models underestimate total AOD by ca. 21 % ± 20 % (as inferred from the ensemble median and interquartile range). Against satellite data, the ensemble AOD biases range from −37 % (MODIS-Terra) to −16 % (MERGED-FMI, a multi-satellite AOD product), which we explain by differences between individual satellites and AERONET measurements themselves. Correlation coefficients (R) between model and observation AOD records are generally high (R>0.75), suggesting that the models are capable of capturing spatio-temporal variations in AOD. We find a much larger underestimate in coarse AODc (∼ −45 % ± 25 %) than in fine AODf (∼ −15 % ± 25 %) with slightly increased inter-model spread compared to total AOD. These results indicate problems in the modelling of DU and SS. The AODc bias is likely due to missing DU over continental land masses (particularly over the United States, SE Asia, and S. America), while marine AERONET sites and the AATSR SU satellite data suggest more moderate oceanic biases in AODc. Column AEs are underestimated by about 10 % ± 16 %. For situations in which measurements show AE > 2, models underestimate AERONET AE by ca. 35 %. In contrast, all models (but one) exhibit large overestimates in AE when coarse aerosol dominates (bias ca. +140 % if observed AE < 0.5). Simulated AE does not span the observed AE variability. These results indicate that models overestimate particle size (or underestimate the fine-mode fraction) for fine-dominated aerosol and underestimate size (or overestimate the fine-mode fraction) for coarse-dominated aerosol. This must have implications for lifetime, water uptake, scattering enhancement, and the aerosol radiative effect, which we can not quantify at this moment. Comparison against Global Atmosphere Watch (GAW) in situ data results in mean bias and inter-model variations of −35 % ± 25 % and −20 % ± 18 % for SCdry and ACdry, respectively. The larger underestimate of SCdry than ACdry suggests the models will simulate an aerosol single scattering albedo that is too low. The larger underestimate of SCdry than ambient air AOD is consistent with recent findings that models overestimate scattering enhancement due to hygroscopic growth. The broadly consistent negative bias in AOD and surface scattering suggests an underestimate of aerosol radiative effects in current global aerosol models. Considerable inter-model diversity in the simulated optical properties is often found in regions that are, unfortunately, not or only sparsely covered by ground-based observations. This includes, for instance, the Sahara, Amazonia, central Australia, and the South Pacific. This highlights the need for a better site coverage in the observations, which would enable us to better assess the models, but also the performance of satellite products in these regions. Using fine-mode AOD as a proxy for present-day aerosol forcing estimates, our results suggest that models underestimate aerosol forcing by ca. −15 %, however, with a considerably large interquartile range, suggesting a spread between −35 % and +10 %.

space-based remote sensing↗

Investigation of a geodesy coexperiment to the Gravity Probe B relativity gyroscope program

Geodesy is the science of measuring the gravitational field of and positions on the Earth. Estimation of the gravitational field via gravitation gradiometry, the measurement of variations in the direction and magnitude of gravitation with respect to position, is this dissertation's focus. Gravity Probe B (GP-B) is a Stanford satellite experiment in gravitational physics. GP-B will measure the precession the rotating Earth causes on the space time around it by observing the precessions of four gyroscopes in a circular, polar, drag-free orbit at 650 km altitude. The gyroscopes are nearly perfect niobium-coated spheres of quartz, operating at 1.8 K to permit observations with extremely low thermal noise. The permissible gyroscope drift rate is miniscule, so the torques on the gyros must be tiny. A drag-free control system, by canceling accelerations caused by nongravitational forces, minimizes the support forces and hence torques. The GP-B system offers two main possibilities for geodesy. One is as a drag-free satellite to be used in trajectory-based estimates of the Earth's gravity field. We described calculations involving that approach in our previous reports, including comparison of laser only, GPS only, and combined tracking and a preliminary estimate of the possibility of estimating relativistic effects on the orbit. The second possibility is gradiometry. This technique has received a more cursory examination in previous reports, so we concentrate on it here. We explore the feasibility of using the residual suspension forces centering the GP-B gyros as gradiometer signals for geodesy. The objective of this work is a statistical prediction of the formal uncertainty in an estimate of the Earth's gravitation field using data from GP-B. We perform an instrument analysis and apply two mathematical techniques to predict uncertainty. One is an analytical approach using a flat-Earth approximation to predict geopotential information quality as a function of spatial wavelength. The second estimates the covariance matrix arising in a least-squares estimate of a spherical harmonic representation of the geopotential using GP-B gradiometer data. The results show that the GP-B data set can be used to create a consistent estimate of the geopotential up to spherical harmonic degree and order 60. The formal uncertainty of all coefficients between degrees 5 and 50 is reduced by factors of up to 30 over current satellite-only estimates and up to 7 over estimates which include surface data. The primary conclusion resulting from this study is that the gravitation gradiometer geodesy coexperiment to GP-B is both feasible and attractive.

Everitt, C. W. F.↗

Experimental Measurement of RCS Jet Interaction Effects on a Capsule Entry Vehicle

An investigation was made in NASA Langley Research Center s 31-Inch Mach 10 Tunnel to determine the effects of reaction-control system (RCS) jet interactions on the aft-body of a capsule entry vehicle. The test focused on demonstrating and improving advanced measurement techniques that would aid in the rapid measurement and visualization of jet interaction effects for the Orion Crew Exploration Vehicle while providing data useful for developing engineering models or validation of computational tools used to assess actual flight environments. Measurements included global surface imaging with pressure and temperature sensitive paints and three-dimensional flow visualization with a scanning planar laser induced fluorescence technique. The wind tunnel model was fabricated with interchangeable parts for two different aft-body configurations. The first, an Apollo-like configuration, was used to focus primarily on the forward facing roll and yaw jet interactions which are known to have significant aft-body heating augmentation. The second, an early Orion Crew Module configuration (4-cluster jets), was tested blowing only out of the most windward yaw jet, which was expected to have the maximum heating augmentation for that configuration. Jet chamber pressures and tunnel flow conditions were chosen to approximate early Apollo wind tunnel test conditions. Maximum heating augmentation values measured for the Apollo-like configuration (>10 for forward facing roll jet and 4 for yaw jet) using temperature sensitive paint were shown to be similar to earlier experimental results (Jones and Hunt, 1965) using a phase change paint technique, but were acquired with much higher surface resolution. Heating results for the windward yaw jet on the Orion configuration had similar augmentation levels, but affected much less surface area. Numerical modeling for the Apollo-like yaw jet configuration with laminar flow and uniform jet outflow conditions showed similar heating patterns, qualitatively, but also showed significant variation with jet exit divergence angle, with as much as 25 percent variation in heat flux intensity for a 10 degree divergence angle versus parallel outflow. These results along with the fabrication methods and advanced measurement techniques developed will be used in the next phase of testing and evaluation for the updated Orion RCS configuration.

Buck, Gregory M.↗

Flight Investigation of Effect of Various Vertical-Tail Modifications on the Directional Stability and Control Characteristics of the P-63A-1 Airplane (AAF No. 42-68889)

Because the results of preliminary flight tests had indicated. the P-63A-1 airplane possessed insufficient directional stability, the NACA and the manufacturer (Bell Aircraft Corporation) suggested three vertical-tail modifications to remedy the deficiencies in the directional characteristics. These modifications included an enlarged vertical tail formed by adding a tip extension to the original vertical tail, a large sharp-edge ventral fin, and a small dorsal fin. The enlarged vertical tail involved only a slight increase in total vertical-tail area from 23.73 to 26.58 square feet but a relatively much larger increase in geometric aspect ratio from 1.24 to 1.73 based on height and area above the horizontal tail. At the request of the Air Material Command, Army Air Forces, flight tests were made to determine the effect of these modifications and of some combinations of these modifications on the directional stability and control characteristics of the airplane, In all, six different vertical-tail. configurations were investigated to determine the lateral and directional oscillation characteristics of the airplane, the sideslip characteristics, the yaw due to ailerons in rudder-fixed rolls from turns and pull-outs, the trim changes due to speed changes; and the trim changes due to power changes. Results of the tests showed that the enlarged vertical tail approximately doubled the directional stability of the airplane and that the pilots considered the directional stability provided by the enlarged vertical tail to be satisfactory. Calculations based on sideslip data obtained at an indicated airspeed of 300 miles per hour showed that the directional stability of the airplane with the original vertical tail corresponded to a value of 0(sub n beta) of -0.00056 whereas for the enlarged vertical tail the estimated va1ue of C(sub n beta) was -0.00130, The ventral fin was found to increase by a moderate amount the directional stability of the airplane with the original vertical tail for smal1 sides1ip angles at low speeds but little consistent change in directional stability was effected by the ventral fin at higher speeds, The effectiveness of the ventral fin was generally much less when used with the enlarged vertical tail than when used with the original vertical tail. The ventral and dorsal fins were found to be very effective in eliminating rudder-force reversals which occurred in low-speed, high-engine-power, sideslipped conditions of flight . Sideslip tests at two altitudes for approximately the sane engine power and indicated airspeed showed that a small decrease in static directional stability occurred with increasing altitude and this decrease in stability was attributed to the increased propeller blade angles required at high altitudes. The variations of rudder pedal force with indicated airspeed using normal rated power and a constant rudder tab setting through the speed range were desirably small for all the configurations tested. The rudder pedal force changed by about 50 pounds for a power change from engine idling power, to normal rated power and this pedal force change was largely independent of airspeed or of vertical-tail configuration for the various configurations tested.

Johnson, Harold I.↗

Designs and Materials for Better Coronagraph Occulting Masks

New designs, and materials appropriate for such designs, are under investigation in an effort to develop coronagraph occulting masks having broad-band spectral characteristics superior to those currently employed. These designs and materials are applicable to all coronagraphs, both ground-based and spaceborne. This effort also offers potential benefits for the development of other optical masks and filters that are required (1) for precisely tailored spatial transmission profiles, (2) to be characterized by optical-density neutrality and phase neutrality (that is, to be characterized by constant optical density and constant phase over broad wavelength ranges), and/or (3) not to exhibit optical- density-dependent phase shifts. The need for this effort arises for the following reasons: Coronagraph occulting masks are required to impose, on beams of light transmitted through them, extremely precise control of amplitude and phase according to carefully designed transmission profiles. In the original application that gave rise to this effort, the concern has been to develop broad-band occulting masks for NASA s Terrestrial Planet Finder coronagraph. Until now, experimental samples of these masks have been made from high-energy-beam-sensitive (HEBS) glass, which becomes locally dark where irradiated with a high-energy electron beam, the amount of darkening depending on the electron-beam energy and dose. Precise mask profiles have been written on HEBS glass blanks by use of electron beams, and the masks have performed satisfactorily in monochromatic light. However, the optical-density and phase profiles of the HEBS masks vary significantly with wavelength; consequently, the HEBS masks perform unsatisfactorily in broad-band light. The key properties of materials to be used in coronagraph occulting masks are their extinction coefficients, their indices of refraction, and the variations of these parameters with wavelength. The effort thus far has included theoretical predictions of performances of masks that would be made from alternative materials chosen because the wavelength dependences of their extinction coefficients and their indices of refraction are such that that the optical-density and phase profiles of masks made from these materials can be expected to vary much less with wavelength than do those of masks made from HEBS glass. The alternative materials considered thus far include some elemental metals such as Pt and Ni, metal alloys such as Inconel, metal nitrides such as TiN, and dielectrics such as SiO2. A mask as now envisioned would include thin metal and dielectric films having stepped or smoothly varying thicknesses (see figure). The thicknesses would be chosen, taking account of the indices of refraction and extinction coefficients, to obtain an acceptably close approximation of the desired spatial transmittance profile with a flat phase profile

Balasubramanian, Kunjithapatham↗

Acceleration Modes and Transitions in Pulsed Plasma Accelerators

Pulsed plasma accelerators typically operate by storing energy in a capacitor bank and then discharging this energy through a gas, ionizing and accelerating it through the Lorentz body force. Two plasma accelerator types employing this general scheme have typically been studied: the gas-fed pulsed plasma thruster and the quasi-steady magnetoplasmadynamic (MPD) accelerator. The gas-fed pulsed plasma accelerator is generally represented as a completely transient device discharging in approximately 1-10 microseconds. When the capacitor bank is discharged through the gas, a current sheet forms at the breech of the thruster and propagates forward under a j (current density) by B (magnetic field) body force, entraining propellant it encounters. This process is sometimes referred to as detonation-mode acceleration because the current sheet representation approximates that of a strong shock propagating through the gas. Acceleration of the initial current sheet ceases when either the current sheet reaches the end of the device and is ejected or when the current in the circuit reverses, striking a new current sheet at the breech and depriving the initial sheet of additional acceleration. In the quasi-steady MPD accelerator, the pulse is lengthened to approximately 1 millisecond or longer and maintained at an approximately constant level during discharge. The time over which the transient phenomena experienced during startup typically occur is short relative to the overall discharge time, which is now long enough for the plasma to assume a relatively steady-state configuration. The ionized gas flows through a stationary current channel in a manner that is sometimes referred to as the deflagration-mode of operation. The plasma experiences electromagnetic acceleration as it flows through the current channel towards the exit of the device. A device that had a short pulse length but appeared to operate in a plasma acceleration regime different from the gas-fed pulsed plasma accelerators was developed by Cheng, et al. The Coaxial High ENerGy (CHENG) thruster operated on the 10-microseconds timescales of pulsed plasma thrusters, but claimed high thrust density, high efficiency and low electrode erosion rates, which are more consistent with the deflagration mode of acceleration. Separate work on gas-fed pulsed plasma thrusters (PPTs) by Ziemer, et al. identified two separate regimes of performance. The regime at higher mass bits (termed Mode I in that work) possessed relatively constant thrust efficiency (ratio of jet kinetic energy to input electrical energy) as a function of mass bit. In the second regime at very low mass bits (termed Mode II), the efficiency increased with decreasing mass bit. Work by Poehlmann et al. and by Sitaraman and Raja sought to understand the performance of the CHENG thruster and the Mode I / Mode II performance in PPTs by modeling the acceleration using the Hugoniot Relation, with the detonation and deflagration modes representing two distinct sets of solutions to the relevant conservation laws. These works studied the proposal that, depending upon the values of the various controllable parameters, the accelerator would operate in either the detonation or deflagration mode. In the present work, we propose a variation on the explanation for the differences in performance between the various pulsed plasma accelerators. Instead of treating the accelerator as if it were only operating in one mode or the other during a pulse, we model the initial stage of the discharge in all cases as an accelerating current sheet (detonation mode). If the current sheet reaches the exit of the accelerator before the discharge is completed, the acceleration mode transitions to the deflagration mode type found in the quasi-steady MPD thrusters. This modeling method is used to demonstrate that standard gas-fed pulsed plasma accelerators, the CHENG thruster, and the quasi-steady MPD accelerator are variations of the same device, with the overall acceleration of the plasma depending upon the behavior of the plasma discharge during initial transient phase and the relative lengths of the detonation and deflagration modes of operation.

Polzin, Kurt A.↗

ISS Ammonia Leak Detection Through X-Ray Fluorescence

Ammonia leaks are a significant concern for the International Space Station (ISS). The ISS has external transport lines that direct liquid ammonia to radiator panels where the ammonia is cooled and then brought back to thermal control units. These transport lines and radiator panels are subject to stress from micrometeorites and temperature variations, and have developed small leaks. The ISS can accommodate these leaks at their present rate, but if the rate increased by a factor of ten, it could potentially deplete the ammonia supply and impact the proper functioning of the ISS thermal control system, causing a serious safety risk. A proposed ISS astrophysics instrument, the Lobster X-Ray Monitor, can be used to detect and localize ISS ammonia leaks. Based on the optical design of the eye of its namesake crustacean, the Lobster detector gives simultaneously large field of view and good position resolution. The leak detection principle is that the nitrogen in the leaking ammonia will be ionized by X-rays from the Sun, and then emit its own characteristic Xray signal. The Lobster instrument, nominally facing zenith for its astrophysics observations, can be periodically pointed towards the ISS radiator panels and some sections of the transport lines to detect and localize the characteristic X-rays from the ammonia leaks. Another possibility is to use the ISS robot arm to grab the Lobster instrument and scan it across the transport lines and radiator panels. In this case the leak detection can be made more sensitive by including a focused 100-microampere electron beam to stimulate X-ray emission from the leaking nitrogen. Laboratory studies have shown that either approach can be used to locate ammonia leaks at the level of 0.1 kg/day, a threshold rate of concern for the ISS. The Lobster instrument uses two main components: (1) a microchannel plate optic (also known as a Lobster optic) that focuses the X-rays and directs them to the focal plane, and (2) a CCD (charge coupled device) focal plane detector that reads out the position and energy of the X-rays, allowing a determination of the leak location. The effective area of the detection system is approximately 2 cm(exp2) at 1 keV. The Lobster astrophysics instrument, designed for monitoring the sky for Xray transients, gives high sensitivity along with large field of view (30×30deg) and good spatial resolution (1 arc min). This offers a significant benefit for detecting ISS ammonia leaks, since the goal is to localize small leaks as efficiently as possible.

Camp, Jordan↗

Investigation of a Laminar Flow Leading Edge

The recent resurgence of interest in utilizing laminar flow on aircraft surfaces for reduction in skin friction drag has generated a considerable amount of research in natural laminar flow (NLF) and hybrid laminar flow control (HLFC) on transonic aircraft wings. This research has focused primarily on airfoil design and understanding transition behavior with little concern for the surface imperfections and manufacturing variations inherent to most production aircraft. In order for laminar flow to find wide-spread use on production aircraft, techniques for constructing the wings must be found such that the large surface imperfections present in the leading edge region of current aircraft do not occur. Toward this end, a modification to existing leading edge construction techniques was devised such that the resulting surface did not contain large gaps and steps as are common on current production aircraft of this class. A lowspeed experiment was first conducted on a simulation of the surface that would result from this construction technique. Preston tube measurements of the boundary layer downstream of the simulated joint and flow visualization using sublimation chemicals validated the literature on the effects of steps on a laminar boundary layer. These results also indicated that the construction technique was indeed compatible with laminar flow. In order to fully validate the compatibility of this construction technique with laminar flow, thus proving that it is possible to build wings that are smooth enough to be used on business jets and light transports in a manner compatible with laminar flow, a flight experiment is being conducted. In this experiment Mach number and Reynolds number will be matched in a real flight environment. The experiment is being conducted using the NASA Dryden F-104 Flight Test Fixture (FTF). The FTF is a low aspect ratio ventral fin mounted beneath an F-104G research aircraft. A new nose shape was designed and constructed for this experiment. This nose shape provides an accelerating pressure gradient in the leading edge region. By flying the aircraft at appropriate Mach numbers and altitudes, this nose shape simulates the leading edge region of a laminar flow wing for a business jet or light transport. Manufactured into the nose shape is a spanwise slot located approximately four inches downstream of the leading edge. The slot, which is an inch wide and one-eighth of an inch deep allows the simulation of surface imperfections, such as gaps and steps at skin joints, which will occur on aircraft using this new construction technique. By placing strips of aluminum of various sizes and shapes in the slot, the effect on the boundary layer of different sizes and shapes of steps and gaps will be examined. It is planned to use five different configurations, differing primarily in the size and number of gaps. Downstream of the slot, the state of the boundary layer is determined using hot film gages and Stanton gages. Agreement between these two very different techniques of measuring boundary layer properties is considered important to being able to state with confidence the effects on the boundary layer of the simulated manufacturing imperfections. To date, the aircraft has not flown. First flights of the aircraft are on schedule to begin October 4, 1993. Low-speed, preliminary experiments at matching Reynolds numbers have been completed.

Drake, Aaron↗

GPS World, Innovation: Autonomous Navigation at High Earth Orbits

Calculating a spacecraft's precise location at high orbital altitudes-22,000 miles (35,800 km) and beyond-is an important and challenging problem. New and exciting opportunities become possible if satellites are able to autonomously determine their own orbits. First, the repetitive task of periodically collecting range measurements from terrestrial antennas to high altitude spacecraft becomes less important-this lessens competition for control facilities and saves money by reducing operational costs. Also, autonomous navigation at high orbital altitudes introduces the possibility of autonomous station keeping. For example, if a geostationary satellite begins to drift outside of its designated slot it can make orbit adjustments without requiring commands from the ground. Finally, precise onboard orbit determination opens the door to satellites flying in formation-an emerging concept for many scientific space applications. The realization of these benefits is not a trivial task. While the navigation signals broadcast by GPS satellites are well suited for orbit and attitude determination at lower altitudes, acquiring and using these signals at geostationary (GEO) and highly elliptical orbits is much more difficult. The light blue trace describes the GPS orbit at approximately 12,550 miles (20,200 km) altitude. GPS satellites were designed to provide navigation signals to terrestrial users-consequently the antenna array points directly toward the earth. GEO and HE0 orbits, however, are well above the operational GPS constellation, making signal reception at these altitudes more challenging. The nominal beamwidth of a Block II/IIA GPS satellite antenna array is approximately 42.6 degrees. At GEO and HE0 altitudes, most of these primary beam transmissions are blocked by the Earth, leaving only a narrow region of nominal signal visibility near opposing limbs of the earth. This region is highlighted in gray. If GPS receivers at GEO and HE0 orbits were designed to use these higher power signals only, precise orbit determination would not be practical. Fortunately, the GPS satellite antenna array also produces side lobe signals at much lower power levels. NASA has designed and tested the Navigator, a new GPS receiver that can acquire and track these weaker signals, thereby dramatically increasing the signal visibility at these altitudes. While using much weaker signals is a fundamental requirement for a high orbital altitude GPS receiver, it is certainly not the only challenge. There are other unique characteristics of this application that must also be considered. For example, Position Dilution of Precision (PDOP) figures are much higher at GEO and HE0 altitudes because visible GPS satellites are concentrated in a much smaller area with respect to the spacecraft antenna. These poor PDOP values contribute considerable error to the point solutions calculated by the spacecraft GPS receiver. Finally, spacecraft GPS receivers must be designed to withstand a variety of extreme environmental conditions. Variations in acceleration between launch and booster separation are extreme. Temperature gradients in the space environment are also severe. Furthermore, radiation effects are a major concern-spacecraft-borne GPS receivers must be designed with radiation-hardened electronics to guard against this phenomenon, otherwise they simply will not work. Perhaps most importantly, there are no opportunities to repair or modify any space-borne GPS receiver after it has been launched. Great care must be taken to ensure all performance characteristics have been analyzed prior to liftoff.

Bamford, William↗

Comparison of removal and spatial mark‐resight models for estimating wild pig density

Density estimation is critical to effectively manage invasive species and elucidate areas of highest concern. For wild pigs (Sus scrofa), the ability to estimate density is complicated because of their variable home range sizes and social structure. Common methods for estimating density (e.g., mark-recapture) may be unsuitable in management applications because additional data needs to be collected before and after management. Removal models offer a suitable alternative to estimate density changes following management and can be applied broadly across areas where management of wild pigs is ongoing. We collected wild pig removal and camera trap data from 25 private properties ranging in size from approximately 0.5 km 2 to 95 km 2 across 3 ecoregions in South Carolina, USA, from 2020–2023. We compared factors affecting consistency and precision of property-level density estimates between removal and spatial mark-resight (SMR) models. In general, excluding 1 large outlier, density estimates from removal models were between 0.60 and 15.85 wild pigs/km 2 (median = 5.34) with a median coefficient of variation (CV) of 0.76 and 95% confidence intervals for the CV between 0.70 and 0.94. Similarly, excluding 1 large outlier, density estimates from SMR were between 0.22 and 30.97 wild pigs/km 2 (median = 5.48) with a median CV of 0.39 and 95% confidence intervals for the CV between 0.38 and 1.20. We found the precision of removal models was affected primarily by the number of wild pigs dispatched in the removal period (3 months) and the ecoregion in which they were removed. None of the covariates, including the number of recaptures (a corresponding measure of sample size), influenced precision of the SMR models, although recaptures did influence the density estimates. At the individual property level, density estimates from our 2 estimators were dissimilar from each other in approximately 80% of instances, although none of the covariates we examined influenced dissimilarity. Our results provide unique insight into how sample size affects density estimates using 2 common methods and into novel SMR models that incorporate both marked and unmarked detections. In addition, the density estimates in this study can be used as a reference for wild pig densities in common land cover types throughout the southeastern United States.

60 APPLIED LIFE SCIENCES↗

Constraints on hydrothermal heat flux through the oceanic lithosphere from global heat flow

A significant discrepancy exists between the heat flow measured at the seafloor and the higher values predicted by thermal models of the cooling lithosphere. This discrepancy is generally interpreted as indicating that the upper oceanic crust is cooled significantly by hydrothermal circulation. The magnitude of this heat flow discrepancy is the primary datum used to estimate the volume of hydrothermal flow, and the variation in the discrepancy with lithospheric age is the primary constraint on how the hydrothermal flux is divided between near-ridge and off-ridge environments. The resulting estimates are important for investigation of both the thermal structure of the lithosphere and the chemistry of the oceans. We reevaluate the magnitude and age variation of the discrepancy using a global heat flow data set substantially larger than in earlier studies, and the GDHI (Global Depth and Heat Flow) model that better predicts the heat flow. We estimate that of the predicted global oceanic heat flux of 32 x 10(exp 12) W, 34% (11 x 10(exp 12) W) occurs by hydrothermal flow. Approximately 30% of the hydrothermal heat flux occurs in crust younger than 1 Ma, so the majority of this flux is off-ridge. These hydrothermal heat flux estimates are upper bounds, because heat flow measurements require sediment at the site and so are made preferentially at topographic lows, where heat flow may be depressed. Because the water temperature for the near-ridge flow exceeds that for the off-ridge flow, the near-ridge water flow will be even a smaller fraction of the total water flow. As a result, in estimating fluxes from geochemical data, use of the high water temperatures appropriate for the ridge axis may significantly overestimate the heat flux for an assumed water flux or underestimate the water flux for an assumed heat flux. Our data also permit improved estimates of the 'sealing' age, defined as the age where the observed heat flow approximately equals that predicted, suggesting that hydrothermal heat transfer has largely ceased. Although earlier studies suggested major differences in sealing ages for different ocean basins, we find that the sealing ages for the Atlantic, Pacific, and Indian oceans are similar and consistent with the sealing age for the entire data set, 65 +/- 10 Ma. The previous inference of a young (approximately 20 Ma) sealing age for the Pacific appears to have biased downward several previous estimates of the global hydrothermal flux. The heat flow data also provide indirect evidence for the mechanism by which the hydrothermal heat flux becomes small, which has often been ascribed to isolation of the igneous crust from seawater due to the hydraulic conductivity of the intervening sediment. We find, however, that even the least sedimented sites show the systematic increase of the ratio of observed to predicted heat flow with age, although the more sedimented sites have a younger sealing age. Moreover, the heat flow discrepancy persists at heavily sedimented sites until approximately 50 Ma. It thus appears that approximately 100-200 m of sediment is neither necessary nor sufficient to stop hydrothermal heat transfer. We therefore conclude that the age of the crust is the primary control on the fraction of heat transported by hydrothermal flow and that sediment thickness has a lesser effect. This inference is consistent with models in which hydrothermal flow decreases with age due to reduced crustal porosity and hence permeability.

Stein, Carol A.↗

Validation of a Manually Oscillating Chair for In-The-Field Assessment of Dynamic Visual Acuity on Crewmembers Within Hours of Returning From Long-Duration Spaceflight

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.

Kreutzberg, G. A.↗