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 55 records · Page 3

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.↗

Convective turbulence with rotation and magnetic fields

A method for calculating mean square velocity fluctuations, mean square temperature fluctuations, and convective flux for a turbulent rotating fluid with externally applied magnetic field is presented. A new spectral model of large scale turbulence is used which requires, as the sole ingredient, the growth rate of the instability generating the turbulence. Results are presented for the convective flux with rotation and magnetic field for a range of parameters of astrophysical interest. This new formula presented here can be viewed as an extension of the mixing length theory to include magnetic fields and rotation.

Canuto, V. M.↗

AVIRIS study of Death Valley evaporite deposits using least-squares band-fitting methods

Minerals found in playa evaporite deposits reflect the chemically diverse origins of ground waters in arid regions. Recently, it was discovered that many playa minerals exhibit diagnostic visible and near-infrared (0.4-2.5 micron) absorption bands that provide a remote sensing basis for observing important compositional details of desert ground water systems. The study of such systems is relevant to understanding solute acquisition, transport, and fractionation processes that are active in the subsurface. Observations of playa evaporites may also be useful for monitoring the hydrologic response of desert basins to changing climatic conditions on regional and global scales. Ongoing work using Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) data to map evaporite minerals in the Death Valley salt pan is described. The AVIRIS data point to differences in inflow water chemistry in different parts of the Death Valley playa system and have led to the discovery of at least two new North American mineral occurrences. Seven segments of AVIRIS data were acquired over Death Valley on 31 July 1990, and were calibrated to reflectance by using the spectrum of a uniform area of alluvium near the salt pan. The calibrated data were subsequently analyzed by using least-squares spectral band-fitting methods, first described by Clark and others. In the band-fitting procedure, AVIRIS spectra are fit compared over selected wavelength intervals to a series of library reference spectra. Output images showing the degree of fit, band depth, and fit times the band depth are generated for each reference spectrum. The reference spectra used in the study included laboratory data for 35 pure evaporite spectra extracted from the AVIRIS image cube. Additional details of the band-fitting technique are provided by Clark and others elsewhere in this volume.

Crowley, J. K.↗

Determining the Products of Inertia for Small Scale UAVs

Moments of inertia and products of inertia often need to be determined for aircraft. As complex bodies, their mass properties need to be determined experimentally for best accuracy. While several moment of inertia experimental techniques have been developed, there are few to determine the products of inertia. Products of inertia can be easily determined mathematically if the angle between the aircraft x body axis and principal x axis is known. This method finds the principal inclination angle by mathematically correlating the measured moments of inertia about a range of axes of the aircraft. This correlation uses a least squares error minimization of a mathematical model that describes the ellipse of inertia in the aircraft's x-z axes plane. Results from a test conducted on a small scale UAV (Unmanned Aerial Vehicle) at NASA Armstrong Flight Research Center is also presented, which is an example of the intended application of this technique.

mass properties↗

Tests of the NACA 0025 and 0035 Airfoils in the Full-Scale Wind Tunnel

This report presents the results of an investigation conducted in the NACA full-scale wind tunnel to determine the aerodynamic characteristics of the 6 by 36-foot rectangular NACA 0025 and 0035 airfoils. The aerodynamic characteristics of the plain airfoils with rounded and square tips were determined by force tests through a complete angle-of-attack range, in addition, the profile drag was determined by the momentum method. The transition points on the airfoils were located by boundary-layer determinations with small total-head and static tubes. Each airfoil was also tested with a 0.20c full-span split flap. Tuft surveys were included to show the progressive breakdown of flow with increasing angles of attack. Previously published data from tests of the NACA 0009, 0012, and 0018 airfoils in the full-scale tunnel have been included in the summary curves.

Bullivant, W Kenneth↗

A Preliminary ZEUS Lightning Location Error Analysis Using a Modified Retrieval Theory

The ZEUS long-range VLF arrival time difference lightning detection network now covers both Europe and Africa, and there are plans for further expansion into the western hemisphere. In order to fully optimize and assess ZEUS lightning location retrieval errors and to determine the best placement of future receivers expected to be added to the network, a software package is being developed jointly between the NASA Marshall Space Flight Center (MSFC) and the University of Nevada Las Vegas (UNLV). The software package, called the ZEUS Error Analysis for Lightning (ZEAL), will be used to obtain global scale lightning location retrieval error maps using both a Monte Carlo approach and chi-squared curvature matrix theory. At the core of ZEAL will be an implementation of an Iterative Oblate (IO) lightning location retrieval method recently developed at MSFC. The IO method will be appropriately modified to account for variable wave propagation speed, and the new retrieval results will be compared with the current ZEUS retrieval algorithm to assess potential improvements. In this preliminary ZEAL work effort, we defined 5000 source locations evenly distributed across the Earth. We then used the existing (as well as potential future ZEUS sites) to simulate arrival time data between source and ZEUS site. A total of 100 sources were considered at each of the 5000 locations, and timing errors were selected from a normal distribution having a mean of 0 seconds and a standard deviation of 20 microseconds. This simulated "noisy" dataset was analyzed using the IO algorithm to estimate source locations. The exact locations were compared with the retrieved locations, and the results are summarized via several color-coded "error maps."

Elander, Valjean↗

Steady induction effects in geomagnetism. Part 1B: Geomagnetic estimation of steady surficial core motions: A non-linear inverse problem

The problem of estimating a steady fluid velocity field near the top of Earth's core which induces the secular variation (SV) indicated by models of the observed geomagnetic field is examined in the source-free mantle/frozen-flux core (SFI/VFFC) approximation. This inverse problem is non-linear because solutions of the forward problem are deterministically chaotic. The SFM/FFC approximation is inexact, and neither the models nor the observations they represent are either complete or perfect. A method is developed for solving the non-linear inverse motional induction problem posed by the hypothesis of (piecewise, statistically) steady core surface flow and the supposition of a complete initial geomagnetic condition. The method features iterative solution of the weighted, linearized least-squares problem and admits optional biases favoring surficially geostrophic flow and/or spatially simple flow. Two types of weights are advanced radial field weights for fitting the evolution of the broad-scale portion of the radial field component near Earth's surface implied by the models, and generalized weights for fitting the evolution of the broad-scale portion of the scalar potential specified by the models.

Voorhies, Coerte V.↗

Backus-Gilbert theory and its application to retrieval of ozone and temperature profiles

The inversion method provides a quantitative evaluation of the trade-off between vertical resolution of a retrieved profile and formal root-mean-square (rms) error due to measurement noise propagation. The problem of retrieving the top-side ozone profile from backscattered ultraviolet (BUV) measurements is considered. For measurements of the type currently being obtained with the Nimbus 4 and AE-E BUV experiments, it is found that a vertical resolution of approximately 0.75 scale height can be achieved for a formal volume mixing ratio profile error of 10%. Other examples include treatments of the retrieval of temperature profiles from measurements in the 15 micron CO2 absorption band for both the terrestrial and Martian atmospheres. Finally, the method is applied to the problem of retrieving temperature profiles of the Jovian planets from measurements in the far infrared pressure induced H2 lines to be obtained from the Mariner Jupiter/Saturn fly-by missions.

Conrath, B. J.↗

Study of the energy spectrum of primary cosmic rays: EAS size fluctuations at a fixed primary energy

During the initial period of the Samarkand EAS array operations the showers were selected on the basis of charged-particle flux density, and during the subsequent periods the showers were selected on the basis of Cerenkov light flux density. This procedure made it possible to measure the shower energy, to estimate the EAS size fluctuations at a fixed primary energy, and to experimentally obtain the scaling factor K(Ne, Eo) from the EAS size spectrum to the primary energy spectrum. Six scintillators of area S = 2 sq m each were added to the array. The fluctuations of EAS sizes in the showers of fixed primary energies and the scaling factors K(Ne, Eo) were inferred from the data obtained. The showers with zenith angles 30 deg were selected. The EAS axis positions were inferred from the amplitude data of the scintillators. The primary energy Eo was determined by the method of least squares for the known EAS axis position using the data of the Cerenkov detector located at 80 to 150 m EAS axis. It is shown that the Cerenkov light flux fluctuations at 100 m from EAS axis, q sub 100, do not exceed 10% at a fixed EAS energy, so the parameter q sub 100 may be used to estimate the EAS-generating primary particle-energy.

Allev, N.↗

Linear regression in astronomy. I

Five methods for obtaining linear regression fits to bivariate data with unknown or insignificant measurement errors are discussed: ordinary least-squares (OLS) regression of Y on X, OLS regression of X on Y, the bisector of the two OLS lines, orthogonal regression, and 'reduced major-axis' regression. These methods have been used by various researchers in observational astronomy, most importantly in cosmic distance scale applications. Formulas for calculating the slope and intercept coefficients and their uncertainties are given for all the methods, including a new general form of the OLS variance estimates. The accuracy of the formulas was confirmed using numerical simulations. The applicability of the procedures is discussed with respect to their mathematical properties, the nature of the astronomical data under consideration, and the scientific purpose of the regression. It is found that, for problems needing symmetrical treatment of the variables, the OLS bisector performs significantly better than orthogonal or reduced major-axis regression.

Isobe, Takashi↗

Estimating Root Mean Square Errors in Remotely Sensed Soil Moisture over Continental Scale Domains

Root Mean Square Errors (RMSE) in the soil moisture anomaly time series obtained from the Advanced Scatterometer (ASCAT) and the Advanced Microwave Scanning Radiometer (AMSR-E; using the Land Parameter Retrieval Model) are estimated over a continental scale domain centered on North America, using two methods: triple colocation (RMSETC ) and error propagation through the soil moisture retrieval models (RMSEEP ). In the absence of an established consensus for the climatology of soil moisture over large domains, presenting a RMSE in soil moisture units requires that it be specified relative to a selected reference data set. To avoid the complications that arise from the use of a reference, the RMSE is presented as a fraction of the time series standard deviation (fRMSE). For both sensors, the fRMSETC and fRMSEEP show similar spatial patterns of relatively highlow errors, and the mean fRMSE for each land cover class is consistent with expectations. Triple colocation is also shown to be surprisingly robust to representativity differences between the soil moisture data sets used, and it is believed to accurately estimate the fRMSE in the remotely sensed soil moisture anomaly time series. Comparing the ASCAT and AMSR-E fRMSETC shows that both data sets have very similar accuracy across a range of land cover classes, although the AMSR-E accuracy is more directly related to vegetation cover. In general, both data sets have good skill up to moderate vegetation conditions.

Remotely sensed soil moisture validation↗

Accuracy and stability of finite element schemes for the duct transmission problem

An investigation is conducted regarding the feasibility of approaches for improving the efficiency and stability of existing finite element method (FEM) schemes, taking into account both analytical and numerical studies. Of the four schemes considered for the 'steady' problem, the Hermitian Galerkin formulation appears to be the most efficient and therefore the most suitable scheme for futher full scale implementation. The Hermitian residual least squares (RLS) scheme although comparable in accuracy for the cases considered exhibits a slight tendency to cumulative errors. The performance of both the Lagrangian element schemes considered compares poorly with that of their Hermitian element counterparts. This is particularly true of the Lagrangian RLS scheme. The presence of internal oscillatory components is an inevitable consequence of all Galerkin schemes irrespective of element type.

Astley, R. J.↗

Implementation and Evaluation of Irrigation Techniques in the Community Land Model

Several previous studies have highlighted the irrigation-induced impacts on the global and regional water cycle, energy budget, and near-surface climate. While land models are widely used to address this question, the implementations of irrigation in these models vary in complexity. Here, we expand the representation of irrigation in Community Land Model to enable six different irrigation methods. We find that using a combination of irrigation methods, including default, sprinkler, flood and paddy techniques performs best as determined by evaluating the simulated irrigation water withdrawals against observations, and therefore select this combination as the new irrigation scheme. Then, the impact of the new irrigation scheme on surface fluxes is evaluated and detected using single-point simulations. Finally, the global and regional irrigation-induced impacts on surface energy and water fluxes are compared using both the original and the new irrigation scheme. The new irrigation scheme substantially reduces the bias and root-mean-square error of simulated irrigation water withdrawal in the USA and other countries, but considerably overestimates withdrawals in Central China. Results of single-point experiments show that different irrigation methods have different effects on surface fluxes, while the magnitudes are small. At the global scale, the new scheme enlarges the irrigation-induced impacts on water and energy variables relative to the original scheme, with varying magnitudes across regions. Overall, our results suggest that this newly developed scheme is a better tool for simulating irrigation-induced impacts on climate, and highlight the added value of incorporating human water management in Earth system models.

Community land models↗

Practical Aspects of Stabilized FEM Discretizations of Nonlinear Conservation Law Systems with Convex Extension

This talk considers simplified finite element discretization techniques for first-order systems of conservation laws equipped with a convex (entropy) extension. Using newly developed techniques in entropy symmetrization theory, simplified forms of the Galerkin least-squares (GLS) and the discontinuous Galerkin (DG) finite element method have been developed and analyzed. The use of symmetrization variables yields numerical schemes which inherit global entropy stability properties of the POE system. Central to the development of the simplified GLS and DG methods is the Degenerative Scaling Theorem which characterizes right symmetrizes of an arbitrary first-order hyperbolic system in terms of scaled eigenvectors of the corresponding flux Jacobean matrices. A constructive proof is provided for the Eigenvalue Scaling Theorem with detailed consideration given to the Euler, Navier-Stokes, and magnetohydrodynamic (MHD) equations. Linear and nonlinear energy stability is proven for the simplified GLS and DG methods. Spatial convergence properties of the simplified GLS and DO methods are numerical evaluated via the computation of Ringleb flow on a sequence of successively refined triangulations. Finally, we consider a posteriori error estimates for the GLS and DG demoralization assuming error functionals related to the integrated lift and drag of a body. Sample calculations in 20 are shown to validate the theory and implementation.

Barth, Timothy↗

Retrieving the Balanced Winds on the Globe as a Generalized Inverse Problem

A generalized inverse technique is applied to retrieve two types of balanced winds that characterize the large-scale dynamics of the atmosphere: rotational winds based upon the linear balance equation, and divergent winds based upon the vorticity budget equation. Both balance equations are singular at or near the equator. The balance equations are transformed in spherical harmonic function space to an under-determined system, for which the scale-weighed least-squares solution consists of a sum of principal and singular components. The principal components represent the response to the source function for the regular eigenmodes, while the singular components are determined by the projection of an independent measurement on the singular eigenmodes. The method was tested with the NCEP/NCAR reanalysis data in which a quasi-balance condition exists. A realistic balanced wind field is retrievable when the singular components are computed based upon the reanalyzed wind data.

Lu, Huei-Iin↗

Reduction, analysis, and properties of electric current systems in solar active regions

The specific attraction and, in large part, the significance of solar magnetograms lie in the fact that they give the most important data on the electric currents and the nonpotentiality of active regions. Using the vector magnetograms from the Marshall Space Flight Center (MSFC), we employ a unique technique in the area of data analysis for resolving the 180 deg ambiguity in order to calculate the spatial structure of the vertical electric current density. The 180 deg ambiguity is resolved by applying concepts from the nonlinear multivariable optimization theory. The technique is shown to be of particular importance in very nonpotential active regions. The characterization of the vertical electric current density for a set of vector magnetograms using this method then gives the spatial scale, locations, and magnitude of these current systems. The method, which employs an intermediate parametric function which covers the magnetogram and which defines the local `preferred' direction, minimizes a specific functional of the observed transverse magnetic field. The specific functional that is successful is the integral of the square of the vertical current density. We find that the vertical electric current densities have common characteristics for the extended bipolar (beta) (gamma) (delta)-regions studied. The largest current systems have j(sub z)'s which maximizes around 30 mA/sq m and have a linear decreasing distribution to a diameter of 30 Mn.

Gary, G. Allen↗

Reduction, Analysis, and Properties of Electric Current Systems in Solar Active Regions

The specific attraction and, in large part, the significance of solar vector magnetograms lie in the fact that they give the most important data on the electric currents and the nonpotentiality of active regions. Using the vector magnetograms from the Marshall Space Flight Center (MSFC), we employ a unique technique in the area of data analysis for resolving the 180 degree ambiguity in order to calculate the spatial structure of the vertical electric current density. The 180 degree ambiguity is resolved by applying concepts from the nonlinear multivariable optimization theory. The technique is shown to be of particular importance in very nonpotential active regions. The characterization of the vertical electric current density for a set of vector magnetograms using this method then gives the spatial scale, locations, and magnitude of these current systems. The method, which employs an intermediate parametric function which covers the magnetogram and which defines the local "preferred" direction, minimizes a specific functional of the observed transverse magnetic field. The specific functional that is successful is the integral of the square of the vertical current density. We find that the vertical electric current densities have common characteristics for the extended bipolar beta gamma delta-regions studied. The largest current systems have j(sub z)'s which maximizes around 30 mA per square meter and have a linear decreasing distribution to a diameter of 30 Mm.

Gary, G. Allen↗

Multilevel methods for transport equations in diffusive regimes

We consider the numerical solution of the single-group, steady state, isotropic transport equation. An analysis by means of the moment equations shows that a discrete ordinate S(sub N) discretization in direction (angle) with a least squares finite element discretization in space does not behave properly in the diffusion limit. A scaling of the S(sub N) equations is introduced so that the least squares discretization has the correct diffusion limit. For the resulting discrete system a full multigrid algorithm was developed.

Manteuffel, Thomas A.↗