Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “motion correction”

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A stochastic solver based on the residence time algorithm for crystal plasticity models

Abstract The deformation of crystalline materials by dislocation motion takes place in discrete amounts determined by the Burgers vector. Dislocations may move individually or in bundles, potentially giving rise to intermittent slip. This confers plastic deformation with a certain degree of variability that can be interpreted as being caused by stochastic fluctuations in dislocation behavior. However, crystal plasticity (CP) models are almost always formulated in a continuum sense, assuming that fluctuations average out over large material volumes and/or cancel out due to multi-slip contributions. Nevertheless, plastic fluctuations are known to be important in confined volumes at or below the micron scale, at high temperatures, and under low strain rate/stress deformation conditions. Here, we develop a stochastic solver for CP models based on the residence-time algorithm that naturally captures plastic fluctuations by sampling among the set of active slip systems in the crystal. The method solves the evolution equations of explicit CP formulations, which are recast as stochastic ordinary differential equations and integrated discretely in time. The stochastic CP model is numerically stable by design and naturally breaks the symmetry of plastic slip by sampling among the active plastic shear rates with the correct probability. This can lead to phenomena such as intermittent slip or plastic localization without adding external symmetry-breaking operations to the model. The method is applied to body-centered cubic tungsten single crystals under a variety of temperatures, loading orientations, and imposed strain rates.

42 ENGINEERING↗
NASA NTRS2014

OCCAMS: Optically Controlled and Corrected Active Meta-Material Space Structures

Telescopes require sub-wavelength figure (shape) error in order to achieve acceptable Strehl ratios. Traditional methods of achieving this require rigid and therefore heavy mirrors and reaction structures as well as proportionally heavy and expensive spacecraft busses and launch vehicles. This effectively limits the diameter and therefore resolution and collecting area of space optics. Large diameter telescopes must either be heavy or actively controlled. Space telescopes of the size proposed for missions such as TPF (Terrestrial Planet Finder), TPI (Terrestrial Planet Imager), etc. will likely require large active primaries and structures and corrective optics to implement downstream wavefront control. We propose a novel viable enabling technology. High resolution space imaging requires the production of lightweight large aperture optics subject to design tradespace constraints based on nanometer physical tolerances, low system aerial density, high control authority, suitable thermal and mechanical properties, deployment capability, launch vehicle volume constraints, as well as production cost and schedule risk mitigation. Segmented systems impose additional requirements on segment alignment, wavefront phasing, and telescope mass with increased edge diffraction. Photons weigh nothing. Why must even small space telescopes "weigh" tons? Our team is working to demonstrate that they do not, by leveraging 15 years of effort and a novel advanced concept.The ultimate research goal is to refine technology to produce a mission capable ultra-lightweight membrane optic whose shape can be remotely controlled using a laser beam. Using our novel optically-controlled molecular actuators allows substitution of optically-induced control for rigidity and mirror mass. In analogy with noise-canceling headphones, this technique also potentially provides a way to excite vibratory modes to couple with and damp out normal modes of oscillation induced by slewing, other motion and thermal changes, thus allowing a robust response to repointing of a large aperture.

Astronomy↗
NASA NTRS2000

Case Study Analyses of the SUCCESS DC-8 Scanning Lidar Database

Under project SUCCESS (Subsonic Aircraft Contrail and Cloud Effects Special Study) funded by the Atmospheric Effects of Aviation Program, SRI International (SRI) developed an angular scanning backscatter lidar for operation on the NASA DC-8 research aircraft and deployed the scanning lidar during the SUCCESS field campaign. The primary purpose of the lidar was to generate real-time video displays of clouds and contrails above, ahead of, and below the DC-8 as a means to help position the aircraft for optimum cloud and contrail sampling by onboard in situ sensors, and to help extend the geometrical domain of the in situ sampling records. A large, relatively complex lidar database was collected and several data examples were processed to illustrate the value of the lidar data for interpreting the other data records collected during SUCCESS. These data examples were used to develop a journal publication for the special SUCCESS Geophysical Research Letters issue. The data examples justified data analyses of a larger part of the DC-8 lidar database and is the objective of the current study. Efficient processing of the SUCCESS DC-8 scanning lidar database required substantial effort to enhance hardware and software components of the data system that was used for the initial analyses. MATLAB instructions are used to generate altitude and distance color-coded lidar displays corrected for effects introduced by aircraft pitch and forward movement during an angular scan time interval. Onboard in situ sensor atmospheric measurements are propagated to distances ahead of the DC-8 using recorded aircraft velocity so that they can be plotted on the lidar displays for comparison with lidar remotely observed aerosol distributions. Resulting lidar and in situ sensor polar scan displays over extended sampling intervals are integrated into a time series movie format for 36 case studies. Contrails and clouds were detected to ranges of 15 km by the forward-viewing angular scanning lidar and were progressively mapped as the aircraft approached and penetrated them. Near aircraft lidar observations were much better correlated with in situ sensor observations than lidar observations at greater distances ahead of the aircraft. The major cause of this difference was thought to be the about 2 deg. offset of the lidar viewing direction from the flight direction. Contrail spatial distributions were not of the quality obtainable from ground-based lidar observations. This results because contrails tend to become horizontally stratified, vertical distance between angular lidar observations increases with increased distance from the aircraft, and erratic aircraft motions during an angular scan. The most useful lidar observations were made with lidar viewing directions of vertically upward or vertically downward. These provided real-time information on aircraft altitudes to achieve optimum in situ cloud and contrail sampling. At sampling altitudes, the forward viewing angular scanning observations were useful for fine-tuning the aircraft altitude for cloud and contrail penetration. Best information on cloud and contrail properties were obtained from vertically directed lidar observations as the aircraft performed a series of upward and downward penetrations of contrails. This operational mode was especially well suited for lidar and radiometric evaluation of cloud and contrail optical and radiative properties. The vertical viewing lidar detected ice crystals thought to be precipitating from an aircraft contrail and their scavenging by a cirrus cloud layer. The lidar display indicates that the crystals are effective for increasing cirrus cloud density. Vertical angular scanning observations can evaluate the sharp decrease in lidar backscatter for small off-vertical viewing directions that result from horizontally aligned ice crystals and perhaps can provide additional information on crystal shapes. The about 2 deg. offset of the lidar viewing direction from the flight direction is thought to have greatly degraded the forward-viewing angular scanning observations and this mode of operation was not fully evaluated. However, the reasoning for this capability remains valid and the angular scan presentations collected during this program justifies modification of the lidar pod for true forward direction lidar viewing during future cloud and contrail studies.

Uthe, Edward E.↗
NASA NTRS2002

Influence of Turbulence on the Restraint of Liquid Jets by Surface Tension in Microgravity Investigated

Microgravity poses many challenges to the designer of spacecraft tanks. Chief among these are the lack of phase separation and the need to supply vapor-free liquid or liquidfree vapor to the spacecraft processes that require fluid. One of the principal problems of phase separation is the creation of liquid jets. A jet can be created by liquid filling, settling of the fluid to one end of the tank, or even closing a valve to stop the liquid flow. Anyone who has seen a fountain knows that jets occur in normal gravity also. However, in normal gravity, the gravity controls and restricts the jet flow. In microgravity, with gravity largely absent, surface tension forces must contain jets. To model this phenomenon, a numerical method that tracks the fluid motion and the surface tension forces is required. Jacqmin has developed a phase model that converts the discrete surface tension force into a barrier function that peaks at the free surface and decays rapidly away. Previous attempts at this formulation were criticized for smearing the interface. This can be overcome by sharpening the phase function, double gridding the fluid function, and using a higher order solution for the fluid function. The solution of this equation can be rewritten as two coupled Poisson equations that also include the velocity. After the code was implemented in axisymmetric form and verified by several test cases at the NASA Glenn Research Center, the drop tower runs of Aydelott were modeled. Work last year with a laminar model was found to overpredict Aydelott's results, except at the lowest Reynolds number conditions of 400. This year, a simple turbulence model was implemented by adding a turbulent viscosity based on the mixing-length hypothesis and empirical measurements of previous works. Predictions made after this change was implemented have been much closer to experimentally observed flow patterns and geyser heights. Two model runs is shown. The first, without any turbulence correction, breaks through the free surface and strikes the far end of the tank. In the second, the turbulence spreads the jet momentum over more of the free surface, enabling the surface tension forces to turn the jet back into the bulk liquid. The model geyser height with the second model is 1.1 cm. This is quite close to the 1.5-cm geyser height measured by Aydelott.

Chato, David J.↗
NASA NTRS2010

Physics of Life: A Model for Non-Newtonian Properties of Living Systems

This innovation proposes the reconciliation of the evolution of life with the second law of thermodynamics via the introduction of the First Principle for modeling behavior of living systems. The structure of the model is quantum-inspired: it acquires the topology of the Madelung equation in which the quantum potential is replaced with the information potential. As a result, the model captures the most fundamental property of life: the progressive evolution; i.e. the ability to evolve from disorder to order without any external interference. The mathematical structure of the model can be obtained from the Newtonian equations of motion (representing the motor dynamics) coupled with the corresponding Liouville equation (representing the mental dynamics) via information forces. All these specific non-Newtonian properties equip the model with the levels of complexity that matches the complexity of life, and that makes the model applicable for description of behaviors of ecological, social, and economical systems. Rather than addressing the six aspects of life (organization, metabolism, growth, adaptation, response to stimuli, and reproduction), this work focuses only on biosignature ; i.e. the mechanical invariants of life, and in particular, the geometry and kinematics of behavior of living things. Living things obey the First Principles of Newtonian mechanics. One main objective of this model is to extend the First Principles of classical physics to include phenomenological behavior on living systems; to develop a new mathematical formalism within the framework of classical dynamics that would allow one to capture the specific properties of natural or artificial living systems such as formation of the collective mind based upon abstract images of the selves and non-selves; exploitation of this collective mind for communications and predictions of future expected characteristics of evolution; and for making decisions and implementing the corresponding corrections if the expected scenario is different from the originally planned one. This approach postulates that even a primitive living species possesses additional, non-Newtonian properties that are not included in the laws of Newtonian or statistical mechanics. These properties follow from a privileged ability of living systems to possess a self-image (a concept introduced in psychology) and to interact with it. The proposed mathematical system is based on the coupling of the classical dynamical system representing the motor dynamics with the corresponding Liouville equation describing the evolution of initial uncertainties in terms of the probability density and representing the mental dynamics. The coupling is implemented by the information-based supervising forces that can be associated with self-awareness. These forces fundamentally change the pattern of the probability evolution, and therefore, lead to a major departure of the behavior of living systems from the patterns of both Newtonian and statistical mechanics. This innovation is meant to capture the signature of life based only on observable behavior, not on any biochemistry. This will not prevent the use of this model for developing artificial living systems, as well as for studying some general properties of behavior of natural, living systems.

Zak, Michail↗
NASA NTRS1994

Steps toward determination of the size and structure of the broad-line region in active galactic nuclei. 8: An intensive HST, IUE, and ground-based study of NGC 5548

We present the data and initial results from a combined HST/IUE/ground-based spectroscopic monitoring campaign on the Seyfert 1 galaxy NGC 5548 that was undertaken in order to address questions that require both higher temporal resolution and higher signal-to-noise ratios than were obtained in our previous multiwavelength monitoring of this galaxy in 1988-89. IUE spectra were obtained once every two days for a period of 74 days beginning on 14 March 1993. During the last 39 days of this campaign, spectroscopic observations were also made with the HST Faint Object Spectrograph (FOS) on a daily basis. Ground-based observations, consisting of 165 optical spectra and 77 photometric observations (both CCD imaging and aperture photometry), are reported for the period 1992 October to 1993 September, although much of the data are concentrated around the time of the satellite-based program. These data constitute a fifth year of intensive optical monitoring of this galaxy. In this contribution, we describe the acquisition and reduction of all of the satellite and ground-based data obtained in this program. We describe in detail various photometric problems with the FOS and explain how we identified and corrected for various anomalies. During the HST portion of the monitoring campaign, the 1350 A continuum flux is found to have varied by nearly a factor of two. In other wavebands, the continuum shows nearly identical behavior, except that the amplitude of variability is larger at shorter wavelengths, and the continuum light curves appear to show more short time-scale variability at shorter wavelengths. The broad emission lines also vary in flux, with amplitudes that are slightly smaller than the UV continuum variations and with a small time delay relative to the UV continuum. On the basis of simple time-series analysis of the UV and optical continuum and emission line light curves, we find (1) that the ultraviolet and optical continuum variations are virtually simultaneous, with any lag between the 1350 A continuum and the 5100 A continuum amounting to less than about one day; (2) that the variations in the highest ionization lines observed, He II lambda 1640 and N V lambda 1240, lag behind the continuum variations by somewhat less than 2 days, and (3) that the velocity field of the C IV-emitting region is not dominated by radial motion. The results on the C IV velocity field are preliminary and quite uncertain, but there are some weak indications that the emission-line (wings absolute value of Delta upsilon is greater than or equal to 3000 km/s) respond to continuum variations slightly more rapidly than does the core. The optical observations show that the variations in the broad H beta line flux follow the continuum variations with a time lag of around two weeks, about twice the lag for Ly alpha and C IV, as in our previous monitoring campaign on this same galaxy. However, the lags measured for Ly alpha, C IV, and H Beta are each slightly smaller than previously determined. We confirm two trends reported earlier, namely (1) that the UV/optical continuum becomes 'harder' as it gets brighter, and (2) that the highest ionization emission lines have the shortest lags, thus indicating radial ionization stratificatin of a broad-line region that spans over an order of magnitude range in radius.

Korista, K.↗
DOE OSTI2020

Stability Analysis of Inline ZFP Compression for Floating-Point Data in Iterative Methods

Currently, the dominating constraint in many high performance computing applications is data capacity and bandwidth, in both internode communications and even moreso in intranode data motion. A new approach to address this limitation is to make use of data compression in the form of a compressed data array. Storing data in a compressed data array and converting to standard IEEE-754 types as needed during a computation can reduce the pressure on bandwidth and storage. However, repeated conversions (lossy compression and decompression) introduce additional approximation errors, which need to be shown to not significantly affect the simulation results. Here, we extend recent work that analyzed the error of a single use of compression and decompression of the ZFP compressed data array representation to the case of time-stepping and iterative schemes, where an advancement operator is repeatedly applied in addition to the conversions. We show that the accumulated error for iterative methods involving fixed-point and time evolving iterations is bounded under standard constraints. An upper bound is established on the number of additional iterations required for the convergence of stationary fixed-point iterations. An additional analysis of traditional forward and backward error of stationary iterative methods using ZFP compressed arrays is also presented. The results of several 1D, 2D, and 3D test problems are provided to demonstrate the correctness of the theoretical bounds.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1992

Pressurized Lunar Rover

The pressurized lunar rover (PLR) consists of a 7 m long, 3 m diameter cylindrical main vehicle and a trailer which houses the power and heat rejection systems. The main vehicle carries the astronauts, life support systems, navigation and communication systems, directional lighting, cameras, and equipment for exploratory experiments. The PLR shell is constructed of a layered carbon-fiber/foam composite. The rover has six 1.5 m diameter wheels on the main body and two 1.5 m diameter wheels on the trailer. The wheels are constructed of composites and flex to increase traction and shock absorption. The wheels are each attached to a double A-arm aluminum suspension, which allows each wheel 1 m of vertical motion. In conjunction with a 0.75 m ground clearance, the suspension aids the rover in negotiating the uneven lunar terrain. The 15 N-m torque brushless electric motors are mounted with harmonic drive units inside each of the wheels. The rover is steered by electrically varying the speeds of the wheels on either side of the rover. The PLR trailer contains a radiosotope thermoelectric generator providing 6.7 kW. A secondary back-up energy storage system for short-term high-power needs is provided by a bank of batteries. The trailer can be detached to facilitate docking of the main body with the lunar base via an airlock located in the rear of the PLR. The airlock is also used for EVA operation during missions. Life support is a partly regenerative system with air and hygiene water being recycled. A layer of water inside the composite shell surrounds the command center. The water absorbs any damaging radiation, allowing the command center to be used as a safe haven during solar flares. Guidance, navigation, and control are supplied by a strapdown inertial measurement unit that works with the on-board computer. Star mappers provide periodic error correction. The PLR is capable of voice, video, and data transmission. It is equipped with two 5 W X-band transponder, allowing simultaneous transmission and reception. An S-band transponder is used to communicate with the crew during EVA. The PLR has a total mass of 6197 kg. It has a nominal speed of 10 km/hr and a top speed of 18 km/hr. The rover is capable of towing 3 metric tons (in addition to the RTG trailer).

Creel, Kenneth↗
NASA NTRS2001

The Transition from Hydrogen Bonding to Ionization in (HCI)n(NH3)n and (HCI)n(H2O)n Clusters: Consequences for Anharmonic Vibrational Spectroscopy

Anharmonic vibrational frequencies and intensities are calculated for 1:1 and 2:2 (HCl)(sub n)(NH3)(sub n) and (HCl)(sub n)(H2O)(sub n) complexes, employing the correlation-corrected vibrational self-consistent field method with ab initio potential surfaces at the MP2/TZP computational level. In this method, the anharmonic coupling between all vibrational modes is included, which is found to be important for the systems studied. For the 4:4 (HCl)(sub n)(H2O)(sub n) complex, the vibrational spectra are calculated at the harmonic level, and anharmonic effects are estimated. Just as the (HCl)(sub n)(NH3)(sub n) structure switches from hydrogen-bonded to ionic for n=2, the (HCl)(sub n)(H2O)(sub n) switches to ionic structure for n=4. For (HCl)2(H2O)2, the lowest energy structure corresponds to the hydrogen-bonded form. However, configurations of the ionic form are separated from this minimum by a barrier of less than an O-H stretching quantum. This suggests the possibility of experiments on ionization dynamics using infrared excitation of the hydrogen-bonded form. The strong cooperative effects on the hydrogen bonding, and concomitant transition to ionic bonding, makes an accurate estimate of the large anharmonicity crucial for understanding the infrared spectra of these systems. The anharmonicity is typically of the order of several hundred wave numbers for the proton stretching motions involved in hydrogen or ionic bonding, and can also be quite large for the intramolecular modes. In addition, the large cooperative effects in the 2:2 and higher order (HCl(sub n)(H2O)(sub n) complexes may have interesting implications for solvation of hydrogen halides at ice surfaces.

Chaban, Galina M.↗
DOE OSTI2023

Computer-controlled finishing via dynamically constraint position-velocity-time scheduler

In a Computer Numerical Controlled (CNC) finishing process, the target material removal from an optical surface is guided by the convolution between the influence function of a machine tool and its dwell time at certain points over the surface. To reduce dynamics stressing and increase machining efficiency, the dwell time must be converted to varying velocities, which are the actual inputs to the machine tool controller. Conventionally, the conversion assumed constant acceleration and relied on linear motion interpolation, which caused discontinuities in velocities. This unsmooth motion affects the material removal distribution, and, thus, the accuracy of the finished surface shape. Many modern CNC machines support the smoother, cubic-polynomial interpolated Position-Velocity-Time (PVT) motion mode; however, the conventional scheduler may fail to provide suitable velocities for the PVT. Here in this study answers this challenge by proposing a novel PVT-based velocity scheduler that achieves smooth motion while considering CNC dynamic limits. Firstly, the principle of the PVT is explained, and the PVT-based velocity scheduler is formulated. Secondly, a quadratic programming is used to optimize the velocities by imposing the CNC dynamic constraints and the C 1 continuities (zeroth and first derivatives are continuous) simultaneously. Thirdly, the smoothness and accuracy of the scheduled velocities are studied on different kinds of tool paths via simulation. Finally, a sub-0.3 nm level surface finishing experiment using ion beam figuring is demonstrated to verify the feasibility of the proposed method. The PVT-based scheduler and simulator code is open-sourced.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Bulk Spray and Individual Plume Characterization of LPG and Iso-Octane Sprays at Engine-Like Conditions

This study presents experimental and numerical examination of directly injected (DI) propane and iso-octane, surrogates for liquified petroleum gas (LPG) and gasoline, respectively, at various engine like conditions with the overall objective to establish the baseline with regards to fuel delivery required for future high efficiency DI-LPG fueled heavy-duty engines. Sprays for both iso-octane and propane were characterized and the results from the optical diagnostic techniques including high-speed Schlieren and planar Mie scattering imaging were applied to differentiate the liquid-phase regions and the bulk spray phenomenon from single plume behaviors. The experimental results, coupled with high-fidelity internal nozzle-flow simulations were then used to define best practices in CFD Lagrangian spray models. Optical imaging revealed that unlike iso-octane, propane’s spray propagation was fed by its flash boiling, spray collapse, and high degree of vaporization, resulting in a direct proportionality of propane’s penetration length to temperature. These unique features of propane and its variation from iso-octane’s spray pattern, contributed to its classification as an unconventional spray. Appropriate corrections to the injection and breakup models were developed to reproduce the under-expanded jet dynamics and to mimic the flash boiling-driven spray development observed with propane sprays. The simulation results were found to be sensitive to cone and inclusion angles of the blob injector. The current work represents a first assessment of the capability of the commonly available models for engine-spray simulations and highlights the fact that, despite the reasonable agreement obtained in the fuel vapor morphology, the representation of the liquid phase lacks accuracy and requires further model development.

02 PETROLEUM↗
DOE OSTI2021

Toward a First-Principles Framework for Predicting Collective Properties of Electrolytes

Conspectus Given the universal importance of electrolyte solutions, it is natural to expect that we have a nearly complete understanding of the fundamental properties of these solutions (e.g., the chemical potential) and that we can therefore explain, predict, and control the phenomena occurring in them. In fact, reality falls short of these expectations. But, recent advances in the simulation and modeling of electrolyte solutions indicate that it should soon be possible to make progress toward these goals. In this Account, we will discuss the use of first-principles interaction potentials based in quantum mechanics (QM) to enhance our understanding of electrolyte solutions. Specifically, we will focus on the use of quantum density functional theory (DFT) combined with molecular dynamics simulation (DFT-MD) as the foundation for our approach. The overarching concept is to understand and accurately reproduce the balance between local or short-ranged (SR) structural details and long-range (LR) correlations, allowing the prediction of the thermodynamics of both single ions in solution as well as the collective interactions characterized by activity/osmotic coefficients. In doing so, relevant collective motions and driving forces characterized by chemical potentials can be determined. Here, in this Account, we will make the case that understanding electrolyte solutions requires a faithful QM representation of the SR nature of the ion–ion, ion–water, and water–water interactions. However, the number of molecules that is required for collective behavior makes the direct application of high-level QM methods that contain the best SR physics untenable, making methods that balance accuracy and efficiency a practical goal. Alternatives such as continuum solvent models (CSMs) and empirically based classical molecular dynamics have been extensively employed to resolve this problem but without yet overcoming the fundamental issue of SR accuracy. We will demonstrate that accurately describing the SR interaction is imperative for predicting both intrinsic properties, namely, at infinite dilution, and collective properties of electrolyte solutions. DFT has played an important role in our understanding of condensed phase systems, e.g., bulk liquid water, the air–water interface, ions in bulk, and at the air–water interface. This approach holds huge promise to provide benchmark calculations of electrolyte solution properties that will allow for the development and improvement of more efficient methods, as well as an enhanced understanding of fundamental phenomena. However, the standard protocol using the generalized gradient approximation with van der Waals (vdW) correction requires improvement in order to achieve a high level of quantitative accuracy. Simply simulating with higher level DFT functionals may not be the best route considering the significant computational cost. Alternative methods of incorporating information from higher levels of QM should be explored; e.g., using force matching techniques on small clusters, where high level benchmark calculations are possible, to develop ideal correction terms to the DFT functional is a promising possibility. We argue that DFT with statistical mechanics is becoming an increasingly useful framework enabling the prediction of collective electrolyte properties.

Duignan, Timothy T.↗
NASA NTRS2015

Psychophysical Calibration of Mobile Touch-Screens for Vision Testing in the Field

The now ubiquitous nature of touch-screen displays in cell phones and tablet computers makes them an attractive option for vision testing outside of the laboratory or clinic. Accurate measurement of parameters such as contrast sensitivity, however, requires precise control of absolute and relative screen luminances. The nonlinearity of the display response (gamma) can be measured or checked using a minimum motion technique similar to that developed by Anstis and Cavanagh (1983) for the determination of isoluminance. While the relative luminances of the color primaries vary between subjects (due to factors such as individual differences in pre-retinal pigment densities), the gamma nonlinearity can be checked in the lab using a photometer. Here we compare results obtained using the psychophysical method with physical measurements for a number of different devices. In addition, we present a novel physical method using the device's built-in front-facing camera in conjunction with a mirror to jointly calibrate the camera and display. A high degree of consistency between devices is found, but some departures from ideal performance are observed. In spite of this, the effects of calibration errors and display artifacts on estimates of contrast sensitivity are found to be small.

gamma correction↗
NASA NTRS2000

Computations of Boiling in Microgravity

The absence (or reduction) of gravity, can lead to major changes in boiling heat transfer. On Earth, convection has a major effect on the heat distribution ahead of an evaporation front, and buoyancy determines the motion of the growing bubbles. In microgravity, convection and buoyancy are absent or greatly reduced and the dynamics of the growing vapor bubbles can change in a fundamental way. In particular, the lack of redistribution of heat can lead to a large superheat and explosive growth of bubbles once they form. While considerable efforts have been devoted to examining boiling experimentally, including the effect of microgravity, theoretical and computational work have been limited. Here, the growth of boiling bubbles is studied by direct numerical simulations where the flow field is fully resolved and the effects of inertia, viscosity, surface deformation, heat conduction and convection, as well as the phase change, are fully accounted for. Boiling involves both fluid flow and heat transfer and thus requires the solution of the Navier-Stokes and the energy equations. The numerical method is based on writing one set of governing transport equations which is valid in both the liquid and vapor phases. This local, single-field formulation incorporates the effect of the interface in the governing equations as source terms acting only at the interface. These sources account for surface tension and latent heat in the equations for conservation of momentum and energy as well as mass transfer across the interface due to phase change. The single-field formulation naturally incorporates the correct mass, momentum and energy balances across the interface. Integration of the conservation equations across the interface directly yields the jump conditions derived in the local instant formulation for two-phase systems. In the numerical implementation, the conservation equations for the whole computational domain (both vapor and liquid) are solved using a stationary grid and the phase boundary is followed by a moving unstructured two-dimensional grid. While two-dimensional simulations have been used for preliminary studies and to examine the resolution requirement, the focus is on fully three-dimensional simulations. The numerical methodology, including the parallelization and grid refinement strategy is discussed, and preliminary results shown. For buoyancy driven flow, the heat transfer is in good agreement with experimental correlations. The changes when gravity is turned off and/or fluid shear is added are discussed, as well as the difference between simulations of a layer freely releasing bubbles versus simulations using only one wavelength initial perturbation. Figure 1 shows the early stages of the formation of a three-dimensional bubble from a thin vapor layer. The boundary conditions are periodic in the x and y direction, the bottom is a hot and the top allows a free outflow. The jagged edge of the surface close to the bottom of the computational domain is due to some of the surface elements being on the other side of the domain and some elements not plotted by our plotting routine. In the second figure, we show the temperature distribution through two perpendicular planes.

Tryggvason, G.↗
NASA NTRS1997

Interaction of Burning Metal Particles

Multiple particle/droplet flames are ubiquitous in practical combustion systems, and thus the flame interaction processes are of great practical importance. This explains the strong current interest in interactive combustion phenomena. This research is aimed at the investigation of combustion parameters of microgravity model aerosols: relatively large uniform metal particles aerosolized in microgravity environment. An experiment consisting of creation and ignition of a metal multiparticle system in microgravity and high-speed video-recording of the combustion events will produce visual records of the development of individual particle flames, their interactions and the particle motion they induce simultaneously with the observation of the entire aerosol combustion process. Frame by frame analysis of the video-images taken using a high-speed movie camera will allow one to determine particle brightness temperatures and the decrease in particle diameter during combustion. Analysis of the experimental results and comparison with the results of single metal particle combustion experiments, conducted under similar microgravity conditions in the framework of a parallel program, will establish the relationship between single and multiple particle burning rates and combustion temperatures, concentrations at which the flame substructure forms rather than individual particle flames, efficiency of radiative heat transfer in metal aerosol combustion, what is the role of electrostatic forces in structuring the flame and the effect of that structure on the flame propagation rate. Although some details of fine particle aerosol clouds, such as the kinetics limited burning rate, radiative heat transfer in a system with a high specific surface, particle induced turbulence, etc., will probably not be very well simulated in the planned experiments, they are relatively well understood and can be accounted for using an adequate individual particle combustion model. On the other hand, the unknown particle interaction phenomena to be studied in this research program are generic and can be readily projected onto a cloud combustion model which would also incorporate a correction for particle size.

Dreizin, Edward L.↗
NASA NTRS1961

A Simulator Study of the Effectiveness of a Pilot's Indicator which Combined Angle of Attack and Rate of Change of Total Pressure as Applied to the Take-Off Rotation and Climbout of a Supersonic Transport

A simulator study has been made to determine the effectiveness of a single instrument presentation as an aid to the pilot in controlling both rotation and climbout path in take-off. The instrument was basically an angle-of-attack indicator, biased with a total-pressure-rate input as a means of suppressing the phugoid oscillation. Linearized six-degree-of-freedom equations of motion were utilized in simulating a hypothetical supersonic transport as the test vehicle. Each of several experienced pilots performed a number of simulated take-offs, using conventional flight instruments and either an angle-of-attack instrument or the combined angle-of-attack and total-pressure-rate instrument. The pilots were able to rotate the airplane, with satisfactory precision, to the 15 deg. angle of attack required for lift-off when using either an angle-of-attack instrument or the instrument which combined total-pressure-rate with angle of attack. At least 4 to 6 second-S appeared to be required for rotation to prevent overshoot, particularly with the latter instrument. The flight paths resulting from take-offs with simulated engine failures were relatively smooth and repeatable within a reasonably narrow band when the combined angle-of-attack and total-pressure-rate instrument presentation was used. Some of the flight paths resulting from take-offs with the same engine-failure conditions were very oscillatory when conventional instruments and an angle-of-attack instrument were used. The pilots considered the combined angle-of-attack and total- pressure-rate instrument a very effective aid. Even though they could, with sufficient practice, perform satisfactory climbouts after simulated engine failure by monitoring the conventional instruments and making correction based on their readings, it was much easier to maintain a smooth flight path with the single combined angle-of-attack and total-pressure-rate instrument.

Hall, Albert W.↗
DOE OSTI2023

Experiment and simulation of high-speed gas jet penetration into a semicircular fluidized bed

This work marks the third in a series of experiments that were in a semi-circular, gas-fluidized bed with side jets. In this work, the particles are (nominally) 1 mm ceramic beads. The bed is operated just at and slightly above and below the minimum fluidization velocity and additional fluidization is provided by two high-speed gas located on the sides of the bed near the flat, front face of the unit. Two primary measurements are taken: high-speed video recording of the front of the bed and bed pressure drop from a tap in the back of the bed. Particle Image Velocimetry (PIV) is used to determine particle motion, characterized as a mean Froude number, from the high-speed video. A CFD-DEM model of the bed is presented using the recently released MFIX-Exa code. Four model subvariants are considered using two methods of representing the jets and two drag models, both of which are calibrated to exactly match the experimentally measured minimum fluidization velocity. Although it is more difficult to determine the jet penetration depths in a straightforward manner as in the previous works using Froude number contours, the CFD-DEM results compare quite well to the PIV measurements, particularly for submodel flow Syam. Unfortunately, the good agreement of the solids-phase is overshadowed by significant disagreement in the gas-phase data. Specifically, the predicted time averaged standard deviation of the pressure drop is found to be over an order of magnitude larger than measured. Due to the low value of the measurements, just 1% of the mean bed pressure drop, it seems possible that the data is in error. On the other hand, the model may not be accurately capturing pressure attenuation through an under-fluidized region in the back of the bed. Without the possibility additional experiments to test the validity of the data, this work is simply being reported “as is” without being able to indicate which, either the simulation or the experiment, is more correct.<br>

Fullmer, William D.↗
NASA NTRS2003

System Identification of X-33 Neural Network

Modern flight control research has improved spacecraft survivability as its goal. To this end we need to have a failure detection system on board. In case the spacecraft is performing imperfectly, reconfiguration of control is needed. For that purpose we need to have parameter identification of spacecraft dynamics. Parameter identification of a system is called system identification. We treat the system as a black box which receives some inputs that lead to some outputs. The question is: what kind of parameters for a particular black box can correlate the observed inputs and outputs? Can these parameters help us to predict the outputs for a new given set of inputs? This is the basic problem of system identification. The X33 was supposed to have the onboard capability of evaluating the current performance and if needed to take the corrective measures to adapt to desired performance. The X33 is comprised of both rocket and aircraft vehicle design characteristics and requires, in general, analytical methods for evaluating its flight performance. Its flight consists of four phases: ascent, transition, entry and TAEM (Terminal Area Energy Management). It spends about 200 seconds in ascent phase, reaching an altitude of about 180,000 feet and a speed of about 10 to 15 Mach. During the transition phase which lasts only about 30 seconds, its altitude may increase to about 190,000 feet but its speed is reduced to about 9 Mach. At the beginning of this phase, the Main Engine is Cut Off (MECO) and the control is reconfigured with the help of aerosurfaces (four elevons, two flaps and two rudders) and reaction control system (RCS). The entry phase brings down the altitude of X33 to about 90,000 feet and its speed to about Mach 3. It spends about 250 seconds in this phase. Main engine is still cut off and the vehicle is controlled by complex maneuvers of aerosurfaces. The last phase TAEM lasts for about 450 seconds and the altitude and speed, both are reduced to zero. The present attempt, as a start, focuses only on the entry phase. Since the main engine remains cut off in this phase, there is no thrust acting on the system. This considerably simplifies the equations of motion. We introduce another simplification by assuming the system to be linear after some non-linearities are removed analytically from our consideration. Under these assumptions, the problem could be solved by Classical Statistics by employing the least sum of squares approach. Instead we chose to use the Neural Network method. This method has many advantages. It is modern, more efficient, can be adapted to work even when the assumptions are diluted. In fact, Neural Networks try to model the human brain and are capable of pattern recognition.

Aggarwal, Shiv↗