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 181 records · Page 10

Computer program for analysis of split-Stirling-cycle cryogenic coolers

A computer program for predicting the detailed thermodynamic performance of split-Stirling-cycle refrigerators has been developed. The mathematical model includes the refrigerator cold head, free-displacer/regenerator, gas transfer line, and provision for modeling a mechanical or thermal compressor. To allow for dynamic processes (such as aerodynamic friction and heat transfer) temperature, pressure, and mass flow rate are varied by sub-dividing the refrigerator into an appropriate number of fluid and structural control volumes. Of special importance to modeling of cryogenic coolers is the inclusion of real gas properties, and allowance for variation of thermo-physical properties such as thermal conductivities, specific heats and viscosities, with temperature and/or pressure. The resulting model, therefore, comprehensively simulates the split-cycle cooler both spatially and temporally by reflecting the effects of dynamic processes and real material properties.

Brown, M. T.↗

Integrated Analysis Capability Program

Integrated Analysis Capability (IAC) software system intended to provide highly effective, interactive analysis tool for integrated design of large structures. Supports needs of engineering analysis groups concerned with interdisciplinary problems. Developed to serve as software interface between computer programs from fields of structures, thermodynamics, controls, and dynamics of systems on one hand and executive software system and data base on other hand to yield highly efficient multi-disciplinary system. Special attention given to such users' requirements as handling data and online assistance with operational features and ability to add new modules of user's choice at future date. Written in FORTRAN 77.

Vos, R. G.↗

Aerothermodynamics of a MESUR Martian entry

Computational fluid dynamics models for the thermodynamics and transport properties used in an equilibrium version of the Langley Aerothermodynamics Upwind Relaxation Algorithm (LAURA) for Martian Atmospheric entries are described. In addition, the physical models used in a nonequilibrium version of LAURA for Martian-entry flows are described. Uncertainties exist in defining constants used in the transport properties for the equilibrium model and in many of the physical models for the nonequilibrium version. Solutions from the two codes using the best available constants are examined at the Martian-entry conditions characteristic of the Mars Environmental Survey (MESUR) mission. While the flowfields are near thermal equilibrium, chemical nonequilibrium effects are present in the entry cases examined. Convective heating at the stagnation point for these flows (assuming fully catalytic wall boundary conditions) is approximately 100 W/sq cm. Radiative heating is negligible.

Mitcheltree, R. A.↗

Zero-G Thermodynamic Venting System (TVS) Performance Prediction Program

This report documents the Zero-g Thermodynamic Venting System (TVS) performance prediction computer program. The zero-g TVS is a device that destratifies and rejects environmentally induced zero-g thermal gradients in the LH2 storage transfer system. A recirculation pump and spray injection manifold recirculates liquid throughout the length of the tank thereby destratifying both the ullage gas and liquid bulk. Heat rejection is accomplished by the opening of the TVS control valve which allows a small flow rate to expand to a low pressure thereby producing a low temperature heat sink which is used to absorb heat from the recirculating liquid flow. The program was written in FORTRAN 77 language on the HP-9000 and IBM PC computers. It can be run on various platforms with a FORTRAN compiler.

Nguyen, Han↗

The early spectral evolution of Nova Cassiopeiae 1993

In this paper we describe the evolution of Nova Cas 1993 over the first two months of its outburst. We present an ultraviolet light curve that covers the period from announcement to just after dust began forming in the ejecta (1994 Feb. 15) and International Ultraviolet Explorer (IUE) spacecraft constraints forced us to halt our observations. We have used spherical, expanding, Non-local Thermodynamic Equilibrium (NLTE) stellar atmospheres to compute synthetic spectra and have compared the results to combined ultraviolet (low-resolution 1200-3400 A and high-resolution 2400-3300 A) spectra. Our fits show that the effective temperature of the ejecta increased from approx. 8000 to about approx. 16 000 K between 1993 Dec. 12 and 1993 Dec. 26. The temperature then increased more slowly to approx. 24 000 on 1994 Jan. 28. A preliminary abundance analysis shows evidence for hydrogen depletion, as we also found for Nova V1974 Cygni; however we find a larger enhancement of carbon, nitrogen, and oxygen. We also show that the principal mechanism for mass ejection in this nova is a radiation pressure driven wind and that mechanical driving is not necessary.

Hauschildt, Peter H.↗

Aerothermodynamic methods for a Mars environmental survey Mars entry

Computational fluid dynamics models for the thermodynamics and transport properties used in an equilibrium version of the Langley aerothermodynamics upwind relaxation algorithm (LAURA) for Mars atmospheric entries are described. In addition, the physical models used in a nonequilibrium version of LAURA for Mars-entry flows are described. Uncertainties exist in defining constants used in the transport properties for the equilibrium model and in many of the physical models for the nonequilibrium version. Solutions from the two codes using the best available constants are examined at the Mars-entry conditions characteristics of the Mars environmental survey mission. While the flowfields are near thermal equilibrium, chemical nonequilibrium effects are present in the entry cases examined. Convective heating at the stagnation point for these flows (assuming fully catalytic wall boundary conditions) is approximately 100 W/cm(exp 2). Radiative heating is negligible.

AEROTHERMODYNAMIC METHODS↗

Risk-Based Probabilistic Approach to Aeropropulsion System Assessment

In an era of shrinking development budgets and resources, where there is also an emphasis on reducing the product development cycle, the role of system assessment, performed in the early stages of an engine development program, becomes very critical to the successful development of new aeropropulsion systems. A reliable system assessment not only helps to identify the best propulsion system concept among several candidates, it can also identify which technologies are worth pursuing. This is particularly important for advanced aeropropulsion technology development programs, which require an enormous amount of resources. In the current practice of deterministic, or point-design, approaches, the uncertainties of design variables are either unaccounted for or accounted for by safety factors. This could often result in an assessment with unknown and unquantifiable reliability. Consequently, it would fail to provide additional insight into the risks associated with the new technologies, which are often needed by decision makers to determine the feasibility and return-on-investment of a new aircraft engine. In this work, an alternative approach based on the probabilistic method was described for a comprehensive assessment of an aeropropulsion system. The statistical approach quantifies the design uncertainties inherent in a new aeropropulsion system and their influences on engine performance. Because of this, it enhances the reliability of a system assessment. A technical assessment of a wave-rotor-enhanced gas turbine engine was performed to demonstrate the methodology. The assessment used probability distributions to account for the uncertainties that occur in component efficiencies and flows and in mechanical design variables. The approach taken in this effort was to integrate the thermodynamic cycle analysis embedded in the computer code NEPP (NASA Engine Performance Program) and the engine weight analysis embedded in the computer code WATE (Weight Analysis of Turbine Engines) with the fast probability integration technique (FPI). FPI was developed by Southwest Research Institute under contract with the NASA Glenn Research Center. The results were plotted in the form of cumulative distribution functions and sensitivity analyses and were compared with results from the traditional deterministic approach. The comparison showed that the probabilistic approach provides a more realistic and systematic way to assess an aeropropulsion system. The current work addressed the application of the probabilistic approach to assess specific fuel consumption, engine thrust, and weight. Similarly, the approach can be used to assess other aspects of aeropropulsion system performance, such as cost, acoustic noise, and emissions. Additional information is included in the original extended abstract.

Tong, Michael T.↗

Computational Tools to Design Molten Salt Electrolytes for High-Temperature Environments

Current Li-ion batteries are designed for a small operating window of 5 °C to 55 °C. Modifications to the electrolyte for operation in Mars atmosphere extended this range on the lower bound to ~ -40 °C. These operating-temperature ranges are far from those presented in high-temperature environments, such as the Venus surface, where temperatures are around 450 °C. Protecting the state-of-the-art Li-ion batteries require insulations that decreases the volumetric capacity and limiting the operational time-window. These unique challenges require a paradigm shift in materials used for designing high-temperature batteries. Molten salt electrolyte-based batteries offer a plausible route to designing high-temperature batteries. ZEBRA batteries are known to be one of the safest energy storage devices operating at 270-350 °C. Nitrate based eutectics, operating at 150 °C have also been used in Li-O2 batteries. In this study, we will examine molten-salt electrolyte transport and electrochemical properties using first-principles computations and benchmark against experiments. Further, we will present thermodynamics-based models for designing and predicting melting point of molten salt mixtures. Based on these computational tools, new molten salt mixtures designed with desired operating temperatures and electrochemical windows will be presented. Implications of these new eutectics in the context of high-temperature environment exploration will also be discussed.

Balachandran Radhakrishnan↗

Addressing the hypersonic simulation problem

The hypersonic simulation problem is solved by the meshing of bench-marked experimental data with results from validated computational fluid dynamic codes. The example presented is for the reentry of the Space Shuttle Orbiter. Mach number effects were assessed by parametrically varying free-stream Mach number and angle of attack in a series of inviscid, perfect-gas computations carried out on a modified Orbiter geometry. Real-gas effects were determined by making calculations at specific points of the reentry trajectory using equilibrium air thermodynamics and comparing with corresponding perfect-gas computations. Viscous computations were also made for both the basic Orbiter body and control surfaces. A methodology is developed that permits the extrapolation of wind tunnel data to flight situations, providing a solution to the hypersonic simulation problem.

Griffith, B. J.↗

Comparisons of CFD Simulations of Icing Wind Tunnel Clouds with Experiments Conducted at the NASA Propulsions Systems Laboratory

This paper evaluates simulation predictions against experimental test data of icing clouds that were produced during 2018 ice crystal icing physics tests conducted at the NASA Propulsion Systems Laboratory icing wind tunnel. Aero-thermal and cloud parameters are set and known upstream at the tunnel inlet and spray system, but change as the cloud and air thermodynamically interact as the flowing masses reach the tunnel test section. Utilizing the ANSYS Fluent Discrete Phase Model function, 3D computational fluid dynamics (CFD) simulations were performed, capturing the thermodynamic interactions between the test parameters, and providing predictions of the aero thermal and cloud conditions at the tunnel test section. Simulation predictions were compared with test data that were measured at the tunnel exit plane. Evaluations focused on the cloud concentration (total water content), humidity content, and air temperature. CFD simulation predictions showed areas of agreement and disagreement. Simulations showed that cloud concentration profiles at the test section are strongly related to the initial spray nozzle pattern used at the tunnel inlet. Experimental data suggest that greater dispersion of the cloud occurred as the simulated cloud predicted areas of high and low cloud concentration compared to test data profiles. Simulations, however, captured the magnitude and location of the change in humidity content and the change in air temperature due to the presence of the cloud reasonably well, when compared to test data. This result would suggest that while the ANSYS simulation did not fully predict the spreading of the cloud as measured during experiment, it did capture evaporation and the molecular movements of air and vapor relatively well.

CFD↗

Comparisons of CFD Simulations of Icing Wind Tunnel Clouds with Experiments Conducted at the NASA Propulsions Systems Laboratory

This paper evaluates simulation predictions against experimental test data of icing clouds that were produced during 2018 ice crystal icing physics tests conducted at the NASA Propulsion Systems Laboratory icing wind tunnel. Aero-thermal and cloud parameters are set and known upstream at the tunnel inlet and spray system, but change as the cloud and air thermodynamically interact as the flowing masses reach the tunnel test section. Utilizing the ANSYS Fluent Discrete Phase Model function, 3D computational fluid dynamics (CFD) simulations were performed, capturing the thermodynamic interactions between the test parameters, and providing predictions of the aero thermal and cloud conditions at the tunnel test section. Simulation predictions were compared with test data that were measured at the tunnel exit plane. Evaluations focused on the cloud concentration (total water content), humidity content, and air temperature. CFD simulation predictions showed areas of agreement and disagreement. Simulations showed that cloud concentration profiles at the test section are strongly related to the initial spray nozzle pattern used at the tunnel inlet. Experimental data suggest that greater dispersion of the cloud occurred as the simulated cloud predicted areas of high and low cloud concentration compared to test data profiles. Simulations, however, captured the magnitude and location of the change in humidity content and the change in air temperature due to the presence of the cloud reasonably well, when compared to test data. This result would suggest that while the ANSYS simulation did not fully predict the spreading of the cloud as measured during experiment, it did capture evaporation and the molecular movements of air and vapor relatively well.

CFD↗

CFD and Thermodynamic Model Predictions of No-Vent Tank Filling in Microgravity

Cryogenic Fluid Management (CFM) of propellant fluids under low gravity conditions present an important challenge in the pursuit of NASA’s future long-duration space missions. Success of these missions relies to a large extent on the ability to store and transfer the cryogenic liquid fuels efficiently and safely. To this end, the Zero-Boil-Off Tank Filling and Transfer Experiment (ZBOT-FT) undertaken by University of Bremen in collaboration with Case Western Reserve University aims at studying two phase flow tank-to-tank transfer operations in microgravity using a simulant fluid, Perfluoro-n-pentane (PNP). As part of this effort, in this work, we present a system level homogeneous thermodynamic model and a two-phase Computational Fluid Dynamics (CFD) model of a no-vent cylindrical tank filling in microgravity. Both 2D axisymmetric and 3D parametric simulations case studies are performed to show the detailed physics of the filling process and the evolution of tank pressure for filling flow rates between 1.0 and 1.65 mL/s. A zero-dimensional homogeneous thermodynamic model is also developed to predict the final steady state tank pressure analytically. Comparison of the 2D and 3D CFD predictions of final pressures to the steady state thermodynamic results displays very good agreement. However, it is shown that the 2D axisymmetric model predicts a maxima in tank pressure before reaching steady state. It is indicted that this pressure maxima is caused by the constraints inherent in the axisymmetric model and is absent in the 3D and thermodynamic pressure predictions.

Two Phase Flow↗

Generalized Fluid System Simulation Program

A general purpose program implemented on a computer analyzes steady state and transient flow in a complex fluid network, modeling phase changes, compressibility, mixture thermodynamics and external body forces such as gravity and centrifugal force. A preprocessor provides for the inter- active development of a fluid network simulation having nodes and branches. Mass, energy, and specie conservation equations are solved at the nodes, and momentum conservation equations are solved in the branches. Contained herein are subroutines for computing "real fluid" thermodynamic and thermophysical properties for 12 fluids, and a number of different source options are provided for model- ing momentum sources or sinks in the branches. The system of equations describing the fluid network is solved by a hybrid numerical method that is a combination of the Newton-Raphson and successive substitution methods. Application and verification of this invention are provided through an example problem, which demonstrates that the predictions of the present invention compare most reasonably with test data.

Alok Kumar Majumdar↗

Progress of Stirling cycle analysis and loss mechanism characterization

An assessment of Stirling engine thermodynamic modeling and design codes shows a general deficiency; this deficiency is due to poor understanding of the fluid flow and heat transfer phenomena that occur in the oscillating flow and pressure level environment within the engines. Stirling engine thermodynamic loss mechanisms are listed. Several experimental and computational research efforts now underway to characterize various loss mechanisms are reviewed. The need for additional experimental rigs and rig upgrades is discussed. Recent developments and current efforts in Stirling engine thermodynamic modeling are also reviewed.

Tew, R. C., Jr.↗

Three-dimensional radiative transfer on a massively parallel computer

We perform 3D radiative transfer calculations in non-local thermodynamic equilibrium (NLTE) in the simple two-level atom approximation on the Mas-Par MP-1, which contains 8192 processors and is a single instruction multiple data (SIMD) machine, an example of the new generation of massively parallel computers. On such a machine, all processors execute the same command at a given time, but on different data. To make radiative transfer calculations efficient, we must re-consider the numerical methods and storage of data. To solve the transfer equation, we adopt the short characteristic method and examine different acceleration methods to obtain the source function. We use the ALI method and test local and non-local operators. Furthermore, we compare the Ng and the orthomin methods of acceleration. We also investigate the use of multi-grid methods to get fast solutions for the NLTE case. In order to test these numerical methods, we apply them to two problems with and without periodic boundary conditions.

Vath, H. M.↗

Statistical equilibrium calculations for silicon in early-type model stellar atmospheres

Line profiles of 36 multiplets of silicon (Si) II, III, and IV were computed for a grid of model atmospheres covering the range from 15,000 to 35,000 K in effective temperature and 2.5 to 4.5 in log (gravity). The computations involved simultaneous solution of the steady-state statistical equilibrium equations for the populations and of the equation of radiative transfer in the lines. The variables were linearized, and successive corrections were computed until a minimal accuracy of 1/1000 in the line intensities was reached. The common assumption of local thermodynamic equilibrium (LTE) was dropped. The model atmospheres used also were computed by non-LTE methods. Some effects that were incorporated into the calculations were the depression of the continuum by free electrons, hydrogen and ionized helium line blocking, and auto-ionization and dielectronic recombination, which later were found to be insignificant. Use of radiation damping and detailed electron (quadratic Stark) damping constants had small but significant effects on the strong resonance lines of Si III and IV. For weak and intermediate-strength lines, large differences with respect to LTE computations, the results of which are also presented, were found in line shapes and strengths. For the strong lines the differences are generally small, except for the models at the hot, low-gravity extreme of our range. These computations should be useful in the interpretation of the spectra of stars in the spectral range B0-B5, luminosity classes III, IV, and V.

Kamp, L. W.↗

GVIPS Models and Software

Two reports discuss, respectively, (1) the generalized viscoplasticity with potential structure (GVIPS) class of mathematical models and (2) the Constitutive Material Parameter Estimator (COMPARE) computer program. GVIPS models are constructed within a thermodynamics- and potential-based theoretical framework, wherein one uses internal state variables and derives constitutive equations for both the reversible (elastic) and the irreversible (viscoplastic) behaviors of materials. Because of the underlying potential structure, GVIPS models not only capture a variety of material behaviors but also are very computationally efficient. COMPARE comprises (1) an analysis core and (2) a C++-language subprogram that implements a Windows-based graphical user interface (GUI) for controlling the core. The GUI relieves the user of the sometimes tedious task of preparing data for the analysis core, freeing the user to concentrate on the task of fitting experimental data and ultimately obtaining a set of material parameters. The analysis core consists of three modules: one for GVIPS material models, an analysis module containing a specialized finite-element solution algorithm, and an optimization module. COMPARE solves the problem of finding GVIPS material parameters in the manner of a design-optimization problem in which the parameters are the design variables.

Arnold, Steven M.↗

A Generalized Fluid System Simulation Program to Model Flow Distribution in Fluid Networks

This paper describes a general purpose computer program for analyzing steady state and transient flow in a complex network. The program is capable of modeling phase changes, compressibility, mixture thermodynamics and external body forces such as gravity and centrifugal. The program's preprocessor allows the user to interactively develop a fluid network simulation consisting of nodes and branches. Mass, energy and specie conservation equations are solved at the nodes; the momentum conservation equations are solved in the branches. The program contains subroutines for computing "real fluid" thermodynamic and thermophysical properties for 33 fluids. The fluids are: helium, methane, neon, nitrogen, carbon monoxide, oxygen, argon, carbon dioxide, fluorine, hydrogen, parahydrogen, water, kerosene (RP-1), isobutane, butane, deuterium, ethane, ethylene, hydrogen sulfide, krypton, propane, xenon, R-11, R-12, R-22, R-32, R-123, R-124, R-125, R-134A, R-152A, nitrogen trifluoride and ammonia. The program also provides the options of using any incompressible fluid with constant density and viscosity or ideal gas. Seventeen different resistance/source options are provided for modeling momentum sources or sinks in the branches. These options include: pipe flow, flow through a restriction, non-circular duct, pipe flow with entrance and/or exit losses, thin sharp orifice, thick orifice, square edge reduction, square edge expansion, rotating annular duct, rotating radial duct, labyrinth seal, parallel plates, common fittings and valves, pump characteristics, pump power, valve with a given loss coefficient, and a Joule-Thompson device. The system of equations describing the fluid network is solved by a hybrid numerical method that is a combination of the Newton-Raphson and successive substitution methods. This paper also illustrates the application and verification of the code by comparison with Hardy Cross method for steady state flow and analytical solution for unsteady flow.

Majumdar, Alok↗