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 685 records · Page 38

Cell prestress. II. Contribution of microtubules

The tensegrity model hypothesizes that cytoskeleton-based microtubules (MTs) carry compression as they balance a portion of cell contractile stress. To test this hypothesis, we used traction force microscopy to measure traction at the interface of adhering human airway smooth muscle cells and a flexible polyacrylamide gel substrate. The prediction is that if MTs balance a portion of contractile stress, then, upon their disruption, the portion of stress balanced by MTs would shift to the substrate, thereby causing an increase in traction. Measurements were done first in maximally activated cells (10 microM histamine) and then again after MTs had been disrupted (1 microM colchicine). We found that after disruption of MTs, traction increased on average by approximately 13%. Because in activated cells colchicine induced neither an increase in intracellular Ca(2+) nor an increase in myosin light chain phosphorylation as shown previously, we concluded that the observed increase in traction was a result of load shift from MTs to the substrate. In addition, energy stored in the flexible substrate was calculated as work done by traction on the deformation of the substrate. This result was then utilized in an energetic analysis. We assumed that cytoskeleton-based MTs are slender elastic rods supported laterally by intermediate filaments and that MTs buckle as the cell contracts. Using the post-buckling equilibrium theory of Euler struts, we found that energy stored during buckling of MTs was quantitatively consistent with the measured increase in substrate energy after disruption of MTs. This is further evidence supporting the idea that MTs are intracellular compression-bearing elements.

Non-NASA Center↗

Ejector Enhanced Pulsejet Based Pressure Gain Combustors: An Old Idea With a New Twist

An experimental investigation of pressure-gain combustion for gas turbine application is described. The test article consists of an off-the-shelf valved pulsejet, and an optimized ejector, both housed within a shroud. The combination forms an effective can combustor across which there is a modest total pressure rise rather than the usual loss found in conventional combustors. Although the concept of using a pulsejet to affect semi-constant volume (i.e., pressure-gain) combustion is not new, that of combining it with a well designed ejector to efficiently mix the bypass flow is. The result is a device which to date has demonstrated an overall pressure rise of approximately 3.5 percent at an overall temperature ratio commensurate with modern gas turbines. This pressure ratio is substantially higher than what has been previously reported in pulsejet-based combustion experiments. Flow non-uniformities in the downstream portion of the device are also shown to be substantially reduced compared to those within the pulsejet itself. The standard deviation of total pressure fluctuations, measured just downstream of the ejector was only 5.0 percent of the mean. This smoothing aspect of the device is critical to turbomachinery applications since turbine performance is, in general, negatively affected by flow non-uniformities and unsteadiness. The experimental rig will be described and details of the performance measurements will be presented. Analyses showing the thermodynamic benefits from this level of pressure-gain performance in a gas turbine will also be assessed for several engine types. Issues regarding practical development of such a device are discussed, as are potential emissions reductions resulting from the rich burning nature of the pulsejet and the rapid mixing (quenching) associated with unsteady ejectors.

Paxson, Daniel E.↗

Wrinkle ridges on Mercury and the Moon within and outside of mascons

Found on all terrestrial planets, wrinkle ridges are anticlines formed by thrust faulting and folding resulting from crustal shortening. The MErcury Surface, Space Environment, Geochemistry, and Ranging (MESSENGER)spacecraft's orbital phase returned high resolution images and topographic data of the previously unimaged northern high latitudes of Mercury where there are large expanses of smooth plains deformed by wrinkle ridges. Concurrently, the Lunar Reconnaissance Orbiter (LRO) is obtaining high resolution images and topographic data covering lunar mare wrinkle ridges. These data allow quantitative comparison of the scale of wrinkle ridges in smooth plains volcanic units on Mercury with mare wrinkle ridges. We evaluate the topographic relief of 300 wrinkle ridges within and outside of mascon basins on the Moon and Mercury. Measured wrinkle ridges range from ~112 to 776 min relief with a mean of ~350m (median=~340m, n=150) on Mercury and from ~47 to 678 min relief with a mean of ~198 m (median=~168m, n=150) on the Moon. Wrinkle ridges on Mercury thus are approximately twice as large in mean relief compared to their counterparts on the Moon. The larger scale of Mercury's wrinkle ridges suggests that their formation can be attributed, in part, to global contraction. As global contraction on the Moon is estimated to be an order of magnitude smaller than on Mercury, the smaller scale of lunar wrinkle ridges suggests they most likely form primarily by load induced subsidence of the mare basalt. Wrinkle ridges located in lunar mascon basins and in the Caloris mascon on Mercury are not statistically significantly different in relief than ridges in non-mascon regions, suggesting comparable levels of contractional strain. The fact that mascon basins do not host wrinkle ridges with greater structural relief relative to non-mascon units may indicate the critical role lithospheric thickness plays in controlling subsidence and contraction of thick volcanic sequences on the Moon and Mercury.

Lisa S. Schleicher↗

Enabling probabilistic learning on manifolds through double diffusion maps

Here, we present a generative learning framework for probabilistic sampling that extends Probabilistic Learning on Manifolds (PLoM), which is designed to generate statistically consistent realizations of a random vector in a finite-dimensional Euclidean space, informed by a (representative) set of observations. In its original form, PLoM constructs a reduced-order probabilistic model by combining three main components: (a) kernel density estimation to approximate the underlying probability measure, (b) Diffusion Maps to characterize the manifold of the data, and (c) a reduced-order Itô Stochastic Differential Equation (ISDE) to sample from the learned distribution. However, its sampling dynamics are posed in the ambient space and the retained number of reduced coordinates is chosen by projection-reconstruction error. In practice, this often (i) requires more coordinates than the data’s intrinsic dimension to achieve stable sampling and (ii) lacks a smooth, basis-independent lifting back to the data domain; moreover, standard Diffusion Maps emphasize harmonic eigenfunctions and can miss non-harmonic latent structure. We address these limitations by decoupling geometry learning from sampling: a first Diffusion Maps pass identifies non-harmonic coordinates on which we formulate a full-order ISDE directly in the latent space, while Double Diffusion Maps captures multiscale geometric features and Geometric Harmonics (GH) learns a smooth lifting map to the ambient variables that is independent of the particular diffusion basis. This hybrid design preserves the system’s dynamical richness with a compact geometric representation and enables principled out-of-sample inference. The effectiveness and robustness of the proposed method are illustrated through two numerical studies: one based on data generated from two-dimensional Hermite polynomial functions and another based on high-fidelity simulations of a detonation wave in a reactive flow.

Double diffusion maps↗

Minor planets and related objects. XV - Asteroid /1620/ Geographos

Analysis of the light curves and color-phase and polarization-phase relations of Geographos according to observations over an eight-month period in 1969. The light curves are relatively smooth with two distinct maxima and minima. The largest amplitude of 2.03 mag is the largest observed on asteroids to date. The average B-V is 0.87 and the average U-B is 0.50. A laboratory model that best reproduced Geographos' light curves was a cylinder having hemispherical ends with a length to width ratio of 2.7; on the basis of this model and a relation of Bowell and Zellner (1973), Geographos was calculated to be approximately 1.50 km wide and 4.0 km long.

Dunlap, J. L.↗

Supersonic aerodynamic characteristics of a series of wrap-around-fin missile configurations

A parametric study of wrap-around-fin missile configurations was conducted at Mach numbers from 1.60 to 2.86 in the Langley Unitary Plan wind tunnel. The fin configurations investigated included variations in chord length, leading edge sweep, thickness ratio, and leading edge shape. The investigation also included a smooth and a stepped-down afterbody required for flush retraction of the wrap-around-fin configuration. The investigation indicated no unusual longitudinal characteristics; however, all the wrap-around-fin configurations tested indicated erratic lateral behavior, particularly in the form of induced roll at zero angle of attack and irregular variations of roll with angle of attack and Mach number. The magnitude of rolling moment at an angle of attack of 0 deg is estimated to represent approximately 0.25 deg or less roll control deflection. The stepped-down afterbody has a marked effect on reducing the induced roll.

Fournier, R. H.↗

The structure of mass-loading shocks

A new two-fluid model which describes mass loading in the solar wind (e.g., the interaction of the solar wind with a cometary coma or the local interstellar medium) is presented. The self-consistent back-reaction of the mass-loaded ions is included through their effective scattering in low-frequency MHD turbulence and the invocation of a diffusive approximation. Such an approximation has the advantage of introducing self-consistent dissipation coefficients into the governing equations, thereby facilitating the investigation of the internal structure of shocks in mass-loading environments. To illustrate the utility of the new model, we consider the structure of cometary shocks in the hypersonic one-dimensional limit, finding that the incoming solar wind is slowed by both mass loading and the development of a large cometary ion pressure gradient. The shock is broadened and smoothed by the cometary ions with a thickness of the order of the cometary ion diffusion scale.

Zank, G. P.↗

Energy density of ionospheric and solar wind origin ions in the near-Earth magnetotail during substorms

Comprehensive energy density studies provide an important measure of the participation of various sources in energization processes and have been relatively rare in the literature. We present a statistical study of the energy density of the near-Earth magnetotail major ions (H(+), O(+), He(++), He(+)) during substorm expansion phase and discuss its implications for the solar wind/magnetosphere/ionosphere coupling. Our aim is to examine the relation between auroral activity and the particle energization during substorms through the correlation between the AE indices and the energy density of the major magnetospheric ions. The data we used here were collected by the charge-energy-mass (CHEM) spectrometer on board the Active Magnetospheric Particle Trace Explorer (AMPTE)/Charge Composition Explorer (CCE) satellite in the near-equatorial nightside magnetosphere, at geocentric distances approximately 7 to 9 R(sub E). CHEM provided the opportunity to conduct the first statistical study of energy density in the near-Earth magnetotail with multispecies particle data extending into the higher energy range (greater than or equal to 20 keV/E). the use of 1-min AE indices in this study should be emphasized, as the use (in previous statistical studies) of the (3-hour) Kp index or of long-time averages of AE indices essentially smoothed out all the information on substorms. Most distinct feature of our study is the excellent correlation of O(+) energy density with the AE index, in contrast with the remarkably poor He(++) energy density - AE index correlation. Furthermore, we examined the relation of the ion energy density to the electrojet activity during substorm growth phase. The O(+) energy density is strongly correlated with the pre-onset AU index, that is the eastward electrojet intensity, which represents the growth phase current system. Our investigation shows that the near-Earth magnetotail is increasingly fed with energetic ionospheric ions during periods of enhanced dissipation of auroral currents. The participation of the ionosphere in the substorm energization processes seems to be closely, although not solely, associated with the solar wind/magnetosphere coupling. That is, the ionosphere influences actively the substorm energization processes by responding to the increased solar wind/magnetosphere coupling as well as to the unloading dissipation of stored energy, with the increased feeding of new material into the magnetosphere.

Daglis, Loannis A.↗

Particle acceleration in the dynamic magnetotail: Orbits in self-consistent three-dimensional MHD fields

The acceleration of protons in a dynamically evolving magnetotail is investigated by tracing particles in the fields obtained from a three-dimensional resistive magnetohydrodynamic (MHD) simulation. The MHD simulation, representing plasmoid formation and ejection through a near-Earth reconnection process, leads to cross-tail electric fields of up to approximately 4 mV/m with integrated voltages across the tail of up to approximately 200 kV. Energization of particles takes place over a wide range along the tail, due to the large spatial extent of the increased electric field together with the finite cross-tail extent of the electric field region. Such accelerated particles appear earthward of the neutral line over a significant portion of the closed field line region inside of the separatrix, not just in the vicinity of the separatrix. Two different acceleration processes are identified: a 'quasi-potential' acceleration, due to particle motion in the direction of the cross-tail electric field, and a 'quasi-betatron' effect, which consists of multiple energy gains from repeated crossings of the acceleration region, mostly on Speiser-type orbits, in the spatially varying induced electric field. The major source region for accelerated particles in the hundreds of keV range is the central plasma sheet at the dawn flank outside the reconnection site. Since this source plasma is already hot and dense, its moderate energization by a factor of approximately 2 may be sufficient to explain the observed increases in the energetic particle fluxes. Particles from the tail are the source of beams at the plasma sheet/lobe boundary. The temporal increase in the energetic particle fluxes, estimated from the increase in energy gain, occurs on a fast timescale of a few minutes, coincident with a strong increase in B(sub z), despite the fact that the inner boundary ('injection boundary') of the distribution of energized particles is fairly smooth.

Birn, Joachim↗

A 3D High-Order Unstructured Finite-Volume Algorithm for Solving Maxwell's Equations

A three-dimensional finite-volume algorithm based on arbitrary basis functions for time-dependent problems on general unstructured grids is developed. The method is applied to the time-domain Maxwell equations. Discrete unknowns are volume integrals or cell averages of the electric and magnetic field variables. Spatial terms are converted to surface integrals using the Gauss curl theorem. Polynomial basis functions are introduced in constructing local representations of the fields and evaluating the volume and surface integrals. Electric and magnetic fields are approximated by linear combinations of these basis functions. Unlike other unstructured formulations used in Computational Fluid Dynamics, the new formulation actually does not reconstruct the field variables at each time step. Instead, the spatial terms are calculated in terms of unknowns by precomputing weights at the beginning of the computation as functions of cell geometry and basis functions to retain efficiency. Since no assumption is made for cell geometry, this new formulation is suitable for arbitrarily defined grids, either smooth or unsmooth. However, to facilitate the volume and surface integrations, arbitrary polyhedral cells with polygonal faces are used in constructing grids. Both centered and upwind schemes are formulated. It is shown that conventional schemes (second order in Cartesian grids) are equivalent to the new schemes using first degree polynomials as the basis functions and the midpoint quadrature for the integrations. In the new formulation, higher orders of accuracy are achieved by using higher degree polynomial basis functions. Furthermore, all the surface and volume integrations are carried out exactly. Several model electromagnetic scattering problems are calculated and compared with analytical solutions. Examples are given for cases based on 0th to 3rd degree polynomial basis functions. In all calculations, a centered scheme is applied in the interior, while an upwind matching scheme is employed at material interfaces and the Engquist-Majda non-reflecting boundary condition is implemented at the numerical outer boundaries. The staggered leapfrog scheme and the Runge-Kutta methods are utilized for the time integration. Excellent agreements are found between the numerical and analytical solutions.

Liu, Yen↗

Fabrication of Spherical Reflectors in Outer Space

A process is proposed for fabrication of lightweight spherical reflectors in outer space for telescopes, radio antennas, and light collectors that would be operated there. The process would obviate the relatively massive substrates and frames needed to support such reflectors in normal Earth gravitation. According to the proposal, fabrication of a reflector would begin with blowing of a bubble to the specified reflector radius. Taking advantage of the outer-space vacuum as a suitable environment for evaporative deposition of metal, a metal-evaporation source would be turned on and moved around the bubble to deposit a reflective metal film over the specified reflector area to a thickness of several microns. Then the source would be moved and aimed to deposit more metal around the edge of the reflector area, increasing the thickness there to approximately equal to 100 micron to form a frame. Then the bubble would be deflated and peeled off the metal, leaving a thin-film spherical mirror having an integral frame. The mirror would then be mounted for use. The feasibility of this technology has been proved by fabricating a prototype at JPL. As shown in the figure, a 2-in. (.5-cm) diameter hemispherical prototype reflector was made from a polymer bubble coated with silver, forming a very smooth surface.

Wang, Yu↗

Computational study of inertial migration of prolate particles in a straight rectangular channel

Inertial migration of spherical particles has been investigated extensively using experiments, theory, and computational modeling. Yet, a systematic investigation of the effect of particle shape on inertial migration is still lacking. Herein, we numerically mapped the migration dynamics of a prolate particle in a straight rectangular microchannel using smoothed particle hydrodynamics at moderate Reynolds number flows. After validation, we applied our model to 2:1 and 3:1 shape aspect ratio particles at multiple confinement ratios. Their effects on the final focusing position, rotational behavior, and transitional dynamics were studied. In addition to the commonly reported tumbling motion, for the first time, we identified a new logrolling behavior of a prolate ellipsoidal particle in the confined channel. This new behavior occurs when the confinement ratio is above an approximate threshold value of K = 0.72. Our microfluidic experiments using cell aggregates with similar shape aspect ratio and confinement ratio confirmed this new predicted logrolling motion. In this study, we also found that the same particle can undergo different rotational modes, including kayaking behavior, depending on its initial cross-sectional position and orientation. Furthermore, we examined the migration speed, angular velocity, and rotation period as well as their dependence on both particle shape aspect ratio and confinement ratio. Our findings are especially relevant to the applications where particle shape and alignment are used for sorting and analysis, such as the use of barcoded particles for biochemical assays through optical reading, or the shape-based enrichment of microalgae, bacteria, and chromosomes.

42 ENGINEERING↗

A method for the measurement and the statistical analysis of atmospheric turbulence

The instantaneous values of output voltages representing the wind velocity vector and the temperature at different elevations of the 250-foot meteorological tower located at NASA Wallops Flight Center are provided with the three dimensional split-film TSI Model 1080 anemometer system. The output voltages are sampled at a rate of one every 5 milliseconds, digitized and stored on digital magnetic tapes for a time period of approximately 40 minutes, with the use of a specially designed data acqusition system. A new calibration procedure permits the conversion of the digital voltages to the respective values of the temperature and the velocity components in a Cartesian coordinate system connected with the TSI probe with considerable accuracy. Power, cross, coincidence and quadrature spectra of the wind components and the temperature are obtained with the use of the fast Fourier transform. The cosine taper data window and ensemble and frequency smoothing techniques are used to provide smooth estimates of the spectral functions.

Tieleman, H. W.↗

Morphologic contrasts between Nirgal and Auqakuh Valles, Mars: Evidence of different crustal properties

Photoclinometric measurements were made of sidewall slopes in Nirgal and Auqakuh Valles and these results were interpreted in terms of the geologic setting and a simple geomorphic model to provide insights into the physical properties of crustal materials in these areas. Nirgal was interpreted to be a runoff channel and Auqakuh to be a fretted channel. Geomorphologic arguments for the sapping origin of Nirgal and Auqakuh Valles were presented. The morphologies of the channels, however, differ greatly: the tributaries of Nirgal end abruptly in theater-headed canyons, whereas the heads of tributaries of Auqakuh shallow gradually. The plateau surface surrounding both channels appears to be covered by smooth materials, presumably lava flows; they are continuous and uneroded in the Nirgal area, but at Auqakuh they are largely eroded and several layers are exposed that total about 200 m in thickness. For Nirgal Valles, the measurements show that sidewalls in the ralatively shallow upper reaches of the channel have average slopes near 30 degrees and, in the lower reaches, sidewall slopes exceed 50 degrees. Auqakuh, on the other hand, has maximum sidewall slopes of 14 degrees and an approximate maximum depth of 1000 m. Faint, horizontal layering in portions of the lower reaches of Nirgal may indicate inhomogeneity in either composition or topography.

Mackinnon, David J.↗

Collisional stripping of Mercury's mantle

A three-dimensional smoothed-particle hydrocode is the basis of the present numerical simulations of conditions under which a giant collision between a proto-Mercury and a planet one-sixth its size would lead to the loss of most of the silicate mantle of Mercury and thereby account for its anomalously high density. A head-on collision at 20 km/sec, and an off-axis impact parameter of half the radius of the proto-Mercury at 35 km/sec, are approximately equal in damage yielded; both will yield a remnant whose characteristics are those of the present Mercury.

Benz, Willy↗

Automatic evaluation of interferograms

A system for the evaluation of interference patterns was developed. For digitizing and processing of the interferograms from classical and holographic interferometers a picture analysis system based upon a computer with a television digitizer was installed. Depending on the quality of the interferograms, four different picture enhancement operations may be used: Signal averaging; spatial smoothing, subtraction of the overlayed intensity function and the removal of distortion-patterns using a spatial filtering technique in the frequency spectrum of the interferograms. The extraction of fringe loci from the digitized interferograms is performed by a foating-threshold method. The fringes are numbered using a special scheme after the removal of any fringe disconnections which appeared if there was insufficient contrast in the holograms. The reconstruction of the object function from the fringe field uses least squares approximation with spline fit. Applications are given.

Becker, F.↗

An Adaptively-Refined, Cartesian, Cell-Based Scheme for the Euler and Navier-Stokes Equations

A Cartesian, cell-based scheme for solving the Euler and Navier-Stokes equations in two dimensions is developed and tested. Grids about geometrically complicated bodies are generated automatically, by recursive subdivision of a single Cartesian cell encompassing the entire flow domain. Where the resulting cells intersect bodies, polygonal 'cut' cells are created. The geometry of the cut cells is computed using polygon-clipping algorithms. The grid is stored in a binary-tree data structure which provides a natural means of obtaining cell-to-cell connectivity and of carrying out solution-adaptive refinement. The Euler and Navier-Stokes equations are solved on the resulting grids using a finite-volume formulation. The convective terms are upwinded, with a limited linear reconstruction of the primitive variables used to provide input states to an approximate Riemann solver for computing the fluxes between neighboring cells. A multi-stage time-stepping scheme is used to reach a steady-state solution. Validation of the Euler solver with benchmark numerical and exact solutions is presented. An assessment of the accuracy of the approach is made by uniform and adaptive grid refinements for a steady, transonic, exact solution to the Euler equations. The error of the approach is directly compared to a structured solver formulation. A non smooth flow is also assessed for grid convergence, comparing uniform and adaptively refined results. Several formulations of the viscous terms are assessed analytically, both for accuracy and positivity. The two best formulations are used to compute adaptively refined solutions of the Navier-Stokes equations. These solutions are compared to each other, to experimental results and/or theory for a series of low and moderate Reynolds numbers flow fields. The most suitable viscous discretization is demonstrated for geometrically-complicated internal flows. For flows at high Reynolds numbers, both an altered grid-generation procedure and a different formulation of the viscous terms are shown to be necessary. A hybrid Cartesian/body-fitted grid generation approach is demonstrated. In addition, a grid-generation procedure based on body-aligned cell cutting coupled with a viscous stensil-construction procedure based on quadratic programming is presented.

Coirier, William John↗

Heat Transfer in a Complex Trailing Edge Passage for a High Pressure Turbine Blade: Simulation Results - Part 2:

A combined experimental and numerical study to investigate the heat transfer distribution in a complex blade trailing edge passage was conducted. The geometry consists of a two pass serpentine passage with taper toward the trailing edge, as well as from hub to tip. The upflow channel has an average aspect ratio of roughly 14:1, while the exit passage aspect ratio is about 5:1. The upflow channel is split in an interrupted way and is smooth on the trailing edge side of the split and turbulated on the other side. A turning vane is placed near the tip of the upflow channel. Reynolds numbers in the range of 31,000 to 61,000, based on inlet conditions, were simulated numerically. The simulation was performed using the Glenn-HT code, a full three-dimensional Navier-Stokes solver using the Wilcox k-omega turbulence model. A structured multi-block grid is used with approximately 4.5 million cells and average y+ values on the order of unity. Pressure and heat transfer distributions are presented with comparison to the experimental data. While there are some regions with discrepancies, in general the agreement is very good for both pressure and heat transfer.

Rigby, David L.↗