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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 397 records · Page 22

All sources NASA NTRS DOE OSTI

NASA NTRS1976

Supersonic Flow of Chemically Reacting Gas-Particle Mixtures. Volume 2: RAMP - A Computer Code for Analysis of Chemically Reacting Gas-Particle Flows

A computer program written in conjunction with the numerical solution of the flow of chemically reacting gas-particle mixtures was documented. The solution to the set of governing equations was obtained by utilizing the method of characteristics. The equations cast in characteristic form were shown to be formally the same for ideal, frozen, chemical equilibrium and chemical non-equilibrium reacting gas mixtures. The characteristic directions for the gas-particle system are found to be the conventional gas Mach lines, the gas streamlines and the particle streamlines. The basic mesh construction for the flow solution is along streamlines and normals to the streamlines for axisymmetric or two-dimensional flow. The analysis gives detailed information of the supersonic flow and provides for a continuous solution of the nozzle and exhaust plume flow fields. Boundary conditions for the flow solution are either the nozzle wall or the exhaust plume boundary.

Penny, M. M.↗

NASA NTRS1978

Structural Dynamics, Stability, and Control of Helicopters

The dynamic synthesis of gyroscopic structures consisting of point-connected substructures is investigated. The objective is to develop a mathematical model capable of an adequate simulation of the modal characteristics of a helicopter using a minimum number of degrees of freedom. The basic approach is to regard the helicopter structure as an assemblage of flexible substructures. The variational equations for the perturbed motion about certain equilibrium solutions are derived. The discretized variational equations can be conveniently exhibited in matrix form, and a great deal of information about the system modal characteristics can be extracted from the coefficient matrices. The derivation of the variational equations requires a monumental amount of algebraic operations. To automate this task a symbolic manipulation program on a digital computer is developed.

Meirovitch, L.↗

NASA NTRS1990

Free-vibration analysis of three-dimensional solids by BEM

This paper discusses the calculation of natural frequencies and mode shapes of three-dimensional solids, using the boundary element method. The metod developed is based on the use of particular integrals of the elastic equilibrium equations, and employs only real arithmetic. It represents an extension to three dimensions of work previously carried out for two-dimensional problems by Ahmad and Banerjee (1986) and Nardini and Brebbia (1982). The method has been incorporated, for multiregion analysis, in the BEST3D (Boundary Element Solution Technology, Three-Dimensional) computer program. Results of the boundary element calculations are compared with analytical, finite element, and experimental results.

Wilson, R. B.↗

DOE OSTI2021

Crystalline Silicotitanate (CST) Ion Exchange Media Performance Evaluations to Support TSCR DSA IX Media Equilibrium Contacts

The primary objective of this work is to calculate the maximum loading expected on the Hanford Tank Side Cesium Removal (TSCR) ion exchange columns. A key consideration in the design of the columns is the amount of 137 Cs that loads onto the Crystalline Silicotitanate (CST) and the heat generated by the loaded column during storage. Per request of Washington River Protection Solutions (WRPS), Savannah River National Laboratory (SRNL) has utilized ZAM, a computer program developed by the research group of Professor Rayford G. Anthony of Texas A&M University, to predict the cesium loading on the CST for a variety of waste compositions expected to be processed by TSCR. The study evaluated cesium loadings for the following waste compositions: 1. Seventeen DFLAW campaign batches to cover projected supernate composition ranges within which TSCR may be expected to operate within the first ten years, 2. Hanford tank AP-105 and AP-107 waste solutions that will be processed by the TSCR system.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2019

Thermal-hydraulic analyses of MIT reactor LEU transition cycles

The Massachusetts Institute of Technology Reactor (MITR) is a research reactor in Cambridge, Massachusetts designed primarily for experiments using neutron beam and in-core irradiation facilities. In the framework of non-proliferation policy, research and test reactors have started a program to convert highly enriched uranium (HEU) fuel to low-enriched uranium (LEU) fuel. A new type of LEU fuel based on a high-density alloy of uranium and molybdenum (U-10Mo) is expected to allow conversion. A transition core plan, from 22 fresh LEU fuel elements to 24 in an equilibrium state, is proposed. This study evaluates the thermal-hydraulic safety margins of each transition cycle and state (i.e. beginning, middle and end). The STAT7 and RELAP5 codes are used in the study. STAT7 provides an integrated platform for statistical propagation of uncertainties in determining steady-state thermal-hydraulic margins. A RELAP5 model was created to verify MITR's STAT7 model, and provide transient analyses. The steady-state analyses determined that, during the entire transition from a fresh to equilibrium core, the Limiting System Safety Settings (LSSS) has a 10.2% margin to the minimum power at which onset of nucleate boiling (ONB) is precluded at 3σ confidence. The loss of flow (LOF) transient, which is the most limiting anticipated transient for the MITR, was simulated. With a conservative initial state of the reactor (operating at μ + 3σ power level and LSSS mass flow rate), the maximum calculated centerline fuel temperature was 115 degrees C, which is a significantly lower value than the blistering limit of 350 degrees C. The cladding wall temperatures did not exceed the ONB temperature throughout the transient. Therefore, during LOF, no nucleate boiling is expected, precluding any critical heat flux trigger. Overall, this study indicates that the proposed LEU transition core specifications have significant margin to thermal-hydraulic limits, during steady-state operation and the LOF transient.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2022

Atmospheric CO 2 estimates for the Miocene to Pleistocene based on foraminiferal $δ^{11}$B at Ocean Drilling Program Sites 806 and 807 in the Western Equatorial Pacific

Constraints on the evolution of atmospheric CO 2 levels throughout Earth's history are foundational to our understanding of past variations in climate. Despite considerable effort, records vary in their temporal and spatial coverage and estimates of past CO 2 levels do not always converge, and therefore new records and proxies are valuable. Here we reconstruct atmospheric CO 2 values across major climate transitions over the past 16 million years using the boron isotopic composition (δ 11 B) of planktic foraminifera from 89 samples obtained from two sites in the West Pacific Warm Pool, Ocean Drilling Program (ODP) Sites 806 and 807, measured using high-precision multi-collector inductively coupled plasma mass spectrometry. We compare our results to published data from ODP Site 872, also in the Western Equatorial Pacific, that goes back to 22 million years ago. These sites are in a region that today is near equilibrium with the atmosphere and are thought to have been in equilibrium with the atmosphere for the interval studied. We show that δ 11 B data from this region are consistent with other boron-based studies. The data show evidence for elevated pCO 2 during the Middle Miocene and Early to Middle Pliocene, and reductions in pCO 2 of ~200 ppm during the Middle Miocene Climate Transition, ~250 ppm during Pliocene Glacial Intensification and ~50 ppm during the Mid-Pleistocene Climate Transition. During the Mid-Pleistocene Transition there is a minimum pCO 2 at marine isotopic stage (MIS) 30. Our results are consistent with a coupling between pCO 2 , temperature and ice sheet expansion from the Miocene to the late Quaternary.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1972

Study of effects of design details on structural response to acoustic excitation

The normal mode vibration characteristics of one-dimensional and two-dimensional panel arrays were investigated analytically and experimentally. The finite element displacement method was used to formulate the structural models. The structural models include a stiffness and consistent mass matrix for thin-walled open-section beams not previously reported, and a modification of a rectangular plate bending element to include a fundamental interior mode for the element as a generalized coordinate. Provision for adding lumped masses to represent accelerometers is also included. For the one-dimensional panel, the normal mode stress resultants are obtained by integrating the equilibrium equations directly. For the two-dimensional panel arrays, the solution for the stress resultants in the cover sheet was attempted by introducing stress functions dependent upon the displacement field of the plate element. This approach, although incomplete, is a natural extension of the one-dimensional analysis for plate elements. The exerimental program which was used to verify the analysis is described and experimental results are compared with the analysis.

Rudder, F. F., Jr.↗

NASA NTRS1975

Importance of helicopter dynamics to the mathematical model of the helicopter

A mathematical model of the helicopter requires appropriate representation of the constituent elements of rotor dynamics. General-purpose programs that model a variety of configurations for a broad range of operating conditions result in varying and incompatible levels of sophistication. Analysis of specific dynamic problems facilitates the identification of configuration parameters which determine system behavior. For the present analysis, the nonlinear equations of a torsionally rigid hingeless rotor are linearized about an equilibrium condition to determine flap-lag stability characteristics in hover. A collocation method was used to obtain the coupled natural frequencies and modes. These modes allow exact treatment of the effect of elastic coupling which more than compensates for the destabilizing inertial coupling. The sensitivity of damping to the number of modes was found to be small, and reasonable accuracy was obtained the first flapwise and edgewise coupled modes. The range of destabilizing precone was found to be small.

White, W. F., Jr.↗

NASA NTRS1981

Forebody and base region real gas flow in severe planetary entry by a factored implicit numerical method. II - Equilibrium reactive gas

The factored-implicit, gas-dynamic algorithm has been adapted to the numerical simulation of equilibrium reactive flows. Changes required in the perfect gas version of the algorithm are developed, and the method of coupling gas-dynamic and chemistry variables is discussed. A flow-field solution that approximates a Jovian entry case was obtained by this method and compared with the same solution obtained by HYVIS, a computer program much used for the study of planetary entry. Comparison of surface pressure distribution and stagnation line shock-layer profiles indicates that the two solutions agree well.

Davy, W. C.↗

NASA NTRS1981

Numerical experiments on galaxy formation. I - Introduction and first results

The large, fully three-dimensional n-body programs designed for numerical experiments on the dynamics of galaxies have been used for investigations of protogalactic collapses. The dynamical aspects of galaxy formation dominate in this approach. The present experiments start from a large rotating gaseous mass that is made up of many gas clouds. The large mass collapses because it is out of equilibrium. Stars form as the collapse proceeds. The problem is formulated, limiting-case examples are described, and results from simple 'first cut' models are given. Results with no star formation are contrasted with results for very rapid star formation. All models with gas formed thin disklike structures; some were circular, while others were oval. The principal result from the simple models is that stars continue to form rapidly even after the gas has been depleted to unacceptably low levels (1% of the galactic mass in gas).

Miller, R. H.↗

DOE OSTI2020

Phosphoproteome Response to Dithiothreitol Reveals Unique Versus Shared Features of Saccharomyces cerevisiae Stress Responses

To cope with sudden changes in the external environment, the budding yeast Saccharomyces cerevisiae orchestrates a multifaceted response that spans many levels of physiology. Several studies have interrogated the transcriptome response to endoplasmic reticulum (ER) stress and the role of regulators such as the Ire1 kinase and Hac1 transcription factors. However, less is known about responses to ER stress at other levels of physiology. Here, we used quantitative phosphoproteomics and computational network inference to uncover the yeast phosphoproteome response to the reducing agent dithiothreitol (DTT) and the upstream signaling network that controls it. We profiled wild-type cells and mutants lacking IRE1 or MAPK kinases MKK1 and MKK2, before and at various times after DTT treatment. In addition to revealing downstream targets of these kinases, our inference approach predicted new regulators in the DTT response, including cell-cycle regulator Cdc28 and osmotic-response kinase Rck2, which we validated computationally. Our results also revealed similarities and surprising differences in responses to different stress conditions, especially in the response of protein kinase A targets. Furthermore, these results have implications for the breadth of signaling programs that can give rise to common stress response signatures.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS2000

Stress Counter-Response Training Via Physiological Self-Regulation During Flight Simulation

This study provided the first evaluation of a new training concept and technology aimed at training pilots to maintain physiological equilibrium during circumstances in an airplane cockpit. Thirty healthy subjects (16 males and 14 females) between the ages of 18 and 35 were randomized into two study groups, A and B. Subjects participated individually in a sequence of four study sessions. In the first visit, subjects were taught to operate a desktop fighter jet flight simulation program. In the three sessions that followed, subjects in group A were trained to minimize their autonomic deviation from baseline values while operating the desktop flight simulation. This was done by making their skin conductance and hand temperature deviations from baseline impair the functionality of the aircraft controls. Subjects also received auditory and visual cues about their autonomic deviation, and were instructed to keep these within pre-set limits to retain full control of the aircraft. Subjects in group B were subjected to periods of impaired aircraft functionality independent of their physiologic activity, and thus served as a control group. No statistically significant group differences were found in the flight performance scores from the three training sessions, and post-training flight performance scores of the two groups were not different. We conclude that this study did not provide clear support for this training methodology in optimizing pilot performance. However, a number of shortcomings in the current status of this training methodology may account for the lack of demonstrable training benefit to the experimental group. Suggested future modifications for research on this training methodology include: Limiting the amount of instrument impairment resulting from physiological deviations; conducting a greater number of physiological training sessions per subject; using pre-post training performance tests which invoke a greater amount of stress in subjects; and developing a more detailed performance scoring system.

Palsson, Olafur S.↗

NASA NTRS1998

Bulk-Flow Analysis of Hybrid Thrust Bearings for Advanced Cryogenic Turbopumps

A bulk-flow analysis and computer program for prediction of the static load performance and dynamic force coefficients of angled injection, orifice-compensated hydrostatic/hydrodynamic thrust bearings have been completed. The product of the research is an efficient computational tool for the design of high-speed thrust bearings for cryogenic fluid turbopumps. The study addresses the needs of a growing technology that requires of reliable fluid film bearings to provide the maximum operating life with optimum controllable rotordynamic characteristics at the lowest cost. The motion of a cryogenic fluid on the thin film lands of a thrust bearing is governed by a set of bulk-flow mass and momentum conservation and energy transport equations. Mass flow conservation and a simple model for momentum transport within the hydrostatic bearing recesses are also accounted for. The bulk-flow model includes flow turbulence with fluid inertia advection, Coriolis and centrifugal acceleration effects on the bearing recesses and film lands. The cryogenic fluid properties are obtained from realistic thermophysical equations of state. Turbulent bulk-flow shear parameters are based on Hirs' model with Moody's friction factor equations allowing a simple simulation for machined bearing surface roughness. A perturbation analysis leads to zeroth-order nonlinear equations governing the fluid flow for the thrust bearing operating at a static equilibrium position, and first-order linear equations describing the perturbed fluid flow for small amplitude shaft motions in the axial direction. Numerical solution to the zeroth-order flow field equations renders the bearing flow rate, thrust load, drag torque and power dissipation. Solution to the first-order equations determines the axial stiffness, damping and inertia force coefficients. The computational method uses well established algorithms and generic subprograms available from prior developments. The Fortran9O computer program hydrothrust runs on a Windows 95/NT personal computer. The program, help files and examples are licensed by Texas A&M University Technology License Office. The study of the static and dynamic performance of two hydrostatic/hydrodynamic bearings demonstrates the importance of centrifugal and advection fluid inertia effects for operation at high rotational speeds. The first example considers a conceptual hydrostatic thrust bearing for an advanced liquid hydrogen turbopump operating at 170,000 rpm. The large axial stiffness and damping coefficients of the bearing should provide accurate control and axial positioning of the turbopump and also allow for unshrouded impellers, therefore increasing the overall pump efficiency. The second bearing uses a refrigerant R134a, and its application in oil-free air conditioning compressors is of great technological importance and commercial value. The computed predictions reveal that the LH2 bearing load capacity and flow rate increase with the recess pressure (i.e. increasing orifice diameters). The bearing axial stiffness has a maximum for a recess pressure rati of approx. 0.55. while the axial damping coefficient decreases as the recess pressure ratio increases. The computer results from three flow models are compared. These models are a) inertialess, b) fluid inertia at recess edges only, and c) full fluid inertia at both recess edges and film lands. The full inertia model shows the lowest flow rates, axial load capacity and stiffness coefficient but on the other hand renders the largest damping coefficients and inertia coefficients. The most important findings are related to the reduction of the outflow through the inner radius and the appearance of subambient pressures. The performance of the refrigerant hybrid thrust bearing is evaluated at two operating speeds and pressure drops. The computed results are presented in dimensionless form to evidence consistent trends in the bearing performance characteristics. As the applied axial load increases, the bearing film thickness and flow rate decrease while the recess pressure increases. The axial stiffness coefficient shows a maximum for a certain intermediate load while the damping coefficient steadily increases. The computed results evidence the paramount of centrifugal fluid inertia at low recess pressures (i.e. low loads), and where there is actually an inflow through the bearing inner diameter, accompanied by subambient pressures just downstream of the bearing recess edge. These results are solely due to centrifugal fluid inertia and advection transport effects. Recommendations include the extension of the computer program to handle flexure pivot tilting pad hybrid bearings and the ability to calculate moment coefficients for shaft angular misalignments.

SanAndres, Luis↗

NASA NTRS1994

Aerothermodynamics of Pyrolizing Surfaces in Hypersonic Rarefied Flows

Direct simulation Monte Carlo (DSMC) calculations of rarefied flows about entry bodies typically employ a fixed surface temperature or a radiative-equilibrium energy balance to compute that temperature. Such boundary conditions neglect any effects of heat capacitance and heat conduction in the spacecraft heat shield and, therefore, provide an upper bound for the surface temperature. Such calculations also neglect pyrolysis from the heat shield which can be significant for a high-energy incident flow at very low densities. Accurate prediction of both heating and aerodynamic forces requires including pyrolysis and surface heat transfer in the models for gas-surface interaction employed in DSMC methods. Although these physical models have long appeared in various continuum flow calculation codes, they have only recently appeared in DSMC codes which are required to simulate rarefied flows during entry at high altitudes. In the current implementation, routines from the widely distributed Charring Material Thermal Response and Ablation (CMA) program are coupled into a DSMC code to calculate the one-dimensional heat transfer into the carbon phenolic heat shield at each point on a vehicle surface. Temperature-dependent material properties, surface re-radiation, and in-depth pyrolysis were included in the calculation, but surface ablation was neglected. Sample calculations for entry of the Galileo probe into the atmosphere of Jupiter demonstrate that including pyrolysis in the model leads to significant differences in predicted aerodynamics. Granted, the drag coefficient does not depend strongly on the surface temperature which can itself be significantly below the radiative equilibrium value during entry. However, the surface mass flux due to pyrolysis of the material is significant once the probe drops to altitudes characterized by transition flow. This leads to a noticeable increase in drag and a decrease in heating compared to a body without pyrolysis.

Haas, Brian L.↗

NASA NTRS1989

Ignition and combustion of metallized propellants

An experimental program was conducted to investigate the ignition and combustion properties of a gelled, metallized fuel, A1/JP20, burned with gaseous oxygen, and to compare the results to the ignition and combustion properties of neat RP1. Both a Cstar efficiency and a specific impulse efficiency were calculated by comparing experimental values with those from a one-dimensional chemical equilibrium computer code. For most operating conditions, the efficiencies of the A1/JP10 were only 2 to 6 percent less than those of the neat RP1.

Galecki, Diane L.↗

NASA NTRS1995

Program Helps To Determine Chemical-Reaction Mechanisms

General Chemical Kinetics and Sensitivity Analysis (LSENS) computer code developed for use in solving complex, homogeneous, gas-phase, chemical-kinetics problems. Provides for efficient and accurate chemical-kinetics computations and provides for sensitivity analysis for variety of problems, including problems involving honisothermal conditions. Incorporates mathematical models for static system, steady one-dimensional inviscid flow, reaction behind incident shock wave (with boundary-layer correction), and perfectly stirred reactor. Computations of equilibrium properties performed for following assigned states: enthalpy and pressure, temperature and pressure, internal energy and volume, and temperature and volume. Written in FORTRAN 77 with exception of NAMELIST extensions used for input.

Bittker, D. A.↗

NASA NTRS1980

Spherical roller bearing analysis. SKF computer program SPHERBEAN. Volume 1: Analysis

The models and associated mathematics used within the SPHERBEAN computer program for prediction of the thermomechanical performance characteristics of high speed lubricated double row spherical roller bearings are presented. The analysis allows six degrees of freedom for each roller and three for each half of an optionally split cage. Roller skew, free lubricant, inertial loads, appropriate elastic and friction forces, and flexible outer ring are considered. Roller quasidynamic equilibrium is calculated for a bearing with up to 30 rollers per row, and distinct roller and flange geometries are specifiable. The user is referred to the material contained here for formulation assumptions and algorithm detail.

Kleckner, R. J.↗

NASA NTRS1999

Development of Curved-Plate Elements for the Exact Buckling Analysis of Composite Plate Assemblies Including Transverse-Shear Effects

The analytical formulation of curved-plate non-linear equilibrium equations including transverse-shear-deformation effects is presented. A unified set of non-linear strains that contains terms from both physical and tensorial strain measures is used. Linearized, perturbed equilibrium equations (stability equations) that describe the response of the plate just after buckling occurs are derived. These equations are then modified to allow the plate reference surface to be located a distance z(sub c) from the centroidal surface. The implementation of the new theory into the VICONOPT exact buckling and vibration analysis and optimum design computer program is described. The terms of the plate stiffness matrix using both classical plate theory (CPT) and first-order shear-deformation plate theory (SDPT) are presented. The effects of in-plane transverse and in-plane shear loads are included in the in-plane stability equations. Numerical results for several example problems with different loading states are presented. Comparisons of analyses using both physical and tensorial strain measures as well as CPT and SDPT are made. The computational effort required by the new analysis is compared to that of the analysis currently in the VICONOPT program. The effects of including terms related to in-plane transverse and in-plane shear loadings in the in-plane stability equations are also examined. Finally, results of a design-optimization study of two different cylindrical shells subject to uniform axial compression are presented.

McGowan, David M.↗