Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermodynamics of computing”

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 307 records · Page 17

Laser Induced Aluminum Surface Breakdown Model

Laser powered propulsion systems involve complex fluid dynamics, thermodynamics and radiative transfer processes. Based on an unstructured grid, pressure-based computational aerothermodynamics; platform, several sub-models describing such underlying physics as laser ray tracing and focusing, thermal non-equilibrium, plasma radiation and air spark ignition have been developed. This proposed work shall extend the numerical platform and existing sub-models to include the aluminum wall surface Inverse Bremsstrahlung (IB) effect from which surface ablation and free-electron generation can be initiated without relying on the air spark ignition sub-model. The following tasks will be performed to accomplish the research objectives.

Chen, Yen-Sen↗

Radiative Losses in the Chromosphere During A C-Class Flare

Context. Solar flares release an enormous amount of energy (~10(exp 32) erg) into the corona. A substantial fraction of this energy is transported to the lower atmosphere, which results in chromospheric heating. The mechanisms that transport energy to the lower solar atmosphere during a flare are still not fully understood. Aims. We aim to estimate the temporal evolution of the radiative losses in the chromosphere at the footpoints of a C-class flare, in order to set observational constraints on the electron beam parameters of a RADYN flare simulation. Methods. We estimated the radiative losses from hydrogen, and singly ionized Ca and Mg using semiempirical model atmospheres, which were inferred from a multiline inversion of observed Stokes profiles obtained with the CRISP and CHROMIS instruments on the Swedish 1-m Solar Telescope. The radiative losses were computed taking into account the effect of partial redistribution and non-local thermodynamic equilibrium. To estimate the integrated radiative losses in the chromosphere, the net cooling rates were integrated between the temperature minimum and the height where the temperature reaches 10 kK. We also compared our time series of radiative losses with those from the RADYN flare simulations. Results. We obtained a high spatial-resolution map of integrated radiative losses around the flare peak time. The stratification of the net cooling rate suggests that the Ca IR triplet lines are responsible for most of the radiative losses in the flaring atmosphere. During the flare peak time, the contribution from Ca II H and K and Mgii h and k lines are strong and comparable to the Ca IR triplet (~32kW m(exp -2)). Since our flare is a relatively weak event, the chromosphere is not heated above 11 kK, which in turn yields a subdued Lyα contribution (~7kW m(exp -2)) in the selected limits of the chromosphere. The temporal evolution of total integrated radiative losses exhibits sharply rising losses (0.4kW m(exp -2) (s(exp -1)) and a relatively slow decay (0.23kW m(exp -2) s(exp -1)). The maximum value of total radiative losses is reached around the flare peak time and can go up to 175kWm􀀀2 for a single pixel located at footpoint. After a small parameter study, we find the best model-data consistency in terms of the amplitude of radiative losses and the overall atmospheric structure with a RADYN flare simulation in the injected energy flux of 5 × 10(exp 10) erg s(exp -1) cm(exp -2).

Chromosphere↗

3D Material Response of the MSL Heatshield Using NuSil-Coated PICA

The Mars Science Laboratory (MSL) was protected during its atmospheric entry by an instrumented heatshield that used NASA's Phenolic Impregnated Carbon Ablator (PICA) material [1]. PICA is a lightweight carbon fiber/polymeric resin material that offers outstanding performance for protecting probes during planetary entry. Data from the Mars Entry Descent and Landing Instrument (MEDLI) suite on MSL offers unique in-flight validation data for models of material response and atmospheric entry. MEDLI recorded, among other things, time-resolved in-depth temperature data of PICA using thermocouple sensors assembled in the MEDLI Integrated Sensor Plugs (MISP) [2]. A space-grade silicone-based coating commercially known as NuSil CV-1144-0 [3] was applied to the entire MSL heatshield, including the MEDLI plugs, to mitigate the spread of dust from PICA. Modeling the thermal response of PICA-NuSil (PICA-N) system is still an open challenge. Ground testing of PICA-N models exhibited surface temperature jumps of the order of 150 K due to oxide scale formation and sub-sequent NuSil burn-off. It is therefore critical to include a validated model for the material response of the coating in engineering codes. A test campaign has been conducted at the NASA’s Langley HyMETS [4] facility to screen the response of PICA-N and gather detailed data on its behavior [5]. A first model of PICA-N thermal response has been developed using the Hy-METS experiments [6]. The objective of this work is to analyze the material response of the latest PICA-N model compared to the engineering model used to simulate the entry of MSL. The environment and material response around the MSL aeroshell during Mars atmospheric entry is simulated using a collection of tools. The Direct Simulation Monte Carlo SPARTA code [7] is used in the rarefied regime, the Data Parallel Line Relaxation (DPLR) code [8] is used in the continuum regime and radiative heating conditions are provided by the Nonequilibrium air radiation (NEQAIR) code [9] to estimate the environmental conditions. The thermal response inside the material is computed using the Porous material Analysis Toolbox based on Open-FOAM (PATO) [10,11,12]. Thermodynamic and chemistry properties are estimated using the Mutation++ library [13]. The approach implemented in PATO as a first cut PICA-N thermal response model is outlined in Figure 1. While the recession is less than the coating thickness, the Surface mass and energy balance Boundary Condition (SBC) uses the NuSil B’ tables. Once the recession removes the coating, the usual PICA B’ tables are used for the SBC. The B’ tables are computed using an equilibrium solver implemented in Mutation++, given the temperature, pressure, blowing rate, composition of the pyrolysis and environment gases, and the condensed species at the surface. Preliminary results of the 3D material response of the MSL heat-shield at the peak heating (80 sec after Entry Interface) are shown in Figure 2. Current NASA’s mission to Mars, Mars 2020, used the spare heatshield of MSL for thermal protection during entry, descent, and landing. In preparation for Mars 2020 post-flight analysis, the PATO high-fidelity material response capability was benchmarked against flight data from MEDLI. This effort represents an important milestone toward the development of validated predictive capabilities for designing thermal protection systems for planetary probes. This bench-marking is awaiting the final release of the MEDLI-2 data.

Aerospace↗

Research study on neutral thermodynamic atmospheric model

The Global Reference Atmospheric Model is used along with the revised perturbation statistics to evaluate and computer graph various atmospheric statistics along a space shuttle reference mission and abort trajectory. The trajectory plots are height vs. ground range, with height from ground level to 155 km and ground range along the reentry trajectory. Cross sectional plots, height vs. latitude or longitude, are also generated for 80 deg longitude, with heights from 30 km to 90 km and latitude from -90 deg to +90 deg, and for 45 deg latitude, with heights from 30 km to 90 km and longitudes from 180 deg E to 180 deg W. The variables plotted are monthly average pressure, density, temperature, wind components, and wind speed and standard deviations and 99th inter-percentile range for each of these variables.

Hargraves, W. R.↗

Computing Properties Of Chemical Mixtures At Equilibrium

Scientists and engineers need data on chemical equilibrium compositions to calculate theoretical thermodynamic properties of chemical systems. Information essential in design and analysis of such equipment as compressors, turbines, nozzles, engines, shock tubes, heat exchangers, and chemical-processing equipment. CET93 is general program that calculates chemical equilibrium compositions and properties of mixtures for any chemical system for which thermodynamic data are available. Includes thermodynamic data for more than 1,300 gaseous and condensed species and thermal-transport data for 151 gases. Written in FORTRAN 77.

Mcbride, B. J.↗

Unstructured Finite Volume Computational Thermo-Fluid Dynamic Method for Multi-Disciplinary Analysis and Design Optimization

This paper describes a finite volume computational thermo-fluid dynamics method to solve for Navier-Stokes equations in conjunction with energy equation and thermodynamic equation of state in an unstructured coordinate system. The system of equations have been solved by a simultaneous Newton-Raphson method and compared with several benchmark solutions. Excellent agreements have been obtained in each case and the method has been found to be significantly faster than conventional Computational Fluid Dynamic(CFD) methods and therefore has the potential for implementation in Multi-Disciplinary analysis and design optimization in fluid and thermal systems. The paper also describes an algorithm of design optimization based on Newton-Raphson method which has been recently tested in a turbomachinery application.

Majumdar, Alok↗

Theoretical kinetic computations in complex reacting systems

Nasa Lewis' studies of complex reacting systems at high temperature are discussed. The changes which occur are the result of many different chemical reactions occurring at the same time. Both an experimental and a theoretical approach are needed to fully understand what happens in these systems. The latter approach is discussed. The differential equations which describe the chemical and thermodynamic changes are given. Their solution by numerical techniques using a detailed chemical mechanism is described. Several different comparisons of computed results with experimental measurements are also given. These include the computation of (1) species concentration profiles in batch and flow reactions, (2) rocket performance in nozzle expansions, and (3) pressure versus time profiles in hydrocarbon ignition processes. The examples illustrate the use of detailed kinetic computations to elucidate a chemical mechanism and to compute practical quantities such as rocket performance, ignition delay times, and ignition lengths in flow processes.

Bittker, David A.↗

A 4-cylinder Stirling engine computer program with dynamic energy equations

A computer program for simulating the steady state and transient performance of a four cylinder Stirling engine is presented. The thermodynamic model includes both continuity and energy equations and linear momentum terms (flow resistance). Each working space between the pistons is broken into seven control volumes. Drive dynamics and vehicle load effects are included. The model contains 70 state variables. Also included in the model are piston rod seal leakage effects. The computer program includes a model of a hydrogen supply system, from which hydrogen may be added to the system to accelerate the engine. Flow charts are provided.

Daniele, C. J.↗

Theoretical and computational dynamics of a compressible flow

An introduction to the theoretical and computational fluid dynamics of a compressible fluid is presented. The general topics addressed include: thermodynamics and physical properties of compressible fluids; 1D flow of an inviscid compressible fluid; shock waves; fundamental equations of the dynamics of a compressible inviscid non-heat-conducting and radiating fluid, method of small perturbations, linearized theory; 2D subsonic steady potential flow; hodograph and rheograph methods, exact solutions of 2D insentropic steady flow equations, 2D steady transonic and hypersonic flows, method of characteristics, linearized theory of 3D potential flow, nonlinear theory of 3D compressibe flow, anisentropic (rotational) flow of inviscid compressible fluid, electromagnetogasdynamics, multiphase flows, flows of a compressible fluid with transport phenomena.

Pai, Shih-I↗

Multispectral Emission of the Sun during the First Whole Sun Month: Magnetohydrodynamic Simulations

We demonstrate that a three-dimensional magnetohydrodynamic (MHD) simulation of the corona can model its global plasma density and temperature structure with sufficient accuracy to reproduce many of the multispectral properties of the corona observed in extreme ultraviolet (EW) and X-ray emission. The key ingredient to this new type of global MHD model is the inclusion of energy transport processes (coronal heating, anisotropic thermal conduction, and radiative losses) in the energy equation. The calculation of these processes has previously been confined to one-dimensional loop models, idealized two-dimensional computations, and three-dimensional active region models. We refer to this as the thermodynamic MHD model, and we apply it to the time period of Carrington rotation 1913 (1996 August 22 to September 18). The form of the coronal heating term strongly affects the plasma density and temperature of the solutions. We perform our calculation for three different empirical heating models: (1) a heating function exponentially decreasing in radius; (2) the model of Schrijver et al.; and (3) a model reproducing the heating properties of the quiet Sun and active regions. We produce synthetic emission images from the density and temperature calculated with these three heating functions and quantitatively compare them with observations from E W Imaging Telescope on the Solar and Heliospheric Observatory and the soft X-ray telescope on Yohkoh. Although none of the heating models provide a perfect match, heating models 2 and 3 provide a reasonable match to the observations.

Lionello, Roberto↗

Simplified curve fits for the thermodynamic properties of equilibrium air

New, improved curve fits for the thermodynamic properties of equilibrium air have been developed. The curve fits are for pressure, speed of sound, temperature, entropy, enthalpy, density, and internal energy. These curve fits can be readily incorporated into new or existing computational fluid dynamics codes if real gas effects are desired. The curve fits are constructed from Grabau-type transition functions to model the thermodynamic surfaces in a piecewise manner. The accuracies and continuity of these curve fits are substantially improved over those of previous curve fits. These improvements are due to the incorporation of a small number of additional terms in the approximating polynomials and careful choices of the transition functions. The ranges of validity of the new curve fits are temperatures up to 25 000 K and densities from 10 to the -7 to 10 to the 3d power amagats.

Srinivasan, S.↗

Nonequilibrium radiation and chemistry models for aerocapture vehicle flowfields, volume 3

The computer programs developed to calculate the shock wave precursor and the method of using them are described. This method calculated the precursor flow field in a nitrogen gas including the effects of emission and absorption of radiation on the energy and composition of gas. The radiative transfer is calculated including the effects of absorption and emission through the line as well as the continuum process in the shock layer and through the continuum processes only in the precursor. The effects of local thermodynamic nonequilibrium in the shock layer and precursor regions are also included in the radiative transfer calculations. Three computer programs utilized by this computational scheme to calculate the precursor flow field solution for a given shock layer flow field are discussed.

Carlson, Leland A.↗

Multi-dimensional modeling of pressurization and expulsion of multi-phase hydrogen propellant

A multidimensional computational model of the pressurization process in a slush hydrogen propellant storage tank was developed and its accuracy evaluated by comparison to experimental data measured for a 5 ft diameter spherical tank. The fluid mechanic, thermodynamic, and heat transfer processes within the ullage are represented by a finite-volume model. the heat and mass fluxes at the ullage boundary were computed in auxiliary analyses and specified as input to the finite-volume model. the model was shown to be in reasonable agreement with the experimental data. A parameter study was undertaken to examine the dependence of the pressurization process on initial ullage temperature distribution and pressurant mass flow rate. It is shown that for a given heat flux rate at the ullage boundary, the pressurization process is nearly independent of initial temperature distribution. The mass flow rate study revealed decreasing pressurant mass requirements with increasing pressurant mass flow rate. Further, significant differences were identified between the ullage temperature and velocity fields predicted for pressurization of slush and those predicted for pressurization of liquid hydrogen. A simplified model of the pressurization process was constructed in search of a dimensionless characterization of the pressurization process. It is shown that the relationship derived from this simplified model collapses all of the pressure history data generated during this study into a single curve.

Source record↗

Calculating Theromodynamic And Transport Properties Of Fluids

Computer program incorporates van der Waals equation and correction tables. FLUID program developed to calculate thermodynamic and transport properties of pure fluids in both liquid and gas phases. Properties calculated by use of simple gas model, empirical corrections, and efficient numerical interpolation scheme. Produces results that agree very well with measured values. Much faster than older, more complex programs developed for same purpose.

Proctor, Margaret P.↗

Solar oscillations and the equation of state

The dependence of frequencies of solar oscillations on the thermodynamic state of the sun is considered. On the basis of an analysis of the equations of stellar structure, as well as the relevant aspects of the properties of the oscillations, it is argued that in the convection zone, information about the equation of state which is relatively unaffected by other uncertainties in the physics of the solar interior can be isolated. The different treatments that have been employed to describe the thermodynamics of stellar plasmas are reviewed. Through application of several of these treatments to the computation of models of the solar envelope, it is demonstrated that the sensitivity of the observed frequencies is in fact sufficient to distinguish even quite subtle features of the physics of solar matter.

Christensen-Dalsgaard, Jorgen↗

Solid Surface Combustion Experiment Completes a Series of Eight Successful Flights

The Solid Surface Combustion Experiment (SSCE) was the first combustion experiment to fly in the space shuttle and the first such experiment in the NASA spaceflight program since Skylab. SSCE was actually a series of experiments designed to begin to characterize flame spreading over solid fuels in microgravity and the differences of this flame spreading from normal gravity behavior. These experiments should lead to a better understanding of the physical processes involved--increasing our understanding of fire behavior, both in space and on Earth. SSCE results will help researchers evaluate spacecraft fire hazards. These experiments were conceived by the principal investigator, Professor Robert A. Altenkirch, Dean of Engineering at Washington State University. In the first five flights, the fuel sample--ashless filter paper instrumented with three thermocouples--was mounted in a sealed chamber filled with a 50-percent or 35-percent mixture of oxygen in nitrogen at pressures of 1.0, 1.5, and 2.0 atm. In the next three flights, a polymethyl methacrylate (plexiglass) fuel was instrumented with three thermocouples and tested in a 70-percent or 50-percent mixture of oxygen and nitrogen at pressures of 1.0 and 2.0 atm. SSCE is a self-contained, battery-operated experiment that can be flown either in the shuttle middeck or in the Spacelab module. More information about the hardware configuration have been published. This past year, the final two of eight flights were completed on STS-64 and STS-63. The NASA Lewis Research Center designed and built the SSCE payload and performed engineering, testing, scientific, and flight operations support. The SSCE project was supported in some way by nearly every major sector of Lewis' organization. Professor Altenkirch developed a numerical simulation of the flame-spreading process from first principles (of fluid mechanics, heat transfer, and reaction kinetics). The spread rates, flame shape, and thermodynamic data from the SSCE flights are being compared directly with the results of the computational model. Results from the eight flights will be used to formulate an improved solid-phase pyrolysis model. In addition, some results of the flights have been published and presented at international combustion symposiums. Additional solid fuel combustion experiments are being investigated for future tests with the existing hardware.

Source record↗