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 433 records · Page 24

Parallel Wavefront Analysis for a 4D Interferometer

This software provides a programming interface for automating data collection with a PhaseCam interferometer from 4D Technology, and distributing the image-processing algorithm across a cluster of general-purpose computers. Multiple instances of 4Sight (4D Technology s proprietary software) run on a networked cluster of computers. Each connects to a single server (the controller) and waits for instructions. The controller directs the interferometer to several images, then assigns each image to a different computer for processing. When the image processing is finished, the server directs one of the computers to collate and combine the processed images, saving the resulting measurement in a file on a disk. The available software captures approximately 100 images and analyzes them immediately. This software separates the capture and analysis processes, so that analysis can be done at a different time and faster by running the algorithm in parallel across several processors. The PhaseCam family of interferometers can measure an optical system in milliseconds, but it takes many seconds to process the data so that it is usable. In characterizing an adaptive optics system, like the next generation of astronomical observatories, thousands of measurements are required, and the processing time quickly becomes excessive. A programming interface distributes data processing for a PhaseCam interferometer across a Windows computing cluster. A scriptable controller program coordinates data acquisition from the interferometer, storage on networked hard disks, and parallel processing. Idle time of the interferometer is minimized. This architecture is implemented in Python and JavaScript, and may be altered to fit a customer s needs.

Rao, Shanti R.↗

Coronal Flux Ropes and their Interplanetary Counterparts

We report on a study comparing coronal flux ropes inferred from eruption data with their interplanetary counterparts constructed from in situ data. The eruption data include the source region magnetic field, post-eruption arcades, and coronal mass ejections (CMEs). Flux ropes were fit to the interplanetary CMEs (ICMEs) considered for the 2011 and 2012 Coordinated Data Analysis Workshops (CDAWs). We computed the total reconnected flux involved in each of the associated solar eruptions and found it to be closely related to flare properties, CME kinematics, and ICME properties. By fitting flux ropes to the white-light coronagraph data, we obtained the geometric properties of the flux ropes and added magnetic properties derived from the reconnected flux. We found that the CME magnetic field in the corona is significantly higher than the ambient magnetic field at a given heliocentric distance. The radial dependence of the flux rope magnetic field strength is faster than that of the ambient magnetic field. The magnetic field strength of the coronal flux rope is also correlated with that in interplanetary flux ropes constructed from in situ data, and with the observed peak magnetic field strength in ICMEs. The physical reason for the observed correlation between the peak field strength in MCs is the higher magnetic field content in faster coronal flux ropes and ultimately the higher reconnected flux in the eruption region. The magnetic flux ropes constructed from the eruption data and coronagraph observations provide a realistic input that can be used by various models to predict the magnetic properties of ICMEs at Earth and other destination in the heliosphere.

Gopalswamy, N.↗

Americium Oxalate: An Experimental and Computational Investigation of Metal–Ligand Bonding

A novel actinide-containing coordination polymer, [Am(C 2 O 4 )(H 2 O) 3 Cl] (Am-1), has been synthesized and structurally characterized. The crystallographic analysis reveals that the structure is two-dimensional and comprised of pseudo-dimeric Am 3+ nodes that are bridged by oxalate ligands to form sheets. Each metal center is nine-coordinate, forming a distorted capped square antiprism geometry with a C 1 symmetry, and features bound oxalate, aqua, and chloro ligands. The Am 3+ –ligand bonds were probed computationally using the quantum theory of atoms in molecules nd natural localized molecular orbital approaches to investigate the underlying mechanisms and hybrid atomic orbital contributions therein. The analyses indicate that the bonds within Am-1 are predominantly ionic and the 5f shell of the Am 3+ metal centers does not add a significant covalent contribution to the bonds. Further, our bonding assessment is supported by measurements on the optical properties of Am-1 using diffuse reflectance and photoluminescence spectroscopies. The position of the principal absorption band at 507 nm ( 5 L 6' ← 7 F 0' ) is notable because it is consistent with previously reported americium oxalate complexes in solution, indicating similarities in the electronic structure and ionic bonding. Compound Am-1 is an active phosphor, featuring strong bright-blue oxalate-based luminescence with no evidence of metal-centered emission.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

How to determine spiral bevel gear tooth geometry for finite element analysis

An analytical method was developed to determine gear tooth surface coordinates of face milled spiral bevel gears. The method combines the basic gear design parameters with the kinematical aspects for spiral bevel gear manufacturing. A computer program was developed to calculate the surface coordinates. From this data a 3-D model for finite element analysis can be determined. Development of the modeling method and an example case are presented.

Handschuh, Robert F.↗

Early orbit determination scheme for the Juno space vehicle

The present status (1959-60) of the ABMA early orbit determination scheme for the Juno space vehicle is described. The design of the scheme was essentially influenced by three factors: its origin centered around postflight vehicle analysis, the flexibility required for accepting varied kinds of data, and the potential of highspeed computers concentrated in ABMA's Computation Laboratory. The scheme is outlined, the tracking stations participating are shortly described, and the physical layout and computational equipment are briefly mentioned. The evaluation techniques used for the various steps of the orbit determination are discussed in detail. Appendices describe briefly the principal coordinate systems used and characteristics of the digital computer programs.

Fletcher Kurtz↗

Spectroscopic characterization of Mn 2+ and Cd 2+ coordination to phosphorothioates in the conserved A9 metal site of the hammerhead ribozyme

Phosphorothioate modifications have widespread use in the field of nucleic acids. As substitution of sulfur for oxygen can alter metal coordination preferences, the phosphorothioate metal-rescue experiment is a powerful method for identifying metal coordination sites that influence specific properties in a large RNAs. The A9/G10.1 metal binding site of the hammerhead ribozyme (HHRz) has previously been shown to be functionally important through phosphorothioate rescue experiments. While an A9-S Rp substitution is inhibitory in Mg 2+ , thiophilic Cd 2+ rescues HHRz activity. Mn 2+ is also often used in phosphorothioate metal-rescue studies but does not support activity for the A9-S Rp HHRz. Here, we use EPR, electron spin-echo envelope modulation (ESEEM), and X-ray absorption spectroscopic methods to directly probe the structural consequences of Mn 2+ and Cd 2+ coordination to R p and S p phosphorothioate modifications at the A9/G10.1 site in the truncated hammerhead ribozyme (tHHRz). The results demonstrate that while Cd 2+ does indeed bind to S in the thio-substituted ligand, Mn 2+ coordinates to the non–sulfur oxo group of this phosphorothioate, regardless of isomer. Computational models demonstrate the energetic preference of Mn—O over Mn—S coordination in metal-dimethylthiophosphate models. In the case of the tHHRz, the resulting Mn 2+ coordination preference of oxygen in either R p or S p A9 phosphorothioates differentially tunes catalytic activity, with Mn—O coordination in the A9-S Rp phosphorothioate enzyme being inhibitory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

TOMCAT - A code for numerical generation of boundary-fitted curvilinear coordinate systems on fields containing any number of arbitrary two-dimensional bodies

A method for automatic generation of boundary-fitted curvilinear coordinate systems, where the transformed coordinates are solutions of an elliptic differential system in the physical plane, and where the coordinate lines are coincident with all boundaries of a general multiply-connected, two-dimensional region containing any number of arbitrarily shaped bodies, and is described along with a suitable computer code for implementing the method. Any partial differential system can be solved on the boundary-fitted coordinate system by appropriate transformations. The transformed equations are approximated by finite differences and solved numerically in the transformed plane. All computations, whether for generating coordinate system or then solving the transformed equations, can be done on a rectangular field with square mesh with no interpolation required on the boundaries. The physical boundaries may even be time-dependent.

Thompson, J. F.↗

Computational fluid dynamics for propulsion technology: Geometric grid visualization in CFD-based propulsion technology research

The coordination is examined of necessary resources, facilities, and special personnel to provide technical integration activities in the area of computational fluid dynamics applied to propulsion technology. Involved is the coordination of CFD activities between government, industry, and universities. Current geometry modeling, grid generation, and graphical methods are established to use in the analysis of CFD design methodologies.

Ziebarth, John P.↗

Constructing field-aligned coordinate systems for gyrokinetic simulations of tokamaks in X-point geometries

Structures in tokamak plasmas are elongated along the direction of the magnetic field and short in the directions perpendicular to the magnetic field. Many tokamak simulation codes take advantage of this by using a field-aligned coordinate system. However, field-aligned coordinate systems have a coordinate singularity at magnetic X-points where the poloidal magnetic field vanishes, which makes it difficult to use field-aligned coordinate systems when simulating the core and scrape-off layer simultaneously. Here, we present an algorithm for grid generation and computing geometric quantities in a standard field-aligned coordinate system that avoids the singularity and allows one to conduct two-dimensional gyrokinetic axisymmetric simulations in X-point geometries. Convergence tests of advection, boundary value problems and geometric quantities all show greater than first-order convergence even in the vicinity of the X-point. We also demonstrate the geometric consistency of our algorithm with an example simulation of the spherical tokamak for energy production, which shows machine-precision particle conservation.

fusion plasma↗

Quasiparticle self-consistent GW band structures and phase transitions of LiAlO 2 in tetrahedrally and octahedrally coordinated structures

In this work, a first-principles computational study is presented of various phases of LiAlO 2 . The relative total energies and equations of state of the α, β, and γ phases are determined after structural relaxation of each phase. The β and γ tetrahedral phases are found to be very close in energy and lattice volume with the γ phase having the lowest energy. The octahedral α phase is a high-pressure phase and the transition pressure from the γ and β phases to α is determined to be about 1 GPa. The electronic band structures of each phase at their own equilibrium volume are determined using the quasiparticle self-consistent (QS) GW method as well as using the 0.8 Σ approach in which the QS GW self-energy is reduced by a factor of 0.8 to correct for the underscreening of W in QS GW. The effective masses of the band edges and the nature of the band gaps are presented. The lowest energy γ phase is found to have a pseudodirect gap of 7.69 eV. The gap is direct at Γ but corresponds to a dipole forbidden transition. The imaginary part of the dielectric function and the absorption coefficient are calculated in the long-wavelength limit and the random phase approximation, without local field or electron-hole interaction effects for each phase, and their anisotropies are discussed. Si doping on the Al site is investigated as a possible n-type dopant in γ-LiAlO 2 using a 128-atom supercell corresponding to 3.125% Si on the Al sublattice in the generalized gradient approximation and a smaller 16-atom cell with 25% Si in the QS GW approximation. The Si is found to significantly perturb the conduction band and lower the gap but a clearly separated deep donor defect level is not found. However, the donor binding energy is still expected to be relatively deep, on the order of a few tenths eV in the hydrogenic effective mass approximation.

36 MATERIALS SCIENCE↗

The Density Matrix of H20 - N2 In the Coordinate Representation: A Monte Carlo Calculation of the Far-Wing Line Shape

The far-wing line shape theory within the binary collision and quasistatic framework has been developed using the coordinate representation. Within this formalism, the main computational task is the evaluation of multidimensional integrals whose variables are the orientational angles needed to specify the initial and final positions of the system during transition processes. Using standard methods, one is able to evaluate the 7-dimensional integrations required for linear molecular systems, or the 7-dimensional integrations for more complicated asymmetric-top (or symmetric-top) molecular systems whose interaction potential contains cyclic coordinates. In order to obviate this latter restriction on the form of the interaction potential, a Monte Carlo method is used to evaluate the 9-dimensional integrations required for systems consisting of one asymmetric-top (or symmetric-top) and one linear molecule, such as H20-N2. Combined with techniques developed previously to deal with sophisticated potential models, one is able to implement realistic potentials for these systems and derive accurate, converged results for the far-wing line shapes and the corresponding absorption coefficients. Conversely, comparison of the far-wing absorption with experimental data can serve as a sensitive diagnostic tool in order to obtain detailed information on the short-range anisotropic dependence of interaction potentials.

Ma, Q.↗

Aerostructural Vortical Flow Interactions with Applications to F/A-18 and F-117 Tail Buffet

The buffet response of the flexible twin-tail configuration of the F/A-18 and a generic F-111 aircraft are computationally simulated and experimentally validated. The problem is a multidisciplinary one which requires the sequential solution of three sets of equations on a multi-block grid structure. The first set is the unsteady, compressible, full Navier-Stokes equations. The second set is the aeroelastic equations for bending and torsional twin-tail responses. The third set is the grid-displacement equations which are used to update the grid coordinates due to the tail deflections. The computational models consist of a 76 deg. swept back, sharp edged delta wing of aspect ratio of one and a swept-back F/A-18 or F-117 twin-tail. The configuration is pitched at 30 deg. angle-of-attack. The problem is solved for the initial flow conditions with the twin tails kept rigid. Next, the aeroelastic equations of the tails are turned on along with the grid-displacement equations to solve for the bending and torsional tails responses due to the unsteady loads produced by the vortex breakdown flow of the leading-edge vortex cores of the delta wing. Several spanwise locations of the twin tails are investigated. The computational results are validated using several existing experimental data.

Kandil, Osama A.↗

A three-dimensional turbulent compressible subsonic duct flow analysis for use with constructed coordinate systems

An approximate analysis, applicable to nonorthogonal coordinate systems having a curved centerline and planar transverse coordinate surfaces normal to the centerline, is presented for computation of three-dimensional subsonic flow in straight and curved diffusers. The formulation is intended to facilitate the use of constructed coordinates in circumstances where it is difficult to maintain smooth behavior in higher derivatives; the use of local Cartesian variables and fluxes leads to governing equations which require only first derivatives of the coordinate transformation. The analysis is applied to a particular family of duct and diffuser geometries having curved centerlines and superelliptic cross sections. Qualitative agreement with experimental measurements is observed with regard to streamwise vortices and distortion of the primary flow.

Levy, R.↗

Computation and Validation of Fluid/Structure Twin Tail Buffet Response

The buffet response of the flexible twin-tail/delta wing configuration, a multidisciplinary problem, is solved using three sets of equations on a multi-block grid structure. The first set is the unsteady, compressible, full Navier-Stokes equations which are used for obtaining the flow-field vector and the aerodynamic loads on the twin tails. The second set is the coupled aeroelastic equations which are used for obtaining the bending and torsional deflections of the twin tails. The third set is the grid-displacement equations which are used for updating the grid coordinates due to the tail deflections. The computational model is similar to the one used by Washburn et. al. which consists of a delta wing of aspect ratio one and twin tails with taper ratios of 0.23. The vortex of the twin tails are located at the wing trailing edge. The configuration is pitched at 30 deg angle of attack, and the freestream Mach number and Reynolds number are 0.3 and 1.25 million, respectively. With the twin tails fixed as rigid surfaces, the problem is solved for the initial flow conditions. Next, the problem is solved for the twin tail response for uncoupled bending and torsional vibrations due to the unsteady loads produced by the vortex breakdown flow of the leading-edge vortex cores. The configuration is investigated for three spanwise positions of the twin tails; inboard, midspan and outboard locations. The computational results are validated and are in very good agreement with the experimental data of Washburn, et. al.

Massey, Steven J.↗

Interactive LORAN-C to geographic and geographic-to-LORAN-C computation

The LORAN program is stored in CMS disk files for use by Avionics Engineering Center terminal users. A CMS EXEC file named LORAN controls program operation. The user types LORAN and the program then prompts for data input and produces output on the terminal. The FORTRAN program refers to a disk file of LORAN master data giving station locations, coding delays, repetition rate and station pair identification letters. For Geographic-to-LORAN conversion, no iterative computations are required; the program is a straightforward coordinate conversion based upon the techniques described by the Navy. For LORAN-to-Geographic conversion, the original Navy program required a dead-reckoned position, near the actual unknown fix, to begin computations. No iteration was performed to obtain the LORAN fix, but internal program errors occurred at execution time if the dead-reckoned fix were displaced from the actual fix by more than a few minutes of latitude or longitude. In order to enhance usefulness of the program for the terminal user, an iterative routine was added which allows a single dead-reckoned position to be entered from the master data file for each LORAN chain. The results compare exactly with the LORAN-C navigation chart, and provide adequate benchmark data for general aviation flight planning and data analysis.

Piecuch, L. M.↗

Approximate and exact numerical computation of supersonic flow over an airfoil

The computation of inviscid supersonic flow over a two-dimensional airfoil is considered. There are two main nonlinear approaches which lead to approximate solutions. Small-amplitude theory gives valid solutions provided the airfoil thickness is not too great and the Mach number is not too high. The second type of approximation, shock expansion theory, employs the fact that even for flows with strong shocks the effect of the down-running characteristics remains small. This leads to an analytic solution at the airfoil. It is pointed out that the approximate solution is accurate enough for many cases of interest. The numerical method furnishes a rapid correction to the solution in those cases where it is not. The characteristic-streamline coordinate system is useful both for the computation of the approximate solution and the corrections, and is also convenient for displaying and interpreting the results.

Lewis, T. S.↗

Polarization Considerations

As light passes through a optical system the reflections and refractions will in general change the polarization state of the light. If we assume that all of the materials in the thin film coatings and substrate are isotropic and homogeneous then calculating the amount of "instrumental" polarization is a relatively straight forward task. In the following sections we will present a of the steps required to perform a 'polarization ray trace' calculation for a single ray and monochromatic and hence polarized light. The thin film portion of the calculation is also shown. The reason for explicitly showing the thin film equations is that there are sign conventions imposed on the boundary value equations by the orientation and handedness of the various coordinate frames which are attached to the geometric rays. The attenuation of light through a optical system, is relatively simple, and requires at the very least a lens (average) reflectivity or transmissivity. Determining the polarization sensitivity of a optical system is still relatively straight forward requiring at least a knowledge of the behavior of the "s" and "p" components at each interface for the chief ray. Determining the thin film induced aberrations of a optical system are somewhat more demanding. Questions about the arithmetic sign of the phase factors and how this relates to the overall "OPD" of a ray are ubiquitous. Many rays are required to construct a wavefront. Thin film codes which modify the OPD's of rays are a requirement for this last mentioned computation. This requires a consistent scheme of coordinate frames and sign conventions and is probably the most demanding task of a polarization ray trace. Only the electric field will used in the discussion. This is not a restriction as the Stokes parameters are functions of the electric field. The following does not attempt to explain, but only to present all of the required concepts and formulas.

Waluschka, Eugene↗