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 163 records · Page 9

Heat Transfer in High-Temperature Fibrous Insulation

The combined radiation/conduction heat transfer in high-porosity, high-temperature fibrous insulations was investigated experimentally and numerically. The effective thermal conductivity of fibrous insulation samples was measured over the temperature range of 300-1300 K and environmental pressure range of 1.33 x 10(exp -5)-101.32 kPa. The fibrous insulation samples tested had nominal densities of 24, 48, and 72 kilograms per cubic meter and thicknesses of 13.3, 26.6 and 39.9 millimeters. Seven samples were tested such that the applied heat flux vector was aligned with local gravity vector to eliminate natural convection as a mode of heat transfer. Two samples were tested with reverse orientation to investigate natural convection effects. It was determined that for the fibrous insulation densities and thicknesses investigated no heat transfer takes place through natural convection. A finite volume numerical model was developed to solve the governing combined radiation and conduction heat transfer equations. Various methods of modeling the gas/solid conduction interaction in fibrous insulations were investigated. The radiation heat transfer was modeled using the modified two-flux approximation assuming anisotropic scattering and gray medium. A genetic-algorithm based parameter estimation technique was utilized with this model to determine the relevant radiative properties of the fibrous insulation over the temperature range of 300-1300 K. The parameter estimation was performed by least square minimization of the difference between measured and predicted values of effective thermal conductivity at a density of 24 kilograms per cubic meters and at nominal pressures of 1.33 x 10(exp -4) and 99.98 kPa. The numerical model was validated by comparison with steady-state effective thermal conductivity measurements at other densities and pressures. The numerical model was also validated by comparison with a transient thermal test simulating reentry aerodynamic heating conditions.

Daryabeigi, Kamran↗

Existence, uniqueness and regularity of a time-periodic probability density distribution arising in a sedimentation-diffusion problem

The sedimentation and diffusion of a nonneutrally buoyant Brownian particle in vertical fluid-filled cylinder of finite length which is instantaneously inverted at regular intervals are investigated analytically. A one-dimensional convective-diffusive equation is derived to describe the temporal and spatial evolution of the probability density; a periodicity condition is formulated; the applicability of Fredholm theory is established; and the parameter-space regions are determined within which the existence and uniqueness of solutions are guaranteed. Numerical results for sample problems are presented graphically and briefly characterized.

Nitsche, Ludwig C.↗

Carbon Nanofiber-Based, High-Frequency, High-Q, Miniaturized Mechanical Resonators

High Q resonators are a critical component of stable, low-noise communication systems, radar, and precise timing applications such as atomic clocks. In electronic resonators based on Si integrated circuits, resistive losses increase as a result of the continued reduction in device dimensions, which decreases their Q values. On the other hand, due to the mechanical construct of bulk acoustic wave (BAW) and surface acoustic wave (SAW) resonators, such loss mechanisms are absent, enabling higher Q-values for both BAW and SAW resonators compared to their electronic counterparts. The other advantages of mechanical resonators are their inherently higher radiation tolerance, a factor that makes them attractive for NASA s extreme environment planetary missions, for example to the Jovian environments where the radiation doses are at hostile levels. Despite these advantages, both BAW and SAW resonators suffer from low resonant frequencies and they are also physically large, which precludes their integration into miniaturized electronic systems. Because there is a need to move the resonant frequency of oscillators to the order of gigahertz, new technologies and materials are being investigated that will make performance at those frequencies attainable. By moving to nanoscale structures, in this case vertically oriented, cantilevered carbon nanotubes (CNTs), that have larger aspect ratios (length/thickness) and extremely high elastic moduli, it is possible to overcome the two disadvantages of both bulk acoustic wave (BAW) and surface acoustic wave (SAW) resonators. Nano-electro-mechanical systems (NEMS) that utilize high aspect ratio nanomaterials exhibiting high elastic moduli (e.g., carbon-based nanomaterials) benefit from high Qs, operate at high frequency, and have small force constants that translate to high responsivity that results in improved sensitivity, lower power consumption, and im - proved tunablity. NEMS resonators have recently been demonstrated using topdown, lithographically fabricated ap - proaches to form cantilever or bridgetype structures. Top-down approaches, however, rely on complicated and expensive e-beam lithography, and often require a release mechanism. Reso - nance effects in structures synthesized using bottom-up approaches have also recently been reported based on carbon nanotubes, but such approaches have relied on a planar two-dimensional (2D) geometry. In this innovation, vertically aligned tubes synthesized using a bottom- up approach have been considered, where the vertical orientation of the tubes has the potential to increase integration density even further. The simulation of a vertically oriented, cantilevered carbon nanotube was performed using COMSOL Multi - physics, a finite element simulation package. All simulations were performed in a 2D geometry that provided consistent results and minimized computational complexity. The simulations assumed a vertically oriented, cantilevered nanotube of uniform density (1.5 g/cu cm). An elastic modulus was assumed to be 600 GPa, relative permittivity of the nanotube was assumed to be 5.0, and Poisson s ratio was assumed to be 0.2. It should be noted that the relative permittivity and Poisson s ratio for the nanotubes of interest are not known accurately. However, as in previous simulations, the relative permittivity and Poisson s ratios were treated as weak variables in the simulation, and no significant changes were recognized when these variables were varied.

Kaul, Anupama B.↗

Effect of microcracks on the thermal expansion of composite laminates

A finite element analysis was used to quantitatively predict the effect of matrix microcracks in the 90 deg plies of graphite/epoxy composites on the coefficient of thermal expansion in the 0 deg direction, alpha (y) (perpendicular to the cracks). Results were generated for (0m/90n)s, (0/+ or 45/90)s and (0/90/ + or - 45)s graphite/epoxy laminate configurations. Analytical predictions were compared with experimental results for the two quasi-isotropic laminate configurations. Both analytical and experimental results showed that microcracks reduced the effective stiffness of the 90 deg plies, thus causing the laminates, thermal response to be more like that of a (0) laminate. The change in alpha(y) was a function of lamina material properties, microcrack density, fiber orientation, and stacking sequence. A combination of classical lamination theory and finite element analysis was used to predict the effect of microcracks in both the 90 deg and 0 deg plies. Analytical results showed that the addition of microcracks in the 0 deg plies do affect alpha(y), but to a lesser extent than those in the 90 deg plies.

Bowles, D. E.↗

Investigation of Transient Gas Phase Column Density Due to Droplet Evaporation

Analytical expressions are developed for estimating column density near a rapidly evaporating droplet along general paths. The influence for instantaneous evaporation is created first as a limiting case, where the peak value occurs at the time it takes a wave of vapor to reach the closest point along the optical path traveling at its most probable thermal speed. Next the case for finite-period evaporation is evaluated for constant conditions. Compared to the instantaneous case, peak column density occurs shortly after droplet extinction at a lower intensity. A new mathematical function is discovered that solves the integrals associated with this case. Finally, ways to account for droplet motion and changes in evaporation rate with size and temperature are discussed.

transient↗

Investigation of Transient Gas Phase Column Density Due to Droplet Evaporation

Analytical expressions are developed for estimating column density near a rapidly evaporating droplet along general paths. The influence for instantaneous evaporation is created first as a limiting case, where the peak value occurs at the time it takes a wave of vapor to reach the closest point along the optical path traveling at its most probable thermal speed. Next the case for finite-period evaporation is evaluated for constant conditions. Compared to the instantaneous case, peak column density occurs shortly after droplet extinction at a lower intensity. A new mathematical function is discovered that solves the integrals associated with this case. Finally, ways to account for droplet motion and changes in evaporation rate with size and temperature are discussed.

Column Density↗

A Graphical Method for Estimating Ion-Rocket Performance

Equations relating the critical temperature and ion current density for surface ionization of cesium on tungsten are derived for the cases of zero and finite electric fields at the ion-emitting surface. These equations are used to obtain a series of graphs that can be used to solve many problems relating to ion-rocket theoretical performance. The effect of operation at less than space-charge-limited current density and the effect of nonuniform propellant flux onto the ion-emitting surface are also treated.

Reynolds, Thaine W.↗

Manual for automatic generation of finite element models of spiral bevel gears in mesh

The goal of this research is to develop computer programs that generate finite element models suitable for doing 3D contact analysis of faced milled spiral bevel gears in mesh. A pinion tooth and a gear tooth are created and put in mesh. There are two programs: Points.f and Pat.f to perform the analysis. Points.f is based on the equation of meshing for spiral bevel gears. It uses machine tool settings to solve for an N x M mesh of points on the four surfaces, pinion concave and convex, and gear concave and convex. Points.f creates the file POINTS.OUT, an ASCI file containing N x M points for each surface. (N is the number of node points along the length of the tooth, and M is nodes along the height.) Pat.f reads POINTS.OUT and creates the file tl.out. Tl.out is a series of PATRAN input commands. In addition to the mesh density on the tooth face, additional user specified variables are the number of finite elements through the thickness, and the number of finite elements along the tooth full fillet. A full fillet is assumed to exist for both the pinion and gear.

Bibel, G. D.↗

Acoustic Radiation Pressure

The theoretical foundation of acoustic radiation pressure in plane wave beams is reexamined. It is shown from finite deformation theory and the Boltzmann-Ehrenfest Adiabatic Principle that the Brillouin stress tensor (BST) is the radiation stress in Lagrangian coordinates (not Eulerian coordinates) and that the terms in the BST are not the momentum flux density and mean excess Eulerian stress but are simply contributions to the variation in the wave oscillation period resulting from changes in path length and true wave velocity, respectively, from virtual variations in the strain. It is shown that the radiation stress in Eulerian coordinates is the mean Cauchy stress (not the momentum flux density, as commonly assumed) and that Langevin's second relation does not yield an assessment of the mean Eulerian pressure, since the enthalpy used in the traditional derivations is a function of the thermodynamic tensions - not the Eulerian pressure. It is shown that the transformation between Lagrangian and Eulerian quantities cannot be obtained from the commonly-used expansion of one of the quantities in terms of the particle displacement, since the expansion provides only the difference between the value of the quantity at two different points in Cartesian space separated by the displacement. The proper transformation is obtained only by employing the transformation coefficients of finite deformation theory, which are defined in terms of the displacement gradients. Finite deformation theory leads to the result that for laterally unconfined, plane waves the Lagrangian and Eulerian radiation pressures are equal with the value (1/4)(2K) along the direction of wave propagation, where (K) is the mean kinetic energy density, and zero in directions normal to the propagation direction. This is contrary to the Langevin result that the Lagrangian radiation pressure in the propagation direction is equal to (2K) and the BST result that the Eulerian radiation pressure in that direction is the momentum flux density.

Finite deformation theory↗

Investigation of Transient Gas Phase Column Density Due to Droplet Evaporation

A NASA robotic refueling mission experiment on board the International Space Station (ISS) was designed to repeatedly transfer simulated cryogenic propellant between two dewars. After each metered transfer, the cryogen was to be vented overboard via sublimation or evaporation. ISS payload providers must conduct analyses to demonstrate that any planned gaseous venting generate no more than a certain level of material that may interfere with optical measurements from other payloads that may be located nearby. This requirement is expressed in terms of a maximum column number density (CND). For the experiment under consideration, fluid droplets likely would accompany any such vapor. Earlier related studies led to development of analytical expressions for column density along general paths for a number of practical cases under steady conditions, including droplet evaporation. These expressions are not appropriate for rapid evaporation of small droplets however, since the droplets cannot sustain mass flow rates consistent with steady conditions. In this work, analytical expressions are developed for estimating column density near a rapidly evaporating droplet along general paths. The influence for instantaneous evaporation is created first as a limiting case, where it is found the peak value occurs at the time it takes a wave of vapor to reach the closest point along the optical path traveling at its most probable thermal speed. Next, the case for finite-period evaporation is evaluated for constant conditions. Compared to the instantaneous case, peak column density occurs shortly after droplet extinction, but at a lower intensity. A new mathematical function is discovered that solves the integrals associated with this case. Finally, ways to account for droplet motion and changes in evaporation rate with size and temperature are discussed in this framework.

Woronowicz, M. S.↗

A design of an all digital FM discriminator.

This paper describes the design of an all digital FM discriminator. The system can be easily constructed using ordinary commercially available, digital IC's. In addition, the output SNR above threshold is computed. It is shown that the output noise consists of two terms: the first being the expected output gaussian noise having a parabolic power spectral density and the second being a quantization noise produced by the finite-level detection process used in the receiver.

Garodnick, J.↗

Theory of type 3b solar radio bursts

During the initial space-time evolution of an electron beam injected into the corona, the strong beam-plasma interaction occurs at the head of the beam, leading to the amplification of a quasi-monochromatic large-amplitude plasma wave that stabilizes by trapping the beam particles. Oscillation of the trapped particles in the wave troughs amplifies sideband electrostatic waves. The sidebands and the main wave subsequently decay to observable transverse electromagnetic waves through the parametric decay instability. This process gives rise to the elementary striation bursts. Owing to velocity dispersion in the beam and the density gradient of the corona, the entire process may repeat at a finite number of discrete plasma levels, producing chains of elementary bursts. All the properties of the type IIIb bursts are accounted for in the context of the theory.

Smith, R. A.↗

Theory of type IIIb solar radio bursts

During the initial space-time evolution of an electron beam injected into the corona, the strong beam-plasma interaction occurs at the head of the beam, leading to the amplification of a quasi-monochromatic large-amplitude plasma wave that stabilizes by trapping the beam particles. Oscillation of the trapped particles in the wave troughs amplifies sideband electrostatic waves. The sidebands and the main wave subsequently decay to observable transverse electromagnetic waves through the parametric decay instability. This process gives rise to the elementary striation bursts. Owing to velocity dispersion in the beam and the density gradient of the corona, the entire process may repeat at a finite number of discrete plasma levels, producing chains of elementary bursts. All the properties of the type IIIb bursts are accounted for in the context of the theory.

Smith, R. A.↗

Algorithm for Unsteady Potential Flow About Airfoils

Implicit finite-difference scheme efficiently computes unsteady potential flow about airfoils. Formulation uses density and velocity potential as dependent variables. Conservation form is retained to assure that shock wave location and speed are computed correctly. Scheme fills need for method to calculate efficiently unsteady potential flow about airfoils and to predict flutter and other unsteady aeroelastic phenomena in transonic flow regimes.

Chipman, R.↗

Systematic generation of nonlinear discretized dynamic equilibrium equations of spinning cantilevers

General nonlinear discretized governing equations of motion of spinning elastic solids and structures are adjusted for the case of a spinning cantilever with initial geometric imperfections. Consideration is given to second degree nonlinearities in the strain-displacement and velocity-displacement relationships. Parameters of the discretization are developed to include the type and number of the coordinate functions used in the admissible trial solution in order to unify the discretization approaches associated with stationarity principles. The coordinate functions comprise both sets of continuous and piecewise continuous functions employed in the Rayleigh-Ritz and the finite element methods, respectively. Coefficient matrices are provided which contain the energy density expressions and which are adaptable to computer programming.

El-Essawi, M.↗

Interpretation of the geoid

The superposition of the first satellite geoid determined by Iszak upon Ootilla's geoid was based on surface gravity determinations. Good agreement was observed except over the Pacific area of the globe. The poor agreement over the Pacific was interpreted as the result of inadequate observations there. Many geoids were determined from satellite observations, including Doppler measurements. It is found that the geoid is the result of density differences in the mantle maintained since the primeval Earth by its finite strength. Various models based on this assumption are developed.

Runcorn, S. K.↗

Involute composite design evaluation using global design sensitivity derivatives

An optimization capability for involute structures has been developed. Its key feature is the use of global material geometry variables which are so chosen that all combinations of design variables within a set of lower and upper bounds correspond to manufacturable designs. A further advantage of global variables is that their number does not increase with increasing mesh density. The accuracy of the sensitivity derivatives has been verified both through finite difference tests and through the successful use of the derivatives by an optimizer. The state of the art in composite design today is still marked by point design algorithms linked together using ad hoc methods not directly related to a manufacturing procedure. The global design sensitivity approach presented here for involutes can be applied to filament wound shells and other composite constructions using material form features peculiar to each construction. The present involute optimization technology is being applied to the Space Shuttle SRM nozzle boot ring redesigns by PDA Engineering.

Hart, J. K.↗