Search NASA⌕ Search

SEARCH · Search NASA

Results for “force vector”

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 19 records
All sourcesNASA NTRSDOE OSTI
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 OSTI2021

Galerkin formulations of isogeometric shell analysis: Alleviating locking with Greville quadratures and higher-order elements

Here we propose new quadrature schemes that asymptotically require only four in-plane points for Reissner–Mindlin shell elements and nine in-plane points for Kirchhoff–Love shell elements in B-spline and NURBS-based isogeometric shell analysis, independent of the polynomial degree p of the elements. The quadrature points are Greville abscissae associated with pth-order B-spline basis functions whose continuities depend on the specific Galerkin formulations, and the quadrature weights are calculated by solving a linear moment fitting problem in each parametric direction. The proposed shell element formulations are shown through numerical studies to be rank sufficient and to be free of spurious modes. The studies reveal comparable accuracy, in terms of both displacement and stress, compared with fully integrated spline-based shell elements, while at the same time reducing storage and computational cost associated with forming element stiffness and mass matrices and force vectors. The high accuracy with low computational cost makes the proposed quadratures along with higher-order spline bases, in particular polynomial orders, p = 5 and 6, good choices for alleviating membrane and shear locking in shells.

42 ENGINEERING↗
DOE OSTI2021

A finite element formulation for deformation twinning induced strain localization in polycrystal magnesium alloys

Deformation twinning induces shear strain localization in hexagonal close-packed crystals and is critical for the material’s ductility and failure. Cracks often occur at twin-twin or twin-grain boundary intersections and propagate along twin bands. However, most crystal plasticity models for deformation twinning are based on a “pseudo-slip” approach and do not capture the localized deformation associated with the formation of each discrete twin band. The few exceptions are discrete twin models that involve very complex numerical algorithms and are often compromised in accuracy due to the numerical convergence. These factors make the discrete twin models hard to adopt. This paper proposes a modification to the conventional finite element weak form, to fully incorporate a twin-induced heterogeneous deformation that does not depend on the “pseudo-slip” assumption. The model starts by splitting the deformation gradient into elastic-slip-twinning components. The twin-induced deformation gradient component is computed separately by solving a microstructural evolution problem and then implemented into finite element weak form by constructing a global “twin-force” vector. The constitutive update (e.g., in the user-defined material subroutine, or UMAT, for ABAQUS) therefore avoids dealing with the twinning and recovers to the form of a regular slip-based crystal plasticity model. The results presented here indicate that the twin-induced strain localization and the associated stress-reversal phenomena near the twin band were naturally captured in the model, which was validated against an in-situ synchrotron X-ray micro-diffraction experiment.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A benchmark dataset for Hydrogen Combustion

The generation of reference data for deep learning models is challenging for reactive systems, and more so for combustion reactions due to the extreme conditions that create radical species and alternative spin states during the combustion process. Here, we extend intrinsic reaction coordinate (IRC) calculations with ab initio MD simulations and normal mode displacement calculations to more extensively cover the potential energy surface for 19 reaction channels for hydrogen combustion. A total of ~290,000 potential energies and ~1,270,000 nuclear force vectors are evaluated with a high quality range-separated hybrid density functional, ωB97X-V, to construct the reference data set, including transition state ensembles, for the deep learning models to study hydrogen combustion reaction.

08 HYDROGEN↗
DOE OSTI2020

Effects of bubbles on the electrochemical behavior of hydrogen-evolving Si microwire arrays oriented against gravity

The size-distribution, coverage, electrochemical impedance, and mass-transport properties of H 2 gas-bubble films were measured for both planar and microwire-array platinized n + -Si cathodes performing the hydrogen-evolution reaction in 0.50 M H 2 SO 4 (aq). Inverted, planar n + -Si/Ti/Pt cathodes produced large, stationary bubbles which contributed to substantial increases in ohmic potential drops. In contrast, regardless of orientation, microwire array n + -Si/Ti/Pt cathodes exhibited a smaller layer of bubbles on the surface, and the formation of bubbles did not substantially increase the steady-state overpotential for H 2 (g) production. Experiments using an electroactive tracer species indicated that even when oriented against gravity, bubbles enhanced mass transport at the electrode surface. Microconvection due to growing and coalescing bubbles dominated effects due to macroconvection of gliding bubbles on Si microwire array cathodes. Electrodes that maintained a large number of small bubbles on the surface simultaneously exhibited low concentrations of dissolved hydrogen and small ohmic potential drops, thus exhibiting the lowest steady-state overpotentials. The results indicate that microstructured electrodes can operate acceptably for unassisted solar-driven water splitting in the absence of external convection and can function regardless of the orientation of the electrode with respect to the gravitational force vector.

25 ENERGY STORAGE↗
DOE OSTI2020

A generalized Boltzmann kinetic theory for strongly magnetized plasmas with application to friction

Coulomb collisions in plasmas are typically modeled using the Boltzmann collision operator, or its variants, which apply to weakly magnetized plasmas in which the typical gyroradius of particles significantly exceeds the Debye length. Conversely, O’Neil has developed a kinetic theory to treat plasmas that are so strongly magnetized that the typical gyroradius of particles is much smaller than the distance of the closest approach in a binary collision. Here, we develop a generalized collision operator that applies across the full range of magnetization strength. Since there is no closed-form solution for the scattering cross section when plasma is strongly magnetized, the input to the collision operator is obtained by numerically calculating particle trajectories. To demonstrate novel physics associated with strong magnetization, it is used to compute the friction force on a massive test charge. In addition to the traditional stopping power component, this is found to exhibit a transverse component that is perpendicular to both the velocity and Lorentz force vectors in the strongly magnetized regime, as was predicted recently using linear response theory. As a result, good agreement is found between the collision theory and linear response theory in the regime in which both apply, but the new collision theory also applies to stronger magnetization strength regimes than the linear response theory is expected to apply in.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

The magnetic states of a van der Waals ferromagnet CrGeTe 3 probed by vector-field magnetic force microscopy

Here, we present a study on the magnetic domain configuration in the van der Waals ferromagnet CrGeTe 3 using a vector-field cryogenic magnetic force microscopy. Our investigation demonstrates the influence of magnetic fields strength and angle on the magnetic domain configuration, showing the coexistence of striped and spike-like magnetic domains. This study contributes to understanding the impact of uniaxial magnetic anisotropy on the domain configuration of van der Waals ferromagnets, offering insights into stabilizing different domain patterns.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Transport Properties of Magnetized High-Energy Density Plasma (Final Project Report)

This report summarizes results of project DE-SC0016159 "Transport Properties of Magnetized High-Energy-Density Plasma,'' which ran from 7/15/2016 - 7/14/2021. The primary objectives of the work, as stated in the original proposal, were to develop a theory that describes transport coefficients in magnetized high energy density plasmas and to test the theory with molecular dynamics (MD) simulations. The science challenge is that strong ion coupling, strong magnetization of electrons, and partial degeneracy of electrons, cause the system to be in a regime that is not well described by current theory. The particular processes that were to be investigated include: ion stopping power, diffusion, electron-ion temperature relaxation, thermal conduction and viscosity. All of the primary research objectives were accomplished during the course of this work. In addition, some unexpected results were discovered along the way that led to new and productive research directions. This work resulted in 17 publications in well-respected peer-reviewed journals (primarily Physics of Plasmas and Physical Review E), and 2 more are being prepared for publication. A few of these were selected as Editor's Picks and one was published as a Rapid Communication. Highlights of the research results include: Transport phase space: The first identification of the parameter space that defines fundamental transport regimes in terms of the Coulomb coupling and magnetization parameters. Mean force kinetic theory: Systematic derivation of a kinetic theory for strongly coupled plasmas based on a new expansion parameter of the BBGKY hierarchy. This provided the derivation of a theory that we had previously posed phenomenologically, and also revealed a new term that captures the equation of state properties at all coupling strengths. Transverse friction force: Discovery that the friction force on a test charge in a strongly magnetized plasma includes a component that is perpendicular to the motion of the test charge in the plane defined by the velocity and magnetic field vectors. Gyrofriction force: Discovery that the friction force on a test charge in a plasma that is both strongly magnetized and strongly coupled includes a component of the force in the direction of the Lorentz force. Kinetic theory for warm dense matter: Extension of the mean force kinetic theory to include partial electron degeneracy, so that it applies to warm dense matter.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Giant Domain Wall Conductivity in Self‐Assembled BiFeO 3 Nanocrystals

Abstract Ever‐increasing demand on electronic devices with ultrahigh‐density non‐volatile data storage has attracted great interest in novel ferroelectric memories based on conductive ferroelectric domain walls. Embedded in an insulating material, ferroelectric domain walls have the capability of being (re)created, displaced, erased, and altered in their spatial configurations and electronic characteristics. However, the domain wall conductivities are in most cases not yet sufficiently high to ensure the current density required to drive read‐out circuits operating at high speeds. In this work, a giant domain wall current (>10 µA) of a single charged domain wall is obtained through conductive atomic force microscopy with a bias field of 4 V. This is achieved in self‐assembled BiFeO 3 nanocrystals grown by sol‐gel method on Nb‐doped SrTiO 3 substrates. Local configurations of domains and domain wall types are studied using vector piezoresponse force microscopy and high‐resolution transmission electronic microscopy. The enhancement of the wall current is shown to be due to the formation of conducting pathways of charged defects accumulated along domain walls and traversing the nanocrystals. The diverse domain walls can be manipulated by electric field in a perpendicular architecture. The perpendicular array structure of BiFeO 3 nanocrystals should have great potentials for developing perpendicular nanoelectronic prototypes.

Liu, Lisha↗
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↗
DOE OSTI2025

Efficient calculation of self magnetic field, self-force, and self-inductance for electromagnetic coils with rectangular cross-section

Abstract For designing high-field electromagnets, the Lorentz force on coils needs to be computed in order to design suitable support structures, and the inductance should be computed to evaluate the stored energy and dynamics. Also, the magnetic field and its variation inside the conductor is of interest for computing stress and strain, and due to superconducting quench limits. For these force, inductance, energy, and internal field calculations, the coils cannot be naively approximated as infinitesimally thin filaments due to divergences when the source and evaluation points coincide, so more computationally demanding calculations are usually required, resolving the finite cross-section of the conductors. Here, we present a new alternative method that enables the internal magnetic field vector, self-force, and self-inductance to be computed rapidly and accurately within a 1D filament model. The method is applicable to coils for which the curve center-line can have general noncircular shape, as long as the conductor width is small compared to the radius of curvature. This paper extends a previous calculation for circular-cross-section conductors (Hurwitz et al 2024 IEEE Trans. Magn. ) to consider the case of rectangular cross-section. The reduced model is derived by rigorous analysis of the singularity, regularizing the filament integrals such that they match the true high-dimensional integrals at high coil aspect ratio. The new filament model exactly recovers analytic results for a circular coil, and is shown to accurately reproduce full finite-cross-section calculations for a non-planar coil of a stellarator magnetic fusion device. Due to the efficiency of the model here, it is well suited for use inside design optimization.

Landreman, Matt (ORCID:000000027233577X)↗
DOE OSTI2019

Transverse force induced by a magnetized wake

The force on a test charge moving through a strongly magnetized plasma is calculated using linear response theory. Strong magnetization is found to generate a component of the force perpendicular to the velocity of the particle in the plane formed by the velocity and magnetic field vectors. Here, this transverse force is generated by an asymmetry with respect to the velocity vector in the induced electrostatic wake potential that is associated with the action of the Lorentz force on the background plasma. The direction depends on the speed of the test charge. If it is faster than a critical speed characteristic of the sound speed, it acts to reduce the component of velocity parallel to the magnetic field and to increase the gyroradius. In contrast, if the speed is below this critical speed, it acts to increase the velocity parallel to the magnetic field and to decrease the gyroradius. Because the transverse force is perpendicular to the velocity, it does not directly influence the total energy of the test charge. Nevertheless, it significantly alters the trajectory on a timescale associated with the Coulomb collision time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A generalized 3D elastic model for nanoscale, self-assembled oxide-metal thin films with pillar-in-matrix configurations

In recent years, functional oxide-metal based vertically aligned nanocomposite (VAN) thin films have gained interest due to their intriguing physical properties and multifunctionalities stemming from the complex interactions between the two phases in the film and the substrate. In this work, we develop a model for studying the energetics of these thin film systems, including the effects of both lattice mismatch and capillary forces due to interface curvature. Each phase is incorporated into the model using a phase indicator function, and we introduce the capillary forces as body forces using a vector density representation of the interface. The model is implemented using the finite element method to study the deformation of the thin film which is composed of Au nanopillars embedded in a La 0.7 Sr 0.3 MnO 3 (LSMO) matrix on an SrTiO 3 (STO) substrate. Further, the results suggest that the total energy is lowest for random configurations of pillars compared to ordered square and hexagonal lattice configurations, consistent with the random distribution of pillars found in experiments. Furthermore, we find that the interfacial energy dominates the total energy of each configuration, suggesting that interfacial energy in the system is an important design parameter for nanocomposite growth, along with the lattice mismatch.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Development of a machine-learning-based ionic-force correction model for quantum molecular dynamic simulations of warm dense matter

In this work Δ learning is used to map orbital-free density functional theory (OF-DFT) ionic forces to the corresponding Kohn-Sham (KS) DFT ionic forces. The development of the approximate force difference in terms of the ion positions is constructed and serves as a stand in for the ground truth force difference. Descriptor vectors for ion configurations are constructed using all distance between ions in conjunction with an indexing based on a nearest neighbor ranking. It is demonstrated that such a scheme of descriptors can uniquely describe an ionic configuration up to a rotation and reflection when no ambiguity in the nearest neighbor ranking exists. How to handle the case when an ambiguity exists in the nearest neighbor ranking is discussed. As a proof of principle, the model is trained and tested on warm dense hydrogen at temperatures between 1 and 15 eV. Once tested, the model was used to perform molecular dynamic simulations of warm dense hydrogen. Furthermore, the resulting energies and pressures are within 1% and 2% of their respective target KS values.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Neural Network Molecular Dynamics at Scale

Neural network molecular dynamics (NNMD) simulations could revolutionize atomistic modeling of materials with quantum-mechanical accuracy at a fraction of computational cost. However, popular NNMD frameworks are generally implemented for a single computing node, and conventional energy-based NN models still suffer from large time-to-solution (T2S), prohibiting the application of NNMD to challenging materials simulations encompassing large spatiotemporal scales. Consequently, no leadership-scale NNMD simulation has thus far been reported. Here, we present a scalable parallel NNMD software (RXMD-NN) based on our scalable reactive molecular dynamics simulation engine named RXMD. RXMD-NN has achieved high scalability up to 786,432 IBM BlueGene/Q cores involving 1.7 billion atoms. Furthermore, we have achieved 4.6-fold reduction of T2S by using a novel network that directly predicts atomic forces from feature vectors. Reduced T2S has for the first time allowed the study of large-scale off-stoichiometry effect in a widely used phase change material, Ge 2 Se 2 Te 5 , thereby resolving its "firstsharp diffraction peak mystery".

fist sharp diffraction peak↗
DOE OSTI2019

Nonlocal Transport and Implied Viscosity and Diffusivity throughout the Boundary Layer in LES of the Southern Ocean with Surface Waves

Abstract Observations from the Southern Ocean Flux Station provide a wide range of wind, buoyancy, and wave (Stokes) forcing for large-eddy simulation (LES) of deep Southern Ocean boundary layers. Almost everywhere there is a nonzero angle Ω between the shear and the stress vectors. Also, with unstable forcing there is usually a depth where there is stable stratification, but zero buoyancy flux and often a number of depths above where there is positive flux, but neutral stratification. These features allow nonlocal transports of buoyancy and of momentum to be diagnosed, using either the Eulerian or Lagrangian shear. The resulting profiles of nonlocal diffusivity and viscosity are quite similar when scaled according to Monin–Obukhov similarity theory in the surface layer, provided the Eulerian shear is used. Therefore, a composite shape function is constructed that may be generally applicable. In contrast, the deeper boundary layer appears to be too decoupled from the Stokes component of the Lagrangian shear. The nonlocal transports can be dominant. The diagnosed across-shear momentum flux is entirely nonlocal and is highly negatively correlated with the across-shear component of the wind stress, just as nonlocal and surface buoyancy fluxes are related. Furthermore, in the convective limit the scaling coefficients become essentially identical, with some consistency with atmospheric experience. The nonlocal contribution to the along-shear momentum flux is proportional to (1 − cosΩ) and is always countergradient, but is unrelated to the aligned wind stress component.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Accelerated Wind-Turbine Wake Recovery Through Actuation of the Tip-Vortex Instability

Advances in wind-plant control have often focused on more effectively balancing power between neighboring turbines. Wake steering is one such method that provides control-based improvements in a quasi-static way, but this does little to fundamentally change the wake recovery process, and thus, it has limited potential. This study investigates use of another control paradigm known as dynamic wake control (DWC) to excite the mutual inductance instability between adjacent tip-vortex structures, thereby accelerating the breakdown of the structures. The current work carries this approach beyond the hypothetical by applying the excitation via turbine control vectors that already exist on all modern wind turbines: blade pitch and rotor speed control. The investigation leverages a free-vortex wake method (FVWM) that allows a thorough exploration of relevant frequencies and amplitudes of harmonic forcing for each control vector (as well as the phase difference between the vectors for a tandem configuration) while still capturing the essential tip-vortex dynamics. The FVWM output feeds into a Fourier stability analysis working to pinpoint candidate DWC strategies suggesting fastest wake recovery. Near-wake length reductions of >80% are demonstrated, although without considering inflow turbulence. Analysis is provided to interpret these predictions considering the presence of turbulence in a real atmospheric inflow.

17 WIND ENERGY↗
DOE OSTI2023

Ground Vehicle Generalized Forces and Moment Governor Design Via Noncertainty-Equivalent Adaptive Prescribed Performance Control

Torque-vectoring technology demonstrates great potential in improving the safety and performance of ground vehicles. In this paper, a novel generalized forces and moment governor for torque vectoring is suggested. The proposed solution strategically combines prescribed performance control, noncertainty equivalent adaptive control design, and a smooth projection operator. The main advantage of the proposed control strategy lies in its capability to guarantee both the transient performance and prompt recovery of the desired deterministic behavior of the closed-loop adaptive system, even in the presence of parametric uncertainties. ASM simulation results are presented to validate the efficacy of the proposed generalized forces and moment governor and to demonstrate its superiority over a baseline solution.

Zhou, Xingyu↗