SEARCH · Engineering Papers

Results for “computable general equilibrium”

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 109 records · Page 6

All sources NASA NTRS DOE OSTI

NASA NTRS1995

A kinematically driven anisotropic viscoelastic constitutive model applied to tires

Aircraft tires are composite structures manufactured with viscoelastic materials such as carbon black filled rubber and nylon cords. When loaded they experience large deflections and moderately large strains. Detailed structural models of tires require the use of either nonlinear shell or nonlinear three dimensional solid finite elements. Computational predictions of the dynamic response of tires must consider the composite viscoelastic material behavior in a realistic fashion. We describe a modification to a nonlinear anisotropic shell finite element so it can be used to model viscoelastic stresses during general deformations. The model is developed by introducing internal variables of the type used to model elastic strain energy. The internal variables are strains, curvatures, and transverse shear angles which are in a one-to-one correspondence with the generalized coordinates used to model the elastic strain energy for nonlinear response. A difference-relaxation equation is used to relate changes in the observable strain field to changes in the internal strain field. The internal stress state is introduced into the equilibrium equations by converting it to nodal loads associated with the element's displacement degrees of freedom. In this form the tangent matrix in the Newton-Raphson solution algorithm is not modified from its form for the nonlinear statics problem. Only the gradient vector is modified and the modification is not computationally costly. The existing finite element model for the Space Shuttle nose gear tire is used to provide examples of the algorithm. In the first example, the tire's rim is displaced at a constant rate up to a fixed value. In the second example, the tire's rim is enforced to follow a saw tooth load and unload curve to generate hysteresis loops.

Johnson, Arthur R.↗

NASA NTRS2003

Simulation of Inviscid Compressible Multi-Phase Flow with Condensation

Condensation of vapours in rapid expansions of compressible gases is investigated. In the case of high temperature gradients the condensation will start at conditions well away from thermodynamic equilibrium of the fluid. In those cases homogeneous condensation is dominant over heterogeneous condensation. The present work is concerned with development of a simulation tool for computation of high speed compressible flows with homogeneous condensation. The resulting ow solver should preferably be accurate and robust to be used for simulation of industrial flows in general geometries.

Kelleners, Philip↗

NASA NTRS1983

Predictions of entry heating for lower surface of shuttle orbiter

A broad base of thermocouple and phase change paint data was assembled and correlated to the nominal design 14414.1 and proposed STS-1 (first flight of the space transportation system) entry trajectories. Averaged data from phase change paint tests compared favorably with thermocouple data for predicting heating rates. Laminar and turbulent radiation equilibrium heating rates were computed on the lower surface of the Shuttle orbiter for both trajectories, and the lower surface center line results were compared both with aerodynamic heating design data and with flight values from the STS-1 and STS-2 trajectories. The peak laminar heating values from the aerodynamic heating design data book were generally 40 to 60 percent higher than the laminar estimates of this study, except at the 55 percent location of maximum span where the design data book values were less than 10 percent higher. Estimates of both laminar and turbulent heating rates compared favorably with flight data.

Edwards, C. L. W.↗

NASA NTRS1990

Symmetrization of conservation laws with entropy for high-temperature hypersonic computations

Results of Hughes, France, and Mallet are generalized to conservation law systems taking into account high-temperature effects. Symmetric forms of different equation sets are derived in terms of entropy variables. First, the case of a general divariant gas is studied; it can be specialized to the usual Navier-Stokes equations, as well as to situations where the gas is vibrationally excited, and undergoes equilibrium chemical reactions. The case of gas in thermochemical nonequilibrium is considered next. Transport phenomena, and in particular mass diffusion, are examined in the framework of symmetric advective-diffusive systems.

Chalot, F.↗

NASA NTRS1984

Computer simulation of surface and film processes

All the investigations which were performed employed in one way or another a computer simulation technique based on atomistic level considerations. In general, three types of simulation methods were used for modeling systems with discrete particles that interact via well defined potential functions: molecular dynamics (a general method for solving the classical equations of motion of a model system); Monte Carlo (the use of Markov chain ensemble averaging technique to model equilibrium properties of a system); and molecular statics (provides properties of a system at T = 0 K). The effects of three-body forces on the vibrational frequencies of triatomic cluster were investigated. The multilayer relaxation phenomena for low index planes of an fcc crystal was analyzed also as a function of the three-body interactions. Various surface properties for Si and SiC system were calculated. Results obtained from static simulation calculations for slip formation were presented. The more elaborate molecular dynamics calculations on the propagation of cracks in two-dimensional systems were outlined.

Tiller, W. A.↗

NASA NTRS1975

Calculations of the radiative and dynamical state of the Venus atmosphere

Results are reported for accurate multiple-scattering calculations to determine the solar-energy deposition profile in the atmosphere of Venus. It is found that most of the absorbed energy is deposited in the main cloud-layer region, located at altitudes above 35 km, and that the ground receives approximately 3% of the energy absorbed in toto by Venus. Using these results, vertical temperature profiles are computed under conditions of pure radiative equilibrium and radiative-convective equilibrium. Since the latter results satisfactorily match the temperature structure determined from various spacecraft observations, it is inferred that the greenhouse effect can account for the high surface temperature. Aerosols make an important contribution to the infrared opacity in these calculations. Preliminary three-dimensional calculations of the general circulation of the atmosphere are discussed which incorporate the results of the radiative calculations.

Pollack, J. B.↗

NASA NTRS1993

LSENS, a general chemical kinetics and sensitivity analysis code for gas-phase reactions: User's guide

A general chemical kinetics and sensitivity analysis code for complex, homogeneous, gas-phase reactions is described. The main features of the code, LSENS, are its flexibility, efficiency and convenience in treating many different chemical reaction models. The models include static system, steady, one-dimensional, inviscid flow, shock initiated reaction, and a perfectly stirred reactor. In addition, equilibrium computations can be performed for several assigned states. An implicit numerical integration method, which works efficiently for the extremes of very fast and very slow reaction, is used for solving the 'stiff' differential equation systems that arise in chemical kinetics. For static reactions, sensitivity coefficients of all dependent variables and their temporal derivatives with respect to the initial values of dependent variables and/or the rate coefficient parameters can be computed. This paper presents descriptions of the code and its usage, and includes several illustrative example problems.

Radhakrishnan, Krishnan↗

NASA NTRS1994

NLTE model stellar atmospheres with line blanketing near the series limits

In this paper we study the influence of line-merging regions at the intermediate long-wavelength side of a continuum threshold on the computed model atmosphere structure and predicted spectrum. In order to model these regions sufficiently accurately, we have developed two concepts. First, we have extended the occupation probability formalism of Hummer and Mihalas to non-local thermodynamic equilibrium (non-LTE) plasmas. Second, in order to treat the very complicated opacity in the line merging region, we have generalized the concept of opacity distribution functions to treat non-LTE situations. All Rydberg states are consistently included within this framework, so that no arbitrary cutoff of high (LTE) levels is made. We have calculated several pure hydrogen models atmospheres for two effective temperatures, T(sub eff) = 20000 and 35000 K, and discussed the differences between models calculated with various treatments of the line merging. In particular, we have shown that the error in the predicted profiles of Balmer lines resulting from the neglect of line merging is typically of the order of 3-4%, while the errors in the far-UV portion of the Balmer continuum reaches 15-35%. The errors generally decrease with increasing effective temperature. At the same time, the internal accuracy of the models is shown to be about or below 0.5% for all predicted spectral features. We conclude that for interpreting current high-accuracy spectrophotometric observations models including the line merging are necessary, and that the formalism developed in this paper is capable of providing a sufficiently accurate and robust modeling technique.

Hubeny, I.↗

NASA NTRS1984

Postbuckling behavior of axially compressed graphite-epoxy cylindrical panels with circular holes

The results of an experimental and analytical study of the effects of circular holes on the postbuckling behavior of graphite-epoxy cylindrical panels loaded in axial compression are presented. The STAGSC-1 general shell analysis computer code is used to determine the buckling and postbuckling response of the panels. The loaded, curved ends of the specimens were clamped by fixtures and the unloaded, straight edges were simply supported by knife-edge restraints. The panels are loaded by uniform end shortening to several times the end shortening at buckling. The unstable equilibrium path of the postbuckling response is obtained analytically by using a method based on controlling an equilibrium-path-arc-length parameter instead of the traditional load parameter. The effects of hole diameter, panel radius, and panel thickness on postbuckling response are considered in the study. Experimental results are compared with the analytical results and the failure characteristics of the graphite-epoxy panels are described.

Knight, N. F., Jr.↗

NASA NTRS2000

Mixing and NOx Emission Calculations of Confined Reacting Jet Flows in Cylindrical and Annular Ducts

Rapid mixing of cold lateral jets with hot cross-stream flows in confined configurations is of practical interest in gas turbine combustors as it strongly affects combustor exit temperature quality, and gaseous emissions in for example rich-lean combustion. It is therefore important to further improve our fundamental understanding of the important processes of dilution jet mixing especially when the injected jet mass flow rate exceeds that of the cross-stream. The results reported in this report describe some of the main flow characteristics which develop in the mixing process in a cylindrical duct. A three-dimensional computational fluid dynamics (CFD) code has been used to predict the mixing flow field characteristics and NOx emission in a quench section of a rich-burn/quick-mix/lean-burn (RQL) combustor. Sixty configurations have been analyzed in both circular and annular geometries in a fully reacting environment simulating the operating condition of an actual RQL gas turbine combustion liner. The evaluation matrix was constructed by varying the number of orifices per row and orifice shape. Other parameters such as J (momentum-flux ratio), MR (mass flowrate ratio), DR (density ratio), and mixer sector orifice ACd (effective orifice area) were maintained constant throughout the entire study. The results indicate that the mixing flow field can be correlated with the NOx production if they are referenced with the stoichiometric equivalence ratio value and not the equilibrium value. The mixing flowfields in both circular and annular mixers are different. The penetration of equal jets in both annular and circular geometries is vastly different which significantly affects the performance of the mixing section. In the computational results with the circular mixer, most of the NOx formation occurred behind the orifice starting at the orifice wake region. General trends have been observed in the NOx production as the number of orifices is changed and this appears to be common for all hole configurations and mixer types (circular or annular). The performance of any orifice shape (in producing minimum NOx) appears to be acceptable if the number of orifices can be freely varied in order to attain the optimum jet penetration.

Oechsle, Victor L.↗

NASA NTRS1960

An Analytical Investigation of Three General Methods of Calculating Chemical-Equilibrium Compositions

The Brinkley, Huff, and White methods for chemical-equilibrium calculations were modified and extended in order to permit an analytical comparison. The extended forms of these methods permit condensed species as reaction products, include temperature as a variable in the iteration, and permit arbitrary estimates for the variables. It is analytically shown that the three extended methods can be placed in a form that is independent of components. In this form the Brinkley iteration is identical computationally to the White method, while the modified Huff method differs only'slightly from these two. The convergence rates of the modified Brinkley and White methods are identical; and, further, all three methods are guaranteed to converge and will ultimately converge quadratically. It is concluded that no one of the three methods offers any significant computational advantages over the other two.

Zeleznik, Frank J.↗

NASA NTRS1975

A study of reacting free and ducted hydrogen/air jets

The mixing and reaction of a supersonic jet of hydrogen in coaxial free and ducted high temperature test gases were investigated. The importance of chemical kinetics on computed results, and the utilization of free-jet theoretical approaches to compute enclosed flow fields were studied. Measured pitot pressure profiles were correlated by use of a parabolic mixing analysis employing an eddy viscosity model. All computations, including free, ducted, reacting, and nonreacting cases, use the same value of the empirical constant in the viscosity model. Equilibrium and finite rate chemistry models were utilized. The finite rate assumption allowed prediction of observed ignition delay, but the equilibrium model gave the best correlations downstream from the ignition location. Ducted calculations were made with finite rate chemistry; correlations were, in general, as good as the free-jet results until problems with the boundary conditions were encountered.

Beach, H. L., Jr.↗

NASA NTRS1982

Hybrid state vector methods for structural dynamic and aeroelastic boundary value problems

A computational technique is developed that is suitable for performing preliminary design aeroelastic and structural dynamic analyses of large aspect ratio lifting surfaces. The method proves to be quite general and can be adapted to solving various two point boundary value problems. The solution method, which is applicable to both fixed and rotating wing configurations, is based upon a formulation of the structural equilibrium equations in terms of a hybrid state vector containing generalized force and displacement variables. A mixed variational formulation is presented that conveniently yields a useful form for these state vector differential equations. Solutions to these equations are obtained by employing an integrating matrix method. The application of an integrating matrix provides a discretization of the differential equations that only requires solutions of standard linear matrix systems. It is demonstrated that matrix partitioning can be used to reduce the order of the required solutions. Results are presented for several example problems in structural dynamics and aeroelasticity to verify the technique and to demonstrate its use. These problems examine various types of loading and boundary conditions and include aeroelastic analyses of lifting surfaces constructed from anisotropic composite materials.

Lehman, L. L.↗

NASA NTRS2009

Data Parallel Line Relaxation (DPLR) Code User Manual: Acadia - Version 4.01.1

Data-Parallel Line Relaxation (DPLR) code is a computational fluid dynamic (CFD) solver that was developed at NASA Ames Research Center to help mission support teams generate high-value predictive solutions for hypersonic flow field problems. The DPLR Code Package is an MPI-based, parallel, full three-dimensional Navier-Stokes CFD solver with generalized models for finite-rate reaction kinetics, thermal and chemical non-equilibrium, accurate high-temperature transport coefficients, and ionized flow physics incorporated into the code. DPLR also includes a large selection of generalized realistic surface boundary conditions and links to enable loose coupling with external thermal protection system (TPS) material response and shock layer radiation codes.

Wright, Michael J.↗

NASA NTRS1997

Gamma-Ray Bursts and Relativistic MHD around Coalescing Neutron Stars

In this proposal we performed the first relativistic hydrodynamic simulations of the binary neutron stars near coalescence. Calculations were performed in three spatial dimensions and were supplemented with magnetohydrodynamic simulations of the magnetic reconnection as the neutron star fluid responds to the relativistic forces of the binary. The ultimate goal of the proposed work was to develop a theoretical model for the temporal and spectral evolution of cosmological gamma-ray bursts produced during the magnetic reconnection. This is the first time such hydrodynamic calculations were performed and some surprising results were obtained. The neutrons stars were observed to compress and heat in the simulations. They could even collapse to black holes many orbits before coalescence. The question then arose as to whether this strange phenomenon could power a gamma-ray burst. The bulk of computing effort in this proposal was therefore devoted to an analysis of the implications of this possible new gamma-ray burst paradigm. The initial focus of the work performed under this proposal was therefore to utilize (3+1) dimensional and spherical numerical general relativistic hydrodynamics to study the origin, evolution, and parametric sensitivity of the collapse instability. We derived modified conditions of hydrostatic equilibrium for stars in the curved space of quasi-static orbits. We followed the magnetic field evolution for approximately one orbit (10 msec). The magnetic field energy was observed to exponentiate with an e-folding time of about 1 msec. Thus, the field should build up very quickly to a magnitude such that reconnection and back reaction of the fluid inhibits further growth.

Mathews, Grant J.↗

NASA NTRS1990

Theoretical investigation of non-equilibrium chemistry and optical radiation in hypersonic flow fields

Future space vehicles returning from distant missions or high earth orbits may enter the upper regions of the atmosphere and use aerodynamic drag to reduce their velocity before they skip out of the atmosphere and enter low earth orbit. The Aeroassist Flight Experiment (AFE) is designed to explore the special problems encountered in such entries. A computer code was developed to calculate the radiative transport along line-or-sight in the general 3-D flow field about an arbitrary entry vehicle, if the temperatures and species concentrations along the line-of-sight are known. The radiative heating calculation at the stagnation point of the AFE vehicle along the entry trajectory was performed, including a detailed line-by-line accounting of the radiative transport in the vacuum ultraviolet (below 200 nm) by the atomic N and O lines. A method was developed for making measurements of the haze particles in the Titan atmosphere above 200 km altitude. Several other tasks of a continuing nature, to improve the technical ability to calculate the nonequilibrium gas dynamic flow field and radiative heating of entry vehicles, were completed or advanced.

Whiting, Ellis E.↗

NASA NTRS1991

Addition of equilibrium air to an upwind Navier-Stokes code and other first steps toward a more generalized flow solver

An upwind three-dimensional volume Navier-Stokes code is modified to facilitate modeling of complex geometries and flow fields represented by proposed National Aerospace Plane concepts. Code enhancements include an equilibrium air model, a generalized equilibrium gas model and several schemes to simplify treatment of complex geometric configurations. The code is also restructured for inclusion of an arbitrary number of independent and dependent variables. This latter capability is intended for eventual use to incorporate nonequilibrium/chemistry gas models, more sophisticated turbulence and transition models, or other physical phenomena which will require inclusion of additional variables and/or governing equations. Comparisons of computed results with experimental data and results obtained using other methods are presented for code validation purposes. Good correlation is obtained for all of the test cases considered, indicating the success of the current effort.

Rosen, Bruce S.↗

NASA NTRS2006

Viking Afterbody Heating Computations and Comparisons to Flight Data

Computational fluid dynamics predictions of Viking Lander 1 entry vehicle afterbody heating are compared to flight data. The analysis includes a derivation of heat flux from temperature data at two base cover locations, as well as a discussion of available reconstructed entry trajectories. Based on the raw temperature-time history data, convective heat flux is derived to be 0.63-1.10 W/sq cm for the aluminum base cover at the time of thermocouple failure. Peak heat flux at the fiberglass base cover thermocouple is estimated to be 0.54-0.76 W/sq cm, occurring 16 seconds after peak stagnation point heat flux. Navier-Stokes computational solutions are obtained with two separate codes using an 8-species Mars gas model in chemical and thermal non-equilibrium. Flowfield solutions using local time-stepping did not result in converged heating at either thermocouple location. A global time-stepping approach improved the computational stability, but steady state heat flux was not reached for either base cover location. Both thermocouple locations lie within a separated flow region of the base cover that is likely unsteady. Heat flux computations averaged over the solution history are generally below the flight data and do not vary smoothly over time for both base cover locations. Possible reasons for the mismatch between flight data and flowfield solutions include underestimated conduction effects and limitations of the computational methods.

Edquist, Karl T.↗