Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Enthalpy”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,333 records · Page 74

Proteus two-dimensional Navier-Stokes computer code, version 2.0. Volume 3: Programmer's reference

A computer code called Proteus 2D was developed to solve the two-dimensional planar or axisymmetric, Reynolds-averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The objective in this effort was to develop a code for aerospace propulsion applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The governing equations are solved in generalized nonorthogonal body-fitted coordinates, by marching in time using a fully-coupled ADI solution procedure. The boundary conditions are treated implicitly. All terms, including the diffusion terms, are linearized using second-order Taylor series expansions. Turbulence is modeled using either an algebraic or two-equation eddy viscosity model. The thin-layer or Euler equations may also be solved. The energy equation may be eliminated by the assumption of constant total enthalpy. Explicit and implicit artificial viscosity may be used. Several time step options are available for convergence acceleration. The documentation is divided into three volumes. The Programmer's Reference contains detailed information useful when modifying the program. The program structure, the Fortran variables stored in common blocks, and the details of each subprogram are described.

Towne, Charles E.↗

Proteus three-dimensional Navier-Stokes computer code, version 1.0. Volume 1: Analysis description

A computer code called Proteus 3D has been developed to solve the three dimensional, Reynolds averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The objective in this effort has been to develop a code for aerospace propulsion applications that is easy to use and easy to modify. Code readability, modularity, and documentation have been emphasized. The governing equations are solved in generalized non-orthogonal body-fitted coordinates by marching in time using a fully-coupled ADI solution procedure. The boundary conditions are treated implicitly. All terms, including the diffusion terms, are linearized using second-order Taylor series expansions. Turbulence is modeled using either an algebraic or two-equation eddy viscosity model. The thin-layer or Euler equations may also be solved. The energy equation may be eliminated by the assumption of constant total enthalpy. Explicit and implicit artificial viscosity may be used. Several time step options are available for convergence acceleration. The documentation is divided into three volumes. This is the Analysis Description, and presents the equations and solution procedure. It describes in detail the governing equations, the turbulence model, the linearization of the equations and boundary conditions, the time and space differencing formulas, the ADI solution procedure, and the artificial viscosity models.

Towne, Charles E.↗

Reliability enhancement of Navier-Stokes codes through convergence enhancement

Reduction of total computing time required by an iterative algorithm for solving Navier-Stokes equations is an important aspect of making the existing and future analysis codes more cost effective. Several attempts have been made to accelerate the convergence of an explicit Runge-Kutta time-stepping algorithm. These acceleration methods are based on local time stepping, implicit residual smoothing, enthalpy damping, and multigrid techniques. Also, an extrapolation procedure based on the power method and the Minimal Residual Method (MRM) were applied to the Jameson's multigrid algorithm. The MRM uses same values of optimal weights for the corrections to every equation in a system and has not been shown to accelerate the scheme without multigriding. Our Distributed Minimal Residual (DMR) method based on our General Nonlinear Minimal Residual (GNLMR) method allows each component of the solution vector in a system of equations to have its own convergence speed. The DMR method was found capable of reducing the computation time by 10-75 percent depending on the test case and grid used. Recently, we have developed and tested a new method termed Sensitivity Based DMR or SBMR method that is easier to implement in different codes and is even more robust and computationally efficient than our DMR method.

Choi, K.-Y.↗

Calibration of HYPULSE for hypervelocity air flows corresponding to flight Mach numbers 13.5, 15, and 17

A series of air calibration tests were performed in GASL's HYPULSE facility in order to more accurately determine test section flow conditions for flows simulating total enthalpies in the Mach 13 to 17 range. Present calibration data supplements previous data and includes direct measurement of test section pitot and static pressure, acceleration tube wall pressure and heat transfer, and primary and secondary incident shock velocities. Useful test core diameters along with the corresponding free-stream conditions and usable testing times were determined. For the M13.5 condition, in-stream static pressure surveys showed the temporal and spacial uniformity of this quantity across the useful test core. In addition, finite fringe interferograms taken of the free-stream flow at the test section did not indicate the presence of any 'strong' wave system for any of the conditions investigated.

Calleja, John↗

Three-dimensional solidification and melting using magnetic field control

A new two-fluid mathematical model for fully three dimensional steady solidification under the influence of an arbitrary acceleration vector and with or without an arbitrary externally applied steady magnetic field have been formulated and integrated numerically. The model includes Joule heating and allows for separate temperature dependent physical properties within the melt and the solid. Latent heat of phase change during melting/solidification was incorporated using an enthalpy method. Mushy region was automatically captured by varying viscosity orders of magnitude between liquidus and solidus temperature. Computational results were obtained for silicon melt solidification in a parallelepiped container cooled from above and from a side. The results confirm that the magnetic field has a profound influence on the solidifying melt flow field thus changing convective heat transfer through the boundaries and the amount and shape of the solid accrued. This suggests that development of a quick-response algorithm for active control of three dimensional solidification is feasible since it would require low strength magnetic fields.

Dulikravich, George S.↗

A theorem regarding roots of the zero-order Bessel function of the first kind

This paper investigates a problem on the steady-state, conduction-convection heat transfer process in cylindrical porous heat exchangers. The governing partial differential equations for the system are obtained using the energy conservation law. Solution of these equations and the concept of enthalpy lead to a new approach to prove a theorem that the sum of inverse squares of all the positive roots of the zero order Bessel function of the first kind equals to one-forth. As a corollary, it is shown that the sum of one over pth power (p greater than or equal to 2) of the roots converges to some constant.

Lin, X.-A.↗

Containerless measurements on liquids at high temperatures

The application of containerless techniques for measurements of the thermophysical properties of high temperature liquids is reviewed. Recent results obtained in the materials research laboratories at Intersonics are also presented. Work to measure high temperature liquid properties is motivated by both the need for reliable property data for modeling of industrial processes involving molten materials and generation of data form basic modeling of materials behavior. The motivation for this work and examples of variations in thermophysical property values from the literature are presented. The variations may be attributed to changes in the specimen properties caused by chemical changes in the specimen and/or to measurement errors. The two methods used to achieve containerless conditions were aeroacoustic levitation and electromagnetic levitation. Their qualities are presented. The accompanying slides show the layout of levitation equipment and present examples of levitated metallic and ceramic specimens. Containerless techniques provide a high degree of control over specimen chemistry, nucleation and allow precise control of liquid composition to be achieved. Effects of minor additions can thus be measured in a systematic way. Operation in reduced gravity enables enhanced control of liquid motion which can allow measurement of liquid transport properties. Examples of nucleation control, the thermodynamics of oxide contamination removal, and control of the chromium content of liquid aluminum oxide by high temperature containerless processes are presented. The feasibility of measuring temperature, emissivity, liquidus temperature, enthalpy, surface tension, density, viscosity, and thermal diffusivity are discussed in the final section of the paper.

Weber, Richard↗

Effects of physical aging on long-term behavior of composites

The HSCT plane, envisioned to have a lifetime of over 60,000 flight hours and to travel at speeds in excess of Mach 2, is the source of intensive study at NASA. In particular, polymer matrix composites are being strongly considered for use in primary and secondary structures due to their high strength to weight ratio and the options of property tailoring. However, an added difficulty in the use of polymer based materials is that their properties change significantly over time, especially at the elevated temperatures that will be experienced during flight, and prediction of properties based on irregular thermal and mechanical loading is extremely difficult. This study focused on one aspect of long-term polymer composite behavior: physical aging. When a polymer is cooled to below its glass transition temperature, the material is not in thermodynamic equilibrium and the free volume and enthalpy evolve over time to approach their equilibrium values. During this time, the mechanical properties change significantly and this change is termed physical aging. This work begins with a review of the concepts of physical aging on a pure polymer system. The effective time theory, which can be used to predict long term behavior based on short term data, is mathematically formalized. The effects of aging to equilibrium are proven and discussed. The theory developed for polymers is then applied first to a unidirectional composite, then to a general laminate. Comparison to experimental data is excellent. It is shown that the effects of aging on the long-term properties of composites can be counter-intuitive, stressing the importance of the development and use of a predictive theory to analyze structures.

Brinson, L. Catherine↗

A numerical study of mixing and combustion in hypervelocity flows through a scramjet combustor model

Interest in high speed, air-breathing propulsion systems such as scramjets has revived in recent years fueled to a large extent by the National Aerospace Plane (NASP) program. These vehicles are expected to fly trans-atmospheric and as a consequence, the Mach number level within the engine/combustor would be rather high (M greater than 5). Ground based testing of such scramjet engines requires a facility that can not only achieve the right Mach number, but also have the proper pressures and temperatures to simulate the combustion processes. At present, only pulse type facilities can provide such high enthalpy flows. The newest of these is the free-piston shock tunnel, T5 located at GALCIT. Recently, a generic combustor model was tested in T5, and the experimental data from that study is analyzed in the present report. The available experimental data from T5 are essentially the static pressures on the injection wall and the one opposite to it. Thus, a principal aim of the present study was to validate the available experimental data by using a proven CFD tool and then investigate the performance characteristics of the combustor model, such as, the mixing efficiency and combustion efficiency. For this purpose, in this study, the code GASP has been used.

Krishamurthy, Ramesh↗

Calculation of the thermoneutral potential of NiCd and NiH2 cells

The thermoneutral potential of a nickel cadmium or nickel hydrogen cell is the potential at which the cell charge or discharge process puts out zero heat, and thus is the potential corresponding to the enthalpy change of the charge/discharge reaction, delta H. A relatively straightforward method for obtaining the thermoneutral potential E(sub tn), is based on the measured potential and temperature derivative of the cell reactions, which are related to the free energy change delta G, and entropy change delta S, respectively. Particularly in the nickel hydrogen cell, the pressure of hydrogen can often vary over an order of magnitude or more during the course of a charge or discharge. In a nickel cadmium cell, although significant changes in oxygen pressure can occur during charge or discharge, since oxygen does not enter into the charge/discharge reaction, these pressure changes are related to the heat generated from oxygen evolution and recombination. However, the entropy changes due to changes in hydrogen pressure relative to the 1 atm standard state must be included to apply this method to the nickel hydrogen cell.

Zimmerman, Albert H.↗

Differential Scanning Calorimetry (DSC) for planetary surface exploration

Differential Scanning Calorimetry (DSC) is the quantitative measurement of the enthalpic response of a material to a systematic change in temperature. In practice, the heat flow into or outward from a sample is measured as the sample is heated or cooled at a carefully controlled rate. DSC superficially resembles, but is not the same as differential thermal analysis (DTA), which is the measurement of temperature differences between a sample and reference material as the pair is heated or cooled. The fundamental properties measured by DSC are enthalpies and temperatures of phase transitions and constant-pressure heat capacities. Depending on instrument design and the nature of the sample, high-quality DSC analyses can be obtained on only a few milligrams of solid materials. DSC requires direct contact with the sample and generally degrades, if not destroys, the sample as a consequence of heating. In laboratory applications, it is common to subject the gaseous effluent from the DSC to analysis by a separate evolved-gas analyzer (EGA).

Gooding, James L.↗

Development of lightweight ceramic ablators and arc-jet test results

Lightweight ceramic ablators (LCA's) were recently developed at Ames to investigate the use of low density fibrous substrates and organic resins as high temperature, high strength ablative heat shields. Unlike the traditional ablators, LCA's use porous ceramic/carbon fiber matrices as substrates for structural support, and polymeric resins as fillers. Several substrates and resins were selected for the initial studies, and the best performing candidates were further characterized. Three arcjet tests were conducted to determine the LCA's thermal performance and ablation characteristics in a high enthalpy, hypersonic flow environment. Mass loss and recession measurements were obtained for each sample at post test, and the recession rates were determined from high speed motion films. Surface temperatures were also obtained from optical pyrometers.

Tran, Huy K.↗

Calorimetric thermometry of meteoritic troilite: Preliminary thermometer relationships

Thermodynamic properties of the alpha/beta phase transformation in terrestrial troilite (FeS), as measured by differential scanning calorimetry (DSC), vary systematically with prior thermal history of the troilite, as imposed under laboratory conditions. Both the transition temperature and enthalpy change for the alpha/beta transformation decrease with increasing maximum temperature of prior heat treatment. DSC measurements on troilite from various meteorites indicate clear differences in the alpha/beta thermodynamic properties that are consistent with differences in the natural thermal histories of the meteorites.

Allton, Judith H.↗

A flux splitting scheme with high-resolution and robustness for discontinuities

A flux splitting scheme is proposed for the general nonequilibrium flow equations with an aim at removing numerical dissipation of Van-Leer-type flux-vector splittings on a contact discontinuity. The scheme obtained is also recognized as an improved Advection Upwind Splitting Method (AUSM) where a slight numerical overshoot immediately behind the shock is eliminated. The proposed scheme has favorable properties: high-resolution for contact discontinuities; conservation of enthalpy for steady flows; numerical efficiency; applicability to chemically reacting flows. In fact, for a single contact discontinuity, even if it is moving, this scheme gives the numerical flux of the exact solution of the Riemann problem. Various numerical experiments including that of a thermo-chemical nonequilibrium flow were performed, which indicate no oscillation and robustness of the scheme for shock/expansion waves. A cure for carbuncle phenomenon is discussed as well.

Wada, Yasuhiro↗

Analysis and design of numerical schemes for gas dynamics. 2: Artificial diffusion and discrete shock structure

The effect of artificial diffusion on discrete shock structures is examined for a family of schemes which includes scalar diffusion, convective upwind and split pressure (CUSP) schemes, and upwind schemes with characteristics splitting. The analysis leads to conditions on the diffusive flux such that stationary discrete shocks can contain a single interior point. The simplest formulation which meets these conditions is a CUSP scheme in which the coefficients of the pressure differences is fully determined by the coefficient of convective diffusion. It is also shown how both the characteristic and CUSP schemes can be modified to preserve constant stagnation enthalpy in steady flow, leading to four variants, the E and H-characteristic schemes, and the E and H-CUSP schemes. Numerical results are presented which confirm the properties of these schemes.

Jameson, Antony↗

Shock-tunnel combustor testing for hypersonic vehicles

Proposed configurations for the next generation of transatmospheric vehicles will rely on air breathing propulsion systems during all or part of their mission. At flight Mach numbers greater than about 7 these engines will operate in the supersonic combustion ramjet mode (scramjet). Ground testing of these engine concepts above Mach 8 requires high pressure, high enthalpy facilities such as shock tunnels and expansion tubes. These impulse, or short duration facilities have test times on the order of a millisecond, requiring high speed instrumentation and data systems. One such facility ideally suited for scramjet testing is the NASA-Ames 16-Inch shock tunnel, which over the last two years has completed a series of tests for the NASP (National Aero-Space Plane) program at simulated flight Mach numbers ranging from 12-16. The focus of the experimental programs consisted of a series of classified tests involving a near-full scale hydrogen fueled scramjet combustor model in the semi-free jet method of engine testing whereby the compressed forebody flow ahead of the cowl inlet is reproduced (see appendix A). The AIMHYE-1 (Ames Integrated Modular Hypersonic Engine) test entry for the NASP program was completed in April 1993, while AIMHYE-2 was completed in May 1994. The test entries were regarded as successful, resulting in some of the first data of its kind on the performance of a near full scale scramjet engine at Mach 12-16. The data was distributed to NASP team members for use in design system verification and development. Due to the classified nature of the hardware and data, the data reports resulting from this work are classified and have been published as part of the NASP literature. However, an unclassified AIAA paper resulted from the work and has been included as appendix A. It contains an overview of the test program and a description of some of the important issues.

Loomis, Mark P.↗

LSENS: A General Chemical Kinetics and Sensitivity Analysis Code for homogeneous gas-phase reactions. Part 1: Theory and numerical solution procedures

LSENS, the Lewis General Chemical Kinetics and Sensitivity Analysis Code, has been developed for solving complex, homogeneous, gas-phase chemical kinetics problems and contains sensitivity analysis for a variety of problems, including nonisothermal situations. This report is part 1 of a series of three reference publications that describe LENS, provide a detailed guide to its usage, and present many example problems. Part 1 derives the governing equations and describes the numerical solution procedures for the types of problems that can be solved. The accuracy and efficiency of LSENS are examined by means of various test problems, and comparisons with other methods and codes are presented. LSENS is a flexible, convenient, accurate, and efficient solver for chemical reaction problems such as static system; steady, one-dimensional, inviscid flow; reaction behind incident shock wave, including boundary layer correction; and perfectly stirred (highly backmixed) reactor. In addition, the chemical equilibrium state can be computed for the following assigned states: temperature and pressure, enthalpy and pressure, temperature and volume, and internal energy and volume. For static problems the code computes the sensitivity coefficients of the dependent variables and their temporal derivatives with respect to the initial values of the dependent variables and/or the three rate coefficient parameters of the chemical reactions.

Radhakrishnan, Krishnan↗

LSENS: A General Chemical Kinetics and Sensitivity Analysis Code for homogeneous gas-phase reactions. Part 3: Illustrative test problems

LSENS, the Lewis General Chemical Kinetics and Sensitivity Analysis Code, has been developed for solving complex, homogeneous, gas-phase chemical kinetics problems and contains sensitivity analysis for a variety of problems, including nonisothermal situations. This report is part 3 of a series of three reference publications that describe LSENS, provide a detailed guide to its usage, and present many example problems. Part 3 explains the kinetics and kinetics-plus-sensitivity analysis problems supplied with LSENS and presents sample results. These problems illustrate the various capabilities of, and reaction models that can be solved by, the code and may provide a convenient starting point for the user to construct the problem data file required to execute LSENS. LSENS is a flexible, convenient, accurate, and efficient solver for chemical reaction problems such as static system; steady, one-dimensional, inviscid flow; reaction behind incident shock wave, including boundary layer correction; and perfectly stirred (highly backmixed) reactor. In addition, the chemical equilibrium state can be computed for the following assigned states: temperature and pressure, enthalpy and pressure, temperature and volume, and internal energy and volume. For static problems the code computes the sensitivity coefficients of the dependent variables and their temporal derivatives with respect to the initial values of the dependent variables and/or the three rate coefficient parameters of the chemical reactions.

Bittker, David A.↗