Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Smooth approximation”

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

Tidal interactions between spherical galaxies

The impulse approximation and full N-body experiments are used to investigate the changes in mass and binding energy which result from gravitational encounters between spherical galaxies. Models with King (1966) and de Vaucouleurs (1948) surface density profiles, and radial, tangential, and isotropic velocity distributions. Results based on the impulse approximation depend on two collision parameters (the impact parameter and the ratio of the perturbing mass to the relative velocity at infinity); they agree well with N-body results. The asymptotic limit, valid for distant encounters, is found to approximate the correct results only for impact parameters larger than 10 effective radii. It is therefore of little use in most situations of interest. No simple 'linear' or 'quadratic' dependence of energy and mass loss on collision strength is found, but rather, a smooth intermediate behavior. Cross sections depend very little on the internal velocity structure of the galaxy, but quite strongly on the surface density profile in its outer parts. The reasons why these results disagree with the conclusions of previous workers are discussed.

Aguilar, L. A.↗

Numerical Analysis of Temperature Gradients and Interface Shape During Directional Solidification of Al and Al-Cu Alloy Under Microgravity Conditions

Numerical modeling was undertaken to analyze the influence of radial thermal gradient on solid/liquid interface shape and convection patterns during solidification of pure Al and Al-4 wt. % Cu alloy. Steady state calculations were performed for different gravity levels and orientations. Furthermore, transient modeling was undertaken to investigate effect of the solidification velocity. The furnace configuration used in this analysis is the proposed International Space Station Furnace. Results from a thermal model of the furnace core were used as initial boundary conditions for solidification modeling. The Solidification model was adopted from previous work and was based on the finite element code FIDAP. Thermocouple data and quenched interface shape from a pure Al sample flown on the Life Sciences and Microgravity Spacelab (LMS) mission, July 1996 was used for model validation. Good agreement was obtained between the predicted interface shape and that measured from the quenched LMS sample. It was found that the imposed temperature boundary condition must be sufficiently smooth for the heat fluxes in the model to be self-consistent. The model predicted that alloy sample Al-Cu is more sensitive to variations in the gravity level compared to pure Al. For Al-Cu alloy, solute diffusivity is approximately 10(exp 4) times smaller than thermal diffusivity of pure Al. Hence, in a microgravity environment the weak convection has no measurable effect on the heat fluxes, but is still strong enough to affect concentration distribution. Since the alloy melting temperature is determined by concentration the interface shape depends on the level of convection. This was found not to be the case for pure metal.

Bune, Andris V.↗

Algebraic Multigrid by Smoothed Aggregation for Second and Fourth Order Elliptic Problems

Multigrid methods are very efficient iterative solvers for system of algebraic equations arising from finite element and finite difference discretization of elliptic boundary value problems. The main principle of multigrid methods is to complement the local exchange of information in point-wise iterative methods by a global one utilizing several related systems, called coarse levels, with a smaller number of variables. The coarse levels are often obtained as a hierarchy of discretizations with different characteristic meshsizes, but this requires that the discretization is controlled by the iterative method. To solve linear systems produced by existing finite element software, one needs to create an artificial hierarchy of coarse problems. The principal issue is then to obtain computational complexity and approximation properties similar to those for nested meshes, using only information in the matrix of the system and as little extra information as possible. Such algebraic multigrid method that uses the system matrix only was developed by Ruge. The prolongations were based on the matrix of the system by partial solution from given values at selected coarse points. The coarse grid points were selected so that each point would be interpolated to via so-called strong connections. Our approach is based on smoothed aggregation introduced recently by Vanek. First the set of nodes is decomposed into small mutually disjoint subsets. A tentative piecewise constant interpolation (in the discrete sense) is then defined on those subsets as piecewise constant for second order problems, and piecewise linear for fourth order problems. The prolongation operator is then obtained by smoothing the output of the tentative prolongation and coarse level operators are defined variationally.

Vanek, Petr↗

Production of polydisperse sprays.

A device which produces a spray with a controllable drop size distribution is described. This system is especially well suited for generating nonsmooth (e.g., bimodal) distributions, although smooth distributions can also be produced. Droplets of uniform size result from the rapid growth of an oscillatory disturbance in a free liquid jet. Simultaneous generation of droplets of various diameters by this method produces the polydisperse spray. The range of drop diameters is approximately 300-3000 microns.

Pierce, T. H.↗

Two-dimensional magnetohydrodynamic turbulence - Cylindrical, non-dissipative model

Incompressible magnetohydrodynamic turbulence is treated in the presence of cylindrical boundaries which are perfectly conducting and rigidly smooth. The model treated is non-dissipative and two-dimensional, the variation of all quantities in the axial direction being ignored. Equilibrium Gibbs ensemble predictions are explored assuming the constraint of constant axial current (appropriate to tokamak operation). No small-amplitude approximations are made. The expectation value of the turbulent kinetic energy is found to approach zero for the state of maximum mean-square vector potential to energy ratio. These are the only states for which large velocity fluctuations are not expected.

Montgomery, D.↗

Determination of wind from NIMBUS 6 satellite sounding data

Objective methods of computing upper level and surface wind fields from NIMBUS 6 satellite sounding data are developed. These methods are evaluated by comparing satellite derived and rawinsonde wind fields on gridded constant pressure charts in four geographical regions. Satellite-derived and hourly observed surface wind fields are compared. Results indicate that the best satellite-derived wind on constant pressure charts is a geostrophic wind derived from highly smoothed fields of geopotential height. Satellite-derived winds computed in this manner and rawinsonde winds show similar circulation patterns except in areas of small height gradients. Magnitudes of the standard deviation of the differences between satellite derived and rawinsonde wind speeds range from approximately 3 to 12 m/sec on constant pressure charts and peak at the jet stream level. Fields of satellite-derived surface wind computed with the logarithmic wind law agree well with fields of observed surface wind in most regions. Magnitudes of the standard deviation of the differences in surface wind speed range from approximately 2 to 4 m/sec, and satellite derived surface winds are able to depict flow across a cold front and around a low pressure center.

Carle, W. E.↗

Theory for computing the field scattered from a smooth inflected surface

A theory is described for computing the reflected or scattered field from a smooth body with inflection points. These inflections occur in certain directions at each surface point for which the total (Gaussian) curvature is zero or negative. For surface illumination in one of these critical directions, the usual reflection formulas obtained by the high-frequency approximation are inapplicable, and a shadow zone exists in the reflected field. Scattering into the shadow zone is treated, as well as specular reflection. This theory should have a variety of applications such as for certain optics problems, computer graphics modeling of three-dimensional shapes, and the design and analysis of specialized microwave reflector antennas.

Barger, R. L.↗

Comparison of Hirs' equation of Moody's equation for determining rotordynamic coefficients of annular pressure seals

The rotordynamic coefficients of an incompressible-flow annular pressure seal were determined using a bulk-flow model in conjunction with two different friction factor relationships. The first, Hirs' equation, assumes the friction factor is a function of Reynolds number only. The second, Moody's equation, approximates Moody's diagram and assumes the friction factor is a function of both Reynolds number and relative roughness. For each value of relative roughness, Hirs' constants were determined so that both equations gave the same magnitude and slope of the friction factor. For smooth seals, both relationships give the same results. For rough seals (e/2 H sub 0 = 0.05) Moody's equation predicts 44% greater direct stiffness, 35% greater cross-coupled stiffness, 19% smaller cross-coupled damping, 59% smaller cross-coupled inertia, and nominally the same direct damping and direct inertia.

Nelson, Clayton C.↗

Comparison of Hirs' equation with Moody's equation for determining rotordynamic coefficients of annular pressure seals

The rotordynamic coefficients of an incompressible-flow annular pressure seal were determined using a bulk-flow model in conjunction with two different friction factor relationships. The first, Hirs' equation, assumes the friction factor is a function of Reynolds number only. The second, Moody's equation, approximates Moody's diagram and assumes the friction factor is a function of both Reynolds number and relative roughness. For each value of relative roughness, Hirs' constants were determined so that both equations gave the same magnitude and slope of the friction factor. For smooth seals, both relationships give the same results. For rough seals, Moody's equation predicts 44 percent greater direct stiffness, 35 percent greater cross-coupled stiffness, 19 percent smaller cross-coupled damping, 59 percent smaller cross-coupled inertia, and nominally the same direct damping and direct inertia.

Nelson, C. C.↗

Three-dimensional multigrid algorithms for the flux-split Euler equations

The Full Approximation Scheme (FAS) multigrid method is applied to several implicit flux-split algorithms for solving the three-dimensional Euler equations in a body fitted coordinate system. Each of the splitting algorithms uses a variation of approximate factorization and is implemented in a finite volume formulation. The algorithms are all vectorizable with little or no scalar computation required. The flux vectors are split into upwind components using both the splittings of Steger-Warming and Van Leer. The stability and smoothing rate of each of the schemes are examined using a Fourier analysis of the complete system of equations. Results are presented for three-dimensional subsonic, transonic, and supersonic flows which demonstrate substantially improved convergence rates with the multigrid algorithm. The influence of using both a V-cycle and a W-cycle on the convergence is examined.

Anderson, W. Kyle↗

Morphology of galaxy mergers at intermediate redshift

Using a combination of numerical simulation and synthesized Hubble Space Telescope Wide Field Camera 2 (HST WFC2) images, we follow the detectability of morphological signatures of disk galaxy mergers at intermediate redshifts z = 0.4 and z = 1.0. Rapid evolution in the surface brightness of tidal tails makes their use as an interaction signature limited; long tidal tails will only be detectable relatively early in the encounter, before the galaxies actually merge. Violent relaxation rapidly smooths isophotal irregularities after the galaxies merge, leaving the very low surface brightness tidal debris surrounding the remnant as the only sign of peculiarity after a few hundred Myr. This debris becomes indistinguishable from the main body of the remnant after less than or approximately equal to 0.2 (z = 1.0) to 1 (z = 0.4) Gyr, leaving a morphologically normal elliptical galaxy behind. Using WFC1, all signatures become undetectable almost immediately once the galaxies merge, due to the aberrated Point Speed Function (PSF). We interpret these results in light of recent HST programs investigating the nature of starburst and E + A galaxies in moderate redshift clusters.

Mihos, J. Christopher↗

Light-Toned, Layered Outcrops of Northern Terra Meridiani Mars: Viking, Phobos 2, and Mars Global Surveyor Observations

Locating outcrops of sedimentary rock on Mars is an important step toward deciphering the planet's geologic and climatologic record. Sedimentary rock representing the earliest martian environments, are of particular interest in this context. This is a report about a vast exposure of material proposed to be martian sedimentary rock. The outcrops cover an area (approximately sq 300,000 km) roughly the size of the Colorado Plateau in North America (approximately 260,000 sq km). The materials occur in northern Terra Meridiani, near of one of the four sites being considered for a 2004 NASA Mars Exploration Rover (MER) landing. The landing ellipse, centered at deg S, deg W, lies in a region exhibiting smooth and rough (at meter scale) dark-toned surfaces, with scattered light-toned patches. Stratigraphically, the dark-toned materials at the MER site lie unconformably on top of a previously-eroded, light-toned surface; the light-toned patches in the landing ellipse are geologic windows down to this lower stratigraphic unit. North of the landing ellipse, the light-toned materials are well-exposed because the darker materials have been removed, stranding outlier remnants in a few locations. The light-toned materials are layered, vertically heterogeneous, and exhibit lateral continuity over hundreds of kilometers. Eroded layers produce cliffs; some outcrops are expressed as mesas, buttes, and spires; and impact craters ranging in diameter from a few meters to tens of kilometers are interbedded with the layers. The purpose of this report is to summarize the results of greater than 6 years of photogeologic investigation into the nature of the light-toned outcrops of northern Terra Meridiani. The work is a 'snapshot' of progress made toward eventual geologic mapping and establishment of the stratigraphic sequence for the materials through 30 September 2002, a day prior to the first release of Mars Odyssey Thermal Emission Imaging System (THEMIS) data to the NASA Planetary Data System (PDS). The main body of data examined were Mars Global Surveyor (MGS) Mars Orbiter Camera (MOC) images acquired through 30 September 2002. The data also 2 include Viking orbiter images, a Phobos 2 Termoscan image, MGS Mars Orbiter Laser Altimeter (MOLA) topographic observations, and the products of published Viking Infrared Thermal Mapper (IRTM) and Mars Global Surveyor (MGS) Thermal Emission Spectrometer (TES) analyses. Through September 2002, over 126,000 MOC images had been acquired, and greater than 600 of the MOC narrow angle (1.5-12 m/pixel) images occur within the portions of Terra Meridiani and southwestern Arabia Terra.

Edgett, Kenneth S.↗

Global Stratigraphy of Venus: Analysis of a Random Sample of Thirty-Six Test Areas

The age relations between 36 impact craters with dark paraboloids and other geologic units and structures at these localities have been studied through photogeologic analysis of Magellan SAR images of the surface of Venus. Geologic settings in all 36 sites, about 1000 x 1000 km each, could be characterized using only 10 different terrain units and six types of structures. These units and structures form a major stratigraphic and geologic sequence (from oldest to youngest): (1) tessera terrain; (2) densely fractured terrains associated with coronae and in the form of remnants among plains; (3) fractured and ridged plains and ridge belts; (4) plains with wrinkle ridges; (5) ridges associated with coronae annulae and ridges of arachnoid annulae which are contemporary with wrinkle ridges of the ridged plains; (6) smooth and lobate plains; (7) fractures of coronae annulae, and fractures not related to coronae annulae, which disrupt ridged and smooth plains; (8) rift-associated fractures; and (9) craters with associated dark paraboloids, which represent the youngest 1O% of the Venus impact crater population (Campbell et al.), and are on top of all volcanic and tectonic units except the youngest episodes of rift-associated fracturing and volcanism; surficial streaks and patches are approximately contemporary with dark-paraboloid craters. Mapping of such units and structures in 36 randomly distributed large regions (each approximately 10(exp 6) sq km) shows evidence for a distinctive regional and global stratigraphic and geologic sequence. On the basis of this sequence we have developed a model that illustrates several major themes in the history of Venus. Most of the history of Venus (that of its first 80% or so) is not preserved in the surface geomorphological record. The major deformation associated with tessera formation in the period sometime between 0.5-1.0 b.y. ago (Ivanov and Basilevsky) is the earliest event detected. In the terminal stages of tessera fon-nation, extensive parallel linear graben swarms representing a change in the style of deformation from shortening to extension were formed on the tessera and on some volcanic plains that were emplaced just after, and perhaps also during the latter stages of the major compressional phase of tessera emplacement. Our stratigraphic analyses suggest that following tessera formation, extensive volcanic flooding resurfaced at least 85% of the planet in the form of the presently-ridged and fractured plains. Several lines of evidence favor a high flux in the post-tessera period but we have no independent evidence for the absolute duration of ridged plains emplacement. During this time, the net state of stress in the lithosphere apparently changed from extensional to compressional, first in the form of extensive ridge belt development, followed by the formation of extensive wrinkle ridges on the flow units. Subsequently, there occurred local emplacement of smooth and lobate plains units which are presently essentially undefortned. The major events in the latest 10% of the presently preserved history of Venus (less than 50 m.y. ago) are continued rifting and some associated volcanism, and the redistribution of eolian material largely derived from impact crater deposits.

Basilevsky, Alexander T.↗

Thermal and nonthermal phenomena in solar flare loops at 20 cm wavelength and in X-rays

X-ray images from the P78-1 satellite are presented for a long-lasting burst at 20 cm wavelength mapped with the Very Large Array on May 19, 1979. The decimeter wave observations were originally interpreted in terms of two models, one invoking thermal electrons radiating at low harmonics of the gyrofrequency, and the other invoking mildly relativistic electrons emitting gyrosynchrotron radiation. If the 20 cm source is thermal, it should also be visible in soft X-rays, while if it is nonthermal, the soft X-ray emission should be weak or spatially or temporally distinct from the 20 cm burst. It is found that only one of the three 20 cm sources was approximately co-spatial with the soft X-ray source, and that it was only partially thermal. The 20 cm burst is therefore primarily decimeter type IV emission from mildly relativistic electrons of the post-flare phase. The over 2 hr lifetime and smooth temporal variation of the burst belie its nonthermal nature and suggest continuous acceleration as well as long-term storage of energetic electrons.

Schmahl, E. J.↗

Squeezed states and path integrals

The continuous-time regularization scheme for defining phase-space path integrals is briefly reviewed as a method to define a quantization procedure that is completely covariant under all smooth canonical coordinate transformations. As an illustration of this method, a limited set of transformations is discussed that have an image in the set of the usual squeezed states. It is noteworthy that even this limited set of transformations offers new possibilities for stationary phase approximations to quantum mechanical propagators.

Daubechies, Ingrid↗

Topology of large-scale structure in seeded hot dark matter models

The topology of the isodensity surfaces in seeded hot dark matter models, in which static seed masses provide the density perturbations in a universe dominated by massive neutrinos is examined. When smoothed with a Gaussian window, the linear initial conditions in these models show no trace of non-Gaussian behavior for r0 equal to or greater than 5 Mpc (h = 1/2), except for very low seed densities, which show a shift toward isolated peaks. An approximate analytic expression is given for the genus curve expected in linear density fields from randomly distributed seed masses. The evolved models have a Gaussian topology for r0 = 10 Mpc, but show a shift toward a cellular topology with r0 = 5 Mpc; Gaussian models with an identical power spectrum show the same behavior.

Beaky, Matthew M.↗

Direct Numerical Simulation of Incompressible Pipe Flow Using a B-Spline Spectral Method

A numerical method based on b-spline polynomials was developed to study incompressible flows in cylindrical geometries. A b-spline method has the advantages of possessing spectral accuracy and the flexibility of standard finite element methods. Using this method it was possible to ensure regularity of the solution near the origin, i.e. smoothness and boundedness. Because b-splines have compact support, it is also possible to remove b-splines near the center to alleviate the constraint placed on the time step by an overly fine grid. Using the natural periodicity in the azimuthal direction and approximating the streamwise direction as periodic, so-called time evolving flow, greatly reduced the cost and complexity of the computations. A direct numerical simulation of pipe flow was carried out using the method described above at a Reynolds number of 5600 based on diameter and bulk velocity. General knowledge of pipe flow and the availability of experimental measurements make pipe flow the ideal test case with which to validate the numerical method. Results indicated that high flatness levels of the radial component of velocity in the near wall region are physical; regions of high radial velocity were detected and appear to be related to high speed streaks in the boundary layer. Budgets of Reynolds stress transport equations showed close similarity with those of channel flow. However contrary to channel flow, the log layer of pipe flow is not homogeneous for the present Reynolds number. A topological method based on a classification of the invariants of the velocity gradient tensor was used. Plotting iso-surfaces of the discriminant of the invariants proved to be a good method for identifying vortical eddies in the flow field.

Loulou, Patrick↗

Heat Transfer Through Turbulent Friction Layers

The "general Prandtl number" Pr(exp 1) - A(sub q)/A Pr, aside from the Reynolds number determines the ratio of turbulent to molecular heat transfer, and the temperature distribution in turbulent friction layers. A(sub q) = exchange coefficient for heat; A = exchange coefficient for momentum transfer. A formula is derived from the equation defining the general Prandtl number which describes the temperature as a function of the velocity. For fully developed thermal boundary layers all questions relating to heat transfer to and from incompressible fluids can be treated in a simple manner if the ratio of the turbulent shear stress to the total stress T(sub t)/T in the layers near the wall is known, and if the A(sub q)/A can be regarded as independent of the distance from the wall. The velocity distribution across a flat smooth channel and deep into the laminar sublayer was measured for isothermal flow to establish the shear stress ratio T(sub t)/T and to extend the universal wall friction law. The values of T(sub t)/T which resulted from these measurements can be approximately represented by a linear function of the velocity in the laminar-turbulent transition zone. The effect of the temperature relationship of the material values on the flow near the wall is briefly analyzed. It was found that the velocity at the laminar boundary (in contrast to the thickness of the laminar layer) is approximately independent of the temperature distribution. The temperature gradient at the wall and the distribution of temperature and heat flow in the turbulent friction layers were calculated on the basis of the data under two equations. The derived formulas and the figures reveal the effects of the Prandtl number, the Reynolds number, the exchange quantities and the temperature relationship of the material values.

Reichardt, H.↗