Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Arbitrary functions”

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 433 records · Page 24

Calculation of power spectrums from digital time series with missing data points

Two algorithms are developed for calculating power spectrums from the autocorrelation function when there are missing data points in the time series. Both methods use an average sampling interval to compute lagged products. One method, the correlation function power spectrum, takes the discrete Fourier transform of the lagged products directly to obtain the spectrum, while the other, the modified Blackman-Tukey power spectrum, takes the Fourier transform of the mean lagged products. Both techniques require fewer calculations than other procedures since only 50% to 80% of the maximum lags need be calculated. The algorithms are compared with the Fourier transform power spectrum and two least squares procedures (all for an arbitrary data spacing). Examples are given showing recovery of frequency components from simulated periodic data where portions of the time series are missing and random noise has been added to both the time points and to values of the function. In addition the methods are compared using real data. All procedures performed equally well in detecting periodicities in the data.

Murray, C. W., Jr.↗

Unsteady aerodynamic modeling for arbitrary motions

A study is presented on the unsteady aerodynamic loads due to arbitrary motions of a thin wing and their adaptation for the calculation of response and true stability of aeroelastic modes. In an Appendix, the use of Laplace transform techniques and the generalized Theodorsen function for two-dimensional incompressible flow is reviewed. New applications of the same approach are shown also to yield airloads valid for quite general small motions. Numerical results are given for the two-dimensional supersonic case. Previously proposed approximate methods, starting from simple harmonic unsteady theory, are evaluated by comparison with exact results obtained by the present approach. The Laplace inversion integral is employed to separate the loads into 'rational' and 'nonrational' parts, of which only the former are involved in aeroelastic stability of the wing. Among other suggestions for further work, it is explained how existing aerodynamic computer programs may be adapted in a fairly straightforward fashion to deal with arbitrary transients.

Edwards, J. W.↗

Theoretical studies of Resonance Enhance Stimulated Raman Scattering (RESRS) of frequency doubled Alexandrite laser wavelengths in cesium vapor

It is well known that the presence of a real atomic level which is nearly resonant with the pump field can greatly enhance the Raman emission cross section. In order to accurately calculate the Raman gain in systems where resonance enhancement plays a dominant role, expressions for the pump and signal susceptibilities must be derived. These expressions should be valid for arbitrary field strengths in order to allow for pump and signal saturation. In addition, the theory should allow for arbitrary longitudinal and transverse relaxation rates. This latter point is extremely vital for three level atomic systems such as the alkali earth metals since they do not have population reservoirs and can have widely varying spontaneous lifetimes on the three pertinent transitions. Moreover, the dephasing rates are strong functions of electron states and are therefore also different for the three coupled pairs of levels. These considerations are not as important when molecular systems are concerned since the large reservoir of rotational states serve to produce essentially equal longitudinal recovery rates for the population of the three levels. The three level system with three arbitrary longitudinal and transverse relaxation rates was solved. There is no need for setting either pair of rates equal and the expressions are valid for arbitrarily strong fields.

Lawandy, N. M.↗

An Analytical Model for the Prediction of a Micro-Dosimeter Response Function

A rapid analytical procedure for the prediction of a micro-dosimeter response function in low Earth orbit (LEO), correlated with the Space Transportation System (STS, shuttle) Tissue Equivalent Proportional Counter (TEPC) measurements is presented. The analytical model takes into consideration the energy loss straggling and chord length distribution of the detector, and is capable of predicting energy deposition fluctuations in a cylindrical micro-volume of arbitrary aspect ratio (height/diameter) by incoming ions through both direct and indirect (ray) events. At any designated (ray traced) target point within the vehicle, the model accepts the differential flux spectrum of Galactic Cosmic Rays (GCR) and/or trapped protons at LEO as input. On a desktop PC, the response function of TEPC for each ion in the GCR/trapped field is computed at the average rate of 30 seconds/ion. The ionizing radiation environment at LEO is represented by O'Neill fs GCR model (2004), covering charged particles in the 1 less than or equal to Z less than or equal to 28. O'Neill's free space GCR model is coupled with the Langley Research Center (LaRC) angular dependent geomagnetic cutoff model to compute the transmission coefficient in LEO. The trapped proton environment is represented by a LaRC developed time dependent procedure which couples the AP8MIN/AP8MAX, Deep River Neutron Monitor (DRNM) and F10.7 solar radio frequency measurements. The albedo neutron environment is represented by the extrapolation of the Atmospheric Ionizing Radiation (AIR) measurements. The charged particle transport calculations correlated with STS 51 and 114 flights are accomplished by using the most recent version (2005) of the LaRC deterministic High charge (Z) and Energy TRaNsport (HZETRN) code. We present the correlations between the TEPC model predictions (response function) and TEPC measured differential/integral spectra in the lineal energy (y) domain for both GCR and trapped protons, with the conclusion that the model correctly accounts for the increase in flux at low y values where energetic ions are the primary contributor. We further discuss that, even with the incorporation of angular dependency in the cutoffs, comparison of the GCR differential/integral flux between STS 51 and 114 TEPC measured data and current calculations indicates that there still exists an underestimation by the simulations at low to mid range y values. This underestimation is partly related the exclusion of the secondary pion particle production from the current version of HZETRN.

Badavi, Francis F.↗

A method for estimating the rolling moment due to spin rate for arbitrary planform wings

The application of aerodynamic theory for estimating the force and moments acting upon spinning airplanes is of interest. For example, strip theory has been used to generate estimates of the aerodynamic characteristics as a function of spin rate for wing-dominated configurations for angles of attack up to 90 degrees. This work, which had been limited to constant chord wings, is extended here to wings comprised of tapered segments. Comparison of the analytical predictions with rotary balance wind tunnel results shows that large discrepancies remain, particularly for those angles-of-attack greater than 40 degrees.

Poppen, W. A., Jr.↗

Synthesis of stiffened conical shells.

The development of a method to effect the automated minimum weight design of ring and stringer stiffened shells is presented. Membrane theory is used for the shell prebuckling analysis. The buckling analysis is based upon an arbitrary shell of revolution computer program. The structural analysis includes both buckling and yielding modes of failure. The synthesis involves the coupling of an exterior penalty function with a method for the unconstrained minimization of a function comprised of a sum of squares. Results of the application of the method to the design of the Viking Aeroshell cone are presented. The least weight Viking Aeroshell appears to be an all magnesium shell with ring stiffeners of hollow circular cross section. Because the method incorporates a general shell of revolution buckling analysis, it can be readily modified and applied to the design of any axisymmetrically loaded uniformly stiffened shell of revolution for which a membrane prebuckling solution exists.

Thornton, W. A.↗

Lattice model calculation of elastic and thermodynamic properties at high pressure and temperature

The elastic constants and the entire frequency spectrum were calculated up to high pressure for the alkali halides in the NaCl lattice, based on an assumed functional form of the inter-atomic potential. The quasiharmonic approximation is used to calculate the vibrational contribution to the pressure and the elastic constants at arbitrary temperature. By explicitly accounting for the effect of thermal and zero point motion, the adjustable parameters in the potential are determined to a high degree of accuracy from the elastic constants and their pressure derivatives measured at zero pressure. The calculated Gruneisen parameter, the elastic constants and their pressure derivatives are in good agreement with experimental results up to about 600 K. The model predicts that for some alkali halides the Grunesen parameter may decrease monotonically with pressure, while for others it may increase with pressure, after an initial decrease.

Demarest, H. H., Jr.↗

Asymptotic stability and instability of large-scale systems

The purpose of this paper is to develop new methods for constructing vector Lyapunov functions and broaden the application of Lyapunov's theory to stability analysis of large-scale dynamic systems. The application, so far limited by the assumption that the large-scale systems are composed of exponentially stable subsystems, is extended via the general concept of comparison functions to systems which can be decomposed into asymptotically stable subsystems. Asymptotic stability of the composite system is tested by a simple algebraic criterion. By redefining interconnection functions among the subsystems according to interconnection matrices, the same mathematical machinery can be used to determine connective asymptotic stability of large-scale systems under arbitrary structural perturbations.

Grujic, L. T.↗

Advanced satellite communication system

The objective of this research program was to develop an innovative advanced satellite receiver/demodulator utilizing surface acoustic wave (SAW) chirp transform processor and coherent BPSK demodulation. The algorithm of this SAW chirp Fourier transformer is of the Convolve - Multiply - Convolve (CMC) type, utilizing off-the-shelf reflective array compressor (RAC) chirp filters. This satellite receiver, if fully developed, was intended to be used as an on-board multichannel communications repeater. The Advanced Communications Receiver consists of four units: (1) CMC processor, (2) single sideband modulator, (3) demodulator, and (4) chirp waveform generator and individual channel processors. The input signal is composed of multiple user transmission frequencies operating independently from remotely located ground terminals. This signal is Fourier transformed by the CMC Processor into a unique time slot for each user frequency. The CMC processor is driven by a waveform generator through a single sideband (SSB) modulator. The output of the coherent demodulator is composed of positive and negative pulses, which are the envelopes of the chirp transform processor output. These pulses correspond to the data symbols. Following the demodulator, a logic circuit reconstructs the pulses into data, which are subsequently differentially decoded to form the transmitted data. The coherent demodulation and detection of BPSK signals derived from a CMC chirp transform processor were experimentally demonstrated and bit error rate (BER) testing was performed. To assess the feasibility of such advanced receiver, the results were compared with the theoretical analysis and plotted for an average BER as a function of signal-to-noise ratio. Another goal of this SBIR program was the development of a commercial product. The commercial product developed was an arbitrary waveform generator. The successful sales have begun with the delivery of the first arbitrary waveform generator.

Staples, Edward J.↗

Aeroelastic stability of wind turbine blade/aileron systems

Aeroelastic stability analyses have been performed for the MOD-5A blade/aileron system. Various configurations having different aileron torsional stiffness, mass unbalance, and control system damping have been investigated. The analysis was conducted using a code recently developed by the General Electric Company - AILSTAB. The code extracts eigenvalues for a three degree of freedom system, consisting of: (1) a blade flapwise mode; (2) a blade torsional mode; and (3) an aileron torsional mode. Mode shapes are supplied as input and the aileron can be specified over an arbitrary length of the blade span. Quasi-steady aerodynamic strip theory is used to compute aerodynamic derivatives of the wing-aileron combination as a function of spanwise position. Equations of motion are summarized herein. The program provides rotating blade stability boundaries for torsional divergence, classical flutter (bending/torsion) and wing/aileron flutter. It has been checked out against fixed-wing results published by Theodorsen and Garrick. The MOD-5A system is stable with respect to divergence and classical flutter for all practical rotor speeds. Aileron torsional stiffness must exceed a minimum critical value to prevent aileron flutter. The nominal control system stiffness greatly exceeds this minimum during normal operation. The basic system, however, is unstable for the case of a free (or floating) aileron. The instability can be removed either by the addition of torsional damping or mass-balancing the ailerons. The MOD-5A design was performed by the General Electric Company, Advanced Energy Program Department under Contract DEN3-153 with NASA Lewis Research Center and sponsored by the Department of Energy.

Strain, J. C.↗

Improved Model of a Mercury Ring Damper

A short document discusses the general problem of mathematical modeling of the three-dimensional rotational dynamics of rigid bodies and of the use of Euler parameters to eliminate the singularities occasioned by the use of Euler angles in such modeling. The document goes on to characterize a Hamiltonian model, developed by the authors, that utilizes the Euler parameters and, hence, is suitable for use in computational simulations that involve arbitrary rotational motion. In this formulation unlike in prior Euler-parameter-based formulations, there are no algebraic constraints. This formulation includes a general potential energy function, incorporates a minimum set of momentum variables, and takes an explicit state-space form convenient for numerical implementation. Practical application of this formulation has been demonstrated by the development of a new and simplified model of the rotational motion of a rigid rotor to which is attached a partially filled mercury ring damper. Models like this one are used in guidance and control of spin-stabilized spacecraft and gyroscope-stabilized seekers in guided missiles.

Fahrenthold, Eric P.↗

An Efficient, Highly Flexible Multi-Channel Digital Downconverter Architecture

In this innovation, a digital downconverter has been created that produces a large (16 or greater) number of output channels of smaller bandwidths. Additionally, this design has the flexibility to tune each channel independently to anywhere in the input bandwidth to cover a wide range of output bandwidths (from 32 MHz down to 1 kHz). Both the flexibility in channel frequency selection and the more than four orders of magnitude range in output bandwidths (decimation rates from 32 to 640,000) presented significant challenges to be solved. The solution involved breaking the digital downconversion process into a two-stage process. The first stage is a 2 oversampled filter bank that divides the whole input bandwidth as a real input signal into seven overlapping, contiguous channels represented with complex samples. Using the symmetry of the sine and cosine functions in a similar way to that of an FFT (fast Fourier transform), this downconversion is very efficient and gives seven channels fixed in frequency. An arbitrary number of smaller bandwidth channels can be formed from second-stage downconverters placed after the first stage of downconversion. Because of the overlapping of the first stage, there is no gap in coverage of the entire input bandwidth. The input to any of the second-stage downconverting channels has a multiplexer that chooses one of the seven wideband channels from the first stage. These second-stage downconverters take up fewer resources because they operate at lower bandwidths than doing the entire downconversion process from the input bandwidth for each independent channel. These second-stage downconverters are each independent with fine frequency control tuning, providing extreme flexibility in positioning the center frequency of a downconverted channel. Finally, these second-stage downconverters have flexible decimation factors over four orders of magnitude The algorithm was developed to run in an FPGA (field programmable gate array) at input data sampling rates of up to 1,280 MHz. The current implementation takes a 1,280-MHz real input, and first breaks it up into seven 160-MHz complex channels, each spaced 80 MHz apart. The eighth channel at baseband was not required for this implementation, and led to more optimization. Afterwards, 16 second stage narrow band channels with independently tunable center frequencies and bandwidth settings are implemented A future implementation in a larger Xilinx FPGA will hold up to 32 independent second-stage channels.

Goodhart, Charles E.↗

SSL Freeform Generator v1.00

The Orion spacecraft utilizes a system of Compact Unique Identifiers (CUIs) to identify telemetry measurements in the process of communicating with the ground. While avionic systems are still in development, simulation of telemetric data must be used for the development of dependent systems. The SuperScript Language (SSL) developed by Honeywell Aerospace provides users a programming language for simulation of this data. The use of tools that provide SSL generation capabilities at a higher level allows for faster and more efficient generation of simulation waveforms. Previous SSL Generator implementations such as Orion Scripted Interface Generator (OrionSIG) relied on software supported routines to generate simulation points for the SSL scripts. This provided relevant simulation capabilities, but was reduced in flexibility by the available software routines available. For this approach, software supported routines have been entirely dropped. Rather, the user is required to enter their own simulation waveform input. This allows arbitrarily complex functions to be utilized simulated with maximum degree of freedom for the user.: The idea for this version of the SSL Generator was a mapping of arbitrary waveforms to Orion CUI telemetry data. The SSL Generator implemented offers a higher degree of freedom to the user than previous versions while attempting to maintain abstraction from the underlying SSL script and the Orion CUI data types. As such, many mappings of waveforms to CUIs will simply not work. Certain implementation decisions were made to minimize the damage caused by these unintended side-effects, such as saturation of values exceeding the maximum or minimum bounds of the data type. Other details, such as signed values being represented in unsigned data types, present problems which cant be fixed to an acceptable degree of intended behavior without user input. Awareness of the limitations of all layers of the software stack will increase the productivity of the end-user.

telemetry↗

Mass Polarization Effect in He-like Systems

Eigenvalues for the ground state S and excited S and P states have been calculated for He-like systems, He, Li(+), Be(+2), and Ne(+8), using Hylleraas-type wave functions. These calculations have been carried out for a number of mass ratios R=mu/M=m(sub e)/(m(sub e)+M), where m(sub e) is the mass of the electron and M is the arbitrary mass of the nucleus. The eigenvalues are fitted to a 5th degree polynomial in R giving the mass polarization term (Delta (sub 1) x Delta (sub 2) and higher order corrections. The mass polarization term obtained from the fitting procedure agrees very well with the first-order result obtained directly. For example, in He we find E=E(sub 0)+Sigma(sup 5)(sub n=1)R(sup n)C(sub n), where E(sub 0)=-5.807448754 Ry and C(sub 1)=0.318138927 which agrees very well with the directly obtained first-order value 0.318138966083 and the result 0.318372 obtained by Yamanaka, using wave functions of the configuration-interaction form. We have carried out a similar calculation for the bound state of H(-).

Bhatia, A. K.↗

Fully unsteady subsonic and supersonic potential aerodynamics for complex aircraft configurations with applications to flutter

A general formulation is presented for the analysis of steady and unsteady, subsonic and supersonic aerodynamics for complex aircraft configurations. The theoretical formulation, the numerical procedure, the description of the program SOUSSA (steady, oscillatory and unsteady, subsonic and supersonic aerodynamics) and numerical results are included. In particular, generalized forces for fully unsteady (complex frequency) aerodynamics for a wing-body configuration, AGARD wing-tail interference in both subsonic and supersonic flows as well as flutter analysis results are included. The theoretical formulation is based upon an integral equation, which includes completely arbitrary motion. Steady and oscillatory aerodynamic flows are considered. Here small-amplitude, fully transient response in the time domain is considered. This yields the aerodynamic transfer function (Laplace transform of the fully unsteady operator) for frequency domain analysis. This is particularly convenient for the linear systems analysis of the whole aircraft.

Tseng, K.↗

Injection efficiency of bound modes

Previous work on efficiency of light injection into the core of a fiber from a thin film and a bulk distribution of sources in the cladding have used the fields of a weakly guiding fiber. This approximation simplifies the analysis of the power efficiency by introducing universal values for the eigenvalues of different fibers with the same V-number, but cannot predict accurately the behavior of the injected light into a fiber with arbitrary differences in indices of refraction. The exact field solution was used in the expressions of the power efficiency, p sub eff, and its behavior as a function of the fiber parameter was analyzed. Weakly guiding results obtained previously are confirmed. However, P sub eff does not always increase with the V-number but with the difference in the indices of refraction, eta sub core-eta sub clad. For the bulk distribution it was found that P sub eff increases with the wavelength, lambda, and decreases with the fiber core radius, a, i.e., it decreases with the V-number. However, for the thin film, the P sub eff remains almost constant with lambda and the fiber core radius.

Egalon, Claudio Oliveira↗

Advance finite element modeling of rotor blade aeroelasticity

An advanced beam finite element has been developed for modeling rotor blade dynamics and aeroelasticity. This element is part of the Element Library of the Second Generation Comprehensive Helicopter Analysis System (2GCHAS). The element allows modeling of arbitrary rotor systems, including bearingless rotors. It accounts for moderately large elastic deflections, anisotropic properties, large frame motion for maneuver simulation, and allows for variable order shape functions. The effects of gravity, mechanically applied and aerodynamic loads are included. All kinematic quantities required to compute airloads are provided. In this paper, the fundamental assumptions and derivation of the element matrices are presented. Numerical results are shown to verify the formulation and illustrate several features of the element.

Straub, F. K.↗