Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “theoretical and computational techniques”

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 109 records · Page 6

Liquefaction dynamics of n-octadecane in cylindrical coordinates

The study presents one method of mathematically modeling the effects of gravity-induced free convection upon the melting phenomena of a finite paraffin slab. The theoretical model is based on the numerical computer solution of the solid phase and liquid phase energy equations, coupled with an imposed ideal-viscous flow approximation. An implicit alternating-direction technique is used as the finite-difference approximation for the basic energy equations in the numerical solution. The theoretical predictions are compared to experimental data from an annular test cell, in which n-octadecane, C18H38, was used as the test material. The results show that the model is a good first approximation to the solution of solid-liquid phase change involving gravity-induced free convection.

Bain, R. L.↗

Real-time neuromorphic algorithms for inverse kinematics of redundant manipulators

The paper presents an efficient neuromorphic formulation to accurately solve the inverse kinematics problem for redundant manipulators. The approach involves a dynamical learning procedure based on a novel formalism in neural network theory: the concept of 'terminal' attractors. Topographically mapped terminal attractors are used to define a neural network whose synaptic elements can rapidly encapture the inverse kinematics transformations, and, subsequently generalize to compute joint-space coordinates required to achieve arbitrary end-effector configurations. Unlike prior neuromorphic implementations, this technique can also systematically exploit redundancy to optimize kinematic criteria, e.g., torque optimization. Simulations on 3-DOF and 7-DOF redundant manipulators, are used to validate the theoretical framework and illustrate its computational efficacy.

Barhen, Jacob↗

Radiative cooling in shock-heated hydrogen-helium plasmas

Axial and off-axis radiative cooling of cylindrical shock-heated hydrogen-helium plasmas is investigated theoretically and experimentally. The coupled fluid dynamic-radiative transfer equations are solved by a combination of approximation techniques aimed at simplifying the computation of the flux divergence term, namely, the quasi-isothermal approximation and the exponential approximation developed for the solid angle integration. The accuracy of the approximation schemes has been assessed and found acceptable for applying the methods to the rapid computation of the radiatively coupled flow problem. Radiative cooling experiments were conducted in a 6-inch annular arc accelerator shock tube (ANAA) for an initial pressure of 1 torr and shock speeds from 35 to 45 Km/sec. The results indicate that the lateral cooling is small compared with the axial cooling, and that better agreement is achieved between the data and the theoretical results by inclusion of the lateral temperature gradient.

Poon, P. T. Y.↗

Studies of flow separation and stalling on one- and two-element airfoils at low speeds

Research has been conducted on the nature of airfoil behavior at pre- and post-separated angles of attack. Detailed wind tunnel studies have been made of boundary layer and wake fields for the GA(W)-1 airfoil, and the airfoil with a 0.3 chord Fowler flap. Experimental data are compared with theoretical predictions from a multi-element viscous flow computer program. Theoretical predictions are reasonably accurate for unseparated flows, but have serious errors when separation is present. Some recent techniques for modeling post-separated flow behavior are discussed in light of the present experiments.

Seetharam, H. C.↗

A standard verification for rainfall estimation from remote platforms

Measurements of reflectivity at horizontal (Zh) and vertical (Zy) polarizations provide adequate information necessary to infer the two parameters of an exponential raindrop size distribution (No, Do) where the distribution is given by N(D) = Noexp(-3.67 D/Do). This distribution enables computation of water content or still air rainfall rates. The physical basis of the radar technique is outlined and illustrated theoretically, and experimental results, comparing radar derived rainfall rates with raingage and disdrometer measurements, are reviewed. The technique is useful for many meteorological and hydrological purposes, including ground truth measurements of rainfall rate over the ocean for comparison with satellite related observations.

Meitin, J. G., Jr.↗

A fast efficient implicit scheme for the gasdynamic equations using a matrix reduction technique

An efficient implicit finite-difference algorithm for the gasdynamic equations utilizing matrix reduction techniques is presented. A significant reduction in arithmetic operations is achieved without loss of the stability characteristics generality found in the Beam and Warming approximate factorization algorithm. Steady-state solutions to the conservative Euler equations in generalized coordinates are obtained for transonic flows and used to show that the method offers computational advantages over the conventional Beam and Warming scheme. Existing Beam and Warming codes can be retrofit with minimal effort. The theoretical extension of the matrix reduction technique to the full Navier-Stokes equations in Cartesian coordinates is presented in detail. Linear stability, using a Fourier stability analysis, is demonstrated and discussed for the one-dimensional Euler equations.

Barth, T. J.↗

A theoretical analysis of the electromagnetic environment of the AS330 super Puma helicopter external and internal coupling

Numerical techniques such as Finite Difference Time Domain (FDTD) computer programs, which were first developed to analyze the external electromagnetic environment of an aircraft during a wave illumination, a lightning event, or any kind of current injection, are now very powerful investigative tools. The program called GORFF-VE, was extended to compute the inner electromagnetic fields that are generated by the penetration of the outer fields through large apertures made in the all metallic body. Then, the internal fields can drive the electrical response of a cable network. The coupling between the inside and the outside of the helicopter is implemented using Huygen's principle. Moreover, the spectacular increase of computer resources, as calculations speed and memory capacity, allows the modellization structures as complex as these of helicopters with accuracy. This numerical model was exploited, first, to analyze the electromagnetic environment of an in-flight helicopter for several injection configurations, and second, to design a coaxial return path to simulate the lightning aircraft interaction with a strong current injection. The E field and current mappings are the result of these calculations.

Flourens, F.↗

Nonrecursive formulations of multibody dynamics and concurrent multiprocessing

Since the late 1980's, research in recursive formulations of multibody dynamics has flourished. Historically, much of this research can be traced to applications of low dimensionality in mechanism and vehicle dynamics. Indeed, there is little doubt that recursive order N methods are the method of choice for this class of systems. This approach has the advantage that a minimal number of coordinates are utilized, parallelism can be induced for certain system topologies, and the method is of order N computational cost for systems of N rigid bodies. Despite the fact that many authors have dismissed redundant coordinate formulations as being of order N(exp 3), and hence less attractive than recursive formulations, we present recent research that demonstrates that at least three distinct classes of redundant, nonrecursive multibody formulations consistently achieve order N computational cost for systems of rigid and/or flexible bodies. These formulations are as follows: (1) the preconditioned range space formulation; (2) penalty methods; and (3) augmented Lagrangian methods for nonlinear multibody dynamics. The first method can be traced to its foundation in equality constrained quadratic optimization, while the last two methods have been studied extensively in the context of coercive variational boundary value problems in computational mechanics. Until recently, however, they have not been investigated in the context of multibody simulation, and present theoretical questions unique to nonlinear dynamics. All of these nonrecursive methods have additional advantages with respect to recursive order N methods: (1) the formalisms retain the highly desirable order N computational cost; (2) the techniques are amenable to concurrent simulation strategies; (3) the approaches do not depend upon system topology to induce concurrency; and (4) the methods can be derived to balance the computational load automatically on concurrent multiprocessors. In addition to the presentation of the fundamental formulations, this paper presents new theoretical results regarding the rate of convergence of order N constraint stabilization schemes associated with the newly introduced class of methods.

Kurdila, Andrew J.↗

Optimal control of quaternion propagation errors in spacecraft navigation

Optimal control techniques are used to drive the numerical error (truncation, roundoff, commutation) in computing the quaternion vector to zero. The normalization of the quaternion is carried out by appropriate choice of a performance index, which can be optimized. The error equations are derived from Friedland's (1978) theoretical development, and a matrix Riccati equation results for the computation of the gain matrix. Simulation results show that a high precision of the order of 10 to the -12th can be obtained using this technique in meeting the q(T)q=1 constraint. The performance of the estimator in the presence of the feedback control that maintains the normalization, is studied.

Vathsal, S.↗

Theoretical Study of the Electric Dipole Moment Function of the CIO Molecule

The potential energy function and electric dipole moment function (EDMF) are computed for CIO Chi(sup 2)Pi using several different techniques to include electron correlation. The EDMF is used to compute Einstein coefficients, vibrational lifetimes, and dipole moments in higher vibrational levels. Remaining questions concerning the position of the maximum of the EDMF may be resolved through experimental measurement of dipole moments of higher vibrational levels. The band strength of the 1-0 fundamental transition is computed to be 12 +/- 2 /sq cm atm in good agreement with three experimental values, but larger than a recent value of 5 /sq cm atm determined from infrared heterodyne spectroscopy. The theoretical methods used include SCF, CASSCF, multireference singles plus doubles configuration interaction (MRCI) and contracted CI, coupled pair functional (CPF), and a modified version of the CPF method. The results obtained using the different methods are critically compared.

Pettersson, Lars G. M.↗

COMOC: Three dimensional boundary region variant, programmer's manual

The three-dimensional boundary region variant of the COMOC computer program system solves the partial differential equation system governing certain three-dimensional flows of a viscous, heat conducting, multiple-species, compressible fluid including combustion. The solution is established in physical variables, using a finite element algorithm for the boundary value portion of the problem description in combination with an explicit marching technique for the initial value character. The computational lattice may be arbitrarily nonregular, and boundary condition constraints are readily applied. The theoretical foundation of the algorithm, a detailed description on the construction and operation of the program, and instructions on utilization of the many features of the code are presented.

Orzechowski, J. A.↗

A semi-implicit dynamic relaxation algorithm for static nonlinear structural analysis

A semiimplicit dynamic relaxation technique for solution of the nonlinear structural equlibrium equation is presented. A previously presented basic transient response analysis algorithm is employed, permitting use of one solution method and one software module for both static and dynamic analyses. A theoretical comparison of the method with explicit dynamic relaxation techniques shows that it offers a substantially improved convergence property without additional computational overhead.

Park, K. C.↗

Detection of elliptical polarization and mode splitting in discrete Schumann resonance excitations

Elliptical polarization and mode splitting have been detected in the magnetic component of discrete, well defined Schumann resonance excitations. These ELF excitations, which are large electromagnetic transients of approximately 1 s duration, are called Q-bursts and typically occur every few minutes. They are believed to be the signature of the impulsive excitation of the earth-ionosphere cavity by ultra-large lightning currents. In this paper the magnetic polarization and spectral characteristics of four large Q-bursts are examined in detail using a new analysis technique. Two events display right-hand polarization and two display left-hand polarization. The theoretical polarization properties of the central and side multiplets of the Schumann resonances are used to define a local orthogonal coordinate system in the measurement frame in which these components may be separated. Maximum entropy spectrums computed separately for what are identified to be the central and side multiplets in this coordinate system show distinctly different eigenfrequencies for the lowest mode near 7.5 Hz. For the limited number of cases examined the magnitude of the line splitting detected using this technique is roughly 1.4-1.8 Hz, larger by nearly a factor of two than theoretical or observed values of the splitting previously reported. The frequencies of the side multiplets may lie either above or below the frequency of the central multiplet.

Sentman, D. D.↗

An integrated radiation physics computer code system.

An integrated computer code system for the semi-automatic and rapid analysis of experimental and analytic problems in gamma photon and fast neutron radiation physics is presented. Such problems as the design of optimum radiation shields and radioisotope power source configurations may be studied. The system codes allow for the unfolding of complex neutron and gamma photon experimental spectra. Monte Carlo and analytic techniques are used for the theoretical prediction of radiation transport. The system includes a multichannel pulse-height analyzer scintillation and semiconductor spectrometer coupled to an on-line digital computer with appropriate peripheral equipment. The system is geometry generalized as well as self-contained with respect to material nuclear cross sections and the determination of the spectrometer response functions. Input data may be either analytic or experimental.

Steyn, J. J.↗

Superpixel-Augmented Endmember Detection for Hyperspectral Images

Superpixels are homogeneous image regions comprised of several contiguous pixels. They are produced by shattering the image into contiguous, homogeneous regions that each cover between 20 and 100 image pixels. The segmentation aims for a many-to-one mapping from superpixels to image features; each image feature could contain several superpixels, but each superpixel occupies no more than one image feature. This conservative segmentation is relatively easy to automate in a robust fashion. Superpixel processing is related to the more general idea of improving hyperspectral analysis through spatial constraints, which can recognize subtle features at or below the level of noise by exploiting the fact that their spectral signatures are found in neighboring pixels. Recent work has explored spatial constraints for endmember extraction, showing significant advantages over techniques that ignore pixels relative positions. Methods such as AMEE (automated morphological endmember extraction) express spatial influence using fixed isometric relationships a local square window or Euclidean distance in pixel coordinates. In other words, two pixels covariances are based on their spatial proximity, but are independent of their absolute location in the scene. These isometric spatial constraints are most appropriate when spectral variation is smooth and constant over the image. Superpixels are simple to implement, efficient to compute, and are empirically effective. They can be used as a preprocessing step with any desired endmember extraction technique. Superpixels also have a solid theoretical basis in the hyperspectral linear mixing model, making them a principled approach for improving endmember extraction. Unlike existing approaches, superpixels can accommodate non-isometric covariance between image pixels (characteristic of discrete image features separated by step discontinuities). These kinds of image features are common in natural scenes. Analysts can substitute superpixels for image pixels during endmember analysis that leverages the spatial contiguity of scene features to enhance subtle spectral features. Superpixels define populations of image pixels that are independent samples from each image feature, permitting robust estimation of spectral properties, and reducing measurement noise in proportion to the area of the superpixel. This permits improved endmember extraction, and enables automated search for novel and constituent minerals in very noisy, hyperspatial images. This innovation begins with a graph-based segmentation based on the work of Felzenszwalb et al., but then expands their approach to the hyperspectral image domain with a Euclidean distance metric. Then, the mean spectrum of each segment is computed, and the resulting data cloud is used as input into sequential maximum angle convex cone (SMACC) endmember extraction.

Thompson, David R.↗

Lightning induced currents in aircraft wiring using low level injection techniques

Various techniques were studied to predict the transient current induced into aircraft wiring bundles as a result of an aircraft lightning strike. A series of aircraft measurements were carried out together with a theoretical analysis using computer modeling. These tests were applied to various aircraft and also to specially constructed cylinders installed within coaxial return conductor systems. Low level swept frequency CW (carrier waves), low level transient and high level transient injection tests were applied to the aircraft and cylinders. Measurements were made to determine the transfer function between the aircraft drive current and the resulting skin currents and currents induced on the internal wiring. The full threat lightning induced transient currents were extrapolated from the low level data using Fourier transform techniques. The aircraft and cylinders used were constructed from both metallic and CFC (carbon fiber composite) materials. The results show the pulse stretching phenomenon which occurs for CFC materials due to the diffusion of the lightning current through carbon fiber materials. Transmission Line Matrix modeling techniques were used to compare theoretical and measured currents.

Stevens, E. G.↗

Macromolecular recognition: Structural aspects of the origin of the genetic system

Theoretical simulation of prebiotic chemical processes is an invaluable tool for probing the phenomenon of evolution of life. Using computational and modeling techniques and guided by analogies from present day systems we, seek to understand the emergence of genetic apparatus, enzymatic catalysis and protein synthesis under prebiotic conditions. In one possible scenario, the RNA enzymatic reaction plays a key role in the emergence of the self-replicating and offers a clue to the onset of enzymatic catalysis prior to the existence of the protein biosynthetic machinery. Our ultimate goal is to propose a simple RNA segment which contains the specificity and catalytic activity of the contemporary RNA enzyme and which could emerge in a primordial chemical environment. To understand the mechanism of ribozyme catalyzed reactions, ab initio and semi-empirical (ZINDO) programs were used to investigate the reaction path of transphosphorylation. A special emphasis was placed on the possible catalytic and structural roles played by the coordinated magnesium cation. Both the inline and adjacent mechanisms of transphosphorylation have been studied. Another important aspect of this reaction is the identity of the functional groups which are essential for the acid base catalysis. The structural characteristics of the target helices, particularly a possible role of G center dot T pair, is under examination by molecular dynamics (MD) simulation technique. Modeling of the ancestral aminoacyl-tRNA synthetases (aRS) may provide important clues to the emergence of the genetic code and the protein synthetic machinery. Assuming that the catalytic function evolved before the elements of specific recognition of a particular amino acid, we are exploring the minimal structural requirements for the catalysis of tRNA aminoacylation. The molecular modeling system SYBYL was used for this study based on the high resolution crystallographic structures of the present day tyrosyl-adenylate:tyrRS and tRNA(Gln): ATP:glnRS complexes. The trinucleotide CCA of the 3'-end tRNA is placed into the active site pocket of tyrRS, based on the scheme of interaction between tRNA(Gln) and glnRS, and upon the stereochemistry of the tyrRS:tRNA:Tyr-AMP transition state. This provides a model of the non-specific recognition of a tRNA's 3'-end by an aRS, which might be similar to that of the ancestral aRS's. In the next step, modeling of the rest of the acceptor stem of tRNA (Tyr) with tyrRS is carried out.

Rein, Robert↗

Rotordynamic Instability Problems in High-Performance Turbomachinery 1996

The first rotordynamics workshop proceedings emphasized a feeling of uncertainty in predicting the stability of characteristics of high-performance turbomachinery. In the second workshop proceedings these uncertainties were reduced through programs established to systematically resolve problems, with emphasis on experimental validation of the forces that influence rotordynamics. In the third proceedings many programs for predicting or measuring forces and force coefficients in high-performance turbomachinery produced results. Data became available for designing new machines with enhanced stability characteristics or for upgrading existing machines. In the fourth proceedings there emerged trends towards a more unified view of rotordynamic instability problems and several encouraging new analytical developments. The fifth workshop supported the continuing trend toward a unified view with several new developments in the design and manufacture of new turbomachineries with enhanced stability characteristics along with new data and associated numerical/theoretical results. The sixth workshop report provided field experience and experimental results, and expanded the use of computational and control techniques with integration of damper, bearing, and eccentric seal operation results. The seventh workshop report provided field experiences, numerical, theoretical, and experimental results and control methods for seals, bearings, and dampers with some attention given to variable thermophysical properties and turbulence measurements, and introduction of two-phase flow results. In the present workshop, active magnetic bearings (AMB's) evolve into a new method of measuring rotordynamic coefficients with discussions on honeycomb seals, drop of magnetically supported rotors, seals, bearings and dampers with new data being reported. The intent of the workshop and this proceedings is to provide a continuing impetus for an understanding and resolution of these problems.

Source record↗