Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “vector fitting”

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 73 records · Page 4

Observations of interstellar helium with a gas absorption cell - Implications for the structure of the local interstellar medium

A photometer sensitive at the 584 A line of He 1, incorporating a helium gas resonance absorption cell, was flown on the Apollo-Soyuz Test Project in July 1975. The instrument observed much of the night-time sky, and returned 42 min of usable data. The data were analyzed by fitting to a model of resonant scattering of solar 584 A flux from nearby interstellar helium. Good model fits were obtained for an interstellar gas bulk velocity vector pointing toward alpha = 72 deg, delta = +15 deg, with speed 20 km/s, with interstellar medium temperatures from 5000 to 20,000 K and with neutral interstellar helium density (8.9 plus or minus 10 to the -3rd/cu cm). In the context of theoretical studies of the interstellar medium by McKee and Ostriker (1977), the results may indicate that the sun lies in the warm, partially ionized periphery of a cold interstellar cloud, surrounded by a high-temperature gas heated by old supernova remnants.

Freeman, J.↗

A New Vehicle-to-Vehicle Communication System: Visual-Enhanced Cooperative Traffic Operations

The advent of Connected and Autonomous Vehicles (CAVs) has highlighted the necessity for robust communication systems between vehicles and their environment. This study introduces a novel vehicle-to-vehicle (V2V) communication system, termed the Visual-Enhanced Cooperative Traffic Operations (VECTOR) system. The VECTOR system addresses the need for robust communication by converting dynamic data (including velocity and yaw angle data) into binary code, which is displayed on an LED panel mounted on the top of the vehicle. Following vehicles detect this panel and decode the information using a camera, implementing a visual-based communication method. VECTOR system employs a comprehensive five-module process. Initially, polynomial fitting techniques are applied to velocity data over fixed time intervals using third-degree polynomials, with validation via R² and MSE metrics. The second module converts velocity and yaw angle data into binary form, thereby enhancing detection and processing efficiency. The third module focuses on improving detection stability across various environmental conditions to enhance traffic safety. The fourth module decodes the binary data back into trajectory information, ensuring the fidelity of velocity and yaw angles. The final module integrates eco-control through the VECTOR system, employing advanced control algorithms to minimize energy consumption in CAVs. Experimental evaluations conducted using a modified CAV test platform based on the Lincoln MKZ demonstrate the feasibility and efficiency of the VECTOR system, achieving a 75% R-squared accuracy rate in replicating original velocity data. This methodology not only highlights potential applications but also underscores significant implications for advancing CAV technology.

Ma, Ke↗

Hybrid NN/SVM Computational System for Optimizing Designs

A computational method and system based on a hybrid of an artificial neural network (NN) and a support vector machine (SVM) (see figure) has been conceived as a means of maximizing or minimizing an objective function, optionally subject to one or more constraints. Such maximization or minimization could be performed, for example, to optimize solve a data-regression or data-classification problem or to optimize a design associated with a response function. A response function can be considered as a subset of a response surface, which is a surface in a vector space of design and performance parameters. A typical example of a design problem that the method and system can be used to solve is that of an airfoil, for which a response function could be the spatial distribution of pressure over the airfoil. In this example, the response surface would describe the pressure distribution as a function of the operating conditions and the geometric parameters of the airfoil. The use of NNs to analyze physical objects in order to optimize their responses under specified physical conditions is well known. NN analysis is suitable for multidimensional interpolation of data that lack structure and enables the representation and optimization of a succession of numerical solutions of increasing complexity or increasing fidelity to the real world. NN analysis is especially useful in helping to satisfy multiple design objectives. Feedforward NNs can be used to make estimates based on nonlinear mathematical models. One difficulty associated with use of a feedforward NN arises from the need for nonlinear optimization to determine connection weights among input, intermediate, and output variables. It can be very expensive to train an NN in cases in which it is necessary to model large amounts of information. Less widely known (in comparison with NNs) are support vector machines (SVMs), which were originally applied in statistical learning theory. In terms that are necessarily oversimplified to fit the scope of this article, an SVM can be characterized as an algorithm that (1) effects a nonlinear mapping of input vectors into a higher-dimensional feature space and (2) involves a dual formulation of governing equations and constraints. One advantageous feature of the SVM approach is that an objective function (which one seeks to minimize to obtain coefficients that define an SVM mathematical model) is convex, so that unlike in the cases of many NN models, any local minimum of an SVM model is also a global minimum.

Rai, Man Mohan↗

Development of iterative techniques for the solution of unsteady compressible viscous flows

During the past two decades, there has been significant progress in the field of numerical simulation of unsteady compressible viscous flows. At present, a variety of solution techniques exist such as the transonic small disturbance analyses (TSD), transonic full potential equation-based methods, unsteady Euler solvers, and unsteady Navier-Stokes solvers. These advances have been made possible by developments in three areas: (1) improved numerical algorithms; (2) automation of body-fitted grid generation schemes; and (3) advanced computer architectures with vector processing and massively parallel processing features. In this work, the GMRES scheme has been considered as a candidate for acceleration of a Newton iteration time marching scheme for unsteady 2-D and 3-D compressible viscous flow calculation; from preliminary calculations, this will provide up to a 65 percent reduction in the computer time requirements over the existing class of explicit and implicit time marching schemes. The proposed method has ben tested on structured grids, but is flexible enough for extension to unstructured grids. The described scheme has been tested only on the current generation of vector processor architecture of the Cray Y/MP class, but should be suitable for adaptation to massively parallel machines.

Hixon, Duane↗

Decoding the B → K ν ν excess at Belle II: Kinematics, operators, and masses

An excess in the branching fraction for B + → K + ν ν recently measured at Belle II may be a hint of new physics. We perform thorough likelihood analyses for different new physics scenarios such as B → K X with a new invisible particle X , or B → K χ χ through a scalar, vector, or tensor current with χ being a new invisible particle or a neutrino. We find that vector-current three-body decay with m X ≃ 0.6 GeV—which may be dark matter—is most favored, while two-body decay with m X ≃ 2 GeV is also competitive. The best-fit branching fractions for the scalar and tensor cases are a few times larger than for the two-body and vector cases. Past measurements provide further discrimination, although the best-fit parameters stay similar. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Mars physical parameters as determined from Mariner 9 observations of the natural satellites and Doppler tracking

Mariner 9 Doppler tracking data and television photographs of Deimos and Phobos were analyzed to determine the gravity field, mass, and spin-axis direction of Mars and the natural satellite orbits. The solutions agree with previously published results. Radio data consisted of an apoapsis state vector for each of revolutions 5-195 obtained from one-revolution fits of Doppler data. Optical data consisted of TV photographs of Phobos and Deimos taken between revolutions 25 and 221. A first-order analytical theory, extended to include dominant second-order resonance effects on the Mariner 9 orbit, was used to calculate the motion of the spacecraft, Deimos, and Phobos. The feasibility of combining radio and optical data in long-arc solutions for accurate determination of orbits and physical parameters is demonstrated. The analytical theory developed for the evolution of a highly eccentric orbit in shallow resonance is accurate to plus or minus 1 km in the apoapsis state vector of Mariner 9 over a period of 200 revolutions.

Born, G. H.↗

Computational methods for the identification of spatially varying stiffness and damping in beams

A numerical approximation scheme for the estimation of functional parameters in Euler-Bernoulli models for the transverse vibration of flexible beams with tip bodies is developed. The method permits the identification of spatially varying flexural stiffness and Voigt-Kelvin viscoelastic damping coefficients which appear in the hybrid system of ordinary and partial differential equations and boundary conditions describing the dynamics of such structures. An inverse problem is formulated as a least squares fit to data subject to constraints in the form of a vector system of abstract first order evolution equations. Spline-based finite element approximations are used to finite dimensionalize the problem. Theoretical convergence results are given and numerical studies carried out on both conventional (serial) and vector computers are discussed.

Banks, H. T.↗

Ginzburg-Landau theory for the solid-liquid interface of bcc elements

Consideration is given to a simple order-parameter theory for the interfacial tension of body-centered-cubic solids in which the principal order parameter is the amplitude of the density wave at the smallest nonzero reciprocal-lattice vector of the solid. The parameters included in the theory are fitted to the measured heat of fusion, melting temperature, and solid-liquid density difference, and to the liquid structure factor and its temperature derivative at freezing. Good agreement is found with experiment for Na and Fe and the calculated anisotropy of the surface tension among different crystal faces is of the order of 2 percent. On the basis of various assumptions about the universal behavior of bcc crystals at melting, the formalism predicts that the surface tension is proportional to the heat of fusion per surface atom.

Shih, W. H.↗

A model for the characterization of the spatial properties in vestibular neurons

Quantitative study of the static and dynamic response properties of some otolith-sensitive neurons has been difficult in the past partly because their responses to different linear acceleration vectors exhibited no "null" plane and a dependence of phase on stimulus orientation. The theoretical formulation of the response ellipse provides a quantitative way to estimate the spatio-temporal properties of such neurons. Its semi-major axis gives the direction of the polarization vector (i.e., direction of maximal sensitivity) and it estimates the neuronal response for stimulation along that direction. In addition, the semi-minor axis of the ellipse provides an estimate of the neuron's maximal sensitivity in the "null" plane. In this paper, extracellular recordings from otolith-sensitive vestibular nuclei neurons in decerebrate rats were used to demonstrate the practical application of the method. The experimentally observed gain and phase dependence on the orientation angle of the acceleration vector in a head-horizontal plane was described and satisfactorily fit by the response ellipse model. In addition, the model satisfactorily fits neuronal responses in three-dimensions and unequivocally demonstrates that the response ellipse formulation is the general approach to describe quantitatively the spatial properties of vestibular neurons.

Non-NASA Center↗

Magnetic Field Vector Measurements Using Doppler-Free Saturation Spectroscopy

Experimentally measuring the magnetic-field vector topology in low-cost fusion-concept devices will provide critical information required to optimize and accelerate the development of those concepts for economical fusion energy. The magnetic-field measurement techniques most commonly used in modern tokamaks, i.e., beam-emission spectroscopy and optical emission spectroscopy, are not suitable below 5000 Gauss due to a lack of sensitivity. To diagnose magnetic field relevant to some of the ARPA-E-supported fusion concepts (50 to 1000 Gauss), a novel approach is required. Doppler-free saturation spectroscopy (DFSS) is a pump/probe laser-based technique and has been successfully demonstrated in the laboratory to provide a magnetic-field measurement accuracy of <10 Gauss. DFSS is implemented to obtain ultra-high-resolution spectra at a localized point in space. The spectral data is then fit to the Schrödinger equation for determination of the magnetic-field vector. To obtain the magnetic-field topology, the measurement spatial location is scanned over the desired region of space by steering the DFSS pump laser beam. The goal of this project is to assemble and commission a portable DFSS diagnostic capable of 2-D magnetic-field vector measurements utilizing the Hα spectral line profile.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Newton algorithm for complex curve fitting

The problem of synthesizing transfer functions from frequency response measurements is considered. Given a complex vector representing the measured frequency response of a physical system, a transfer function of specified order is determined that minimizes the sum of the magnitude-squared of the frequency response errors. This nonlinear least squares minimization problem is solved by an iterative global descent algorithm of the Newton type which converges quadratically near the minimum. The unknown transfer function is expressed as a sum of second order rational polynomials, a parameterization that facilitates a numerically robust computer implementation. The algorithm is developed for single-input, single-output, causal, stable transfer functions.

Spanos, J. T.↗

A model for the motion of the Philippine Sea plate consistent with NUVEL-1 and geological data

We investigate angular velocity vectors of the Philippine Sea (PH) plate relative to the adjacent major plates, Eurasia (EU) and Pacific (PA), and the smaller Caroline (CR) plate. Earthquake slip vector data along the Philippine Sea plate are inverted, subject to the constraint that EU-PA motion equals that predicted by the global relative plate model NUVEL-1. The resulting solution fails to satisfy geological constraints along the Caroline-Pacific boundary: convergence along the Mussau Trench and divergence along the Sorol Trough. We then seek solutions satisfying both the CR-PA boundary conditions and the Philippine Sea slip vector data, by adjusting the PA-PH and EU-PH best fitting poles within their error ellipses. We also consider northern Honshu to be part of the North American plate and impose the constraint that the Philippine Sea plate subducts beneath northern Honshu along the Sagmi Trough in a NNW-NW direction. Of the solutions satisfying these conditions, we select the best EU-PH as 48.2 deg N, 157.0 deg E, 1.09 deg/my, corresponding to a pole far from Japan and south of Kamchatka, and PA-PH, 1.2 deg N, 134.2 deg E, 1.00 deg/my. Predicted NA-PH and EU-PH convergence rates in central Honshu are consistent with estimated seismic slip rates. Previous estimates of the EU-PH pole close to central Honshu are inconsistent with extension within the Bonin backarc implied by earthquake slip vectors and NNW-NW convergence of the Bonin forearc at the Sagami Trough.

Seno, Tetsuzo↗

Algorithms for l2 and l-infinity transfer function curve fitting

In this paper algorithms for fitting transfer functions to frequency response data are developed. Given a complex vector representing the measured frequency response of a physical system, a transfer function of specified order is determined that minimizes either of the following criteria: (1) the sum of the magnitude-squared of the frequency response errors, and (2) the magnitude of the maximum error. Both of these criteria are nonlinear in the coefficients of the unknown transfer function, and iterative minimization algorithms are proposed. A numerical example demonstrates the effectiveness of the proposed algorithms.

Spanos, John T.↗

Light Detecting and Ranging (LIDAR) for in-situ heliostat optical error assessment

A system and method for optical assessment of a heliostat includes obtaining a point cloud data representing an image of the heliostat; isolating the data; filtering and fitting the filtered heliostat data to a bounding box; translating the heliostat data to a plane to aid in segmentation; segmenting a plurality of facets of the heliostat fitting each of the segmented facets to a respective plane; generating normal vectors characterizing each of the plurality of facets; and calculating a canting angle associated with each respective facet of the plurality of facets. A heliostat with mirrored facets and a scanner are provided. The scanner captures point cloud data representing the heliostat, which is segmented for each facet. Normal vectors characterize the facets and a canting angle is calculated for the respective facet.

Small, Daniel E.↗

Backus Effect on a Perpendicular Errors in Harmonic Models of Real vs. Synthetic Data

Measurements of geomagnetic scalar intensity on a thin spherical shell alone are not enough to separate internal from external source fields; moreover, such scalar data are not enough for accurate modeling of the vector field from internal sources because of unmodeled fields and small data errors. Spherical harmonic models of the geomagnetic potential fitted to scalar data alone therefore suffer from well-understood Backus effect and perpendicular errors. Curiously, errors in some models of simulated 'data' are very much less than those in models of real data. We analyze select Magsat vector and scalar measurements separately to illustrate Backus effect and perpendicular errors in models of real scalar data. By using a model to synthesize 'data' at the observation points, and by adding various types of 'noise', we illustrate such errors in models of synthetic 'data'. Perpendicular errors prove quite sensitive to the maximum degree in the spherical harmonic expansion of the potential field model fitted to the scalar data. Small errors in models of synthetic 'data' are found to be an artifact of matched truncation levels. For example, consider scalar synthetic 'data' computed from a degree 14 model. A degree 14 model fitted to such synthetic 'data' yields negligible error, but amplifies 4 nT (rmss) added noise into a 60 nT error (rmss); however, a degree 12 model fitted to the noisy 'data' suffers a 492 nT error (rmms through degree 12). Geomagnetic measurements remain unaware of model truncation, so the small errors indicated by some simulations cannot be realized in practice. Errors in models fitted to scalar data alone approach 1000 nT (rmss) and several thousand nT (maximum).

Voorhies, C. V.↗

Precision measurements of Higgs hadronic decay modes at the FCC-ee

The expected precision at the FCC-ee on the product $σ$ x $\mathcal{B} \large{(}H → b\bar{b}, c\bar{c}, s\bar{s}, gg\large{)}$ of Higgs boson production cross sections times branching ratios of hadronic decays is presented. This study provides the first comprehensive determination of all major hadronic Higgs decay modes in a combined fit at future $e^+e^−$ colliders, using both Higgs-strahlung $(ZH)$ and Vector boson fusion $(v\bar{v}H)$ production processes, with a full treatment of interference effects in the $v\bar{v}jj$ final state. It assumes four identical IDEA detectors collecting $e^+e^−$ collisions at $\sqrt{s}$ = 240 and 365 GeV. The combination of all channels across both energies, with full covariance between production and decay modes, yields a production cross-section times branching-ratio precision at the percent to per-mil level for the dominant hadronic final states $(b\bar{b}, c\bar{c}, gg)$. These results provide a comprehensive input to the determination of Higgs coupling projections at the FCC-ee, and they establish for the first time sensitivity to the rare decay $H → s\bar{s}$, demonstrating that FCC-ee has the potential to provide evidence of the strange-quark Yukawa coupling.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Computer solutions for studying correlations between solar magnetic fields and Skylab X-ray observations

A method is described which correlates the NASA-Marshall Space Flight Center (MSFC) Image Data Processing System (IDAPS) and MSFC magnetograph data to X-ray and H-alpha observations from the Skylab mission. Solutions of Laplace's equation in three dimensions, based on the magnetograph data, are convolved with observed X-ray and H-alpha regions. Matched filtering (template matching) provides a best fit of the observed X-ray regions to the computed total magnetic vector magnitude between 10,000 and 15,000 km above the photosphere.

Teuber, D.↗

Extension of state-vector splitting to the Navier-Stokes equations

A multi-dimensional upwind scheme for solving the Euler equations has been extended to solve the Navier-Stokes equations. These methods share the property of not requiring body-fitted computational grids. A time advance scheme which is second-order accurate in time and space has been developed for use in conjunction with state-vector splitting to accurately resolve regions with significant gradients in flow quantities. The scheme has been tested on several simple viscous flow problems for which analytic solutions are known; the agreement of computation with theory is excellent.

Gooch, Carl F.↗