Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite density”

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

Use of shape-preserving interpolation methods in surface modeling

In many large-scale scientific computations, it is necessary to use surface models based on information provided at only a finite number of points (rather than determined everywhere via an analytic formula). As an example, an equation of state (EOS) table may provide values of pressure as a function of temperature and density for a particular material. These values, while known quite accurately, are typically known only on a rectangular (but generally quite nonuniform) mesh in (T,d)-space. Thus interpolation methods are necessary to completely determine the EOS surface. The most primitive EOS interpolation scheme is bilinear interpolation. This has the advantages of depending only on local information, so that changes in data remote from a mesh element have no effect on the surface over the element, and of preserving shape information, such as monotonicity. Most scientific calculations, however, require greater smoothness. Standard higher-order interpolation schemes, such as Coons patches or bicubic splines, while providing the requisite smoothness, tend to produce surfaces that are not physically reasonable. This means that the interpolant may have bumps or wiggles that are not supported by the data. The mathematical quantification of ideas such as physically reasonable and visually pleasing is examined.

Ftitsch, F. N.↗

Calculation of real-gas effects on blunt-body trim angles

The effect of vibrational excitation and dissociation at high temperatures on the trim angle of attack of a blunt lifting body is calculated for a nonequilibrium flow regime in air using a CFD technique. The vibrational-electronic temperature and the species densities are calculated assuming the flow to be in a nonequilibrium state. The forebody flow of a two-dimensional blunt body of the shape of the Apollo Command Module at a finite angle of attack is calculated. The results show that the pitching moment around a reference point is larger and the trim angle of attack is smaller for a reacting gas than for a perfect gas. The calculated shift in the trim angle due to the real-gas effect is of the same order as that seen during the Apollo flights.

Park, Chul↗

Calculations of the flow properties of a confined diffusion flame

A finite element algorithm for the computation of confined, axisymmetric, turbulent diffusion flames is developed. The mean mixture properties were obtained by three methods based on diffusion flame concept: without using a probability density function (PDF), with a double-delta PDF, and with a beta PDF. A comparison is made for the combustion models, and the effect of turbulence on combustion are discussed.

Kim, Yongmo↗

The study of PDF turbulence models in combustion

In combustion computations, it is known that the predictions of chemical reaction rates are poor if conventional turbulence models are used. The probability density function (pdf) method seems to be the only alternative that uses local instantaneous values of the temperature, density, etc., in predicting chemical reaction rates, and thus is the only viable approach for more accurate turbulent combustion calculations. The fact that the pdf equation has a very large dimensionality renders finite difference schemes extremely demanding on computer memories and thus impractical. A logical alternative is the Monte Carlo scheme. Since CFD has a certain maturity as well as acceptance, it seems that the use of a combined CFD and Monte Carlo scheme is more beneficial. Therefore, a scheme is chosen that uses a conventional CFD flow solver in calculating the flow field properties such as velocity, pressure, etc., while the chemical reaction part is solved using a Monte Carlo scheme. The discharge of a heated turbulent plane jet into quiescent air was studied. Experimental data for this problem shows that when the temperature difference between the jet and the surrounding air is small, buoyancy effect can be neglected and the temperature can be treated as a passive scalar. The fact that jet flows have a self-similar solution lends convenience in the modeling study. Futhermore, the existence of experimental data for turbulent shear stress and temperature variance make the case ideal for the testing of pdf models wherein these values can be directly evaluated.

Hsu, Andrew T.↗

Probability density function approach for compressible turbulent reacting flows

The objective of the present work is to extend the probability density function (PDF) tubulence model to compressible reacting flows. The proability density function of the species mass fractions and enthalpy are obtained by solving a PDF evolution equation using a Monte Carlo scheme. The PDF solution procedure is coupled with a compression finite-volume flow solver which provides the velocity and pressure fields. A modeled PDF equation for compressible flows, capable of treating flows with shock waves and suitable to the present coupling scheme, is proposed and tested. Convergence of the combined finite-volume Monte Carlo solution procedure is discussed. Two super sonic diffusion flames are studied using the proposed PDF model and the results are compared with experimental data; marked improvements over solutions without PDF are observed.

Hsu, A. T.↗

Informatics and High Throughput Screening of Thermophysical Properties

The combination of computer-aided experiments with computational modeling enables a new class of powerful tools for materials research. A non-contact method for measuring density, thermal expansion, and creep of undercooled and high-temperature materials has been developed, using electrostatic levitation and optical diagnostics, including digital video. These experiments were designed to take advantage of the large volume of data (many gigabytes/experiment, terabytes/campaign) to gain additional information about the samples. For example, using sub-pixel interpolation to measure about 1000 vectors per image of the sample's surface allows the density of an axisymmetric sample to be determined to an accuracy of about 200 ppm (0.02%). A similar analysis applied to the surface shape of a rapidly rotating sample is combined with finite element modeling to determine the stress-dependence of creep in the sample in a single test. Details of the methods for both the computer-aided experiments and computational models will be discussed.

Hyers, Robert W.↗

Experimental Determination of in Situ Utilization of Lunar Regolith for Thermal Energy Storage

A Lunar Thermal Energy from Regolith (LUTHER) experiment has been designed and fabricated at the NASA Lewis Research Center to determine the feasibility of using lunar soil as thermal energy storage media. The experimental apparatus includes an alumina ceramic canister (25.4 cm diameter by 45.7 cm length) which contains simulated lunar regolith, a heater (either radiative or conductive), 9 heat shields, a heat transfer cold jacket, and 19 type B platinum rhodium thermocouples. The simulated lunar regolith is a basalt, mined and processed by the University of Minnesota, that closely resembles the lunar basalt returned to earth by the Apollo missions. The experiment will test the effects of vacuum, particle size, and density on the thermophysical properties of the regolith. The properties include melt temperature (range), specific heat, thermal conductivity, and latent heat of storage. Two separate tests, using two different heaters, will be performed to study the effect of heating the system using radiative and conductive heat transfer. The physical characteristics of the melt pattern, material compatibility of the molten regolith, and the volatile gas emission will be investigated by heating a portion of the lunar regolith to its melting temperature (1435 K) in a 10(exp -4) pascal vacuum chamber, equipped with a gas spectrum analyzer. A finite differencing SINDA model was developed at NASA Lewis Research Center to predict the performance of the LUTHER experiment. The analytical results of the code will be compared with the experimental data generated by the LUTHER experiment. The code will predict the effects of vacuum, particle size, and density has on the heat transfer to the simulated regolith.

Richter, Scott W.↗

Effects of a magnetic field gradient on the lower hybrid drift instability

The effect of a magnetic field gradient on the lower hybrid drift instability is examined in the low drift velocity regime (i.e., of the order of the ion thermal velocity). The theoretical model includes density and temperature gradients, and considers only waves propagating perpendicular to the ambient magnetic field with wavelengths of the order of the electron Larmor radius. It is found, for the case of a finite ratio of the electron pressure to the magnetic field pressure, that a magnetic gradient drift-wave resonance can be important and is stabilizing or destabilizing depending upon the physical conditions. Substantial plasma heating and a large anomalous resistivity can also occur.

Huba, J. D.↗

An instability of finite amplitude circularly polarized Alfven waves

A demonstration is presented that a finite-amplitude circularly-polarized Alfven wave is generally unstable in a MHD fluid. The wave decays by a four-wave coupling process in which the daughter waves are forward propagating random density and magnetic fluctuations and a backward-propagating magnetic wave. For parameters typical of the solar corona and the solar wind (thermal to magnetic energy density ratios between 0.1 and 1, and values between 0.1 and 0.9 for the ratio of magnetic energy density of the initial Alfven wave to that of the background magnetic field), large decay rates are found.

Goldstein, M. L.↗

On the varying slope of velocity spectra

Spectra of zonal, meridional and vertical wind velocity, measured during a 24 hour period with the spaced-antenna technique indicate quite a variable slope as a function of height. It is found that the spectral slope (1h to 24h) of all three components correlates with the mean horizontal wind velocity. A possible conclusion is that the frequency dependence of power density of horizontal and vertical fluctuation component apparently depends on the mean wind velocity. However, the vertical spectra at periods larger than about 1 hour can also be influenced by spillover (due to finite radar antenna beam width) from the horizontal fluctuation component or by a Doppler shift.

Rottger, J.↗

Upwind Navier-Stokes solutions for leading-edge vortex flows

An incompressible Navier-Stokes solver using an upwind finite-difference algorithm is employed to investigate low-speed, three-dimensional, laminar, leading-edge vortex flows over three round-edged low-aspect-ratio wings. The effects of grid density, Reynolds number, and wing planform on the flowfield structures and integral values are studied. Computed results show good qualitative and quantitative agreement with the available experimental data.

Hsu, C.-H.↗

Entropic uncertainty relation at finite temperature

We discussed how much information is lost when a particle is in equilibrium with the thermal reservoir of temperature T = 1/beta. The universal temperature correction to the r.h.s. of U(X,P:psi) greater than or = 1 + ln(pi) is determined. For this purpose, it is convenient to employ the framework of thermo-field dynamics (TFD). This formulation of finite-temperature (T not = 0) quantum theory utilizes the doubled Hilbert space, the normal operator (A) acting on the objective space, and its corresponding tildian operator on the fictitious space. The physical probability density associated with the measurement of the normal operator, A, is given, and the information entropy at T not = 0 is defined. The results describe how the thermal disturbance effects in S sub X or S sub P (delta X or delta P) can be suppressed by squeezing with keeping U = S sub X + S sub P (delta X x delta P) its minimum value.

Abe, Sumiyoshi↗

Solution strategies and heat transfer calculations for three-dimensional configurations at hypersonic speeds

A procedure which reduces the memory requirements for computing the viscous flow over a modified Orbiter geometry at a hypersonic flight condition is presented. The Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) code which incorporates a thermochemical nonequilibrium chemistry model, a finite rate catalytic wall boundary condition and wall temperature distribution based on radiation equilibrium is used in this study. In addition, the effect of choice of 'min mod' function, eigenvalue limiter and grid density on surface heating is investigated. The surface heating from a flowfield calculation at Mach number 22, altitude of 230,000 ft and 40 deg angle of attack is compared with flight data from three Orbiter flights.

Weilmuenster, K. J.↗

Impact Cratering Calculations

Understanding the physical processes of impact cratering on planetary surfaces and atmospheres as well as collisions of finite-size self-gravitating objects is vitally important to planetary science. The observation has often been made that craters are the most ubiquitous landform on the solid planets and the satellites. The density of craters is used to date surfaces on planets and satellites. For large ringed basin craters (e.g. Chicxulub), the issue of identification of exactly what 'diameter' transient crater is associated with this structure is exemplified by the arguments of Sharpton et al. (1993) versus those of Hildebrand et al. (1995). The size of a transient crater, such as the K/T extinction crater at Yucatan, Mexico, which is thought to be the source of SO,-induced sulfuric acid aerosol that globally acidified surface waters as the result of massive vaporization of CASO, in the target rock, is addressed by our present project. The impact process excavates samples of planetary interiors. The degree to which this occurs (e.g. how deeply does excavation occur for a given crater diameter) has been of interest, both with regard to exposing mantle rocks in crater floors, as well as launching samples into space which become part of the terrestrial meteorite collection (e.g. lunar meteorites, SNC's from Mars). Only in the case of the Earth can we test calculations in the laboratory and field. Previous calculations predict, independent of diameter, that the depth of excavation, normalized by crater diameter, is d(sub ex)/D = 0.085 (O'Keefe and Ahrens, 1993). For Comet Shoemaker-Levy 9 (SL9) fragments impacting Jupiter, predicted excavation depths of different gas-rich layers in the atmosphere, were much larger. The trajectory and fate of highly shocked material from a large impact on the Earth, such as the K/T bolide is of interest. Melosh et al. (1990) proposed that the condensed material from the impact upon reentering the Earth's atmosphere induced. radiative heating, and producing global firestorms. The observed reentry splash of the SL-9 impact-induced plumes that reimpact Jupiter (Boslough et al., 1994) supported Melosh's K/T model. The fate of early primitive planetary atmospheres during the latter stages of planetary accretion, resulting from impactors in the 100 to 103 km diameter require modeling, e.g. Newman et al. (1997). Ahrens (1990; 1993) and Chen and Ahrens (1997) found that upon delivery of most of the impact energy to the solid planet, very large ground motions arise, which couple sufficient kinetic energy to the atmosphere to cause substantial atmospheric escape. The trade-off of this model with that of Cameron (1997) who suggests that atmospheric blow-off occurs as a result of the massive impact-induced heating of the atmosphere and Pepin (1997) who uses this heating event to model differential hydrodynamic loss of lighter atmospheric gases, requires further research.

Ahrens, Thomas J.↗

Analysis of Static Spacecraft Floating Potential at Low Earth Orbit (LEO)

Spacecraft floating potential is the charge on the external surfaces of orbiting spacecraft relative to the space. Charging is caused by unequal negative and positive currents to spacecraft surfaces. The charging process continues until the accelerated particles can be collected rapidly enough to balance the currents at which point the spacecraft has reached its equilibrium or floating potential. In low inclination. Low Earth Orbit (LEO), the collection of positive ion and negative electrons. in a particular direction. are typically not equal. The level of charging required for equilibrium to be established is influenced by the characteristics of the ambient plasma environment. by the spacecraft motion, and by the geometry of the spacecraft. Using the kinetic theory, a statistical approach for studying the interaction is developed. The approach used to study the spacecraft floating potential depends on which phenomena are being applied. and on the properties of the plasma. especially the density and temperature. The results from kinetic theory derivation are applied to determine the charging level and the electric potential distribution at an infinite flat plate perpendicular to a streaming plasma using finite-difference scheme.

Herr, Joel L.↗

Multiscale Modeling of Grain Boundaries in ZrB2: Structure, Energetics, and Thermal Resistance

A combination of ab initio, atomistic and finite element methods (FEM) were used to investigate the structures, energetics and lattice thermal conductance of grain boundaries for the ultra high temperature ceramic ZrB2. Atomic models of idealized boundaries were relaxed using density functional theory. Information about bonding across the interfaces was determined from the electron localization function. The Kapitza conductance of larger scale versions of the boundary models were computed using non-equilibrium molecular dynamics. The interfacial thermal parameters together with single crystal thermal conductivities were used as parameters in microstructural computations. FEM meshes were constructed on top of microstructural images. From these computations, the effective thermal conductivity of the polycrystalline structure was determined.

Lawson, John W.↗

Probing atmospheric water vapor profiles via multiple scattering of electromagnetic waves

A theoretical analysis on the multiple scattering of electromagnetic waves propagating in a finite inhomogeneous medium is presented and applied to the study of wave propagation in a clear atmosphere (fine weather conditions). It is shown that the analysis offers a method of synthesizing the water vapor density profile in a clear atmosphere by measuring the resultant reflections from the density profile at several different frequencies. It is also shown that the resultant reflection emerges as the consequence of multiple scattering of partial reflections from various parts of the inhomogeneous medium. The solutions of the multiple scattering approach are shown to be more accurate than those of the WKB approach, which neglects the multiple scattering effects.

Tang, C. C. H.↗

Calculations of transonic potential flow over cascades

Transonic flow through a cascade was studied by using the full potential equation and the finite volume method of Jameson and Caughey. The C-type computational grid is generated by an electrostatic analogy and simple shearing transformation. The solution algorithm includes an option of using either an artificial density or an artificial viscosity formulation of the dissipative term. Using the developed code, flows through a cascade of NACA 0012 airfoils and flows through a cascade of shockless blades were computed. It is found that the designed flow through the shockless blade is accurately predicted, the artificial density formulation shows more tolerance to the mesh irregularity, and the C-type mesh does not extend very far upstream for a small pitch-cord ratio.

Jou, W. H.↗