SEARCH · Engineering Papers
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.
By use of rational approximation, a procedure was found which allows calculation of the z-transfer function from a restricted number of time sequence values. The order of the system need not be known. Transformation to the s-domain is easily provided.
Description of an analytically obtained optimum adaptive delta modulator-demodulator configuration. The device utilizes two past samples to obtain a step size which minimizes the mean square error for a Markov-Gaussian source. The optimum system is compared, using computer simulations, with a linear delta modulator and an enhanced Abate delta modulator. In addition, the performance is compared to the rate distortion bound for a Markov source. It is shown that the optimum delta modulator is neither quantization nor slope-overload limited. The highly nonlinear equations obtained for the optimum transmitter and receiver are approximated by piecewise-linear equations in order to obtain system equations which can be transformed into hardware. The derivation of the experimental system is presented.
An algorithm and ILLIAC computer program, developed for the simulation of homogeneous incompressible turbulence in the presence of an applied mean strain, are described. The turbulence field is represented spatially by a truncated triple Fourier series (spectral method) and followed in time using a fourth-order Runge-Kutta algorithm. These include: (1) transformation of variables suggested by Taylor's sudden-distortion theory; (2) implicit viscous diffusion by use of an integrating factor; (3) implicit pressure calculation suggested by Taylor's sudden-distortion theory, and (4) inexpensive control of aliasing by random and phased coordinate shifts.
The derivation of a Hamiltonian field theory for nonlinear density waves in Saturn's rings is discussed. Starting with a Hamiltonian for a discrete system of gravitating streamlines, an averaged Hamiltonian is obtained by successive applications of Lie transforms. The transformation may be carried out to any desired order in q, where q is the nonlinearity parameter defined in the work of Shu, et al (1985) and Borderies et al (1985). Subsequent application of the Wentzel-Kramer-Brillouin Method approximation yields an asymptotic field Hamiltonian. Both the nonlinear dispersion relation and the wave action transport equation are easily derived from the corresponding Lagrangian by the standard variational principle.
The phenomenon of forced unsteady separation and eruption of boundary-layer vorticity is a highly-complex, high-Reynolds number flow phenomenon, which abruptly leads to the formation of a dynamic stall vortex as demonstrated earlier by the authors for a NACA 0015 airfoil undergoing constant rate pitch-up motion. This, as well as the results of other researchers, have convincingly demonstrated a complex vortical structure within the state of unsteady separation prior to the evolution of dynamic stall. This phenomenon of vortex eruption, although observed in studying dynamic stall phenomena, is also associated with transition from laminar to turbulence flow and its generic nature has been stressed by many researchers including the present investigators. An unsteady Navier-Stokes (NS) analysis is developed for arbitrarily maneuvering bodies using velocity-vorticity variables; this formulation is nearly form-invariant under a generalized non-inertial coordinate transformation. A fully-implicit uniformly second-order accurate method is used, with the nonlinear convective terms approximated using a biased third-order upwind differencing scheme to be able to simulate higher-Re flows. No explicit artificial dissipation is added. The numerical method is fully vectorized and currently achieves a computational index of 7 micro-seconds per time step per mesh point, using a single processor on a CRAY Y-MP. The simulation results show that the energetic free shear from the leading edge is responsible for the wall viscous layer to abruptly erupt near the center of the counterclockwise rotating eddy in the unsteady boundary layer. Primary, secondary, tertiary and quaternary vortices have been observed before the dynamic stall vortex evolves and gathers its maximum strength. This study will discuss the simulation results of Reynolds number up to Re = 45,000 and will also discuss the efforts of initial acceleration in a specific maneuver, on the evolution of the stall vortex.
Conduct of space exploration is undergoing a significant transformation. Initial reconnaissance missions are giving way to long duration observations with data-intensive instruments, in situ investigations and complex operations. To keep pace, a transformation in the Deep Space Network is in order.
Fast Fourier transform (FFT) techniques are applied to the problem of finding the flow due to source or vortex distributions in the field outside an airfoil or other two-dimensional body. Either the complex potential or the complex velocity may be obtained to a high order of accuracy, with computational effort similar to that required by second-order fast Poisson solvers. These techniques are applicable to general flow problems with compressibility and rotation. An example is given of their use for inviscid compressible flow.
An algorithm is presented for efficient p-iterative solution of the Lambert/Gauss orbit-determination problem using second-order Newton iteration. The algorithm is based on a universal transformation of Kepler's time-of-flight equation and approximate inverse solutions of this equation for short-way and long-way flight paths. The approximate solutions provide both good starting values for iteration and simplified computation of the second-order term in the iteration formula. Numerical results are presented which indicate that in many cases of practical significance (except those having collinear position vectors) the algorithm produces at least eight significant digits of accuracy with just two or three steps of iteration.
Several gamma-ray bursts in the BATSE data have sufficiently long durations and complex temporal structures with pulses that appear to be spaced quasi-periodically. In order to test and quantify these periods we have applied fast Fourier transformations (FFT) to all these events. We have also performed cross spectral analyses of the FFT of the two extreme (high-low) energy bands in each case to determine the lead/lag of the pulses in different energies.
Frozen orbits of the Hill problem are determined in the double averaged problem, where short and long period terms are removed by means of Lie transforms. The computation of initial conditions of corresponding quasi periodic solutions in the non-averaged problem is straightforward for the perturbation method used provides the explicit equations of the transformation that connects the averaged and non-averaged models. A fourth order analytical theory reveals necessary for the accurate computation of quasi periodic, frozen orbits.
This paper deals with the dynamic response of anisotropic laminated composite flat panels exposed to sonic boom and explosive blast-type loadings. The pertinent governing equations incorporating transverse shear deformation, transverse normal stress, the higher order effects as well as the viscous structural damping are solved by using the integral-transform technique. The obtained results are compared with their counterparts obtained within the framework of the first order transverse shear deformation and the classical plate theories and some conclusions concerning their range of applicability are outlined. The paper also contains a detailed analysis of the influence played by the various parameters characterizing the considered pressure pulses as well as the material and geometry of the plate.
Registration of test and reference images is a key component of nearly all PSP data reduction techniques. This is done to ensure that a test image pixel viewing a particular point on the model is ratioed by the reference image pixel which views the same point. Typically registration is needed to account for model motion due to differing airloads when the wind-off and wind-on images are taken. Registration is also necessary when two cameras are used for simultaneous acquisition of data from a dual-frequency paint. This presentation will discuss the advantages and disadvantages of several different image registration techniques. In order to do so, it is necessary to propose both an accuracy requirement for image registration and a means for measuring the accuracy of a particular technique. High contrast regions in the unregistered images are most sensitive to registration errors, and it is proposed that these regions be used to establish the error limits for registration. Once this is done, the actual registration error can be determined by locating corresponding points on the test and reference images, and determining how well a particular registration technique matches them. An example of this procedure is shown for three transforms used to register images of a semispan model. Thirty control points were located on the model. A subset of the points were used to determine the coefficients of each registration transform, and the error with which each transform aligned the remaining points was determined. The results indicate the general superiority of a third-order polynomial over other candidate transforms, as well as showing how registration accuracy varies with number of control points. Finally, it is proposed that image registration may eventually be done away with completely. As more accurate image resection techniques and more detailed model surface grids become available, it will be possible to map raw image data onto the model surface accurately. Intensity ratio data can then be obtained by a "model surface ratio," rather than an image ratio. The problems and advantages of this technique will be discussed.
Many physical processes appear to exhibit fractional order behavior that may vary with time or space. The continuum of order in the fractional calculus allows the order of the fractional operator to be considered as a variable. This paper develops the concept of variable and distributed order fractional operators. Definitions based on the Riemann-Liouville definitions are introduced and behavior of the operators is studied. Several time domain definitions that assign different arguments to the order q in the Riemann-Liouville definition are introduced. For each of these definitions various characteristics are determined. These include: time invariance of the operator, operator initialization, physical realization, linearity, operational transforms. and memory characteristics of the defining kernels. A measure (m2) for memory retentiveness of the order history is introduced. A generalized linear argument for the order q allows the concept of "tailored" variable order fractional operators whose a, memory may be chosen for a particular application. Memory retentiveness (m2) and order dynamic behavior are investigated and applications are shown. The concept of distributed order operators where the order of the time based operator depends on an additional independent (spatial) variable is also forwarded. Several definitions and their Laplace transforms are developed, analysis methods with these operators are demonstrated, and examples shown. Finally operators of multivariable and distributed order are defined in their various applications are outlined.
The thermodynamics of organic chemistry under mild aqueous conditions was examined in order to begin to understand its influence on the structure and operation of metabolism and its antecedents. Free energies were estimated for four types reactions of biochemical importance carbon-carbon bond cleavage and synthesis, hydrogen transfer between carbon groups, dehydration of alcohol groups, and aldo-keto isomerization. The energies were calculated for mainly aliphatic groups composed of carbon, hydrogen, and oxygen. The energy values showed that (1) when carbon-carbon bond cleavage involves two different types of functional groups, transfer of the shared electron-pair to the more reduced carbon group is energetically favored over transfer to the more oxidized carbon group, and (2) the energy of carbon-carbon bond transformation is strongly dependent on the type of functional group that donates the shared electron-pair during cleavage, and the group that accepts the shared electron-pair during synthesis, and (3) the energetics of C-C bond transformation is determined primarily by the half-reaction energies of the couples: carbonyl/carboxylic acid, carboxylic acid/carbon dioxide, alcohol/carbonyl, and hydrocarbon/alcohol. The energy of hydrogen-transfer between carbon groups was found to depend on the functional group class of both the hydrogen-donor and hydrogen-acceptor. From these and other observations we concluded that the chemistry of the origin of metabolism (and to a lesser degree modem metabolism) is strongly constrained by the (1) limited disproportionation energy of organic substrates that can be dissipated in a few irreversible reactions, (2) the energy-dominance of few half-reaction couples in carbon-carbon bond transformation that establishes whether a chemical reaction is energetically irreversible, reversible or unfeasible, and (3) the dependence of the transformation-energy on the oxidation state of carbon groups (functional group type) which is contingent on prior reactions in the synthetic pathway.
Fractional-Order systems, or systems containing fractional derivatives and integrals, have been studied by many in the engineering area. Additionally, very readable discussions, devoted specifically to the subject, are presented by Oldham and Spanier, Miller and Ross, and Pudlubny (1999a). It should be noted that there are a growing number of physical systems whose behavior can be compactly described using fractional system theory. Of specific interest to electrical engineers are long lines, electrochemical processes, dielectric polarization, colored noise, viscoelastic materials, and chaos. With the growing number of applications, it is important to establish a theory of control for these fractional-order systems, and for the potential use of fractional-order systems as feedback compensators. This topic is addressed in this paper. The first section discusses the control of fractional-order systems using a vector space representation, where initialization is included in the discussion. It should be noted that Bagley and Calico and Padovan and Sawicki both present a fractional state-space representation, which do not include the important historic effects. Incorporation of these effects based on the initialized fractional calculus is presented . The control methods presented in this paper are based on the initialized fractional order system theory. The second section presents an input-output approach. Some of the problems encountered in these sections are: a) the need to introduce a new complex plane to study the dynamics of fractional-order systems, b) the need to properly define the Laplace transform of the fractional derivative, and c) the proper inclusion of the initialization response in the system and control formulation. Following this, the next section generalizes the proportional-plus-integral-control (PI-control) and PID-control (PI-plus- derivative) concepts using fractional integrals. This is then further generalized using general fractional- order compensators. Finally the compensator concept is generalized by the use of a continuum of fractions in the compensator via the concept of order-distributions. The last section introduces fractional feedback in discrete-time.
An efficient, direct, second-order solver for the discrete solution of a class of two-dimensional separable elliptic equations on the sphere (which generally arise in implicit and semi-implicit atmospheric models) is presented. The method involves a Fourier transformation in longitude and a direct solution of the resulting coupled second-order finite-difference equations in latitude. The solver is made efficient by vectorizing over longitudinal wave-number and by using a vectorized fast Fourier transform routine. It is evaluated using a prescribed solution method and compared with a multigrid solver and the standard direct solver from FISHPAK.
We present techniques for discovering and exploiting regularity in large curvilinear data sets. The data can be based on a single mesh or a mesh composed of multiple submeshes (also known as zones). Multi-zone data are typical to Computational Fluid Dynamics (CFD) simulations. Regularities include axis-aligned rectilinear and cylindrical meshes as well as cases where one zone is equivalent to a rigid-body transformation of another. Our algorithms can also discover rigid-body motion of meshes in time-series data. Next, we describe a data model where we can utilize the results from the discovery process in order to accelerate large data visualizations. Where possible, we replace general curvilinear zones with rectilinear or cylindrical zones. In rigid-body motion cases we replace a time-series of meshes with a transformed mesh object where a reference mesh is dynamically transformed based on a given time value in order to satisfy geometry requests, on demand. The data model enables us to make these substitutions and dynamic transformations transparently with respect to the visualization algorithms. We present results with large data sets where we combine our mesh replacement and transformation techniques with out-of-core paging in order to achieve significant speed-ups in analysis.
Proposed digital image processor improved version of Programmable Remapper, which performs geometric and radiometric transformations on digital images. Features include overlapping and variably sized preimages. Overcomes some of limitations of image-warping circuit boards implementing only those geometric tranformations expressible in terms of polynomials of limited order. Also overcomes limitations of existing Programmable Remapper and made to perform transformations at video rate.