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 325 records · Page 18

Isospin-breaking effects in the three-pion contribution to hadronic vacuum polarization

Isospin-breaking (IB) effects are required for an evaluation of hadronic vacuum polarization at subpercent precision. While the dominant contributions arise from the e + e - → π + π - channel, also IB in the subleading channels can become relevant for a detailed understanding, e.g., of the comparison to lattice QCD. Here, we provide such an analysis for e + e - → 3π by extending our dispersive description of the process, including estimates of final-state radiation (FSR) and ρ–ω mixing. In particular, we develop a formalism to capture the leading infrared-enhanced effects in terms of a correction factor η3π that generalizes the analog treatment of virtual and final-state photons in the 2π case. The global fit to the e+e- → 3π data base, subject to constraints from analyticity, unitarity, and the chiral anomaly, gives for the total 3π contribution to the anomalous magnetic moment of the muon, of which ${a}^{FSR}_\mu [3\pi]$ = 0.51(1) x 10 -10 and ${a}^{\rho - \omega}_\mu [3\pi]$ = -2.68(70) x 10 -10 can be ascribed to IB. We argue that the resulting cancellation with ρ–ω mixing in e + e - → 2π can be understood from a narrow-resonance picture, and provide updated values for the vacuum-polarization-subtracted vector-meson parameters M ω = 782.70(3) MeV, M Φ = 1019.21(2) MeV, Γ ω = 8.71(3) MeV, and Γ Φ = 4.27(1) MeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Conservation equations and physical models for hypersonic air flows in thermal and chemical nonequilibrium

The conservation equations for simulating hypersonic flows in thermal and chemical nonequilibrium and details of the associated physical models are presented. These details include the curve fits used for defining thermodynamic properties of the 11 species air model, curve fits for collision cross sections, expressions for transport properties, the chemical kinetics models, and the vibrational and electronic energy relaxation models. The expressions are formulated in the context of either a two or three temperature model. Greater emphasis is placed on the two temperature model in which it is assumed that the translational and rotational energy models are in equilibrium at the translational temperature, T, and the vibrational, electronic, and electron translational energy modes are in equilibrium at the vibrational temperature, T sub v. The eigenvalues and eigenvectors associated with the Jacobian of the flux vector are also presented in order to accommodate the upwind based numerical solutions of the complete equation set.

Peter A. Gnoffo↗

Dipolar Spin Ice Regime Proximate to an All-In-All-Out Néel Ground State in the Dipolar-Octupolar Pyrochlore Ce 2 Sn 2 O 7

The dipolar-octupolar (DO) pyrochlores, R 2 M 2 O 7 ( R = Ce , Sm , Nd ), are key players in the search for realizable novel quantum spin liquid (QSL) states as a large parameter space within the DO pyrochlore phase diagram is theorized to host QSL states of both dipolar and octupolar nature. New single crystals and powders of Ce 2 Sn 2 O 7 , synthesized by hydrothermal techniques, present an opportunity for a new characterization of the exchange parameters in Ce 2 Sn 2 O 7 using the near-neighbor X Y Z model Hamiltonian associated with DO pyrochlores. Utilizing quantum numerical linked cluster expansion fits to heat capacity and magnetic susceptibility measurements, and classical Monte Carlo calculations to the diffuse neutron diffraction of the new hydrothermally grown Ce 2 Sn 2 O 7 samples, we place Ce 2 Sn 2 O 7 ’s ground state within the ordered dipolar all-in-all-out (AIAO) Néel phase, with quantum Monte Carlo calculations showing a transition to long-range order at temperatures below those accessed experimentally. Indeed, our new neutron diffraction measurements on the hydrothermally grown Ce 2 Sn 2 O 7 powders show a broad signal at low scattering wave vectors, reminiscent of a spin ice, in striking contrast from previous powder neutron diffraction on samples grown from solid-state synthesis, which found diffuse scattering at high scattering wave vectors associated with magnetic and suggested an octupolar quantum spin ice state. We conclude that new hydrothermally grown Ce 2 Sn 2 O 7 samples host a finite-temperature proximate dipolar spin ice phase, above the expected transition to AIAO Néel order. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Amplitude Analysis of $\omega\pi^0$ Photoproduction at GlueX

Excited meson states can often lie hidden within mass spectra beneath more dominant resonances, making it difficult to extract their physical properties. We can unveil these states through amplitude analysis, which disentangles the overlapping states via fits to their unique production and decay angular distributions. Understanding the light-meson spectrum is essential for confirming Lattice QCD predictions, especially in regards to the search for possible hybrid mesons. The Gluonic Excitation (GlueX) experiment at Jefferson Lab aids this search by studying the production of excited light mesons in $\gamma p$ interactions in the 8.2 - 8.8 GeV photon-beam energy range. We show in these proceedings preliminary results of an amplitude analysis of a large $\omega\pi^0$ dataset, concentrating on measuring the interference between the $b_1(1235)$ and an excited $J^{PC}=1^{--}$ vector resonance.

Scheuer, Kevin↗

GPS-Based Reduced Dynamic Orbit Determination Using Accelerometer Data

Currently two gravity field satellite missions, CHAMP and GRACE, are equipped with high sensitivity electrostatic accelerometers, measuring the non-conservative forces acting on the spacecraft in three orthogonal directions. During the gravity field recovery these measurements help to separate gravitational and non-gravitational contributions in the observed orbit perturbations. For precise orbit determination purposes all these missions have a dual-frequency GPS receiver on board. The reduced dynamic technique combines the dense and accurate GPS observations with physical models of the forces acting on the spacecraft, complemented by empirical accelerations, which are stochastic parameters adjusted in the orbit determination process. When the spacecraft carries an accelerometer, these measured accelerations can be used to replace the models of the non-conservative forces, such as air drag and solar radiation pressure. This approach is implemented in a batch least-squares estimator of the GPS High Precision Orbit Determination Software Tools (GHOST), developed at DLR/GSOC and DEOS. It is extensively tested with data of the CHAMP and GRACE satellites. As accelerometer observations typically can be affected by an unknown scale factor and bias in each measurement direction, they require calibration during processing. Therefore the estimated state vector is augmented with six parameters: a scale and bias factor for the three axes. In order to converge efficiently to a good solution, reasonable a priori values for the bias factor are necessary. These are calculated by combining the mean value of the accelerometer observations with the mean value of the non-conservative force models and empirical accelerations, estimated when using these models. When replacing the non-conservative force models with accelerometer observations and still estimating empirical accelerations, a good orbit precision is achieved. 100 days of GRACE B data processing results in a mean orbit fit of a few centimeters with respect to high-quality JPL reference orbits. This shows a slightly better consistency compared to the case when using force models. A purely dynamic orbit, without estimating empirical accelerations thus only adjusting six state parameters and the bias and scale factors, gives an orbit fit for the GRACE B test case below the decimeter level. The in orbit calibrated accelerometer observations can be used to validate the modelled accelerations and estimated empirical accelerations computed with the GHOST tools. In along track direction they show the best resemblance, with a mean correlation coefficient of 93% for the same period. In radial and normal direction the correlation is smaller. During days of high solar activity the benefit of using accelerometer observations is clearly visible. The observations during these days show fluctuations which the modelled and empirical accelerations can not follow.

VanHelleputte, Tom↗

Nonlinear Matrix Approximation with Radial Basis Function Components

We introduce and investigate matrix approximation by decomposition into a sum of radial basis function (RBF) components. An RBF component is a generalization of the outer product between a pair of vectors, where an RBF function replaces the scalar multiplication between individual vector elements. Even though the RBF functions are positive definite, the summation across components is not restricted to convex combinations and allows us to compute the decomposition for any real matrix that is not necessarily symmetric or positive definite. We formulate the problem of seeking such a decomposition as an optimization problem with a nonlinear and non-convex loss function. Several modern versions of the gradient descent method, including their scalable stochastic counterparts, are used to solve this problem. We provide extensive empirical evidence of the effectiveness of the RBF decomposition and that of the gradient-based fitting algorithm. While being conceptually motivated by singular value decomposition (SVD), our proposed nonlinear counterpart outperforms SVD by drastically reducing the memory required to approximate a data matrix with the same L2 error for a wide range of matrix types. For example, it leads to 2 to 6 times memory save for Gaussian noise, graph adjacency matrices, and kernel matrices. Moreover, this proximity-based decomposition can offer additional interpretability in applications that involve, e.g., capturing the inner low-dimensional structure of the data, retaining graph connectivity structure, and preserving the acutance of images.

Rebrova, Elizaveta↗

Data-driven Stellar Models

We developed a data-driven model to map stellar parameters ( T eff , log g , and [Fe/H]) accurately and precisely to broadband stellar photometry. This model must, and does, simultaneously constrain the passband-specific dust reddening vector in the Milky Way, R. The model uses a neural network to learn the (de-reddened) absolute magnitude in one band and colors across many bands, given stellar parameters from spectroscopic surveys and parallax constraints from Gaia. To demonstrate the effectiveness of this approach, we train our model on a data set with spectroscopic parameters from LAMOST, APOGEE, and GALAH, Gaia parallaxes, and optical and near-infrared photometry from Gaia, Pan-STARRS 1, Two Micron All Sky Survey and Wide-field Infrared Survey Explorer. Testing the model on these data sets leads to an excellent fit and a precise—and by construction—accurate prediction of the color–magnitude diagrams in many bands. This flexible approach rigorously links spectroscopic and photometric surveys, and also results in an improved, T eff -dependent R. As such, it provides a simple and accurate method for predicting photometry in stellar evolutionary models. Our model will form a basis to infer stellar properties, distances, and dust extinction from photometric data, which should be of great use in 3D mapping of the Milky Way. Our trained model can be obtained at doi:10.5281/zenodo.3902382.

79 ASTRONOMY AND ASTROPHYSICS↗

On the implications of positive W mass shift

We investigate the phenomenological implications of the recent W mass measurement by the CDF collaboration, which exhibits tension with the standard model (SM) electroweak fit. Performing the fit to the electroweak observables within the SM effective field theory, we find that the new physics that contributes either to the determination of the electroweak vacuum expectation value, or to the oblique parameters, can improve the agreement with data. The best description is obtained from a fit where flavor universality is not required in the new physics operators, with 2 to 3 σ indications for several nonzero Wilson coefficients. We point out that top partners with order TeV masses could lead to the observed shift in the W mass.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Remanent magnetization and 3-dimensional density model of the Kentucky anomaly region

A three-dimensional model of the Kentucky body was developed to fit surface gravity and long wavelength aeromagnetic data. Magnetization and density parameters for the model are much like those of Mayhew et al (1982). The magnetic anomaly due to the model at satellite altitude is shown to be much too small by itself to account for the anomaly measured by Magsat. It is demonstrated that the source region for the satellite anomaly is considerably more extensive than the Kentucky body sensu stricto. The extended source region is modeled first using prismatic model sources and then using dipole array sources. Magnetization directions for the source region found by inversion of various combinations of scalar and vector data are found to be close to the main field direction, implying the lack of a strong remanent component. It is shown by simulation that in a case (such as this) where the geometry of the source is known, if a strong remanent component is present its direction is readily detectable, but by scalar data as readily as vector data.

Mayhew, M. A.↗

Computations of Unsteady Viscous Compressible Flows Using Adaptive Mesh Refinement in Curvilinear Body-fitted Grid Systems

A methodology for accurate and efficient simulation of unsteady, compressible flows is presented. The cornerstones of the methodology are a special discretization of the Navier-Stokes equations on structured body-fitted grid systems and an efficient solution-adaptive mesh refinement technique for structured grids. The discretization employs an explicit multidimensional upwind scheme for the inviscid fluxes and an implicit treatment of the viscous terms. The mesh refinement technique is based on the AMR algorithm of Berger and Colella. In this approach, cells on each level of refinement are organized into a small number of topologically rectangular blocks, each containing several thousand cells. The small number of blocks leads to small overhead in managing data, while their size and regular topology means that a high degree of optimization can be achieved on computers with vector processors.

Steinthorsson, E.↗

Photoionization of the hydrogen atom in strong magnetic fields

The photoionization of the hydrogen atom in magnetic fields B about 10 exp 11 - 10 exp 13 G typical of the surface layers of neutron stars is investigated analytically and numerically. We consider the photoionization from various tightly bound and hydrogen-like states of the atom for photons with arbitrary polarizations and wave-vector directions. It is shown that the length form of the interaction matrix elements is more appropriate in the adiabatic approximation than the velocity form, at least in the most important frequency range omega much less than omega(B), where omega(B) is the electron cyclotron frequency. Use of the length form yields nonzero cross sections for photon polarizations perpendicular to the magnetic field at omega less than omega(B); these cross sections are the ones that most strongly affect the properties of the radiation escaping from an optically thick medium, e.g., from the atmosphere of a neutron star. The results of the numerical calculations are fitted by simple analytical formulas.

Potekhin, Aleksandr IU.↗

Time-domain analysis of planar microstrip devices using a generalized Yee-algorithm based on unstructured grids

The generalized Yee-algorithm is presented for the temporal full-wave analysis of planar microstrip devices. This algorithm has the significant advantage over the traditional Yee-algorithm in that it is based on unstructured and irregular grids. The robustness of the generalized Yee-algorithm is that structures that contain curved conductors or complex three-dimensional geometries can be more accurately, and much more conveniently modeled using standard automatic grid generation techniques. This generalized Yee-algorithm is based on the the time-marching solution of the discrete form of Maxwell's equations in their integral form. To this end, the electric and magnetic fields are discretized over a dual, irregular, and unstructured grid. The primary grid is assumed to be composed of general fitted polyhedra distributed throughout the volume. The secondary grid (or dual grid) is built up of the closed polyhedra whose edges connect the centroid's of adjacent primary cells, penetrating shared faces. Faraday's law and Ampere's law are used to update the fields normal to the primary and secondary grid faces, respectively. Subsequently, a correction scheme is introduced to project the normal fields onto the grid edges. It is shown that this scheme is stable, maintains second-order accuracy, and preserves the divergenceless nature of the flux densities. Finally, for computational efficiency the algorithm is structured as a series of sparse matrix-vector multiplications. Based on this scheme, the generalized Yee-algorithm has been implemented on vector and parallel high performance computers in a highly efficient manner.

Gedney, Stephen D.↗

Modal vector estimation for closely spaced frequency modes

Techniques for obtaining improved modal vector estimates for systems with closely spaced frequency modes are discussed. In describing the dynamical behavior of a complex structure modal parameters are often analyzed: undamped natural frequency, mode shape, modal mass, modal stiffness and modal damping. From both an analytical standpoint and an experimental standpoint, identification of modal parameters is more difficult if the system has repeated frequencies or even closely spaced frequencies. The more complex the structure, the more likely it is to have closely spaced frequencies. This makes it difficult to determine valid mode shapes using single shaker test methods. By employing band selectable analysis (zoom) techniques and by employing Kennedy-Pancu circle fitting or some multiple degree of freedom (MDOF) curve fit procedure, the usefulness of the single shaker approach can be extended.

Craig, R. R., Jr.↗

High-latitude E Region Ionosphere-thermosphere Coupling: A Comparative Study Using in Situ and Incoherent Scatter Radar Observations

We present in situ and ground-based measurements of the ratio k of ion cyclotronangular frequency to ion-neutral momentum transfer collision frequency to investigateionosphere-thermosphere (IT) coupling in the auroral E region. In situ observations were obtained by NASA sounding rocket 36.234, which was launched into the nightsideE region ionosphere at 1229 UT on 19 January 2007 from Poker Flat, AK. The payload carried instrumentation to determine ion drift angle and electric field vectors. Neutral winds were measured by triangulating a chemical tracer released from rocket 41.064 launched two minutes later. k is calculated from the rotation of the ion drift angle relative to the E-cross-B drift direction in a frame co-rotating with the payload. Between the altitudes of 118 km and 130 km k increases exponentially with a scale height of 9.3 +/- 0.7 km, deviating from an exponential above 130 km. k = 1 at an altitude z(sub0) of 119.9 +/- 0.5 km. The ratio was also estimated from Poker Flat Incoherent Scatter Radar (PFISR) measurements using the rotation of ion velocity with altitude. Exponential fits to the PFISR measurements made during the flight of 41.064 yield z(sub0) 115.9 +/- 1.2 km and a scale height of 9.1 +/- 1.0 km. Differences between in situ and ground-based measurements show that the E region atmospheric densities were structured vertically and/or horizontally on scales of 1 km to 10 km. There were no signs of ionospheric structure in ion density or ion temperature below scales of 1 km. The observations demonstrate the accuracy with which the in situ and PFISR data may be used as probes of IT coupling.

Burchill, J. K.↗

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↗

Data-driven variational multiscale reduced order models

We propose a new data-driven reduced order model (ROM) framework that centers around the hierarchical structure of the variational multiscale (VMS) methodology and utilizes data to increase the ROM accuracy at a modest computational cost. The VMS methodology is a natural fit for the hierarchical structure of the ROM basis: In the first step, we use the ROM projection to separate the scales into three categories: (i) resolved large scales, (ii) resolved small scales, and (iii) unresolved scales. In the second step, we explicitly identify the VMS–ROM closure terms, i.e., the terms representing the interactions among the three types of scales. In the third step, we use available data to model the VMS–ROM closure terms. Thus, instead of phenomenological models used in VMS for standard numerical discretizations (e.g., eddy viscosity models), we utilize available data to construct new structural VMS–ROM closure models. Specifically, we build ROM operators (vectors, matrices, and tensors) that are closest to the true ROM closure terms evaluated with the available data. We test the new data-driven VMS–ROM in the numerical simulation of four test cases: (i) the 1D Burgers equation with viscosity coefficient v = 10 -3 ; (ii) a 2D flow past a circular cylinder at Reynolds numbers, Re = 100, Re = 500, and Re= 1000; (iii) the quasi-geostrophic equations at Reynolds number Re = 450 and Rossby number Ro = 0.0036; and (iv) a 2D flow over a backward facing step at Reynolds number Re = 1000. The numerical results show that the data-driven VMS–ROM is significantly more accurate than standard ROMs.

42 ENGINEERING↗

Wind dynamic range video camera

A television camera apparatus is disclosed in which bright objects are attenuated to fit within the dynamic range of the system, while dim objects are not. The apparatus receives linearly polarized light from an object scene, the light being passed by a beam splitter and focused on the output plane of a liquid crystal light valve. The light valve is oriented such that, with no excitation from the cathode ray tube, all light is rotated 90 deg and focused on the input plane of the video sensor. The light is then converted to an electrical signal, which is amplified and used to excite the CRT. The resulting image is collected and focused by a lens onto the light valve which rotates the polarization vector of the light to an extent proportional to the light intensity from the CRT. The overall effect is to selectively attenuate the image pattern focused on the sensor.

Craig, G. D.↗

The X17 with Chiral Couplings

In recent years, the ATOMKI collaboration has performed a series of measurements of excited nuclei, observing a resonant excess of electron-positron pairs at large opening angles compared to the Standard Model prediction. The excess has been hypothesized to be due to the production of a new spin-1 or spin-0 particle, X17, with a mass of about 17 MeV. Recently, the PADME experiment has reported an excess in the $e^+e^-$ cross section at center-of-mass energies near 17 MeV, perhaps further hinting at the existence of a new state. Studies of the spin-1 case have hitherto focused on either vector {\em or} axial-vector couplings to quarks and leptons, whereas UV theories more naturally produce {\em both} vector and axial-vector (\textit{i.e.} chiral) couplings, analogous to the Standard Model weak interactions. We consider the ATOMKI anomalies in the context of an $X$ with chiral couplings to quarks and explore the parameter space that can explain the ATOMKI anomalies, contrasting them with experimental constraints. We find that it is possible to accommodate the reported ATOMKI signals. However, the $99\%$ CL region is in tension with null results from searches for atomic parity violation and direct searches for new low mass physics coupled to electrons. This tension is found to be driven by the magnitude of the reported excess in the transition of $^{12}{\rm C}(17.23)$, which drives the best-fit region towards excluded couplings.

Fieg, Max H. [UC, Irvine; UC, Irvine (main); Fermi↗