Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “conditionally specified models”

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 379 records · Page 21

Lander Lighting Solution to Reduce Pilot & Autonomous Approach Errors

The south pole lighting environment will have harsh low inclination sunlight, making overhead judgement of surface features difficult. Autonomous solutions are great, however, the need for visual monitoring and independent go/no-go decisions remain. Our project proposes that lunar landing systems will be better served by including a powerful landing light system that improves visibility of surfaces from overhead by illuminating the ground at critical distances for the crew to make last minute decisions regarding an approach. The project utilized computer-based optical modeling software to predict requirements for a potential landing light system. The analysis based the lamp prediction from commercially available LED chip sets and lamp optics. The goal was to illustrate a method to raise the surface contrast of a landing site within an acceptable contrast threshold for most camera systems and human observers to recognize hazards that would not be noticed with low inclination sunlight alone. The Apollo lunar landings benefitted from overhead sun or dark conditions. The surface lighting at the Lunar South Pole is a harsh environment where surfaces are lit from a low inclination angle by the sun (from the side). This change in lighting condition precipitates a need for updated lunar landing systems that facilitate improved recognition of landing sites, and thereby increase pilot awareness of landing hazards. The reliance on LIDAR and other autonomous mechanisms alone is risky given the known usage of visual monitoring for operator concurrence on current spacecraft programs and present-day autonomous land-based vehicles. Visual monitoring via cameras or windows requires the surface contrast to be within 3 orders of magnitude for reasonable recognition of objects. Artificial overhead illumination, when sufficiently sized, provides a means to even out contrast problems created by low inclination sunlight, potentially reducing piloting errors. Current vehicle requirements do not specify this type of guidance for the purpose of increasing mission success. An optical ray-trace simulation model was developed in Zemax Optics Studio to predict the best combination of LED power, LED optics, lamp quantity, and lamp location to raise the surface contrast to within 2 orders of magnitude from 3 orders required to further increased visibility and reduce risk. The project considered the following design constraints: potential base diameter of lander, approach distance(s) for a go-no-go decision point (200 meter), solar inclination angle (2-7), lunar surface reflectance (10%), LED chip sets, LED focusing optics, LED power, lamp quantity, lamp locations, and illumination diameter of lunar surface landing zone. The results can be used to establish minimum design constraints for vehicle landing light systems. With a solar inclination angle ranging from 2-7 degrees, the horizontal illumination of the lunar surface is attenuated by about 10% when compared to overhead illumination from the Sun. This modifies the sun's maximum of 130,000 lux to 13,000 lux horizontal illuminance. The artificial lighting system was designed to provide an 18-meter-wide illumination zone, to create viewing clearances around a 6-meter-wide lander. The system provides an average illuminance of 300 lux, meeting the 2 orders of magnitude criteria. The solution utilized modern Chip On Board LEDs, that each utilized 17 watts, with focusing Total Internal Reflection (TIR) lenses. A lighting system of 300 LEDs was arrayed along the "bottom" of a “lander”. With 17 watts per LED, the system is estimated to require 5100 watts. This is a large amount of power, but it would only be needed during critical phases during the landing. LED lighting systems can be dimmed, and it is assumed that as the lander arrives closer to the landing site, the lighting system power can be adjusted as needed to produce the necessary surface illuminance. The designed reduction of contrast improves reliability of safety assessments using real time visible light camera systems and out the window viewing by the crew.

T A Clark↗

Manifold traversing as a model for learning control of autonomous robots

This paper describes a recipe for the construction of control systems that support complex machines such as multi-limbed/multi-fingered robots. The robot has to execute a task under varying environmental conditions and it has to react reasonably when previously unknown conditions are encountered. Its behavior should be learned and/or trained as opposed to being programmed. The paper describes one possible method for organizing the data that the robot has learned by various means. This framework can accept useful operator input even if it does not fully specify what to do, and can combine knowledge from autonomous, operator assisted and programmed experiences.

Szakaly, Zoltan F.↗

A hydromagnetic model of corotating conductive solar wind streams

Under the assumption of quasi-azimuthal symmetry the governing equations of a steady hydromagnetic flow in a thermally conductive flux tube possess six invariants. Four of them represent constancy of mass efflux, energy efflux, angular momentum efflux and magnetic flux. Based on the entropy equation we obtain useful approximation in explicit expressions for the two remaining invariants. One of them provides the constraint which determines the compatible heat flux to ensure a vanishing pressure at infinity. Thus, the admissible solution that represents a corotating solar wind stream in terms of specified interplanetary condition can be calculated by an algebraic method, without the necessity of numerical integration. A two-point relationship is then derived, which correlates the solar wind properties at two separated interplanetary sites measured at two properly separated instants. This relationship may be applied to observational data from spacecraft and earth-bound satellites to discern the corotation feature in the solar wind.

Yeh, T.↗

Numerical simulation of thermally induced near-surface flows over Martian terrain

Numerical simulations of the Martian near-surface wind regime using a mesoscale atmospheric model have shown that the thermally induced near-surface winds are analogous to terrestrial circulations. In particular, katabatic wind displays a striking similarity to flow observed over Antarctica. Introduction of solar radiation strongly perturbs the slope flows; anabatic conditions develop in middle to high latitudes during the daytime hours due to the solar heating of the sloping terrain. There appears to be a rapid transition from the katabatic to the anabatic flow regimes, emphasizing the primary importance of radiative exchanges at the surface in specifying the horizontal pressure gradient force.

Parish, T. R.↗

Impact of Variations on 1-D Flow in Gas Turbine Engines via Monte Carlo Simulations

The unsteady compressible inviscid flow is characterized by the conservations of mass, momentum, and energy; or simply the Euler equations. In this paper, a study of the subsonic one-dimensional Euler equations with local preconditioning is presented using a modal analysis approach. Specifically, this study investigates the behavior of airflow in a gas turbine engine using the specified conditions at the inflow and outflow boundaries of the compressor, combustion chamber, and turbine, to determine the impact of variations in pressure, velocity, temperature, and density at low Mach numbers. Two main questions motivate this research: 1) Is there any aerodynamic problem with the existing gas turbine engines that could impact aircraft performance? 2) If yes, what aspect of a gas turbine engine could be improved via design to alleviate that impact and to optimize aircraft performance? This paper presents an initial attempt to model the flow behavior in terms of their eigenfrequencies subject to the assumption of the uncertainty or variation (perturbation). The flow behavior is explored using simulation outputs from a customer-deck model obtained from Pratt & Whitney. Variations of the main variables (i.e., pressure, temperature, velocity, density) about their mean states at the inflow and outflow boundaries of the compressor, combustion chamber, and turbine are modeled. Flow behavior is analyzed for the high-pressure compressor and combustion chamber utilizing the conditions on their left and right boundaries. In the same fashion, similar analyses are carried out for the high-pressure and low-pressure turbines. In each case, the eigenfrequencies that are obtained for different boundary conditions are examined closely based on their probabilistic distributions, a result of a Monte Carlo 10,000 sample simulation. Furthermore, the characteristic waves and wave response are analyzed and contrasted among different cases, with and without preconditioners. The results reveal the existence of flow instabilities due to the combined effect of variations and excessive pressures in the case of the combustion chamber and high-pressure turbine. Finally, a discussion is presented on potential impacts of the instabilities and what can be improved via design to alleviate them for a better aircraft performance.

Ngo, Khiem Viet↗

NDARC - NASA Design and Analysis of Rotorcraft Validation and Demonstration

Validation and demonstration results from the development of the conceptual design tool NDARC (NASA Design and Analysis of Rotorcraft) are presented. The principal tasks of NDARC are to design a rotorcraft to satisfy specified design conditions and missions, and then analyze the performance of the aircraft for a set of off-design missions and point operating conditions. The aircraft chosen as NDARC development test cases are the UH-60A single main-rotor and tail-rotor helicopter, the CH-47D tandem helicopter, the XH-59A coaxial lift-offset helicopter, and the XV-15 tiltrotor. These aircraft were selected because flight performance data, a weight statement, detailed geometry information, and a correlated comprehensive analysis model are available for each. Validation consists of developing the NDARC models for these aircraft by using geometry and weight information, airframe wind tunnel test data, engine decks, rotor performance tests, and comprehensive analysis results; and then comparing the NDARC results for aircraft and component performance with flight test data. Based on the calibrated models, the capability of the code to size rotorcraft is explored.

Johnson, Wayne↗

Method and system for fault accommodation of machines

A method for multi-objective fault accommodation using predictive modeling is disclosed. The method includes using a simulated machine that simulates a faulted actual machine, and using a simulated controller that simulates an actual controller. A multi-objective optimization process is performed, based on specified control settings for the simulated controller and specified operational scenarios for the simulated machine controlled by the simulated controller, to generate a Pareto frontier-based solution space relating performance of the simulated machine to settings of the simulated controller, including adjustment to the operational scenarios to represent a fault condition of the simulated machine. Control settings of the actual controller are adjusted, represented by the simulated controller, for controlling the actual machine, represented by the simulated machine, in response to a fault condition of the actual machine, based on the Pareto frontier-based solution space, to maximize desirable operational conditions and minimize undesirable operational conditions while operating the actual machine in a region of the solution space defined by the Pareto frontier.

Goebel, Kai Frank↗

On the Development and Demise of the Fall 2019 Southeast U. S. Flash Drought: Links to an Extreme Positive IOD

Much of the southeast United States experienced record dry conditions during September of 2019, with the area in abnormally dry to exceptional drought conditions growing from 25% at the beginning of the month to 80% by the end of the month. The drought ended just as abruptlydue to above normal rain that fell during the second half of October. In this study we employed MERRA-2 and the GEOS-5 AGCM to diagnose the underlying causes of the drought’s onset, maintenance, and demise. The basic approach involves performing a series of AGCM simulations in which the model is constrained to remain close to MERRA-2 over pre-specified areas that are external to the drought region. The start of the drought appears to have been forced by anomalous heating in the central/western tropical Pacific that resulted in low level anti-cyclonic flow and a tendency for descending motion over much of the southeast. An anomalous ridge associated with a Rossby wave train (emanating from the Indian Ocean region) is found to be the main source of the most intense temperature and precipitation anomalies that develop over the southeast during the last week of September. A second Rossby wave train (emanating from the same region) is responsible for the substantial rain that fell during the second half of October to end the drought. The links to the Indian Ocean Dipole (with record positive values) as well a sa waning El Nino allow some speculation as to the likelihood of similar events occurring in the future.

Siegfried D Schubert↗

SSME Turbopump Turbine Computations

A two-dimensional viscous code was developed to be used in the prediction of the flow in the SSME high-pressure turbopump blade passages. The rotor viscous code (RVC) employs a four-step Runge-Kutta scheme to solve the two-dimensional, thin-layer Navier-Stokes equations. The Baldwin-Lomax eddy-viscosity model is used for these turbulent flow calculations. A viable method was developed to use the relative exit conditions from an upstream blade row as the inlet conditions to the next blade row. The blade loading diagrams are compared with the meridional values obtained from an in-house quasithree-dimensional inviscid code. Periodic boundary conditions are imposed on a body-fitted C-grid computed by using the GRAPE GRids about Airfoils using Poisson's Equation (GRAPE) code. Total pressure, total temperature, and flow angle are specified at the inlet. The upstream-running Riemann invariant is extrapolated from the interior. Static pressure is specified at the exit such that mass flow is conserved from blade row to blade row, and the conservative variables are extrapolated from the interior. For viscous flows the noslip condition is imposed at the wall. The normal momentum equation gives the pressure at the wall. The density at the wall is obtained from the wall total temperature.

Jorgenson, P. G. E.↗

Impact of Improvements in Volcanic Implementation on Atmospheric Chemistry and Climate in the GISS-E2 Model

The representation of volcanic eruptions in climate models introduces some of the largest errors when evaluating historical simulations, partly due to the crude model parameterizations. We will show preliminary results from the Goddard Institute for Space Studies (GISS)-E2 model comparing traditional highly parameterized volcanic implementation (specified Aerosol Optical Depth, Effective Radius) to deploying the full aerosol microphysics module MATRIX and directly emitting SO2 allowing us the prognosically determine the chemistry and climate impact. We show a reasonable match in aerosol optical depth, effective radius, and forcing between the full aerosol implementation and reconstructions/observations of the Mt. Pinatubo 1991 eruption, with a few areas as targets for future improvement. This allows us to investigate not only the climate impact of the injection of volcanic aerosols, but also influences on regional water vapor, O3, and OH distributions. With the skill of the MATRIX volcano implementation established, we explore (1) how the height of the injection column of SO2 influence atmospheric chemistry and climate response, (2) how the initial condition of the atmosphere influences the climate and chemistry impact of the eruption with a particular focus on how ENSO and QBO and (3) how the coupled chemistry could mitigate the climate signal for much larger eruptions (i.e. the 1258 eruption, reconstructed to be approximately 10x Pinatubo). During each sensitivity experiment we assess the impact on profiles of water vapor, O3, and OH, and assess how the eruption impacts the budget of each.

aerosols↗

On the numerical solution of the dynamically loaded hydrodynamic lubrication of the point contact problem

The transient analysis of hydrodynamic lubrication of a point-contact is presented. A body-fitted coordinate system is introduced to transform the physical domain to a rectangular computational domain, enabling the use of the Newton-Raphson method for determining pressures and locating the cavitation boundary, where the Reynolds boundary condition is specified. In order to obtain the transient solution, an explicit Euler method is used to effect a time march. The transient dynamic load is a sinusoidal function of time with frequency, fractional loading, and mean load as parameters. Results include the variation of the minimum film thickness and phase-lag with time as functions of excitation frequency. The results are compared with the analytic solution to the transient step bearing problem with the same dynamic loading function. The similarities of the results suggest an approximate model of the point contact minimum film thickness solution.

Lim, Sang G.↗

Theoretical analysis of an augmentor wing for a VTOL fighter

A method based on potential flow theory was developed for predicting forces and moments acting on augmentor wings for prescribed ejector jet characteristics. A three dimensional nonplanar vortex lattice is laid out on the chordal planes of the augmentor wing components. Jet induced effects are included in the boundary condition from which the horseshoe vortex strengths are obtained. The jet within the diffusor is made to expand from the primary nozzles to the diffusor exit and is represented by a distribution of vorticity on the jet boundary to provide proper entrainment. The jet downstream of the diffusor exit is modeled by a vorticity distribution and blockage panels and its centerline location and spreading rate are taken from experimental data. The vortex lattice and jet models are used in an iterative manner until the predicted diffusor exit velocity matches the specified one. Some comparisons with available data show good agreement at lower power settings.

Dillenius, M. F. E.↗

Thermal modeling of a Ni-H2 battery cell

The nickel-hydrogen secondary battery has many desirable features which make it attractive for satellite power systems. It can provide a significant improvement over the energy density of present spacecraft nickel-cadnium batteries, combined with longer life, tolerance to overcharge and possibility of state-of-charge indication. However, to realize these advantages, accurate thermal modeling of nickel-hydrogen cells is required in order to properly design the battery pack so that it operates within a specified temperature range during the operation. Maintenance of a low operating temperature and a uniform temperature profile within the cell will yield better reliability, improved cycle life and better charge/discharge efficiencies. This research has the objective of developing and testing a thermal model which can be used to characterize battery operation. Primarily, temperature distribution with the heat generation rates as a function of position and time will be evaluated for a Ni-H2 cell in the three operating modes: (1) charge cycle, (2) discharge cycle, and (3) overcharge condition, if applicable. Variables to be examined include charging current, discharge rates, state of charge, pressure and temperature. Once the thermal model has been developed, this resulting model will predict the actual operating temperature and temperature gradient for the specific cell geometry to be used.

Ryu, Si-Ok↗

Prediction of Turbulence-Generated Noise in Unheated Jets: JeNo Users' Manual (Version 1.0) - Part 2

JeNo (Version 1.0) is a Fortran90 computer code that calculates the far-field sound spectral density produced by axisymmetric, unheated jets at a user specified observer location and frequency range. The user must provide a structured computational grid and a mean flow solution from a Reynolds-Averaged Navier Stokes (RANS) code as input. Turbulence kinetic energy and its dissipation rate from a k-epsilon or k-omega turbulence model must also be provided. JeNo is a research code, and as such, its development is ongoing. The goal is to create a code that is able to accurately compute far-field sound pressure levels for jets at all observer angles and all operating conditions. In order to achieve this goal, current theories must be combined with the best practices in numerical modeling, all of which must be validated by experiment. Since the acoustic predictions from JeNo are based on the mean flow solutions from a RANS code, quality predictions depend on accurate aerodynamic input.This is why acoustic source modeling, turbulence modeling, together with the development of advanced measurement systems are the leading areas of research in jet noise research at NASA Glenn Research Center.

Khavaran, Abbas↗

Understanding and Utilizing PBL Height Data from Multiple Observing Systems in the GEOS System

The accuracy of PBL height simulation is a key issue in many applications including forecasting near surface meteorology and air quality, however, it is a very challenging problem due to the lack of not only comprehensive, global Planetary Boundary Layer (PBL) observations but also a strategy and infrastructure to utilize PBL height data from a variety of sensors. Following the designation of PBL as an incubation class observable in the 2017 Decadal Survey, the PBL Incubation Study Team Report [14] made clear that “a future global PBL observing system requires modeling and data assimilation as essential components.” There is an urgent need for global modeling development in order to utilize Program of Record (POR) observations, assess their impacts, and identify gaps to be filled by future PBL missions. Our overall objective is to develop PBL data assimilation capabilities in the NASA Global Earth Observing System (GEOS), focusing on PBL height from multiple observing systems, to support the assessment and use of future PBL observations. The NASA GEOS system is composed of the GEOS global atmospheric general circulation model (AGCM) and the atmospheric data assimilation system (ADAS). The PBL parameterizations include the “Lock” K-profile scheme driven by surface and cloud-top buoyancy fluxes ([4]), and the “Louis” local scheme for stable conditions based on the Richardson number ([5]). Above the mixed layer defined by the Lock surface plume, shallow cumulus convection is represented by the mass flux scheme of [9]. Additional parameterizations are summarized in [1]. The ADAS employs the hybrid 4D Ensemble- Variational (EnVar) configuration ([15]), with the ensemble providing flow-dependent background error covariance information. The resultant analysis increments are fed back to the forecast model through the 4D incremental analysis update (IAU) approach ([11]). In this study, PBL height data are being or have been generated from radiosondes, GNSS RO, satellite (CATS, CALIPSO and ICESat-2) and ground-based (MPLNET) lidars, and wind profiler. Investigations have been conducted to specify quality marks for PBL height retrievals for the data assimilation purpose. These PBL height data have different strengths and weaknesses ([2], [3], [6], [7], [8], [10]), and the satellite PBL height data provide better global coverage and complement in-situ PBL height data. Radiosondes offer high accuracy and in situ measurement of temperature and humidity profiles, but with poor spatio-temporal sampling. The in-situ observing systems like MPLNET and wind profiler provide long history of PBL height records at each station. The GNSS RO based PBL height is retrieved based on the sharp gradients in refractivity profile that represent the fine vertical structure of temperature and moisture changes above the PBL. However, not all RO refractivity profiles reach the surface depending on location and regime, and RO refractivity retrievals can be negatively biased below 2km. The PBL height data from satellite lidars provide high resolution along track PBL height retrievals, but over land they are affected by previous day convective PBL aerosol and strongly associated with mixing layer and retrievals cannot be made below thick, attenuating clouds. A successful assimilation of PBL height data requires a thorough understanding of the observing method and the retrieval algorithm for each observing system in order to use the PBL height data from multiple observing systems properly. Due to the sensitivity of PBL height data to the observing method and choice of algorithm, it is important to use a model definition appropriate for each observation type to compute differences between PBL height data and model PBL height (OmFs). The GEOS model currently includes two PBL height definitions suitable for direct comparison with observed PBL height, and additional definitions are being added in this study. Evaluation of different model PBL height definitions is underway. Meanwhile, efforts have been made in the GEOS data assimilation system to develop PBL height data assimilation capability. PBL height data can be assimilated using two different approaches. The traditional approach is to construct an observation operator and its tangent linear and adjoint, which link control variables to PBL height data from each observing system. This observation operator can be very complicated, e.g., the lidar-based PBL height observation operator includes the backscatter lidar forward observation operator, the algorithm to derive PBL height from attenuated total backscatter, interpolation, and calculations handling the mismatch between observed and model scales. The other approach is to augment PBL height to the control variable vector, and it is adopted in this study. The latter approach was also used in previous studies, e.g., the assimilation of PBL height data from radiosonde and aircraft in the Real Time Mesoscale Analysis (RTMA) system for a dispersion modelling study ([13]); the PBL height assimilation study using lidar PBL height data at Greensburg, Kansas for a field campaign ([12]). The PBL height assimilation from multiple observing systems in this study allows us to take advantage of the diverse PBL height data that provide much better global coverage collectively under different meteorological conditions and with different temporal and spatial scales. As all the PBL heights are tightly coupled with the PBL thermodynamic variables, the strong correlations, which are provided by the 4D ensemble forecast, enable PBL height data from various sources to interact and combine coherently and provide additional information for PBL temperature and moisture fields. The results of comparisons among PBL height data from different sources and the evaluation of the model PBL height definitions with the PBL height data will be presented, and the PBL height data synergy strategies and preliminary results will also be discussed at the conference.

Y. Zhu↗

Development of an Accelerated Test Design for Predicting the Service Life of the Solar Array at Mead, Nebraska

An accelerated life test is described which was developed to predict the life of the 25 kW photovoltaic array installed near Mead, Nebraska. A quantitative model for accelerating testing using multiple environmental stresses was used to develop the test design. The model accounts for the effects of thermal stress by a relation of the Arrhenius form. This relation was then corrected for the effects of nonthermal environmental stresses, such as relative humidity, atmospheric pollutants, and ultraviolet radiation. The correction factors for the nonthermal stresses included temperature-dependent exponents to account for the effects of interactions between thermal and nonthermal stresses on the rate of degradation of power output. The test conditions, measurements, and data analyses for the accelerated tests are presented. Constant-temperature, cyclic-temperature, and UV types of tests are specified, incorporating selected levels of relative humidity and chemical contamination and an imposed forward-bias current and static electric field.

Gaines, G. B.↗

Mass transport phenomena between bubbles and dissolved gases in liquids under reduced gravity conditions

The experimental and analytical work that was done to establish justification and feasibility for a shuttle middeck experiment involving mass transfer between a gas bubble and a liquid is described. The experiment involves the observation and measurement of the dissolution of an isolated immobile gas bubble of specified size and composition in a thermostatted solvent liquid of known concentration in the reduced gravity environment of earth orbit. Methods to generate and deploy the bubble were successful both in normal gravity using mutually buoyant fluids and under reduced gravity conditions in the NASA Lear Jet. Initialization of the experiment with a bubble of a prescribed size and composition in a liquid of known concentration was accomplished using the concept of unstable equilibrium. Subsequent bubble dissolution or growth is obtained by a step increase or decrease in the liquid pressure. A numerical model was developed which simulates the bubble dynamics and can be used to determine molecular parameters by comparison with the experimental data. The primary objective of the experiment is the elimination of convective effects that occur in normal gravity.

Dewitt, Kenneth J.↗

Calculation of the Pressure Distribution on Bodies of Revolution in the Subsonic Flow of a Gas: Axially Symmetrical Flow - Part 1

The present report concerns a method of computing the velocity and pressure distributions on bodies of revolution in axially symmetrical flow in the subsonic range. The differential equation for the velocity potential Phi of a compressible fluid motion is linearized tn the conventional manner, and then put in the form Delta(Phi) = 0 by affine transformation. The quantity Phi represents the velocity potential of a fictitious incompressible flow, for which a constant superposition of sources by sections is secured by a method patterned after von Karman which must comply with the boundary condition delta(phi)/delta(n) = 0 at the originally specified contour. This requirement yields for the "pseudo-stream function" psi a differential equation which must be fulfilled for as many points on the contour as source lengths are assumed. In this manner, the problem of defining the still unknown source intensities is reduced to the solution of an inhomogeneous equation system. The pressure distribution is then determined with the aid of Bernoulli's equation and adiabatic equation of state. Lastly, the pressure distributions in compressible and incompressible medium are compared on a model problem.

Bilharz, Herbert↗