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 19 records

Data Assimilation in the Ocean and in the Atmosphere: What Should be Next?

As part of the International Symposium on Assimilation of Observation in Meteorology and Oceanography, a panel discussion was held on the evening of 15 March 1995. The purpose of this panel discussion was focus on what the next major areas of research in data assimilation should be. The panelists had five minutes each for short presentations (Kalman filters, representers, etc.) and this was followed by an open discussion. This preprocessing will require a good understanding of the fine-scale phenomena. Least square methods such as Kalman filters and variational schemes are inefficient estimators of non-Gaussian field such as chemical traces. Regardless of the modeling technique employed (Lagrangian methods seem best), a least squares assimilation scheme will smear fine structure. The estimator of maximum likelihood must be sought, by examination of tracer probability distributions.

Kalnay, Eugenia↗

Large-scale computation of incompressible viscous flow by least-squares finite element method

The least-squares finite element method (LSFEM) based on the velocity-pressure-vorticity formulation is applied to large-scale/three-dimensional steady incompressible Navier-Stokes problems. This method can accommodate equal-order interpolations and results in symmetric, positive definite algebraic system which can be solved effectively by simple iterative methods. The first-order velocity-Bernoulli function-vorticity formulation for incompressible viscous flows is also tested. For three-dimensional cases, an additional compatibility equation, i.e., the divergence of the vorticity vector should be zero, is included to make the first-order system elliptic. The simple substitution of the Newton's method is employed to linearize the partial differential equations, the LSFEM is used to obtain discretized equations, and the system of algebraic equations is solved using the Jacobi preconditioned conjugate gradient method which avoids formation of either element or global matrices (matrix-free) to achieve high efficiency. To show the validity of this scheme for large-scale computation, we give numerical results for 2D driven cavity problem at Re = 10000 with 408 x 400 bilinear elements. The flow in a 3D cavity is calculated at Re = 100, 400, and 1,000 with 50 x 50 x 50 trilinear elements. The Taylor-Goertler-like vortices are observed for Re = 1,000.

Jiang, Bo-Nan↗

Development of Jet Noise Power Spectral Laws Using SHJAR Data

High quality jet noise spectral data measured at the Aeroacoustic Propulsion Laboratory at the NASA Glenn Research Center is used to examine a number of jet noise scaling laws. Configurations considered in the present study consist of convergent and convergent-divergent axisymmetric nozzles. Following the work of Viswanathan, velocity power factors are estimated using a least squares fit on spectral power density as a function of jet temperature and observer angle. The regression parameters are scrutinized for their uncertainty within the desired confidence margins. As an immediate application of the velocity power laws, spectral density in supersonic jets are decomposed into their respective components attributed to the jet mixing noise and broadband shock associated noise. Subsequent application of the least squares method on the shock power intensity shows that the latter also scales with some power of the shock parameter. A modified shock parameter is defined in order to reduce the dependency of the regression factors on the nozzle design point within the uncertainty margins of the least squares method.

Khavaran, Abbas↗

Computational Issues in Damping Identification for Large Scale Problems

Two damping identification methods are tested for efficiency in large-scale applications. One is an iterative routine, and the other a least squares method. Numerical simulations have been performed on multiple degree-of-freedom models to test the effectiveness of the algorithm and the usefulness of parallel computation for the problems. High Performance Fortran is used to parallelize the algorithm. Tests were performed using the IBM-SP2 at NASA Ames Research Center. The least squares method tested incurs high communication costs, which reduces the benefit of high performance computing. This method's memory requirement grows at a very rapid rate meaning that larger problems can quickly exceed available computer memory. The iterative method's memory requirement grows at a much slower pace and is able to handle problems with 500+ degrees of freedom on a single processor. This method benefits from parallelization, and significant speedup can he seen for problems of 100+ degrees-of-freedom.

Pilkey, Deborah L.↗

Body frame close coupling wave packet approach to gas phase atom-rigid rotor inelastic collisions

The close coupling wave packet (CCWP) method is formulated in a body-fixed representation for atom-rigid rotor inelastic scattering. For J greater than j-max (where J is the total angular momentum and j is the rotational quantum number), the computational cost of propagating the coupled channel wave packets in the body frame is shown to scale approximately as N exp 3/2, where N is the total number of channels. For large numbers of channels, this will be much more efficient than the space frame CCWP method previously developed which scales approximately as N-squared under the same conditions.

Sun, Y.↗

Signal-to-noise enhancement in ground-based intensity observations of solar p modes

Intensity observations of solar p modes are needed to form a complete picture of wave propagation in the photosphere. Ground-based intensity observations are severely hampered by terrestrial atmospheric noise. Partial cancellation of the noise power can be achieved if two spectra having disparate signal-to- noise ratios, and based on time series acquired simultaneously at the same site, are combined. A method of combining the spectra is suggested in which one amplitude is scaled and subtracted from the other. The result is squared yielding a positive-definite power density. To test the method, the intensity of light scattered by the Earth's atmnosphere was recorded at fifteen- second intervals in two narrow bands centered on 0.5 microns and 1.6 microns. When the two resulting spectra were combined, the noise power was attenuated by a factor of 2.7. The scale factor was varied about its optimum value, revealing that noise peaks have a different siganture than signal peaks, and opening up the possibility of a new tool in discrimination against noise peaks. Maxima at symmetry-allowed frequencies and minima at symmetry- forbidden frequencies indicate that the possibility that these results are obtained by chance is only 6.1 x 10(exp -4). The positions of these maxima and minima also support the solar-cycle dependent frequency shifts found by Palle, Regulo, and Roca Cortes.

Germain, Marvin E.↗

Astrometric observations of the satellites of the outer planets. I - The Galilean satellites in 1977

Astrometric observations of the Galilean satellites of Jupiter performed around the opposition of 1977 with a 67-cm photovisual refractor are reported. The measurements are reduced to absolute and intersatellite positions by two different methods: a linear least-squares plate-constant solution to the AGK3 reference-star frame and the trail-scale method of Pascu (1977). Results of a plate-constant reduction are also presented for two observations of JV (Amalthea). Comparison of the data with Liske's (1978) theoretical predictions for the intersatellite positions indicates no systematic bias in the data; a probable scatter of about 15 arcsec in both right ascension and declination is estimated.

Ianna, P. A.↗

A renormalization group analysis of two-dimensional magnetohydrodynamic turbulence

The renormalization group (RNG) method is used to study the physics of two-dimensional (2D) magnetohydrodynamic (MHD) turbulence. It is shown that, for a turbulent magnetofluid in two dimensions, no RNG transformation fixed point exists on account of the coexistence of energy transfer to small scales and mean-square magnetic flux transfer to large scales. The absence of a fixed point renders the RNG method incapable of describing the 2D MHD system. A similar conclusion is reached for 2D hydrodynamics, where enstrophy flows to small scales and energy to large scales. These analyses suggest that the applicability of the RNG method to turbulent systems is intrinsically limited, especially in the case of systems with dual-direction transfer.

Liang, Wenli Z.↗

Rapid estimation of frequency response functions by close-range photogrammetry

The accuracy of a rapid method which estimates the frequency response function from stereoscopic dynamic data is computed. It is shown that reversal of the order of the operations of coordinate transformation and Fourier transformation, which provides a significant increase in computational speed, introduces error. A portion of the error, proportional to the perturbation components normal to the camera focal planes, cannot be eliminated. The remaining error may be eliminated by proper scaling of frequency data prior to coordinate transformation. Methods are developed for least squares estimation of the full 3x3 frequency response matrix for a three dimensional structure.

Tripp, J. S.↗

Calculation of Weibull strength parameters and Batdorf flow-density constants for volume- and surface-flaw-induced fracture in ceramics

The calculation of shape and scale parameters of the two-parameter Weibull distribution is described using the least-squares analysis and maximum likelihood methods for volume- and surface-flaw-induced fracture in ceramics with complete and censored samples. Detailed procedures are given for evaluating 90 percent confidence intervals for maximum likelihood estimates of shape and scale parameters, the unbiased estimates of the shape parameters, and the Weibull mean values and corresponding standard deviations. Furthermore, the necessary steps are described 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 distribution. It also shows how to calculate the Batdorf flaw-density constants by uing the Weibull distribution statistical parameters. The techniques described were verified with several example problems, from the open literature, and were coded. The techniques described were verified with several example problems from the open literature, and were coded in the Structural Ceramics Analysis and Reliability Evaluation (SCARE) design program.

Pai, Shantaram S.↗

Calculation of Weibull strength parameters and Batdorf flow-density constants for volume- and surface-flaw-induced fracture in ceramics

The calculation of shape and scale parametes of the two-parameter Weibull distribution is described using the least-squares analysis and maximum likelihood methods for volume- and surface-flaw-induced fracture in ceramics with complete and censored samples. Detailed procedures are given for evaluating 90 percent confidence intervals for maximum likelihood estimates of shape and scale parameters, the unbiased estimates of the shape parameters, and the Weibull mean values and corresponding standard deviations. Furthermore, the necessary steps are described 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 distribution. It also shows how to calculate the Batdorf flaw-density constants by using the Weibull distribution statistical parameters. The techniques described were verified with several example problems, from the open literature, and were coded in the Structural Ceramics Analysis and Reliability Evaluation (SCARE) design program.

Shantaram, S. Pai↗

Structural Requirements and Scaling Analysis of a Fluidic Mirror Space Telescope Support Structure

The NASA FLUTE project proposes large-scale (50m) fluidic telescopes for astronomy applications. To continue to explore the universe, astronomers require larger and larger telescope apertures. The highest priority astrophysics targets such as exoplanets and early galaxies are extremely faint, motivating larger telescope apertures. However, mission costs depend on aperture diameter, and scaling apertures beyond 10-m apertures faces economic and technological viability challenges. An unsegmented primary mirror made in space via fluidic microgravity shaping would provide a scalable and cost-effective method to scale apertures to 50-m scale while achieving sub-nanometer (root mean square) surface quality. Such microgravity fluidic shaping has been demonstrated in laboratory neutral buoyance environments, parabolic microgravity experiments, as well aboard the International Space Station. One of the main components of a fluidic observatory is the mirror frame. The frame must provide a stable bounding circular ring which the edges of the fluid mirror surface can wet. The frame can optionally provide a ‘floor’ surface on the interior of the ring to provide additional fluid support and reduce required fluid volume. In this work, we evaluate several classes of structural frame architectures potentially suitable for a fluidic telescope support structure. We start by estimating stability requirements, orbital, station keeping, and slew loads based on a notional CONOPS. The scaling of overall fluid mass required for each architecture is evaluated. Preliminary results elucidate the importance of a support floor for overall mission viability above 10-m diameter. We then investigate the scaling of a tetrahedral truss frame support structure. We show that segmented solid shell support surfaces can provide sufficient stability at modest mass fractions. We estimate that the total fluid and frame mass for a 50-m telescope could be on the order of 15,000 kg. Finally, implementation considerations are discussed, including deployment/assembly methodologies. The results of this study establish feasibility of a large-scale fluidic telescope and will guide further architecture development and detailed structural design.

Christine Gregg↗

Steady surficial core motions - An alternate method

A new method for deriving steady surficial core motions from geomagnetic field models is presented and applied. The original method determines the steady velocity field at the top of a frozen-flux core which best fits, in both the spatial and temporal linear least squares sense, a model of the large scale geomagnetic secular variation (SV) at the base of a source-free mantle (Voorhies, 1986). New method is based upon the same physical assumptions, but fits SV at earth's surface instead of the core-mantle boundary (CMB). The newly derived flow differs somewhat from prior solutions, and shows but slight westward drift; however, other key global properties are virtually unchanged. Over a 15 year interval the new method provides a better fit to the time-varying geomagnetic field at earth's surface than does the old.

Voorhies, Coerte V.↗

Estimating the Properties of Hard X-Ray Solar Flares by Constraining Model Parameters

We wish to better constrain the properties of solar flares by exploring how parameterized models of solar flares interact with uncertainty estimation methods. We compare four different methods of calculating uncertainty estimates in fitting parameterized models to Ramaty High Energy Solar Spectroscopic Imager X-ray spectra, considering only statistical sources of error. Three of the four methods are based on estimating the scale-size of the minimum in a hypersurface formed by the weighted sum of the squares of the differences between the model fit and the data as a function of the fit parameters, and are implemented as commonly practiced. The fourth method is also based on the difference between the data and the model, but instead uses Bayesian data analysis and Markov chain Monte Carlo (MCMC) techniques to calculate an uncertainty estimate. Two flare spectra are modeled: one from the Geostationary Operational Environmental Satellite X1.3 class flare of 2005 January 19, and the other from the X4.8 flare of 2002 July 23.We find that the four methods give approximately the same uncertainty estimates for the 2005 January 19 spectral fit parameters, but lead to very different uncertainty estimates for the 2002 July 23 spectral fit. This is because each method implements different analyses of the hypersurface, yielding method-dependent results that can differ greatly depending on the shape of the hypersurface. The hypersurface arising from the 2005 January 19 analysis is consistent with a normal distribution; therefore, the assumptions behind the three non- Bayesian uncertainty estimation methods are satisfied and similar estimates are found. The 2002 July 23 analysis shows that the hypersurface is not consistent with a normal distribution, indicating that the assumptions behind the three non-Bayesian uncertainty estimation methods are not satisfied, leading to differing estimates of the uncertainty. We find that the shape of the hypersurface is crucial in understanding the output from each uncertainty estimation technique, and that a crucial factor determining the shape of hypersurface is the location of the low-energy cutoff relative to energies where the thermal emission dominates. The Bayesian/MCMC approach also allows us to provide detailed information on probable values of the low-energy cutoff, Ec, a crucial parameter in defining the energy content of the flare-accelerated electrons. We show that for the 2002 July 23 flare data, there is a 95% probability that Ec lies below approximately 40 keV, and a 68% probability that it lies in the range 7-36 keV. Further, the low-energy cutoff is more likely to be in the range 25-35 keV than in any other 10 keV wide energy range. The low-energy cutoff for the 2005 January 19 flare is more tightly constrained to 107 +/- 4 keV with 68% probability.

X-rays↗

Techniques for control of long-term reliability of complex integrated circuits. I - Reliability assurance by test vehicle qualification.

Development of an alternate approach to the conventional methods of reliability assurance for large-scale integrated circuits. The product treated is a large-scale T squared L array designed for space applications. The concept used is that of qualification of product by evaluation of the basic processing used in fabricating the product, providing an insight into its potential reliability. Test vehicles are described which enable evaluation of device characteristics, surface condition, and various parameters of the two-level metallization system used. Evaluation of these test vehicles is performed on a lot qualification basis, with the lot consisting of one wafer. Assembled test vehicles are evaluated by high temperature stress at 300 C for short time durations. Stressing at these temperatures provides a rapid method of evaluation and permits a go/no go decision to be made on the wafer lot in a timely fashion.

Van Vonno, N. W.↗

Do large rate coefficients for ion-polar neutral reactions have a serious effect on chemical models of dense clouds?

In order to incorporate large ion-polar neutral rate coefficients into existing gas phase reaction networks, it is necessary to utilize simplified theoretical treatments because of the significant number of rate coefficients needed. The authors have used two simple theoretical treatments: the locked dipole approach of Moran and Hamill for linear polar neutrals and the trajectory scaling approach of Su and Chesnavich for nonlinear polar neutrals. The former approach is suitable for linear species because in the interstellar medium these are rotationally relaxed to a large extent and the incoming charged reactants can lock their dipoles into the lowest energy configuration. The latter approach is a better approximation for nonlinear neutral species, in which rotational relaxation is normally less severe and the incoming charged reactants are not as effective at locking the dipoles. The treatments are in reasonable agreement with more detailed long range theories and predict an inverse square root dependence on kinetic temperature for the rate coefficient. Compared with the locked dipole method, the trajectory scaling approach results in rate coefficients smaller by a factor of approximately 2.5.

Herbst, E.↗

Separation of core noise and jet noise

A method of identification and measurement of core noise and jet noise separately has been developed based on cross-correlation of signals from microphones located at widely separated angles in the far field of a jet. The different coherent properties of core noise and jet noise are used in this method to achieve this separation. Experimental data obtained in a small scale facility is analyzed to demonstrate that this method can be used successfully to separate the mean square pressures of core noise and jet noise.

Parthasarathy, S. P.↗

Silicon Ingot Casting: Heat Exchanger Method. Multi-wire Slicing: Fixed Abrasine Slicing Technique, Phase 3

Ingot casting was scaled up to 16 cm by 16 cm square cross section size and ingots weighing up to 8.1 kg were cast. The high degree of crystallinity was maintained in the large ingot. For large sizes, the nonuniformity of heat treatment causes chipping of the surface of the ingot. Progress was made in the development of a uniform graded structure in the silica crucibles. The high speed slicer blade-head weight was reduced to 37 pounds, allowing surface speeds of up to 500 feet per minute. Slicing of 10 cm diameter workpieces at these speeds increased the through-put of the machine to 0.145 mm/min.

Schmid, F.↗