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 1,099 records · Page 61

Processing laser velocimeter high-speed burst counter data

The statistical calculations of the mean and potential bias errors are considered with respect to their applicability to the randomly sampled velocity data obtained from a laser velocimeter. Techniques for displaying the resulting data in the form of arrow plots, flow streamline plots, and contour maps are also discussed. A description is presented of the operation of a new hardware interface between the signal processing electronics and the computer. The interface has been constructed to make it possible to use recently developed data processing techniques for obtaining turbulent power spectra from laser velocimeter data. Additional capabilities which were built into the interface allow the measurement of the velocity vector for each particle passing through the laser velocimeter sample window. By processing the data with extended versions of the turbulent power spectra techniques, cross spectra, interactions between velocity magnitude fluctuations, and flow angle fluctuations may be investigated.

Meyers, J. F.↗

Seismic Models for Near-Surface Explosion Yield Estimation in Alluvium and Sedimentary Rock

Seismic ground-motion data provide valuable constraints on explosion characteristics, such as yield and height-of-burst/depth-of-burial (HOB/DOB). In this study, we investigate a range of seismic amplitude features and their efficacy in minimizing errors associated with yield estimation. Using a set of explosion recordings from experiments conducted in alluvium and sedimentary rock geologies, we investigated the effectiveness of three different seismic feature types over a range of different frequencies. Using both velocity and displacement data, we investigated the zero-to-peak (ZTP) amplitude of the first arriving P wave, the peak-to-peak (PTP) amplitude of the first arriving P wave, and the P rms amplitude over various time windows starting with the P wave. These three basic features were measured on both vertical component only data and the vertical–radial vector sum. In total, there were 56 different combinations investigated. Our seismic models combine the effects of scaled range and scaled HOB/DOB on the observed seismic features. The results show that the vertical–radial vector sum of the ZTP amplitudes measured on displacement seismograms most consistently produce models with the smallest difference between predicted amplitude values and observed amplitude values when accounting for both alluvium and sedimentary rock lithology. The improvement in fit to the data when incorporating the difference in source lithology is significant. Lastly, if an alluvium lithology is assumed for a sedimentary rock lithology, the difference in yield could be up to ~3.0 times the appropriate yield in the far field.

58 GEOSCIENCES↗

A comprehensive method of estimating electric fields from vector magnetic field and Doppler measurements

Photospheric electric fields, estimated from sequences of vector magnetic field and Doppler measurements, can be used to estimate the flux of magnetic energy (the Poynting flux) into the corona and as time-dependent boundary conditions for dynamic models of the coronal magnetic field. We have modified and extended an existing method to estimate photospheric electric fields that combines a poloidal-toroidal decomposition (PTD) of the evolving magnetic field vector with Doppler and horizontal plasma velocities. Our current, more comprehensive method, which we dub the 'PTD-Doppler-FLCT Ideal' (PDFI) technique, can now incorporate Doppler velocities from non-normal viewing angles. It uses the FISHPACK software package to solve several two-dimensional Poisson equations, a faster and more robust approach than our previous implementations. Here, we describe systematic, quantitative tests of the accuracy and robustness of the PDFI technique using synthetic data from anelastic MHD (ANMHD) simulations, which have been used in similar tests in the past. We find that the PDFI method has less than 1% error in the total Poynting flux and a 10% error in the helicity flux rate at a normal viewing angle (θ = 0) and less than 25% and 10% errors, respectively, at large viewing angles (θ < 60°). We compare our results with other inversion methods at zero viewing angle and find that our method's estimates of the fluxes of magnetic energy and helicity are comparable to or more accurate than other methods. We also discuss the limitations of the PDFI method and its uncertainties.

79 ASTRONOMY AND ASTROPHYSICS↗

Dark energy survey internal consistency tests of the joint cosmological probes analysis with posterior predictive distributions

Beyond ΛCDM, physics or systematic errors may cause subsets of a cosmological data set to appear inconsistent when analysed assuming ΛCDM. Here, we present an application of internal consistency tests to measurements from the Dark Energy Survey Year 1 (DES Y1) joint probes analysis. Our analysis relies on computing the posterior predictive distribution (PPD) for these data under the assumption of ΛCDM. We find that the DES Y1 data have an acceptable goodness of fit to ΛCDM, with a probability of finding a worse fit by random chance of p = 0.046. Using numerical PPD tests, supplemented by graphical checks, we show that most of the data vector appears completely consistent with expectations, although we observe a small tension between large- and small-scale measurements. A small part (roughly 1.5 per cent) of the data vector shows an unusually large departure from expectations; excluding this part of the data has negligible impact on cosmological constraints, but does significantly improve the p-value to 0.10. The methodology developed here will be applied to test the consistency of DES Year 3 joint probes data sets.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Comparison of earth rotation as inferred from radio interferometric, laser ranging and astrometric observations

Results of measurements of the earth's rotation vector for a 400-day period from late September 1980 to December 1981, for which date from VLBI, satellite laser ranging (SLR), and lunar laser ranging (LLR) were available, are compared. The acquisition of the data and their evaluation are described. VLBI, SLR, and classical astrometric determinations of the X-parameter required to describe the location of the rotation pole on the earth's surface are shown, and VLBI, LLR, and classical astrometric determinations of the angle of rotation about this pole (UT1) are presented. The results indicate that VLBI and SLR, at their present stages of development, yield standard errors under 20 cm in the determinations of X, about twofold smaller than obtained from classical measurements, and that VLBI and LLR yield determination of UT1 with standard errors less than 40 cm, somewhat smaller than that of the corresponding determinations from classical observations. Methods for improving these types of intercomparisons are suggested.

Robertson, D. S.↗

Determination of convergence rates across the Ventura Basin, Southern California, using GPS and historical triangulation

Comparison of angles from historical triangulation observations dating as far back as 1932 with Global Positions System (GPS) measurements taken in 1987 indicates that rapid convergence may be taking place on decade timescales in the central and eastern part of the Ventura basin, an east-west trending trough bounded by thrust faults. Changes in angles over this time were analyzed using Prescott's modified Frank's method and in terms of a model which assumes that the regions to the north and south of the basin are rigid blocks undergoing relative motion. For the two block model, inversion of the observed angle changes over the last 28 years for the relative motion vector leads to north-south convergence across the basin of 30 + or - 5 mm/yr, with a left lateral component of 10 + or - 1 mm/yr in the Fillmore-Santa Paula area in the central part of the basin. The modified Frank's method yields strain rates of approximately 2 microrad/yr in both the east and central parts of the basin for measurements spanning the 1971 San Fernando earthquake. Assuming no east-west strain yeilds north-south compression of approximately 3.5 + or - .2 cm/yr. Comparison of triangulation data prior to the earthquake shows no strain outside the margin of error. The convergence rates determined by geodetic techniques are consistent with geologic observations in the area. Such large geodetic deformation rates, with no apparent near-surface creep on the major thrust, can be understood if these faults become subhorizontal at relatively shallow depths and if the subhorizontal portions of the faults are creeping. An alternative explanation of the large displacement rates might be that the pumping of oil in the vicinity of the benchmarks caused large horizontal motions, although it is unlikely that meter scale horizontal motions are due to oil withdrawal. These and other hypotheses are evaluated to better constrain the tectonics of this active region.

Donnellan, Andrea↗

Tele-Autonomous control involving contact

Object localization and its application in tele-autonomous systems are studied. Two object localization algorithms are presented together with the methods of extracting several important types of object features. The first algorithm is based on line-segment to line-segment matching. Line range sensors are used to extract line-segment features from an object. The extracted features are matched to corresponding model features to compute the location of the object. The inputs of the second algorithm are not limited only to the line features. Featured points (point to point matching) and featured unit direction vectors (vector to vector matching) can also be used as the inputs of the algorithm, and there is no upper limit on the number of the features inputed. The algorithm will allow the use of redundant features to find a better solution. The algorithm uses dual number quaternions to represent the position and orientation of an object and uses the least squares optimization method to find an optimal solution for the object's location. The advantage of using this representation is that the method solves for the location estimation by minimizing a single cost function associated with the sum of the orientation and position errors and thus has a better performance on the estimation, both in accuracy and speed, than that of other similar algorithms. The difficulties when the operator is controlling a remote robot to perform manipulation tasks are also discussed. The main problems facing the operator are time delays on the signal transmission and the uncertainties of the remote environment. How object localization techniques can be used together with other techniques such as predictor display and time desynchronization to help to overcome these difficulties are then discussed.

Shao, Lejun↗

Perturbations, TTVs & the (Un)reliability of Ephemerides of Kepler Planets

Fractional uncertainties quoted for orbital periods of the vast majority of planet candidates listed in recent Kepler catalogs are $<10^{-5}$, with values $\sim 10^{-6}$ (corresponding to 2 minutes per 4 years) being typical. These small uncertainties suggest ephemeris predictions for most Kepler planets are robust for decades to come. Tabulated periods list formal uncertainties in the mean times between midpoints of successive transits during the time interval in which transits were observed. However, transit timing variations (TTVs) produce errors in estimates of some planets’ orbital periods that need to be accounted for in certain dynamical investigations and ephemeris predictions. Periodic TTVs with timescales short compared to the interval of Kepler observations largely average out and do not produce significant errors in estimates of orbital periods. TTVs with timescales comparable to the interval of Kepler observations have been fit for dozens of Kepler planet candidates to estimate long-term average orbital periods (Holczer et al. 2016, ApJS 225, 9), and more detailed dynamical models have been used to estimate long-term average periods of a small number of well-studied planets. Most Kepler planets that show large TTVs are near mean-motion orbital resonances with other planets. The largest effect for planets moderately close to two-body resonances are due to rotation of the forced eccentricity vector by resonant perturbations; the timescale for this precession is typically short, so the variations tend to average out. Libration of planets locked in resonances typically occurs on timescales longer than the Kepler baseline, but few Kepler planets appear to be resonantly-locked. The more general but smaller (during the era of Kepler observations) effect is caused by secular precession of the planets' free eccentricities usually has a period much longer than the baseline of the Kepler observations, so it is not accounted for in estimates of mean periods or uncertainties. This precession causes a discrepancy between Kepler era mean orbital period and long-term mean orbital period can exist even for planets having TTVs that are too small to be observable during the epoch of Kepler observations.

Kepler↗

A nine-dimensional ab initio global potential energy surface for the H{sub 2}O{sup +} + H{sub 2} → H{sub 3}O{sup +} + H reaction

An accurate full-dimensional global potential energy surface (PES) is developed for the title reaction. While the long-range interactions in the reactant asymptote are represented by an analytical expression, the interaction region of the PES is fit to more than 81 000 of ab initio points at the UCCSD(T)-F12b/AVTZ level using the permutation invariant polynomial neural network approach. Fully symmetric with respect to permutation of all four hydrogen atoms, the PES provides a faithful representation of the ab initio points, with a root mean square error of 1.8 meV or 15 cm{sup −1}. The reaction path for this exoergic reaction features an attractive and barrierless entrance channel, a submerged saddle point, a shallow H{sub 4}O{sup +} well, and a barrierless exit channel. The rate coefficients for the title reaction and kinetic isotope effect have been determined on this PES using quasi-classical trajectories, and they are in good agreement with available experimental data. It is further shown that the H{sub 2}O{sup +} rotational enhancement of reactivity observed experimentally can be traced to the submerged saddle point. Using our recently proposed Sudden Vector Projection model, we demonstrate that a rotational degree of freedom of the H{sub 2}O{sup +} reactant is strongly coupled with the reaction coordinate at this saddle point, thus unraveling the origin of the pronounced mode specificity in this reaction.

36 MATERIALS SCIENCE↗

Unobtrusive Sensor-Based Occupancy Facing Direction Detection and Tracking Using Advanced Machine Learning Algorithms

Facing direction detection plays a critical role in human-computer interaction and has a wide application in surveillance systems, driving awareness recognition, smart home appliances, computer games, and so on. Current detection methods are mainly focused on extracting specific patterns from user's optical images, which raises concerns on privacy invasion and these detection techniques do not usually work in a dark environment. To address these concerns, here this paper proposes an activity recognition system guided by an unobtrusive sensor (ARGUS). By using a low pixel infrared thermopile array sensor, ARGUS is capable of identifying five facing directions (left 45°/90°, right 45°/90°, and front) through the support vector machine classifier. Also two feature extraction methods are compared. One is manually-defined and the other is based on a pre-trained convolutional neural network (CNN) model. The facing direction detection accuracy resulting from manually-defined features reaches 85.3%, 90.6%, and 85.2% at detection distance of 0.6, 1.2, and 1.8 m, respectively. The level of accuracy resulted from using pre-trained CNN features demonstrates a much more reliable performance (89.1%, 95.3%, and 95.1% at distance of 0.6, 1.2, and 1.8 m, respectively). In addition, ARGUS has been successfully applied for occupancy tracking with a root mean square error of 0.19 m.

42 ENGINEERING↗

Inversion of data from diffraction-limited multiwavelength remote sensors. I - Linear case

The remote sensing inverse problem is considered in which the sensor does not necessarily view the same area on the earth at each wavelength, and spatial correlations of the geophysical parameters may be present. Under the conditions of linearity and stationary statistics, the minimum mean-square error solution to the problem of inverting such data is a spatial filter of the Wiener-Kolmogorov class. The resulting remote-sensing system can be characterized by an impulse response matrix in ordinary space or by a transfer matrix in frequency space. A signal-to-noise matrix for the geophysical parameters to be sensed is also defined; this matrix depends on the postulated a priori statistics and on the characteristics of the remote-sensing instrument. The system transfer matrix and the signal-to-noise matrix are simultaneously diagonalizable. The optimum transfer matrix filters out of the estimate vector those eigenparameters for which eigenvalues of the signal-to-noise matrix are less than unity.

Rosenkranz, P. W.↗

On linear structure and phase rotation invariant properties of block 2(sup l)-PSK modulation codes

Two important structural properties of block 2(l)-ary PSK (phase shift keying) modulation codes, linear structure and phase symmetry, are investigated. For an additive white Gaussian noise (AWGN) channel, the error performance of a modulation code depends on its squared Euclidean distance distribution. Linear structure of a code makes the error performance analysis much easier. Phase symmetry of a code is important in resolving carrier phase ambiguity and ensuring rapid carrier phase resynchronization after temporary loss of synchronization. It is desirable for a code to have as many phase symmetries as possible. A 2(l)-ary modulation code is represented here as a code with symbols from the integer group. S sub 2(l) PSK = (0,1,2,...,2(l)-1), under the modulo-2(l) addition. The linear structure of block 2(l)-ary PSK modulation codes over S sub 2(l)-ary PSK with respect to the modulo-2(l) vector addition is defined, and conditions under which a block 2(l)-ary PSK modulation code is linear are derived. Once the linear structure is developed, phase symmetry of a block 2(l)-ary PSK modulation code is studied. It is a necessary and sufficient condition for a block 2(l)-PSK modulation code, which is linear as a binary code, to be invariant under 180 deg/2(l-h) phase rotation, for 1 is less than or equal to h is less than or equal to l. A list of short 8-PSK and 16-PSK modulation codes is given, together with their linear structure and the smallest phase rotation for which a code is invariant.

Lin, Shu↗

Mass Inferencing Model Creation and Deployment to the RASSOR Lunar Excavation Robot

The Regolith Advanced Surface Systems Operations Robot (RASSOR) Excavator is a teleoperated mobile robotic platform with a unique space regolith excavation capability. The Intelligent Capabilities Enhanced RASSOR research project developed functionality for inferencing regolith mass ingested during RASSOR operation, enhancing RASSOR’s ability to successfully complete ISRU missions. To teleoperate or run autonomously, it is crucial for the quantity of regolith mass ingested by RASSOR to be available as a system state for efficient operation. For example, during autonomous operation, RASSOR should navigate and move to a processing plant to offload the collected regolith when the drums are full; without knowledge of how much mass is in the drums, this type of high-level planning is not possible. Four distinct modeling approaches were employed in developing a mass inferencing approach that could work on RASSOR. All take in system states, such as arm/drum positions, velocities, currents, voltages, and robot pose, and output a mass prediction for each set of the robot’s bucket drums.1) A neural network model that takes a vector of normalized system states; 2) A model that uses the integrated power consumption of an arm-raise (normalized by velocity); 3) A model that uses average drum current over a variable length interval of the drum disengaged from the surface; and 4) A real-time estimation model that aggregates excavation drum current. The developed models run in real time, outputting predictions for the front and rear drums, timestamp of the last prediction, and total mass in RASSOR’s drums. Further testing is required to validate the arm-raise model (2), though initial tests indicate reasonable performance (<10% mean error) on the hardware. The linear fit of average drum-current model (3) had a front value of r^2=0.99 and a rear value of r^2=0.98 on the validation dataset. This model currently has the best performance on unseen data. The real time model (4) is still in development, though initial results on a small subset of the training data show that it has high accuracy in predicting the increase in mass during excavation. Though work remains to be done with deploying a high-fidelity model to the physical system that makes predictions with error below the desired threshold, the modular architecture for model development allows quick adjustment of parameters to increase model fidelity. This architecture can also be adapted to use lunar excavation data to create models that are reflective of RASSOR’s dynamics when operating on the lunar surface. The results are promising as it has been shown that models can be developed that accurately estimate excavated regolith mass.

rassor↗

Three-dimensional traction microscopy accounting for cell-induced matrix degradation

Tractions exerted by cells on the extracellular matrix (ECM) are critical in many important physiological and pathological processes such as embryonic morphogenesis, wound healing, and cancer metastasis. Three-dimensional Traction Microscopy (3DTM) is a tool to quantify cellular tractions by first measuring the displacement field in the ECM in response to these tractions, and then using this measurement to infer tractions. Most applications of 3DTM have assumed that the ECM has spatially-uniform mechanical properties, but cells secrete enzymes that can locally degrade the ECM. In this work, a novel computational method is developed to quantify both cellular tractions and ECM degradation. In particular, the ECM is modeled as a hyperelastic, Neo-Hookean solid, whose material parameters are corrupted by a single degradation parameter. The feasibility of determining both the traction and the degradation parameter is first demonstrated by showing the existence and uniqueness of the solution. An inverse problem is then formulated to determine the nodal values of the traction vector and the degradation parameter, with the objective of minimizing the difference between a predicted and measured displacement field, under the constraint that the predicted displacement field satisfies the equation of equilibrium. The inverse problem is solved by means of a gradient-based optimization approach, and the gradient is computed efficiently using appropriately derived adjoint fields. The computational method is validated in-silico using a geometrically realistic neuronal cell model and synthetic traction and degradation fields. It is found that the method accurately recovers both the traction and degradation fields. Moreover, it is found that neglecting ECM degradation can yield significant errors in traction measurements. Finally, our method can extend the range of context where tractions can be appropriately measured.

42 ENGINEERING↗

POLARIZED SYNCHROTRON EMISSIVITIES AND ABSORPTIVITIES FOR RELATIVISTIC THERMAL, POWER-LAW, AND KAPPA DISTRIBUTION FUNCTIONS

Synchrotron emission and absorption determine the observational appearances of many astronomical systems. In this paper, we describe a numerical scheme for calculating synchrotron emissivities and absorptivities in all four Stokes parameters for arbitrary gyrotropic electron distribution functions, building on earlier work by Leung, Gammie, and Noble. We use this technique to evaluate the emissivities and the absorptivities for a thermal (Maxwell–Jüttner), isotropic power-law, and an isotropic kappa distribution function. The latter contains a power-law tail at high particle energies that smoothly merges with a thermal core at low energies, as is characteristic of observed particle spectra in collisionless plasmas. We provide fitting formulae and error bounds on the fitting formulae for use in codes that solve the radiative transfer equation. The numerical method and the fitting formulae are implemented in a compact C library called symphony. We find that the kappa distribution has a source function that is indistinguishable from a thermal spectrum at low frequency and transitions to the characteristic self-absorbed synchrotron spectrum, ∝ν{sup 5/2}, at high frequency; the linear polarization fraction for a thermal spectrum is near unity at high frequency; and all distributions produce O(10%) circular polarization at low frequency for lines of sight sufficiently close to the magnetic field vector.

79 ASTRONOMY AND ASTROPHYSICS↗

ON ESTIMATING FORCE-FREENESS BASED ON OBSERVED MAGNETOGRAMS

It is a common practice in the solar physics community to test whether or not measured photospheric or chromospheric vector magnetograms are force-free, using the Maxwell stress as a measure. Some previous studies have suggested that magnetic fields of active regions in the solar chromosphere are close to being force-free whereas there is no consistency among previous studies on whether magnetic fields of active regions in the solar photosphere are force-free or not. Here we use three kinds of representative magnetic fields (analytical force-free solutions, modeled solar-like force-free fields, and observed non-force-free fields) to discuss how measurement issues such as limited field of view (FOV), instrument sensitivity, and measurement error could affect the estimation of force-freeness based on observed magnetograms. Unlike previous studies that focus on discussing the effect of limited FOV or instrument sensitivity, our calculation shows that just measurement error alone can significantly influence the results of estimates of force-freeness, due to the fact that measurement errors in horizontal magnetic fields are usually ten times larger than those in vertical fields. This property of measurement errors, interacting with the particular form of a formula for estimating force-freeness, would result in wrong judgments of the force-freeness: a truly force-free field may be mistakenly estimated as being non-force-free and a truly non-force-free field may be estimated as being force-free. Our analysis calls for caution when interpreting estimates of force-freeness based on measured magnetograms, and also suggests that the true photospheric magnetic field may be further away from being force-free than it currently appears to be.

79 ASTRONOMY AND ASTROPHYSICS↗

Accurate solutions, parameter studies and comparisons for the Euler and potential flow equations

Parameter studies are conducted using the Euler and potential flow equation models for unsteady and steady flows in both two and three dimensions. The Euler code is an implicit, upwind, finite volume code which uses the Van Leer method of flux-vector-splitting which has been recently extended for use on dynamic meshes and maintain all the properties of the original splitting. The potential flow code is an implicit, finite difference method for solving the transonic small disturbance equations and incorporates both entropy and vorticity corrections into the solution procedures thereby extending its applicability into regimes where shock strength normally precludes its use. Parameter studies resulting in benchmark type calculations include the effects of spatial and temporal refinement, spatial order of accuracy, far field boundary conditions for steady flow, frequency of oscillation, and the use of subiterations at each time step to reduce linearization and factorization errors. Comparisons between Euler and potential flows results are made as well as with experimental data where available.

Anderson, W. Kyle↗

Accurate solutions, parameter studies and comparisons for the Euler and potential flow equations

Parameter studies are conducted using the Euler and potential flow equation models for steady and unsteady flows in both two and three dimensions. The Euler code is an implicit, upwind, finite volume code which uses the Van Leer method of flux vector splitting which has been recently extended for use on dynamic meshes and maintain all the properties of the original splitting. The potential flow code is an implicit, finite difference method for solving the transonic small disturbance equations and incorporates both entropy and vorticity corrections into the solution procedures thereby extending its applicability into regimes where shock strength normally precludes its use. Parameter studies resulting in benchmark type calculations include the effects of spatial and temporal refinement, spatial order of accuracy, far field boundary conditions for steady flow, frequency of oscillation, and the use of subiterations at each time step to reduce linearization and factorization errors. Comparisons between Euler and potential flow results are made, as well as with experimental data where available.

Anderson, W. Kyle↗