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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS2011

On Calculation Methods and Results for Straight Cylindrical Roller Bearing Deflection, Stiffness, and Stress

The purpose of this study was to assess some calculation methods for quantifying the relationships of bearing geometry, material properties, load, deflection, stiffness, and stress. The scope of the work was limited to two-dimensional modeling of straight cylindrical roller bearings. Preparations for studies of dynamic response of bearings with damaged surfaces motivated this work. Studies were selected to exercise and build confidence in the numerical tools. Three calculation methods were used in this work. Two of the methods were numerical solutions of the Hertz contact approach. The third method used was a combined finite element surface integral method. Example calculations were done for a single roller loaded between an inner and outer raceway for code verification. Next, a bearing with 13 rollers and all-steel construction was used as an example to do additional code verification, including an assessment of the leading order of accuracy of the finite element and surface integral method. Results from that study show that the method is at least first-order accurate. Those results also show that the contact grid refinement has a more significant influence on precision as compared to the finite element grid refinement. To explore the influence of material properties, the 13-roller bearing was modeled as made from Nitinol 60, a material with very different properties from steel and showing some potential for bearing applications. The codes were exercised to compare contact areas and stress levels for steel and Nitinol 60 bearings operating at equivalent power density. As a step toward modeling the dynamic response of bearings having surface damage, static analyses were completed to simulate a bearing with a spall or similar damage.

Krantz, Timothy L.↗
NASA NTRS2011

High Precision Thermal, Structural and Optical Analysis of an External Occulter Using a Common Model and the General Purpose Multi-Physics Analysis Tool Cielo

The efficient simulation of multidisciplinary thermo-opto-mechanical effects in precision deployable systems has for years been limited by numerical toolsets that do not necessarily share the same finite element basis, level of mesh discretization, data formats, or compute platforms. Cielo, a general purpose integrated modeling tool funded by the Jet Propulsion Laboratory and the Exoplanet Exploration Program, addresses shortcomings in the current state of the art via features that enable the use of a single, common model for thermal, structural and optical aberration analysis, producing results of greater accuracy, without the need for results interpolation or mapping. This paper will highlight some of these advances, and will demonstrate them within the context of detailed external occulter analyses, focusing on in-plane deformations of the petal edges for both steady-state and transient conditions, with subsequent optical performance metrics including intensity distributions at the pupil and image plane.

integrated modeling↗
NASA NTRS2000

Dynamic Response of a Planetary Gear System Using a Finite Element/Contact Mechanics Model

The dynamic response of a helicopter planetary gear system is examined over a wide range of operating speeds and torques. The analysis tool is a unique, semianalytical finite element formulation that admits precise representation of the tooth geometry and contact forces that are crucial in gear dynamics. Importantly, no a priori specification of static transmission error excitation or mesh frequency variation is required; the dynamic contact forces are evaluated internally at each time step. The calculated response shows classical resonances when a harmonic of mesh frequency coincides with a natural frequency. However, peculiar behavior occurs where resonances expected to be excited at a given speed are absent. This absence of particular modes is explained by analytical relationships that depend on the planetary configuration and mesh frequency harmonic. The torque sensitivity of the dynamic response is examined and compared to static analyses. Rotation mode response is shown to be more sensitive to input torque than translational mode response.

Parker, Robert G.↗
NASA NTRS1983

On the off-axis tension test for unidirectional composites

The off axis tensile test was examined experimentally to obtain actual displacement fields over the surface of graphite polyimide coupon specimens. The experimental results were compared with approximate analytical solutions and generated finite element results. An optical method of high sensitivity moireinterferometry was used to determine the actual displacements to high precision. The approximate analytical solution and the finite element results compare very favorably with the measured centerline displacements in the test section, and the finite element displacement fields provide excellent agreement with the moiredisplacements throughout the specimen. A 15 degree fiber orientation and coupon aspect ratios of 5 and 15 are presented. Previously announced in STAR as N82-23553

Herakovich, C. T.↗
NASA NTRS1982

On the off-axis tensile test for unidirectional composites

The off axis tensile test was examined experimentally to obtain actual displacement fields over the surface of graphite polyimide coupon specimens. The experimental results were compared with approximate analytical solutions and generated finite element results. An optical method of high sensitivity moire interferometry was used to determine the actual displacements to high precision. The approximate analytical solution and the finite element results compare very favorably with the measured centerline displacements in the test section, and the finite element displacement fields provide excellent agreement with the moire displacements throughout the specimen. A 15 degree fiber orientation and coupon aspect ratios of 5 and 15 are presented.

Nemeth, M. P.↗
DOE OSTI2023

Precisely computing phonons via irreducible derivatives

Computing phonons from first principles is typically considered a solved problem, yet inadequacies in existing techniques continue to yield deficient results in systems with sensitive phonons. Here, in this study, we circumvent this issue using the lone irreducible derivative (LID) and bundled irreducible derivative (BID) approaches to computing phonons via finite displacements, where the former optimizes precision via energy derivatives and the latter provides the most efficient algorithm using force derivatives. A condition number optimized basis for BID is derived which guarantees the minimum amplification of error. Additionally, a hybrid LID-BID approach is formulated, in which select irreducible derivatives computed using LID replace BID results. We illustrate our approach on two prototypical systems with sensitive phonons: the shape memory alloy AuZn and metallic lithium. Comparing our resulting phonons in the aforementioned crystals to calculations in the literature reveals nontrivial inaccuracies. Our approaches can be fully automated, making them well suited for both niche systems of interest and high-throughput approaches.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Explicit block encodings of boundary value problems for many-body elliptic operators

Simulation of physical systems is one of the most promising use cases of future digital quantum computers. In this work we systematically analyze the quantum circuit complexities of block encoding the discretized elliptic operators that arise extensively in numerical simulations for partial differential equations, including high-dimensional instances for many-body simulations. When restricted to rectangular domains with separable boundary conditions, we provide explicit circuits to block encode the many-body Laplacian with separable periodic, Dirichlet, Neumann, and Robin boundary conditions, using standard discretization techniques from low-order finite difference methods. To obtain high-precision, we introduce a scheme based on periodic extensions to solve Dirichlet and Neumann boundary value problems using a high-order finite difference method, with only a constant increase in total circuit depth and subnormalization factor. We then present a scheme to implement block encodings of differential operators acting on more arbitrary domains, inspired by Cartesian immersed boundary methods. We then block encode the many-body convective operator, which describes interacting particles experiencing a force generated by a pair-wise potential given as an inverse power law of the interparticle distance. This work provides concrete recipes that are readily translated into quantum circuits, with depth logarithmic in the total Hilbert space dimension, that block encode operators arising broadly in applications involving the quantum simulation of quantum and classical many-body mechanics.

Kharazi, Tyler [University of California, Berkeley↗
DOE OSTI2024

Dynamics of the O ( 4 ) critical point in QCD: Critical pions and diffusion in model G

We present a detailed study of the finite momentum dynamics of the O ( 4 ) critical point of QCD, which lies in the dynamic universality class of “model G.” The critical scaling of the model is analyzed in multiple dynamical channels. For instance, the finite momentum analysis allows us to precisely extract the pion dispersion curve below the critical point. The pion velocity is in striking agreement with the predictions relation and static universality. The pion damping rate and velocity are both consistent with the dynamical critical exponent ζ = 3 / 2 of model G. Similarly, although the critical amplitude for the diffusion coefficient of the conserved O ( 4 ) charges is small, it is clearly visible both in the restored phase and with finite explicit symmetry breaking, and its dynamical scaling is again consistent with ζ = 3 / 2 . We determine a new set of universal dynamical critical amplitude ratios relating the diffusion coefficient to a suitably defined order parameter relaxation time. We also show that in a finite volume simulation, the chiral condensate diffuses on the coset manifold in a manner consistent with dynamical scaling, and with a diffusion coefficient that is determined by the transport coefficients of hydrodynamic pions. Finally, the amplitude ratios (together with other nonuniversal amplitudes also reported here) compile all relevant information for further studies of model G both in and out of equilibrium. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1987

Application of optimization techniques to spacecaft fuel usage minimization in deep space navigation

Mathematical analysis of the minimization of spacecraft fuel usage for both impulsive and finite motor burns is presented. A high precision integrated trajectory search program (SEPV) and several optimization software libraries are used to solve minimum fuel usage problems. The SEPV program has the capacity to vary either the initial spacecraft state or the finite burn parameters to acquire a specified set of target values. Several test examples for the Voyager 2 Uranus Encounter and the Galileo Jupiter Orbiter are presented to show that spacecraft fuel consumption can be minimized in targeting maneuver strategies. The fuel savings achieved by the optimum solution can be significant.

Wang, Tseng-Chan↗
DOE OSTI2024

Adaptive finite differencing in high accuracy electronic structure calculations

Abstract A multi-order Adaptive Finite Differencing (AFD) method is developed for the kinetic energy operator in real-space, grid-based electronic structure codes. It uses atomic pseudo orbitals produced by the corresponding pseudopotential codes to optimize the standard finite difference (SFD) operators for improved precision. Results are presented for a variety of test systems and Bravais lattice types, including the well-known Δ test for 71 elements in the periodic table, the Mott insulator NiO, and borax decahydrate, which contains covalent, ionic, and hydrogen bonds. The tests show that an 8th-order AFD operator leads to the same average Δ value as that achieved by plane-wave codes and is typically far more accurate and has a much lower computational cost than a 12th-order SFD operator. The scalability of real-space electronic calculations is demonstrated for a 2016-atom NiO cell, for which the computational time decreases nearly linearly when scaled from 18 to 144 CPU-GPU nodes.

Briggs, E. L. (ORCID:0000000343983492)↗
NASA NTRS1990

Nonlinear viscoelastic response of highly filled elastomers under multiaxial finite deformation

A biaxial tester was used to obtain precise biaxial stress responses of highly filled, high strain capability elastomers. Stress-relaxation experiments show that the time-dependent part of the relaxation response can be reasonably approximated by a function which is strain and biaxiality independent. Thus, isochronal data from the stress-relaxation curves can be used to determine the stored energy density function. The complex behavior of the elastomers under biaxial deformation may be caused by dewetting.

Peng, Steven T. J.↗
NASA NTRS1979

Recent finite element studies in plasticity and fracture mechanics

The paper reviews recent work on fundamentals of elastic-plastic finite-element analysis and its applications to the mechanics of crack opening and growth in ductile solids. The presentation begins with a precise formulation of incremental equilibrium equations and their finite-element forms in a manner valid for deformations of arbitrary magnitude. Special features of computational procedures are outlined for accuracy in view of the near-incompressibility of elastic-plastic response. Applications to crack mechanics include the analysis of large plastic deformations at a progressively opening crack tip, the determination of J integral values and of limitations to J characterizations of the intensity of the crack tip field, and the determination of crack tip fields in stable crack growth.

Rice, J. R.↗
NASA NTRS1993

A comparison of viscoelastic damping models

Modern finite element methods (FEM's) enable the precise modeling of mass and stiffness properties in what were in the past overwhelmingly large and complex structures. These models allow the accurate determination of natural frequencies and mode shapes. However, adequate methods for modeling highly damped and high frequency dependent structures did not exist until recently. The most commonly used method, Modal Strain Energy, does not correctly predict complex mode shapes since it is based on the assumption that the mode shapes of a structure are real. Recently, many techniques have been developed which allow the modeling of frequency dependent damping properties of materials in a finite element compatible form. Two of these methods, the Golla-Hughes-McTavish method and the Lesieutre-Mingori method, model the frequency dependent effects by adding coordinates to the existing system thus maintaining the linearity of the model. The third model, proposed by Bagley and Torvik, is based on the Fractional Calculus method and requires fewer empirical parameters to model the frequency dependence at the expense of linearity of the governing equations. This work examines the Modal Strain Energy, Golla-Hughes-McTavish and Bagley and Torvik models and compares them to determine the plausibility of using them for modeling viscoelastic damping in large structures.

Slater, Joseph C.↗
NASA NTRS1979

Finite-strain large-deflection elastic-viscoplastic finite-element transient response analysis of structures

A method of analysis for thin structures that incorporates finite strain, elastic-plastic, strain hardening, time dependent material behavior implemented with respect to a fixed configuration and is consistently valid for finite strains and finite rotations is developed. The theory is formulated systematically in a body fixed system of convected coordinates with materially embedded vectors that deform in common with continuum. Tensors are considered as linear vector functions and use is made of the dyadic representation. The kinematics of a deformable continuum is treated in detail, carefully defining precisely all quantities necessary for the analysis. The finite strain theory developed gives much better predictions and agreement with experiment than does the traditional small strain theory, and at practically no additional cost. This represents a very significant advance in the capability for the reliable prediction of nonlinear transient structural responses, including the reliable prediction of strains large enough to produce ductile metal rupture.

Rodal, J. J. A.↗
DOE OSTI2021

Interactions of two and three mesons including higher partial waves from lattice QCD

We study two- and three-meson systems composed either of pions or kaons at maximal isospin using Monte Carlo simulations of lattice QCD. Utilizing the stochastic LapH method, we are able to determine hundreds of two- and three-particle energy levels, in nine different momentum frames, with high precision. We fit these levels using the relativistic finite-volume formalism based on a generic effective field theory in order to determine the parameters of the two- and three-particle K-matrices. We find that the statistical precision of our spectra is sufficient to probe not only the dominant s-wave interactions, but also those in d waves. In particular, we determine for the first time a term in the three-particle K-matrix that contains two-particle d waves. We use three N f = 2 + 1 CLS ensembles with pion masses of 200, 280, and 340 MeV. This allows us to study the chiral dependence of the scattering observables, and compare to the expectations of chiral perturbation theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1985

A General-applications Direct Global Matrix Algorithm for Rapid Seismo-acoustic Wavefield Computations

A new matrix method for rapid wave propagation modeling in generalized stratified media, which has recently been applied to numerical simulations in diverse areas of underwater acoustics, solid earth seismology, and nondestructive ultrasonic scattering is explained and illustrated. A portion of recent efforts jointly undertaken at NATOSACLANT and NORDA Numerical Modeling groups in developing, implementing, and testing a new fast general-applications wave propagation algorithm, SAFARI, formulated at SACLANT is summarized. The present general-applications SAFARI program uses a Direct Global Matrix Approach to multilayer Green's function calculation. A rapid and unconditionally stable solution is readily obtained via simple Gaussian ellimination on the resulting sparsely banded block system, precisely analogous to that arising in the Finite Element Method. The resulting gains in accuracy and computational speed allow consideration of much larger multilayered air/ocean/Earth/engineering material media models, for many more source-receiver configurations than previously possible. The validity and versatility of the SAFARI-DGM method is demonstrated by reviewing three practical examples of engineering interest, drawn from ocean acoustics, engineering seismology and ultrasonic scattering.

Schmidt, H.↗
NASA NTRS1988

Piezoelectric polymer-based isolation mount for articulated pointing systems on large flexible spacecraft

When a pointing system is attached to a large flexible body that is subject to continuous excitation, a mechanical isolator becomes a key element in achieving high performance. An active softmount based on the piezoelectric polymer poly (vinylidene fluoride) is considered here for use in a precision pointing system for Space Station. A finite element model of the softmount is developed and added to simple Space Station and payload models for performance analysis. Both linear frequency domain and nonlinear time domain simulations are carried out in order to demonstrate the wideband disturbance rejection capabilities of the design. Nonlinear geometric effects of large nodal deflections are also considered.

Sirlin, Samuel W.↗