Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Equation of State”

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 379 records · Page 21

Acoustic transmission matrix of a variable area duct or nozzle carrying a compressible subsonic flow

The differential equations governing the propagation of sound in a variable area duct or nozzle carrying a one dimensional subsonic compressible fluid flow are derived and put in state variable form using acoustic pressure and particle velocity as the state variables. The duct or nozzle is divided into a number of regions. The region size is selected so that in each region the Mach number can be assumed constant and the area variation can be approximated by an exponential area variation. Consequently, the state variable equation in each region has constant coefficients. The transmission matrix for each region is obtained by solving the constant coefficient acoustic state variable differential equation. The transmission matrix for the duct or nozzle is the product of the individual transmission matrices of each region. Solutions are presented for several geometries with and without mean flow.

Miles, J. H.↗

Steady-state and transient analysis of a squeeze film damper bearing for rotor stability

A study of the steady-state and transient response of the squeeze film damper bearing is presented. Both the steady-state and transient equations for the hydrodynamic bearing forces are derived. The bearing equivalent stiffness and damping coefficients are determined by steady-state equations. These coefficients are used to find the bearing configuration which will provide the optimum support characteristics based on a stability analysis of the rotor-bearing system. The transient analysis of rotor-bearing systems is performed by coupling the bearing and journal equations and integrating forward in time. The effects of unbalance, cavitation, and retainer springs are included in the analysis. Methods of determining the stability of a rotor-bearing system under the influence of aerodynamic forces and internal shaft friction are discussed with emphasis on solving the system characteristic frequency equation and on producing stability maps. It is shown that for optimum stability and low force transmissability the squeeze bearing should operate at an eccentricity ratio epsilon 0.4.

Barrett, L. E.↗

Acoustic transmission matrix of a variable area duct or nozzle carrying a compressible subsonic flow

The differential equations governing the propagation of sound in a variable area duct or nozzle carrying a one-dimensional subsonic compressible fluid flow are derived and put in state variable form using acoustic pressure and particle velocity as the state variables. The duct or nozzle is divided into a number of regions. The region size is selected so that in each region the Mach number can be assumed constant and the area variation can be approximated by an exponential area variation. Consequently, the state variable equation in each region has constant coefficients. The transmission matrix for each region is obtained by solving the constant coefficient acoustic state variable differential equation. The transmission matrix for the duct or nozzle is the product of the individual transmission matrices of each region. Solutions are presented for several geometries with and without mean flow.

Miles, J. H.↗

Surfactant Facilitated Spreading of Aqueous Drops on Hydrophobic Surfaces

Microgravity technologies often require aqueous phases to spread over nonwetting hydrophobic solid surfaces. Surfactants facilitate the wetting of water on hydrophobic surfaces by adsorbing on the water/air and hydrophobic solid/water interfaces and lowering the surface tensions of these interfaces. The tension reductions decrease the contact angle, which increases the equilibrium wetted area. Hydrocarbon surfactants; (i.e., amphiphiles with a hydrophobic moiety consisting of an extended chain of (aliphatic) methylene -CH2- groups attached to a large polar group to give aqueous solubility) are capable of reducing the contact angles on surfaces which are not very hydrophobic, but do not reduce significantly the contact angles of the very hydrophobic surfaces such as parafilm, polyethylene or self assembled monolayers. Trisiloxane surfactants (amphiphiles with a hydrophobe consisting of methyl groups linked to a trisiloxane backbone in the form of a disk ((CH3)3-Si-O-Si-O-Si(CH3)3) and an extended ethoxylate (-(OCH2CH2)a-) polar group in the form of a chain with four or eight units) can significantly reduce the contact angle of water on a very hydrophobic surface and cause rapid and complete (or nearly complete) spreading (termed superspreading). The overall goal of the research described in this proposal is to establish and verify a theory for how trisiloxanes cause superspreading, and then use this knowledge as a guide to developing more general hydrocarbon based surfactant systems which superspread. We propose that the trisiloxane surfactants superspread because their structure allows them to strongly lower the high hydrophobic solid/aqueous tension when they adsorb to the solid surface. When the siloxane adsorbs, the hydrophobic disk parts of the molecule adsorb onto the surface removing the surface water. Since the cross-sectional area of the disk is larger than that of the extended ethoxylate chain, the disks can form a space-filling mat on the surface which removes a significant amount of the surface water. In this presentation, we report the results of measurements of the molecular packing and rates of kinetic exchange of the trisiloxane surfactants at the air/water interface in order to confirm our picture of trisiloxane packing, and provide additional insight into the superspreading process. We used the pendant bubble technique as a Langmuir trough to measure the trisiloxane equation of state which relates the tension to the surface concentration. From these measurements we obtain accurate values for the maximum packing density. We find that trisiloxanes with 4 and 8 ethoxylate groups have the same maximum packing concentration, indicating that the maximum packing is controlled by the cross section of the head group. For trisiloxanes with larger than eight ethoxylates, the maximum packing increases with ethoxylate number, indicating that the disposition of the ethoxylate chain (i.e., its effective size) is controlling. This supports our picture of superspreading: The superspreading ability of trisiloxanes decreases considerably for trisiloxanes with larger than eight ethoxylates; the packing measurements indicate that for the higher ethoxylate number trisiloxanes, the compact nonpolar head groups are pushed apart by the ethoxylate chain. They leave spaces of surface water on adsorption and do not lower the solid tension as much as their lower chain analogues. Finally the report measurements of the dynamic tension reduction accompanying the adsorption of trisiloxanes onto an initially clean interface using the pendant bubble technique, and we obtain from these relaxations, the equation of state and a mass transfer model, the rate constants for kinetic exchange. We find that the rate constants for desorption of trisiloxanes are generally much slower than for analogous aliphatic polyethoxylate surfactants with identical ethoxylate chain lengths. When an aqueous drop of a superspreader solution is placed on a hydrophobic surface and begins to spread, the lower desorption rates allows the tension at the drop center to remain reduced relative to the tension of the expanding periphery, thus strengthening Marangoni forces which can assist the spreading. Marangoni forces can be especially significant in the case of superspreaders because their maximum equilibrium reduction in tension is several dynes/cm lower than for aliphatic surfactants.

Kumar, Nitin↗

Tethers as Debris: Simulating Impacts of Kevlar Tethers on Shuttle Tiles

In a previous paper I examined the effects of impacts of polymer tethers on aluminum plates using the SPHC hydrodynamic code. In this paper I apply tether models to a new target - models of Space Shuttle tiles developed during the STS 107 accident investigation. In this three-dimensional simulation, a short tether fragment strikes a single tile supported on an aluminum backing plate. A tile of the LI-900 material is modeled. Penetration and damage to the tile and the backwall are characterized for three normal impact velocities. The tether is modeled as a bundle of eight 1-mm strands, with the bundle having dimensions 2-mm x 4-mm x 20-cm. The bulk material properties used are those of Kevlar(TradeMark) 49, for which a Mie-Gruneisen multiphase equation of state (eos) is used. In addition, the strength model is applied in a linear sense, such that tensile loads along the strand length are supported, but there is no strength in the lateral directions. Tile models include the various layers making up the tile structure. The outermost layer is a relatively dense borosilicate glass, known as RCG, 0.5-mm thick. The RCG layer is present on the top and four sides of the tile. Below this coating is the bulk of the tile, 1.8- in thick, made of LI-900, a product consisting of rigidized fiberous silica with a density of 9 lWft3. Below the main insulating layer is a bottom layer of the same material that has been treated to increase its density by approximately 69% to improve its strength. This densified layer is bonded to a Strain Isolation Pad (SIP), modeled as a refractory felt fabric. The SIP is bonded to an aluminum 2024 wall 0.1-in thick. The tile and backwall materials use a Me-Gruneisen multiphase eos, with the exception of the SIP felt, which uses a fabric equation of state. Fabrics must be crushed to the full bulk material density before bulk material properties and a Mie-Gruneisen eos are applied. Tether fragment impact speeds of 3,7, and 10 km/s are simulated, with impact velocities normal to the tile face. Damage results are presented in tabular format.

Evans, Steven W.↗

A finite element model of conduction, convection, and phase change near a solid/melt interface

Detailed understanding of heat transfer and fluid flow is required for many aerospace thermal systems. These systems often include phase change and operate over a range of accelerations or effective gravitational fields. An approach to analyzing such systems is presented which requires the simultaneous solution of the conservation laws of energy, momentum, and mass, as well as an equation of state. The variable property form of the governing equations are developed in two-dimensional Cartesian coordinates for a Newtonian fluid. A numerical procedure for solving the governing equations is presented and implemented in a computer program. The Galerkin form of the finite element method is used to solve the spatial variation of the field variables, along with the implicit Crank-Nicolson time marching algorithm. Quadratic Langrangian elements are used for the internal energy and the two components of velocity. Linear Lagrangian elements are used for the pressure. The location of the solid/liquid interface as well as the temperatures are determined form the calculated internal energy and pressure. This approach is quite general in that it can describe heat transfer without phase change, phase change with a sharp interface, and phase change without an interface. Analytical results from this model are compared to those of other researchers studying transient conduction, convection, and phase change and are found to be in good agreement. The numerical procedure presented requires significant computer resources, but this is not unusual when compared to similar studies by other researchers. Several methods are suggested to reduce the computational times.

Viterna, Larry A.↗

The computation of standard solar models

Procedures for calculating standard solar models with the usual simplifying approximations of spherical symmetry, no mixing except in the surface convection zone, no mass loss or gain during the solar lifetime, and no separation of elements by diffusion are described. The standard network of nuclear reactions among the light elements is discussed including rates, energy production and abundance changes. Several of the equation of state and opacity formulations required for the basic equations of mass, momentum and energy conservation are presented. The usual mixing-length convection theory is used for these results. Numerical procedures for calculating the solar evolution, and current evolution and oscillation frequency results for the present sun by some recent authors are given.

Ulrich, Roger K.↗

Optimization methods, flux conserving methods for steady state Navier-Stokes equation

Navier-Stokes equation as discretized by new flux conserving method proposed by Chang and Scott results in the system: vector F(vector x) = 0, where F is a vector valued function. The Optimization method we use is based on Quasi-Newton methods: given a nonlinear function vector F(vector x) = 0, we solve, Delta(vector x) = -BF(vector x), where Delta(vector x) is the correction term and B is the inverse Jacobian of F(x). Then, iteratively, vector(x(sub (i+1))) = vector(x (sub i)) + alpha.Delta(vector x(sub i)), where alpha is a line search correction term determined by a line search routine. We use the BFCG's update the Jacobian matrix B(sub k) at each iteration. It is well known that B(sub k) approaches B(*) at the solution X(*). This algorithm has several advantages over the Newton-Raphson method. For example, we do not need to calculate the Jacobian matrix at each iteration which is computationally very expensive.

Adeyeye, John↗

A Second Law Based Unstructured Finite Volume Procedure for Generalized Flow Simulation

An unstructured finite volume procedure has been developed for steady and transient thermo-fluid dynamic analysis of fluid systems and components. The procedure is applicable for a flow network consisting of pipes and various fittings where flow is assumed to be one dimensional. It can also be used to simulate flow in a component by modeling a multi-dimensional flow using the same numerical scheme. The flow domain is discretized into a number of interconnected control volumes located arbitrarily in space. The conservation equations for each control volume account for the transport of mass, momentum and entropy from the neighboring control volumes. In addition, they also include the sources of each conserved variable and time dependent terms. The source term of entropy equation contains entropy generation due to heat transfer and fluid friction. Thermodynamic properties are computed from the equation of state of a real fluid. The system of equations is solved by a hybrid numerical method which is a combination of simultaneous Newton-Raphson and successive substitution schemes. The paper also describes the application and verification of the procedure by comparing its predictions with the analytical and numerical solution of several benchmark problems.

Majumdar, Alok↗

Quasi-Newton methods for parameter estimation in functional differential equations

A state-space approach to parameter estimation in linear functional differential equations is developed using the theory of linear evolution equations. A locally convergent quasi-Newton type algorithm is applied to distributed systems with particular emphasis on parameters that induce unbounded perturbations of the state. The algorithm is computationally implemented on several functional differential equations, including coefficient and delay estimation in linear delay-differential equations.

Brewer, Dennis W.↗

Real-gas effects 1: Simulation of ideal gas flow by cryogenic nitrogen and other selected gases

The thermodynamic properties of nitrogen gas do not thermodynamically approximate an ideal, diatomic gas at cryogenic temperatures. Choice of a suitable equation of state to model its behavior is discussed and the equation of Beattie and Bridgeman is selected as best meeting the needs for cryogenic wind tunnel use. The real gas behavior of nitrogen gas is compared to an ideal, diatomic gas for the following flow processes: isentropic expansion; normal shocks; boundary layers; and shock wave boundary layer interactions. The only differences in predicted pressure ratio between nitrogen and an ideal gas that may limit the minimum operating temperatures of transonic cryogenic wind tunnels seem to occur at total pressures approaching 9atmospheres and total temperatures 10 K below the corresponding saturation temperature, where the differences approach 1 percent for both isentropic expansions and normal shocks. Several alternative cryogenic test gases - air, helium, and hydrogen - are also analyzed. Differences in air from an ideal, diatomic gas are similar in magnitude to those of nitrogen. Differences for helium and hydrogen are over an order of magnitude greater than those for nitrogen or air. Helium and hydrogen do not approximate the compressible flow of an ideal, diatomic gas.

Hall, R. M.↗

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↗

Analytic Steady-State Accuracy of a Spacecraft Attitude Estimator

This paper extends Farrenkopf's analysis of a single-axis spacecraft attitude estimator using gyro and angle sensor data to include the angle output white noise of a rate-integrating gyro. Analytic expressions are derived for the steady-state pre-update and post-update angle and drift bias variances and for the state update equations. It is shown that only part of the state update resulting from the angle sensor measurement is propagated to future times.

Markley, F. Landis↗

Computation techniques and computer programs to analyze Stirling cycle engines using characteristic dynamic energy equations

The basic equations that are used to describe the physical phenomena in a Stirling cycle engine are the general energy equations and equations for the conservation of mass and conversion of momentum. These equations, together with the equation of state, an analytical expression for the gas velocity, and an equation for mesh temperature are used in this computer study of Stirling cycle characteristics. The partial differential equations describing the physical phenomena that occurs in a Stirling cycle engine are of the hyperbolic type. The hyperbolic equations have real characteristic lines. By utilizing appropriate points along these curved lines the partial differential equations can be reduced to ordinary differential equations. These equations are solved numerically using a fourth-fifth order Runge-Kutta integration technique.

Larson, V. H.↗

An optimal momentum management controller for the Space Station

This paper presents a new sequential design procedure for determining an optimal CMG (control moment gyro) momentum management and attitude control system for the Space Station Freedom. First, the Space Station equations of motion are linearized and uncoupled, and the associated state space equations are defined. Next, a new sequential procedure is used for the development of a continuous quadratic regulator with eigenvalue placement in a specified region of the complex plane. The regional pole assignment method is utilized since it is best suited for tradeoffs between eigenvalue locations and robustness with respect to parameter variations, sensor failures, implementation accuracies and gain reductions. The matrix sign function is used for solving the Riccati equations which appear in the design procedure. Simulation results are given which show that the resultant design provides desired system performance.

Sunkel, J. W.↗