Search NASA⌕ Search

SEARCH · Search NASA

Results for “inverse kinematics”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Direct reactions with the AT-TPC

Direct reactions are crucial tools for accessing properties of the atomic nucleus. Fundamental and exotic phenomena such as collective modes, pairing, weakbinding effects and evolution of single-particles energies can be investigated in peripheral collisions between a heavy nucleus and a light target. The necessity of using inverse kinematics to reveal how these structural properties change with isospin imbalance renders direct reactions a challenging technique when using the missing mass method. In this scenario, Active Target Time Projection Chambers (AT-TPC) have demonstrated an outstanding performance in enabling these types of reactions even under conditions of very low beam intensities. The AT-TPC of the Facility for Rare Isotope Beams (FRIB) is a next generation multipurpose Active Target. When operated inside a solenoidal magnet, direct reactions benefit from the measurement of the magnetic rigidity that enables particle identification and the determination of the excitation energy with high resolution without the need of auxiliary detectors. Additionally, the AT-TPC can be coupled to a magnetic spectrometer improving even further its spectroscopic investigation capability. In this contribution, we discuss inelastic scattering and transfer reaction data obtained via the AT-TPC and compare them to theory. In particular, we present the results for the 14 C(p,p′) and 12 Be (p,d) 11 Be reactions. For 14 C, we compare the experimental excitation energy of the first 1 – excited state with coupled-cluster calculationsbased on nuclear interactions from chiral effective field theory and with available shell-model predictions. For 12 Be, we determine the theoretical spectroscopic factors of the 12 Be (p,d) 11 Be transfer reaction in the shell modeland compare them to the experimental excitation spectrum from a qualitative standpoint.

active target↗
DOE OSTI2025

A Novel Approach for Computing Rigid Body Motion Using Linear Accelerations

Here, a novel approach is presented for computing general rigid body motion based on a few known linear accelerations. This method utilizes linear acceleration data obtained from three distinct points on the body, all within a body-fixed reference frame. The only requirement is that the three chosen points must not be collinear. A system of differential-algebraic equations is derived, combining principles of rigid body kinematics with theory of the rotation group SO(3). These equations provide a framework for numerically computing various motion parameters, including angular velocity, angular acceleration, body orientation, velocity field, acceleration field, and displacement field. By numerically solving this system of equations, we can fully characterize rigid body motion in three-dimensional space. A numerical example is provided to demonstrate the practical implementation and efficacy of the proposed technique, illustrating its potential for accurate motion computation in various applications.

42 ENGINEERING↗
DOE OSTI2025

Signature of 0 + excited state and shape coexistence in 94 Kr through 93 Kr(d,p) 94 Kr reaction

A measurement of the excitation spectrum of 94 Kr via one-neutron transfer to the ground state of 93 Kr using the 93 Kr(d,p) 94 Kr reaction at 8A MeV, observing the outgoing protons, performed with the IRIS facility at TRIUMF is reported. Two states in 94 Kr, at 1.50±0.14 MeV and 2.20±0.14 MeV, were observed. An adiabatic-wave approximation analysis of the differential cross sections leads us to identify the lower energy state as being populated with neutron transfer to the 3 s1/2 orbital. This leads to the first observation of the lowest 0 + excited state in 94 Kr, hence signaling shape co-existence. Theoretical calculations performed within the in-medium similarity renormalization group framework are presented that also show the existence of a low-energy 0 + state, aligning qualitatively with the observation.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Comparison of forward and inverse cam generation methods for the design of cam-linkage mechanisms

This paper presents a comparison between the forward and inverse cam generation methods applied to the design of a multi-lobed cam-rocker slider mechanism for function generation. In the forward method, the cam profile is specified directly and then the linkage output is determined using kinematics. In the inverse method, the trajectory of the output link is specified first and then the cam profile is determined using inverse linkage kinematics. First, the cam-linkage mechanism is introduced, and relevant kinematic equations are presented. Then, the cam generation methods are thoroughly explained. Here, the two cam generation methods are applied to a case study of a hydraulic motor using a multi-lobed cam and multiple linkage mechanisms. While the forward method has direct control of the cam shape, the inverse method has direct control of the linkage output and offers an advantage in calculation speed. For the mechanism selected, the cam designed with the inverse method achieved a reduction of at least 50.97% in the torque ripple when compared to the cam designed with the forward method.

42 ENGINEERING↗
DOE OSTI2022

A mixed, unified forward/inverse framework for earthquake problems: fault implementation and coseismic slip estimate

SUMMARY We introduce a new finite-element (FE) based computational framework to solve forward and inverse elastic deformation problems for earthquake faulting via the adjoint method. Based on two advanced computational libraries, FEniCS and hIPPYlib for the forward and inverse problems, respectively, this framework is flexible, transparent and easily extensible. We represent a fault discontinuity through a mixed FE elasticity formulation, which approximates the stress with higher order accuracy and exposes the prescribed slip explicitly in the variational form without using conventional split node and decomposition discrete approaches. This also allows the first order optimality condition, that is the vanishing of the gradient, to be expressed in continuous form, which leads to consistent discretizations of all field variables, including the slip. We show comparisons with the standard, pure displacement formulation and a model containing an in-plane mode II crack, whose slip is prescribed via the split node technique. We demonstrate the potential of this new computational framework by performing a linear coseismic slip inversion through adjoint-based optimization methods, without requiring computation of elastic Green’s functions. Specifically, we consider a penalized least squares formulation, which in a Bayesian setting—under the assumption of Gaussian noise and prior—reflects the negative log of the posterior distribution. The comparison of the inversion results with a standard, linear inverse theory approach based on Okada’s solutions shows analogous results. Preliminary uncertainties are estimated via eigenvalue analysis of the Hessian of the penalized least squares objective function. Our implementation is fully open-source and Jupyter notebooks to reproduce our results are provided. The extension to a fully Bayesian framework for detailed uncertainty quantification and non-linear inversions, including for heterogeneous media earthquake problems, will be analysed in a forthcoming paper.

58 GEOSCIENCES↗
DOE OSTI2024

Bayesian calibration and uncertainty quantification of a rate-dependent cohesive zone model for polymer interfaces

In this work we present a rate-dependent cohesive zone model for the fracture of polymeric interfaces and performs a Bayesian calibration, an uncertainty quantification, and a sensitivity analysis for the model. The proposed cohesive zone model accounts for both reversible elastic and irreversible rate-dependent separation sliding deformation at the interface. The viscous dissipation due to the irreversible opening at the interface is modeled using elastic-viscoplastic kinematics that incorporates the effects of strain rate. Inverse calibration of parameters for such complex models through trial and error is challenging due to the large number of parameters of the model. Moreover, the calibrated parameter values are often non-unique and uncertain when the available experimental data is limited. To tackle this challenge, we employ a Bayesian calibration approach to identify parameters from experimental data, the resulting parameters significantly enhance the accuracy of the model. To quantify the uncertainty associated with the inverse parameter estimation, a modular Bayesian approach is employed to calibrate the unknown model parameters, accounting for the parameter uncertainty of the cohesive zone model. The advantages of the Bayesian calibration over a deterministic parameter fit are demonstrated. Further, to quantify the model uncertainties, such as incorrect assumptions or missing physics, a discrepancy function is introduced, which significantly improves the model’s prediction. Finally, the total uncertainty of the model is quantified in a predictive setting. A sensitivity analysis is performed to assess how changes in the input variables of the model affect the peak load, facilitating the identification of a concise set of highly influential parameters. The present approach can be used for calibration and uncertainty quantification for other complex computational mechanics models. It should also facilitate the designing of interface materials under uncertainty.

42 ENGINEERING↗
DOE OSTI2020

Direct CP violation and the ΔI=1/2 rule in K→ππ decay from the standard model

We present a lattice QCD calculation of the Δ I = 1 / 2 , K → π π decay amplitude A 0 and ϵ ' , the measure of direct C P violation in K → π π decay, improving our 2015 calculation [1] of these quantities. Both calculations were performed with physical kinematics on a 32 3 × 64 lattice with an inverse lattice spacing of a - 1 = 1.3784 ( 68 ) GeV . However, the current calculation includes nearly 4 times the statistics and numerous technical improvements allowing us to more reliably isolate the π π ground state and more accurately relate the lattice operators to those defined in the standard model. We find Re ( A 0 ) = 2.99 ( 0.32 ) ( 0.59 ) × 10 - 7 GeV and Im ( A 0 ) = - 6.98 ( 0.62 ) ( 1.44 ) × 10 - 11 GeV , where the errors are statistical and systematic, respectively. The former agrees well with the experimental result Re ( A 0 ) = 3.3201 ( 18 ) × 10 - 7 GeV . These results for A 0 can be combined with our earlier lattice calculation of A 2 [2] to obtain Re ( ϵ ' / ϵ ) = 21.7 ( 2.6 ) ( 6.2 ) ( 5.0 ) × 10 - 4 , where the third error represents omitted isospin breaking effects, and Re ( A 0 ) / Re ( A 2 ) = 19.9 ( 2.3 ) ( 4.4 ) . The first agrees well with the experimental result of Re ( ϵ ' / ϵ ) = 16.6 ( 2.3 ) × 10 - 4 . A comparison of the second with the observed ratio Re ( A 0 ) / Re ( A 2 ) = 22.45 ( 6 ) , demonstrates the standard model origin of this “ Δ I = 1 / 2 rule” enhancement.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Imaging Polarity in Two Dimensional Materials by Breaking Friedel's Law

Friedel's law guarantees an inversion-symmetric diffraction pattern for thin, light materials where a kinematic approximation or a single-scattering model holds. Typically, breaking Friedel symmetry is ascribed to multiple scattering events within thick, non-centrosymmetric crystals. However, two-dimensional (2D) materials such as a single monolayer of MoS 2 can also violate Friedel's law, with unexpected contrast between conjugate Bragg peaks. Here, we show analytically that retaining higher order terms in the power series expansion of the scattered wavefunction can describe the anomalous contrast between hkl and $\bar{hkl}$ peaks that occurs in 2D crystals with broken in-plane inversion symmetry. These higher-order terms describe multiple scattering paths starting from the same atom in an atomically thin material. Furthermore, 2D materials containing heavy elements, such as WS 2 , always act as strong phase objects, violating Friedel's law no matter how high the energy of the incident electron beam. Experimentally, this understanding can enhance diffraction-based techniques to provide rapid imaging of polarity, twin domains, in-plane rotations, or other polar textures in 2D materials.

2D materials↗
DOE OSTI2026

Machine learning for seismic low-frequency extrapolation

The cycle-skipping problem that plagues full waveform inversion (FWI) can be at least partially mitigated if low frequencies (which encode the kinematics of wave propagation in seismic data) are recorded. However, seismic sources and receivers are band-limited, so seismic data does not generally include signals down to 0 Hz. To improve our ability to solve the seismic inverse problem, one can synthesize this missing low-frequency (LF) content from the recorded high-frequency (HF) data using machine learning (ML) models. Deep learning models such as convolutional neural networks (CNNs) demonstrate impressive ability to perform low frequency extrapolation. However, such models require powerful hardware (GPU machines) and careful training. We assess the extrapolation capabilities of three different ML models that do not require GPU machines, namely, random forest, Gaussian process regression and gradient boosting, on both synthetic and real data. Experimental results on two synthetic data sets (generated from a low velocity lens embedded in a homogeneous medium, and the Marmousi model) demonstrate that FWI applied to the extrapolated data consistently improves inversion accuracy relative to FWI applied to the original data sets that do not contain low frequencies. Application of low-frequency extrapolation to real data from the Northwest Shelf of Australia demonstrates that tree-based ML models such as gradient boosting can outperform CNNs in terms of both accuracy and computational cost on non-GPU architectures.

58 GEOSCIENCES↗
DOE OSTI2022

Neutral to charged kaon yield fluctuations in Pb – Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV

We present the first measurement of event-by-event fluctuations in the kaon sector in Pb – Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV with the ALICE detector at the LHC. The robust fluctuation correlator νdyn is used to evaluate the magnitude of fluctuations of the relative yields of neutral and charged kaons, as well as the relative yields of charged kaons, as a function of collision centrality and selected kinematic ranges. While the correlator $ν_{dyn}[K^+, K^–]$ exhibits a scaling approximately in inverse proportion of the charged particle multiplicity, $ν_{dyn}[K^0_S, K^±]$ features a significant deviation from such scaling. Within uncertainties, the value of $ν_{dyn}[K^0_S, K^±]$ is independent of the selected transverse momentum interval, while it exhibits a pseudorapidity dependence. The results are compared with HIJING, AMPT and EPOS–LHC predictions, and are further discussed in the context of the possible production of disoriented chiral condensates in central Pb – Pb collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Neutrino signal from proto-neutron star evolution: Effects of opacities from charged-current–neutrino interactions and inverse neutron decay

In this work, we investigate the impact of charged-current–neutrino processes on the formation and evolution of neutrino spectra during the deleptonization of proto-neutron stars. To this end we develop the full kinematics of these reaction rates consistent with the nuclear equation of state, including weak magnetism contributions. This allows us to systematically study the impact of inelastic contributions and weak magnetism on the $ν_e$ and $\bar{ν}_e$ luminosities and average energies. Furthermore, we explore the role of the inverse neutron decay, also known as the direct Urca process, on the emitted spectra of $\bar{ν}_e$. This process is commonly considered in the cooling scenario of cold neutron stars but has so far been neglected in the evolution of hot proto-neutron stars. We find that the inverse neutron decay becomes the dominating opacity source for low-energy $\bar{ν}_e$. Accurate three-flavor Boltzmann neutrino transport enables us to relate the magnitude of neutrino fluxes and spectra to details of the treatment of weak processes. This allows us to quantify the corresponding impact on the conditions relevant for the nucleosynthesis in the neutrino-driven wind, which is ejected from the proto-neutron star surface during the deleptonization phase.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Search for quantum gravity using astrophysical neutrino flavour with IceCube

Along their long propagation from production to detection, neutrinos undergo flavour conversions that convert their types or flavours. High-energy astrophysical neutrinos propagate unperturbed over a billion light years in vacuum and are sensitive to small effects caused by new physics. Effects of quantum gravity are expected to appear at the Planck energy scale. Such a high-energy universe would have existed only immediately after the Big Bang and is inaccessible by human technologies. On the other hand, quantum gravity effects may exist in our low-energy vacuum, but are suppressed by inverse powers of the Planck energy. Measuring the coupling of particles to such small effects is difficult via kinematic observables, but could be observable through flavour conversions. Furthermore, we report a search with the IceCube Neutrino Observatory, using astrophysical neutrino flavours to search for new space–time structure. We did not find any evidence of anomalous flavour conversion in the IceCube astrophysical neutrino flavour data. We apply the most stringent limits of any known technologies, down to 10 –42 GeV –2 with Bayes factor greater than 10 on the dimension-six operators that parameterize the space–time defects. We thus unambiguously reach the parameter space of quantum-gravity-motivated physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Refinements to the Graves–Pitarka Kinematic Rupture Generator, Including a Dynamically Consistent Slip-Rate Function, Applied to the 2019 Mw 7.1 Ridgecrest Earthquake

The main objective of this study is to develop physics–based constraints on the spatiotemporal variation of the slip–rate function using a simplified dynamic rupture model. First, we performed dynamic rupture modeling of the 2019 M w 7.1 Ridgecrest, California, earthquake, to analyze the effects of depth–dependent stress and material friction on slip rate. Then, we used our modeling results to guide refinements to the slip–rate function that were implemented in the Graves–Pitarka kinematic rupture generation technique. The dynamic ruptures were computed on a surface–rupturing, planar strike–slip fault that includes a weak (negative to low–stress–drop) zone in the upper 4 km of the crust. Below the weak zone, we placed high–stress–drop patches designed to mirror the large–slip areas seen in various rupture model inversions of the event. The locations of the high–stress–drop patches and the hypocenter were varied in multiple realizations to investigate how changing the dynamic conditions affected the resulting rupture kinematics, in particular, the slip rate. From these simulations, we observed a systematic change in the shape of the slip–rate function from Kostrov type below the weak zone to a predominantly symmetric shape within the weak zone, along with a depth–dependent reduction of peak slip rate. We generalized these shallow rupture features into a depth–dependent parametric variation of the slip–rate function and implemented it in the Graves–Pitarka kinematic rupture model generator. The performance of the updated kinematic approach was then verified in 0–4 Hz simulations of the M w 7.1 Ridgecrest earthquake, which showed that incorporating the depth–dependent variation in the shape of the slip–rate function improves the fit to the observed near–fault ground motions in the 0.5–3 s period range.

58 GEOSCIENCES↗
DOE OSTI2023

Compensating attenuation effects in full-waveform inversion with dissipation-dispersion decoupling

Seismic attenuation poses challenges to velocity model building from acoustic/elastic full-waveform inversion (FWI). In particular, when constructing the FWI gradient for velocity inversion, recent studies have indicated that the high-attenuation structure could distort the gradient by damping the amplitudes and shifting the kinematic phases, resulting in an imbalanced update and thus unreliable velocity model. These Q effects are particularly significant in the reflection acquisition geometry due to the “double-damping” issue. Here, we develop a Q-compensated FWI algorithm for constructing a Q-free FWI gradient. By using a recently developed viscoacoustic wave propagator, this compensation can be done conveniently by keeping the dispersion term and flipping the dissipation term in the wave equation when we simulate the forward and time-reversed adjoint wavefields. The resultant gradient obtained by interacting these two wavefields has correct kinematics and Q-free amplitudes. This Q-compensated FWI can balance the update between Q- and no- Q-areas, which we determine using synthetic examples. In addition, we illustrate how to take advantage of the dissipation-dispersion decoupling to determine the anomalous Q value in the Q model building workflow via Q-compensated reverse time migration.

Geochemistry & Geophysics↗
DOE OSTI2024

S -matrix positivity without Lorentz invariance: a case study

We investigate the analytic structure of scattering amplitudes in theories in which Lorentz invariance is spontaneously broken. We do so by computing and studying the S-matrix for a simple example: a superfluid described by a complex scalar with quartic interactions. The computation is confined to tree-level, for there are no absolutely stable single-particle states, though the lifetime can be made long by lowering the chemical potential. For the 2 → 2 amplitude in center-of-mass configurations, not only is crossing symmetry violated, there appears a tree level branch cut for unphysical kinematics. Its appearance is a consequence of non-analyticity in the dispersion relation. The branch point defines a new scale in the problem, which scales inversely with the chemical potential. In this example, even derivatives of the forward amplitude are positive while odd derivatives are negative. This pattern can be understood in a general way in the limit of a small chemical potential, or weak Lorentz breaking.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Effective interactions of the open bosonic string via field theory

Abstract We describe a method to extract an effective Lagrangian description for open bosonic strings, at zero transcendentality. The method relies on a particular formulation of its scattering amplitudes derived from color-kinematics duality. More precisely, starting from a (DF) 2 + YM quantum field theory, we integrate out all the massive degrees of freedom to generate an expansion in the inverse string tensionα′. We explicitly compute the Lagrangian terms through$$ \mathcal{O} $$ O (α′ 4 ), and target the sector of operators proportional toF 4 to all orders inα′.

Physics↗
DOE OSTI2026

Computing neutrino cross sections from Euclidian responses

Energy integrated neutrino cross sections are integrals of nuclear responses weighted with kinematic prefactors. We decompose the prefactors into a limited set of functions of energy transfer and show the relevant integrals are the moments of the responses, and integrals weighted with $1/(a+ω)^n$ with $n\leq 2$. These can be directly obtained from the Euclidean response, avoiding the need for inversion of the Laplace transform. As a proof of concept we study the procedure with toy-model responses for the quasielastic peak. We show that the different contributions can be straightforwardly organized in terms of relative importance, and how flux-averaged cross sections can be obtained. Using a realistic model for the response and numerical uncertainty we show that it is feasible to obtain the required integrals from the Euclidean response, with large uncertainties only for the third moment. Due to kinematic restrictions, the integrals contain contributions from the unphysical region for neutrino scattering, coming from high-momentum nucleons. We show that (in the absence of two-body currents) robust corrections for this contamination are obtained from the single-nucleon momentum distribution. These results present an opportunity to compute certain neutrino cross sections with ab-initio methods with controlled uncertainties.

Nikolakopoulos, A. [Washington U., Seattle]↗