Engineering Papers⌕ Search

SEARCH · Engineering Papers

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

How numerical simulations helped to achieve breakeven on the NIF

The inertial confinement fusion program relies upon detailed simulations with inertial confinement fusion (ICF) codes to design targets and to interpret the experimental results. These simulations treat as much physics from essential principles as is practical, including laser deposition, cross beam energy transfer, x-ray production and transport, nonlocal thermal equilibrium kinetics, thermal transport, hydrodynamic instabilities, thermonuclear burn, and transport of reaction products. Improvements in radiation hydrodynamic code capabilities and vast increases in computing power have enabled more realistic, accurate 3D simulations that treat all known asymmetry sources. We describe how numerical simulations helped to guide the program, assess the impediments to breakeven, and optimize every aspect of target design. A preshot simulation of the first National Ignition Facility experiment that surpassed breakeven predicted an increased yield that matches the experimental result, within the preshot predicted uncertainty, with a target gain of 1.5. We will cover the key developments in Lawrence Livermore National Laboratory ICF codes that enabled these simulations and give specific examples of how they helped to guide the program.

Marinak, M. M. (ORCID:0009000331127696)↗

A vectorized algorithm for 3D dynamics of a tethered satellite

Equations of motion characterizing the three dimensional motion of a tethered satellite during the retrieval phase are studied. The mathematical model involves an arbitrary number of point masses connected by weightless cords. Motion occurs in a gravity gradient field. The formulation presented accounts for general functions describing support point motion, rate of tether retrieval, and arbitrary forces applied to the point masses. The matrix oriented program language MATLAB is used to produce an efficient vectorized formulation for computing natural frequencies and mode shapes for small oscillations about the static equilibrium configuration; and for integrating the nonlinear differential equations governing large amplitude motions. An example of time response pertaining to the skip rope effect is investigated.

Wilson, Howard B.↗

Data-based and secure switched cyber–physical systems

In this work, we develop a completely model-free moving target defense framework for the detection and mitigation of sensor and/or actuator attacks in cyber–physical systems with dynamics that evolve in discrete-time. We incorporate an intrusion detection mechanism based on an approximate dynamic programming technique that learns the policies for optimal regulation and optimal tracking while simultaneously defending against actuator and sensor attacks in a model-free fashion. Switching rules are leveraged to force proactive and reactive defense mechanisms as well as, guarantee the stability of the equilibrium point. Finally, as a case study, we apply the proposed moving target defense framework to a DC–DC converter that is used in electric vehicles.

42 ENGINEERING↗

Multivariable control altitude demonstration on the F100 turbofan engine

The F100 Multivariable control synthesis (MVCS) program, was aimed at demonstrating the benefits of LGR synthesis theory in the design of a multivariable engine control system for operation throughout the flight envelope. The advantages of such procedures include: (1) enhanced performance from cross-coupled controls, (2) maximum use of engine variable geometry, and (3) a systematic design procedure that can be applied efficiently to new engine systems. The control system designed, under the MVCS program, for the Pratt & Whitney F100 turbofan engine is described. Basic components of the control include: (1) a reference value generator for deriving a desired equilibrium state and an approximate control vector, (2) a transition model to produce compatible reference point trajectories during gross transients, (3) gain schedules for producing feedback terms appropriate to the flight condition, and (4) integral switching logic to produce acceptable steady-state performance without engine operating limit exceedance.

Lehtinen, B.↗

Potential theory of radiation

A theoretical method is being developed by which the structure of a radiation field can be predicted by a radiation potential theory, similar to a classical potential theory. The introduction of a scalar potential is justified on the grounds that the spectral intensity vector is irrotational. The vector is also solenoidal in the limits of a radiation field in complete radiative equilibrium or in a vacuum. This method provides an exact, elliptic type equation that will upgrade the accuracy and the efficiency of the current CFD programs required for the prediction of radiation and flow fields. A number of interesting results emerge from the present study. First, a steady state radiation field exhibits an optically modulated inverse square law distribution character. Secondly, the unsteady radiation field is structured with two conjugate scalar potentials. Each is governed by a Klein-Gordon equation with a frictional force and a restoring force. This steady potential field structure and the propagation of radiation potentials are consistent with the well known results of classical electromagnetic theory. The extension of the radiation potential theory for spray combustion and hypersonic flow is also recommended.

Chiu, Huei-Huang↗

Computer code for single-point thermodynamic analysis of hydrogen/oxygen expander-cycle rocket engines

This analysis and this computer code apply to full, split, and dual expander cycles. Heat regeneration from the turbine exhaust to the pump exhaust is allowed. The combustion process is modeled as one of chemical equilibrium in an infinite-area or a finite-area combustor. Gas composition in the nozzle may be either equilibrium or frozen during expansion. This report, which serves as a users guide for the computer code, describes the system, the analysis methodology, and the program input and output. Sample calculations are included to show effects of key variables such as nozzle area ratio and oxidizer-to-fuel mass ratio.

Glassman, Arthur J.↗

Rarefied-continuum gas dynamics transition for SUMS project

This program is to develop an analytic method for reducing SUMS data for the determination of the undisturbed atmosphere conditions ahead of the shuttle along its descending trajectory. It is divided into an internal flow problem, an external flow problem and their matching conditions. Since the existing method of Direct Simulation Monte Carlo (DSMC) failed completely for the internal flow problem, the emphasis is on the internal flow of a highly non-equilibrium, rarefied air through a short tube of a diameter much less than the gaseous mean free path. A two fluid model analysis of this internal flow problem has been developed and studied with typical results illustrated. A computer program for such an analysis and a technical paper published in Lecture Notes in Physics No. 323 (1989) are included as Appendices 3 and 4. A proposal for in situ determination of the surface accommodation coefficients sigma sub t and sigma e is included in Appendix 5 because of their importance in quantitative data reduction. A two fluid formulation for the external flow problem is included as Appendix 6 and a review article for AIAA on Hypersonic propulsion, much dependent on ambient atmospheric density, is also included as Appendix 7.

Cheng, Sin-I↗

Atomic Spectroscopic Databases at NIST

We describe recent work at NIST to develop and maintain databases for spectra, transition probabilities, and energy levels of atoms that are astrophysically important. Our programs to critically compile these data as well as to develop a new database to compare plasma calculations for atoms that are not in local thermodynamic equilibrium are also summarized.

Reader, J.↗

A Fuzzy Approach of the Competition on the Air Transport Market

The aim of this communication is to study with a new scope the conditions of the equilibrium in an air transport market where two competitive airlines are operating. Each airline is supposed to adopt a strategy maximizing its profit while its estimation of the demand has a fuzzy nature. This leads each company to optimize a program of its proposed services (frequency of the flights and ticket prices) characterized by some fuzzy parameters. The case of monopoly is being taken as a benchmark. Classical convex optimization can be used to solve this decision problem. This approach provides the airline with a new decision tool where uncertainty can be taken into account explicitly. The confrontation of the strategies of the companies, in the ease of duopoly, leads to the definition of a fuzzy equilibrium. This concept of fuzzy equilibrium is more general and can be applied to several other domains. The formulation of the optimization problem and the methodological consideration adopted for its resolution are presented in their general theoretical aspect. In the case of air transportation, where the conditions of management of operations are critical, this approach should offer to the manager elements needed to the consolidation of its decisions depending on the circumstances (ordinary, exceptional events,..) and to be prepared to face all possibilities. Keywords: air transportation, competition equilibrium, convex optimization , fuzzy modeling,

Charfeddine, Souhir↗

Theory Manual for the Fuel Cycle Analysis Code REBUS

REBUS is a system of programs designed for the fuel cycle analysis of fast reactors and accelerator-driven systems. The first version of REBUS was developed in 1970 for calculating the equilibrium operating conditions of liquid metal fast breeder reactors. During the past five decades, various extensions and refinements have been made in the simulation of material motions in fuel cycle facilities and the solution of the neutronics and isotopic transmutation equations that govern the reactor behavior. The latest version complies with the standard code practices and interface dataset specifications of the committee on computer code coordination (CCCC). This report provides a detailed description of the physical processes and their mathematical models and solution methods of the current production version of REBUS. The models of the facilities and physical processes involved in fast reactor fuel cycle are first discussed. Then the mathematical formulation of these physical models is presented, followed by detailed discussions of solution methods of the resulting set of equations. Finally the output edit quantities are discussed along with their definitions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Solute partitioning under continuous cooling conditions as a cooling rate indicator

A model of solute partitioning in a finite body under conditions of continuous cooling is developed for the determination of cooling rates from concentration profile data, and applied to the partitioning of zirconium between ilmenite and ulvospinel in the Apollo 15 Elbow Crater rocks. Partitioning in a layered composite solid is described numerically in terms of concentration profiles and diffusion coefficients which are functions of time and temperature, respectively; a program based on the model can be used to calculate concentration profiles for various assumed cooling rates given the diffusion coefficients in the two phases and the equilibrium partitioning ratio over a range of temperatures. In the case of the Elbow Rock gabbros, the cooling rates are calculated from measured concentration ratios 10 microns from the interphase boundaries under the assumptions of uniform and equilibrium initial conditions at various starting temperatures. It is shown that the specimens could not have had uniform concentrations profiles at the previously suggested initial temperature of 1350 K. It is concluded that even under conditions where the initial temperature, grain sizes and solute diffusion coefficients are not well characterized, the model can be used to estimate the cooling rate of a grain assemblage to within an order of magnitude.

Onorato, P. I. K.↗

A review of turbulent-boundary-layer heat transfer research at Stanford, 1958-1983

For the past 25 years, there has existed in the Thermosciences Laboratory of the Mechanical Engineering Department of Stanford University a research program, primarily experimental, concerned with heat transfer through turbulent boundary layers. In the early phases of the program, the topics considered were the simple zero-pressure-gradient turbulent boundary layer with constant and with varying surface temperature, and the accelerated boundary layer. Later equilibrium boundary layers were considered along with factors affecting the boundary layer, taking into account transpired flows, flows with axial pressure gradients, transpiration, acceleration, deceleration, roughness, full-coverage film cooling, surface curvature, free convection, and mixed convection. A description is provided of the apparatus and techniques used, giving attention to the smooth plate rig, the rough plate rig, the full-coverage film cooling rig, the curvature rig, the concave wall rig, the mixed convection tunnel, and aspects of data reduction and uncertainty analysis.

Moffat, R. J.↗

Adsorption and recombination of hydrogen atoms on a model graphite surface

The adsorption and recombination of atomic hydrogen on a model graphite grain have been examined in a series of calculations in which a modified, iterative, extended Hueckel program was used. The hydrogen atom is found to be chemisorbed at a site with a zero-point binding energy of 0.7 eV and at an equilibrium distance of 2.25 A above the site. Despite a barrier of about 0.4 eV between adjacent sites, calculations suggest that at temperatures as low as 10 K, an H atom will tunnel through to adjacent sites in less than one nanosecond. However, a potential barrier to the recombination of two hydrogen atoms has been found which displays high sensitivity to the mutual arrangement of the two hydrogen atoms with respect to the graphite surface. Results show that at very low temperatures, recombinations can occur only by tunneling. Consistent with experiment, the region in which H2 begins to form exhibits a repulsive potential with respect to possible chemisorption of the incipient H2 entity.

Aronowitz, S.↗

Aeroelastic tailoring for oblique wing lateral trim

Composite material aeroelastic tailoring is presently explored as a means for the correction of the roll trim imbalance of oblique-wing aircraft configurations. The concept is demonstrated through the analysis of a realistic oblique wing by a static aeroelastic computational procedure encompassing the full potential transonic aerodynamic code FLO22 and a Ritz structural plate program that models the stiffness due to symmetrical-but-unbalanced composite wing skins. Results indicate that asymetric composite tailoring reduces the aileron deflection needed for roll equilibrium, and reduces control surface hinge moment and drag. Wing skin stresses are, however, very high.

Bohlmann, Jonathan D.↗

Development of AFE aerobrake aerothermodynamic data book

The computation method developed for the NASA Aeroassist Flight Experiment (AFE) data book generates a design reference for the AFE's aerothermodynamic environment using an optimized technology for a 4100-lb vehicle. This environment is defined by convective, radiative, and total heating rates, radiation equilibrium temperatures, and local surface pressures along the AFE pitch-plane and associated off-pitch planes. The Boundary Layer Integral Matrix Procedure is the major program code used in this analysis; a partially catalytic wall was assumed on the basis of measured recombination rates.

Ting, Paul C.↗

Regenerable biocide delivery unit, volume 2

Source code for programs dealing with the following topics are presented: (1) life cycle test stand-parametric test stand control (in BASIC); (2) simultaneous aqueous iodine equilibria-true equilibrium (in C); (3) simultaneous aqueous iodine equilibria-pseudo-equilibrium (in C); (4) pseudo-(fast)-equilibrium with iodide initially present (in C); (5) solution of simultaneous iodine rate expressions (Mathematica); (6) 2nd order kinetics of I2-formic acid in humidity condensate (Mathematica); (7) prototype RMCV onboard microcontroller (CAMBASIC); (8) prototype RAM data dump to PC (in BASIC); and (9) prototype real time data transfer to PC (in BASIC).

Atwater, James E.↗

High-Flux, High-Temperature Thermal Vacuum Qualification Testing of a Solar Receiver Aperture Shield

As part of the International Space Station (ISS) Phase 1 program, NASA Lewis Research Center (LERC) and the Russian Space Agency (RSA) teamed together to design, build and flight test the world's first orbital Solar Dynamic Power System (SDPS) on the Russian space station Mir. The Solar Dynamic Flight Demonstration (SDFD) program was to operate a nominal 2 kWe SDPS on Mir for a period up to 1-year starting in late 1997. Unfortunately, the SDFD mission was demanifested from the ISS phase 1 shuttle program in early 1996. However, substantial flight hardware and prototypical flight hardware was built including a heat receiver and aperture shield. The aperture shield comprises the front face of the cylindrical cavity heat receiver and is located at the focal plane of the solar concentrator. It is constructed of a stainless steel plate with a 1-m outside diameter, a 0.24-m inside diameter and covered with high-temperature, refractory metal Multi-Foil Insulation (MFI). The aperture shield must minimize heat loss from the receiver cavity, provide a stiff, high strength structure to accommodate shuttle launch loads and protect receiver structures from highly concentrated solar fluxes during concentrator off-pointing events. To satisfy Mir operational safety protocols, the aperture shield was required to accommodate direct impingement of the intensely concentrated solar image for a 1-hour period. To verify thermal-structural durability under the anticipated high-flux, high-temperature loading, an aperture shield test article was constructed and underwent a series of two tests in a large thermal vacuum chamber configured with a reflective, point-focus solar concentrator and a solar simulator. The test article was positioned near the focal plane and exposed to concentrated solar flux for a period of 1-hour. In the first test, a near equilibrium temperature of 1862 K was attained in the center of the shield hot spot. In the second test, with increased incident flux, a near equilibrium temperature of 2072 K was achieved. The aperture shield sustained no visible damage as a result of the exposures. This paper describes the aperture shield thermal-vacuum qualification test program including the test article, test facility, procedures, data collection, test success criteria, results and conclusions.

Kerslake, Thomas W.↗

Master equation study of hydrogen abstraction from HCHO by OH via a chemically activated intermediate

The abstraction reaction of hydrogen from formaldehyde by OH radical plays an important role in formaldehyde oxidation. The reaction involves a bimolecular association to form a chemically activated hydrogen-bonded reaction complex followed by a unimolecular reaction of the complex to generate the products. The reaction rate is usually considered to be pressure-independent by assuming equilibrium between the reactants and the complex. However, our nonequilibrium calculations based on the chemically significant eigenmode of the master equation, carried out with our recently developed TUMME program, indicate that the reaction complex makes the rate constant dependent on pressure at low temperatures (T < 200 K). The calculations include anharmonicity, variational effects, and multi-dimensional tunneling. We find that the reaction rate constant reaches a low-pressure limit at pressures below 10 Torr over the whole investigated temperature range (20–1800 K), which explains why the available low-temperature experiments, which are for pressures below 2 Torr, did not observe the pressure dependence. A new extension of the TUMME master-equation program is used to explore the time evolutions of the concentrations of the OH radical and the complex under pseudo-first-order conditions. The time-dependent evolution of the concentrations of the complex at a low temperature provide direct evidence for the stabilization of the reaction complex at high pressures, and it shows the negligible role of the stabilized reaction complex at low pressures. Finally, the picture that emerges is qualitatively consistent with our previous study of the reaction of methanol with OH in that the tunneling in the unimolecular step from the complex to the products affects the phenomenological reaction rate constants differently at high and low pressures and leads to a significant pressure effect.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗