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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
NASA NTRS2009

Estimating Relative Positions of Outer-Space Structures

A computer program estimates the relative position and orientation of two structures from measurements, made by use of electronic cameras and laser range finders on one structure, of distances and angular positions of fiducial objects on the other structure. The program was written specifically for use in determining errors in the alignment of large structures deployed in outer space from a space shuttle. The program is based partly on equations for transformations among the various coordinate systems involved in the measurements and on equations that account for errors in the transformation operators. It computes a least-squares estimate of the relative position and orientation. Sequential least-squares estimates, acquired at a measurement rate of 4 Hz, are averaged by passing them through a fourth-order Butterworth filter. The program is executed in a computer aboard the space shuttle, and its position and orientation estimates are displayed to astronauts on a graphical user interface.

Balian, Harry↗
NASA NTRS1993

On the generation of secondary motion in circular to rectangular transition ducts

A computer code has been developed for the calculation of three-dimensional turbulent flow in circular to rectangular transition ducts. The governing equations are written in a general curvilinear coordinate system so that all singly-connected geometrical configurations can be accommodated. Turbulence closure is through the standard k-epsilon model or a Reynolds stress closure model. Computed results show qualitative agreement with experimental data. Pressure-induced secondary flow is the dominant mode through the transition, but turbulence-induced secondary flow becomes prominent some distance downstream of the transition.

Demuren, A. O.↗
NASA NTRS1976

Computation of turbulent boundary layers on curved surfaces, 1 June 1975 - 31 January 1976

An accurate method was developed for predicting effects of streamline curvature and coordinate system rotation on turbulent boundary layers. A new two-equation model of turbulence was developed which serves as the basis of the study. In developing the new model, physical reasoning is combined with singular perturbation methods to develop a rational, physically-based set of equations which are, on the one hand, as accurate as mixing-length theory for equilibrium boundary layers and, on the other hand, suitable for computing effects of curvature and rotation. The equations are solved numerically for several boundary layer flows over plane and curved surfaces. For incompressible boundary layers, results of the computations are generally within 10% of corresponding experimental data. Somewhat larger discrepancies are noted for compressible applications.

Wilcox, D. C.↗
NASA NTRS1986

Implicit methods for computing chemically reacting flow

The backward Euler scheme was used to solve a large system of inviscid flow and chemical rate equations in three spatial coordinates. The flow equations were integrated simultaneously in time by a conventional ADI factorization technique, then the species equations were solved by either simultaneous or successive techniques. The methods were evaluated in their efficiency and robustness for a hypersonic flow problem involving an aerobrake configuration. It was found that both implicit methods can effectively reduce the stiffness associated with the chemical production term and that the successive solution for the species was as stable as the simultaneous solution. The latter method is more economical because the computation time varies linearly with the number of species.

Li, C. P.↗
DOE OSTI2020

General formulas for adiabatic invariants in nearly periodic Hamiltonian systems

While it is well known that every nearly periodic Hamiltonian system possesses an adiabatic invariant, extant methods for computing terms in the adiabatic invariant series are inefficient. The most popular method involves the heavy intermediate calculation of a non-unique near-identity coordinate transformation, even though the adiabatic invariant itself is a uniquely defined scalar. Here, a less well-known method, developed by S. Omohundro, avoids calculating intermediate sequences of coordinate transformations but is also inefficient as it involves its own sequence of complex intermediate calculations. In order to improve the efficiency of future calculations of adiabatic invariants, we derive generally applicable, readily computable formulas for the first several terms in the adiabatic invariant series. To demonstrate the utility of these formulas, we apply them to charged-particle dynamics in a strong magnetic field and magnetic field-line dynamics when the field lines are nearly closed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1979

Numerical solution of the Navier-Stokes equations for laminar, transonic flows

An implicit finite difference solution of the Navier-Stokes equations yielded time histories of the transonic laminar flow development about a circular cylinder and NACA-0018 airfoil. Reynolds numbers ranged from those corresponding to purely laminar flow to those corresponding to significant turbulence in the boundary layer. Body thermal conditions of an adiabatic wall and a specified body temperature were considered. Versatility in treating arbitrary bodies was incorporated by using numerically generated, body-fitted coordinate transformations. Solution of the simultaneous difference equations for the dependent variables was obtained using an accelerated Gauss-Seidel iterative scheme. Computational results are presented in the form of velocity vector fields, Mach number contours, aerodynamic coefficients, heat transfer rates at the body surface, and body temperature distributions. Truncation analyses of first and second derivative difference approximations resulted in general criteria for numerically generated coordinates so that flow near a body is more accurately represented.

Turner, L., III↗
DOE OSTI2023

The five-point bootstrap

We study five-point correlation functions of scalar operators in d-dimensional conformal field theories. We develop a new approach to computing the five-point conformal blocks for exchanged primary operators of arbitrary spin by introducing a generalization of radial coordinates, using an appropriate ansatz, and perturbatively solving two quadratic Casimir differential equations. We then study five-point correlators 〈σσϵσσ〉 in the critical 3d Ising model. We truncate the operator product expansions (OPEs) in the correlator by including a finite number of primary operators with conformal dimension below a cutoff ∆ ⩽ ∆ cutoff . We then compute several OPE coefficients involving ϵ and two spinning operators by demanding that the truncated correlator approximately satisfies the crossing relation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Nonlinear Interface Reduction for Time-Domain Analysis of Hurty/Craig-Bampton Superelements with Frictional Contact

Virtual prototyping in engineering design relies today on modern numerical models of contacting structures with accurate resolution of interface mechanics, which strongly affect the system-level stiffness and energy dissipation due to frictional losses. High-fidelity modeling within the localized interfaces is required to resolve local quantities of interest that may drive design decisions. The high-resolution finite element meshes necessary to resolve inter-component stresses tend to be computationally expensive, particularly when the analyst is interested in response time histories. The Hurty/Craig-Bampton (HCB) transformation is a widely used method in structural dynamics for reducing the interior portion of a finite element model while having the ability to retain all nonlinear contact degrees of freedom (DOF) in physical coordinates. These models may still require many DOF to adequately resolve the kinematics of the interface, leading to inadequate reduction and computational savings. This study proposes a novel interface reduction method to overcome these challenges by means of system-level characteristic constraint (SCC) modes and properly orthogonal interface modal derivatives (POIMDs) for transient dynamic analyses. Both SCC modes and POIMDs are computed using the reduced HCB mass and stiffness matrices, which can be directly computed from many commercial finite element analysis software. Comparison of time history responses to an impulse-type load in a mechanical beam assembly indicate that the interface-reduced model correlates well with the HCB truth model. Localized features like slip and contact area are well-represented in the time domain when the beam assembly is loaded with a broadband excitation. The proposed method also yields reduced-order models with greater critical timestep lengths for explicit integration schemes.

42 ENGINEERING↗
NASA NTRS1986

Deformable Subreflector Computed by Geometric Optics

Distorted antenna surfaces forced to produce a uniform wave front. SUBFORMING employs geometric optics in determining subreflector coordinates to match main reflector surface with known distortions. Antenna with distorted paraboloidal reflecting surface forced to produce uniform wave front by using a Cassegrainian geometry with path-length-compensating subreflector. Program written in FORTRAN V for batch execution.

Williams, W. F.↗
NASA NTRS1986

Numerical Solution of the Incompressible Navier-Stokes Equations in Three-dimensional Generalized Curvilinear Coordinates

Numerically solving the incompressible Navier-Stokes equations is known to be time consuming and expensive. Testing of the INS3D computers code, which solves these equations with the use of the pseudocompressibility method, shows this method to be an efficient way to obtain the steady state solution. The effects of the waves introduced by the pseudocompressibility method are analyzed and criteria are set and tested for the choice of the pseudocompressibility parameter which governs the artificial sound speed. The code is tested using laminar flow over a two dimensional backward-facing step, and laminar flow over a two dimensional circular cylinder. The results of the computations over the backward-facing step are in excellent agreement with experimental results. The transient solution of the flow over the cylinder impulsively started from rest is in good agreement with experimental results. However, the computed frequency of periodic shedding of vortices behind the cylinder is not in agreement with the experimental value. For a three dimensional test case, computations were conducted for a cylinder end wall junction. The saddle point separation and horseshoe vortex system appear in the computed field. The solution also shows secondary vortex filaments which wrap around the cylinder and spiral up in the wake.

Rogers, S. E.↗
NASA NTRS2009

Lunar Pole Illumination and Communications Maps Computed from GSSR Elevation Data

A Digital Elevation Model of the lunar south pole was produced using Goldstone Solar System RADAR (GSSR) data obtained in 2006.12 This model has 40-meter horizontal resolution and about 5-meter relative vertical accuracy. This Digital Elevation Model was used to compute average solar illumination and Earth visibility with 100 kilometers of the lunar south pole. The elevation data were converted into local terrain horizon masks, then converted into lunar-centric latitude and longitude coordinates. The horizon masks were compared to latitude, longitude regions bounding the maximum Sun and Earth motions relative to the moon. Estimates of Earth visibility were computed by integrating the area of the region bounding the Earth's motion that was below the horizon mask. Solar illumination and other metrics were computed similarly. Proposed lunar south pole base sites were examined in detail, with the best site showing yearly solar power availability of 92 percent and Direct-To-Earth (DTE) communication availability of about 50 percent. Similar analysis of the lunar south pole used an older GSSR Digital Elevation Model with 600-meter horizontal resolution. The paper also explores using a heliostat to reduce the photovoltaic power system mass and complexity.

lunar base↗
NASA NTRS2010

Stereoscopic Machine-Vision System Using Projected Circles

A machine-vision system capable of detecting obstacles large enough to damage or trap a robotic vehicle is undergoing development. The system includes (1) a pattern generator that projects concentric circles of laser light forward onto the terrain, (2) a stereoscopic pair of cameras that are aimed forward to acquire images of the circles, (3) a frame grabber and digitizer for acquiring image data from the cameras, and (4) a single-board computer that processes the data. The system is being developed as a prototype of machine- vision systems to enable robotic vehicles ( rovers ) on remote planets to avoid craters, large rocks, and other terrain features that could capture or damage the vehicles. Potential terrestrial applications of systems like this one could include terrain mapping, collision avoidance, navigation of robotic vehicles, mining, and robotic rescue. This system is based partly on the same principles as those of a prior stereoscopic machine-vision system in which the cameras acquire images of a single stripe of laser light that is swept forward across the terrain. However, this system is designed to afford improvements over some of the undesirable features of the prior system, including the need for a pan-and-tilt mechanism to aim the laser to generate the swept stripe, ambiguities in interpretation of the single-stripe image, the time needed to sweep the stripe across the terrain and process the data from many images acquired during that time, and difficulty of calibration because of the narrowness of the stripe. In this system, the pattern generator does not contain any moving parts and need not be mounted on a pan-and-tilt mechanism: the pattern of concentric circles is projected steadily in the forward direction. The system calibrates itself by use of data acquired during projection of the concentric-circle pattern onto a known target representing flat ground. The calibration- target image data are stored in the computer memory for use as a template in processing terrain images. During operation on terrain, the images acquired by the left and right cameras are analyzed. The analysis includes (1) computation of the horizontal and vertical dimensions and the aspect ratios of rectangles that bound the circle images and (2) comparison of these aspect ratios with those of the template. Coordinates of distortions of the circles are used to identify and locate objects. If the analysis leads to identification of an object of significant size, then stereoscopicvision algorithms are used to estimate the distance to the object. The time taken in performing this analysis on a single pair of images acquired by the left and right cameras in this system is a fraction of the time taken in processing the many pairs of images acquired in a sweep of the laser stripe across the field of view in the prior system. The results of the analysis include data on sizes and shapes of, and distances and directions to, objects. Coordinates of objects are updated as the vehicle moves so that intelligent decisions regarding speed and direction can be made. The results of the analysis are utilized in a computational decision-making process that generates obstacle-avoidance data and feeds those data to the control system of the robotic vehicle.

Mackey, Jeffrey R.↗
NASA NTRS1986

An exact solution for orbit view-periods from a station on a tri-axial ellipsoidal planet

This paper presents the concise exact solution for predicting view-periods to be observed from a masked or unmasked tracking station on a tri-axial ellipsoidal surface. The new exact approach expresses the azimuth and elevation angles of a spacecraft in terms of the station-centered geodetic topocentric coordinates in an elegantly concise manner. A simple and efficient algorithm is developed to avoid costly repetitive computations in searching for neighborhoods near the rise and set times of each satellite orbit for each station. Only one search for each orbit is necessary for each station. Sample results indicate that the use of an assumed spherical earth instead of an 'actual' tri-axial ellipsoidal earth could introduce an error up to a few minutes in a view-period prediction for circular orbits of low or medium altitude. For an elliptical orbit of high eccentricity and long period, the maximum error could be even larger. The analytic treatment and the efficient algorithm are designed for geocentric orbits, but they should be applicable to interplanetary trajectories by an appropriate coordinates transformation at each view-period calculation. This analysis can be accomplished only by not using the classical orbital elements.

Tang, C. C. H.↗
DOE OSTI2024

Improving the five-point bootstrap

We present a new algorithm for the numerical evaluation of five-point conformal blocks in d-dimensions, greatly improving the efficiency of their computation. To do this we use an appropriate ansatz for the blocks as a series expansion in radial coordinates, derive a set of recursion relations for the unknown coefficients in the ansatz, and evaluate the series using a Padé approximant to accelerate its convergence. We then study the 〈σσϵσσ〉 correlator in the 3d critical Ising model by truncating the operator product expansion (OPE) and only including operators with conformal dimension below a cutoff ∆ ⩽ ∆cutoff. We approximate the contributions of the operators above the cutoff by the corresponding contributions in a suitable disconnected five-point correlator. Using this approach, we compute a number of OPE coefficients with greater accuracy than previous methods.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1984

3-D inelastic analysis methods for hot section components (base program)

A 3-D inelastic analysis methods program consists of a series of computer codes embodying a progression of mathematical models (mechanics of materials, special finite element, boundary element) for streamlined analysis of combustor liners, turbine blades, and turbine vanes. These models address the effects of high temperatures and thermal/mechanical loadings on the local (stress/strain) and global (dynamics, buckling) structural behavior of the three selected components. These models are used to solve 3-D inelastic problems using linear approximations in the sense that stresses/strains and temperatures in generic modeling regions are linear functions of the spatial coordinates, and solution increments for load, temperature and/or time are extrapolated linearly from previous information. Three linear formulation computer codes, referred to as MOMM (Mechanics of Materials Model), MHOST (MARC-Hot Section Technology), and BEST (Boundary Element Stress Technology), were developed and are described.

Wilson, R. B.↗
NASA NTRS1992

SKYMAP: Exploring the Universe in software

SKYMAP is a computer program which produces maps of arbitrary portions of the sky in a variety of projections and coordinate systems. Over the past ten years it has been used to produce finder charts for occultations by planets, display scan and image data from the Spacelab 2 Infrared Telescope, and make maps of fields for astronomical observations at X-ray, optical, infrared, and radio wavelengths. It can display multiple source catalogs, including the HST Guide Star Catalog, as well as solar system objects with astrometric accuracy using the JPL DE-130 ephemeris or tabulated positions. SKYMAP can be tuned to a specific task using an ASCII parameter file which controls how information is displayed on any Tektronix-compatible graphics display. The program contains a variety of interactive graphic and image processing features and has been ported to a variety of computer systems. A recent project visually demonstrates source density variation in various commonly-used all-sky catalogs.

Mink, D. J.↗
DOE OSTI2022

$\mathrm{BPS}$ coherent states and localization

We introduce coherent states averaged over a gauge group action to study correlators of half BPS states in $\mathcal{N}$= 4 SYM theory. The overlaps of these averaged coherent states are a generating function of correlators and can be written in terms of the Harish-Chandra-Itzykzon-Zuber (HCIZ) integral. We show that this formula immediately leads to a computation of the normalization of two point functions in terms of characters obtained originally in the work of Corley, Jevicki and Ramgoolam. We also find various generalizations for A n–1 quivers that follow directly from other solvable integrals over unitary groups. All of these can be computed using localization methods. When we promote the parameters of the generating function to collective coordinates, there is a dominant saddle that controls the effective action of these coherent states in the regime where they describe single AdS giant gravitons. We also discuss how to add open strings to this formulation. These will produce calculations that rely on correlators of matrix components of unitaries in the ensemble that is determined by the HCIZ integral to determine anomalous dimensions. We also discuss how sphere giants arise from Grassman integrals, how one gets a dominant saddle and how open strings are added in that case. The fact that there is a dominant saddle helps to understand how a 1/N expansion arises for open strings. We generalize the coherent state idea to study 1/4 and 1/8 BPS states as more general integrals over unitary groups.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Computing material volume fractions on a superimposed mesh as applied to Monte Carlo particle transport simulations

Here, we present a newly implemented ray tracing algorithm in OpenMC for efficiently computing material volume fractions on superimposed meshes in complex geometries. By firing rays along each coordinate direction through the geometry, the approach accumulates track-length data in each mesh element, thereby determining the fractional composition of each material. Scaling studies on three different models—a random tetrahedra configuration, the Frascati Neutron Generator ITER dose rate benchmark, and a stellarator design—show excellent parallel performance, with nearly linear speedup on modern multi-threaded and distributed-memory systems. An analysis of the residual error relative to high-resolution reference solutions demonstrated that under optimal conditions it decreases as 1/R, where R is the number of rays fired, making it straightforward to achieve user-prescribed accuracy. This new functionality enables practical, mesh-based approaches for detailed nuclear analyses in production Monte Carlo workflows without resorting to expensive, fully conformal or unstructured meshing.

Monte Carlo↗