Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite precision”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Correlated Electromagnetic Levitation Actuator: A Reaction Sphere-Based Attitude Control System

To address problems experienced by current reaction wheels and control moment gyroscopebased attitude control systems (ACS), researchers at National Aeronautics and Space Administration’s (NASA’s) Marshall Space Flight Center (MSFC) have begun developing a reaction sphere actuator based on correlated electromagnetic levitation that will be immune to destructive bearing friction, momentum saturation, and gimbal lock. The Correlated Electromagnetic Levitation Actuator (CELA) advances the state of the art of reaction sphere ACSs by employing the concept of correlated magnetics. It is a frictionless, direct-drive reaction sphere that harnesses a unique technology with an array of applications across multiple disciplines. Correlated electromagnets function in a manner that is analogous to a matched filter; the convolution of two signals is peaked at the index representing the greatest match. For CELA, the signals are the patterns of magnetic flux density as a function of position. The magnitude of the convolution equates to an attractive or repulsive force, and these forces can be azimuthal or radial. The development of CELA is based in four distinct disciplines: Advanced Manufacturing, Prototype Development, Electromagnetic Modeling, and Controls. We are developing novel manufacturing techniques required to build arrays of permanent and electromagnet dipoles on curved surfaces. To print the permanent magnetic array, we have developed a probe with pyramidal magnets that will reside on a robotic arm to induce localized magnetic fields on a surface. We have also used high temperature ovens to erase dipole patterns from a permanent magnet by heating the surface to its Curie temperature. A number of test articles and prototypes have been developed using additive manufacturing methods. These prototypes have included hemispherical motors to test the drive algorithm, and a levitation test bed that demonstrates a magnetic bearing method based on attractive magnetic forces and ratiometric Hall effect sensors. We developed an array of electromagnetic dipoles on a printed circuit board (PCB) with individual H-bridges controlling each coil. This device created various flux density patterns and we measured their magnetic fields using a custom Hall effect 3-D probe and a LabVIEW virtual instrument. These data will serve as a benchmark for characterizing the accuracy of future models. Current work is focused on modeling the magnetic fields of our prototype arrays using COMSOL finite element analysis (FEA) and verifying the model against our test data. Accurate modeling will allow us to quickly test new patterns of electromagnets and their macro behavior. Eventually, the magnetic field models will be implemented in our controls simulations to facilitate precise control of the reaction sphere. Initial model results agree with field measurements to within 1 G (5% of measured flux density). Currently, we are testing different material properties of the electromagnets and their magnetic fields and thermal effects. These results will be used to refine the design of the electromagnetic dipoles. Our control efforts have centered on developing commutation, levitation, and field pattern shaping hardware in the form of breadboards and PCBs with software running on a local microcontroller. In addition, our partners developed MATLAB Simulink models to demonstrate a PID controller that mitigates disturbance forces resulting from the interaction of drive and levitation magnetics. Finally, we have designed a three-axis test stand that will be used in future work to demonstrate CELA’s orientation control capability.

reaction sphere↗
DOE OSTI2021

Mixed Precision Fermi-Operator Expansion on Tensor Cores from a Machine Learning Perspective

Here we present a second-order recursive Fermi-operator expansion scheme using mixed precision floating point operations to perform electronic structure calculations using tensor core units. A performance of over 100 teraFLOPs is achieved for half-precision floating point operations on Nvidia’s A100 tensor core units. The second-order recursive Fermi-operator scheme is formulated in terms of a generalized, differentiable deep neural network structure, which solves the quantum mechanical electronic structure problem. We demonstrate how this network can be accelerated by optimizing the weight and bias values to substantially reduce the number of layers required for convergence. We also show how this machine learning approach can be used to optimize the coefficients of the recursive Fermi-operator expansion to accurately represent the fractional occupation numbers of the electronic states at finite temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Observation and Modeling of the Thermal Runaway of High-Capacity Pouch Cells Due to an Internal Short Induced by an Indenter

Owing to the complexity of coupling mechanical and electrical solvers for finite element analysis (FEA) of mechanical abuse of a battery, the simulation should be performed by incorporating a local short circuit model activated based on the status of battery elements. However, the approach should also consider various other conditions of the battery simulated, such as the initial state of charge (SOC), for precise simulations of the internal short. This study proposes an approach developed with a mechanical model and single-particle battery model, along with a reaction kinetics model for thermal abuse using LS-DYNA. Several key parameters were selected, such as the activation criteria, kinetic abuse and element integration points during the internal short circuit induced by the indenter. The impact of these key parameters in the FEA on thermal runaway response for various initial SOCs from 25 to 100 % were analyzed. The FEA results obtained were compared against experimental results with various parameters such as activation criterion of internal short, kinetic abuse based on accelerating rate calorimetry (ARC) data. Finally, the root causes of the discrepancy between the FEA and experimental results are discussed regarding integration point what solver can make difference during calculation.

ENERGY STORAGE↗
DOE OSTI2025

AL-CRADA-2025-01 Final Report: American-Made Wind Turbine Materials Recycling Prize Accelerate! Phase 2

This report summarizes the Phase 2 achievements of the Wind Turbine Materials Recycling Prize project led by Critical Materials Recycling (CMR) in partnership with Ames National Laboratory. The team demonstrated the technical feasibility and economic viability of recovering and reusing rare-earth Nd-Fe-B magnets from end-of-life (EOL) wind turbine generators. Through mechanical disassembly, demagnetization, and precision processing, magnets were harvested and characterized to assess structural, compositional, and magnetic properties. Results showed that recovered magnets performed on par with commercial-grade counterparts, enabling their integration into actuator and motor prototypes without redesign. System-level modeling and finite-element simulations validated operational performance, while Life Cycle Cost (LCC) analysis revealed up to 20% material cost savings and 2–3% reductions in overall motor system costs. Life Cycle Assessment (LCA) further confirmed substantial environmental benefits, including a ~94% reduction in global warming potential relative to virgin magnet production. The project’s success was enabled by a tightly integrated partnership between industry and national lab experts, illustrating a viable path for domestic magnet-to-magnet recycling. These findings support a scalable, circular solution for rare-earth material recovery, reducing dependence on imported critical minerals and advancing sustainable clean energy technologies.

Lograsso, Thomas [Ames Laboratory (AMES), Ames, IA↗
NASA NTRS2009

On the Quality of the Nimbus 7 LIMS Version 6 Water Vapor Profiles and Distributions

This report describes the quality of the Nimbus 7 Limb Infrared Monitor of the Stratosphere (LIMS) water vapor (H2O) profiles of 1978/79 that were processed with a Version 6 (V6) algorithm and archived in 2002. The V6 profiles incorporate a better knowledge of the instrument attitude for the LIMS measurements along its orbits, leading to improvements for its temperature profiles and for the registration of its water vapor radiances with pressure. As a result, the LIMS V6 zonal-mean distributions of H2O exhibit better hemispheric symmetry than was the case from the original Version 5 (V5) dataset that was archived in 1982. Estimates of the precision and accuracy of the V6 H2O profiles are developed and provided. Individual profiles have a precision of order 5% and an estimated accuracy of about 19% at 3 hPa, 14% at 10 hPa, and 26% at 50 hPa. Profile segments within about 2 km of the tropopause are often affected by emissions from clouds that appear in the finite field-of-view of the detector for the LIMS H2O channel. Zonally-averaged distributions of the LIMS V6 H2O are compared with those from the more recent Microwave Limb Sounder (MLS) satellite experiment for November, February, and May of 2004/2005. The patterns and values of their respective distributions are similar in many respects. Effects of a strengthened Brewer-Dobson circulation are indicated in the MLS distributions of the recent decade versus those of LIMS from 1978/79. A tropical tape recorder signal is present in the 7-month time series of LIMS V6 H2O with lowest values in February 1979, and the estimated, annually-averaged "entry-level" H2O is 3.5 to 3.8 ppmv. It is judged that this historic LIMS water vapor dataset is of good quality for studies of the near global-scale chemistry and transport for pressure levels from 3 hPa to about 70 to 100 hPa.

Remsberg, E. E.↗
NASA NTRS2011

Porting a Hall MHD Code to a Graphic Processing Unit

We present our experience porting a Hall MHD code to a Graphics Processing Unit (GPU). The code is a 2nd order accurate MUSCL-Hancock scheme which makes use of an HLL Riemann solver to compute numerical fluxes and second-order finite differences to compute the Hall contribution to the electric field. The divergence of the magnetic field is controlled with Dedner?s hyperbolic divergence cleaning method. Preliminary benchmark tests indicate a speedup (relative to a single Nehalem core) of 58x for a double precision calculation. We discuss scaling issues which arise when distributing work across multiple GPUs in a CPU-GPU cluster.

Dorelli, John C.↗
NASA NTRS1980

Composite materials research and education program: The NASA-Virginia Tech composites program

Major areas of study include: (1) edge effects in finite width laminated composites subjected to mechanical, thermal and hygroscopic loading with temperature dependent material properties and the influence of edge effects on the initiation of failure; (2) shear and compression testing of composite materials at room and elevated temperatures; (3) optical techniques for precise measurement of coefficients of thermal expansion of composites; (4) models for the nonlinear behavior of composites including material nonlinearity and damage accumulation and verification of the models under biaxial loading; (5) compressive failure of graphite/epoxy plates with circular holes and the buckling of composite cylinders under combined compression and torsion; (6) nonlinear mechanical properties of borsic/aluminum, graphite/polyimide and boron/aluminum; (7) the strength characteristics of spliced sandwich panels; and (8) curved graphite/epoxy panels subjected to internal pressure.

Herakovich, C. T.↗
DOE OSTI2021

Equations of Motion for the Vertical Rigid-Body Rotor: Linear and Nonlinear Cases

Centuries ago, the prolific mathematician Leonhard Euler (1707–1783) wrote down the equations of motion (EOM) for the heavy symmetrical top with one point fixed. The resulting set of equations turned out to be nonlinear and had a limited number of closed-form solutions.Today, tools such as transfer matrix and finite elements enable the calculation of the rotordynamic properties for rotor-bearing systems. Some of these tools rely on the “linearized” version of the EOM to calculate the eigenvalues, unbalance response, or transients in these systems.In fact, industry standards mandate that rotors be precisely balanced to have safe operational characteristics. However, in some cases, the nonlinear aspect of the EOM should be considered.The purpose of this chapter is to show examples of how the linear vs. nonlinear formulations differ. This chapter also shows how excessive unbalance is capable of dramatically altering the behavior of the system and can produce chaotic motions associated with the “jump” phenomenon.

Laos, Hector↗
DOE OSTI2020

A least-squares finite element method based on the Helmholtz decomposition for hyperbolic balance laws

In this paper, a least-squares finite element method for scalar nonlinear hyperbolic balance laws is proposed and studied. The approach is based on a formulation that utilizes an appropriate Helmholtz decomposition of the flux vector and is related to the standard notion of a weak solution. This relationship, together with a corresponding connection to negative-norm least-squares, is described in detail. As a consequence, an important numerical conservation theorem is obtained, similar to the famous Lax–Wendroff theorem. The numerical conservation properties of the method in this paper do not fall precisely in the framework introduced by Lax and Wendroff, but they are similar in spirit as they guarantee that when L 2 convergence holds, the resulting approximations approach a weak solution to the hyperbolic problem. The least-squares functional is continuous and coercive in an H -1 -type norm, but not L 2 -coercive. Nevertheless, the L 2 convergence properties of the method are discussed. Convergence can be obtained either by an explicit regularization of the functional, that provides control of the L 2 norm, or by properly choosing the finite element spaces, providing implicit control of the L 2 norm. Numerical results for the inviscid Burgers equation with discontinuous source terms are shown, demonstrating the L 2 convergence of the obtained approximations to the physically admissible solution. The numerical method utilizes a least-squares functional, minimized on finite element spaces, and a Gauss–Newton technique with nested iteration. Finally, we believe that the linear systems encountered with this formulation are amenable to multigrid techniques and combining the method with adaptive mesh refinement would make this approach an efficient tool for solving balance laws (this is the focus of a future study).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Contact in the Unitary Fermi Gas across the Superfluid Phase Transition

A quantity known as the contact is a fundamental thermodynamic property of quantum many-body systems with short-range interactions. Determination of the temperature dependence of the contact for the unitary Fermi gas of infinite scattering length has been a major challenge, with different calculations yielding qualitatively different results. Here we use finite-temperature auxiliary-field quantum Monte Carlo (AFMC) methods on the lattice within the canonical ensemble to calculate the temperature dependence of the contact for the homogeneous spin-balanced unitary Fermi gas. We extrapolate to the continuum limit for 40, 66, and 114 particles, eliminating systematic errors due to finite-range effects. We observe a dramatic decrease in the contact as the superfluid critical temperature is approached from below, followed by a gradual weak decrease as the temperature increases in the normal phase. Our theoretical results are in excellent agreement with the most recent precision ultracold atomic gas experiments. Here, we also present results for the energy as a function of temperature in the continuum limit.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Improving tide-estuary representation in MPAS-Ocean (Final Report)

The overarching goal of this research is to produce a skill-assessed tide-estuary representation, including BGC (biogeochemistry) components of the ecosystem, into the MPAS-O global climate modeling system. Specifically, in the original proposal we proposed one-way nesting, into global MPAS-O, of a high-resolution model for TES (tidal and estuary systems) including the coupling of the BGC components. We have completed the proposed tasks during the past 4 years. Specifically, after extensive consultation with the MPAS-O group during the first year, we have decided to go beyond the original proposal and to two-way nest SCHISM into MPAS-O. As a result, SCHISM-CoSiNE (with CoSiNE being the BGC model) is now fully integrated into MPAS-O as a subroutine (more precisely, an analysis member) and can be invoked as MPAS-O is running with two-way exchange of information between the TES and global models (with proper spatial interpolation done on the fly). As originally proposed, we have also developed the new finite-volume (FV) model MPAS-OI (aka ELCIRC-sub) for TES. We have conducted several tests using the coupled MPAS-SCHISM- CoSiNE model and also MPAS-OI for Columbia River, San Francisco Bay and Chesapeake Bay. Four journal publications have been produced so far as the result of the project and 1 is under review (cf. Publication section). With these achievements, we will be ready to conduct fully coupled simulations of TES and global ocean anywhere, as soon as some missing capabilities in MPAS-O (e.g. tides, improved ability to handle the boundary conditions, and streamlining of the global BGC model) are worked out. The significance of the outcome of this project cannot be overstated. It marks the first time, as far as we know, that a global climate model is tightly coupled to a bona fide nearshore and terrestrial model (SCHISM system); in fact SCHISM is the first bona fide compound flooding model that seamlessly covers the coastal ocean, estuaries, rivers and watersheds (Ye et al. 2020; Zhang et al. 2020). The coupled model will reduce the errors and uncertainties of the MPAS-O global climate model system. It also complements and extends on-going E3SM project, as the latter seeks in the long-term to include components required to simulate impacts of sea-level change and storm surge on coastal regions, including wave models and focusing resolution in coastal and storm-track regions, as an important component of the Cryosphere-ocean system. Note that a spectral wave model (the unstructured-grid version of Wave Watch III, aka Wind Wave Model) is already incorporated inside the SCHISM system and is therefore directly accessible in the coupled system.

58 GEOSCIENCES↗
DOE OSTI2023

Delamination imaging in composites using cross-correlation method by non-contact air-coupled Lamb waves

Abstract This paper presents a cross-correlation damage detection technique using damage scattered waves of all directions cross-correlate with the incident waves for delamination imaging in carbon fiber reinforced polymer (CFRP) composite plate. The presented imaging can not only detect the delamination but also precisely quantify the damage with locations and sizes. First, a CFRP composite sample with interlaminar delamination by inserting the Teflon layer is designed and manufactured in house using the hot-press machine. And a three-dimensional model of composite sample is simulated using the finite element method. Next, the cross-correlation imaging algorithm is introduced in detail, and the cross-correlation imaging proof-of-concept study is conducted with the simulated Lamb waves in the composite sample. Finally, a fully non-contact air-coupled transducer and scanning laser Doppler vibrometer system with a single-mode Lamb wave method is established to actuate and sense the interrogating Lamb waves in the structure. The imaging method is experimentally implemented for the one delamination and two delaminations imaging and quantification.

Instruments & Instrumentation↗
DOE OSTI2022

Dynamic Response of an Electron Gas: Towards the Exact Exchange-Correlation Kernel

Precise calculations of dynamics in the homogeneous electron gas (jellium model) are of fundamental importance for design and characterization of new materials. In this work, we introduce a diagrammatic Monte Carlo technique based on algorithmic Matsubara integration that allows us to compute frequency and momentum resolved finite temperature response directly in the real frequency domain using series of connected Feynman diagrams. The data for charge response at moderate electron density are used to extract the frequency dependence of the exchange-correlation kernel at finite momenta and temperature. These results are as important for development of the time-dependent density functional theory for materials dynamics as ground state energies are for the density functional theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Equations of Motion for the Vertical Rigid-Body Rotor: Linear and Nonlinear Cases

Centuries ago, the prolific mathematician Leonhard Euler (1707–1783) wrote down the equations of motion (EOM) for the heavy symmetrical top with one point fixed. The resulting set of equations turned out to be nonlinear and had a limited number of closed-form solutions. Today, tools such as transfer matrix and finite elements enable the calculation of the rotor dynamic properties for rotor-bearing systems. Some of these tools rely on the “linearized” version of the EOM to calculate the eigenvalues, unbalance response, or transients in these systems. In fact, industry standards mandate that rotors be precisely balanced to have safe operational characteristics. However, in some cases, the nonlinear aspect of the EOM should be considered. The purpose of this paper is to show examples of how the linear vs. nonlinear formulations differ. This paper will also show how excessive unbalance is capable of dramatically altering the behavior of the system and can produce chaotic motions associated with the “jump” phenomenon.

42 ENGINEERING↗
DOE OSTI2026

Calculation of nuclear magnetic form factors [Slides]

The knowledge of how electron and neutrinos interact in nuclei is critical to an understanding of finite nucle and of astrophysical environments like neutron stars and supernovae. Within this project we are studying nuclear magnetic form factors that reflect how electrons interact with nuclei at low momenta. We implemented realistic treatments of many-nucleon correlations and currents to enable high-precision studies of electron physics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1989

Qualitative simulation for process modeling and control

A qualitative model is developed for a first-order system with a proportional-integral controller without precise knowledge of the process or controller parameters. Simulation of the qualitative model yields all of the solutions to the system equations. In developing the qualitative model, a necessary condition for the occurrence of oscillatory behavior is identified. Initializations that cannot exhibit oscillatory behavior produce a finite set of behaviors. When the phase-space behavior of the oscillatory behavior is properly constrained, these initializations produce an infinite but comprehensible set of asymptotically stable behaviors. While the predictions include all possible behaviors of the real system, a class of spurious behaviors has been identified. When limited numerical information is included in the model, the number of predictions is significantly reduced.

Dalle Molle, D. T.↗
NASA NTRS2002

Designing Crane Controls with Applied Mechanical and Electrical Safety Features

The use of overhead traveling bridge cranes in many varied applications is common practice. In particular, the use of cranes in the nuclear, military, commercial, aerospace, and other industries can involve safety critical situations. Considerations for Human Injury or Casualty, Loss of Assets, Endangering the Environment, or Economic Reduction must be addressed. Traditionally, in order to achieve additional safety in these applications, mechanical systems have been augmented with a variety of devices. These devices assure that a mechanical component failure shall reduce the risk of a catastrophic loss of the correct and/or safe load carrying capability. ASME NOG-1-1998, (Rules for Construction of Overhead and Gantry Cranes, Top Running Bridge, and Multiple Girder), provides design standards for cranes in safety critical areas. Over and above the minimum safety requirements of todays design standards, users struggle with obtaining a higher degree of reliability through more precise functional specifications while attempting to provide "smart" safety systems. Electrical control systems also may be equipped with protective devices similar to the mechanical design features. Demands for improvement of the cranes "control system" is often recognized, but difficult to quantify for this traditionally "mechanically" oriented market. Finite details for each operation must be examined and understood. As an example, load drift (or small motions) at close tolerances can be unacceptable (and considered critical). To meet these high functional demands encoders and other devices are independently added to control systems to provide motion and velocity feedback to the control drive. This paper will examine the implementation of Programmable Electronic Systems (PES). PES is a term this paper will use to describe any control system utilizing any programmable electronic device such as Programmable Logic Controllers (PLC), or an Adjustable Frequency Drive (AID) 'smart' programmable motion controller. Therefore the use of the term Programmable Electronic Systems (PES) is an encompassing description for a large spectrum of programmable electronic control devices.

Lytle, Bradford P.↗
NASA NTRS2018

Telescope Development for LISA

A laser interferometer space-based gravitational wave observatory requires an optical telescope to efficiently transfer laser light between pairs of widely-separated sciencecraft. The application is precision interferometric metrology, and therefore requires the telescope to have high optical pathlength stability, and low scattered light performance. We summarize recent activity on the development of a telescope to meet requirements for the LISA mission, including optimization of the optical design, development of a CAD and finite element model, and thinking about how to live with a high level of particulate contamination on the optics and still meet challenging scattered light requirements.

Livas, Jeff↗