Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “error vector”

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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A three-dimensional model of terrain-induced updrafts for movement ecology studies

Abstract Background Spatially explicit simulation models of animal movements through the atmosphere necessarily require a representation of the spatial and temporal variation of atmospheric conditions. In particular, for movements of soaring birds that rely extensively on vertical updrafts to avoid flapping flight, accurate and reliable estimation of the vertical component of wind is critical. The interaction between wind and complex terrain shapes both the horizontal and vertical wind fields, highlighting the need to model the coupling between local terrain features and atmospheric conditions at scales relevant to animal movement. Methods In this work, we propose a new empirical model for estimating the orographic updraft field. The model is developed using computational fluid dynamics simulations of canonical atmospheric conditions over moderately complex terrain. To isolate buoyancy and thermal effects, and focus on terrain-induced effects, we use only simulations of a neutrally stratified atmosphere to develop the model. The model, which we name Engineering Vertical Velocity Estimator (EVVE), is simple to implement and is a function of the underlying terrain elevation map, the desired height above ground level (AGL), and wind conditions at a reference height (80 m). We validate the model with data from the Alaiz mountain (Spain) field campaign. Results Compared to observations, the proposed improved model estimates the updrafts at 120 m AGL with a mean error of 0.11 m/s ( $$\sigma =0.28$$ σ = 0.28 m/s), compared to 0.85 m/s ( $$\sigma =0.58$$ σ = 0.58 m/s) for its baseline. For typical land-based wind turbine hub heights of 80 m AGL, the proposed model has a mean error of 0.04 m/s ( $$\sigma =0.25$$ σ = 0.25 m/s), compared to baseline 0.54 m/s ( $$\sigma =0.45$$ σ = 0.45 m/s) estimations. We illustrate an application of the model in movement ecology by comparing simulated tracks and presence maps of golden eagles ( Aquila chrysaetos ) moving across two distinct landscapes. The tracks and presence maps are obtained using a simple heuristic-based movement model, with the updraft field given by the proposed model and a wind vector-based estimation approach that is currently in wide use in movement ecology studies of raptors and other soaring birds. Conclusions We highlight that movement model results can be sensitive to the underlying orographic updraft model, especially in studies of fine-scale movements in regions of complex topography. We suggest adopting the proposed model rather than the wind vector estimation method for studies of soaring bird movements.

17 WIND ENERGY↗
NASA NTRS1987

Extension and applications of flux-vector splitting to unsteady calculations on dynamic meshes

The Van Leer method of flux-vector splitting for the Euler equations is extended for use on moving meshes and all the properties of the original splittings are maintained. The solution is advanced in time with an implicit, approximately factored algorithm. The use of multiple grids to reduce the computer time is investigated. A substantial reduction in computer time to resolve a pitching cycle is easily obtained with virtually no loss in accuracy. A subiterative procedure to eliminate factorization and linearization errors so that larger time steps can be used is also investigated. Subsequent computations show good agreement with experimental data for transonic and supersonic airfoils and wings undergoing forced pitching oscillation.

Anderson, W. Kyle↗
DOE OSTI2019

The Buzzard Flock: Dark Energy Survey Synthetic Sky Catalogs

We present a suite of 18 synthetic sky catalogs designed to support science analysis of galaxies in the Dark Energy Survey Year 1 (DES Y1) data. For each catalog, we use a computationally efficient empirical approach, ADDGALS, to embed galaxies within light-cone outputs of three dark matter simulations that resolve halos with masses above ~5x10^12 h^-1 m_sun at z <= 0.32 and 10^13 h^-1 m_sun at z~2. The embedding method is tuned to match the observed evolution of galaxy counts at different luminosities as well as the spatial clustering of the galaxy population. Galaxies are lensed by matter along the line of sight --- including magnification, shear, and multiple images --- using CALCLENS, an algorithm that calculates shear with 0.42 arcmin resolution at galaxy positions in the full catalog. The catalogs presented here, each with the same LCDM cosmology (denoted Buzzard), contain on average 820 million galaxies over an area of 1120 square degrees with positions, magnitudes, shapes, photometric errors, and photometric redshift estimates. We show that the weak-lensing shear catalog, redMaGiC galaxy catalogs and redMaPPer cluster catalogs provide plausible realizations of the same catalogs in the DES Y1 data by comparing their magnitude, color and redshift distributions, angular clustering, and mass-observable relations, making them useful for testing analyses that use these samples. We make public the galaxy samples appropriate for the DES Y1 data, as well as the data vectors used for cosmology analyses on these simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Reducing computation in an i-vector speaker recognition system using a tree-structured universal background model

The majority of state-of-the-art speaker recognition systems (SR) utilize speaker models that are derived from an adapted universal background model (UBM) in the form of a Gaussian mixture model (GMM). This is true for GMM supervector systems, joint factor analysis systems, and most recently i-vector systems. In all of the identified systems, the posterior probabilities and sufficient statistics calculations represent a computational bottleneck in both enrollment and testing. We propose a multi-layered hash system, employing a tree-structured GMM–UBM which uses Runnalls’ Gaussian mixture reduction technique, in order to reduce the number of these calculations. Moreover, with this tree-structured hash, we can trade-off reduction in computation with a corresponding degradation of equal error rate (EER). As an example, we also reduce this computation by a factor of 15× while incurring less than 10% relative degradation of EER (or 0.3% absolute EER) when evaluated with NIST 2010 speaker recognition evaluation (SRE) telephone data.

42 ENGINEERING↗
NASA NTRS1987

A vector-based failure detection and isolation algorithm for a dual fail-operational redundant strapdown inertial measurement unit

A vector-based failure detection and isolation technique for a skewed array of two degree-of-freedom inertial sensors is developed. Failure detection is based on comparison of parity equations with a threshold, and isolation is based on comparison of logic variables which are keyed to pass/fail results of the parity test. A multi-level approach to failure detection is used to ensure adequate coverage for the flight control, display, and navigation avionics functions. Sensor error models are introduced to expose the susceptibility of the parity equations to sensor errors and physical separation effects. The algorithm is evaluated in a simulation of a commercial transport operating in a range of light to severe turbulence environments. A bias-jump failure level of 0.2 deg/hr was detected and isolated properly in the light and moderate turbulence environments, but not detected in the extreme turbulence environment. An accelerometer bias-jump failure level of 1.5 milli-g was detected over all turbulence environments. For both types of inertial sensor, hard-over, and null type failures were detected in all environments without incident. The algorithm functioned without false alarm or isolation over all turbulence environments for the runs tested.

Morrell, Frederick R.↗
NASA NTRS1991

Experimental study of trajectory planning and control of a high precision robot manipulator

The kinematic and trajectory planning is presented for a 6 DOF end-effector whose design was based on the Stewart Platform mechanism. The end-effector was used as a testbed for studying robotic assembly of NASA hardware with passive compliance. Vector analysis was employed to derive a closed-form solution for the end-effector inverse kinematic transformation. A computationally efficient numerical solution was obtained for the end-effector forward kinematic transformation using Newton-Raphson method. Three trajectory planning schemes, two for fine motion and one for gross motion, were developed for the end-effector. Experiments conducted to evaluate the performance of the trajectory planning schemes showed excellent tracking quality with minimal errors. Current activities focus on implementing the developed trajectory planning schemes on mating and demating space-rated connectors and using the compliant platform to acquire forces/torques applied on the end-effector during the assembly task.

Nguyen, Charles C.↗
NASA NTRS1996

Evaluation of IGS Orbits with Satellite Laser Ranging

The accuracy with which orbits for the Global Positioning System (GPS) spacecraft, can be computed directly affects the accuracy of the resulting site coordinates and polar motion. Several groups routinely analyze GPS ground tracking data to compute precise orbits and terrestrial reference frame solutions. In this paper, we infer the accuracy of the orbits of two of the GPS satellites by comparing to independent laser ranges of subcentimeter accuracy obtained by a small but reasonably well distributed network of tracking sites. We find that all seven International GPS Service for Geodynamics (IGS) analysis centers achieve range residual root mean square (rms) errors at or below the 100 mm level. The best orbit solutions, from JPL, CODE, and the IGS combined product, yield a residual rms of about 50 mm. These residuals are consistent with three dimensional orbit errors of less than 150 mm. Estimating yaw rates for the spacecraft during shadow events, and using these estimates to compute the laser residual, significantly improves the fit. A small mean residual value of -15 to -30 mm seems to exist for most centers and laser sites which is not fully explained at present, but may be due to uncertainties in the corrections to the laser data, such as the reflector to spacecraft center of mass vector or small reference frame differences between the SLR sites and the GPS orbits.

Watkins, M. M.↗
NASA NTRS2004

In-Flight Assessment of a Pursuit Guidance Display Format for Manually Flown Precision Instrument Approaches

In-flight evaluations of a pursuit guidance display system for manually flown precision instrument approaches were performed. The guidance system was integrated into the RASCAL JUH-60A Black Hawk helicopter. The applicability of the pursuit guidance disp1aFs to the operation of Runway Independent Aircraft (RIA) is made evident because the displays allow the pilot to fly a complex, multi-segment, descending, decelerating approach trajectory. The complex trajectory chosen for this in-flight assessment began from a downwind abeam position at 110 knots and was hand-flown to a 50 ft decision altitude at 40 knots using a rate-command/attitude-hold plus turn-coordination control system. The elements of the pursuit guidance format displayed on a 10-inch liquid crystal display (LCD) flat panel consisted of a flightpath vector and a "leader" aircraft as the pursuit guidance element. Approach guidance was based primarily on carrier-phase differential Global Positioning System (GPS) navigation, and secondarily on both medium accuracy inertial navigation unit states and air data computer states. Required Navigation Performance (RNP) concepts were applied to the construction of display elements such as lateral/vertical deviation indicators and a tunnel that indicated to the pilot, in real-time, the performance with respect to RNP error bounds. The results of the flight evaluations of the guidance display show that precise path control for operating within tight RNP boundaries (RNP 0.007NM/24ft for initial approach, RNP 0.008NM/19ft for intermediate approach, and RNP 0.002NM/9ft for final approach) is attainable with minimal to moderate pilot workload.

Moralez, Ernesto, III↗
NASA NTRS2011

The Extreme Spin of the Black Hole in Cygnus X-1

The compact primary in the X-ray binary Cygnus X-1 was the first black hole to be established via dynamical observations. We have recently determined accurate values for its mass and distance, and for the orbital inclination angle of the binary. Building on these results, which are based on our favored (asynchronous) dynamical model, we have measured the radius of the inner edge of the black hole s accretion disk by fitting its thermal continuum spectrum to a fully relativistic model of a thin accretion disk. Assuming that the spin axis of the black hole is aligned with the orbital angular momentum vector, we have determined that Cygnus X-1 contains a near-extreme Kerr black hole with a spin parameter a* > 0.95 (3(sigma)). For a less probable (synchronous) dynamical model, we find a. > 0.92 (3 ). In our analysis, we include the uncertainties in black hole mass, orbital inclination angle, and distance, and we also include the uncertainty in the calibration of the absolute flux via the Crab. These four sources of uncertainty totally dominate the error budget. The uncertainties introduced by the thin-disk model we employ are particularly small in this case given the extreme spin of the black hole and the disk s low luminosity.

Gou, Lijun↗
NASA NTRS1978

Roll tracking effects of G-vector tilt and various types of motion washout

In a dogfight scenario, the task was to follow the target's roll angle while suppressing gust disturbances. All subjects adopted the same behavioral strategies in following the target while suppressing the gusts, and the MFP-fitted math model response was generally within one data symbol width. The results include the following: (1) comparisons of full roll motion (both with and without the spurious gravity tilt cue) with the static case. These motion cues help suppress disturbances with little net effect on the visual performance. Tilt cues were clearly used by the pilots but gave only small improvement in tracking errors. (2) The optimum washout (in terms of performance close to real world, similar behavioral parameters, significant motion attenuation (60 percent), and acceptable motion fidelity) was the combined attenuation and first-order washout. (3) Various trends in parameters across the motion conditions were apparent, and are discussed with respect to a comprehensive model for predicting adaptation to various roll motion cues.

Jex, H. R.↗
DOE OSTI2011

Adaptive sparse polynomial chaos expansion based on least angle regression

Polynomial chaos (PC) expansions are used in stochastic finite element analysis to represent the random model response by a set of coefficients in a suitable (so-called polynomial chaos) basis. The number of terms to be computed grows dramatically with the size of the input random vector, which makes the computational cost of classical solution schemes (may it be intrusive (i.e. of Galerkin type) or non intrusive) unaffordable when the deterministic finite element model is expensive to evaluate. To address such problems, the paper describes a non intrusive method that builds a sparse PC expansion. First, an original strategy for truncating the PC expansions, based on hyperbolic index sets, is proposed. Then an adaptive algorithm based on least angle regression (LAR) is devised for automatically detecting the significant coefficients of the PC expansion. Beside the sparsity of the basis, the experimental design used at each step of the algorithm is systematically complemented in order to avoid the overfitting phenomenon. The accuracy of the PC metamodel is checked using an estimate inspired by statistical learning theory, namely the corrected leave-one-out error. As a consequence, a rather small number of PC terms are eventually retained (sparse representation), which may be obtained at a reduced computational cost compared to the classical 'full' PC approximation. The convergence of the algorithm is shown on an analytical function. Then the method is illustrated on three stochastic finite element problems. The first model features 10 input random variables, whereas the two others involve an input random field, which is discretized into 38 and 30 - 500 random variables, respectively.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1986

Computing Relative Joint Positions of Robot Arms

Vector-algebra method developed for extracting Denavit-Hartenberg parameters for any assembled robot arm. Method for extracting relative joint geometry of robot arms useful to researchers who need data for existing robot arms for either validation of mathematical models or for studies involving actual control of these devices. Method, does not require robot arm to be disassembled, also useful in recalibration of misalined or bent robot arm and becomes useful industrial procedure. Merit of method is errors not propagated.

Barker, L. K.↗
NASA NTRS1992

Plate tectonics from VLBI and SLR global data

This study is based on data derived from fifteen years of observations of the SLR (side-looking radar) network and six years of the VLBI (very long baseline interferometry) network. In order to use all available information VLBI and SLR global data sets were combined in a least squares fashion to calculate station horizontal velocities. All significant data pertaining to a single site contribute to the station horizontal motion. The only constraint on the solution is that no vertical motion is allowed. This restriction does not greatly affect the precision of the overall solution given the fact that the expected vertical motion for most stations, even those experiencing post glacial uplift, is well under 1 cm/yr. Since the average baseline is under 4,000 km, only a small fraction of the station vertical velocity is translated into baseline rates so that the error introduced in the solution by restricting up-down station movement is minimal. As a reference, station velocities were then compared to the ones predicted by the NUVEL-1 geological model of DeMets et al. (1990). The focus of the study is on analyzing these discrepancies for global plate tectonics as well as regional tectonic settings. The method used also allows us not only to derive horizontal motion for individual stations but also to calculate Euler vectors for those plates that have enough stations located on the stable interior like North America, Pacific, Eurasia, and Australia.

Harrison, Christopher G. A.↗
NASA NTRS2018

Estimation of Lower-Body Kinetics from Loading Profile and Kinematics Alone, Without Measured Ground Reaction Forces

Biomechanical models of human motion can estimate kinetic outcomes, such as joint moments, joint forces and muscle forces. Typically, one performs an inverse dynamics (ID) analysis to compute joint moments from joint angles and measured external forces. Sometimes it is impractical to measure ground reaction forces and moments (GRF&M). We devised an empirical method for performing ID analysis of resistance exercises without measured GRF&M. The method solves the multibody dynamics equations of motion with four key assumptions about the GRF&M that reduce the number of unknowns. The assumptions are 1) negligible ground reaction moments, 2) fixed lateral/medial location of the center of pressure (COP), 3) equal fore/aft location of the COP between the feet, and 4) constant angle of the GRF vector relative to the vertical axis in the frontal plane. We used evaluation trials from a spaceflight countermeasure resistance training device to test this approach. Four participants performed squat and deadlift exercises at various loads. We compared results from traditional ID analysis to results without measured GRF&M using our method. We found that joint moment trajectories in the sagittal plane were qualitatively similar in shape between the two methods, and the amount of root mean squared error (RMSE), measured by difference in joint moment impulse, was typically under 10 percent. Non-sagittal joint moment trajectories, which are much lower in overall magnitude, were not qualitatively similar in shape between the two methods. Non-sagittal moments displayed much higher RMSE, with typical values well over 50 percent. These findings were further supported by validation metrics (Sprague and Geers' P and M metrics, Pearson's r correlation coefficient). Based on these findings, we concluded that useful kinetic results are obtained from ID analysis of squat and deadlift exercises, even when GRF&M are not measured, as long as the outcomes of interest lie in the sagittal plane.

computational modeling↗
DOE OSTI2011

Searches for Physics beyond the Standard Model

The Thomas Jefferson National Accelerator Laboratory has the demonstrated ability to test the fundamental symmetries of nature to very great precision and thereby probe for new physics beyond the Standard Model (SM). In the following objectives and descriptions will be given of three Jefferson Laboratory (JLab) experiments: Qweak, MOLLER, and PVDIS. The Qweak experiment is to measure the weak charge of the proton (via the vector coupling of the Z0 boson to the proton). The MOLLER experiment is to measure the weak charge of the electron. The PV-DIS experiment will measure combinations of vector couplings and axial-vector couplings to the quarks of the nucleons. These three experiments follow from the advances made in precision parity-violating electron scattering measurements at the CEBAF of JLab. The Standard Model makes accurate predictions of the 'running' of the electroweak mixing angle or sin2(θw) from the Z0 pole down to low energies and therefore of the weak charges of the proton and electron. The Qweak experiment will make the first precision determination of the weak charge of the proton, Qpw = 1 – 4 sin2 (θw), from a measurement of the parity-violating asymmetry in the elastic scattering of longitudinally polarized electrons from the protons in a liquid hydrogen target at very low momentum transfer. The projected result will determine the proton's weak charge with a 4.1% total error and consequently sin2 (θw) with a 0.3% error. The Qweak experiment is at present three months into its commissioning run. The MOLLER experiment is to measure the parity-violating asymmetry in the scattering of 11 GeV longitudinally polarized electrons from the atomic electrons in a liquid hydrogen target. The longitudinal analyzing power Az is predicted to be 35.6 ppb at the kinematics of the experiment and is to be determined with a precision of 0.73 ppb, which would make the MOLLER experiment the most precise parity-violation experiment ever undertaken. The result would yield a measurement of the weak charge of the electron to 2.3% at an average Q2 value of 0.0056 (GeV/c)2 and in turn a determination of the electroweak mixing angle sin2(θw) with an uncertainty of ±0.00026 (stat) ±0.00013 (syst), comparable to the accuracy of the two best determinations at the Z0 pole. The PVDIS experiment is to measure the parity-violating asymmetry in deep inelastic scattering of longitudinally polarized electrons from an unpolarized deuterium target. The longitudinal analyzing power can be expressed in terms of the quark distribution functions of the deuterium target and the couplings C1q (axial electron x vector quark) and C2q (vector electron x axial quark), which in the Standard Model can be expressed in terms of sin2(θw).

08 HYDROGEN↗
DOE OSTI2022

Entropic lens on stabilizer states

The n-qubit stabilizer states are those left invariant by a 2 n -element subset of the Pauli group. The Clifford group is the group of unitaries which take stabilizer states to stabilizer states; a physically motivated generating set, the Hadamard, phase, and controlled-not (cnot) gates which comprise the Clifford gates, impose a graph structure on the set of stabilizers. We explicitly construct these structures, the “reachability graphs,” at n ≤ 5. When we consider only a subset of the Clifford gates, the reachability graphs separate into multiple, often complicated, connected components. Seeking an understanding of the entropic structure of the stabilizer states, which is ultimately built up by cnot gate applications on two qubits, we are motivated to consider the restricted subgraphs built from the Hadamard and cnot gates acting on only two of the n qubits. We show how the two subgraphs already present at two qubits are embedded into more complicated subgraphs at three and four qubits. We argue that no additional types of subgraph appear beyond four qubits, but that the entropic structures within the subgraphs can grow progressively more complicated as the qubit number increases. Starting at four qubits, some of the stabilizer states have entropy vectors which are not allowed by holographic entropy inequalities. Here, we comment on the nature of the transition between holographic and nonholographic states within the stabilizer reachability graphs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Generating moment matching scenarios using optimization techniques

An optimization based method is proposed to generate moment matching scenarios for numerical integration and its use in stochastic programming. The main advantage of the method is its flexibility: it can generate scenarios matching any prescribed set of moments of the underlying distribution rather than matching all moments up to a certain order, and the distribution can be defined over an arbitrary set. This allows for a reduction in the number of scenarios and allows the scenarios to be better tailored to the problem at hand. The method is based on a semi-infinite linear programming formulation of the problem that is shown to be solvable with polynomial iteration complexity. A practical column generation method is implemented. The column generation subproblems are polynomial optimization problems; however, they need not be solved to optimality. It is found that the columns in the column generation approach can be efficiently generated by random sampling. The number of scenarios generated matches a lower bound of Tchakaloff's. The rate of convergence of the approximation error is established for continuous integrands, and an improved bound is given for smooth integrands. Extensive numerical experiments are presented in which variants of the proposed method are compared to Monte Carlo and quasi-Monte Carlo methods on both numerical integration problems and stochastic optimization problems. The benefits of being able to match any prescribed set of moments, rather than all moments up to a certain order, is also demonstrated using optimization problems with 100-dimensional random vectors. Here, empirical results show that the proposed approach outperforms Monte Carlo and quasi-Monte Carlo based approaches on the tested problems.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1971

Computer simulation of supersonic rarefied gas flow in the transition region, about a spherical probe; a Monte Carlo approach with application to rocket-borne ion probe experiments

This report describes a Monte Carlo simulation of transition flow around a sphere. Conditions for the simulation correspond to neutral monatomic molecules at two altitudes (70 and 75 km) in the D region of the ionosphere. Results are presented in the form of density contours, velocity vector plots and density, velocity and temperature profiles for the two altitudes. Contours and density profiles are related to independent Monte Carlo and experimental studies, and drag coefficients are calculated and compared with available experimental data. The small computer used is a PDP-15 with 16 K of core, and a typical run for 75 km requires five iterations, each taking five hours. The results are recorded on DECTAPE to be printed when required, and the program provides error estimates for any flow field parameter.

Horton, B. E.↗