Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computing methods”

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 1,711 records · Page 95

Electron Impact Ionization: A New Parameterization for 100 eV to 1 MeV Electrons

Low, medium and high energy electrons can penetrate to the thermosphere (90-400 km; 55-240 miles) and mesosphere (50-90 km; 30-55 miles). These precipitating electrons ionize that region of the atmosphere, creating positively charged atoms and molecules and knocking off other negatively charged electrons. The precipitating electrons also create nitrogen-containing compounds along with other constituents. Since the electron precipitation amounts change within minutes, it is necessary to have a rapid method of computing the ionization and production of nitrogen-containing compounds for inclusion in computationally-demanding global models. A new methodology has been developed, which has parameterized a more detailed model computation of the ionizing impact of precipitating electrons over the very large range of 100 eV up to 1,000,000 eV. This new parameterization method is more accurate than a previous parameterization scheme, when compared with the more detailed model computation. Global models at the National Center for Atmospheric Research will use this new parameterization method in the near future.

Fang, Xiaohua↗

Using NASTRAN to solve symmetric structures with nonsymmetric loads

A method for computation of reflective dihedral symmetry in symmetrical structures under nonsymmetric loads is described. The method makes it possible to confine the analysis to a half, a quarter, or an octagonal segment. The symmetry of elastic deformation is discussed, and antisymmetrical deformation is distinguished from nonsymmetrical deformation. Modes of deformation considered are axial, bending, membrane, and torsional deformation. Examples of one and two dimensional elements are presented and extended to three dimensional elements. The method of setting up a problem within NASTRAN is discussed. The technique is applied to a thick structure having quarter symmetry which was modeled with polyhedra and subjected to five distinct loads having varying degrees of symmetry.

Butler, T. G.↗

The NCOREL computer program for 3D nonlinear supersonic potential flow computations

An innovative computational technique (NCOREL) was established for the treatment of three dimensional supersonic flows. The method is nonlinear in that it solves the nonconservative finite difference analog of the full potential equation and can predict the formation of supercritical cross flow regions, embedded and bow shocks. The method implicitly computes a conical flow at the apex (R = 0) of a spherical coordinate system and uses a fully implicit marching technique to obtain three dimensional cross flow solutions. This implies that the radial Mach number must remain supersonic. The cross flow solutions are obtained by using type dependent transonic relaxation techniques with the type dependency linked to the character of the cross flow velocity (i.e., subsonic/supersonic). The spherical coordinate system and marching on spherical surfaces is ideally suited to the computation of wing flows at low supersonic Mach numbers due to the elimination of the subsonic axial Mach number problems that exist in other marching codes that utilize Cartesian transverse marching planes.

Siclari, M. J.↗

Performance of FFT methods in local gravity field modelling

Fast Fourier transform (FFT) methods provide a fast and efficient means of processing large amounts of gravity or geoid data in local gravity field modelling. The FFT methods, however, has a number of theoretical and practical limitations, especially the use of flat-earth approximation, and the requirements for gridded data. In spite of this the method often yields excellent results in practice when compared to other more rigorous (and computationally expensive) methods, such as least-squares collocation. The good performance of the FFT methods illustrate that the theoretical approximations are offset by the capability of taking into account more data in larger areas, especially important for geoid predictions. For best results good data gridding algorithms are essential. In practice truncated collocation approaches may be used. For large areas at high latitudes the gridding must be done using suitable map projections such as UTM, to avoid trivial errors caused by the meridian convergence. The FFT methods are compared to ground truth data in New Mexico (xi, eta from delta g), Scandinavia (N from delta g, the geoid fits to 15 cm over 2000 km), and areas of the Atlantic (delta g from satellite altimetry using Wiener filtering). In all cases the FFT methods yields results comparable or superior to other methods.

Forsberg, Rene↗

Theoretical study of the bonding of the first- and second-row transition-metal positive ions to methylene

The geometries of the molecules formed by the interaction of the first- and second-row transition-metal cations with methylene are optimized at the modified coupled-pair functional (MCPF) level of theory using large Gaussian basis sets, and their dissociation energies are computed employing both the MCPF and internally contracted averaged coupled-pair functional (ICACPF) methods. The computed binding energies are generally in good agreement with the available experimental results, although the calculations indicate that the experimental values for ScCH2(+), TiCH2(+), and NbCH2(+) are probably too large. The nature of the bonding in each case and trends in the bonding patterns across the transition-metal rows are discussed.

Bauschlicher, Charles W., Jr.↗

Improved Gaussian Beam-Scattering Algorithm

The localized model of the beam-shape coefficients for Gaussian beam-scattering theory by a spherical particle provides a great simplification in the numerical implementation of the theory. We derive an alternative form for the localized coefficients that is more convenient for computer computations and that provides physical insight into the details of the scattering process. We construct a FORTRAN program for Gaussian beam scattering with the localized model and compare its computer run time on a personal computer with that of a traditional Mie scattering program and with three other published methods for computing Gaussian beam scattering. We show that the analytical form of the beam-shape coefficients makes evident the fact that the excitation rate of morphology-dependent resonances is greatly enhanced for far off-axis incidence of the Gaussian beam.

Lock, James A.↗

A Spherical Earth Solution for TOA Lightning Location Retrieval

With the advent of high-speed digital computers, the method of chi square minimization is a highly practical means for analyzing a wide variety of (otherwise intractable) nonlinear inversion problems in applied mathematical physics. Little thought or effort is required to apply the chi square method to obtain quick and reasonable estimates of a solution, and the method offers a means to assess retrieval errors. Because the method is simple and practical it is sometimes hastily applied to problems that can be solved by formal analytic or quasi-analytic means. Presently, Global Atmospherics Inc. (GAI) finds the minimum of a chi square function to analyze time-of-arrival (TOA) and magnetic bearing data derived from the National Lightning Detection Network (NLDN); ellipsoidal Earth geometry is assumed. An analytic solution to this problem has not yet been attained, but the consideration and solving of less general problem statements might eventually lead to a final solution. In the present study, the problem of retrieving lightning ground-strike location on a spherical Earth surface using a network of 4 or more time-of-arrival (TOA) sensors is considered. It is shown that this problem has an analytic solution and therefore does not require the use of nonlinear estimation theory (such as the chi square method mentioned above). The mathematical robustness of the analytic solution is tested using computer-generated lightning sources and simulated TOA measurement errors. A quasi-analytic extension of the spherical Earth solution for an oblate spheroidal Earth geometry is considered in a related study. The incorporation of magnetic bearing information into these analytic solutions would lead to a general and elegant analytic retrieval scheme that would most likely replace the chi square estimation theory currently employed by Global Atmospherics Inc. (GAI).

Koshak, William J,↗

Compensating For GPS Ephemeris Error

Method of computing position of user station receiving signals from Global Positioning System (GPS) of navigational satellites compensates for most of GPS ephemeris error. Present method enables user station to reduce error in its computed position substantially. User station must have access to two or more reference stations at precisely known positions several hundred kilometers apart and must be in neighborhood of reference stations. Based on fact that when GPS data used to compute baseline between reference station and user station, vector error in computed baseline is proportional ephemeris error and length of baseline.

Wu, Jiun-Tsong↗

The FEM-2 design method

The FEM-2 parallel computer is designed using methods differing from those ordinarily employed in parallel computer design. The major distinguishing aspects are: (1) a top-down rather than bottom-up design process; (2) the design considers the entire system structure in terms of layers of virtual machines; and (3) each layer of virtual machine is defined formally during the design process. The result is a complete hardware/software system design. The basic design method is discussed and the advantages of the method are considered. A status report on the FEM-2 design is included.

Pratt, T. W.↗

Calculation of potential flow past airship bodies in yaw

An outline of Von Karman's method of computing the potential flow of airships in yaw by means of partially constant dipolar superposition on the axis of the body is followed by several considerations for beginning and end of the superposition. Then this method is improved by postulating a continuous, in part linearly variable dipolar superposition on the axis. The second main part of the report brings the calculation of the potential flow by means of sources and sinks, arranged on the surface of the airship body. The integral equation which must satisfy this surface superposition is posed, and the core reduced to functions developed from whole elliptical normal integrals. The functions are shown diagrammatically. The integration is resolvable by iteration. The consequence of the method is good. The formulas for computing the velocity on the surface and of the potential for any point conclude the report.

Lotz, I↗

Method to Measure Total Noise Temperature of a Wireless Receiver During Operation

An electromagnetic signal receiver and methods for determining the noise level and signal power in a signal of interest while the receiver is operating. In some embodiments, the signal of interest is a GPS signal. The receiver includes a noise source that provides a noise signal of known power during intervals while the signal of interest is observed. By measuring a signal-to-noise ratio for the signal of interest and the noise power in the signal of interest, the noise level and signal power of the signal of interest can be computed. Various methods of making the measurements and computing the power of the signal of interest are described. Applications of the system and method are described.

Young, Lawrence E.↗

Method of Characteristics Calculations and Computer Code for Materials with Arbitrary Equations of State and Using Orthogonal Polynomial Least Square Surface Fits

A numerical scheme using the method of characteristics to calculate the flow properties and pressures behind decaying shock waves for materials under hypervelocity impact is developed. Time-consuming double interpolation subroutines are replaced by a technique based on orthogonal polynomial least square surface fits. Typical calculated results are given and compared with the double interpolation results. The complete computer program is included.

Chang, T. S.↗

Recent Enhancements to the Development of CFD-Based Aeroelastic Reduced-Order Models

Recent enhancements to the development of CFD-based unsteady aerodynamic and aeroelastic reduced-order models (ROMs) are presented. These enhancements include the simultaneous application of structural modes as CFD input, static aeroelastic analysis using a ROM, and matched-point solutions using a ROM. The simultaneous application of structural modes as CFD input enables the computation of the unsteady aerodynamic state-space matrices with a single CFD execution, independent of the number of structural modes. The responses obtained from a simultaneous excitation of the CFD-based unsteady aerodynamic system are processed using system identification techniques in order to generate an unsteady aerodynamic state-space ROM. Once the unsteady aerodynamic state-space ROM is generated, a method for computing the static aeroelastic response using this unsteady aerodynamic ROM and a state-space model of the structure, is presented. Finally, a method is presented that enables the computation of matchedpoint solutions using a single ROM that is applicable over a range of dynamic pressures and velocities for a given Mach number. These enhancements represent a significant advancement of unsteady aerodynamic and aeroelastic ROM technology.

Silva, Walter A.↗

Root-cubing and general root-powering methods for finding the zeros of polynomials

Mathematical analysis technique generalizes a root squaring and root cubing method into a general root powering method. The introduction of partitioned polynomials into this general root powering method simplifies the coding of the polynomial transformations into input data suitable for processing by computer. The method includes analytic functions.

Bareiss, E. H.↗

Stability of matrix equations arising in electromagnetics

Many problems in electromagnetics can be formulated as a Fredholm equation of the first kind. In these problems it is desired to determine some function f(x) from its integral transform g(y). Except for special cases for which an explicit inverse operator is known, the equation is usually discretized by the method of moments into a matrix equation Af = g, which is then solved on a computer. A method of calculating the sensitivity of the solution of Af = g due to errors in either the vector g or the matrix A is developed. The usefulness of this method is demonstrated with problems from antenna theory, scattering theory, and remote sensing.

Klein, C.↗

Simulation of multicorrelated random processes using the FFT algorithm

A technique for the digital simulation of multicorrelated Gaussian random processes is described. This technique is based upon generating discrete frequency functions which correspond to the Fourier transform of the desired random processes, and then using the fast Fourier transform (FFT) algorithm to obtain the actual random processes. The main advantage of this method of simulation over other methods is computation time; it appears to be more than an order of magnitude faster than present methods of simulation. One of the main uses of multicorrelated simulated random processes is in solving nonlinear random vibration problems by numerical integration of the governing differential equations. The response of a nonlinear string to a distributed noise input is presented as an example.

Wittig, L. E.↗

A velocity-command controller for a VTOL aircraft

A technique is presented for calculating feedback and feedforward gain matrices that enable a VTOL aircraft to track input commands of forward and vertical velocity while maintaining acceptable responses to pilot inputs. Leverrier's algorithm is used for determining a set of state-variable, feedback gains that force the closed-loop poles and zeros of one pilot-input transfer function to pre-selected positions in the s-plane. This set of feedback gains is then used to calculate the feedback and feedforward gains for the velocity-command controller. The method is computationally attractive since the gains are determined by solving systems of linear, simultaneous equations. The method is used in a digital simulation of the CH-47 helicopter to control longitudinal dynamics.

Reid, G. F.↗