Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reference 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 1,297 records · Page 72

Development and Implementation of an Empirical Ionosphere Variability Model

Spacecraft designers and operations support personnel involved in space environment analysis for low Earth orbit missions require ionospheric specification and forecast models that provide not only average ionospheric plasma parameters for a given set of geophysical conditions but the statistical variations about the mean as well. This presentation describes the development of a prototype empirical model intended for use with the International Reference Ionosphere (IRI) to provide ionospheric Ne and Te variability. We first describe the database of on-orbit observations from a variety of spacecraft and ground based radars over a wide range of latitudes and altitudes used to obtain estimates of the environment variability. Next, comparison of the observations with the IRI model provide estimates of the deviations from the average model as well as the range of possible values that may correspond to a given IRI output. Options for implementation of the statistical variations in software that can be run with the IRI model are described. Finally, we provide example applications including thrust estimates for tethered satellites and specification of sunrise Ne, Te conditions required to support spacecraft charging issues for satellites with high voltage solar arrays.

Minow, Joesph I.↗

An extended-source spatial acquisition process based on maximum likelihood criterion for planetary optical communications

This paper describes an extended-source spatial acquisition process based on the maximum likelihood criterion for interplanetary optical communications. The objective is to use the sun-lit Earth image as a receiver beacon and point the transmitter laser to the Earth-based receiver to establish a communication path. The process assumes the existence of a reference image. The uncertainties between the reference image and the received image are modeled as additive white Gaussian disturbances. It has been shown that the optimal spatial acquisition requires solving two nonlinear equations to estimate the coordinates of the transceiver from the received camera image in the transformed domain. The optimal solution can be obtained iteratively by solving two linear equations. Numerical results using a sample sun-lit Earth as a reference image demonstrate that sub-pixel resolutions can be achieved in a high disturbance environment. Spatial resolution is quantified by Cramer-Rao lower bounds.

Yan, Tsun-Yee↗

Extracting Damping Ratio from Dynamic Data and Numerical Solutions

There are many ways to extract damping parameters from data or models. This Technical Memorandum provides a quick reference for some of the more common approaches used in dynamics analysis. Described are six methods of extracting damping from data: the half-power method, logarithmic decrement (decay rate) method, an autocorrelation/power spectral density fitting method, a frequency response fitting method, a random decrement fitting method, and a newly developed half-quadratic gain method. Additionally, state-space models and finite element method modeling tools, such as COMSOL Multiphysics (COMSOL), provide a theoretical damping via complex frequency. Each method has its advantages which are briefly noted. There are also likely many other advanced techniques in extracting damping within the operational modal analysis discipline, where an input excitation is unknown; however, these approaches discussed here are objective, direct, and can be implemented in a consistent manner.

Casiano, M. J.↗

Air-water analogy and the study of hydraulic models

The author first sets forth some observations about the theory of models. Then he established certain general criteria for the construction of dynamically similar models in water and in air, through reference to the perfect fluid equations and to the ones pertaining to viscous flow. It is, in addition, pointed out that there are more cases in which the analogy is possible than is commonly supposed.

SIMILITUDE↗

A seal test facility for the measurement of isotropic and anisotropic linear rotordynamic characteristics

A new seal test facility for measuring high-pressure seal rotor-dynamic characteristics has recently been made operational at Case Western Reserve University (CWRU). This work is being sponsored by the Electric Power Research Institute (EPRI). The fundamental concept embodied in this test apparatus is a double-spool-shaft spindle which permits independent control over the spin speed and the frequency of an adjustable circular vibration orbit for both forward and backward whirl. Also, the static eccentricity between the rotating and non-rotating test seal parts is easily adjustable to desired values. By accurately measuring both dynamic radial displacement and dynamic radial force signals, over a wide range of circular orbit frequency, one is able to solve for the full linear-anisotropic model's 12 coefficients rather than the 6 coefficients of the more restrictive isotropic linear model. Of course, one may also impose the isotropic assumption in reducing test data, thereby providing a valid qualification of which seal configurations are well represented by the isotropic model and which are not. In fact, as argued in reference (1), the requirement for maintaining a symmetric total system mass matrix means that the resulting isotropic model needs 5 coefficients and the anisotropic model needs 11 coefficients.

Adams, M. L.↗

Transmission Telescope Optical Dynamic Alignment

The Integrated Radio Optical Communication System (iROC) is designed to transmit data between Mars and Earth by means of radio waves at 32.67 GHz (Ka Band) and laser beam at 1550 nm, both transmitted via a combined teletenna. The iROC terminal will provide ‘beaconless’ operations to allow full function from the outer planets. In order to point without the aid of an uplink beacon, the proof of concept presented here is addressing the need for an accurate determination and control of the relative position of the laser beam with respect to a reference star. The experiment presented simulates a surrogate transmission telescope system in laboratory setting and presents the model used in the correction of the outgoing beam. The results of the model show a nonlinear dependence between the outgoing and the reference beam, indicating the necessity of a minimum of two metrology instruments placed along the optical system for increased pointing precision.

optical communication↗

Four-D global reference atmosphere technical description, part 1

An empirical atmospheric model was developed which generates values for pressure, density, temperature, and winds from surface levels to orbital altitudes. The output parameters consist of components for: (1) latitude, longitude, and altitude dependent monthly and annual means; (2) quasi-biennial oscillations; and (3) random perturbations to simulate partially the variability due to synoptic, diurnal, planetary wave, and gravity wave variations. Quasi-biennial and random variation perturbations are computed from parameters determined from various empirical studies and are added to the monthly mean values. This model has been developed as a computer program called PROFILE which can be used to generate altitude profiles of atmospheric parameters along any simulated trajectory through the atmosphere. The PROFILE program was developed for design applications in the space shuttle program. Other applications of the model are discussed, such as for global circulation and diffusion studies, and for generating profiles for comparison with other atmospheric measurement techniques, (e.g. satellite measured temperature profiles).

Justus, C. G.↗

On identified predictive control

Self-tuning control algorithms are potential successors to manually tuned PID controllers traditionally used in process control applications. A very attractive design method for self-tuning controllers, which has been developed over recent years, is the long-range predictive control (LRPC). The success of LRPC is due to its effectiveness with plants of unknown order and dead-time which may be simultaneously nonminimum phase and unstable or have multiple lightly damped poles (as in the case of flexible structures or flexible robot arms). LRPC is a receding horizon strategy and can be, in general terms, summarized as follows. Using assumed long-range (or multi-step) cost function the optimal control law is found in terms of unknown parameters of the predictor model of the process, current input-output sequence, and future reference signal sequence. The common approach is to assume that the input-output process model is known or separately identified and then to find the parameters of the predictor model. Once these are known, the optimal control law determines control signal at the current time t which is applied at the process input and the whole procedure is repeated at the next time instant. Most of the recent research in this field is apparently centered around the LRPC formulation developed by Clarke et al., known as generalized predictive control (GPC). GPC uses ARIMAX/CARIMA model of the process in its input-output formulation. In this paper, the GPC formulation is used but the process predictor model is derived from the state space formulation of the ARIMAX model and is directly identified over the receding horizon, i.e., using current input-output sequence. The underlying technique in the design of identified predictive control (IPC) algorithm is the identification algorithm of observer/Kalman filter Markov parameters developed by Juang et al. at NASA Langley Research Center and successfully applied to identification of flexible structures.

Bialasiewicz, Jan T.↗

Numerical Determination of Critical Conditions for Thermal Ignition

The determination of ignition or thermal explosion in an oxidizing porous body of material, as described by a dimensionless reaction-diffusion equation of the form .tu = .2u + .e-1/u over the bounded region O, is critically reexamined from a modern perspective using numerical methodologies. First, the classic stationary model is revisited to establish the proper reference frame for the steady-state solution space, and it is demonstrated how the resulting nonlinear two-point boundary value problem can be reexpressed as an initial value problem for a system of first-order differential equations, which may be readily solved using standard algorithms. Then, the numerical procedure is implemented and thoroughly validated against previous computational results based on sophisticated path-following techniques. Next, the transient nonstationary model is attacked, and the full nonlinear form of the reaction-diffusion equation, including a generalized convective boundary condition, is discretized and expressed as a system of linear algebraic equations. The numerical methodology is implemented as a computer algorithm, and validation computations are carried out as a prelude to a broad-ranging evaluation of the assembly problem and identification of the watershed critical initial temperature conditions for thermal ignition. This numerical methodology is then used as the basis for studying the relationship between the shape of the critical initial temperature distribution and the corresponding spatial moments of its energy content integral and an attempt to forge a fundamental conjecture governing this relation. Finally, the effects of dynamic boundary conditions on the classic storage problem are investigated and the groundwork is laid for the development of an approximate solution methodology based on adaptation of the standard stationary model.

Luo, W.↗

Methodology for analysis and simulation of large multidisciplinary problems

The Integrated Structural Modeling (ISM) program is being developed for the Air Force Weapons Laboratory and will be available for Air Force work. Its goal is to provide a design, analysis, and simulation tool intended primarily for directed energy weapons (DEW), kinetic energy weapons (KEW), and surveillance applications. The code is designed to run on DEC (VMS and UNIX), IRIS, Alliant, and Cray hosts. Several technical disciplines are included in ISM, namely structures, controls, optics, thermal, and dynamics. Four topics from the broad ISM goal are discussed. The first is project configuration management and includes two major areas: the software and database arrangement and the system model control. The second is interdisciplinary data transfer and refers to exchange of data between various disciplines such as structures and thermal. Third is a discussion of the integration of component models into one system model, i.e., multiple discipline model synthesis. Last is a presentation of work on a distributed processing computing environment.

Russell, William C.↗

High-Fidelity Modeling Approach for Coupling Thermal, Fluid, and Neutronic Analysis in an NTP Engine

Nuclear thermal propulsion (NTP), with its increased efficiency compared to traditional rocket engines, is a potentially enabling technology for space missions of higher complexity, such as delivering heavy machinery or crew to the moon or other planets. Thus, exploring modeling approaches that can effectively couple the thermal, fluid, and neutronic behaviors of NTP reactors is essential to mitigate cost and schedule challenges with hardware development. This paper details one such approach which utilizes Cardinal to unify reactor physics analysis performed using OpenMC, heat transfer analysis performed using the MOOSE framework, and thermal hydraulic analysis performed using the MOOSE Thermal Hydraulic Module. This approach is demonstrated with a proof of concept that focuses on a legacy reactor geometry from the historic NERVA program. This preliminary model found an overprediction in, when compared to the reference data, axial power peaking and the temperatures of both solids and fluids by roughly 20%. It is clear, however, that the model conserved a number of broad trends, and a few model refinements should reduce the observed error. This approach has definite merit and should enable the production of informative results for more modern NTP systems. This work is a first step in creating a transient reactor model that can be updated across a sequence of time steps based on continually fluctuating, time-dependent boundary conditions.

Nuclear Thermal Propulsion↗

Waterhammer Testing and Modeling of the Ares I Upper Stage Reaction Control System

NASA's Ares I rocket is the agency's first step in completing the goals of the Constellation Program, which plans to deliver a new generation of space explorers into low earth orbit for future missions to the International Space Station, the moon, and other destinations within the solar system. Ares I is a two-stage rocket topped by the Orion crew capsule and its service module. The launch vehicle's First Stage is a single, five-segment reusable solid rocket booster (RSRB), derived from the Space Shuttle Program's four segment RSRB. The vehicle's Upper Stage, being designed at Marshall Space Flight Center (MSFC), is propelled by a single J-2X Main Engine fueled with liquid oxygen and liquid hydrogen. During active Upper Stage flight of the Ares I launch vehicle, the Upper Stage Reaction Control System (US ReCS) will perform attitude control operations for the vehicle. The US ReCS will provide three-axis attitude control capability (roll, pitch, and yaw) for the Upper Stage while the J-2X is not firing and roll control capability while the engine is firing. Because of the requirements imposed upon the system, the design must accommodate rapid pulsing of multiple thrusters simultaneously to maintain attitude control. In support of these design activities and in preparation for Critical Design Review, analytical models of the US ReCS propellant feed system have been developed using the Thermal Hydraulic Library of MSC.EASY5 v.2008, herein referred to as EASY5. EASY5 is a commercially available fluid system modeling package with significant history of modeling space propulsion systems. In Fall 2009, a series of development tests were conducted at MSFC on a cold-flow test article for the US ReCS, herein referred to as System Development Test Article (SDTA). A subset of those tests performed were aimed at examining the effects of waterhammer on a flight-representative system and to ensure that those effects could be quantified with analytical models and incorporated into the design of the flight system. This paper presents an overview of the test article and the test approach, along with a discussion of the analytical modeling methodology. In addition, the results of that subset of development tests, along with analytical model pre-test predictions and post-test model correlations, will also be discussed in detail.

Williams, J. Hunter↗

On the amplitudes reached by baroclinically unstable disturbances

The maximum eddy energies reached by baroclinically unstable disturbances in some simplified numerical models (Gall, 1976; Simmons and Hoskins, 1978) are compared to a scale-dependent measure of the energy available to the disturbances. It is found that, except for the longest wavelengths, the ratio of the maximum eddy energies to this available energy (the eddy efficiency) tends to remain constant relative to the large variations in the maximum energies. The roles of the growth rate of the initial disturbance and the beta-effect are discussed with reference to results from a maximally simplified model. It is shown that dimensional arguments paralleling Charney (1971) lead to the maximum energy versus scale relationships found in the numerical results.

Schneider, E. K.↗

Predicting NonInertial Effects with Algebraic Stress Models which Account for Dissipation Rate Anisotropies

Three types of turbulence models which account for rotational effects in noninertial frames of reference are evaluated for the case of incompressible, fully developed rotating turbulent channel flow. The different types of models are a Coriolis-modified eddy-viscosity model, a realizable algebraic stress model, and an algebraic stress model which accounts for dissipation rate anisotropies. A direct numerical simulation of a rotating channel flow is used for the turbulent model validation. This simulation differs from previous studies in that significantly higher rotation numbers are investigated. Flows at these higher rotation numbers are characterized by a relaminarization on the cyclonic or suction side of the channel, and a linear velocity profile on the anticyclonic or pressure side of the channel. The predictive performance of the three types of models are examined in detail, and formulation deficiencies are identified which cause poor predictive performance for some of the models. Criteria are identified which allow for accurate prediction of such flows by algebraic stress models and their corresponding Reynolds stress formulations.

Jongen, T.↗

Production of Lightning NO(x) and its Vertical Distribution Calculated from 3-D Cloud-scale Chemical Transport Model Simulations

A 3-D cloud scale chemical transport model that includes a parameterized source of lightning NO(x), based on observed flash rates has been used to simulate six midlatitude and subtropical thunderstorms observed during four field projects. Production per intracloud (P(sub IC) and cloud-to-ground (P(sub CG)) flash is estimated by assuming various values of P(sub IC) and P(sub CG) for each storm and determining which production scenario yields NO(x) mixing ratios that compare most favorably with in-cloud aircraft observations. We obtain a mean P(sub CG) value of 500 moles NO (7 kg N) per flash. The results of this analysis also suggest that on average, P(sub IC) may be nearly equal to P(sub CG), which is contrary to the common assumption that intracloud flashes are significantly less productive of NO than are cloud-to-ground flashes. This study also presents vertical profiles of the mass of lightning NO(x), after convection based on 3-D cloud-scale model simulations. The results suggest that following convection, a large percentage of lightning NO(x), remains in the middle and upper troposphere where it originated, while only a small percentage is found near the surface. The results of this work differ from profiles calculated from 2-D cloud-scale model simulations with a simpler lightning parameterization that were peaked near the surface and in the upper troposphere (referred to as a "C-shaped" profile). The new model results (a backward C-shaped profile) suggest that chemical transport models that assume a C-shaped vertical profile of lightning NO(x) mass may place too much mass neat the surface and too little in the middle troposphere.

Ott, Lesley↗

Global mean sea surface computation using GEOS 3 altimeter data

A mean sea surface map has been determined for the global ocean areas between +62 deg and -62 deg latitude using GEOS 3 altimeter data. A grid of laser reference orbits computed using the GEM 10B gravity model has been used to orient the altimeter data in a center of mass coordinate system. The density of the altimeter tracks has enabled the computation of the sea surface heights above the reference ellipsoid on 1 deg x 1 deg grid in most of the oceanic areas. In the northwest Atlantic the dense coverage has enabled computations on a 0.25 deg x 0.25 deg grid. Comparisons of the global surface with an independently computed mean sea surface based upon SEASAT altimeter data indicate an rms agreement of a little over a meter. Comparisons of the regional solution in the northwest Atlantic with SEASAT profiles indicate a precision of a few decimeters in this surface. An analysis of the global crossover differences has indicated the possibility of a timing error in the altimeter data. Timing bias values of within 0.5 msec of 9.2 msec for 1975 and within 0.6 msec of 18.7 msec for 1976 have been recovered from an analysis of the altimeter data.

Marsh, J. G.↗

Buckling Testing and Analysis of Honeycomb Sandwich Panel Arc Segments of a Full-Scale Fairing Barrel Part 1: 8-Ply In-Autoclave Facesheets: 8-Ply In-Autoclave Facesheets - Part 1

Four honeycomb sandwich panels, representing 1/16th arc segments of a 10-m diameter barrel section of the heavy lift launch vehicle, were manufactured under the NASA Composites for Exploration program and the NASA Space Launch Systems program. Two configurations were chosen for the panels: 6-ply facesheets with 1.125 in. honeycomb core and 8-ply facesheets with 1.000 in. honeycomb core. Additionally, two separate carbon fiber/epoxy material systems were chosen for the facesheets: inautoclave IM7/977-3 and out-of-autoclave T40-800b/5320-1. Smaller 3.00- by 5.00-ft panels were cut from the 1/16th barrel sections. These panels were tested under compressive loading at the NASA Langley Research Center. Furthermore, linear eigenvalue and geometrically nonlinear finite element analysis was performed to predict the compressive response of the 3.00- by 5.00-ft panels. This manuscript summarizes the experimental and analytical modeling efforts pertaining to the panel composed of 8-ply, IM7/977-3 facesheets (referred to Panel A). To improve the robustness of the geometrically nonlinear finite element model, measured surface imperfections were included in the geometry of the model. Both the linear and nonlinear models yield good qualitative and quantitative predictions. Additionally, it was predicted correctly that the panel would fail in buckling prior to failing in strength. Furthermore, several imperfection studies were performed to investigate the influence of geometric imperfections, fiber misalignments, and three-dimensional (3 D) effects on the compressive response of the panel.

Myers, David E.↗