SEARCH · Engineering Papers
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.
A method is presented for calculating the global distribution of diabatic heating rate. Preliminary results of global heating rate evaluated from the European center for Medium Range Weather Forecasts Level IIIb analysis data is also presented.
In this paper we describe a powerful method for mapping the distribution of dust through a molecular cloud using data obtained in large-scale, multiwavelength, infrared imaging surveys. This method combines direct measurements of near-infrared color excess and certain techniques of star counting to derive mean extinctions and map the dust column density distribution through a cloud at higher angular resolutions and greater optical depths than those achieved previously by optical star counting. We report the initial results of the application of this method to a dark cloud complex near the cluster IC 5146, where we have performed coordinated, near-infrared, JHK imaging and (13)CO, C(18)O, and CS millimeter-wave, molecular-line surveys of a large portion of the complex. More than 4000 stars were detected in our JHK survey of the cloud. Of these, all but about a dozen appear to be field stars not associated with the cloud. Star count maps at J band show a striking and detailed anticorrelation between the surface density of J-band sources and CO and CS molecular-line emission. We used the (H-K) colors and positions of nearly 1300 sources to directly measure and map the extinction and thus trace the dust column density through the cloud at an effective angular resolution of 1 min .5. We report an interesting correlation between the measured dispersion in our extinction determinations and the extinction. Modeling this relation indicates that effects of small-scale cloud structure dominate the uncertainties in our measurements. Moreover, we demonstrate that such observations can be used to place constraints on the nature of the spatial distribution of extinction on scales smaller than our resolution. In particular, we show that models in which the dust is distributed uniformly or in discrete high-extinction clumps on scales smaller than (1 min .5) are inconsistent with the observations. We have derived extinctions at the same positions and at the same angular resolution (1 min .7) as our molecular-line observations. This enabled a direct comparison of (13)CO, C(18)O, and CS integrated intensities and column densities with A(sub V) for more than 500 positions in the cloud, corresponding to a range in A(sub V) between 0 to 32 mag of extinction. We found the integrated intensities of (13)CO, C(18)O, and CS to be roughly linearly correlated with extinction over different ranges of extinction. However, for all three molecules we find the scatter in the observed relations to be larger than can be accounted for by instrumental error, suggesting that there are large intrinsic variations in the abundances or excitation of the molecules through the cloud. Mean abundances for all the molecules relative to hydrogen were directly derived from the data. The ratio of (13)CO to C(18)O abundances was found to be significantly higher than the terrestrial ratio in regions where extinction is less than 10 mag. In the same region, the dispersion in the abundance ratio is also found to be very large, suggesting that the abundances of one or both molecules are very unstable even at relatively large cloud optical depths. Beyond 10 mag of extinction the abundances of both species appear very stable with their ratio close to the terrestrial value.
Two different singularity methods have been utilized to calculate the potential flow past a three dimensional non-lifting body. Two separate FORTRAN computer programs have been developed to implement these theoretical models, which will in the future allow inclusion of the fuselage effect in a pair of existing subcritical wing design computer programs. The first method uses higher order axial singularity distributions to model axisymmetric bodies of revolution in an either axial or inclined uniform potential flow. Use of inset of the singularity line away from the body for blunt noses, and cosine-type element distributions have been applied to obtain the optimal results. Excellent agreement to five significant figures with the exact solution pressure coefficient value has been found for a series of ellipsoids at different angles of attack. Solutions obtained for other axisymmetric bodies compare well with available experimental data. The second method utilizes distributions of singularities on the body surface, in the form of a discrete vortex lattice. This program is capable of modeling arbitrary three dimensional non-lifting bodies. Much effort has been devoted to finding the optimal method of calculating the tangential velocity on the body surface, extending techniques previously developed by other workers.
This project, supported by the NASA New Investigator Program, has primarily been funding a graduate student, Darren McKague. Since August 1999 Darren has been working part time at Raytheon, while continuing his PhD research. Darren is planning to finish his thesis work in May 2001, thus some of the work described here is ongoing. The proposed research was to use GOES visible and infrared imager data and SSM/I microwave data to obtain joint distributions of cirrus cloud ice mass and precipitation for a study region in the Eastern Tropical Pacific. These joint distributions of cirrus cloud and rainfall were to be compared to those from the CSU general circulation model to evaluate the cloud microphysical amd cumulus parameterizations in the GCM. Existing algorithms were to be used for the retrieval of cloud ice water path from GOES (Minnis) and rainfall from SSM/I (Wilheit). A theoretical study using radiative transfer models and realistic variations in cloud and precipitation profiles was to be used to estimate the retrieval errors. Due to the unavailability of the GOES satellite cloud retrieval algorithm from Dr. Minnis (a co-PI), there was a change in the approach and emphasis of the project. The new approach was to develop a completely new type of remote sensing algorithm - one to directly retrieve joint probability density functions (pdf's) of cloud properties from multi-dimensional histograms of satellite radiances. The usual approach is to retrieve individual pixels of variables (i.e. cloud optical depth), and then aggregate the information. Only statistical information is actually needed, however, and so a more direct method is desirable. We developed forward radiative transfer models for the SSM/I and GOES channels, originally for testing the retrieval algorithms. The visible and near infrared ice scattering information is obtained from geometric ray tracing of fractal ice crystals (Andreas Macke), while the mid-infrared and microwave scattering is computed with Mie scattering. The radiative transfer is performed with the Spherical Harmonic Discrete Ordinate Method (developed by the PI), and infrared molecular absorption is included with the correlated k-distribution method. The SHDOM radiances have been validated by comparison to version 2 of DISORT (the community "standard" discrete-ordinates radiative transfer model), however we use SHDOM since it is computationally more efficient.
The GMRES method is parallelized, and combined with local preconditioning to construct an implicit parallel solver to obtain steady-state solutions for the Navier-Stokes equations of fluid flow on distributed-memory machines. The new implicit parallel solver is designed to preserve the convergence rate of the equivalent 'serial' solver. A static domain-decomposition is used to partition the computational domain amongst the available processing nodes of the parallel machine. The SPMD (Single-Program Multiple-Data) programming model is combined with message-passing tools to develop the parallel code on a 32-node Intel Hypercube and a 512-node Intel Delta machine. The implicit parallel solver is validated for internal and external flow problems, and is found to compare identically with flow solutions obtained on a Cray Y-MP/8. A peak computational speed of 2300 MFlops/sec has been achieved on 512 nodes of the Intel Delta machine,k for a problem size of 1024 K equations (256 K grid points).
The Modified Garabedian McFadden (MGM) design procedure was incorporated into an existing 2-D multigrid Navier-Stokes airfoil analysis method. The resulting design method is an iterative procedure based on a residual correction algorithm and permits the automated design of airfoil sections with prescribed surface pressure distributions. The new design method, Multigrid Modified Garabedian McFadden (MG-MGM), is demonstrated for several different transonic pressure distributions obtained from both symmetric and cambered airfoil shapes. The airfoil profiles generated with the MG-MGM code are compared to the original configurations to assess the capabilities of the inverse design method.
Observations of precipitation over Greenland are limited. Direct precipitation measurements for the whole ice sheet are impractical, and those in the coastal region have substantial uncertainty but may be correctable with some effort. However, the analyzed wind, geopotential height and moisture fields are available for recent years, and the precipitation is retrievable from these fields by a dynamic method. Based on recent Greenland precipitation from dynamic studies, several deficiencies in the precipitation spatial distributions from these dynamic methods were evaluated by Bromwich et al.
The statistical treatment of univariate censored data is discussed. A heuristic derivation of the Kaplan-Meier maximum-likelihood estimator from first principles is presented which results in an expression amenable to analytic error analysis. Methods for comparing two or more censored samples are given along with simple computational examples, stressing the fact that most astronomical problems involve upper limits while the standard mathematical methods require lower limits. The application of univariate survival analysis to six data sets in the recent astrophysical literature is described, and various aspects of the use of survival analysis in astronomy, such as the limitations of various two-sample tests and the role of parametric modelling, are discussed.
This work formulates a method for the modeling of material damping characteristics in distributed parameter models which may be easily applied to models such as rod, plate, and beam equations. The general linear boundary value vibration equation is modified to incorporate hysteresis effects represented by complex stiffness using the transfer function approach proposed by Golla and Hughes. The governing characteristic equations are decoupled through separation of variables yielding solutions similar to those of undamped classical theory, allowing solution of the steady state as well as transient response. Example problems and solutions are provided demonstrating the similarity of the solutions to those of the classical theories and transient responses of nonviscous systems.
A method is presented for the optimization of the lift distribution across the wing of an aircraft in formation flight. The usual elliptical distribution is no longer optimal for the trailing wing in the formation due to the asymmetric nature of the encountered flow field. Control surfaces along the trailing edge of the wing can be configured to obtain a non-elliptical profile that is more optimal in terms of minimum drag. Due to the difficult-to-predict nature of formation flight aerodynamics, a Newton-Raphson peak-seeking controller is used to identify in real time the best aileron and flap deployment scheme for minimum total drag. Simulation results show that the peak-seeking controller correctly identifies an optimal trim configuration that provides additional drag savings above those achieved with conventional anti-symmetric aileron trim.
Presented is a statistical approach to measuring fault latency in a digital processor. The method relies on the use of physical fault injection where the duration of the fault injection can be controlled. Although a specific fault's latency period is never directly measured, the method indirectly determines the distribution of fault latency.
Researchers propose a method for recognizing two-dimensional shapes (hand-drawn characters, for example) with an associative memory. The method consists of two stages: first, the image is preprocessed to extract tangents to the contour of the shape; second, the set of tangents is converted to a long bit string for recognition with sparse distributed memory (SDM). SDM provides a simple, massively parallel architecture for an associative memory. Long bit vectors (256 to 1000 bits, for example) serve as both data and addresses to the memory, and patterns are grouped or classified according to similarity in Hamming distance. At the moment, tangents are extracted in a simple manner by progressively blurring the image and then using a Canny-type edge detector (Canny, 1986) to find edges at each stage of blurring. This results in a grid of tangents. While the technique used for obtaining the tangents is at present rather ad hoc, researchers plan to adopt an existing framework for extracting edge orientation information over a variety of resolutions, such as suggested by Watson (1987, 1983), Marr and Hildreth (1980), or Canny (1986).
To enable development of a vector wind gust model suitable for orbital flight test operations and trade studies, hypotheses concerning the distributions of gust component variables were verified. Methods for verification of hypotheses that observed gust variables, including gust component magnitude, gust length, u range, and L range, are gamma distributed and presented. Observed gust modulus has been drawn from a bivariate gamma distribution that can be approximated with a Weibull distribution. Zonal and meridional gust components are bivariate gamma distributed. An analytical method for testing for bivariate gamma distributed variables is presented. Two distributions for gust modulus are described and the results of extensive hypothesis testing of one of the distributions are presented. The validity of the gamma distribution for representation of gust component variables is established.
A method is described of rapidly estimating the position of the laminar separation point from the given pressure distribution along a body; the method is applicable to a fairly wide variety of cases. The laminar separation point is found by the von Karman-Millikan method for a series of velocity distributions along a flat plate, which consist of a region of uniform velocity followed by a region of uniform decreased velocity. It is shown that such a velocity distribution can frequently replace the actual velocity distribution along a body insofar as the effects on laminar separation are concerned. An example of the application of the method is given by using it to calculate the position of the laminar separation point on the NACA 0012 airfoil section at zero lift. The agreement between the position of the separation point calculated according to the present method and that found from more elaborate computations is very good.
Two techniques for extending the range of applicability of the basic vortex-lattice method are discussed. The first improves the computation of aerodynamic forces on thin, low-aspect-ratio wings of arbitrary planforms at subsonic Mach numbers by including the effects of leading-edge and tip vortex separation, characteristic of this type wing, through use of the well-known suction-analogy method of E. C. Polhamus. Comparisons with experimental data for a variety of planforms are presented. The second consists of the use of the vortex-lattice method to predict pressure distributions over thick multi-element wings (wings with leading- and trailing-edge devices). A method of laying out the lattice is described which gives accurate pressures on the top and part of the bottom surface of the wing. Limited comparisons between the result predicted by this method, the conventional lattice arrangement method, experimental data, and 2-D potential flow analysis techniques are presented.
Theoretical and numerical results are presented for inverse problems involving estimation of spatially varying parameters such as stiffness and damping in distributed models for elastic structures such as Euler-Bernoulli beams. An outline of algorithms used and a summary of computational experiences are presented.
A model that is currently used to predict the priming and performance limitations of a monogroove heat pipe is expanded to include the boiling limitation and the cross-sectional temperature distribution as determined from a multidimensional finite element analysis technique. The improved model is verified experimentally and shown to accurately predict the cross-sectional temperature distribution when the heat flux distribution is known. The model provides a way to estimate the level at which nucleate boiling and the associated dryout of the capillary wick occurs.
Explore the source record for details and available documents.