SEARCH · Engineering Papers

Results for “method of multipliers”

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.

SourceResearch connectionFrom yearThrough year

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

At least 487 records · Page 27

All sources NASA NTRS DOE OSTI

NASA NTRS2023

A Finite Element Approach for Simplified 2D Nonlinear Dynamic Contact/Impact Analysis

In this paper, a simplified numerical approach for finite element dynamic analysis of an inelastic solid structure subjected to solid object impact is presented. The approach approximates the impacting solid as the selected multiple nodes, for which mass of the impactor is distributed. The node-to-segment contact formulation with the penalty constraint technique incorporated is employed to impose contact conditions between the nodes and the surface of the receiver structure. The node-to-segment algorithm is integrated into Newton–Raphson time integration scheme and the Lagrange multiplier technique is applied to enforce the identical displacements for the selected nodes throughout the analysis process. The approach is verified using two-dimensional plane strain models considering elastic-perfectly-plastic material behavior. The results obtained using the proposed approach are in a good agreement with those simulated using a commercial finite element code, ABAQUS dynamic/implicit, in terms of displacements and stress distribution fields. The proposed approach is shown to be computationally superior to general finite element method-based contact/impact analysis without significantly sacrificing the accuracy.

Seungwook Seok↗

NASA NTRS1993

Spur-reduced digital sinusoid synthesis

This article presents and analyzes a technique for reducing the spurious signal content in digital sinusoid synthesis. Spurious-harmonic (spur) reduction is accomplished through dithering both amplitude and phase values prior to word-length reduction. The analytical approach developed for analog quantization is used to produce new bounds on spur performance in these dithered systems. Amplitude dithering allows output word-length reduction without introducing additional spurs. Effects of periodic dither similar to those produced by a pseudonoise (PN) generator are analyzed. This phase-dithering method provides a spur reduction of 6(M plus one) dB per phase bit when the dither consists of M uniform variates. While the spur reduction is at the expense of an increase in system noise, the noise power can be made white, making the power spectral density small. This technique permits the use of a smaller number of phase bits addressing sinusoid lookup tables, resulting in an exponential decrease in system complexity. Amplitude dithering allows the use of less complicated multipliers and narrower data paths in purely digital applications, as well as the use of coarse resolution, highly linear digital to analog converters (DAC's) to obtain spur performance limited by the DAC linearity rather than its resolution.

Flanagan, M. J.↗

NASA NTRS2017

A Fuzzy Logic Approach for Separation Assurance and Collision Avoidance for Unmanned Aerial Systems

In the coming years, operations in low altitude airspace will vastly increase as the capabilities and applications of small Unmanned Aerial Systems (sUAS) continue to multiply. Therefore, solutions to managing sUAS in highly congested airspace must be explored. In this study, a Fuzzy Logic based approach was used to help mitigate the risk of collisions between aircraft using separation assurance and collision avoidance techniques. The system was evaluated for its effectiveness to mitigate the risk of mid-air collisions between aircraft. This system utilizes only current state information and can resolve potential conflicts without knowledge of intruder intent. The avoidance logic was verified using formal methods and shown to select the correct action in all instances. Additionally, the Fuzzy Logic Controllers were shown to always turn the vehicles in the correct direction. Numerical testing demonstrated that the avoidance system was able to prevent a mid-air collision between two sUAS in all tested cases. Simulations were also performed in a three-dimensional environment with a heterogenous fleet of sUAS performing a variety of realistic missions. Simulations showed that the system was 99.98 effective at preventing mid-air collisions when separation assurance was disabled (unmitigated case) and 100 effective when enabled (mitigated case).

Verification & Validation of a Fuzzy Logic Control↗

NASA NTRS1995

Method and apparatus for spur-reduced digital sinusoid synthesis

A technique for reducing the spurious signal content in digital sinusoid synthesis is presented. Spur reduction is accomplished through dithering both amplitude and phase values prior to word-length reduction. The analytical approach developed for analog quantization is used to produce new bounds on spur performance in these dithered systems. Amplitude dithering allows output word-length reduction without introducing additional spurs. Effects of periodic dither similar to that produced by a pseudo-noise (PN) generator are analyzed. This phase dithering method provides a spur reduction of 6(M + 1) dB per phase bit when the dither consists of M uniform variates. While the spur reduction is at the expense of an increase in system noise, the noise power can be made white, making the power spectral density small. This technique permits the use of a smaller number of phase bits addressing sinusoid look-up tables, resulting in an exponential decrease in system complexity. Amplitude dithering allows the use of less complicated multipliers and narrower data paths in purely digital applications, as well as the use of coarse-resolution, highly-linear digital-to-analog converters (DAC's) to obtain spur performance limited by the DAC linearity rather than its resolution.

Zimmerman, George A.↗

NASA NTRS2003

Mechanical Amplifier for a Piezoelectric Transducer

A mechanical amplifier has been devised to multiply the stroke of a piezoelectric transducer (PZT) intended for use at liquid helium temperatures. Interferometry holds the key to high angular resolution imaging and astrometry in space. Future space missions that will detect planets around other solar systems and perform detailed studies of the evolution of stars and galaxies will use new interferometers that observe at mid- and far-infrared wavelengths. Phase-measurement interferometry is key to many aspects of astronomical interferometry, and PZTs are ideal modulators for most methods of phase measurement, but primarily at visible wavelengths. At far infrared wavelengths of 150 to 300 m, background noise is a severe problem and all optics must be cooled to about 4 K. Under these conditions, piezos are ill-suited as modulators, because their throw is reduced by as much as a factor of 2, and even a wavelength or two of modulation is beyond their capability. The largest commercially available piezo stacks are about 5 in. (12.7 cm) long and have a throw of about 180 m at room temperature and only 90 m at 4 K. It would seem difficult or impossible to use PZTs for phase measurements in the far infrared were it not for the new mechanical amplifier that was designed and built.

Moore, James↗

DOE OSTI2024

Structure preserving and energy dissipative contact approaches for implicit dynamics

In this work, several structure preserving and energy dissipative contact approaches are proposed and evaluated. The time integration schemes considered are general with regard to the version of constraint type, but here the emphasis was on mortar contact. The proposed mortar contact approach conserves both linear and angular momentum for mortar contact in a novel way. The proposed time integration scheme can conserve energy or provide strict contact dissipation. In addition, the proposed scheme enforces both gap constraints and gap velocity constraints (i.e., persistency). Using a midstep time integrator often causes energy dissipation during initial impact, here this energy can be recovered in a novel way. Enforcing the gap velocity constraint mitigates the contact chatter of the contact pressure and nodes in many problems. Whereas some approaches enforce the persistency condition and gap constraints simultaneously during the solution of the equations of motion (EOM) requiring multipliers for both constraints included in the equation set, here the gap constraint is solved through the equations of motion and the persistency condition is satisfied in the time integration scheme by the velocity update after the equations of motion. It is shown that this approach is strictly dissipative in that a plastic contact condition can be achieved. Analogous to a coefficient of restitution for rigid bodies, any dissipated energy can then be returned upon release if energy conservation is desired. Structure preserving methods are good for long-time dynamics simulations and energy conserving and strictly dissipative methods can overcome stability issues associated with standard time integration algorithms such as the Newmark method.

42 ENGINEERING↗

DOE OSTI2026

Sap Velocity Data for East River Watershed Sites (2023-2025)

This dataset includes sap velocity measurements for aspen (populus tremuloides), fir (abies lasiocarpa), spruce (Engelmann spruce) and lodgepole pine (pinus contorta) trees at nine sites in the East River Watershed near Gothic, CO. This dataset was generated following a similar method as a previous dataset (Dataset. doi:10.15485/1647654) but was conducted at different sites in the area and now includes lodgepole pine. Site selection was done to explicitly improve understanding of topographic controls of tree water use. The data collection began in June 2023 and we provide data until December 2025 - though data collection is ongoing. The sap flux data were collected using ICT SFM1 sensors and are presented in both units of cm h^-1 and as kg h^-1 by multiplying the sap flux by the sapwood area of the tree. All sap flow data has been been corrected using estimates of wounding diameter, water content of wood and sap wood depth. We also provide a normalized sap velocity estimate by subtracting each measurement from that trees' annual minimum and dividing by that trees' annual maximum. This provides data for each tree and year on a 0-1 scale. This data entry contains one CSV file that includes all available sap flow data. The timestamps are provided in local time as year, day of year and hour and each measurement contains an associated latitude, longitude and site number which can be used to identify trees in a given stand. Each species is given a numeric value as: 1= populus tremuloides, 2=Engelmann spruce, 3=abies lasiocarpa and 4=pinus contorta.

abies↗

NASA NTRS1994

Buckling and vibration analysis of laminated panels using VICONOPT

The analysis aspects of the 23,000-line FORTRAN program VIPASA with constraints and optimization (VICONOPT) are described. Overall stiffness matrices assembled from the earlier exact vibration and instability of plate assemblies including shear and anisotropy (VIPASA) flat plate stiffnesses are optionally coupled by Langrangian multipliers to find critical buckling loads, or natural frequencies of undamped vibration, of prismatic assemblies of anisotropic flat plates with arbitrarily located point supports or simple transverse supporting frames. The longitudunal continuity of typical wing and fuselage panels is closely approximated because the solutions are for the infinitely long structure obtained by repeating a bay and its supports longitudinally. Any longitudinally invariant in-plane plate stresses are permitted, and very rapid solutions are guaranteed by numerous refinements, including multilevel substructuring and a method for repetitive cross sections that is exact for regular polygons used to represent cylinders. Modal displacements and stresses in or between plies of laminated plates are calculated and plotted, with values being recovered at all nodes of substructures. Comparison with usual approximate finite-element methods confirms that, for comparably converged solutions, VICONOPT is typically between 100 and 10(exp 4) times faster.

Kennedy, David↗

NASA NTRS2006

Estimating Random Errors Due to Shot Noise in Backscatter Lidar Observations

In this paper, we discuss the estimation of random errors due to shot noise in backscatter lidar observations that use either photomultiplier tube (PMT) or avalanche photodiode (APD) detectors. The statistical characteristics of photodetection are reviewed, and photon count distributions of solar background signals and laser backscatter signals are examined using airborne lidar observations at 532 nm using a photon-counting mode APD. Both distributions appear to be Poisson, indicating that the arrival at the photodetector of photons for these signals is a Poisson stochastic process. For Poisson-distributed signals, a proportional, one-to-one relationship is known to exist between the mean of a distribution and its variance. Although the multiplied photocurrent no longer follows a strict Poisson distribution in analog-mode APD and PMT detectors, the proportionality still exists between the mean and the variance of the multiplied photocurrent. We make use of this relationship by introducing the noise scale factor (NSF), which quantifies the constant of proportionality that exists between the root-mean-square of the random noise in a measurement and the square root of the mean signal. Using the NSF to estimate random errors in lidar measurements due to shot noise provides a significant advantage over the conventional error estimation techniques, in that with the NSF uncertainties can be reliably calculated from/for a single data sample. Methods for evaluating the NSF are presented. Algorithms to compute the NSF are developed for the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) lidar and tested using data from the Lidar In-space Technology Experiment (LITE). OCIS Codes:

Liu, Zhaoyan↗

DOE OSTI2025

Multiplexing and Demultiplexing Signals for Radiography Application Using the Discrete Fourier Transform

Our goal is to develop an X-ray phase-contrast imaging system that can provide excellent soft tissue contrast of phase, attenuation, and small-angle scatter. We propose to replace the common system of G0, G1, and G2 gradings with a biprism array to replace the G1 grading and introduce a novel X-ray tube designed to replace the motion of the phase stepping grading G2. The proposed X-ray tube uses temporal multiplexing to provide simultaneous virtual “electronic phase stepping.” In this work the discrete Fourier transform is used to separate from the composite measurement individual X-ray phase contrast measurements sampled at different frequencies. The method performs a discrete Fourier transform of a composite refence sequence to obtain using the frequency amplitudes calibration factors needed to extract the X-ray phase contrast measurement amplitudes from the composite image. The composite reference sequence is the sum of the individual sequences, at different frequencies, with amplitudes of one. The method takes the discrete Fourier transform of this composite reference sequence; whereby, the amplitude of each frequency component is compared with the total sum of its stand-alone sequence amplitude. A calibration factor is determined so that the amplitude of this composite reference frequency times the calibration factor must equal the total sum of the sequence amplitude—the zero-frequency amplitude of the discrete Fourier transform of its stand-alone sequence. To demultiplex the composite measured signal these calibration factors are multiplied by the amplitudes of the frequency components of the discrete Fourier transform of the composite X-phase-contrast measurement to obtain the amplitude of each frequency encoded measurement. Using these calibration factors, we demonstrate with the discrete Fourier transform in Mathematica the extraction of individual images from a composite image that one would expect obtaining from our proposed new X-ray phase contrast imaging system. We then demonstrate as an example how using images from X-ray phase contrast data one can calculate phase, attenuation and the dark field images using grading phase step data supplied to use from Microworks, GmbH in Karlsruhe, Germany.

42 ENGINEERING↗

DOE OSTI2021

Method and apparatus for shape-based energy analysis of solids

A computer-readable medium stores instructions including storing a model of a physical structure and defining a mesh for the model. The mesh includes finite elements, each defined by a respective set of edges. The instructions include, for each finite element, identifying a governing differential equation and a set of complementary functions that exactly satisfy the governing differential equation. The instructions include determining an applied physical stimulus for the physical structure. The instructions include generating an energy optimization model that minimizes a difference between internal energy and external energy of the finite elements in response to the applied physical stimulus. The model includes a matrix of respective scalar multipliers for the complementary functions for each finite element. The instructions include transforming the matrix and calculating a physical parameter of interest. The instructions include, in response to the physical parameter not satisfying a design parameter, updating the model of the physical structure.

Spears, Robert E.↗

DOE OSTI2025

YSO implantation detector for beta-delayed neutron spectroscopy

A segmented-scintillator-based implantation detector was developed to study the energy distribution of β-delayed neutrons emitted from exotic isotopes. The detector comprises a 34 × 34 YSO scintillator coupled to an 8 × 8 anode Position-Sensitive Photo-Multiplier Tube (PSPMT) via a tapered light guide. The detector was used at RIBF, RIKEN for time-of-flight-based neutron spectroscopy measurement in the 78Ni region. The detector provides the position and timing resolution necessary for ion-beta correlations and time of flight measurements. The detector provides a high ~80% beta-detection efficiency and a sub-nanosecond timing resolution. This contribution discusses the details of the design, operation, implementation, and analysis developed to obtain neutron time-of-flight spectrum and the analysis methods in the context of neutron-rich nuclei in the 78Ni region.

Implantation detector↗

DOE OSTI2023

DEEPEN 3D PFA Favorability Models and 2D Favorability Maps at Newberry Volcano

DEEPEN stands for DE-risking Exploration of geothermal Plays in magmatic ENvironments. Part of the DEEPEN project involved developing and testing a methodology for a 3D play fairway analysis (PFA) for multiple play types (conventional hydrothermal, superhot EGS, and supercritical). This was tested using new and existing geoscientific exploration datasets at Newberry Volcano. This GDR submission includes images, data, and models related to the 3D favorability and uncertainty models and the 2D favorability and uncertainty maps. The DEEPEN PFA Methodology is based on the method proposed by Poux et al. (2020), which uses the Leapfrog Geothermal software with the Edge extension to conduct PFA in 3D. This method uses all available data to build a 3D geodata model which can be broken down into smaller blocks and analyzed with advanced geostatistical methods. Each data set is imported into a 3D model in Leapfrog and divided into smaller blocks. Conditional queries can then be used to assign each block an index value which conditionally ranks each block's favorability, from 0-5 with 5 being most favorable, for each model (e.g., lithologic, seismic, magnetic, structural). The values between 0-5 assigned to each block are referred to as index values. The final step of the process is to combine all the index models to create a favorability index. This involves multiplying each index model by a given weight and then summing the resulting values. The DEEPEN PFA Methodology follows this approach, but split up by the specific geologic components of each play type. These components are defined as follows for each magmatic play type: 1. Conventional hydrothermal plays in magmatic environments: Heat, fluid, and permeability 2. Superhot EGS plays: Heat, thermal insulation, and producibility (the ability to create and sustain fractures suitable for and EGS reservoir) 3. Supercritical plays: Heat, supercritical fluid, pressure seal, and producibility (the proper permeability and pressure conditions to allow production of supercritical fluid) More information on these components and their development can be found in Kolker et al., 2022. For the purposes of subsurface imaging, it is easier to detect a permeable fluid-filled reservoir than it is to detect separate fluid and permeability components. Therefore, in this analysis, we combine fluid and permeability for conventional hydrothermal plays, and supercritical fluid and producibility for supercritical plays. More information on this process is described in the following sections. We also project the 3D favorability volumes onto 2D surfaces for simplified joint interpretation, and we incorporate an uncertainty component. Uncertainty was modeled using the best approach for the dataset in question, for the datasets where we had enough information to do so. Identifying which subsurface parameters are the least resolved can help qualify current PFA results and focus future efforts in data collection. Where possible, the resulting uncertainty models/indices were weighted using the same weights applied to the respective datasets, and summed, following the PFA methodology above, but for uncertainty. There are two different versions of the Leapfrog model and associated favorability models: - v1.0: The first release in June 2023 - v2.1: The second release, with improvements made to the earthquake catalog (included additional identified events, removed duplicate events), to the temperature model (fixed a deep BHT), and to the index models (updated the seismicity-heat source index models for supercritical and EGS, and the resistivity-insulation index models for all three play types). Also uses the jet color map rather than the magma color map for improved interpretability. - v2.1.1: Updated to include v2.0 uncertainty results (see below for uncertainty model versions) There are two different versions of the associated uncertainty models: - v1.0: The first release in June 2023 - v2.0: The se...

15 GEOTHERMAL ENERGY↗

DOE OSTI2024

DEEPEN: Final 3D PFA Favorability Models and 2D Favorability Maps at Newberry Volcano

Part of the DEEPEN (DE-risking Exploration of geothermal Plays in magmatic ENvironments) project involved developing and testing a methodology for a 3D play fairway analysis (PFA) for multiple play types (conventional hydrothermal, superhot EGS, and supercritical). This was tested using new and existing geoscientific exploration datasets at Newberry Volcano. This GDR submission includes images, data, and models related to the 3D favorability and uncertainty models and the 2D favorability and uncertainty maps. The DEEPEN PFA Methodology, detailed in the journal article below, is based on the method proposed by Poux & O'brien (2020), which uses the Leapfrog Geothermal software with the Edge extension to conduct PFA in 3D. This method uses all available data to build a 3D geodata model which can be broken down into smaller blocks and analyzed with advanced geostatistical methods. Each data set is imported into a 3D model in Leapfrog and divided into smaller blocks. Conditional queries can then be used to assign each block an index value which conditionally ranks each block's favorability, from 0-5 with 5 being most favorable, for each model (e.g., lithologic, seismic, magnetic, structural). The values between 0-5 assigned to each block are referred to as index values. The final step of the process is to combine all the index models to create a favorability index. This involves multiplying each index model by a given weight and then summing the resulting values. The DEEPEN PFA Methodology follows this approach, but split up by the specific geologic components of each play type. These components are defined as follows for each magmatic play type: 1. Conventional hydrothermal plays in magmatic environments: Heat, fluid, and permeability 2. Superhot EGS plays: Heat, thermal insulation, and producibility (the ability to create and sustain fractures suitable for and EGS reservoir) 3. Supercritical plays: Heat, supercritical fluid, pressure seal, and producibility (the proper permeability and pressure conditions to allow production of supercritical fluid) More information on these components and their development can be found in Kolker et al., (2022). For the purposes of subsurface imaging, it is easier to detect a permeable fluid-filled reservoir than it is to detect separate fluid and permeability components. Therefore, in this analysis, we combine fluid and permeability for conventional hydrothermal plays, and supercritical fluid and producibility for supercritical plays. We also project the 3D favorability volumes onto 2D surfaces for simplified joint interpretation, and we incorporate an uncertainty component. Uncertainty was modeled using the best approach for the dataset in question, for the datasets where we had enough information to do so. Identifying which subsurface parameters are the least resolved can help qualify current PFA results and focus future efforts in data collection. Where possible, the resulting uncertainty models/indices were weighted using the same weights applied to the respective datasets, and summed, following the PFA methodology above, but for uncertainty.

15 GEOTHERMAL ENERGY↗

NASA NTRS1992

A multiple scales approach to sound generation by vibrating bodies

The problem of determining the acoustic field in an inviscid, isentropic fluid generated by a solid body whose surface executes prescribed vibrations is formulated and solved as a multiple scales perturbation problem, using the Mach number M based on the maximum surface velocity as the perturbation parameter. Following the idea of multiple scales, new 'slow' spacial scales are introduced, which are defined as the usual physical spacial scale multiplied by powers of M. The governing nonlinear differential equations lead to a sequence of linear problems for the perturbation coefficient functions. However, it is shown that the higher order perturbation functions obtained in this manner will dominate the lower order solutions unless their dependence on the slow spacial scales is chosen in a certain manner. In particular, it is shown that the perturbation functions must satisfy an equation similar to Burgers' equation, with a slow spacial scale playing the role of the time-like variable. The method is illustrated by a simple one-dimenstional example, as well as by three different cases of a vibrating sphere. The results are compared with solutions obtained by purely numerical methods and some insights provided by the perturbation approach are discussed.

Geer, James F.↗

NASA NTRS1983

Modified jet noise source model for twin-jet shielding analysis

An analytical method to estimate the influence that a jet of heated flow has on the noise emission from a parallel jet is presented. The shielding jet is modelled as a cylinder of constant cross-section in which the flow speed and temperature are uniform throughout. The jet noise emission is modelled by a point source with directivity imposed. The directivity term consists of: a self-noise term, a shear-noise term, and a convection factor. The self- and shear-noise terms each contain a basic directivity factor multiplying a spectral shape function. The various components are evaluated based on comparison with isothermal jet radiation experimental data. The modified source term is incorporated into the jet shielding model and compared to heated twin jet shielding data. The estimated spectra agree well except further downstream of the nozzle where peak of the noise spectrum estimated by the model lies approximately one octave below the experimental peak. The noise reduction estimated by the model agrees favorably with experiment in the near downstream region. This discrepancy is explained in terms of the shielding mechanism which is dominant far downstream.

Gerhold, C. H.↗

NASA NTRS1985

New results on the PLL false lock phenomenon

Conditions, based on the multipliers of a linear differential equation with periodic coefficients, for the occurrence of false lock in a phase lock loop (PLL) are presented. Two new and powerful algorithms are presented for computing the frequency difference between the reference and voltage controlled oscillator when the loop is false locked. The first algorithm expresses omega f sup -1 as a power series in the closed loop gain. The second algorithm expresses omega sup -1 as part of a solution to a nonlinear algebraic system of equations which can be solved by numerical methods. The power of both algorithms is illustrated by application and comparison.

Stensby, J. L.↗

NASA NTRS1998

Analysis of Physical Properties of Dust Suspended in the Mars Atmosphere

Methods for iteratively determining the infrared optical constants for dust suspended in the Mars atmosphere are described. High quality spectra for wavenumbers from 200 to 2000 1/cm were obtained over a wide range of view angles by the Mariner 9 spacecraft, when it observed a global Martian dust storm in 1971-2. In this research, theoretical spectra of the emergent intensity from Martian dust clouds are generated using a 2-stream source-function radiative transfer code. The code computes the radiation field in a plane-parallel, vertically homogeneous, multiply scattering atmosphere. Calculated intensity spectra are compared with the actual spacecraft data to iteratively retrieve the optical properties and opacity of the dust, as well as the surface temperature of Mars at the time and location of each measurement. Many different particle size distributions a-re investigated to determine the best fit to the data. The particles are assumed spherical and the temperature profile was obtained from the CO2 band shape. Given a reasonable initial guess for the indices of refraction, the searches converge in a well-behaved fashion, producing a fit with error of less than 1.2 K (rms) to the observed brightness spectra. The particle size distribution corresponding to the best fit was a lognormal distribution with a mean particle radius, r(sub m) 0.66 pm, and variance, omega(sup 2) = 0.412 (r(sub eff) = 1.85 microns, v(sub eff) =.51), in close agreement with the size distribution found to be the best fit in the visible wavelengths in recent studies. The optical properties and the associated single scattering properties are shown to be a significant improvement over those used in existing models by demonstrating the effects of the new properties both on heating rates of the Mars atmosphere and in example spectral retrieval of surface characteristics from emission spectra.

Snook, Kelly↗