Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “grid operators”

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 217 records · Page 12

Domain Decomposition By the Advancing-Partition Method for Parallel Unstructured Grid Generation

A new method of domain decomposition has been developed for generating unstructured grids in subdomains either sequentially or using multiple computers in parallel. Domain decomposition is a crucial and challenging step for parallel grid generation. Prior methods are generally based on auxiliary, complex, and computationally intensive operations for defining partition interfaces and usually produce grids of lower quality than those generated in single domains. The new technique, referred to as "Advancing Partition," is based on the Advancing-Front method, which partitions a domain as part of the volume mesh generation in a consistent and "natural" way. The benefits of this approach are: 1) the process of domain decomposition is highly automated, 2) partitioning of domain does not compromise the quality of the generated grids, and 3) the computational overhead for domain decomposition is minimal. The new method has been implemented in NASA's unstructured grid generation code VGRID.

Pirzadeh, Shahyar Z.↗

The BANDIT computer program for the reduction of matrix bandwidth for NASTRAN

A computer program, BANDIT, is described which was developed for use as a preprocessor to NASTRAN to automatically resequence the grid point numbers for reduced matrix bandwidth. The BANDIT program accepts a standard NASTRAN data deck as input, resequences the grid point numbers for reduced bandwidth, and generates a set of SEQGP cards for insertion into the NASTRAN bulk data deck. BANDIT uses the Cuthill-McKee strategy for resequencing grid points. The program is operational on the CDC 6000 series of computers, the IBM 360, and the Univac 1108.

Everstine, G. C.↗

Spatial land-use inventory, modeling, and projection/Denver metropolitan area, with inputs from existing maps, airphotos, and LANDSAT imagery

A landscape model was constructed with 34 land-use, physiographic, socioeconomic, and transportation maps. A simple Markov land-use trend model was constructed from observed rates of change and nonchange from photointerpreted 1963 and 1970 airphotos. Seven multivariate land-use projection models predicting 1970 spatial land-use changes achieved accuracies from 42 to 57 percent. A final modeling strategy was designed, which combines both Markov trend and multivariate spatial projection processes. Landsat-1 image preprocessing included geometric rectification/resampling, spectral-band, and band/insolation ratioing operations. A new, systematic grid-sampled point training-set approach proved to be useful when tested on the four orginal MSS bands, ten image bands and ratios, and all 48 image and map variables (less land use). Ten variable accuracy was raised over 15 percentage points from 38.4 to 53.9 percent, with the use of the 31 ancillary variables. A land-use classification map was produced with an optimal ten-channel subset of four image bands and six ancillary map variables. Point-by-point verification of 331,776 points against a 1972/1973 U.S. Geological Survey (UGSG) land-use map prepared with airphotos and the same classification scheme showed average first-, second-, and third-order accuracies of 76.3, 58.4, and 33.0 percent, respectively.

Tom, C.↗

Performance characteristics of ring-cusp thrusters with xenon propellant

The performance characteristics and operating envelope of several 30-cm ring-cusp ion thrusters with xenon propellant were investigated. Results indicate a strong performance dependence on the discharge chamber boundary magnetic fields and resultant distribution of electron currents. Significant improvements in discharge performance over J-series divergent-field thrusters were achieved for large throttling ranges, which translate into reduced cathode emission currents and reduced power dissipation which should be of significant benefit for operation at thruster power levels in excess of 10 kW. Mass spectrometry of the ion beam was documented for both the ring-cusp and J-series thrusters with xenon propellant for determination of overall thruster efficiency, and lifetime. Based on the lower centerline values of doubly charged ions in the ion beam and the lower operating discharge voltage, the screen grid erosion rate of the ring-cusp thruster is expected to be lower than the divergent-field J-series thruster by a factor of 2.

Patterson, M. J.↗

A linearized Euler analysis of unsteady flows in turbomachinery

A method for calculating unsteady flows in cascades is presented. The model, which is based on the linearized unsteady Euler equations, accounts for blade loading shock motion, wake motion, and blade geometry. The mean flow through the cascade is determined by solving the full nonlinear Euler equations. Assuming the unsteadiness in the flow is small, then the Euler equations are linearized about the mean flow to obtain a set of linear variable coefficient equations which describe the small amplitude, harmonic motion of the flow. These equations are discretized on a computational grid via a finite volume operator and solved directly subject to an appropriate set of linearized boundary conditions. The steady flow, which is calculated prior to the unsteady flow, is found via a Newton iteration procedure. An important feature of the analysis is the use of shock fitting to model steady and unsteady shocks. Use of the Euler equations with the unsteady Rankine-Hugoniot shock jump conditions correctly models the generation of steady and unsteady entropy and vorticity at shocks. In particular, the low frequency shock displacement is correctly predicted. Results of this method are presented for a variety of test cases. Predicted unsteady transonic flows in channels are compared to full nonlinear Euler solutions obtained using time-accurate, time-marching methods. The agreement between the two methods is excellent for small to moderate levels of flow unsteadiness. The method is also used to predict unsteady flows in cascades due to blade motion (flutter problem) and incoming disturbances (gust response problem).

Hall, Kenneth C.↗

Performance characteristics of ring-cusp thrusters with xenon propellant

The performance characteristics and operating envelope of several 30-cm ring-cusp ion thrusters with xenon propellant were investigated. Results indicate a strong performance dependence on the discharge chamber boundary magnetic fields and resultant distribution of electron currents. Significant improvements in discharge performance over J-series divergent-field thrusters were achieved for large throttling ranges, which translate into reduced cathode emission currents and reduced power dissipation which should be of significant benefit for operation at thruster power levels in excess of 10 kW. Mass spectrometer of the ion beam was documented for both the ring-cusp and J-series thrusters with xenon propellant for determination of overall thruster efficiency, and lifetime. Based on the lower centerline values of doubly charged ions in the ion beam and the lower operating discharge voltage, the screen grid erosion rate of the ring-cusp thruster is expected to be lower than the divergent-field J-series thruster by a factor of 2.

Patterson, M. J.↗

Calculation of unsteady flows in turbomachinery using the linearized Euler equations

A method for calculating unsteady flows in cascades is presented. The model, which is based on the linearized unsteady Euler equations, accounts for blade loading shock motion, wake motion, and blade geometry. The mean flow through the cascade is determined by solving the full nonlinear Euler equations. Assuming the unsteadiness in the flow is small, then the Euler equations are linearized about the mean flow to obtain a set of linear variable coefficient equations which describe the small amplitude, harmonic motion of the flow. These equations are discretized on a computational grid via a finite volume operator and solved directly subject to an appropriate set of linearized boundary conditions. The steady flow, which is calculated prior to the unsteady flow, is found via a Newton iteration procedure. An important feature of the analysis is the use of shock fitting to model steady and unsteady shocks. Use of the Euler equations with the unsteady Rankine-Hugoniot shock jump conditions correctly models the generation of steady and unsteady entropy and vorticity at shocks. In particular, the low frequency shock displacement is correctly predicted. Results of this method are presented for a variety of test cases. Predicted unsteady transonic flows in channels are compared to full nonlinear Euler solutions obtained using time-accurate, time-marching methods. The agreement between the two methods is excellent for small to moderate levels of flow unsteadiness. The method is also used to predict unsteady flows in cascades due to blade motion (flutter problem) and incoming disturbances (gust response problem).

Hall, Kenneth C.↗

The large discretization step method for time-dependent partial differential equations

A new method for the acceleration of linear and nonlinear time dependent calculations is presented. It is based on the Large Discretization Step (LDS) approximation, defined in this work, which employs an extended system of low accuracy schemes to approximate a high accuracy discrete approximation to a time dependent differential operator. Error bounds on such approximations are derived. These approximations are efficiently implemented in the LDS methods for linear and nonlinear hyperbolic equations, presented here. In these algorithms the high and low accuracy schemes are interpreted as the same discretization of a time dependent operator on fine and coarse grids, respectively. Thus, a system of correction terms and corresponding equations are derived and solved on the coarse grid to yield the fine grid accuracy. These terms are initialized by visiting the fine grid once in many coarse grid time steps. The resulting methods are very general, simple to implement and may be used to accelerate many existing time marching schemes.

Haras, Zigo↗

High-order Two-way Artificial Boundary Conditions for Nonlinear Wave Propagation with Backscattering

When solving linear scattering problems, one typically first solves for the impinging wave in the absence of obstacles. Then, by linear superposition, the original problem is reduced to one that involves only the scattered waves driven by the values of the impinging field at the surface of the obstacles. In addition, when the original domain is unbounded, special artificial boundary conditions (ABCs) that would guarantee the reflectionless propagation of waves have to be set at the outer boundary of the finite computational domain. The situation becomes conceptually different when the propagation equation is nonlinear. In this case the impinging and scattered waves can no longer be separated, and the problem has to be solved in its entirety. In particular, the boundary on which the incoming field values are prescribed, should transmit the given incoming waves in one direction and simultaneously be transparent to all the outgoing waves that travel in the opposite direction. We call this type of boundary conditions two-way ABCs. In the paper, we construct the two-way ABCs for the nonlinear Helmholtz equation that models the laser beam propagation in a medium with nonlinear index of refraction. In this case, the forward propagation is accompanied by backscattering, i.e., generation of waves in the direction opposite to that of the incoming signal. Our two-way ABCs generate no reflection of the backscattered waves and at the same time impose the correct values of the incoming wave. The ABCs are obtained for a fourth-order accurate discretization to the Helmholtz operator; the fourth-order grid convergence is corroborated experimentally by solving linear model problems. We also present solutions in the nonlinear case using the two-way ABC which, unlike the traditional Dirichlet boundary condition, allows for direct calculation of the magnitude of backscattering.

Fibich, Gadi↗

ICEG2D (v2.0) - An Integrated Software Package for Automated Prediction of Flow Fields for Single-Element Airfoils With Ice Accretion

An integrated geometry/grid/simulation software package, ICEG2D, is being developed to automate computational fluid dynamics (CFD) simulations for single- and multi-element airfoils with ice accretions. The current version, ICEG213 (v2.0), was designed to automatically perform four primary functions: (1) generate a grid-ready surface definition based on the geometrical characteristics of the iced airfoil surface, (2) generate high-quality structured and generalized grids starting from a defined surface definition, (3) generate the input and restart files needed to run the structured grid CFD solver NPARC or the generalized grid CFD solver HYBFL2D, and (4) using the flow solutions, generate solution-adaptive grids. ICEG2D (v2.0) can be operated in either a batch mode using a script file or in an interactive mode by entering directives from a command line within a Unix shell. This report summarizes activities completed in the first two years of a three-year research and development program to address automation issues related to CFD simulations for airfoils with ice accretions. As well as describing the technology employed in the software, this document serves as a users manual providing installation and operating instructions. An evaluation of the software is also presented.

Thompson David S.↗

Translation Optics for 30 cm Ion Engine Thrust Vector Control

Data were obtained from a 30 cm xenon ion thruster in which the accelerator grid was translated in the radial plane. The thruster was operated at three different throttle power levels, and the accelerator grid was incrementally translated in the X, Y, and azimuthal directions. Plume data was obtained downstream from the thruster using a Faraday probe mounted to a positioning system. Successive probe sweeps revealed variations in the plume direction. Thruster perveance, electron backstreaming limit, accelerator current, and plume deflection angle were taken at each power level, and for each accelerator grid position. Results showed that the thruster plume could easily be deflected up to six degrees without a prohibitive increase in accelerator impingement current. Results were similar in both X and Y direction.

Haag, Thomas↗

Study of Near-Surface Models in Large-Eddy Simulations of a Neutrally Stratified Atmospheric Boundary Layer

Large-eddy simulation (LES) is a widely used technique in armospheric modeling research. In LES, large, unsteady, three dimensional structures are resolved and small structures that are not resolved on the computational grid are modeled. A filtering operation is applied to distinguish between resolved and unresolved scales. We present two near-surface models that have found use in atmospheric modeling. We also suggest a simpler eddy viscosity model that adopts Prandtl's mixing length model (Prandtl 1925) in the vicinity of the surface and blends with the dynamic Smagotinsky model (Germano et al, 1991) away from the surface. We evaluate the performance of these surface models by simulating a neutraly stratified atmospheric boundary layer.

Senocak, I.↗

5,120 Superconducting Bolometers for the PIPER Balloon-Borne CMB Polarization Experiment

We are constructing the Primordial Inflation Polarization Explorer (PIPER) to measure the polarization of the cosmic microwave background (CMB) and search for the imprint of gravity waves produced during an inflationary epoch in the early universe. The signal is faint and lies behind confusing foregrounds, both astrophysical and cosmological, and so many detectors are required to complete the measurement in a limited time. We will use four of our matured 1,280 pixel, high-filling-factor backshort-under-grid bolometer arrays for efficient operation at the PIPER CMB wavelengths. All four arrays observe at a common wavelength set by passband filters in the optical path. PIPER will fly four times to observe at wavelengths of 1500, 1100, 850, and 500 microns in order to separate CMB from foreground emission. The arrays employ leg-isolated superconducting transition edge sensor bolometers operated at 145 mK; tuned resonant backshorts for efficient optical coupling; and a second-generation superconducting quantum interference device multiplexer readout. We describe the design, development, and performance of PIPER bolometer array technology to achieve background-limited sensitivity for a cryogenic balloon-borne telescope.

Benford, Dominic J.↗

5,120 Superconducting Bolometers for the PIPER Balloon-Borne CMB Polarization Experiment

We are constructing the Primordial Inflation Polarization Explorer (PIPER) to measure the polarization o[ the cosmic microwave background (CMB) and search for the imprint of gravity waves produced during an inflationary epoch in the early universe. The signal is faint and lies behind confusing foregrounds, both astrophysical and cosmological, and so many detectors are required to complete the measurement in a limited time. We will use four of our matured 1,280 pixel, high-filling-factor backshort-under-grid bolometer arrays for efficient operation at the PIPER CMB wavelengths. All four arrays observe at a common wavelength set by passband filters in the optical path. PIPER will fly four times to observe at wavelengths of 1500, 1100, 850, and 500 microns in order to separate CMB from foreground emission. The arrays employ leg-isolated superconducting transition edge sensor bolometers operated at 128mK; tuned resonant backshorts for efficient optical coupling; and a second-generation superconducting quantum interference device (SQUID) multiplexer readout. We describe the design, development, and performance of PIPER bo|ometer array technology to achieve background-limited sensitivity for a cryogenic balloon-borne telescope.

Benford, Dominic J.↗

Advanced Hall Electric Propulsion for Future In-space Transportation

The Hall thruster is an electric propulsion device used for multiple in-space applications including orbit raising, on-orbit maneuvers, and de-orbit functions. These in-space propulsion functions are currently performed by toxic hydrazine monopropellant or hydrazine derivative/nitrogen tetroxide bi-propellant thrusters. The Hall thruster operates nominally in the 1500 sec specific impulse regime. It provides greater thrust to power than conventional gridded ion engines, thus reducing trip times and operational life when compared to that technology in Earth orbit applications. The technology in the far term, by adding a second acceleration stage, has shown promise of providing over 4000s Isp, the regime of the gridded ion engine and necessary for deep space applications. The Hall thruster system consists of three parts, the thruster, the power processor, and the propellant system. The technology is operational and commercially available at the 1.5 kW power level and 5 kW application is underway. NASA is looking toward 10 kW and eventually 50 kW-class engines for ambitious space transportation applications. The former allows launch vehicle step-down for GEO missions and demanding planetary missions such as Europa Lander, while the latter allows quick all-electric propulsion LEO to GEO transfers and non-nuclear transportation human Mars missions.

Oleson, Steven R.↗

Benchmarking Memory Performance with the Data Cube Operator

Data movement across a computer memory hierarchy and across computational grids is known to be a limiting factor for applications processing large data sets. We use the Data Cube Operator on an Arithmetic Data Set, called ADC, to benchmark capabilities of computers and of computational grids to handle large distributed data sets. We present a prototype implementation of a parallel algorithm for computation of the operatol: The algorithm follows a known approach for computing views from the smallest parent. The ADC stresses all levels of grid memory and storage by producing some of 2d views of an Arithmetic Data Set of d-tuples described by a small number of integers. We control data intensity of the ADC by selecting the tuple parameters, the sizes of the views, and the number of realized views. Benchmarking results of memory performance of a number of computer architectures and of a small computational grid are presented.

Frumkin, Michael A.↗

Controlling Electron Backstreaming Phenomena Through the Use of a Transverse Magnetic Field

DEEP-SPACE mission propulsion requirements can be satisfied by the use of high specific impulse systems such as ion thrusters. For such missions. however. the ion thruster will be required to provide thrust for long periods of time. To meet the long operation time and high-propellant throughput requirements, thruster lifetime must be increased. In general, potential ion thruster failure mechanisms associated with long-duration thrusting can be grouped into four areas: (1) ion optics failure; (2) discharge cathode failure; (3) neutralizer failure; and (4) electron backstreaming caused by accelerator grid aperture enlargement brought on by accelerator grid erosion. The work presented here focuses on electron backstreaming. which occurs when the potential at the center of an accelerator grid aperture is insufficient to prevent the backflow of electrons into the ion thruster. The likelihood of this occurring depends on ion source operation time. plasma density, and grid voltages, as accelerator grid apertures enlarge as a result of erosion. Electrons that enter the gap between the high-voltage screen and accelerator grids are accelerated to the energies approximately equal to the beam voltage. This energetic electron beam (typically higher than 1 kV) can damage not only the ion source discharge cathode assembly. but also any of the discharge surfaces upstream of the ion acceleration optics that the electrons happen to impact. Indeed. past backstreaming studies have shown that near the backstreaming limit, which corresponds to the absolute value of the accelerator grid voltage below which electrons can backflow into the thruster, there is a rather sharp rise in temperature at structures such as the cathode keeper electrode. In this respect operation at accelerator grid voltages near the backstreaming limit is avoided. Generally speaking, electron backstreaming is prevented by operating the accelerator grid at a sufficiently negative voltage to ensure a sufficiently negative aperture center potential. This approach can provide the necessary margin assuming an expected aperture enlargement. Operation at very negative accelerator grid voltages, however, enhances ion charge-exchange and direct impingement erosion of the accelerator grid. The focus of the work presented here is the mitigation of electron backstreaming by the use of a magnetic field. The presence of a magnetic field oriented perpendicular to the thruster axis can significantly decrease the magnitude of the backflowing electron current by significantly reducing the electron diffusion coefficient. Negative ion sources utilize this principle to reduce the fraction of electrons in the negative ion beam. The focus of these efforts has been on the attenuation of electron current diffusing from the discharge plasma into the negative ion extraction optics by placing the transverse magnetic field upstream of the extraction electrodes. In contrast. in the case of positive ion sources such as ion thrusters, the approach taken in the work presented here is to apply the transverse field downstream of the ion extraction system so as to prevent electrons from flowing back into the source. It was found in the work presented here that the magnetic field also reduces the absolute value of the electron backstreaming limit voltage. In this respect. the applied transverse magnetic field provides two mechanisms for electron backstreaming mitigation: (1) electron current attenuation and (2) backstreaming limit voltage shift. Such a shift to less negative voltages can lead to reduced accelerator grid erosion rates.

Foster, John E.↗

A Dark Target research aerosol algorithm for MODIS observations over eastern China: increasing coverage while maintaining accuracy at high aerosol loading

Satellite aerosol products such as the Dark Target (DT) produced from the MODerate resolution Imaging Spectroradiometer (MODIS) are useful for monitoring the progress of air pollution. Unfortunately, the DT often fails to retrieve during the heaviest aerosol events as well as the more moderate events in winter. Some of the literature at-tributes this lack of retrieval to the cloud mask. However, we found this lack of retrieval is mainly traced to thresholds used for masking of inland water and snow. Modifications to these two masks greatly increase 50 % of the retrievals of aerosol optical depth at 0.55 μm (AOD) greater than 1.0. The “extra”-high-AOD retrievals tend to be biased when com-pared with a ground-based sun photometer (AErosol RObotic NETwork, AERONET). Reducing bias in new retrievals re-quires two additional steps. One is an update to the assumed aerosol optical properties (aerosol model); the haze in this region is both less absorbing and lower in altitude than what is assumed in the global algorithm. The second is account-ing for the scale height of the aerosol, specifically that the heavy-aerosol events in the region are much closer to the surface than what is assumed by the global DT algorithm. The resulting combination of modified masking thresholds, new aerosol model, and lower aerosol layer scale height was applied to 3 months of MODIS observations (January–March2013) over eastern China. After these two additional steps are implemented, the significant increase in new retrievals introduces no overall bias at a high-AOD regime but does degrade other overall validation statistics. We also find that the research algorithm is able to identify additional pollution events that AERONET instruments may not due to different spatial sampling. Mean AOD retrieved from the re-search algorithm increases from 0.11 to 0.18 compared to values calculated from the operational DT algorithm during January to March of 2013 over the study area. But near Beijing, where the severe pollution occurs, the new algorithm increases AOD by as much as 3.0 for each 0.5°grid box over the previous operational-algorithm values.

Dark Target↗