Search NASA⌕ Search

SEARCH · Search NASA

Results for “reference frames”

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

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

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
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

High-Resolution Sampling of a River Plume Front with Uncrewed Underwater and Aerial Vehicles

Sampling fast-propagating oceanic features is inherently challenging and demands versatile instrumentation and innovative strategies. This paper introduces a novel sampling strategy designed to capture such phenomena, exemplified by a river plume front. Our method revolves around modifying the preprogrammed pathway of an uncrewed underwater vehicle (UUV) to dynamically track and three-dimensionally sample the evolution of the front. To enable the UUV to follow the feature, we adapt the use of a drifting gateway buoy to be positioned and trapped at the front’s convergence zone, allowing underway navigation relative to the buoy. In our demonstration, we showcase the effectiveness of this strategy by successfully conducting over 30 crossings of a river plume front within a 6-h window. The UUV sensors allowed a comprehensive assessment of key front characteristics, including density, velocity, and turbulence. Supplemental drone footage contributed to the overall picture and facilitated the transformation of the dataset into a front-following reference frame. This article provides an in-depth description of the deployment strategy and required postcollection data processing, including frontal crossing detection, the assessment of the frontal orientation from drone footage, and defining the plume bottom boundaries using backscatter intensity contours.

autonomous observations↗
DOE OSTI2025

SSAPy - Space Situational Awareness for Python

SSAPy is a fast and flexible orbit modeling and analysis tool for orbits spanning from low-Earth into the cislunar regime. Orbits can be flexibly specified from common input formats such as Keplerian elements or two-line element (TLE) data files. SSAPy allows users to model satellites and specify parameters such as satellite area, mass, and drag coefficients. SSAPy includes a customizable force-propagation with a range of Earth, Lunar, radiation, atmospheric, and maneuvering models. SSAPy makes use of various community integration methods and can calculate time-evolved orbital quantities, including satellite magnitudes and state vectors. Users can specify various space- and ground-based observation models with support for multiple coordinate and reference frames. SSAPy also supports orbit analysis and propagation methods such as multiple hypothesis tracking and has built-in uncertainty quantification. The majority of SSAPy’s methods are vectorized and parallelizable, allowing for effective use of high-performance computer (HPC) systems. Finally, SSAPy has plotting functionality, allowing users to visualize orbits and trajectories. Examples are shown in Figure 1 and Figure 2.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

High-Speed and High-Quality Field Welding Repair Based on Advanced Non-Destructive Evaluation and Numerical Modeling

Creep strength-enhanced ferritic (CSEF) steels such as Grade 91 (9Cr-1Mo-V) and Grade 92 (Fe-9Cr-2W-0.5Mo) steels are widely used in the fossil-fuel-fired and nuclear power plants. The weld integrity of these steels is crucial for power plants' safe and reliable operations. Due to harsh service conditions, the steel weld can become susceptible to environmental degradation. Field welding repair is used to restore the degraded weld’s performance where a controlled temper-bead welding technique is commonly used to temper the freshly formed martensite during welding. However, knowledge of weld repairability is limited and experimental trial and error optimization to achieve desired microstructure and joint properties is expensive and time-consuming. Many existing computational models, e.g., finite element models, are limited to solving heat conduction equation and ignoring convective heat transfer due to molten metal flow. These models can result in over-prediction of peak temperatures of weld pool and heat-affected zone (HAZ), which in turn can affect the accuracy of tempering prediction. Moreover, these finite element models require an input of the deposit profiles in advance and thus limits the usability of these models. Here, a molten pool-based, multi-pass multi-layer model has been developed based on computational fluid dynamics (CFD) approach with the Volume of Fluid (VOF) method. The model calculates the bead formation, thereby eliminating the need for pre-determined bead profiles required by finite element models. For computational efficiency, a coordinate system attached to the moving heat source is utilized. A subroutine is developed to convert the temperature profiles in the reference frame stationary to the heat source to that stationary to the workpiece. The converted thermal cycles are then imported into a microstructure model to compute the tempering kinetics and resultant hardness using a Johnson-Mehl-Avrami-Kolmogorov (JMAK), and modified Grange-Baughman parameter. The modeling approach is first developed and validated on single- and multi-pass deposition of stainless steel filler metal onto a SA-533 high strength steel substrate. The models are then applied to a multi-pass V-groove repair weld of Grade 91 steel plate as well as directed energy deposition of Grade 92 steel. Non-destructive characterization of microstructures was performed on Grade 91 and 92 steel welds. Two welding processes, cold metal transfer (CMT) and flux-cored arc welding (FCAW), were investigated for the Grade 91 steel weld samples. For the Grade 92 weld samples, three different heat inputs (low, medium, and high) of gas tungsten arc welding (GTAW) were utilized to replicate traditional field welding processes. The non-destructive evaluation (NDE) method used for this research was immersion ultrasonic testing (UT) using a micro-resolution ultrasonic imaging methodology specifically designed to operate in the through-transmission configuration operating at 20 MHz of frequency. The system used a focused ultrasonic beam spot size diameter between 250-300 μm, and a 6 μm laser vibrometer spot size for detection, to produce highly defined images with longitudinal and mode-converted shear waves. From the micro-resolution ultrasonic C-scan images, three microstructural regions, i.e., weld metal (WM), HAZ, and base metal (BM), were clearly identifiable. Various levels of ultrasonic amplitudes distributed over the three regions were correlated with electron beam backscattered diffraction (EBSD) images using grain size, grain boundaries, and dislocation densities. The results showed that areas with relatively higher ultrasonic amplitude levels were associated with smaller grains and higher dislocation densities, while areas with lower amplitude levels were associated with larger grains and lower dislocation densities. In addition, ultrasonic velocity data obtained across the three different weld microstructural regions of Grade 91 test samples were correlated with optical metallographic images and hardness measurements. The results showed distinctive decreases in ultrasonic velocity and hardness over the HAZ region, where weld failures often occur during service.

36 MATERIALS SCIENCE↗
DOE OSTI2025

DROP DURABILITY ASSESSMENT OF ELECTRONIC ASSEMBLIES UNDER OFF-AXIS LOADING WITH SKEWED FIXTURES

This thesis studies drop durability of electronic assemblies when the acceleration vector is oriented at 45° to the out-of-plane direction of the circuit card. The off-axis drop tests are accomplished with a skewed fixture and are conducted as a proxy for multiaxial drop testing. Advanced shock testing and vibration test methods have been developed over the last few decades to better represent real-world field environments during ground-based laboratory testing. However, many of these test methods require expensive and specialized equipment not available in most laboratories. An alternative approach for approximating simultaneous loading along multiple axes on conventional equipment utilizes skewed fixtures which have seen use in off-axis random vibration and drop impact testing. These methods generally rely on the conversion of a uniaxial input load from the test equipment (using a uniaxial drop tower or shaker) into a multiaxial load when resolved in the reference frame of the test article (mounted on a skewed fixture). Skewed fixture design is presented and recommendations for conducting skewed angle drop testing are introduced based on local measurements along the skewed face of the fixture to accurately monitor the impact event. Characterization tests were performed with a skewed fixture, at simultaneous acceleration loads from 500 to 3,000 g in two (in-plane and out-of-plane) directions, while meeting standard time domain tolerances. Upon experimental characterization, drop shock durability tests were conducted on a printed circuit assembly (PCA). Mean drops-to-failure were measured and quantified with Weibull statistics. Dominant solder joint failure modes were identified via failure analysis. Prior work on inclined angle impact testing is limited, and the majority of solder joint interconnect level fatigue studies are conducted considering perpendicular loading normal the circuit card. Low-cycle fatigue curves are generated based on plastic strain and plastic work density within the solder joint. A multiscale nonlinear finite element model is used to relate board-level flexure to solder joint interconnect level plastic strain. A high strain rate solder constitutive model allows for accurate modeling of solder plasticity resulting from high-impact drop shock. Fatigue parameters are computed from the Coffin-Manson relation and Palmgren-Miner damage accumulation. This work serves to apply established low-cycle fatigue methods for conventional drop shock loading (impact normal to circuit card) to non-perpendicular loading with a skewed fixture.

Hower, Jonathan [Kansas City National Security Cam↗
DOE OSTI2025

Enabling Grid-Forming Control Under Unbalanced Conditions

Standalone microgrids often experience unbalanced loading and faults, which can cause grid-forming control designed for balanced conditions to produce oscillatory responses. To address this issue, a compact time-domain transformation appropriate for inverter control is proposed, allowing the conversion of unbalanced three-phase signals to positive and negative synchronous reference frames. This transformation supports the development of a grid-forming control with fault ride-through, featuring frequency and voltage droop controllers and nested current and voltage control loops that seamlessly integrate an enhanced current limiter. The effectiveness of the proposed control and transformation is demonstrated through analytical results and electromagnetic transient simulation.

24 - POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

An Optimized Parameterization of Sub‐Grid Scale Advection for Convection Permitting Models

Convection‐permitting models (CPMs) explicitly resolve deep convection yet under‐resolve the organized lateral exchanges among drafts and their environment that control entrainment/detrainment, precipitation efficiency, and mesoscale structure. In this work, we introduce the Optimized Advection Scheme (OAS), which introduces a small rotation of the Cartesian frame of reference for the horizontal winds relative to other variables used in advection that induces cross‐gradient transport to mimic under‐resolved convective mixing. The rotation angle is selected to minimize the Kullback–Leibler divergence between the simulated and satellite observed precipitation intensity distributions, yielding a physically consistent perturbation that is computationally inexpensive and portable. Optimized Advection Scheme is implemented in WRF and evaluated over Amazon (April 2014). It shifts precipitation–precipitable‐water joint distributions toward lighter rain, reduces overly intense rates, and improves mesoscale convective system (MCS) lifetime and propagation. Mechanistically, the added cross‐gradient transport promotes convective detrainment and environmental mixing, which cools and moistens the mid‐troposphere, weakens downward momentum transport, alleviates excessive downwelling shortwave biases, and warms the surface temperature. The optimized rotation angle yields comparable improvements at 4‐km and 1‐km grid spacing, demonstrating resolution‐independent benefits across the CPM gray zone. By targeting the dynamical root of under‐mixed convective circulations, rather than tuning model microphysics or closures, OAS delivers robust, scale‐aware improvements in precipitation statistics, cloud vertical structure, and characteristics of MCS (MCSs), offering a practical pathway to more reliable CPM simulations for weather and climate applications.

CPM↗
DOE OSTI2024

Explicit simulation of the Brownian rotation of arbitrary shaped aerosol particles using quaternions

The shape of an aerosol particle strongly influences its mass and momentum transfer cross-sections, charging properties, and other physical properties. Here, we present an explicit time-stepping procedure to simulate the rotational Brownian motion of arbitrary shaped aerosol particles by solving Euler’s equation of rotation. A Langevin formulation of the rotation equations is used, wherein Brownian motion due to thermal collisions between a particle and background gas molecules is represented using a stochastic fluctuating torque and fluid resistance is included as a drag torque. To avoid singularities associated with describing the orientation of a shape with Euler angles, we employ a quaternion formulation that leads to first-order stochastic differential equations to describe the evolution of the angular position and angular velocity of a rigid body. We perform all the rotational dynamics calculations in the body-fixed frame of reference attached to the rotating shape whose basis vectors are the normalized eigenvectors of the inertia tensor of the particle. Numerical solutions to rotation under torque-free conditions, damped rotation without Brownian motion, and stochastic rotation for arbitrary shapes are presented and discussed. The presented method enables time-resolved simulation of Brownian rotation for direct comparison with experimentally measured trajectories or statistical measures. The second order accuracy of the used time-stepping procedure places a severe restriction on the timestep that can be used for obtaining accurate results. Animations of presented simulations are included for visualizing rotational motion at various gas pressures. To aid implementation, MATLAB ® codes are also provided. Extension to include translation Brownian motion is straightforward.

Roy, Mrittika↗
DOE OSTI2025

Generalized parton distributions from lattice QCD with asymmetric momentum transfer: Unpolarized quarks at nonzero skewness

We extend the formalism of asymmetric frames of reference for generalized parton distributions (GPDs) to the case of nonzero skewness, i.e., including longitudinal momentum transfer. The framework, based on Lorentz-invariant amplitudes and previously developed and numerically implemented for unpolarized, helicity and transversity GPDs at zero skewness, gives efficient access to a broad range of kinematics, making full mapping of GPDs from the lattice realistic. The general-skewness formalism is tested using lattice data with both transverse and longitudinal or only longitudinal momentum transfer, the latter being a special case with a reduced number of independent amplitudes. We extract the amplitudes in coordinate space and express the GPDs 𝐻 and 𝐸 in terms of these amplitudes. This is followed by reconstruction of quasidistributions and their matching to the light cone. We further identify and discuss the principal challenges for nonzero skewness GPDs.

Lattice QCD↗
DOE OSTI2026

Comprehensive framework for assessing and optimizing existing research networks

Conservation, monitoring, and research networks, or collections of ecological research sites unified under a common mission of data collection or a research mission, are essential infrastructure for understanding large landscapes. However, most networks developed opportunistically over decades rather than through systematic design, creating potential limitations in the ability to address conservation challenges across entire regions. We developed a framework to evaluate how well an existing research network represents the environmental conditions its members study and devised an approach to rank sites of priority for strategic expansion. Our approach measures performance through environmental representativeness, geographic coverage, and adequacy for scientific inference and thus optimizes limited monitoring resources to maximize scientific impact. We demonstrated this approach with the U.S. Department of Agriculture (USDA) Forest Service Experimental Forests and Ranges Network (EFRN), a 79‐site network across the United States that grew opportunistically over a century. At the national scale, the network effectively captured high‐biomass forests important for carbon cycle research; 82% of forest biomass was in well‐represented areas. Some areas in Texas, Florida, the Rocky Mountains, and the West Coast had no relevant EFRN sites, which limits the ability to make regional inferences. A fundamental challenge for the EFRN was that sites improving regional extent coverage sometimes provided minimal national benefits, which can create conflicts between local and global priorities. Adding the highest‐ranked candidate site provided a relevant site for 17% of currently poorly represented 1‐km pixel cells nationally, but regional and national site rankings varied considerably due to nested spatial inference. This framework provides quantitative tools for strategic infrastructure decision‐making, ensures that limited monitoring resources maximize conservation impact, and can be applied broadly to address the widespread challenge of optimizing conservation and monitoring networks worldwide.

additional site↗
DOE OSTI2026

Finite deformation implementation of a mixed-mode single-integral type cohesive zone with reorienting surfaces of separation

To model material ductile failure and crack propagation, cohesive zone elements can be embedded along potential fracture paths in a finite element simulation. When damage criteria are met, elements in the mesh decohere, simulating the formation and propagation of a crack. In this paper, we present a novel computational algorithm based on finite deformation theory, essential to modeling crack initiation and growth in solids undergoing large deformations. This new algorithm was formulated within a Lagrangian frame of reference to extend previous cohesive zone algorithms to include modeling crack growth in finite deformation contexts. The local coordinate system, necessary for defining an embedded cohesive zone, is constructed based upon the current configuration and is updated within the nonlinear iteration process, thereby resulting in the convergence of the solution for a growing crack in a large deformation quasi-static setting. The model’s accuracy was demonstrated by comparing finite element model simulation results with the analytic case of a constant surface separation, as shown in the verification examples. The power and efficacy of the algorithm to capture large deformations during crack growth were then demonstrated with a double cantilever beam example case. It indicates that the model can be applied to a variety of physical circumstances for predicting crack initiation and growth with delamination and fracture.

42 ENGINEERING↗
DOE OSTI2026

From Raw to Curated Data: A Lakehouse Approach for Scientific Workflows

This report provides a technical overview of how to go from raw to curated data in three stages using a lakehouse approach. We focus on the application of open source tools in scientific use cases (while noting parallels to enterprise and commercial alternatives). Our goal is to provide scientific data managers and infrastructure providers with a common frame of reference for understanding and applying modern lakehouse technologies and approaches.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2026

A State-Space Model for Stability Boundary Analysis of Grid-Following Voltage Source Converters Considering Grid Conditions

With the growing significance of renewable energy resources and energy storage systems, the number of grid-connected inverters has been rising at an increasingly rapid pace. Generally, these inverters are directly integrated with the distribution network by synchronizing with the grid voltage at the point of common coupling. However, the low grid strength and varying R/X ratios, as the common characteristics of most distribution networks or weak grids, can lead to dynamic interactions that comprise stability and limit the power transfer capacity of grid-connected inverters. To ensure stable operation of the inverters, researchers must determine the stability boundary, described as the maximum power transfer capacity of grid-connected inverters under the premise of maintaining system small-signal stability. For this purpose, we propose to formulate a state-space model of the system in the synchronously rotating dq-frame of reference and perform eigenvalue analysis to determine the stability boundary. With a detailed model of the control structure and parameters of the grid-connected inverters, the stability boundary is identified as a surface with respect to different grid strengths and R/X ratios. Case study results of proposed eigenvalue analysis are compared with those of admittance model-based stability analysis as well as time-domain simulation using a switching model in Matlab/Simulink, validating the effectiveness and accuracy of the proposed eigenvalue analysis for stability boundary identification.

grid-connected inverters↗
DOE OSTI2026

On Transforming the Specific Force Vector for an ECEF Inertial Navigation System

An inertial measurement unit accumulates the instantaneous specific force applied to the sensor over a small window of time. Importantly, this accumulation occurs in the sensor’s frame, yet the specific force must be mathematically integrated in a different navigation frame, and these two frames may be rotating with respect to each other. Consequently, a correction is needed to not only account for the frame difference but also for the time variation in the frame difference. Many references provide derivations of this correction [1], but the derivation is often abbreviated, and important assumptions are sometimes left unstated. This memo addresses this gap and provides a clear derivation of the specific force transform with all assumptions specifically stated. This derivation is specific to the Earth-centered, Earth-fixed navigation frame, but, hopefully, the derivation is sufficiently clear that readers can apply it to other navigation frames. Generally, notation follows the same format as described in both [1, 2].

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Coriolis forces modify magnetostatic ponderomotive potentials

It is possible to produce a ponderomotive effect in a plasma system without time-varying fields, if the plasma flows over spatial oscillations in the field. This can be achieved by superimposing a spatially oscillatory perturbation on a guide field, then setting up an electric field perpendicular to the guide field to drive flow over the perturbation. However, subtle distinctions in the structure of the resulting electric field can entirely change the behavior of the resulting ponderomotive force. Previous work has shown that, in slab models, these distinctions can be explained in terms of the polarization of the effective wave that appears in the co-moving frame. Here, we consider what happens to this picture in a cylindrical system, where the transformation to the co-moving (rotating) frame is not inertial. It turns out that the non-inertial nature of this frame transformation can lead to counterintuitive behavior, partly due to the appearance of parallel (magnetic-field-aligned) electric fields in the rotating frame even in cases where none existed in the laboratory frame. Apart from the academic interest of this study, the practical impact lies in being better able to anticipate the antenna configuration on the plasma periphery of a cylindrical plasma that will lead to optimal ponderomotive barrier formation in the interior plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Proton Improvement Plan II (PIP-II) Clock and Timing

The PIP-II timing system is planned to be a two-part system consisting of a global timing system (referred to as ACLK) that provides high level, event-based timing for the whole Fermilab accelerator complex while the second part is a RF synchronized clock system unique to the PIP-II Linac itself (referred to as LCLK). The ACLK System will make use of an external 10 MHz GPS based signal source as a reference for its 650 MHz phase lock as it is the reference frequency for the TCLK output of the system (needed to support legacy hardware around complex). The LCLK System will use a PIP-II Linac RF reference (162.5 MHz) from the Linac LLRF system to allow beam synchronized event placement. Both the ACLK and LCLK systems will have a clock output with a data frame of 16 event bits + 32 data bits with frames broadcast at 650 MHz, phase locked to the 10MHz reference in the case of ACLK and the PIP-II LLRF sourced 162.5 MHz reference in the case of LCLK.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

PRISMA: PARALLEL REFINEMENT AND INTEGRATION SYSTEM FOR MULTI-AZIMUTHAL ANALYSIS

The Parallel Refinement and Integration System for Multi-azimuthal Analysis (PRISMA, version 1.1.0) is a Python application for processing X-ray diffraction (XRD) image data. PRISMA wraps GSAS-II to perform azimuthally-binned peak refinement, computes per-frame strain and d-spacing from those fits, and provides three PyQt5 graphical interfaces: (1) a Recipe Builder for selecting GSAS-II control (.imctrl) files, optional mask (.immask) files or threshold-ased masking, reference and experiment image sets, peaks, zimuthal range and bin size, and an optional ceria-based auto-calibration; (2) a Batch Processor that uses Dask on local workstations and pure MPI (mpi4py.futures.MPICommExecutor) on HPC to distribute GSAS-II refinement across cores or compute nodes and write results to a 4-dimensional (peaks x frames x azimuths x measurements) Zarr dataset; and (3) a Data Analyzer that renders heatmaps of fit parameters, strain, frame-to-frame deltas, and percent-change-vs-reference, and exports user-defined subsections to CSV or Excel. The peak-refinement algorithm is deterministic. Benchmark on ALCF Crux: a 20,000-image set, single-peak fit in frame mode with 44 azimuthal bins on 128 nodes x 128 workers, 48 seconds total wall time.

Lorenzo Martin, Maria De La Cinta [Argonne Nationa↗
DOE OSTI2024

Breakthrough Conductivity Enhancement in Deep Eutectic Solvents via Grotthuss–Type Proton Transport

There is an increasing demand for the development of ion-conducting electrolytes for energy storage systems. Much attention is directed toward deep eutectic solvents as potential candidates. In the search for highly conductive systems, the possibility of designing deep eutectic solvents with Grotthuss-type proton transport is widely overlooked. Herein, ethaline, a mixture of choline chloride and ethylene glycol is used in a 1:2 molar ratio, to induce a significant conductivity increase with the addition of water and sulfuric acid (H 2 SO 4 ). The achieved breakthrough conductivity is analyzed experimentally and simulated with ab initio molecular dynamics (AIMD). At sufficient water content, an H-bonding network is formed that leads to a significant breakthrough conductivity based on H 2 SO 4 -derived proton transfer following the long-established Grotthuss proton transport mechanism. This result is substantiated by the positive deviation from the ideal KCl line in the Walden plot. Specifically, the data series positioned above the reference line indicates a Grotthuss mechanism in action. The AIMD simulations demonstrate proton transfer between water and ethylene glycol, supported by simulation frames captured at various times.

36 MATERIALS SCIENCE↗