Search NASA⌕ Search

SEARCH · Search NASA

Results for “Scaling and squaring method”

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 37 records · Page 2

Calculation of Weibull strength parameters, Batdorf flaw density constants and related statistical quantities using PC-CARES

This manual describes the operation and theory of the PC-CARES (Personal Computer-Ceramic Analysis and Reliability Evaluation of Structures) computer program for the IBM PC and compatibles running PC-DOS/MS-DOR OR IBM/MS-OS/2 (version 1.1 or higher) operating systems. The primary purpose of this code is to estimate Weibull material strength parameters, the Batdorf crack density coefficient, and other related statistical quantities. Included in the manual is the description of the calculation of shape and scale parameters of the two-parameter Weibull distribution using the least-squares analysis and maximum likelihood methods for volume- and surface-flaw-induced fracture in ceramics with complete and censored samples. The methods for detecting outliers and for calculating the Kolmogorov-Smirnov and the Anderson-Darling goodness-of-fit statistics and 90 percent confidence bands about the Weibull line, as well as the techniques for calculating the Batdorf flaw-density constants are also described.

Szatmary, Steven A.↗

Characterization, parameter estimation, and aircraft response statistics of atmospheric turbulence

A nonGaussian three component model of atmospheric turbulence is postulated that accounts for readily observable features of turbulence velocity records, their autocorrelation functions, and their spectra. Methods for computing probability density functions and mean exceedance rates of a generic aircraft response variable are developed using nonGaussian turbulence characterizations readily extracted from velocity recordings. A maximum likelihood method is developed for optimal estimation of the integral scale and intensity of records possessing von Karman transverse of longitudinal spectra. Formulas for the variances of such parameter estimates are developed. The maximum likelihood and least-square approaches are combined to yield a method for estimating the autocorrelation function parameters of a two component model for turbulence.

Mark, W. D.↗

Simulation of ocean SAR images via phase history generation

A method for simulating an ocean synthetic-aperture radar (SAR) image is illustrated for a simple internal wave current pattern. The method calculates both the amplitude image and the radar signal history. The simulation model used consists of six stages. In the first stage, the full wave spectrum in two spatial coordinates is calculated from the wind speed and direction using the action spectral density equation. In the second stage, the pixel size is selected and used as the basis for dividing the spectrum into large and small scale motions. Realizations for the large-scale ocean surface height and velocity are then calculated. In the third stage, the sensor wavelength and geometry are used to calculate the small scale statistics (radar cross section, coherence time, root-mean square (RMS), radial velocity and RMS slope). In the fourth stage, an autoregressive method is used to generate a realization of the surface reflectivity history that is consistent with the radial velocity and the radial velocity variance. In the fifth stage, the signal history is generated by summing the reflectivities at the proper times with the appropriate antenna weighting function for the SAR. As a consequence of this process, speckle is automatically included in the signal. In the final stage, the SAR image can be created by using any of the traditional ways to process the signal history, including variable focusing and multilook processing.

Bennett, John R.↗

Spatial Variability in Column CO2 Inferred from High Resolution GEOS-5 Global Model Simulations: Implications for Remote Sensing and Inversions

Column CO2 observations from current and future remote sensing missions represent a major advancement in our understanding of the carbon cycle and are expected to help constrain source and sink distributions. However, data assimilation and inversion methods are challenged by the difference in scale of models and observations. OCO-2 footprints represent an area of several square kilometers while NASA s future ASCENDS lidar mission is likely to have an even smaller footprint. In contrast, the resolution of models used in global inversions are typically hundreds of kilometers wide and often cover areas that include combinations of land, ocean and coastal areas and areas of significant topographic, land cover, and population density variations. To improve understanding of scales of atmospheric CO2 variability and representativeness of satellite observations, we will present results from a global, 10-km simulation of meteorology and atmospheric CO2 distributions performed using NASA s GEOS-5 general circulation model. This resolution, typical of mesoscale atmospheric models, represents an order of magnitude increase in resolution over typical global simulations of atmospheric composition allowing new insight into small scale CO2 variations across a wide range of surface flux and meteorological conditions. The simulation includes high resolution flux datasets provided by NASA s Carbon Monitoring System Flux Pilot Project at half degree resolution that have been down-scaled to 10-km using remote sensing datasets. Probability distribution functions are calculated over larger areas more typical of global models (100-400 km) to characterize subgrid-scale variability in these models. Particular emphasis is placed on coastal regions and regions containing megacities and fires to evaluate the ability of coarse resolution models to represent these small scale features. Additionally, model output are sampled using averaging kernels characteristic of OCO-2 and ASCENDS measurement concepts to create realistic pseudo-datasets. Pseudo-data are averaged over coarse model grid cell areas to better understand the ability of measurements to characterize CO2 distributions and spatial gradients on both short (daily to weekly) and long (monthly to seasonal) time scales

Ott, L.↗

Metafitting: Weight optimization for least-squares fitting of PTTI data

For precise time intercomparisons between a master frequency standard and a slave time scale, we have found it useful to quantitatively compare different fitting strategies by examining the standard uncertainty in time or average frequency. It is particularly useful when designing procedures which use intermittent intercomparisons, with some parameterized fit used to interpolate or extrapolate from the calibrating intercomparisons. We use the term 'metafitting' for the choices that are made before a fitting procedure is operationally adopted. We present methods for calculating the standard uncertainty for general, weighted least-squares fits and a method for optimizing these weights for a general noise model suitable for many PTTI applications. We present the results of the metafitting of procedures for the use of a regular schedule of (hypothetical) high-accuracy frequency calibration of a maser time scale. We have identified a cumulative series of improvements that give a significant reduction of the expected standard uncertainty, compared to the simplest procedure of resetting the maser synthesizer after each calibration. The metafitting improvements presented include the optimum choice of weights for the calibration runs, optimized over a period of a week or 10 days.

Douglas, Rob J.↗

The magnetic field of Mercury

The paper examines the magnetic field observations and their analyses relating to the determination of the Mercury magnetic field. Methods of analyzing data included: (1) comparison of bow shock and magnetopause relative positions at Mercury to the earth, (2) direct spherical harmonic analysis, (3) magnetosphere modeling by an image dipole, and (4) scaling of a mathematical model for the terrestrial magnetosphere. Dipole moments were determined using partial quadrupole and octupole terms to improve the least-square fit of models to observations; analyses by method (2) yield a convergent series of dipole moments values considered to best represent the intrinsic planetary field. Finally, it is suggested that the origin of the magnetic field of Mercury cannot be uniquely determined, but the sources of convective energy may be radiogenic decay and heat release, gravitational settling, and differentiation of processional torques.

Ness, N. F.↗

Ionospheric gravity wave measurements with the USU dynasonde

A method for the measurement of ionospheric Gravity Wave (GW) using the USU Dynasonde is outlined. This method consists of a series of individual procedures, which includes functions for data acquisition, adaptive scaling, polarization discrimination, interpolation and extrapolation, digital filtering, windowing, spectrum analysis, GW detection, and graphics display. Concepts of system theory are applied to treat the ionosphere as a system. An adaptive ionogram scaling method was developed for automatically extracting ionogram echo traces from noisy raw sounding data. The method uses the well known Least Mean Square (LMS) algorithm to form a stochastic optimal estimate of the echo trace which is then used to control a moving window. The window tracks the echo trace, simultaneously eliminating the noise and interference. Experimental results show that the proposed method functions as designed. Case studies which extract GW from ionosonde measurements were carried out using the techniques described. Geophysically significant events were detected and the resultant processed results are illustrated graphically. This method was also developed for real time implementation in mind.

Berkey, Frank T.↗

Clouding tracing: Visualization of the mixing of fluid elements in convection-diffusion systems

This paper describes a highly interactive method for computer visualization of the basic physical process of dispersion and mixing of fluid elements in convection-diffusion systems. It is based on transforming the vector field from a traditionally Eulerian reference frame into a Lagrangian reference frame. Fluid elements are traced through the vector field for the mean path as well as the statistical dispersion of the fluid elements about the mean position by using added scalar information about the root mean square value of the vector field and its Lagrangian time scale. In this way, clouds of fluid elements are traced and are not just mean paths. We have used this method to visualize the simulation of an industrial incinerator to help identify mechanisms for poor mixing.

Ma, Kwan-Liu↗

Evidence for the Kinematic Sunyaev-Zel'dovich Effect with the Atacama Cosmology Telescope and Velocity Reconstruction from the Baryon Oscillation Spectroscopic Survey

We use microwave temperature maps from two seasons of data from the Atacama Cosmology Telescope at 146 GHz, together with the "Constant Mass" CMASS galaxy sample from the Baryon Oscillation Spectroscopic Survey to measure the kinematic Sunyaev-Zel'dovich (kSZ) effect over the redshift range z1/4 0.4-0.7. We use galaxy positions and the continuity equation to obtain a reconstruction of the line-of-sight velocity field. We stack the microwave temperature at the location of each halo, weighted by the corresponding reconstructed velocity. We vary the size of the aperture photometry filter used, thus probing the free electron profile of these halos from within the virial radius out to three virial radii, on the scales relevant for investigating the missing baryons problem. The resulting best fit kSZ model is preferred over the no-kSZ hypothesis at 3.3 and 2.9 sigma for two independent velocity reconstruction methods, using 25,537 galaxies over 660 square degrees. The data suggest that the baryon profile is shallower than the dark matter in the inner regions of the halos probed here, potentially due to energy injection from active galactic nucleus or supernovae. Thus, by constraining the gas profile on a wide range of scales, this technique will be useful for understanding the role of feedback in galaxy groups and clusters. The effect of foregrounds that are uncorrelated with the galaxy velocities is expected to be well below our signal, and residual thermal Sunyaev-Zel'dovich contamination is controlled by masking the most massive clusters. Finally, we discuss the systematics involved in converting our measurement of the kSZ amplitude into the mean free electron fraction of the halos in our sample.

Schaan, Emmanuel S.↗

Subpixel resolution from multiple images

Multiple images taken from similar locations and under similar lighting conditions contain similar, but not identical, information. Slight differences in instrument orientation and position produces mismatches between the projected pixel grids. These mismatches ensure that any point on the ground is sampled differently in each image. If all the images can be registered with respect to each other to a small fraction of a pixel accuracy, then the information from the multiple images can be combined to increase linear resolution by roughly the square root of the number of images. In addition, the gray-scale resolution of the composite image is also improved. We describe methods for multiple image registration and combination, and discuss some of the problems encountered in developing and extending them. We display test results with 8:1 resolution enhancement, and Viking Orbiter imagery with 2:1 and 4:1 enhancements.

Cheeseman, Peter↗

On the wavelength dependence of the effects of turbulence on average refraction angles in occultations by planetary atmospheres

The dependence of the effects of planetary atmospheric turbulence on radio or optical wavelength in occultation experiments is discussed, and the analysis of Hubbard and Jokipii (1977) is criticized. It is argued that in deriving a necessary condition for the applicability of their method, Hubbard and Jokipii neglect a factor proportional to the square of the ratio of atmospheric or local Fresnel zone radius and the inner scale of turbulence, and fail to establish sufficient conditions, thereby omitting the square of the ratio of atmospheric scale height and the local Fresnel zone radius. The total discrepancy is said to mean that the results correspond to geometrical optics instead of wave optics, as claimed, thus being inapplicable in a dicussion of wavelength dependence. Calculations based on geometrical optics show that the bias in the average bending angle depends on the wavelength in the same way as does the bias in phase path caused by turbulence in a homogeneous atmosphere. Hubbard and Jokipii comment that the criterion of Haugstad and Eshleman is incorrect and show that there is a large wave optical domain where the results are independent of wavelength.

Haugstad, B. S.↗

Statistics of some atmospheric turbulence records relevant to aircraft response calculations

Methods for characterizing atmospheric turbulence are described. The methods illustrated include maximum likelihood estimation of the integral scale and intensity of records obeying the von Karman transverse power spectral form, constrained least-squares estimation of the parameters of a parametric representation of autocorrelation functions, estimation of the power spectra density of the instantaneous variance of a record with temporally fluctuating variance, and estimation of the probability density functions of various turbulence components. Descriptions of the computer programs used in the computations are given, and a full listing of these programs is included.

Mark, W. D.↗

Computer programs for smoothing and scaling airfoil coordinates

Detailed descriptions are given of the theoretical methods and associated computer codes of a program to smooth and a program to scale arbitrary airfoil coordinates. The smoothing program utilizes both least-squares polynomial and least-squares cubic spline techniques to smooth interatively the second derivatives of the y-axis airfoil coordinates with respect to a transformed x-axis system which unwraps the airfoil and stretches the nose and trailing-edge regions. The corresponding smooth airfoil coordinates are then determined by solving a tridiagonal matrix of simultaneous cubic-spline equations relating the y-axis coordinates and their corresponding second derivatives. A technique for computing the camber and thickness distribution of the smoothed airfoil is also discussed. The scaling program can then be used to scale the thickness distribution generated by the smoothing program to a specific maximum thickness which is then combined with the camber distribution to obtain the final scaled airfoil contour. Computer listings of the smoothing and scaling programs are included.

Morgan, H. L., Jr.↗

Analysis of images obtained by the Phobos 2 spacecraft

The reflection spectrum and density of the Martian satellite, Phobos, suggests that it may be a captured C-type asteroid and hence may contain relatively unaltered organic material from the early history of the solar system. Since there is evidence from isotopic anomalies that meteorites presumably derived from such asteroids contain preserved interstellar material, it is conceivable that evidence for interstellar matter might also be present on Phobos. The present research concerns interpretation of the chemical composition of Phobos' surface and analysis of images of Phobos taken by the FREGAT camera equipment aboard the Phobos 2 Spacecraft, including both disk-integrated and disk-resolved photometry. Because the surface of Phobos is a regolith, spectral information on the surface composition will be convolved with information on the surface texture. A primary goal was to untangle these effects. Unfortunately, the Phobos 1 Spacecraft lost contact with Earth during its transit to Mars, and Phobos 2 likewise lost contact after two months in orbit around Mars. Thus, the hovering phase of very close approach and the landing of instruments on the surface of Phobos did not take place, so no data on the elemental composition of the surface is available. However, 37 images were obtained of the satellite. The results include the first such observations at near infrared wavelengths, for which disk-integrated results show that the corresponding geometric albedo is quite small and very similar to the value in the blue. The resolved surface photometry was analyzed by developing new methods of radiative transfer in rough surface layers and a procedure referred to as statistical photoclinometry. The results provide root mean square surface slopes on scales from the order of 1 millimeter to 250 meters, albedo fluctuations, and corresponding correlation lengths.

Irvine, W. M.↗

Least-squares finite element solution of 3D incompressible Navier-Stokes problems

Although significant progress has been made in the finite element solution of incompressible viscous flow problems. Development of more efficient methods is still needed before large-scale computation of 3D problems becomes feasible. This paper presents such a development. The most popular finite element method for the solution of incompressible Navier-Stokes equations is the classic Galerkin mixed method based on the velocity-pressure formulation. The mixed method requires the use of different elements to interpolate the velocity and the pressure in order to satisfy the Ladyzhenskaya-Babuska-Brezzi (LBB) condition for the existence of the solution. On the other hand, due to the lack of symmetry and positive definiteness of the linear equations arising from the mixed method, iterative methods for the solution of linear systems have been hard to come by. Therefore, direct Gaussian elimination has been considered the only viable method for solving the systems. But, for three-dimensional problems, the computer resources required by a direct method become prohibitively large. In order to overcome these difficulties, a least-squares finite element method (LSFEM) has been developed. This method is based on the first-order velocity-pressure-vorticity formulation. In this paper the LSFEM is extended for the solution of three-dimensional incompressible Navier-Stokes equations written in the following first-order quasi-linear velocity-pressure-vorticity formulation.

Jiang, Bo-Nan↗

The thermodynamic properties of nitrogen from 65 to 2000 K with pressures to 10,000 atmospheres

An equation of state is presented for liquid and gaseous nitrogen for temperatures from 65 degrees K to 2000 degrees K and pressures to 10,000 atmospheres. All the pressure-density-temperature data available from published literature have been corrected and applied to bring experimental temperatures into accord with the International Practical Temperature Scale of 1968. The coefficients of the equation of state were determined by a weighted least squares fit to selected published pressure-density-temperature data. The methods of weighting the various data for simultaneous fitting are presented and discussed.

Jacobsen, R. T.↗

Aspects of the determination of winds by means of scatterometry and of the utilization of vector wind data for meteorological forecasts

The present paper provides a description of four aspects of scatterometer winds and their uses. The theory of wave generation by the wind is considered along with an analysis of the properties of superobservations, and studies of intermittent versus continuous data assimilation methods for numerical weather predictions which use remotely sensed data. A comparison of the sum of squares versus the maximum likelihood method for recovering the vector winds is also conducted. Questions regarding wind speed, friction velocity, or normal stress are discussed and synoptic scale fields from Seasat-SASS data are examined.

Pierson, W. J., Jr.↗

Large-eddy simulation of a backward facing step flow using a least-squares spectral element method

We report preliminary results obtained from the large eddy simulation of a backward facing step at a Reynolds number of 5100. The numerical platform is based on a high order Legendre spectral element spatial discretization and a least squares time integration scheme. A non-reflective outflow boundary condition is in place to minimize the effect of downstream influence. Smagorinsky model with Van Driest near wall damping is used for sub-grid scale modeling. Comparisons of mean velocity profiles and wall pressure show good agreement with benchmark data. More studies are needed to evaluate the sensitivity of this method on numerical parameters before it is applied to complex engineering problems.

Chan, Daniel C.↗