Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coordinate computing”

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,261 records · Page 70

Energy conservative SBP discretizations of the acoustic wave equation in covariant form on staggered curvilinear grids

In this work, we develop a numerical method for solving the acoustic wave equation in covariant form on staggered curvilinear grids in an energy conserving manner. The use of a covariant basis decomposition leads to a rotationally invariant scheme that outperforms a Cartesian basis decomposition on rotated grids. The discretization is based on high order Summation-By-Parts (SBP) operators and preserves both symmetry and positive definiteness of the contravariant metric tensor. To improve accuracy and decrease computational cost, we also derive a modified discretization of the metric tensor that leads to a conditionally stable discretization. Bounds are derived that yield a point-wise condition that can be evaluated to check for stability of the modified discretization. This condition shows that the interpolation operators should be constructed such that their norm is close to one.

97 MATHEMATICS AND COMPUTING↗

Crystallographic Evidence for Formation of M 6 O 8 /M 6 O 9 f-Element Clusters in Hydrolysis Reactions

In this contribution, we have conducted an examination of crystallographic data which includes 184 4f- and 5f-structures containing the M 6 O 8 or M 6 O 9 moiety. Since much of the work in this field has actually been done to explore the formation and properties of metal organic frameworks (MOFs), we hope that this collection can be used by actinide scientists to further investigate structural and mechanistic trends and gain insight into the formation of these compounds. Here we discuss the general trends in coordination, ligands used, and bonding within the M 6 O 8 or M 6 O 9 core in the hope that further computational work can benefit.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Two-State Reactivity in Iron-Catalyzed Alkene Isomerization Confers σ-Base Resistance

A low-coordinate, high spin (S = 3/2) organometallic iron(I) complex is a catalyst for the isomerization of alkenes. A combination of experimental and computational mechanistic studies supports a mechanism in which alkene isomerization occurs by the allyl mechanism. Importantly, while substrate binding occurs on the S = 3/2 surface, oxidative addition to an η 1 -allyl intermediate only occurs on the S = 1/2 surface. In conclusion, since this spin state change is only possible when the alkene substrate is bound, the catalyst has high immunity to typical σ-base poisons due to the antibonding interactions of the high spin state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Three dimensional nozzle-exhaust flow field analysis by a reference plane technique.

A numerical method based on reference plane characteristics has been developed for the calculation of highly complex supersonic nozzle-exhaust flow fields. The difference equations have been developed for three coordinate systems. Local reference plane orientations are employed using the three coordinate systems concurrently thus catering to a wide class of flow geometries. Discontinuities such as the underexpansion shock and contact surfaces are computed explicitly for nonuniform vehicle external flows. The nozzles considered may have irregular cross-sections with swept throats and may be stacked in modules using the vehicle undersurface for additional expansion. Results are presented for several nozzle configurations.

Dash, S. M.↗

Large-deformation modal coordinates for nonrigid vehicle dynamics

The derivation of minimum-dimension sets of discrete-coordinate and hybrid-coordinate equations of motion of a system consisting of an arbitrary number of hinge-connected rigid bodies assembled in tree topology is presented. These equations are useful for the simulation of dynamical systems that can be idealized as tree-like arrangements of substructures, with each substructure consisting of either a rigid body or a collection of elastically interconnected rigid bodies restricted to small relative rotations at each connection. Thus, some of the substructures represent elastic bodies subjected to small strains or local deformations, but possibly large gross deformations, in the hybrid formulation, distributed coordinates referred to herein as large-deformation modal coordinates, are used for the deformations of these substructures. The equations are in a form suitable for incorporation into one or more computer programs to be used as multipurpose tools in the simulation of spacecraft and other complex electromechanical systems.

Likins, P. W.↗

SICLOPS - A system of computer programs for rectified mapping of airborne scanner imagery

This paper describes SICLOPS, a system of computer programs, which has been developed to perform geometrically corrected line-by-line mapping of airborne multispectral scanner data to ground coordinates and to estimate ground areas. As inputs, SICLOPS accepts aircraft attitude and positional information furnished by ancillary aircraft equipment, plus ground truth information giving the location, in map grid coordinates, of a number of recognizable landmarks on the ground. The geometric correction and mapping procedure locates the scan lines, or the pixels on each line, in terms of map grid coordinates. The area estimation procedure gives ground area for each pixel or for a predesignated parcel specified in map grid coordinates. The results of exercising the system with simulated data showed the uncorrected video and corrected imagery and produced area estimates accurate to better than 99.7%. Real input data, though prepared, have not yet been processed by the system.

Spencer, M. M.↗

A general method for calculating three-dimensional compressible laminar and turbulent boundary layers on arbitrary wings

The method described utilizes a nonorthogonal coordinate system for boundary-layer calculations. It includes a geometry program that represents the wing analytically, and a velocity program that computes the external velocity components from a given experimental pressure distribution when the external velocity distribution is not computed theoretically. The boundary layer method is general, however, and can also be used for an external velocity distribution computed theoretically. Several test cases were computed by this method and the results were checked with other numerical calculations and with experiments when available. A typical computation time (CPU) on an IBM 370/165 computer for one surface of a wing which roughly consist of 30 spanwise stations and 25 streamwise stations, with 30 points across the boundary layer is less than 30 seconds for an incompressible flow and a little more for a compressible flow.

Cebeci, T.↗

A finite element algorithm for sound propagation in axisymmetric ducts containing compressible mean flow

An accurate mathematical model for sound propagation in axisymmetric aircraft engine ducts with compressible mean flow is reported. The model is based on the usual perturbation of the basic fluid mechanics equations for small motions. Mean flow parameters are derived in the absence of fluctuating quantities and are then substituted into the equations for the acoustic quantities which were linearized by eliminating higher order terms. Mean swirl is assumed to be zero from the restriction of axisymmetry. A linear rectangular serendipity element is formulated from these equations using a Galerkin procedure and assembled in a special purpose computer program in which the matrix map for a rectangular mesh was specifically coded. Representations of the fluctuating quantities, mean quantities and coordinate transformations are isoparametric. The global matrix is solved by foreward and back substitution following an L-U decomposition with pivoting restricted internally to the blocks. Results from the model were compared with results from several alternative analyses and yielded satisfactory agreement.

Abrahamson, A. L.↗

A high Reynolds number numerical solution of the Navier-Stokes equations in stream function-vorticity form

Problems encoutered in investigations of high Reynolds number, incompressible flow are reviewed. A numerical solution computer program was modified to solve the stream function-vorticity form of the Navier-Stokes equations. Using a body fitted coordinate system with a U-shaped outer boundary, a simulation of incompressible flow at a Reynolds number of one million and a body angle of attack of zero was achieved.

Bearden, J. H.↗

A calculation procedure for viscous flow in turbomachines, volume 3

A method for analyzing the nonadiabatic viscous flow through turbomachine blade passages was developed. The field analysis is based upon the numerical integration of the full incompressible Navier-Stokes equations, together with the energy equation on the blade-to-blade surface. A FORTRAN IV computer program was written based on this method. The numerical code used to solve the governing equations employs a nonorthogonal boundary fitted coordinate system. The flow may be axial, radial or mixed and there may be a change in stream channel thickness in the through-flow direction. The inputs required for two FORTRAN IV programs are presented. The first program considers laminar flows and the second can handle turbulent flows. Numerical examples are included to illustrate the use of the program, and to show the results that are obtained.

Khalil, I.↗

Grid and metric generation on the assembly of locally bi-quadratic coordinate transformations

The generation of metric coefficients of the coordinate transformation from a generally curved-sided domain boundary to the unit square (cube) is required for efficient solution algorithms in computational fluid mechanics. An algebraic procedure is presented for establishment of these data on the union of arbitrarily selected sub-domains of the global solution domain. A uniformly smooth progression of grid refinement is readily generated, including multiple specification of refined grids for a given macro-element domain discretization. The procedure is illustrated as generally applicable to non-simply connected domains in two and three dimensions.

Baker, A. J.↗

Study of the time-differentiated particle flux density at various distances from EAS axis

The EAS time structure is studied using the enlarged EAS array of the Moscow State University. The time measurements are made using 22 scintillators which form 13 rectanges of 180x190 sq m size covering the entire array area. The array is triggered by a signal of 4-fold coincidences of the pulses from the detectors forming each of the rectangles. The data were obtained during 2200 hours of the array operation in 1984. A total of 816 showers, to which at least 14 of 22 scintillator detectors responded, were selected among all those detected. The coordinates of the EAS axis in the observation plane and the EAS sizes were determined by the maximum likelihood method using a computer on the assumption that the electron LDF is the NKG form. A total of 492 showers in the interval of EAS size Ne = 5x10 to the 6th power - 2x10 to the 8th power (N bar e = 1.7x 10 to the 7th power) with zenith angles theta or = 45 deg and axes within the array are analyzed.

Atrashkevich, V. G.↗

Digital fringe reduction techniques applied to the measurement of three-dimensional transonic flow fields

A system for digital inteferogram evaluation based on an image-processing system connected to a host computer has been implemented. The system supports one- and two-dimensional interferogram evaluations. Interferograms are digitized, enhanced, and segmented. The fringe coordinates are extracted, and the fringes are represented as polygon data structures. Fringe numbering and fringe interpolation modules are implemented. The system supports editing and interactive features as well as graphic visualization. An application of the system to the evaluation of double exposure interferograms from the transonic flow field around a helicopter blade and the reconstruction of the 3-D flow field are given.

Becker, F.↗

Microwave backscattering from an anisotropic soybean canopy

Electromagnetic backscattering from a soybean canopy is modeled in the L band region of the spectrum. Mature soybean plants are taken as an ensemble of leaves and stems which are represented by lossy dielectric disks and rods respectively. Field data indicated that leaves and stems are not distributed uniformly in the azimuth coordinate. The plant has a tendency to grow out into the area between the rows. The effects on backscattered radar waves was computed by the distorted Born approximation. Results for look directions along the rows and perpendicular to the rows show that only a modest difference occurs in the L band frequency range. The use of another nonuniform distribution, different from those observed experimentally, results in a significant effect due to vegetation asymmetry.

Lang, R. H.↗

Application of digital interferogram evaluation techniques to the measurement of 3-D flow fields

A system for digitally evaluating interferograms, based on an image processing system connected to a host computer, was implemented. The system supports one- and two-dimensional interferogram evaluations. Interferograms are digitized, enhanced, and then segmented. The fringe coordinates are extracted, and the fringes are represented as polygonal data structures. Fringe numbering and fringe interpolation modules are implemented. The system supports editing and interactive features, as well as graphic visualization. An application of the system to the evaluation of double exposure interferograms from the transonic flow field around a helicopter blade and the reconstruction of the three dimensional flow field is given.

Becker, Friedhelm↗

A review of modelling techniques for the open and closed-loop dynamics of large space systems

This paper reviews the steps in the development of mathematical models that can be used to simulate the in-orbit dynamic behavior of large flexible systems. A general continuum formulation is compared with the hybrid coordinate formulation and also a finite element representation of the total system. A review of structural analysis routines emphasizes the use of computer generated graphics to help understand the different modal elastic shape functions of complex systems. Numerical techniques employed to synthesize shape and attitude control laws are summarized. Finally, the modeling of environmental disturbance torques due to the interaction of solar radiation pressure on vibrating and thermally deflected systems is discussed.

Bainum, Peter M.↗

A direct algorithm for solution of incompressible three-dimensional unsteady Navier-Stokes equations

A direct, implicit, numerical solution algorithm, with second-order accuracy in space and time, is constructed for the three-dimensional unsteady incompressible Navier-Stokes equations formulated in terms of velocity and vorticity, using generalized orthogonal coordinates to achieve the accurate solution of complex viscous flow configurations. A numerically stable, efficient, direct inversion procedure is developed for the computationally intensive divergence-curl elliptic velocity problem. This overdetermined partial differential operator is first formulated as a uniquely determined, nonsingular matrix-vector problem; this aspect of the procedure is a unique feature of the present analysis. The three-dimensional vorticity-transport equation is solved by a modified factorization technique which completely eliminates the need for any block-matrix inversions and only scalar tridiagonal matrices need to be inverted. The method is applied to the test problem of the three-dimensional flow within a shear-driven cubical box. Coherent streamwise vortex structures are observed within the steady-state flow at Re = 100.

Osswald, G. A.↗

Algorithm For Hypersonic Flow In Chemical Equilibrium

Implicit, finite-difference, shock-capturing algorithm calculates inviscid, hypersonic flows in chemical equilibrium. Implicit formulation chosen because overcomes limitation on mathematical stability encountered in explicit formulations. For dynamical portion of problem, Euler equations written in conservation-law form in Cartesian coordinate system for two-dimensional or axisymmetric flow. For chemical portion of problem, equilibrium state of gas at each point in computational grid determined by minimizing local Gibbs free energy, subject to local conservation of molecules, atoms, ions, and total enthalpy. Major advantage: resulting algorithm naturally stable and captures strong shocks without help of artificial-dissipation terms to damp out spurious numerical oscillations.

Palmer, Grant↗