Search NASA⌕ Search

SEARCH · Search NASA

Results for “Equilibrium programming”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS1975

Computer program system for dynamic simulation and stability analysis of passive and actively controlled spacecraft. Volume 1. Theory

A theoretical development and associated digital computer program system is presented. The dynamic system (spacecraft) is modeled as an assembly of rigid and/or flexible bodies not necessarily in a topological tree configuration. The computer program system may be used to investigate total system dynamic characteristics including interaction effects between rigid and/or flexible bodies, control systems, and a wide range of environmental loadings. Additionally, the program system may be used for design of attitude control systems and for evaluation of total dynamic system performance including time domain response and frequency domain stability analyses. Volume 1 presents the theoretical developments including a description of the physical system, the equations of dynamic equilibrium, discussion of kinematics and system topology, a complete treatment of momentum wheel coupling, and a discussion of gravity gradient and environmental effects. Volume 2, is a program users' guide and includes a description of the overall digital program code, individual subroutines and a description of required program input and generated program output. Volume 3 presents the results of selected demonstration problems that illustrate all program system capabilities.

Bodley, C. S.↗
NASA NTRS1979

Radiative transfer in spherical dust shells using a generalized two-stream Eddington approximation

Application of a generalized two-stream Eddington approximation to the problem of radiative transfer in extended, spherically symmetric dust shells is presented. It is assumed that the radiation field can be characterized by the mean intensity, the flux, and a positionally dependent direction cosine specifying the division into two solid-angle ranges. The direction cosine is not specified a priori and is a function of the geometry, opacity, and emissivity in the dust shell. A multiple-grain-size multiple-temperature-distribution dust shell is postulated in which isotropic and anisotropic scattering as well as absorption and thermal reemission are allowed. A program has been developed that solves for the multiple temperature distributions by applying the constraint of radiative equilibrium to each grain size, and then calculates emergent fluxes. Results of one such calculation are presented for a model dust shell having a maximum optical depth (approximately 41) in the visible, clearly showing large optical extinction and a moderate infrared excess.

Haisch, B. M.↗
NASA NTRS1994

The metallicity of M4: Accurate spectroscopic fundamental parameters for four giants

High-quality spectra, covering the wavelength range 5480 to 7080 A, have been obtained for four giant stars in the intermediate-metallicity CN-bimodal globular cluster M4 (NGC 6121). We have employed a model atmosphere analysis that is entirely independent from cluster parameters, such as distance, age, and reddening, in order to derive accurate values for the stellar parameters effective temperature, surface gravity, and microturbulence, and for the abundance of iron relative to the Sun, (Fe/H), and of calcium, Ca/H, for each of the four stars. Detailed radiative transfer and statistical equilibrium calculations carried out for iron and calcium suggest that departures from local thermodynamic equilibrium are not significant for the purposes of our analysis. The spectroscopically derived effective temperatures for our program stars are hotter by about 200 K than existing photometric calibrations suggest. We conclude that this is due partly to the uncertain reddening of M4 and to the existing photometric temperature calibration for red giants being too cool by about 100 K. Comparison of our spectroscopic and existing photometric temperatures supports the prognosis of a significant east-west gradient in the reddening across M4. Our derived iron abundances are slightly higher than previous high-resolution studies suggested; the differences are most probably due to the different temperature scale and choice of microturbulent velocities adopted by earlier workers. The resulting value for the metallicity of M4 is (Fe/H )(sub M4) = -1.05 + or - 0.15. Based on this result, we suggest that metallicities derived in previous high-dispersion globular cluster abundance analyses could be too low by 0.2 to 0.3 dex. Our calcium abundances suggest an enhancement of calcium, an alpha element, over iron, relative to the Sun, in M4 of (Ca/H) = 0.23.

Drake, J. J.↗
NASA NTRS1977

General instability analysis

HOLBOAT computer program provides instability analysis of inhomogeneous, anisotropic, right-circular cylinder or segment under combined loading. Program is based on Kirchhoff-Love hypothesis, general anisotropic constitutive equations, and Flugge's differential equations of equilibrium.

Bianca, C. J.↗
NASA NTRS1992

Swept shock/boundary layer interaction experiments in support of CFD code validation

Research on the topic of shock wave/turbulent boundary-layer interaction was carried out during the past three years at the Penn State Gas Dynamics Laboratory. This report describes the experimental research program which provides basic knowledge and establishes new data on heat transfer in swept shock wave/boundary-layer interactions. An equilibrium turbulent boundary-layer on a flat plate is subjected to impingement by swept planar shock waves generated by a sharp fin. Five different interactions with fin angle ranging from 10 deg to 20 deg at freestream Mach numbers of 3.0 and 4.0 produce a variety of interaction strengths from weak to very strong. A foil heater generates a uniform heat flux over the flat plate surface, and miniature thin-film-resistance sensors mounted on it are used to measure the local surface temperature. The heat convection equation is then solved for the heat transfer distribution within an interaction, yielding a total uncertainty of about +/- 10 percent. These experimental data are compared with the results of numerical Navier-Stokes solutions which employ a k-epsilon turbulence model. Finally, a simplified form of the peak heat transfer correlation for fin interactions is suggested.

Settles, G. S.↗
NASA NTRS1993

On the preferred source location for the convective amplification of ion cyclotron waves

The propagation, growth, and absorption of electromagnetic ion cyclotron waves (EMICs) in the Pc1 range are studied using the HOTRAY ray tracing program for a realistic distribution of thermal plasma (H+, He+, O+), which is assumed to be in diffusive equilibrium inside the plasmasphere and collisionless in the low-density region outside the plasmapause. It is demonstrated that there are two principal source regions for the growth of EMIC waves. The theoretical results are related to the most recent satellite and ground-based observations.

Horne, Richard B.↗
NASA NTRS1994

Heat-transfer measurements and computations of swept-shock-wave boundary-layer interactions

An experimental and computational research program providing new knowledge of the heat transfer in swept-shock-wave/boundary-layer interactions is described. An equilibrium turbulent boundary layer on a flat plate is subjected to impingement by a swept planar shock wave generated by a sharp fin. Five different interactions with fin angles ranging from 10 to 20 deg at freestream Mach numbers of 3 and 4 produce a variety of interaction strengths ranging from weak to very strong. A foil heater generates a uniform heat flux over the flat plate surface, and miniature thin-film-resistance sensors are used to measure the local surface temperature. The heat convection equation is then solved for the heat transfer distribution within an interaction, yielding an uncertainty of about +/- 10%. These data are compared with numerical Navier-Stokes solutions that employ a k-epsilon turbulence model. A simple peak heat transfer correlation for fin interactions is suggested.

Lee, Y.↗
NASA NTRS1991

The Navy/NASA Engine Program (NNEP89): A user's manual

An engine simulation computer code called NNEP89 was written to perform 1-D steady state thermodynamic analysis of turbine engine cycles. By using a very flexible method of input, a set of standard components are connected at execution time to simulate almost any turbine engine configuration that the user could imagine. The code was used to simulate a wide range of engine cycles from turboshafts and turboprops to air turborockets and supersonic cruise variable cycle engines. Off design performance is calculated through the use of component performance maps. A chemical equilibrium model is incorporated to adequately predict chemical dissociation as well as model virtually any fuel. NNEP89 is written in standard FORTRAN77 with clear structured programming and extensive internal documentation. The standard FORTRAN77 programming allows it to be installed onto most mainframe computers and workstations without modification. The NNEP89 code was derived from the Navy/NASA Engine program (NNEP). NNEP89 provides many improvements and enhancements to the original NNEP code and incorporates features which make it easier to use for the novice user. This is a comprehensive user's guide for the NNEP89 code.

Plencner, Robert M.↗
NASA NTRS1988

An integrated approach to the synthesis of geometrically non-linear structures

An integrated approach to the minimum weight design of geometrically nonlinear three-dimensional truss structures with geometric imperfections, subject to inequality constraints on static displacements, stresses, local buckling and cross sectional areas, is investigated. The integrated structural synthesis problem involves design and response quantities as independent variables and equilibrium equations, describing the finite element model, as equality constraints. The nonlinear structural analysis and the optimization are thus merged together into a single process. A computer program developed to compute the constraint values and analytical gradients is coupled with a generalized reduced gradient algorithm to solve the integrated problem. Numerical results for a geometrically nonlinear shallow dome example problem are presented for various types of imperfections. Furthermore, it is found that the algorithm is capable of detecting and guarding against system as well as element elastic instability using equilibrium information only, that is, without imposing system and local buckling inequality constraints.

Smaoui, H.↗
NASA NTRS1996

A Parametric Study of a Plug Nozzle, Using the Liquid Propellant Program (LPP) Code

The Liquid Propellant Program (LPP) computer code is a super-set of the industry standard Two Dimensional Kinetics (TDK) computer code. The TDK code uses a two dimensional method of characteristics solution with fully coupled finite rate kinetics for axially symmetric nozzles. The chemical reactions are modeled with a generalized reaction package that include three dimensional body efficiencies and four reaction rate forms. The code performs optional solutions for frozen or equilibrium flow. TDK evaluates discrete shocks, both attached or induced. The Transonic module models variable mixture ratio profiles from the combustion chamber injector. The Mass Addition Boundary Layer module (MABL) calculates the boundary parameters with the same chemistry options, and includes transpiration or tangential slot injection of gas at the wall. The LPP upgrades include: planar nozzle, scarfed nozzles, plug nozzles, and scramjet nozzle configurations. The code evaluates both upper and lower wall flow simulation, and includes the interaction with the external flow. The MABL module evaluates equilibrium radiation heat transfer for both upper and lower walls. In addition, LPP code models combustion effects due to injector inefficiencies with the Spray Combustion Analysis Program (SCAP) module. The LPP package provides extensive post plotting capabilities for flow visualization. The LPP is sufficiently fast and robust to provide performance predictions for extensive parametric studies and sufficiently accurate to provide flow field and performance solutions for detailed studies.

Dunn, Stuart S.↗
NASA NTRS1983

Postglacial rebound observed by Lageos and the effective viscosity of the lower mantle

Sixty-four observations of the orbital node made by the Lageos satellite over a five year time interval reveal an acceleration of (-8.1 + or 1.8) x 10 to the -8 power arcseconds day/2 due to a source which is not presently modeled in the GEODYN orbit determination computer program. This acceleration cannot be explained by the ocean tide with 18.6 year period, assuming it to be an equilibrium tide. Instead it seems to be due to postglacial rebound, which changes the J(2) coefficient in the spherical harmonic expansion of the earth's gravitational field at the rate of (-8.2 + or - 18) x 10 to the -19th power/s; this in turn accelerates the node. This rate does not agree with the -32 x 10 to the -19th power/s predicted by Wu and Peltier's (1982) L2 model, which has upper and lower mantle effective viscosities of 10 to the 21st and 22nd powers Pa's, respectively. it does agree well with their L1 model, which gives about 10 x 10 to the 19th power/s. Since the effective viscosity is 10 to the 21st power Pa s throughout the entire mantle in the L1 model, the results support the contentions that the effective viscosity is near 10 to the 21st power Pa s everywhere in the mantle, and this relatively low value for the effective viscosity may have permitted several degrees of polar wander due to glaciation during the Quaternary Ice Age.

Rubincam, D. P.↗
NASA NTRS1984

Postglacial rebound observed by Lageos and the effective viscosity of the lower mantle

Sixty-four observations of the orbital node made by the Lageos satellite over a five year time interval reveal an acceleration of (-8.1 + or - 1.8) x 10 to the -8 power arcseconds day/2 due to a source which is not presently modeled in the GEODYN orbit determination computer program. This acceleration cannot be explained by the ocean tide with 18.6 year period, assuming it to be an equilibrium tide. Instead it seems to be due to postglacial rebound, which changes the J(2) coefficient in the spherical harmonic expansion of the earth's gravitational field at the rate of (-8.2 + or - 18) x 10 to the -19th power/s; this in turn accelerates the node. This rate does not agree with the -32 x 10 to the -19th power/s predicted by Wu and Peltier's (1982) L2 model, which has upper and lower mantle effective viscosities of 10 to the 21st and 22nd powers Pa's, respectively. It does agree well with their L1 model, which gives about 10 x 10 to the 19th power/s. Since the effective viscosity is 10 to the 21st power Pa s throughout the entire mantle in the L1 model, the results support the contentions that the efective viscosity is near 10 to the 21st power Pa s everyhere in the mantle, and this relatively low value for the effective viscosity may have permitted several degrees of polar wander due to glaciation during the Quaternary Ice Age. Previously announced in STAR as N84-13705

Rubincam, D. P.↗
NASA NTRS1988

Fluid-thermal-structural study of aerodynamically heated leading edges

A finite element approach for integrated fluid-thermal-structural analysis of aerodynamically heated leading edges is presented. The Navier-Stokes equations for high speed compressible flow, the energy equation, and the quasi-static equilibrium equations for the leading edge are solved using a single finite element approach in one integrated, vectorized computer program called LIFTS. The fluid-thermal-structural coupling is studied for Mach 6.47 flow over a 3-in diam cylinder for which the flow behavior and the aerothermal loads are calibrated by experimental data. Issues of the thermal-structural response are studied for hydrogen-cooled, super thermal conducting leading edges subjected to intense aerodynamic heating.

Deuchamphai, Pramote↗
NASA NTRS1988

Flow-thermal-structural study of aerodynamically heated leading edges

A finite element approach for integrated fluid-thermal-structural analysis of aerodynamically heated leading edges is presented. The Navier-Stokes equations for high speed compressible flow, the energy equation, and the quasi-static equilibrium equations for the leading edge are solved using a single finite element approach in one integrated, vectorized computer program called LIFTS. The fluid-thermal-structural coupling is studied for Mach 6.47 flow over a 3-inch diameter cylinder for which the flow behavior and the aerothermal loads are calibrated by experimental data. Issues of the thermal-structural response are studied for hydrogen cooled, super thermal conducting leading edges subjected to intense aerodynamic heating.

Dechaumphai, Pramote↗
NASA NTRS1999

Frequency-Dependent Viscosity of Xenon Near the Critical Point

We used a novel, overdamped oscillator aboard the Space Shuttle to measure the viscosity eta of xenon near its critical density rho(sub c), and temperature T(sub c). In microgravity, useful data were obtained within 0.1 mK of T(sub c), corresponding to a reduced temperature t = (T -T(sub c))/T(sub c) = 3 x 10(exp -7). The data extend two decades closer to T(sub c) than the best ground measurements, and they directly reveal the expected power-law behavior eta proportional to t(sup -(nu)z(sub eta)). Here nu is the correlation length exponent, and our result for the small viscosity exponent is z(sub eta) = 0.0690 +/- 0.0006. (All uncertainties are one standard uncertainty.) Our value for z(sub eta) depends only weakly on the form of the viscosity crossover function, and it agrees with the value 0.067 +/- 0.002 obtained from a recent two-loop perturbation expansion. The measurements spanned the frequency range 2 Hz less than or equal to f less than or equal to 12 Hz and revealed viscoelasticity when t less than or equal to 10(exp -1), further from T(sub c) than predicted. The viscoelasticity scales as Af(tau), where tau is the fluctuation-decay time. The fitted value of the viscoelastic time-scale parameter A is 2.0 +/- 0.3 times the result of a one-loop perturbation calculation. Near T(sub c), the xenon's calculated time constant for thermal diffusion exceeded days. Nevertheless, the viscosity results were independent of the xenon's temperature history, indicating that the density was kept near rho(sub c), by judicious choices of the temperature vs. time program. Deliberately bad choices led to large density inhomogeneities. At t greater than 10(exp -5), the xenon approached equilibrium much faster than expected, suggesting that convection driven by microgravity and by electric fields slowly stirred the sample.

Berg, Robert F.↗
NASA NTRS1992

Dynamics and energetics of the solar corona

The primary objective of this research program is to improve our understanding of the dynamics and energetics of the solar corona both in the quiescent dynamic equilibrium state when coronal structure is dominated by the equatorial streamer belt and in the eruptive state when coronal plasma is ejected into the interplanetary medium. Numerical solutions of the time-dependent magnetohydrodynamic (MHD) equations and comparisons of the computed results with observations form the core of the approach to achieving this objective. Some of the specific topics that have been studied are: (1) quiescent coronal streamers in an atmosphere dominated by a dipole magnetic field at large radii, (2) the formation of coronal mass ejections (CMEs) in quiescent streamers due to the emergence of new magnetic flux and due to photospheric shear motion, (3) MHD shock formation near the leading edge of CMEs, (4) coronal magnetic arcade eruption as a result of applied photospheric shear motion, and (5) the three-dimensional structure of CMEs.

Steinolfson, R. S.↗
NASA NTRS1992

Heat transfer measurements and CFD comparison of swept shock wave/boundary-layer interactions

An experimental research program providing basic knowledge and establishing new data on the heat transfer in swept shock wave/boundary-layer interactions is described. An equilibrium turbulent boundary-layer on a flat plate is subjected to impingement by swept planar shock waves generated by a sharp fin. Five different interactions with fin angles ranging from 10 to 20 deg at freestream Mach numbers of 3.0 and 4.0 produce a variety of interaction strengths from weak to very strong. A foil heater generates a uniform heat flux over the flat plate surface and miniature thin-film-resistance sensors mounted on it are used to measure the local surface temperature. The heat convection equation is then solved for the heat transfer distribution within an interaction, yielding a total uncertainty of about +/- 10 percent. These experimental data are compared with the results of numerical Navier-Stokes solutions which employ a kappa-epsilon turbulence model. Finally, a simplified form of the peak heat transfer correlation for fin interactions is suggested.

Lee, Y.↗
NASA NTRS1986

Assessment of two-temperature kinetic model for dissociating and weakly-ionizing nitrogen

The validity of the author's two-temperature, chemical/kinetic model which the author has recently improved is assessed by comparing the calculated results with the existing experimental data for nitrogen in the dissociating and weakly ionizing regime produced behind a normal shock wave. The computer program Shock Tube Radiation Program (STRAP) based on the two-temperature model is used in calculating the flow properties behind the shock wave and the Nonequilibrium Air Radiation (NEQAIR) program, in determining the radiative characteristics of the flow. Both programs were developed earlier. Comparison is made between the calculated and the existing shock tube data on (1) spectra in the equilibrium region, (2) rotational temperature of the N2(+) B state, (3) vibrational temperature of the N2(+) B state, (4) electronic excitation temperature of the N2 B state, (5) the shape of time-variation of radiation intensities, (6) the times to reach the peak in radiation intensity and equilibrium, and (7) the ratio of nonequilibrium to equilibrium radiative heat fluxes. Good agreement is seen between the experimental data and the present calculation except for the vibrational temperature. A possible reason for the discrepancy is given.

Park, C.↗