Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “orthogonal weights”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

Creating a Test-Validated Finite-Element Model of the X-56A Aircraft Structure

Small modeling errors in a finite-element model will eventually induce errors in the structural flexibility and mass, thus propagating into unpredictable errors in the unsteady aerodynamics and the control law design. One of the primary objectives of the X-56A Multi-Utility Technology Testbed aircraft is the flight demonstration of active flutter suppression and, therefore, in this study, the identification of the primary and secondary modes for the structural model tuning based on the flutter analysis of the X-56A aircraft. The ground-vibration test-validated structural dynamic finite-element model of the X-56A aircraft is created in this study. The structural dynamic finite-element model of the X-56A aircraft is improved using a model-tuning tool. In this study, two different weight configurations of the X-56A aircraft have been improved in a single optimization run. Frequency and the cross-orthogonality (mode shape) matrix were the primary focus for improvement, whereas other properties such as c.g. location, total weight, and off-diagonal terms of the mass orthogonality matrix were used as constraints. The end result was an improved structural dynamic finite-element model configuration for the X-56A aircraft. Improved frequencies and mode shapes in this study increased average flutter speeds of the X-56A aircraft by 7.6% compared to the baseline model.

ground vibration test↗

Ultra-Stable Beacon Source for Laboratory Testing of Optical Tracking

The ultra-stable beacon source (USBS) provides a laser-beam output with a very low angular jitter and can be used as an absolute angular reference to simulate a beacon in the laboratory. The laser is mounted on the top of a very short (approximately equal to 1 m) inverted pendulum (IP) with its optical axis parallel to the carbon fiber pendulum leg. The 85-cm, carbon fiber rods making up the leg are very lightweight and rigid, and are supported by a flex-joint at the bottom (see figure). The gimbal-mounted laser is a weight-adjustable load of about 1.5 kg with its center of rotation co-located with the center of percussion of the inverted pendulum. This reduces the coupling of transverse motion at the base of the pendulum to angular motion of the laser at the top. The inverted pendulum is mounted on a gimbal with its center of rotation coinciding with the pivot position of the inverted pendulum flexure joint. This reduces coupling of ground tilt at the inverted pendulum base to motion of the laser mounted at the top. The mass of the top gimbal is adjusted to give the pendulum a very low resonant frequency (approximately equal to 10 mHz) that filters transverse seismic disturbances from the ground where the base is attached. The motion of the IP is monitored by an optical-lever sensor. The laser light is reflected by the mirror on the IP, and then is detected by a quadrant photo-detector (QPD). The position of the beam spot on the QPD corresponds to the tilt of the IP. Damping of this motion is provided by two coil and magnet pairs. The bottom gimbal mount consists of two plates. The IP is mounted on the second plate. The first plate is supported by two posts through needles and can be rotated about the axis connecting the tips of the needles. The second plate hangs from the first plate and can be rotated about the axis perpendicular to the first plate. As a result, the second plate acts as a two-axis rotation stage. Its center of rotation is located at the effective bending point of the flex-joint. The second plate is pressed against two screw actuators by the weight of the IP. The screw actuators are orthogonal to each other and are used to adjust the inclination of the second plate. The actuators are driven by stepper motors. The whole IP system is housed in a box made of Lexan plastic plates to provide isolation from air currents and temperature variations. The signals from the sensors are processed and recorded with a PC using the xPC Target realtime environment of Math- Works. The control algorithms are written using the Simulink package from The MathWorks.

Aso, Yoichi↗

GENIE - Generation of computational geometry-grids for internal-external flow configurations

Progress realized in the development of a master geometry-grid generation code GENIE is presented. The grid refinement process is enhanced by developing strategies to utilize bezier curves/surfaces and splines along with weighted transfinite interpolation technique and by formulating new forcing function for the elliptic solver based on the minimization of a non-orthogonality functional. A two step grid adaptation procedure is developed by optimally blending adaptive weightings with weighted transfinite interpolation technique. Examples of 2D-3D grids are provided to illustrate the success of these methods.

Soni, B. K.↗

Polar decomposition for attitude determination from vector observations

This work treats the problem of weighted least squares fitting of a 3D Euclidean-coordinate transformation matrix to a set of unit vectors measured in the reference and transformed coordinates. A closed-form analytic solution to the problem is re-derived. The fact that the solution is the closest orthogonal matrix to some matrix defined on the measured vectors and their weights is clearly demonstrated. Several known algorithms for computing the analytic closed form solution are considered. An algorithm is discussed which is based on the polar decomposition of matrices into the closest unitary matrix to the decomposed matrix and a Hermitian matrix. A somewhat longer improved algorithm is suggested too. A comparison of several algorithms is carried out using simulated data as well as real data from the Upper Atmosphere Research Satellite. The comparison is based on accuracy and time consumption. It is concluded that the algorithms based on polar decomposition yield a simple although somewhat less accurate solution. The precision of the latter algorithms increase with the number of the measured vectors and with the accuracy of their measurement.

Bar-Itzhack, Itzhack Y.↗

Polar Plate Theory for Orthogonal Anisotropy

Laminated fiber-reinforced (or filamentary) composites are used today for their high strength-to-weight and stiffness-to-weight ratios. However, because of the anisotropic behavior of composites, determining the response on a macroscopic scale is challenging. This is particularly evident in the evaluation of the governing differential equations of a circular disk with the fibers of the lamina oriented with rectilinear orthogonality. This includes any situation involving a composite plate of circular geometry in which out-of-plane displacements due to load are desired, such as fastener pull through loading of a composite plate. Current analysis techniques use numerical methods with rectilinear coordinate systems to solve problems with circular geometry. These analyses over predict plate stiffness by 20% and underpredict failure by 70%. Consequently, there is a need to transform classical composite plate theory to a polar coordinate system. In order to better analyze structures with circular geometries the classical composite plate equations are transformed into the plate equations for a rectilinearly anisotropic composite in polar coordinates. A composite plate is typically a laminate of fibers in rectilinear directions. Subsequent to the lay-tip the necessary geometry is cut out of a rectangular plate. In a similar manner, the derivation of the plate equation starts with the fundamental definitions of strain, displacement and curvature and incorporates the material property angular dependence into the equilibrium equations for a differential polar element. In the transformed state, the stiffness coefficients are no longer constant, adding to the complexity of the governing differential equations. This paper discusses the new derivation and evaluation of the plate equations for a circular composite disk with orthogonal rectilinear anisotropy. The resultant new three partial differential equations, which describe the circular anisotropic plate, can be used to evaluate out-of-plane displacements for given load conditions in design of composite Structures. Without the formation of these unique plate equations, design of such structures is not a precise engineering accomplishment due to the lack of precise design tools. Such structures are over designed to compensate and costly tests need to be performed after production to validate safety. Though impossible to present the mathematical solution within the limits of this paper, the application of the newly derived plate equation in its orthotropic form for a balanced symmetric laminate compared favorably with test results and surpassed the results of standard numerical methods. Within the limitations of the problem foundation the predictive model provides a mechanism for establishing out-of-plane deflection levels for a circular composite plate. Future work in this area should include incorporating the circular composite plate equations in Finite Element Models. The new anisotropic polar plate equations call be utilized as a design tool so that the design more accurately meets the requirements, thus taking full advantage of the weight savings of composites. These accomplishments have not previously been presented to the engineering community nor ever proposed to the IAF. Should the proposed paper be accepted for presentation, attendance is assured.

Bailey, Michelle D.↗

Model-free distributed learning

Model-free learning for synchronous and asynchronous quasi-static networks is presented. The network weights are continuously perturbed, while the time-varying performance index is measured and correlated with the perturbation signals; the correlation output determines the changes in the weights. The perturbation may be either via noise sources or orthogonal signals. The invariance to detailed network structure mitigates large variability between supposedly identical networks as well as implementation defects. This local, regular, and completely distributed mechanism requires no central control and involves only a few global signals. Thus it allows for integrated on-chip learning in large analog and optical networks.

Dembo, Amir↗

Global coordinates with centimeter accuracy in the International Terrestrial Reference Frame using GPS

Using 21 days of GPS data from 21 globally distributed receivers operating during early 1991, a 7-parameter transformation between a GPS free-network solution and coordinates of 12 stations listed in the International Terrestrial Reference Frame (ITRF) is solved. Standard errors of GPS coordinates are derived by applying an orthogonal projection operator to the free-network covariance. The weighted rms difference between 33 transformed GPS and ITRF coordinates is 12 mm in the Northern Hemisphere. Best results are obtained by mapping ITRF coordinates to the epoch of this experiment assuming no vertical site motions. Fixing selected sites in the GPS solution to ITRF '90 does not improve the agreement. It is concluded that the use of fiducial constraints is unnecessary for global networks.

Blewitt, Geoffrey↗

Innovative mechanism for measuring the mass properties of an object

The Kennedy Space Center Robotics Group recently completed development and testing on a novel approach to measure the mass properties of a rigid body. This unique design can measure the payload's weight, mass center location, and moments of inertia about three orthogonal axes. Furthermore, these measurements only require a single torque sensor and a single angular position sensor.

Wolcott, Kedron R.↗

Kravchuk functions for the finite oscillator approximation

Kravchuk orthogonal functions - Kravchuk polynomials multiplied by the square root of the weight function - simplify the inversion algorithm for the analysis of discrete, finite signals in harmonic oscillator components. They can be regarded as the best approximation set. As the number of sampling points increases, the Kravchuk expansion becomes the standard oscillator expansion.

Atakishiyev, Natig M.↗

Estimation of Non-Seasonal Oceanic Mass Redistribution and its Geodynamic Effects from the Topex/Poseidon Observation

Mass redistribution in terms of angular momentum will change the Earth's rotation under the conservation of angular momentum. It also changes the external gravitational field according to Newton's gravitational law. These small geodynamic variations can be measured very accurately by modem space geodetic techniques. We report on new estimates of the oceanic tidal angular momentum and hence the tidal variations in the Earth's rotation rate and polar motion. These are deduced from the new tidal solutions GOT99 and TPXO.3, plus M2 and Mf estimates (with error bars) from TPXO.4a. Comparisons with observed diurnal/semidiurnal tidal variations are generally similar to previous model estimates. Comparisons of models and observations of Mf polar motion are confused, with neither models nor observations agreeing among themselves or with each other. Part of this is caused by differences in Mf currents; model TPXO.4 shows current-intensive Rossby modes along western boundaries of major basins, which are indeed expected on theoretical grounds. Unfortunately, there are no reliable in situ observations of Mf currents useful for comparing models. On non-tidal oceanic mass redistribution we compute the principal components of the multi-year T/P sea surface height (SSH) field using EOF (empirical orthogonal functions) taking into consideration of the latitude-dependent area weighting. Emphasis is placed upon non-seasonal and interannual oscillations on regional scales, including ENSO and the North Atlantic Oscillation. Steric corrections to the observed SSH due to sea surface temperature changes (which has no geodynamic effects) is now under study using POCM ocean model output.

Chao, Benjamin F.↗

In search of the elusive flow field

A specific set of problems must be overcome before an optical flow field can be derived from an image sequence. These problems are discussed and an outline of a scheme for overcoming these obstacles is presented. A two-dimensional motion sensor is presented which is capable of measuring edge velocity in twelve separate directions at any point in the image. However, because of the aperture problem, the sensors can only determine the velocity component in the direction orthogonal to the edge. A second stage which uses a cosine-weighted voting scheme is used to derive the true flow field from the edge velocities. The whole system is based on a parallel network of connections and weights and is potentially a very fast method of obtaining the two-dimensional flow field.

Perrone, John A.↗

A Posteriori Error Estimation for Finite Volume and Finite Element Approximations Using Broken Space Approximation

We consider a posteriori error estimates for finite volume and finite element methods on arbitrary meshes subject to prescribed error functionals. Error estimates of this type are useful in a number of computational settings: (1) quantitative prediction of the numerical solution error, (2) adaptive meshing, and (3) load balancing of work on parallel computing architectures. Our analysis recasts the class of Godunov finite volumes schemes as a particular form of discontinuous Galerkin method utilizing broken space approximation obtained via reconstruction of cell-averaged data. In this general framework, weighted residual error bounds are readily obtained using duality arguments and Galerkin orthogonality. Additional consideration is given to issues such as nonlinearity, efficiency, and the relationship to other existing methods. Numerical examples are given throughout the talk to demonstrate the sharpness of the estimates and efficiency of the techniques. Additional information is contained in the original.

Barth, Timothy J.↗

Programs for computing abscissas and weights for classical and nonclassical Gaussian quadrature formulas

Computer programs for computing Gaussian quadrature abscissas and weights are described. For the classical case the programs use Laguerre iteration to compute abscissas as zeros of orthogonal polynomials. The polynomials are evaluated from known recursion coefficients. The nonclassical case is handled similarly except that the recursion coefficients are computed by numerical integration. A sample problem, with input and output, is presented to illustrate the use of the programs. It computes the quadrature abscissas and weights associated with the weight function over the interval (0,1) for quadrature orders from 16 to 96 in increments of 8.

Desmarais, R. N.↗

Testing the Cosmological Principle: Astrometric Limits on Systemic Motion of Quasars at Different Cosmological Epochs

A sample of 60,410 bona fide optical quasars with astrometric proper motions in Gaia Early Data Release 3 and spectroscopic redshifts above 0.5 in an oval 8400 square degree area of the sky is constructed. Using orthogonal Zernike functions of polar coordinates, the proper motion fields are fitted in a weighted least-squares adjustment of the entire sample and of six equal bins of sorted redshifts. The overall fit with 37 Zernike functions reveals a statistically significant pattern, which is likely to be of instrumental origin. The main feature of this pattern is a chain of peaks and dips mostly in the R.A. component with an amplitude of 25 μ as yr -1 . This field is subtracted from each of the six analogous fits for quasars grouped by redshifts covering the range 0.5 through 7.03, with median values of 0.72, 1.00, 1.25, 1.52, 1.83, 2.34. The resulting residual patterns are noisier, with formal uncertainties up to 8 μ as yr -1 in the central part of the area. We detect a single high-confidence Zernike term for the R.A. proper motion components of quasars with redshifts around 1.52 representing a general gradient of 30 μ as yr -1 over 150° on the sky. We do not find any small- or medium-scale systemic variations of the residual proper motion field as functions of redshift above the 2.5 σ significance level.

79 ASTRONOMY AND ASTROPHYSICS↗

Oblique convergence and deformation along the Kuril and Japan trenches

The hypothesis that present-day deformation within the southern Kuril forearc is driven by oblique subduction of the Pacific plate is tested using 397 horizontal slip directions derived from shallow-thrust earthquakes from the Kuril and Japan trenches for the period 1963-1991. A simple two-plate model fits the 397 slip vectors significantly worse than a model that permits strike-slip motion of the southern Kuril forearc relative to the overlying plate. Weighted, mean slip directions along the southern Kuril trench are systematically rotated toward the direction orthogonal to the trench, which implies that the net convergence is partitioned into less oblique subduction and trench-parallel displacement of the southern Kuril forearc. The angular discrepancy between the observed slip direction and the direction predicted by the NUVEL-1 Pacific-North America Euler vector implies that the southern Kuril forearc translates 6-11 mm/yr to the southwest relative to the overlying North American plate. These results are consistent with geologically, geodetically, and seismologically observed convergence at the leading edge of the forearc sliver in southern Kokkaido and with previously inferred extension at the trailing edge of the sliver, which is located at the Bussol Strait at 46 deg N.

Demets, Charles↗

Non-Tidal Non-Seasonal Oceanic Mass Redistribution Estimated from the TOPEX/Poseidon Observation

Topex/Poseidon altimetry data are used to estimate the non-tidal mass redistribution as a function of space-time. The goal is to study the contribution of ocean circulations in the geodynamic effects including Earth's rotational and gravitational variations. We examine the non-seasonal anomalies at monthly sampling rate over the T/P span of eight years, concentrating especially on interannual variabilities. Since the sea-surface height data obtained from altimetry is the combined effect of steric change (primarily thermal effect) and the mass flux, and because the former has zero contribution to the geodynamic effects, one needs to do a so-called steric correction by removing from the altimetry data the steric contributions. We achieve it using multiyear monthly sea-surface temperature maps, together with monthly "climatology" for mixed-layer depth maps for mean-months of the year. We analyze both sets of the altimetry data and the steric correction data using the empirical orthogonal function/principal component analysis (in which we take care of issues associated with the area-weighting and non-zero mean), and examine the most important modes, either globally or regionally. In particular, the ENSO in the tropical Pacific/Indian Oceans exhibits the most prominent pattern. The net mass transport after the steric correction can then be compared with: (1) ocean general circulation model outputs for the same period of time (such as from POCM-4B); (2) non-atmospheric Earth rotation variations obtained from space geodesy data and atmospheric angular momentum data; (3) non-atmospheric low-degree gravitational variations from satellite-laser-ranging observations and global atmospheric data.

Chao, Benjamin F.↗

Proper orthogonal descriptors for efficient and accurate interatomic potentials

Here, we present the proper orthogonal descriptors for efficient and accuracy representation of the potential energy surface. The potential energy surface is represented as a many-body expansion of parametrized potentials in which the potentials are functions of atom positions and parameters. The proper orthogonal decomposition is employed to decompose the parametrized potentials into a set of proper orthogonal descriptors (PODs). Because of the rapid convergence of the proper orthogonal decomposition, relevant snapshots can be sampled exhaustively to represent the atomic neighborhood environment accurately with a small number of descriptors. The proper orthogonal descriptors are used to develop interatomic potentials by using a linear expansion of the descriptors and determining the expansion coefficients from a weighted least-squares regression against a density functional theory (DFT) training set. We present a comprehensive evaluation of the POD potentials on previously published DFT data sets comprising Li, Mo, Cu, Ni, Si, Ge, and Ta elements. The data sets represent a diverse pool of metals, transition metals, and semiconductors. The accuracy of the POD potentials are comparable to that of state-of-the-art machine learning potentials such as the spectral neighbor analysis potential (SNAP) and the atomic cluster expansion (ACE).

97 MATHEMATICS AND COMPUTING↗

Orthogonal Polynomials Defined by Self-Similar Measures with Overlaps

Here, we study orthogonal polynomials with respect to self-similar measures, focusing on the class of infinite Bernoulli convolutions, which are defined by iterated function systems with overlaps, especially those defined by the Pisot, Garsia, and Salem numbers. By using an algorithm of Mantica, we obtain graphs of the coefficients of the 3-term recursion relation defining the orthogonal polynomials. We use these graphs to predict whether the singular infinite Bernoulli convolutions belong to the Nevai class. Based on our numerical results, we conjecture that all infinite Bernoulli convolutions with contraction ratios greater than or equal to 1/2 belong to Nevai’s class, regardless of the probability weights assigned to the self-similar measures.

97 MATHEMATICS AND COMPUTING↗