Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Bounded motion”

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 235 records · Page 13

Use of satellite data in a diagnostic parameterization of convective heating

Heating estimates derived from a diagnostic technique using observed rainfall and GOES IR digital imagery were completed and evaluated for accuracy. A sensitivity analysis was done to examine assumptions regarding shape of the normalized mass flux profile, cloud precipitation efficiency, and existence of convective scale downdrafts. The results, which were derived using what are felt to be bounding limits of the assumptions, indicate that the heating estimates are reliable for use in diagnostic available potential energy (APE) budgets. Comparison to heating estimates derived as residuals in the thermodynamic equation show the level of maximum heating (near 300 mb) to be the same on a time averaged basis. Heating estimates were used to study the response of the large scale environment to the cumulus scale thermodynamic forcing. An analysis of the thermally forced component of vertical motion through the omega equation showed that a significant fraction of the total grid scale upward motion results from the heating provided by the condensation and vertical eddy heat transport in the convective cells.

Robertson, F. R.↗

Convective effects in float-zone and Czochralski melts

The hydrodynamics of crystal-growth melts is a relatively new research area. Numerical modeling of these processes is necessary. The work discussed herein is in two parts: numerical simulations of the flow in a Czochralski melt, and also of that in a float zone. In addition, for the float-zone case, energy stability theory will be used to determine stability bounds for the onset of oscillatory thermo-capillary flow. Convective effects in crystal-growth melts arise from a variety of mechanisms. Temperature gradients both in the direction of gravity and normal to it give rise to convection due to buoyancy effects. Rotation of the crucible and/or crystal causes a forced convection which may augment or oppose the buoyancy-driven flow. Finally, thermo-capillary forces (due to the variation of surface tension with temperature) drive surface motions which in turn generate convection in the bulk fluid. All of these mechanisms are present in either Czochralski or float-zone growth. The objective of the Czochralski modeling is to develop an accurate numerical simulation of the flow in a Czochralski silicon melt and to investigate the effects of various parameters on the flow properties. Like some earlier investigations, the intent is to simulate the effects of buoyancy, forced and thermo-capillary convection, including unsteady effects. Unlike earlier work, the aim is to include the effects of a variable free surface and freezing interface and, possibly incorporate nonaxisymmetric effects.

Neitzel, G. P.↗

Collision detection for spacecraft proximity operations

Collision Detection for Spacecraft Proximity Operations This thesis describes the development of a new collision detection algorithm to be used when two spacecraft are operating in the same vicinity. The two spacecraft are modelled as unions of convex polyhedra, where the polyhedron resulting from the union may be either convex or nonconvex. The relative motion of the two spacecraft is assumed to be such that one vehicle is moving with constant linear and angular velocity with respect to the other. The algorithm determines if a collision is possible and, if so, predicts the time when the collision will take place. The theoretical basis for the new collision detection algorithm is the C-function formulation of the configuration space approach recently introduced by researchers in robotics. Three different types of C-functions are defined that model the contacts between the vertices, edges, and faces of the polyhedra representing the two spacecraft. These C-functions are used to formulate three "collision" conditions. The first of these conditions limits the points representing potential collisions to the zeros of the C-functions. The new algorithm is fundamentally a search for the smallest zero of any C-function that satisfies the second and third collision conditions. The C-functions are shown to be transcendental functions of time for the assumed trajectory of the moving spacecraft. The zeros of these functions cannot be expressed in dosed form. Therefore, numerical search procedures are developed to find aLl of the zeros of a C-function in specified bounded intervals of time. These bounded intervals of time are found by examining the second and third collision conditions. The capabilities of the new algorithm are demonstrated for several example cases. These include examples of collisions determined by zeros of each of the three different types of C-functions. In addition to predicting the time of first contact of the polyhedra, the algorithm identifies the features of the two polyhedra that are touching at this time. The new collision detection algorithm is the first such algorithm that is capable of solving the collision detection problem exactly for the case where the moving object has constant linear and angular velocities. This is a significant improvement on previous collision detection algorithms described in the literature. In particular, the ability to handle constant angular velocity represents a more realistic type of rotational motion than those which have been used in other algorithms.

Robin M Vaughan↗

Some new results on spherical dynamos

Small-scale motions in the earth's liquid core are likely to be highly anisotropic because of the effects of rotation. Guided by physical considerations, the 'alpha-effect' described by an anisotropic tensor alpha sub ik is formulated and the corresponding boundary value problem for a sphere is solved for a variety of boundary conditions. A converged solution has been obtained only in the case of the Fermi condition of an infinitely conducting exterior of the sphere. Some remarks are made on the hypothetical upper bound on magnetic field strengths in planetary cores originally proposed by Busse (1976).

Busse, F. H.↗

Trompe L'Oeil 'binary' pulsars

A freely precessing pulsar produces pulse phase residuals which can mimic those of a pulsar in a binary orbit. In particular, discrete sets of phase residuals due to precessional motion of an isolated pulsar are sampled; it is shown that this data is well fit by residuals from a binary pulsar in a sufficiently tight orbit. Analytic and numerical relationships between the projected orbital size, a(p) sin i, and the orbital eccentricity, e, of a misidentified binary pulsar; are found the observations that would distinguish between these models are discussed. Regardless of the mechanism that causes the precession, the maximum amplitude of the phase residual is pi/2: consequently, a(p)sin i is (approximately) bounded by cP(puls)/4. The newly discovered 'binary' millisecond pulsars in the globular cluster 47 Tuc is discussed, and it is shown that the periodic frequency modulation reported cannot be explained by free precession.

Nelson, Robert W.↗

Time Dependence of Collision Probabilities During Satellite Conjunctions

The NASA Conjunction Assessment Risk Analysis (CARA) team has recently implemented updated software to calculate the probability of collision (P (sub c)) for Earth-orbiting satellites. The algorithm can employ complex dynamical models for orbital motion, and account for the effects of non-linear trajectories as well as both position and velocity uncertainties. This “3D P (sub c)” method entails computing a 3-dimensional numerical integral for each estimated probability. Our analysis indicates that the 3D method provides several new insights over the traditional “2D P (sub c)” method, even when approximating the orbital motion using the relatively simple Keplerian two-body dynamical model. First, the formulation provides the means to estimate variations in the time derivative of the collision probability, or the probability rate, R (sub c). For close-proximity satellites, such as those orbiting in formations or clusters, R (sub c) variations can show multiple peaks that repeat or blend with one another, providing insight into the ongoing temporal distribution of risk. For single, isolated conjunctions, R (sub c) analysis provides the means to identify and bound the times of peak collision risk. Additionally, analysis of multiple actual archived conjunctions demonstrates that the commonly used “2D P (sub c)” approximation can occasionally provide inaccurate estimates. These include cases in which the 2D method yields negligibly small probabilities (e.g., P (sub c)) is greater than 10 (sup -10)), but the 3D estimates are sufficiently large to prompt increased monitoring or collision mitigation (e.g., P (sub c) is greater than or equal to 10 (sup -5)). Finally, the archive analysis indicates that a relatively efficient calculation can be used to identify which conjunctions will have negligibly small probabilities. This small-P (sub c) screening test can significantly speed the overall risk analysis computation for large numbers of conjunctions.

Hall, Doyle T.↗

Study of Magnetic Motions in the Solar Photosphere and their Implications for Heating the Solar Atmosphere

We continued our program of CO observations with the McMath-Pierce facility at Kitt Peak National Solar Observatory. Uitenbroek has developed a two-and three dimensional radiative transfer code that now includes chemical equilibrium calculations. This code allows us to compute a CO spectrum from for instance a snapshot of a solar granulation simulation (e.g. Stein & Nordlund 1989, Apj 342, L95) and compare these theoretical spectra with our spatially resolved CO spectroscopy. Van Ballegooijen and Uitenbroek have started calculations of two-dimensional fluxtube models that account consistently for hydrogen ionization in the calculation of the electron density. To this end we solve radiative transfer for hydrogen (bound-bound and bound-free transitions) in the two-dimensional models, including the effect of partial frequency redistribution (PRD) in the Lyman (alpha) and (beta) lines. From our internally consistent models we will calculate emergent spectra and the way these vary with location of some well-known spectral diagnostics and compare our results with observed line profiles. We can readily compare theoretical CO profiles from out models with our spatially resolved CO observations. Also we can compare with spatially resolved Ca II (ground based) and Mg II (we have observations done with the UVSP/SMM instrument), and Lyman (alpha) observations that should be available from SUMER/SOHO.

Noyes, Robert W.↗

A Candidate Wake Vortex Strength Definition for Application to the NASA Aircraft Vortex Spacing System (AVOSS)

A significant effort is underway at NASA Langley to develop a system to provide dynamical aircraft wake vortex spacing criteria to Air Traffic Control (ATC). The system under development, the Aircraft Vortex Spacing System (AVOSS), combines the inputs of multiple subsystems to provide separation matrices with sufficient stability for use by ATC and sufficient monitoring to ensure safety. The subsystems include a meteorological subsystem, a wake behavior prediction subsystem, a wake sensor subsystem, and system integration and ATC interfaces. The proposed AVOSS is capable of using two factors, singly or in combination, for reducing in-trail spacing. These factors are wake vortex motion out of a predefined approach corridor and wake decay below a strength that is acceptable for encounter. Although basic research into the wake phenomena has historically used wake total circulation as a strength parameter, there is a requirement for a more specific strength definition that may be applied across multiple disciplines and teams to produce a real-time, automated system. This paper presents some of the limitations of previous applications of circulation to aircraft wake observations and describes the results of a preliminary effort to bound a spacing system strength definition.

Hinton, David A.↗

The Effect of a Rotating Magnetic Field on Rayleigh-Benard Convection

A destabilizing vertical temperature gradient and a rotating magnetic field have been applied to a cylindrical column of liquid gallium. The convective flows which arise as a function of these driving parameters are identified. For small magnetic field strengths, a region of stationary flow is observed. This region is bounded by critical values of the Rayleigh and Hartmann numbers. As the rotating magnetic field is increased, the critical Rayleigh number can increase by more than a factor of 10. The rotating magnetic field itself induces an instability at a critical value of the Hartmann number independent of the Rayleigh number. The nature of the bifurcations (whether subcritical or supercritical) and the convective flows occurring at the critical Rayleigh numbers are dependent upon the Hartmann number. In the stationary flow region, the rotating magnetic field can induce fluid motion with velocities on the order of cm/s. Thus, a rotating magnetic field can be used to achieve the benefits of forced convection during crystal growth processes without triggering deleterious instabilities.

Volz, Martin P.↗

On Target to Mars

This viewgraph presentation reviews the use of Descent Image Motion Estimation System (DIMES) for the descent of a spacecraft onto the surface of Mars. In the past this system was used to assist in the landing of the MER spacecraft. The overall algorithm is reviewed, and views of the hardware, and views from Spirit's descent are shown. On Spirit, had DIMES not been used, the impact velocity would have been at the limit of the airbag capability and Spirit may have bounced into Endurance Crater. By using DIMES, the velocity was reduced to well within the bounds of the airbag performance and Spirit arrived safely at Mars. Views from Oppurtunity's descent are also shown. The system to avoid and detect hazards is reviewed next. Landmark Based Spacecraft Pinpoint Landing is also reviewed. A cartoon version of a pinpoint landing and the various points is shown. Mars s surface has a large amount of craters, which are ideal landmarks . According to literatures on Martian cratering, 60 % of Martian surface is heavily cratered. The ideal (craters) landmarks for pinpoint landing will be between 1000 to 50 meters in diagonal The ideal altitude for position estimation should greater than 2 km above the ground. The algorithms used to detect and match craters are reviewed.

Mars↗

Large eddy simulation of incompressible turbulent channel flow

The three-dimensional, time-dependent primitive equations of motion were numerically integrated for the case of turbulent channel flow. A partially implicit numerical method was developed. An important feature of this scheme is that the equation of continuity is solved directly. The residual field motions were simulated through an eddy viscosity model, while the large-scale field was obtained directly from the solution of the governing equations. An important portion of the initial velocity field was obtained from the solution of the linearized Navier-Stokes equations. The pseudospectral method was used for numerical differentiation in the horizontal directions, and second-order finite-difference schemes were used in the direction normal to the walls. The large eddy simulation technique is capable of reproducing some of the important features of wall-bounded turbulent flows. The resolvable portions of the root-mean square wall pressure fluctuations, pressure velocity-gradient correlations, and velocity pressure-gradient correlations are documented.

Moin, P.↗

Investigation of thermospheric winds relative to space station orbital altitudes

An investigation of thermospheric winds, relative to the space station orbital altitudes, was made in order to provide information that is useful in an environmental disturbance assessment. Current plans are for this low Earth orbiting facility to orbit at an inclination of 28.5 deg. The orbital altitudes were not yet defined due to the evolutionary configuration of the Space Station. The upper and lower bounds of the orbital altitudes will be based on constraints set by the drag and expected orbital decay and delivery altitude capability of the Shuttle. The orbital altitude will be estimated on the order of 500 km. Neutral winds in the region from about 80 to 600 km which were derived from satellite drag data, Fabry-Perot interferometers, sounding rockets, ground-based optical Doppler techniques, incoherent scatter radar measurements from Millstone Hill combined with the mass spectrometer and lithium trail neutral wind measurements are examined. The equations of motion of the low Earth orbiting facility are also discussed.

Susko, M.↗

Buoyancy Driven Shear Flows of Bubble Suspensions

In this work the gas volume fraction and the root-mean-squared fluid velocity are measured in buoyancy driven shear flows of bubble suspensions in a tall, inclined, rectangular channel. The experiments are performed under conditions where We << 1a nd Re >> 1, for which comparisons are made with kinetic theory and numerical simulations. Here Re = gamma(a(exp 2)/nu is the Reynolds number and We = rho(gamma(exp 2))a(exp 3)/sigma is the Weber number; gamma is the shear rate, a is the bubble radius, nu is the kinematic viscosity of the liquid, rho is the density of the liquid, and sigma is the surface tension of the gas/liquid interface. Kang et al. calculated the bubble phase pressure and velocity variance of sheared bubble suspensions under conditions where the bubbles are spherical and the liquid phase velocity field can be approximated using potential flow theory, i.e. We= 0 and Re >> 1. Such conditions can be achieved in an experiment using gas bubbles, with a radius of O(0.5mm), in water. The theory requires that there be no average relative motion of the gas and liquid phases, hence the motivation for an experimental program in microgravity. The necessity of performing preliminary, Earth based experiments, however, requires performing experiments where the gas phase rises in the liquid, which significantly complicates the comparison of experiments with theory. Rather than comparing experimental results with theory for a uniform, homogeneous shear flow, experiments can be compared directly with solutions of the averaged equations of motion for bubble suspensions. This requires accounting for the significant lift force acting on the gas phase when the bubbles rise parallel to the average velocity of the sheared suspension. Shear flows can be produced in which the bubble phase pressure gradient, arising from shear induced collisions amongst the bubbles, balances a body force (centrifugal or gravitational) on the gas phase. A steady, non-uniform gas volume fraction can be measured, from which the bubble phase pressure gradient can be obtained and compared to theory and numerical simulations. The presence of bounding walls further complicates the experiments, since the detailed interactions of the bubbles with bounding walls is not well understood, especially in the presence of gravity, where the momentum and energy exchange depends on the inclination of the wall.

Koch, D. L.↗

Computer Modeling of the Earliest Cellular Structures and Functions

In the absence of extinct or extant record of protocells (the earliest ancestors of contemporary cells). the most direct way to test our understanding of the origin of cellular life is to construct laboratory models of protocells. Such efforts are currently underway in the NASA Astrobiology Program. They are accompanied by computational studies aimed at explaining self-organization of simple molecules into ordered structures and developing designs for molecules that perform proto-cellular functions. Many of these functions, such as import of nutrients, capture and storage of energy. and response to changes in the environment are carried out by proteins bound to membrane< We will discuss a series of large-scale, molecular-level computer simulations which demonstrate (a) how small proteins (peptides) organize themselves into ordered structures at water-membrane interfaces and insert into membranes, (b) how these peptides aggregate to form membrane-spanning structures (eg. channels), and (c) by what mechanisms such aggregates perform essential proto-cellular functions, such as proton transport of protons across cell walls, a key step in cellular bioenergetics. The simulations were performed using the molecular dynamics method, in which Newton's equations of motion for each item in the system are solved iteratively. The problems of interest required simulations on multi-nanosecond time scales, which corresponded to 10(exp 6)-10(exp 8) time steps.

Pohorille, Andrew↗

Aircraft agility maneuvers

A new dynamic model for aircraft motions is presented. This model can be viewed as intermediate between a point-mass model, in which the body attitude angles are control-like, and a rigid-body model, in which the body-attitude angles evolve according to Newton's Laws. Specifically, consideration is given to the case of symmetric flight, and a model is constructed in which the body roll-rate and the body pitch-rate are the controls. In terms of this body-rate model a minimum-time heading change maneuver is formulated. When the bounds on the body-rates are large the results are similar to the point-mass model in that the model can very quickly change the applied forces and produce an acceleration to turn the vehicle. With finite bounds on these rates, the forces change in a smooth way. This leads to a measurable effect of agility.

Cliff, Eugene M.↗

Quasi-relativistic electron precipitation due to interactions with coherent VLF waves in the magnetosphere

The equations of motion for the cyclotron resonance interaction between coherent whistler mode waves and energetic particles are rederived with the inclusion of relativistic effects. The pitch angle scattering of the near-loss-cone quasi-relativistic electrons trapped in the magnetosphere is studied using a test particle method employing these relativistic equations, and the precipitated energy spectrum due to the wave-induced perturbations of a full distribution of particles is computed. Results show that the full width at half maximum peak width of the rms scattering pattern of the near-loss-cone particles would give an upper bound to the peak width of the associated precipitated energy spectrum under the conditions of moderate wave intensities in the low L shell region. In addition, it is found that the peak widths are within the upper limit values measured by recent satellite experiments. It is concluded that interactions of inner radiation belt particles with monochromatic waves could produce precipitated fluxes with relatively sharp spectral widths, and that therefore the L-dependent narrow peaks observed by low altitude satellite particle detectors could be caused by such interactions.

Chang, H. C.↗

Evaluation of an elastic curve technique for finding the auroral oval from satellite images automatically

An automated algorithm for finding the inner boundary based on recently proposed computer vision technique is described. The algorithm is analogous to solving the equations of motion for an elastic curve, where the forces are provided by the image. The resulting equilibrium position of the elastic curve provides an automated method for finding the shape and location of the inner boundary of the auroral oval. Two methods for the evaluation of the automated algorithm, both based on the comparisons with manual measurements, are developed. The first method compares the areas within the automated and the manual boundaries. The second method measures the overlap between the interiors of the two boundaries. The expected variation between two sets of manual measurements is used to set an upper bound to the allowed discrepancy between the automated results and a single set of manual measurements. The algorithm, when tested with 71 satellite images, is found to perform best for those images without overlap between the aurora and the dayside hemisphere.

Samadani, Ramin↗

Transmission-zero bounds for large space structures, with applications

Many large space structure control problems lead quite naturally to the application of an optimal regulator, so the transmission zeros of the open-loop system give fundamental information about the speed of response achievable by the closed-loop system. Despite the importance of this and other well-known zeros properties, little attention has been given to the transmission zeros of large space structures, except for the special case of a rigid spacecraft with flexible appendages. The object of this paper is to remedy this deficiency. In particular, it is proved that the zeros of a structure with colocated sensors and actuators must lie in a region of the complex plane that is defined by its natural frequencies and damping ratios. This generic result, a consequence of the special form of the equations of motion of structural dynamics, admits a very simple graphical interpretation: it is the generalization of the classical pole-zero interlacing property of undamped single-input/single-output structures. The number of sensor/actuator pairs, and their locations, specify where in the permissible region transmission zeros actually lie, thus quantifying the effect of sensor/actuator placement on closed-loop system performance. These points are illustrated by simple examples.

Williams, Trevor↗