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 415 records · Page 23

The thermal catastrophe model of substorms

Equations are presented for the physical mechanisms involved in the resonant absorption of Alfven waves in the plasma sheet boundary layer. It is shown that energy absorbed by the plasma-sheet particles is a function of the central plasma sheet temperature. The heating curve, when coupled with convective transport, yielded an equation of state for the steady state plasma sheet, whose solution has the form of a mathematical catastrophe. The master equation includes dynamic terms describing the transition across the thermal catastrophe, making it possible to evaluate the time scale for the catastrophe to occur.

Goertz, C. K.↗

Laboratory Simulations of Ammonia-Rich Oceans in Icy Worlds

Improved equations of state for ammonia-water solutions are important for properly understanding the interiors of large icy satellites hosting deep interior oceans. Titan is the primary example of such a world, but water-rich dwarf planets Measurements of solution density are now possible at relevant pressures (above approx. 250 megapascals) using the Simulator for Icy World Interiors. Analysis of sound velocity measurements in aqueous magnesium sulfate obtained in our laboratory, shows a correction on the order of 5% to 700 megapascals (7 kilobar) from -20 to 100 C and to 3 m (approx. 30 percentage by weight) concentration. Accurate prediction of density as a function of pressure, temperature, and ammonia concentration are needed for interpretation of remote observations to address questions of interior liquid layer depth, composition, and fluid dynamics.

equations of state↗

Energy-state formulation of lumped volume dynamic equations with application to a simplified free piston Stirling engine

Lumped volume dynamic equations are derived using an energy-state formulation. This technique requires that kinetic and potential energy state functions be written for the physical system being investigated. To account for losses in the system, a Rayleigh dissipation function is also formed. Using these functions, a Lagrangian is formed and using Lagrange's equation, the equations of motion for the system are derived. The results of the application of this technique to a lumped volume are used to derive a model for the free-piston Stirling engine. The model was simplified and programmed on an analog computer. Results are given comparing the model response with experimental data.

Daniele, C. J.↗

Dynamic simulation of coronal mass ejections

A model is developed for the formation and propagation through the lower corona of the loop-like coronal transients in which mass is ejected from near the solar surface to the outer corona. It is assumed that the initial state for the transient is a coronal streamer. The initial state for the streamer is a polytropic, hydrodynamic solution to the steady-state radial equation of motion coupled with a force-free dipole magnetic field. The numerical solution of the complete time-dependent equations then gradually approaches a stationary coronal streamer configuration. The streamer configuration becomes the initial state for the coronal transient. The streamer and transient simulations are performed completely independent of each other. The transient is created by a sudden increase in the pressure at the base of the closed-field region in the streamer configuration. Both coronal streamers and coronal transients are calculated for values of the plasma beta (the ratio of thermal to magnetic pressure) varying from 0.1 to 100.

Steinolfson, R. S.↗

Dynamical effects on the core-mantle boundary from depth-dependent thermodynamical properties of the lower mantle

A common assumption in modeling dynamical processes in the lower mantle is that both the thermal expansivity and thermal conductivity are reasonably constant. Recent work from seismic equation of state leads to substantially higher values for the thermal conductivity and much lower thermal expansivity values in the deep mantle. The dynamical consequences of incorporating depth-dependent thermodynamic properties on the thermal-mechanical state of the lower mantle are examined with the spherical-shell mean-field equations. It is found that the thermal structure of the seismically resolved anomalous zone at the base of the mantle is strongly influenced by these variable properties and, in particular, that the convective distortion of the core-mantle boundary (CMB) is reduced with the decreasing thermal expansivity. Such a reduction of the dynamically induced topography from pure thermal convection would suggest that some other dynamical mechanism must be operating at the CMB.

Zhang, Shuxia↗

Going Beyond Einstein with the Constellation-X Mission

The Constellation-X mission will address the questions: "What happens to matter close to a black hole?" and "What is Dark Energy?" These questions are central to the NASA Beyond Einstein Program, where Constellation-X plays a central role. The mission will address these questions by using high throughput X-ray spectroscopy to observe the effects of strong gravity close to the event horizon of black holes, and to observe the formation and evolution of clusters of galaxies in order to precisely determine Cosmological parameters. To achieve these primary science goals requires a factor of 25-100 increase in sensitivity for high resolution X-ray spectroscopy.'The mission will also perform routine high-resolution X-ray spectroscopy of faint 2nd extended X-ray source populations. This will provide diagnostic information such as density, elemental abundances, velocity; and ionization state for a wide range of astrophysical problems, including new constraints on the Neutron Star equation of state.

White, Nicholas↗

Relative Shock Effects in Mixed Powders of Calcite, Gypsum, and Quartz: A Calibration Scheme from Shock Experiments

The shock behavior of calcite and gypsum is important in understanding the Cretaceous/Tertiary event and other terrestrial impacts that contain evaporite sediments in their targets. Most interest focuses on issues of devolatilization to quantify the production of CO2 or SO2 to better understand their role in generating a temporary atmosphere and its effects on climate and biota [e.g., papers in 1,2,3,4]. Devolatilization of carbonate is also important because the dispersion and fragmentation of ejecta is strongly controlled by the expansion of large volumes of gas during the impact process as well [5,6]. Shock recovery experiments for calcite yield seemingly conflicting results: early experimental devolatilization studies [7,8,9] suggested that calcite was substantially outgassed at 30 GPa (> 50%). However, the recent petrographic work of [10,11,12] presented evidence that essentially intact calcite is recovered from 60 GPa experiments. [13] reported results of shock experiments on anhydrite, gypsum, and mixtures of those phases with silica. Their observations indicate little or no devolatilization of anhydrite shocked to 42 GPa and that the fraction of sulfur, by mass, that degassed is approx.10(exp -2) of theoretical prediction. In another (preliminary) report of shock experiments on calcite, anhydrite, and gypsum, [14] observe calcite recrystallization when shock loaded at 61 GPa, only intensive plastic deformation in anhydrite shock loaded at 63 GPa, and gypsum converted to anhydrite when shock loaded at 56 GPa. [15] shock loaded anhydrite and quartz to a peak pressure of 60 GPa. All of the quartz grains were trans-formed to glass and the platy anhydrite grains were completely pseudomorphed by small crystallized anhydrite grains. However, no evidence of interaction between the two phases could be observed and they suggest that recrystallization of anhydrite grains is the result of a solid state transformation. [16] reanalyzed the calcite and anhydrite shock wave experiments of [17] using improved equations of state of porous materials and vaporized products. They determined the pressures for incipient and complete vaporization to be 32.5 and 122 GPa for anhydrite and 17.8 and 54.1 GPa for calcite, respectively, a factor of 2 to 3 lower than reported earlier by [17].

Bell, Mary S.↗

Nodal Modeling of Liquid Propellant Feed and Pressurization System

Nodal or network flow modeling plays an important role in the design and operation of the feed and pressurization system of a liquid rocket engine. Model development and execution time is relatively short for nodal codes in comparison to the Navier-Stokes based CFD codes. Nodal models also allow the inclusion of several components into one model to predict the behavior of a larger system. Unlike CFD models, the nodal models do not need very fine discretization of the flow field because they use empirical correlations to model fluid friction and heat transfer. This paper presents several applications of nodal modeling of liquid propellant feed and pressurization systems using the Generalized Fluid System Simulation Program (GFSSP), a nodal code developed at NASA/ Marshall Space Flight Center. GFSSP discretizes the flow field into nodes which are connected by branches. The mass and energy conservation equations and the equation of state are solved to calculate pressure, temperature, and resident mass at the nodes. The momentum equations are solved at the branches to calculate flow rates. Applications include a) tank pressurization by inert gas as well as autogenous pressurization by gaseous propellant, b) submerged pressurization by helium, c) self-pressurization due to boil-off of cryogenic propellant, d) chilldown of a transfer line of a cryogenic propellant feed system, and e) chilldown and filling of a cryogenic tank. Each of the above-mentioned models was verified and validated by comparing with test data.

Nodal Model↗

Nodal Modeling of Liquid Propellant Feed and Pressurization System

Nodal or network flow modeling plays an important role in the design and operation of the feed and pressurization system of a liquid rocket engine. Model development and execution time is relatively short for nodal codes in comparison to the Navier-Stokes based CFD codes. Nodal models also allow the inclusion of several components into one model to predict the behavior of a larger system. Unlike CFD models, the nodal models do not need very fine discretization of the flow field because they use empirical correlations to model fluid friction and heat transfer. This paper presents several applications of nodal modeling of liquid propellant feed and pressurization systems using the Generalized Fluid System Simulation Program (GFSSP), a nodal code developed at NASA/ Marshall Space Flight Center. GFSSP discretizes the flow field into nodes which are connected by branches. The mass and energy conservation equations and the equation of state are solved to calculate pressure, temperature, and resident mass at the nodes. The momentum equations are solved at the branches to calculate flow rates. Applications include a) tank pressurization by inert gas as well as autogenous pressurization by gaseous propellant, b) submerged pressurization by helium, c) self-pressurization due to boil-off of cryogenic propellant, d) chilldown of a transfer line of a cryogenic propellant feed system, and e) chilldown and filling of a cryogenic tank. Each of the above-mentioned models was verified and validated by comparing with test data.

Nodal Model↗

Microscopic Simulation and Macroscopic Modeling for Thermal and Chemical Non-Equilibrium

This paper deals with the accurate microscopic simulation and macroscopic modeling of extreme non-equilibrium phenomena, such as encountered during hypersonic entry into a planetary atmosphere. The state-to-state microscopic equations involving internal excitation, de-excitation, dissociation, and recombination of nitrogen molecules due to collisions with nitrogen atoms are solved time-accurately. Strategies to increase the numerical efficiency are discussed. The problem is then modeled using a few macroscopic variables. The model is based on reconstructions of the state distribution function using the maximum entropy principle. The internal energy space is subdivided into multiple groups in order to better describe the non-equilibrium gases. The method of weighted residuals is applied to the microscopic equations to obtain macroscopic moment equations and rate coefficients. The modeling is completely physics-based, and its accuracy depends only on the assumed expression of the state distribution function and the number of groups used. The model makes no assumption at the microscopic level, and all possible collisional and radiative processes are allowed. The model is applicable to both atoms and molecules and their ions. Several limiting cases are presented to show that the model recovers the classical twotemperature models if all states are in one group and the model reduces to the microscopic equations if each group contains only one state. Numerical examples and model validations are carried out for both the uniform and linear distributions. Results show that the original over nine thousand microscopic equations can be reduced to 2 macroscopic equations using 1 to 5 groups with excellent agreement. The computer time is decreased from 18 hours to less than 1 second.

Modeling↗

Collisionless solar wind. II - Variable electron temperature.

Examination of a two-component ?model' for the solar wind, in which the protons become collisionless beyond a distance from the solar center equal to or greater than 10 solar radii, where they are already highly supersonic. The proton temperatures are found from the double adiabatic equation of state. The electrons are highly subsonic, and their temperature profile is prescribed ad hoc. The momentum equations for the electrons and protons are solved subject to the conditions of quasi-neutrality and zero charge efflux from the sun. Some of the principal results are: (1) the proton thermal anisotropy is substantially reduced when solar rotation is considered; (2) solar rotation leads to significantly lower mean proton temperatures; and (3) the electron temperature profile in the supersonic region is the primary parameter determining flow acceleration there.

Hollweg, J. V.↗

Buckling of cylindrical shells with random imperfections

Buckling of long cylindrical shells with random imperfections subjected to axial loads is treated using the method of truncated hierarchy. A system of homogeneous variational equations is set up to examine the existence of bifurcation in the neighborhood of the equilibrium state. These equations are linear and involve stationary random coefficients. The truncated-hierarchy method is applied, and characteristic equations are obtained. Various exponential cosine correlation functions associated with asymmetric imperfections are examined numerically. Qualitatively, the results obtained are as anticipated.

Fersht, R. S.↗

Polynomial elimination theory and non-linear stability analysis for the Euler equations

Numerical methods are presented that exploit the polynomial properties of discretizations of the Euler equations. It is noted that most finite difference or finite volume discretizations of the steady-state Euler equations produce a polynomial system of equations to be solved. These equations are solved using classical polynomial elimination theory, with some innovative modifications. This paper also presents some preliminary results of a new non-linear stability analysis technique. This technique is applicable to determining the stability of polynomial iterative schemes. Results are presented for applying the elimination technique to a one-dimensional test case. For this test case, the exact solution is computed in three iterations. The non-linear stability analysis is applied to determine the optimal time step for solving Burgers' equation using the MacCormack scheme. The estimated optimal time step is very close to the time step that arises from a linear stability analysis.

Kennon, S. R.↗

Hydrodynamics of coalescing binary neutron stars: Ellipsoidal treatment

We employ an approximate treatment of dissipative hydrodynamics in three dimensions to study the coalescence of binary neutron stars driven by the emission of gravitational waves. The stars are modeled as compressible ellipsoids obeying a polytropic equation of state; all internal fluid velocities are assumed to be linear functions of the coordinates. The hydrodynamics equations then reduce to a set of coupled ordinary differential equations for the evolution of the principal axes of the ellipsoids, the internal velocity parameters, and the binary orbital parameters. Gravitational radiation reaction and viscous dissipation are both incorporated. We set up exact initial binary equilibrium configurations and follow the transition from the quasi-static, secular decay of the orbit at large separation to the rapid dynamical evolution of the configurations just prior to contact. A hydrodynamical instability resulting from tidal interactions significantly accelerates the coalescence at small separation, leading to appreciable radial infall velocity and tidal lag angles near contact. This behavior is reflected in the gravitational waveforms and may be observable by gravitational wave detectors under construction. In cases where the neutron stars have spins which are not aligned with the orbital angular momentum, the spin-induced quadrupole moment can lead to precession of the orbital plane and therefore modulation of the gravitational wave amplitude even at large orbital radius. However, the amplitude of the modulation is small for typical neutron star binaries with spins much smaller than the orbital angular momentum.

Lai, Dong↗

Two-Dimensional Network Flow Modeling of No-Vent Tank Filling of A Cryogenic Tank With Thermo-Dynamic Vent System Assisted Injector

A finite volume-based network flow modeling tool, Generalized Fluid System Simulation Program (GFSSP), was used to simulate the no-vent filling of a cryogenic tank with a thermodynamic vent system assisted injector. The tank was discretized into nodes in the radial and axial directions to calculate two-dimensional axisymmetric flow with heat transfer between the solid wall and fluid. The nodes are connected by branches where momentum equations are solved to calculate the flowrates. The mass and energy conservation equations and the equation of state are solved to calculate pressure, temperature, and resident mass at the nodes. The energy conservation equations at the solid nodes were solved to calculate the temperature of the solid wall. The system of equations was solved by a combination of successive substitution and simultaneous Newton-Raphson method. The numerical model accounts for a) different regimes of pool boiling heat transfer, b) condensation of vapor around the spray droplets, c) condensation of vapor at the interface of the cooled injector and ullage, and d) condensation at the liquid-vapor interphase. The numerical predictions of tank pressure, filling rate, and wall temperatures were compared with the test data. The paper also presents the effect of nodal discretization by comparing the results of the 1D and 2D models and examines the effect of the droplet diameter on tank filling.

Nodal Model↗

Two-Dimensional Network Flow Modeling of No-Vent Tank Filling of A Cryogenic Tank With Thermo-Dynamic Vent System Assisted Injector

A finite volume-based network flow modeling tool, Generalized Fluid System Simulation Program (GFSSP), was used to simulate the no-vent filling of a cryogenic tank with a thermodynamic vent system assisted injector. The tank was discretized into nodes in the radial and axial directions to calculate two-dimensional axisymmetric flow with heat transfer between the solid wall and fluid. The nodes are connected by branches where momentum equations are solved to calculate the flowrates. The mass and energy conservation equations and the equation of state are solved to calculate pressure, temperature, and resident mass at the nodes. The energy conservation equations at the solid nodes were solved to calculate the temperature of the solid wall. The system of equations was solved by a combination of successive substitution and simultaneous Newton-Raphson method. The numerical model accounts for a) different regimes of pool boiling heat transfer, b) condensation of vapor around the spray droplets, c) condensation of vapor at the interface of the cooled injector and ullage, and d) condensation at the liquid-vapor interphase. The numerical predictions of tank pressure, filling rate, and wall temperatures were compared with the test data. The paper also presents the effect of nodal discretization by comparing the results of the 1D and 2D models and examines the effect of the droplet diameter on tank filling.

Nodal Model↗

An approximate factorization solution of the Navier-Stokes equations for transonic flow using body-fitted coordinates with application to NACA 64A010 airfoils

The implementation of the approximate factorization algorithm and its ability to efficiently and accurately describe transonic flow about an NACA 64A010 airfoil section is examined. The approximate factorization algorithm is developed from the nondimensional, conservative, vectorized Navier-Stokes equations expressed in curvilinear coordinates. Equations of state and transport coefficient relations appropriate to atmospheric air are appended to close the system of partial differential equations. An algebraic turbulence model is also incorporated into the equation set. This algorithm was verified by investigating the flow about an NACA 64A010 airfoil at 0, 2, and 3.5 deg angle of attack for free-stream conditions of 2,000,000 Reynolds number and 0.8 Mach number. Overall results were in good qualitative agreement with wind tunnel data sets. However, while nondimensional times of six were attained, numerical difficulties prevented any case from reaching a true steady state.

Copper, G. K.↗