Search NASA⌕ Search

SEARCH · Search NASA

Results for “Approximations and Expansions”

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 361 records · Page 20

An exact derivation of contact resistance to planar devices

A mixed boundary-value problem is formulated for an imperfect rectangular contact to a semiconducting sheet (layer). An exact model for derivation of contact resistance is developed and its solution by conformal mapping and eigenfunction expansion is presented. The expression for contact resistance is similar to that for a perfect (lossless) contact but includes an additional term containing the lowest-order coefficient in the infinite series expansion of the complex potential function. The calculated contact resistances are compared with those obtained from three approximate models: lossless contact, transmission line, and extended transmission line models. The extended transmission line model appears to be a very satisfactory approximation provided the ratio of contact length to sheet thickness is no less than 0.5.

Schuldt, S. B.↗

New approach to the solution of large, full matrix equations

A new approach to the solution of matrix equations resulting from integral equations is presented and applied to the solution of two-dimensional Neumann problems describing the inviscid, incompressible flow past an airfoil. The problem is reformulated in terms of a preselected set of mode functions giving an equivalent matrix equation to be solved for the mode-function expansion coefficients. Because of the inherent smoothness of the original problem, the coefficient problem can be solved approximately without significantly affecting the accuracy of the final solution. Very promising two-dimensional results are obtained and the extension of the method to three-dimensional problems is investigated. On the basis of these results it is shown that the computing time for the matrix solution for a large three-dimensional panel method calculation could be reduced by an order of magnitude compared with that required for a direct solution.

Clark, R. W.↗

An approximate atmospheric guidance law for aeroassisted plane change maneuvers

An approximate optimal guidance law for the aeroassisted plane change problem is presented which is based upon an expansion of the Hamilton-Jacobi-Bellman equation with respect to the small parameter of Breakwell et al. (1985). The present law maximizes the final velocity of the reentry vehicle while meeting terminal constraints on altitude, flight path angle, and heading angle. The integrable zeroth-order solution found when the small parameter is set to zero corresponds to a solution of the problem where the aerodynamic forces dominate the inertial forces. Higher order solutions in the expansion are obtained from the solution of linear partial differential equations requiring only quadrature integration.

Speyer, Jason L.↗

A linear shock cell model for non-circular jets using conformal mapping with a pseudo-spectral hybrid scheme

The shock structure in non-circular supersonic jets is predicted using a linear model. This model includes the effects of the finite thickness of the mixing layer and the turbulence in the jet shear layer. A numerical solution is obtained using a conformal mapping grid generation scheme with a hybrid pseudo-spectral discretization method. The uniform pressure perturbation at the jet exit is approximated by a Fourier-Mathieu series. The pressure at downstream locations is obtained from an eigenfunction expansion that is matched to the pressure perturbation at the jet exit. Results are presented for a circular jet and for an elliptic jet of aspect ratio 2.0. Comparisons are made with experimental data.

Bhat, Thonse R. S.↗

Segmented multigrid domain decomposition procedure for incompressible viscous flows

In this paper, the Navier-Stokes (NS) equations are approximated with a reduced NS system, that represents the lowest-order terms in an asymptotic Re expansion. This system allows for simplified boundary conditions, more generality in the location of the outflow boundary, and insures mass conservation in all subdomain grid interfaces, and also at the outflow boundary. The higher-order NS diffusion terms are included through a deferred corrector, in chosen subdomains, when required.

Srinivasan, Kumar↗

Quantifying shapes of volcanoes on Venus

A large population of discrete volcanic edifices on Venus has been identified and cataloged by means of Magellan SAR images, and an extensive database describing thousands of such features is in final preparation. Those volcanoes categorized as Intermediate to Large in scale, while relatively small in number (approx. 400), nonetheless constitute a significant volumetric component (approx. 13 x 10(exp 6) cu km) of the total apparent crustal volume of Venus. For this reason, we have focused attention on the morphometry of a representative suite of the larger edifices on Venus and, in particular, on ways of constraining the eruptive histories of these possibly geologically youthful landforms. Our approach has been to determine a series of reproducible morphometric parameters for as many of the discrete volcanoes on Venus that have an obvious expression within the global altimetry data acquired by Magellan. In addition, we have attempted to objectively and systematically define the mathematical essence of the shapes of these larger volcanoes using a polynomial cross-section approximation involving only parameters easily measured from digital topography, as well as with simple surface cylindrical harmonic expansions. The goal is to reduce the topological complexities of the larger edifices to a few simple parameters which can then be related to similar expressions for well-studied terrestrial and martian features.

Garvin, J. B.↗

Radiative signatures from impact of Comet Shoemaker-Levy-9 on Jupiter

The visible optical power emitted from the expansion plumes from 0.4 and 2 km diameter fragments of Shoemaker-Levy are expected to be, approximately 25% and comparable to, the visible solar flux reflected from Jupiter, respectively, for several minutes, and could be easily observed by sensors on the Galileo spacecraft. Earth-based observers can detect these plumes as these expand over the SW limb of Jupiter and come into earth view some minutes after impact!

Ahrens, Thomas J.↗

Accuracy of the Generalized Self-Consistent Method in Modelling the Elastic Behaviour of Periodic Composites

Local stress and strain fields in the unit cell of an infinite, two-dimensional, periodic fibrous lattice have been determined by an integral equation approach. The effect of the fibres is assimilated to an infinite two-dimensional array of fictitious body forces in the matrix constituent phase of the unit cell. By subtracting a volume averaged strain polarization term from the integral equation we effectively embed a finite number of unit cells in a homogenized medium in which the overall stress and strain correspond to the volume averaged stress and strain of the constrained unit cell. This paper demonstrates that the zeroth term in the governing integral equation expansion, which embeds one unit cell in the homogenized medium, corresponds to the generalized self-consistent approximation. By comparing the zeroth term approximation with higher order approximations to the integral equation summation, both the accuracy of the generalized self-consistent composite model and the rate of convergence of the integral summation can be assessed. Two example composites are studied. For a tungsten/copper elastic fibrous composite the generalized self-consistent model is shown to provide accurate, effective, elastic moduli and local field representations. The local elastic transverse stress field within the representative volume element of the generalized self-consistent method is shown to be in error by much larger amounts for a composite with periodically distributed voids, but homogenization leads to a cancelling of errors, and the effective transverse Young's modulus of the voided composite is shown to be in error by only 23% at a void volume fraction of 75%.

Walker, Kevin P.↗

Theoretical Performance of Hydrogen-Oxygen Rocket Thrust Chambers

Data are presented for liquid-hydrogen-liquid-oxygen thrust chambers at chamber pressures from 15 to 1200 pounds per square inch absolute, area ratios to approximately 300, and percent fuel from about 8 to 34 for both equilibrium and frozen composition during expansion. Specific impulse in vacuum, specific impulse, combustion-chamber temperature, nozzle-exit temperature, characteristic velocity, and the ratio of chamber-to-nozzle-exit pressure are included. The data are presented in convenient graphical forms to allow quick calculation of theoretical nozzle performance with over- or underexpansion, flow separation, and introduction of the propellants at various initial conditions or heat loss from the combustion chamber.

Sievers, Gilbert K.↗

The First in situ Observation of Kelvin-Helmholtz Waves at High-Latitude Magnetopause during Strongly Dawnward Interplanetary Magnetic Field Conditions

We report the first in situ observation of high-latitude magnetopause (near the northern duskward cusp) Kelvin-Helmholtz waves (KHW) by Cluster on January 12, 2003, under strongly dawnward interplanetary magnetic field (IMF) conditions. The fluctuations unstable to Kelvin-Helmholtz instability (KHI) are found to propagate mostly tailward, i.e., along the direction almost 90 deg. to both the magnetosheath and geomagnetic fields, which lowers the threshold of the KHI. The magnetic configuration across the boundary layer near the northern duskward cusp region during dawnward IMF is similar to that in the low-latitude boundary layer under northward IMF, in that (1) both magnetosheath and magnetospheric fields across the local boundary layer constitute the lowest magnetic shear and (2) the tailward propagation of the KHW is perpendicular to both fields. Approximately 3-hour-long periods of the KHW during dawnward IMF are followed by the rapid expansion of the dayside magnetosphere associated with the passage of an IMF discontinuity that characterizes an abrupt change in IMF cone angle, Phi = acos (B(sub x) / absolute value of Beta), from approx. 90 to approx. 10. Cluster, which was on its outbound trajectory, continued observing the boundary waves at the northern evening-side magnetopause during sunward IMF conditions following the passage of the IMF discontinuity. By comparing the signatures of boundary fluctuations before and after the IMF discontinuity, we report that the frequencies of the most unstable KH modes increased after the discontinuity passed. This result demonstrates that differences in IMF orientations (especially in f) are associated with the properties of KHW at the high-latitude magnetopause due to variations in thickness of the boundary layer, and/or width of the KH-unstable band on the surface of the dayside magnetopause.

solar wind↗

Electromagnetic and scalar diffraction by a right-angled wedge with a uniform surface impedance

The diffraction of an electromagnetic wave by a perfectly-conducting right-angled wedge with one surface covered by a dielectric slab or absorber is considered. The effect of the coated surface is approximated by a uniform surface impedance. The solution of the normally incident electromagnetic problem is facilitated by introducing two scalar fields which satisfy a mixed boundary condition on one surface of the wedge and a Neumann of Dirichlet boundary condition on the other. A functional transformation is employed to simplify the boundary conditions so that eigenfunction expansions can be obtained for the resulting Green's functions. The eigenfunction expansions are transformed into the integral representations which then are evaluated asymptotically by the modified Pauli-Clemmow method of steepest descent. A far zone approximation is made to obtain the scattered field from which the diffraction coefficient is found for scalar plane, cylindrical or sperical wave incident on the edge. With the introduction of a ray-fixed coordinate system, the dyadic diffraction coefficient for plane or cylindrical EM waves normally indicent on the edge is reduced to the sum of two dyads which can be written alternatively as a 2 X 2 diagonal matrix.

Hwang, Y. M.↗

Plasma waves associated with the first AMPTE magnetotail barium release

Plasma waves observed during the March 21, 1985, AMPTE magnetotail barium release are described. Electron plasma oscillations provided local measurements of the plasma density during both the expansion and decay phases. Immediately after the explosion, the electron density reached a peak of about 400,000/cu cm, and then started decreasing approximately as t to the -2.4 as the cloud expanded. About 6 minutes after the explosion, the electron density suddenly began to increase, reached a secondary peak of about 240/cu cm, and then slowly decayed down to the preevent level over a period of about 15 minutes. The density increase is believed to be caused by the collapse of the ion cloud into the diamagnetic cavity created by the initial expansion. The plasma wave intensities observed during the entire event were quite low. In the diamagnetic cavity, electrostatic emissions were observed near the barium ion plasma frequency, and in another band at lower frequencies. A broadband burst of electrostatic noise was also observed at the boundary of the diamagnetic cavity. Except for electron plasma oscillations, no significant wave activity was observed outside of the diamagnetic cavity.

Gurnett, D. A.↗

Computational Analyses in Support of Sub-scale Diffuser Testing for the A-3 Facility: Steady Predictions - Part 1

Simulation technology can play an important role in rocket engine test facility design and development by assessing risks, providing analysis of dynamic pressure and thermal loads, identifying failure modes and predicting anomalous behavior of critical systems. Advanced numerical tools assume greater significance in supporting testing and design of high altitude testing facilities and plume induced testing environments of high thrust engines because of the greater inter-dependence and synergy in the functioning of the different sub-systems. This is especially true for facilities such as the proposed A-3 facility at NASA SSC because of a challenging operating envelope linked to variable throttle conditions at relatively low chamber pressures. Facility designs in this case will require a complex network of diffuser ducts, steam ejector trains, fast operating valves, cooling water systems and flow diverters that need to be characterized for steady state performance. In this paper, we will demonstrate with the use of CFD analyses s advanced capability to evaluate supersonic diffuser and steam ejector performance in a sub-scale A-3 facility at NASA Stennis Space Center (SSC) where extensive testing was performed. Furthermore, the focus in this paper relates to modeling of critical sub-systems and components used in facilities such as the A-3 facility. The work here will address deficiencies in empirical models and current CFD analyses that are used for design of supersonic diffusers/turning vanes/ejectors as well as analyses for confined plumes and venting processes. The primary areas that will be addressed are: (1) supersonic diffuser performance including analyses of thermal loads (2) accurate shock capturing in the diffuser duct; (3) effect of turning duct on the performance of the facility (4) prediction of mass flow rates and performance classification for steam ejectors (5) comparisons with test data from sub-scale diffuser testing and assessment of confidence levels in CFD based flowpath modeling of the facility. The analyses tools used here expand on the multi-element unstructured CFD which has been tailored and validated for impingement dynamics of dry plumes, complex valve/feed systems, and high pressure propellant delivery systems used in engine and component test stands at NASA SSC. The analyses performed in the evaluation of the sub-scale diffuser facility explored several important factors that influence modeling and understanding of facility operation such as (a) importance of modeling the facility with Real Gas approximation, (b) approximating the cluster of steam ejector nozzles as a single annular nozzle, (c) existence of mixed subsonic/supersonic flow downstream of the turning duct, and (d) inadequacy of two-equation turbulence models in predicting the correct pressurization in the turning duct and expansion of the second stage steam ejectors. The procedure used for modeling the facility was as follows: (i) The engine, test cell and first stage ejectors were simulated with an axisymmetric approximation (ii) the turning duct, second stage ejectors and the piping downstream of the second stage ejectors were analyzed with a three-dimensional simulation utilizing a half-plane symmetry approximation. The solution i.e. primitive variables such as pressure, velocity components, temperature and turbulence quantities were passed from the first computational domain and specified as a supersonic boundary condition for the second simulation. (iii) The third domain comprised of the exit diffuser and the region in the vicinity of the facility (primary included to get the correct shock structure at the exit of the facility and entrainment characteristics). The first set of simulations comprising the engine, test cell and first stage ejectors was carried out both as a turbulent real gas calculation as well as a turbulent perfect gas calculation. A comparison for the two cases (Real Turbulent and Perfect gas turbulent) of the Ma Number distribution and temperature distributions are shown in Figures 1 and 2 respectively. The Mach Number distribution shows small yet distinct differences between the two cases such as locations of shocks/shock reflections and a slightly different impingement point on the wall of the diffuser from the expansion at the exit of the nozzle. Similarly the temperature distribution indicates different flow recirculation patterns in the test cell. Both cases capture all the essential flow phenomena such as the shock-boundary layer interaction, plume expansion, expansion of the first stage ejectors, mixing between the engine plume and the first stage ejector flow and pressurization due to the first stage ejectors. The final paper will discuss thermal loads on the walls of the diffuser and cooling mechanisms investigated.

Allgood, Daniel C.↗

Melting of iron by significant structure theory

From Eyring's method of significant structures, a partition function is derived for liquid iron. The solid at high temperature is described in the Einstein approximation. Magnetic and electronic contributions to the thermodynamic properties of both the liquid and solid phases are considered. The model is compatible with properties (thermal expansion, compressibility, heat capacity, entropy of melting, and volume change on melting) at one atm. The melting temperature at high pressure is found by satisfying the requirement that the Gibbs free energies of the liquid and solid phases are equal at the melting temperature. Under conditions at the earth's core-mantle boundary, the melting temperature of iron is greater than approximately 5000 K, and under inner-outer core conditions the melting temperature is greater than approximately 7000 K. These estimates are consistent with the Lindemann melting law, but not with the Kraut-Kennedy melting law.

Leppaluoto, D. A.↗

Computational Analyses in Support of Sub-scale Diffuser Testing for the A-3 Facility: Steady Predictions - Part 1

Simulation technology can play an important role in rocket engine test facility design and development by assessing risks, providing analysis of dynamic pressure and thermal loads, identifying failure modes and predicting anomalous behavior of critical systems. Advanced numerical tools assume greater significance in supporting testing and design of high altitude testing facilities and plume induced testing environments of high thrust engines because of the greater inter-dependence and synergy in the functioning of the different sub-systems. This is especially true for facilities such as the proposed A-3 facility at NASA SSC because of a challenging operating envelope linked to variable throttle conditions at relatively low chamber pressures. Facility designs in this case will require a complex network of diffuser ducts, steam ejector trains, fast operating valves, cooling water systems and flow diverters that need to be characterized for steady state performance. In this paper, we will demonstrate with the use of CFD analyses s advanced capability to evaluate supersonic diffuser and steam ejector performance in a sub-scale A-3 facility at NASA Stennis Space Center (SSC) where extensive testing was performed. Furthermore, the focus in this paper relates to modeling of critical sub-systems and components used in facilities such as the A-3 facility. The work here will address deficiencies in empirical models and current CFD analyses that are used for design of supersonic diffusers/turning vanes/ejectors as well as analyses for confined plumes and venting processes. The primary areas that will be addressed are: (1) supersonic diffuser performance including analyses of thermal loads (2) accurate shock capturing in the diffuser duct; (3) effect of turning duct on the performance of the facility (4) prediction of mass flow rates and performance classification for steam ejectors (5) comparisons with test data from sub-scale diffuser testing and assessment of confidence levels in CFD based flowpath modeling of the facility. The analyses tools used here expand on the multi-element unstructured CFD which has been tailored and validated for impingement dynamics of dry plumes, complex valve/feed systems, and high pressure propellant delivery systems used in engine and component test stands at NASA SSC. The analyses performed in the evaluation of the sub-scale diffuser facility explored several important factors that influence modeling and understanding of facility operation such as (a) importance of modeling the facility with Real Gas approximation, (b) approximating the cluster of steam ejector nozzles as a single annular nozzle, (c) existence of mixed subsonic/supersonic flow downstream of the turning duct, and (d) inadequacy of two-equation turbulence models in predicting the correct pressurization in the turning duct and expansion of the second stage steam ejectors. The procedure used for modeling the facility was as follows: (i) The engine, test cell and first stage ejectors were simulated with an axisymmetric approximation (ii) the turning duct, second stage ejectors and the piping downstream of the second stage ejectors were analyzed with a three-dimensional simulation utilizing a half-plane symmetry approximation. The solution i.e. primitive variables such as pressure, velocity components, temperature and turbulence quantities were passed from the first computational domain and specified as a supersonic boundary condition for the second simulation. (iii) The third domain comprised of the exit diffuser and the region in the vicinity of the facility (primary included to get the correct shock structure at the exit of the facility and entrainment characteristics). The first set of simulations comprising the engine, test cell and first stage ejectors was carried out both as a turbulent real gas calculation as well as a turbulent perfect gas calculation. A comparison for the two cases (Real Turbulent and Perfect gas turbulent) of the Ma Number distribution and temperature distributions are shown in Figures 1 and 2 respectively.

Allgood, Daniel C.↗

Growth-phase thinning of the near-Earth current sheet during the CDAW 6 substorm

The thinning of the near-Earth current sheet during the growth phase of the Coordinated Data Analysis Workshop (CDAW) 6 magnetospheric substorm is studied. The expansion onset of the substorm occurred at 1054 UT, March 22, 1979. During the growth phase, two spacecraft, International Sun Earth Explorer (ISEE) 1 and ISEE 2, were within the current sheet approximately 13 R(sub E) from the Earth and obtained simultaneous high-resolution magnetic data at two points in the current sheet. Plasma data were also provided by the ISEE spacecraft and solar wind data by IMP 8. To facilitate the analysis, the GSM magnetic field data are transformed to a 'neutral sheet coordinate system' in which the new x axis is parallel to the average magnetic field above and below the neutral sheet and the new y axis lies in the GSM equatorial plane. A model based on the assumption that the current sheet is a time-invariant structure fails to predict neutral sheet crossing times. Consequently, the Harris sheet model, which allows one to remove the restriction of time invariancy, is used instead. It is found that during the growth phase, a model parameter corresponding to the thickness of the current sheet decreased exponentially from about 5 R(sub E) to 1 R(sub E) with a time constant of about 14 min. In addition, the ISEE 1 and ISEE 2 neutral sheet crossings after expansion onset indicate that the neutral sheet was moving upward at 7 km/s relative to the spacecraft. Since both crossings occurred in approximately 80 s, the current sheet thickness is estimated to be about 500 km. These results demonstrate that the near-Earth current sheet undergoes dramatic thinning during the substorm growth phase and expansion onset.

Sanny, Jeff↗

Field dependence of gaseous-ion mobility - Theoretical tests of approximate formulas.

The approximate formulas considered include relations based on the Wannier free-flight theory, the Kihara medium-field expansion, and the Frost-Patterson interpolation formulas. A few accurate theoretical results are available for testing the foregoing formulas. Cases concerning high fields, intermediate fields, and resonant charge transfer are examined. It is found that of the formulas tested, the one based on the Wannier free-flight theory is the most flexible, since it can be used for all fields and all ion-neutral force laws and mass ratios.

Hahn, H.↗