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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS1981

An alternating-direction implicit algorithm for unsteady potential flow

An implicit finite-difference scheme is presented for the efficient computation of unsteady potential flow about airfoils. The formulation uses density and the velocity potential as dependent variables and is cast in conservation form to assure the theoretically correct determination of shockwave location and speed. To enable boundary conditions to be imposed directly on the airfoil surface, a time-varying sheared-rectilinear coordinate transformation is employed. Calculated time-history solutions on a pulsating airfoil are compared with the results of other unsteady transonic codes, including a previous method of the authors. The present method is demonstrated to be unconditionally stable and to give accurate solutions with sharply resolved shocks.

Chipman, R.↗
NASA NTRS1982

Transonic potential flow and coordinate generation for bodies in a wind tunnel

An accurate method has been developed for computation of transonic potential flow about a 2-D lifting airfoil or an axisymmetric body in a wind tunnel. The computational mesh is nearly orthogonal everywhere and is generated by a sequence of Schwarz-Christoffel transformations and shearings to obtain an '0' grid near the body. A conservative finite-volume scheme for the full-potential equation and exact boundary conditions is used together with 'retarded' density to solve efficiently transonic flow with embedded shocks and large regions of supersonic flow. Supersonic free-stream flows are also solved with captured bow shocks and embedded subsonic regions.

Doria, M. L.↗
NASA NTRS1984

Crack tip field and fatigue crack growth in general yielding and low cycle fatigue

Fatigue life consists of crack nucleation and crack propagation periods. Fatigue crack nucleation period is shorter relative to the propagation period at higher stresses. Crack nucleation period of low cycle fatigue might even be shortened by material and fabrication defects and by environmental attack. In these cases, fatigue life is largely crack propagation period. The characteristic crack tip field was studied by the finite element method, and the crack tip field is related to the far field parameters: the deformation work density, and the product of applied stress and applied strain. The cyclic carck growth rates in specimens in general yielding as measured by Solomon are analyzed in terms of J-integral. A generalized crack behavior in terms of delta is developed. The relations between J and the far field parameters and the relation for the general cyclic crack growth behavior are used to analyze fatigue lives of specimens under general-yielding cyclic-load. Fatigue life is related to the applied stress and strain ranges, the deformation work density, crack nucleus size, fracture toughness, fatigue crack growth threshold, Young's modulus, and the cyclic yield stress and strain. The fatigue lives of two aluminum alloys correlate well with the deformation work density as depicted by the derived theory. The general relation is reduced to Coffin-Manson low cycle fatigue law in the high strain region.

Minzhong, Z.↗
NASA NTRS1984

The structures and energetics of fully nonlinear symmetric baroclinic waves

Rotating baroclinic flow for Richardson number lower than about 1 is studied by means of a finite difference Navier-Stokes model assuming no variations except in the vertical plane that completely contains the density gradient vector. The horizontally infinite channel to which attention is given further assumes periodic boundary conditions at the vertical computational boundaries and no-slip conducting horizontal boundaries. Two configurations are considered. Symmetric baroclinic waves developed in the flows in a manner consistent with linear theory, and it is noted that the structures and energetics of the fully developed waves were dependent on the Prandtl number Pr. For Pr greater than 1, the conversion from potential energy to wave kinetic energy was direct, via temperature and vertical motion correlation, while for Pr of less than 1, conversion proceeded from potential energy to average kinetic energy by means of an induced meridional flow, and then to wave kinetic energy.

Miller, T. L.↗
NASA NTRS1984

Effect of reabsorbed recombination radiation on the saturation current of direct gap p-n junctions

The application of the radiative transfer theory for semiconductors to p-n homojunctions subject to low level injection conditions is discussed. By virtue of the interaction of the radiation field with free carriers across the depletion layer, the saturation current density in Shockley's expression for the diode current is reduced at high doping levels. The reduction, due to self-induced photon generation, is noticeable for n-type material owing to the small electron effective mass in direct band-gap III-V compounds. The effect is insignificant in p-type material. At an equilibrium electron concentration of 2 x 10 to the 18th/cu cm in GaAs, a reduction of the saturation current density by 15 percent is predicted. It is concluded that realistic GaAs p-n junctions possess a finite thickness.

Von Roos, O.↗
NASA NTRS1990

Nonaxisymmetric instabilities in a slender torus - Two- and three-dimensional simulations

Nonaxisymmetric instabilities in accretion disks are investigated in terms of the slender, non-self-gravitating, ideal fluid torus model. Nonlinear simulations are carried out in two and tree dimensions, using Cartesian-grid finite differencing. The fastest growing instability is the principal mode. This mode saturates with the formation of ellipsoidal density distributions, designated as 'planets'. The simulations provide evidence for nonlinear mode-mode coupling. When several wavenumbers m are unstable, multiple 'planets' form that subsequently merge. Sources of numerical error are examined, and the effects are contrasted with those of the physical instability. The same qualitative evolution is seen in both the two- and the three-dimensional simulations, even though strict vertical hydrostatic equilibrium is no longer rigorously maintained after mode saturation.

Hawley, John F.↗
NASA NTRS1991

Studies on the dynamics of Venus ionosphere - Effects of varying the ionopause height and the onset of turbulence

The conservation equations of plasma dynamics in the upper ionosphere of Venus have been solved by using a spectral method in the horizontal and finite differencing in the vertical direction. The effect of varying the ionopause height on the computed nightside ion densities is investigated. These ion densities show a sharp decrease as the ionopause altitude is reduced to 300 km. The effect of viscous forces on the horizontal plasma flow is investigated for a wide range of values of the coefficient of viscosity. The Reynolds numbers characteristics of the flow are calculated and the conditions for the onset of turbulence discussed. It is found that the Reynolds number can be large (greater than 1000) in the subsolar region for a coefficient of viscosity of up to 1.6 x 10 to the -10th g/cm s. The influence of magnetic fields on viscosity is also discussed.

Singhal, R. P.↗
NASA NTRS1992

Numerical models for the collapse and fragmentation of centrally condensed molecular cloud cores

The gravitational collapse and fragmentation of centrally condensed molecular cloud cores are investigated using a new hydrodynamical code. The numerical scheme is second-order accurate and uses explicit finite difference methods to advance the fluid variables on a 3D Cartesian grid. Two initial power-law density profiles, rho varies as r exp -1 and rho varies as r exp -2, are considered, as well as two initial density perturbations in the azimuthal coordinate theta, rho-i - rho(1 + a cos 2theta) where a = 0.1 and 0.5. Fragmentation is found to be possible in these centrally condensed cores if the initial conditions also include differential rotation. Models which collapse with initial uniform rotation do not produce fragments. If molecular cloud cores are indeed centrally condensed, as suggested by observations of star-forming regions and by studies of ambipolar diffusion, then differential rotation may be a mechanism for producing binary protostars during gravitational collapse.

Myhill, Elizabeth A.↗
NASA NTRS1992

The prediction of noise radiation from supersonic elliptic jets

This paper describes the prediction of noise radiation from supersonic elliptic jets. The noise is associated with the large scale structures in the jet mixing layer. These structures are described as instability waves. The local characteristics of the instability waves are determined from a compressible, linear, analysis. The jet mean velocity and density are described in elliptic cylindrical coordinates. The local eigensolution for the instability waves is determined from a finite difference solution of the non-separable boundary value problem. This inner solution which is formulated in terms of the method of multiple scales is matched with the radiated field using the method of matched asymptotic expansions. The form of the far-field directivity is derived. Predictions are presented for the noise radiation by the several modes of instability in the elliptic jet. The radiated field is not axisymmetric and certain modes radiate strongly in the directions of the major and minor axes of the jet. The extension of the present work to other geometries and flow fields is discussed.

Morris, Philip J.↗
NASA NTRS1990

Amplification of a high-frequency electromagnetic wave by a relativistic plasma

The amplification of a high-frequency transverse electromagnetic wave by a relativistic plasma component, via the synchrotron maser process, is studied. The background plasma that supports the transverse wave is considered to be cold, and the energetic component whose density is much smaller than that of the background component has a loss-cone feature in the perpendicular momentum space and a finite field-aligned drift speed. The ratio of the background plasma frequency squared to the electron gyrofrequency squared is taken to be sufficiently larger than unity. Such a parameter regime is relevant to many space and astrophysical situations. A detailed study of the amplification process is carried out over a wide range of physical parameters including the loss-cone index, the ratio of the electron mass energy to the temperature of the energetic component, the field-aligned drift speed, the normalized density, and the wave propagation angle.

Yoon, Peter H.↗
NASA NTRS1987

Power flow as a complement to statistical energy analysis and finite element analysis

Present methods of analysis of the structural response and the structure-borne transmission of vibrational energy use either finite element (FE) techniques or statistical energy analysis (SEA) methods. The FE methods are a very useful tool at low frequencies where the number of resonances involved in the analysis is rather small. On the other hand SEA methods can predict with acceptable accuracy the response and energy transmission between coupled structures at relatively high frequencies where the structural modal density is high and a statistical approach is the appropriate solution. In the mid-frequency range, a relatively large number of resonances exist which make finite element method too costly. On the other hand SEA methods can only predict an average level form. In this mid-frequency range a possible alternative is to use power flow techniques, where the input and flow of vibrational energy to excited and coupled structural components can be expressed in terms of input and transfer mobilities. This power flow technique can be extended from low to high frequencies and this can be integrated with established FE models at low frequencies and SEA models at high frequencies to form a verification of the method. This method of structural analysis using power flo and mobility methods, and its integration with SEA and FE analysis is applied to the case of two thin beams joined together at right angles.

Cuschieri, J. M.↗
NASA NTRS1992

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 which contains simulated lunar regolith, a heater, nine heat shields, a heat transfer cold jacket, and 19 type-B platinum rhodium thermocouples. The simulated lunar regolith is a basalt 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, which include melt temperature, 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. A finite differencing SINDA model was developed at NASA Lewis Research Center to predict the performance of 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.↗
NASA NTRS1985

Contributions to the finite element solution of the fan noise radiation problem

The radiation of fan generated noise to the far field from a nacelle of realistic geometry is investigated using the finite element method. Several innovations have been introduced to minimize the computational requirements and create a highly efficient numerical scheme. The innovations include: (1) formulation of the problem in terms of velocity potential and density in such a way that no inlet mean flow velocity derivatives are required in the field equations, (2) the use of 'wave envelope' elements in an outer region permitting a grid much coarser than would be used for conventional finite elements, (3) the use of a mesh which deforms with an increase of forward flight speed so that mesh lines are always lines of constant phase and rays for a point source, permitting the use of wave envelope elements and simple boundary conditions for any case of forward velocity, (4) an efficient scheme for introducing the noise source via modal amplitude coefficients, and (5) the use of a frontal solution technique which for physically realistic problems drastically reduces the active storage requirements. The finite element scheme is outlined, as are the specific details of the innovations. Results are given for cases where comparable experimental data are available.

Eversman, W.↗
NASA NTRS2022

A Novel High-Performance Mission-Enabling Multi-Purpose Radioisotope Heat Source

Recent studies indicate science mission concepts targeting access to the sub-surface oceans of icy moons require ice-penetrating cryobots powered by advanced Radioisotope Power Systems (RPS). These systems would deliver waste heat for ice-melting in the range of 10 kW. Minimizing the transit time through the kilometers-thick ice shells to just a few years requires these RPS to utilize heat sources having a higher thermal energy volumetric density than the existing flight-qualified General-Purpose Heat Source (GPHS). A Compact Heat Source (CPHS)has been conceptualized in which the graphite impact shells(GIS) of the existing GPHS are rearranged in a hexagonal aeroshell containing seven GIS per module, as opposed to the standard two per module; offering a thermal energy density of 0.57 W/cm3versus 0.29 W/cm3 offered by the GPHS simply from the repackaging of Technology Readiness Level (TRL)9 subassemblies. Preliminary thermal modeling of the CPHS integrated into a notional radioisotope thermoelectric generator(RTG) structure further suggests that centerline temperatures are well within allowable limits during nominal operation. Given the need for the CPHS for a subset of missions, it is worth exploring the applicability of the CPHS for more general RTG purposes. We discuss herein how the CPHS may be implemented with either heritage or in-development thermoelectric converter technologies into a Next-Generation RTG concept. Due to a higher energy density, the legacy heat rejection fin arrangement must be modified to permit a sufficiently low cold-side temperature. Preliminary finite element analysis suggests fin-root temperatures can be kept as low as 520 K while allowing the generator to fit within the usable dimensions of currently available United States Department of Energy shipping containers. Such temperatures would certainly be compatible with the use of high temperature thermoelectric converter technologies.. A prime candidate is the heritage silicon-germanium (SiGe) unicouple, whose design could be adapted by approximately halving the leg-length, but without changes in hot and cold junction interfaces, which are features critical to the proven performance and reliability of these devices. The estimated Beginning of Life power for a SiGe-based CPHS-RTG using 12 CPHS for a thermal inventory of 10.5 kW is greater than 600 W under deep space operating conditions. Using higher performance segmented couples currently in development that are based on skutterudite,La3−xTe4and 14-1-11 Zintl thermoelectric materials in lieu of the SiGe unicouples would increase the power level to more than1 kW. The high specific power (We/kg) attribute of CPHS-RTGs found in this study could potentially enable Radioisotope Electric Propulsion (REP) mission concepts. Past NASA REP mission concept studies identified specific power needs in excess of 6to 8 We/kg. Based on a GPHS-RTG-like system configuration, we show that at fin root temperatures between 530 K and 570K (deep space environment), specific powers exceeding 10 We/kg are achievable using high performance segmented thermoelectric converters. The compact sizing and power density of the CPHS-RTG would constitute a significant step upgrade in specific power when compared to heritage GPHS-RTG (approximately5.1 We/kg) and off-the-shelf Multi-Mission RTG (approximately 2.6 We/kg).

Nesmith, Bill J.↗
NASA NTRS1986

On the electromagnetic scattering from infinite rectangular grids with finite conductivity

A variety of methods can be used in constructing solutions to the problem of mesh scattering. However, each of these methods has certain drawbacks. The present paper is concerned with a new technique which is valid for all spacings. The new method involved, called the fast Fourier transform-conjugate gradient method (FFT-CGM), represents an iterative technique which employs the conjugate gradient method to improve upon each iterate, utilizing the fast Fourier transform. The FFT-CGM method provides a new accurate model which can be extended and applied to the more difficult problems of woven mesh surfaces. The formulation of the FFT-conjugate gradient method for aperture fields and current densities for a planar periodic structure is considered along with singular operators, the formulation of the FFT-CG method for thin wires with finite conductivity, and reflection coefficients.

Christodoulou, C. G.↗
NASA NTRS1991

Electric fields and current densities under small Florida thunderstorms

Results are presented of measurements of the electric field E and Maxwell current density that were performed simultaneously under and near small Florida thunderstorms. It is shown that the amplitude of JM is of the order of 1 nA/sq cm or less in the absence of precipitation and that there are regular time variations in JM during the intervals between lightning discharges that tend to have the same shapes after different discharges in different storms. It is argued that the major causes of time variations in JM between lightning discharges are currents that flow in the finitely conducting atmosphere in response to the field changes rather than rapid time variations in the strength of cloud current sources. The displacement current densities that are computed from the E records dominate JM except when there is precipitation, when E is large and steady, or when E is unusually noisy.

Deaver, Lance E.↗
NASA NTRS1987

Probabilistic finite elements

In the Probabilistic Finite Element Method (PFEM), finite element methods have been efficiently combined with second-order perturbation techniques to provide an effective method for informing the designer of the range of response which is likely in a given problem. The designer must provide as input the statistical character of the input variables, such as yield strength, load magnitude, and Young's modulus, by specifying their mean values and their variances. The output then consists of the mean response and the variance in the response. Thus the designer is given a much broader picture of the predicted performance than with simply a single response curve. These methods are applicable to a wide class of problems, provided that the scale of randomness is not too large and the probabilistic density functions possess decaying tails. By incorporating the computational techniques we have developed in the past 3 years for efficiency, the probabilistic finite element methods are capable of handling large systems with many sources of uncertainties. Sample results for an elastic-plastic ten-bar structure and an elastic-plastic plane continuum with a circular hole subject to cyclic loadings with the yield stress on the random field are given.

Belytschko, Ted↗
NASA NTRS1972

Convection in the tanks of a rotating spacecraft

A general study of convection and mixing of a stratified fluid in a rotating container is described, with application to the special problem of fluid heating and convection in a spacecraft tank. The analysis was based on a set of approximate equations for the Navier-Stokes description of fluid convection with small density variations in a rotating system, including effects of body forces due to temperature stratification (caused by a heater) and arbitrary time-dependent rotation of the tank about a noncentral axis. An efficient numerical finite difference scheme and computational method are described for the convection of vorticity and energy in a two-dimensional tank. Special procedures were developed for analysis of the thermodynamic states resulting from the approximate flow equations derived for small density variations. Results of the numerical simulation are presented for studying the effectiveness of rotation maneuvers in mixing stratified oxygen in the tanks of an Apollo spacecraft. Significant effects of the rotation maneuvers are discussed, including the reduction of the potential for pressure collapse.

Source record↗