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.

SourceResearch connectionFrom yearThrough year

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

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

NASA NTRS2003

Coulomb Stress Accumulation along the San Andreas Fault System

Stress accumulation rates along the primary segments of the San Andreas Fault system are computed using a three-dimensional (3-D) elastic half-space model with realistic fault geometry. The model is developed in the Fourier domain by solving for the response of an elastic half-space due to a point vector body force and analytically integrating the force from a locking depth to infinite depth. This approach is then applied to the San Andreas Fault system using published slip rates along 18 major fault strands of the fault zone. GPS-derived horizontal velocity measurements spanning the entire 1700 x 200 km region are then used to solve for apparent locking depth along each primary fault segment. This simple model fits remarkably well (2.43 mm/yr RMS misfit), although some discrepancies occur in the Eastern California Shear Zone. The model also predicts vertical uplift and subsidence rates that are in agreement with independent geologic and geodetic estimates. In addition, shear and normal stresses along the major fault strands are used to compute Coulomb stress accumulation rate. As a result, we find earthquake recurrence intervals along the San Andreas Fault system to be inversely proportional to Coulomb stress accumulation rate, in agreement with typical coseismic stress drops of 1 - 10 MPa. This 3-D deformation model can ultimately be extended to include both time-dependent forcing and viscoelastic response.

Smith, Bridget↗

DOE OSTI2022

Immunogenicity and protective efficacy of a rhesus adenoviral vaccine targeting conserved COVID-19 replication transcription complex

The COVID-19 pandemic marks the third coronavirus pandemic this century (SARS-CoV-1, MERS, SARS-CoV-2), emphasizing the need to identify and evaluate conserved immunogens for a pan-sarbecovirus vaccine. Here we investigate the potential utility of a T-cell vaccine strategy targeting conserved regions of the sarbecovirus proteome. We identified the most conserved regions of the sarbecovirus proteome as portions of the RNA-dependent RNA polymerase (RdRp) and Helicase proteins, both of which are part of the coronavirus replication transcription complex (RTC). Fitness constraints suggest that as SARS-CoV-2 continues to evolve these regions may better preserve cross-reactive potential of T-cell responses than Spike, Nucleocapsid, or Membrane proteins. We sought to determine if vaccine-elicited T-cell responses to the highly conserved regions of the RTC would reduce viral loads following challenge with SARS-CoV-2 in mice using a rhesus adenovirus serotype 52 (RhAd52) vector. The RhAd52.CoV.Consv vaccine generated robust cellular immunity in mice and led to significant reductions in viral loads in the nasal turbinates following challenge with a mouse-adapted SARS-CoV-2. These data suggest the potential utility of T-cell targeting of conserved regions for a pan-sarbecovirus vaccine.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2021

A diagrammatic approach for automatically deriving analytical gradients of tensor hyper-contracted electronic structure methods

In this work, we introduce a diagrammatic approach to facilitate the automatic derivation of analytical nuclear gradients for tensor hyper-contraction (THC) based electronic structure methods. The automatically derived gradients are guaranteed to have the same scaling in terms of both operation count and memory footprint as the underlying energy calculations, and the computation of a gradient is roughly three times as costly as the underlying energy. The new diagrammatic approach enables the first cubic scaling implementation of nuclear derivatives for THC tensors fitted in molecular orbital basis (MO-THC). Furthermore, application of this new approach to THC-MP2 analytical gradients leads to an implementation, which is at least four times faster than the previously reported, manually derived implementation. Finally, we apply the new approach to the 14 tensor contraction patterns appearing in the supporting subspace formulation of multireference perturbation theory, laying the foundation for developments of analytical nuclear gradients and nonadiabatic coupling vectors for multi-state CASPT2.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Ab-initio simulation of spin-vibronic spectra of methoxy radical

Despite the fact that experimental and theoretical work on the spectrum of methoxy has stretched from the microwave to the ultraviolet and proceeded for nearly 50 years, parts of the spectrum have remained a challenge to simulate theoretically and make reliable line-by-line assignments. The spectral complexity arises because the radical has a non-zero electron spin and significant vibronic coupling between the two elec- tronic components of the ground state due to the presence of a conical intersection. This work describes a completely ab initio effort to understand and assign the spin- vibronic levels of the X 2E state from 0 to above 3000 cm−1, a region that includes the fundamental transitions of the C-H symmetric and asymmetric stretches that have not previously been identified uniquely. A potential energy surface for methoxy was calculated at the EOM-CCSDT/ANO1 level of theory. Subsequently this potential energy surface was fit to a quartic power series expansion of all nine vibrational nor- mal coordinates (as determined at the minimum of the conical intersection) by the use of a machine-learning-based algorithm. After the addition of spin-orbit coupling, the spin-vibronic problem was solved using both the Krylov-Schur and Lanczos algorithms with the SOCJT3 software to converge eigenvalues up to 3500 cm−1 and their eigen- vectors. The latter were used, in conjunction with the calculated dipole moment and its derivatives (calculated using finite differences at EOM-CCSDT/ANO1 level), to determine spectral intensities for the spin-vibronic spectra. The calculated transition frequencies and intensities were used to simulate and assign the observed transitions of the spin-vibronic spectra of the radical. The credibility of the assignments and their significance is discussed in detail.

Sharma, Ketan [University of Florida, Gainesville,↗

NASA NTRS1978

A time dependent model of the Jovian current sheet

Analysis of vector helium magnetometer measurements during the Pioneer 10 outbound pass through the Jovian magnetosphere reveals that the average location of the near-equatorial current sheet is a surface whose shape is determined by the velocity at which the magnetic perturbations produced by a rotating tilted dipole propagate to large radial distances. A model is presented which is a surface encircling the planet and rotating rigidly with it. Based on a linear fit to the times of the current sheet crossings, a formula is derived for the surface representing the average position of the near-equatorial current sheet in the Jovian magnetosphere.

Kivelson, M. G.↗

NASA NTRS1985

Magnetization models for the source of the 'Kentucky anomaly' observed by Magsat

Both the aeromagnetic data and magnetic anomaly data obtained by Magsat indicate the presence of a very magnetic source region within the crust beneath Kentucky and Tennessee. A source model was previously developed to fit surface gravity and long-wavelength aeromagnetic data, using limited seismic constraint. For the present study the model was further developed, and 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 using both prismatic model sources and 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 determinable, but by scalar data as readily as vector data. Magnetization magnitude for the extended source region is about 3 A/m if the vertical extent of the source includes the whole of the crust.

Mayhew, M. A.↗

DOE OSTI2025

Efficient lattice QCD computation of radiative-leptonic-decay form factors at multiple positive and negative photon virtualities

In previous work [D. Giusti, Methods for high-precision determinations of radiative-leptonic decay form factors using lattice QCD, Phys. Rev. D 107, 074507 (2023)], we showed that form factors for radiative leptonic decays of pseudoscalar mesons can be determined efficiently and with high precision from lattice QCD using the “three-dimensional (3D) method,” in which three-point functions are computed for all values of the current insertion time and the time integral is performed at the data-analysis stage. Here, we demonstrate another benefit of the 3D method: the form factors can be extracted for any number of nonzero photon virtualities from the same three-point functions at no extra cost. We present results for the $D_s → ℓνγ*$ vector form factor as a function of photon energy and photon virtuality, for both positive and negative virtuality, for a single ensemble with 340 MeV pion mass and 0.11 fm lattice spacing. In our analysis, we separately consider the two different time orderings and the different quark flavors in the electromagnetic current. We discuss in detail the behavior of the unwanted exponentials contributing to the three-point functions, as well as the choice of fit models and fit ranges used to remove them for various values of the virtuality. While positive photon virtuality is relevant for decays to multiple charged leptons, negative photon virtuality suppresses soft contributions and is of interest in QCD-factorization studies of the form factors.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Dynamical model of J/ψ photoproduction on the nucleon

Here, a dynamical model based on a phenomenological charm quark-nucleon (c – N) potential v cN and the Pomeronexchange mechanism is constructed to investigate the J/ψ photoproduction on the nucleon from threshold to invariant mass W = 300 GeV. The J/ψ – N potential, V J/ψN (r), is constructed by folding v cN into the wave function φ J/ψ ($c\bar{c}$) of J/ψ within a constituent quark model (CQM) of Segovia et al. [Int. J. Mod. Phys. E 22, 1330026 (2013)]. A photoproduction amplitude is also generated by v cN by a $c\bar{c}$–loop integration over the γ → $c\bar{c}$ vertex function and φ J/ψ ($c\bar{c}$). No commonly used vector meson dominance assumption is used to define this photoproduction amplitude which is needed to describe the data near the threshold. The c – N potential v cN (r) is parameterized in a form such that the predicted V J/ψN (r) at large distances has the same Yukawa potential form extracted from a lattice QCD (LQCD) calculation of Kawanai and Sasaki, [Phys. Rev. D 82, 091501(R) (2010)]. The parameters of vcN are determined by fitting the total cross-section data of Jefferson Laboratory (JLab) by performing calculations that include J/ψ – N final-state interactions (FSI). The resulting differential cross sections dσ /dt are found in good agreements with the data. It is shown that the FSI effects dominate the cross section in the very near-threshold region, allowing for sensitive testing of the predicted J/ψ – N scattering amplitudes. By imposing the constraints of J/ψ – N potential extracted from the LQCD calculation of Kawanai and Sasaki, [Phys. Rev. D 82, 091501(R) (2010)], we have obtained three J/ψ – N potentials which fit the JLab data equally well. The resulting J/ψ – N scattering lengths are in the range of a = [-0.05, -0.25] fm. With the determined v cN (r) and the wave functions generated from the same CQM, the constructed model is used to predict the cross sections of photoproduction of η c (1S) and ψ(2S) mesons for future experimental tests.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1995

Ultraviolet spectropolarimetry of high-redshift quasars with the Hubble Space Telescope

Ultraviolet spectropolarimetry of three bright high-redshift low-polarization quasars (LPQs) was obtained with the Faint Object Spectrograph of the Hubble Space Telescope (HST). Two of the quasars, PG 1634+706 and PG 2302+029, had polarizations p approximately = 0.5%-1.0% throughout the ultraviolet, and showed no significant variation of polarization amplitude or position angle with wavelength. PG 2302+029 was also marginally (2.4 sigma) circularly polarized in the optical continuum. For the highest redshift quasar, PG 1222+228 (Ton 1530), the polarization was measured down to rest wavelengths below 800 A. Although the continuum of PG 1222+228 was weakened by Lyman limit absorption from an intergalactic gas cloud, the polarization increased sharply from 1% to about 4.5%, a change of 4 sigma significance. This abrupt rise in polarization does not appear attributable to any known instrumental artifact. These UV polarizations were only slightly less than those previously observed for these same objects in the optical. The polarization spectra were flat with a typical slope of the polarized flux pF(sub nu) proportional to nu(exp -0.8 +/- 0.5). Unlike the case of several high luminosity Seyfert 1 nuclei studied previously, polarization caused by scattering from dust grains does not provide the best fit to the polarization spectra of these luminous quasars. These observed spectra are consistent with a wavelength-independent polarization proportional to the total nonstellar light or, possibly, to the contribution of the blue thermal component. The polarization spectra have insufficient signal-to-noise to locate the scatterers with respect to the continuum source and the much larger broad line region. A decrease in amplitude and rotation of the position angle of the polarization vector at the shortest wavelengths, which could result from general relativistic effects near a spinning black hole, was not observed. In fact, in PG 1222+228, the polarization was observed to increase at the shortest wavelengths. The rise in polarization with frequency is so sharp that it cannot be due to any wavelength-independent polarizing mechanism at any radius in an accretion disk. Such a rise could be attributable, for example, to a relative increase in scattering opacity over absorption at higher frequencies.

Impey, C. D.↗

NASA NTRS2025

Improving Adhesive Bondline Time of Flight Predictions During Autoclave Cure Utilizing Machine Learning

Composite materials are increasingly being used in aerospace applications due to their superior strength-to-weight ratio compared to commonly used metals. A current limitation to widespread adoption is the certification of adhesively bonded joints. One approach to improving adhesive bonding in composites is accurately measuring the thickness of adhesive bondlines in composite laminates. Precise bondline thickness control is essential for aerospace applications where adhesive layer thickness directly affects joint fracture properties and structural performance. This study focused on implementing machine learning techniques to determine the ultrasonic time of flight (directly correlated to thickness) in adhesive bondlines throughout autoclave cure cycles. A high-temperature (use up to 180°C) ultrasonic scanning system was deployed in an autoclave to provide time of flight data through composite panels. Three experiments were conducted on the curing of 305 mm × 305 mm unidirectional composite panels. In the first experiment, a piecewise function was fit for the temperature correction factor to account for changing autoclave temperatures. Due to deficiencies in the first calibration experiment, a second experiment was run, and the results were used to train a machine learning model. The revised experiment, in combination with the machine learning model, significantly increased the accuracy of the bondline time of flight predictions (~14% error reduced to <1%). Data was processed using the Regression Learner Application in MATLAB®, with a Support Vector Machine selected for the model. The result was a machine learning algorithm capable of reliably quantifying ultrasonic time of flight through adhesive bondlines. The third experiment provided independent test data for the machine learning model, demonstrating that the model produces accurate predictions from data beyond its training set.

Machine Learning↗

DOE OSTI2023

Dark photon and dark Z mediated B meson decays

We study flavor changing neutral current decays of B and K mesons in the dark U(1) D model, with the dark photon/dark Z mass between 10 MeV and 2 GeV. Although the model provides an improved fit (compared to the standard model) to the differential decay distributions of B → K ($\ast$) ℓ + ℓ – , with ℓ = μ, e, and B s → $\phi$μ + μ – , the allowed parameter space is ruled out by measurements of atomic parity violation, K + → μ + + invisible decay, and ${B}_s-{\overline{B}}_s$ mixing, among others. To evade constraints from low energy data, we extend the model to allow for (1) additional invisible Z D decay, (2) a direct vector coupling of Z D to muons, and (3) a direct coupling of Z D to both muons and electrons, with the electron coupling fine-tuned to cancel the Z D coupling to electrons via mixing. We find that only the latter case survives all constraints.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Search for leptoquark pair production decaying into $te^- \bar{t}e^+$ or $t\mu ^- \bar{t}\mu ^+$ in multi-lepton final states in pp collisions at $\sqrt{s} = 13\,\textrm{TeV}$ with the ATLAS detector

A search for leptoquark pair production decaying into $te^- \bar{t}e^+$ or $t\mu ^- \bar{t}\mu ^+$ in final states with multiple leptons is presented. The search is based on a dataset of pp collisions at $\sqrt{s}=13~\text {TeV}$ recorded with the ATLAS detector during Run 2 of the Large Hadron Collider, corresponding to an integrated luminosity of 139 fb -1 . Four signal regions, with the requirement of at least three light leptons (electron or muon) and at least two jets out of which at least one jet is identified as coming from a b-hadron, are considered based on the number of leptons of a given flavour. The main background processes are estimated using dedicated control regions in a simultaneous fit with the signal regions to data. No excess above the Standard Model background prediction is observed and 95% confidence level limits on the production cross section times branching ratio are derived as a function of the leptoquark mass. Under the assumption of exclusive decays into $te^{-}$ ($t\mu ^{-}$ ), the corresponding lower limit on the scalar mixed-generation leptoquark mass $m_{\textrm{LQ}_{\textrm{mix}}^{\textrm{d}}}$ is at 1.58 (1.59) TeV and on the vector leptoquark mass $m_{{\tilde{U}}_1}$ at 1.67 (1.67) TeV in the minimal coupling scenario and at 1.95 (1.95) TeV in the Yang–Mills scenario.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS2007

Aeroelastic Model Structure Computation for Envelope Expansion

Structure detection is a procedure for selecting a subset of candidate terms, from a full model description, that best describes the observed output. This is a necessary procedure to compute an efficient system description which may afford greater insight into the functionality of the system or a simpler controller design. Structure computation as a tool for black-box modelling may be of critical importance in the development of robust, parsimonious models for the flight-test community. Moreover, this approach may lead to efficient strategies for rapid envelope expansion which may save significant development time and costs. In this study, a least absolute shrinkage and selection operator (LASSO) technique is investigated for computing efficient model descriptions of nonlinear aeroelastic systems. The LASSO minimises the residual sum of squares by the addition of an l(sub 1) penalty term on the parameter vector of the traditional 2 minimisation problem. Its use for structure detection is a natural extension of this constrained minimisation approach to pseudolinear regression problems which produces some model parameters that are exactly zero and, therefore, yields a parsimonious system description. Applicability of this technique for model structure computation for the F/A-18 Active Aeroelastic Wing using flight test data is shown for several flight conditions (Mach numbers) by identifying a parsimonious system description with a high percent fit for cross-validated data.

Kukreja, Sunil L.↗

DOE OSTI2025

Emergent topological polarization textures in relaxor ferroelectrics

Relaxor ferroelectrics underpin high-performance actuators and sensors, yet the nature of polar heterogeneities driving their broadband dielectric response remains debated. Using a unified, multimodal structural refinement framework— simultaneously fitting complementary X-ray and neutron total scattering, X-ray absorption spectra, and diffuse scattering—we reconstruct 3D mesoscale polarization maps in the classic relaxor system PbMg 1/3 Nb 2/3 O 3 –PbTiO 3 . We uncover self-organized swirling polarization textures with half-skyrmion (meron) vortices, challenging models of independent polar nanoregions. These textures, characterized by smooth changes in the polarization direction, originate from overlapping volumes in which the projections of locally correlated polarization vectors onto each volume’s long axis share the same sign. Vortex cores correlate strongly with local charge and strain gradients imposed by compositional heterogeneities. In this work, our results suggest that chemical disorder, acting via depolarizing and strain fields, stabilizes topological vortex textures of the polarization field, offering a route for engineering new dielectric and ferroelectric functionalities.

Eremenko, Maksim [National Institute of Standards ↗

NASA NTRS1988

Application of advanced grid generation techniques for flow field computations about complex configurations

In the computation of flowfields about complex configurations, it is very difficult to construct a boundary-fitted coordinate system. An alternative approach is to use several grids at once, each of which is generated independently. This procedure is called the multiple grids or zonal grids approach; its applications are investigated. The method conservative providing conservation of fluxes at grid interfaces. The Euler equations are solved numerically on such grids for various configurations. The numerical scheme used is the finite-volume technique with a three-stage Runge-Kutta time integration. The code is vectorized and programmed to run on the CDC VPS-32 computer. Steady state solutions of the Euler equations are presented and discussed. The solutions include: low speed flow over a sphere, high speed flow over a slender body, supersonic flow through a duct, and supersonic internal/external flow interaction for an aircraft configuration at various angles of attack. The results demonstrate that the multiple grids approach along with the conservative interfacing is capable of computing the flows about the complex configurations where the use of a single grid system is not possible.

Kathong, Monchai↗

DOE OSTI2023

Analysis of Three Body Decays in Quasi-Real Photoproduction

This thesis presents preliminary measurements of the moments of angular distribution for the K?K+ mesonic final state using a quasi-real, linearly polarised photon. Two of the main points of focus for this analysis were the development of the formalism for photoproduced vector-pseudoscalar decay and the tools by which the analysis was carried out. The formalism described how the moments of angular distribution could be extracted from the angular decay variables, and how they related to resonance decays in terms of partial waves. In extracting the moments, it was found that adapted Markov Chain Monte Carlo (MCMC) methods proved to be more effective in terms of extracting results from the data compared to the sole use of gradient descent based fitting algorithms. These measurements are complementary to other decay channels and production mechanisms, with the global aim of probing the existence of exotic mesons. Data from this analysis was taken by the CLAS collaboration at the Thomas Jefferson National Accelerator Facility (JLab) by scattering from a liquid hydrogen target using a 10.6 GeV electron. This work was done within the MesonEx program at CLAS12 that aims to map the spectrum of mesons, and in doing so, gain a greater understanding of Quantum ChromoDynamics (QCD) and the strong force.

Wishart, Robert↗

DOE OSTI2025

Light and Strange Vector Resonances from Lattice QCD at Physical Quark Masses

We present the first ab initio calculation at physical quark masses of scattering amplitudes describing the lightest pseudoscalar mesons interacting via the strong force in the vector channel. Using lattice quantum chromodynamics, we postdict the defining parameters for two short-lived resonances, the ρ(770) and K*(892) which manifest as complex energy poles in ππ and Kπ scattering amplitudes, respectively. The calculation proceeds by first computing the finite-volume energy spectrum of the two-hadron systems and then determining the amplitudes from the energies using the Lüscher formalism. The error budget includes a data-driven systematic error, obtained by scanning possible fit ranges and fit models to extract the spectrum from Euclidean correlators, as well as the scattering amplitudes from the latter. The final results, obtained by analytically continuing multiple parametrizations into the complex energy plane, are M ρ = 796(5)(50)MeV, Γ ρ = 192(10)(31) MeV, M K* = 893(2)(54) MeV, and Γ K = 51(2)(11) MeV, where the subscript indicates the resonance and M and Γ stand for the mass and width, respectively, and where the first bracket indicates the statistical and the second bracket the systematic uncertainty.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1975

Neutral hydrogen flux measured at 100- to 200-km altitude in an electron aurora

Neutral hydrogen fluxes were measured at altitudes of 120-200 km by a rocket payload that also measured electron and proton fluxes and vector magnetic fields. An intense electron arc was crossed, while an upper limit to the flux of 0.5- to 20-keV protons was 1,000,000 per sq cm s sr keV. A neutral flux of 50,000,000 per sq cm s sr was observed, assuming hydrogen with greater than 1-keV energy, with greater north-south extent than the electron flux. Its pitch angle distribution was peaked toward 90 deg, tending toward isotropy in the center. This is fitted to a model describing spreading of an initial proton arc above 500 km.

Iglesias, G. E.↗