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 253 records · Page 14

Impact of the Assimilation of Hyperspectral Infrared Retrieved Profiles on Advanced Weather and Research Model Simulations of a Non-Convective Wind Event

Non-convective wind events commonly occur with passing extratropical cyclones and have significant societal and economic impacts. Since non-convective winds often occur in the absence of specific phenomena such as a thunderstorm, tornado, or hurricane, the public are less likely to heed high wind warnings and continue daily activities. Thus non-convective wind events result in as many fatalities as straight line thunderstorm winds. One physical explanation for non-convective winds includes tropopause folds. Improved model representation of stratospheric air and associated non-convective wind events could improve non-convective wind forecasts and associated warnings. In recent years, satellite data assimilation has improved skill in forecasting extratropical cyclones; however errors still remain in forecasting the position and strength of extratropical cyclones as well as the tropopause folding process. The goal of this study is to determine the impact of assimilating satellite temperature and moisture retrieved profiles from hyperspectral infrared (IR) sounders (i.e. Atmospheric Infrared Sounder (AIRS), Cross-track Infrared and Microwave Sounding Suite (CrIMSS), and Infrared Atmospheric Sounding Interferometer (IASI)) on the model representation of the tropopause fold and an associated high wind event that impacted the Northeast United States on 09 February 2013. Model simulations using the Advanced Research Weather Research and Forecasting Model (ARW) were conducted on a 12-km grid with cycled data assimilation mimicking the operational North American Model (NAM). The results from the satellite assimilation run are compared to a control experiment (without hyperspectral IR retrievals), 32-km North American Regional Reanalysis (NARR) interpolated to a 12-km grid, and 13-km Rapid Refresh analyses.

Berndt, E. B.↗

Global Discrete Artificial Boundary Conditions for Time-Dependent Wave Propagation

We construct global artificial boundary conditions (ABCs) for the numerical simulation of wave processes on unbounded domains using a special non-deteriorating algorithm that has been developed previously for the long-term computation of wave-radiation solutions. The ABCs are obtained directly for the discrete formulation of the problem; in so doing, neither a rational approximation of 'non-reflecting kernels,' nor discretization of the continuous boundary conditions is required. The extent of temporal nonlocality of the new ABCs appears fixed and limited; in addition, the ABCs can handle artificial boundaries of irregular shape on regular grids with no fitting/adaptation needed and no accuracy loss induced. The non-deteriorating algorithm, which is the core of the new ABCs is inherently three-dimensional, it guarantees temporally uniform grid convergence of the solution driven by a continuously operating source on arbitrarily long time intervals, and provides unimprovable linear computational complexity with respect to the grid dimension. The algorithm is based on the presence of lacunae, i.e., aft fronts of the waves, in wave-type solutions in odd-dimension spaces, It can, in fact, be built as a modification on top of any consistent and stable finite-difference scheme, making its grid convergence uniform in time and at the same time keeping the rate of convergence the same as that of the non-modified scheme. In the paper, we delineate the construction of the global lacunae-based ABCs in the framework of a discretized wave equation. The ABCs are obtained for the most general formulation of the problem that involves radiation of waves by moving sources (e.g., radiation of acoustic waves by a maneuvering aircraft). We also present systematic numerical results that corroborate the theoretical design properties of the ABCs' algorithm.

Ryabenkii, V. S.↗

Design and Expected Performance of GISMO-2, a Two Color Millimeter Camera for the IRAM 30 m Telescope

We present the main design features for the GISMO-2 bolometer camera, which we build for background-limited operation at the IRAM 30 m telescope on Pico Veleta, Spain. GISMO-2 will operate simultaneously in the 1 and 2 mm atmospherical windows. The 1 mm channel uses a 32 × 40 TES-based backshort under grid (BUG) bolometer array, the 2 mm channel operates with a 16 × 16 BUG array. The camera utilizes almost the entire full field of view provided by the telescope. The optical design of GISMO-2 was strongly influenced by our experience with the GISMO 2mm bolometer camera, which is successfully operating at the 30 m telescope. GISMO is accessible to the astronomical community through the regularIRAMcall for proposals.

Two-Color mm camera↗

Drive Train Normal Modes Analysis for the ERDA/NASA 100-Kilowatt Wind Turbine Generator

Natural frequencies, as a function of power were determined using a finite element model. Operating conditions investigated were operation with a resistive electrical load and operation synchronized to an electrical utility grid. The influence of certain drive train components on frequencies and mode shapes is shown. An approximate method for obtaining drive train natural frequencies is presented.

Sullivan, T. L.↗

Phase 1 of the First Solar Small Power System Experiment (experimental System No. 1). Volume 1: Technical Studies for Solar Point-focusing, Distributed Collector System, with Energy Conversion at the Collector, Category C

The technical and economic feasibility of a solar electric power plant for a small community is evaluated and specific system designs for development and demonstration are selected. All systems investigated are defined as point focusing, distributed receiver concepts, with energy conversion at the collector. The preferred system is comprised of multiple parabolic dish concentrators employing Stirling cycle engines for power conversion. The engine, AC generator, cavity receiver, and integral sodium pool boiler/heat transport system are combined in a single package and mounted at the focus of each concentrator. The output of each concentrator is collected by a conventional electrical distribution system which permits grid-connected or stand-alone operation, depending on the storage system selected.

Clark, T. B.↗

Low-level circulations induced by an upper-level jet in a fine mesh numerical model

Attention is given to varying the initial conditions in the Penn State three-dimensional numerical model in the development of the low-level jet over the Ohio Valley during the 10 May 1973 severe weather event. The results of two forecasts are summarized: (1) the control forecast which interpolated the National Meteorological Center's operational analysis to the mesoscale model grid; and (2) modifying the control forecast by inserting artificial soundings and surface observations. Enhancements are made on the basis of subsynoptic scale analyses by Uccellini and surface data. The results indicate that the low-level jet was in part induced by pressure changes associated with deceleration and horizontal divergence of air in the exit region of the upper-level jet. The enhancement of the initial moisture field using surface reports of clouds and precipitation has a substantial effect on precipitation.

Burkhart, R. P.↗

Applications of the AVE-Sesame data sets to mesoscale studies

Data collected by the lightning data concentrator are available for research. The Mark 3 McIDAS capability provides greater flexibility for the Marshall user community and serves as a model of future UW McIDAS to remote computer links. Techniques were investigated for the display of dynamic 3-D data sets. To date the most promising display technology is a polarized two CRT perspective display which allows both dynamic 3-D images and graphics presentations with full color capability. Algorithms were for the preparation and display of conventional and satellite based weather data in 3-D. These include gridding, contouring, and streamlining processors which operate on both real time and case study data bases. An upper air trajectory model was implemented which creates a display of air parcel trajectories in perspective 3-D. A subsystem for the generation of 3-D solid surface display with shading and hidden surface display with shading and hidden surface removal was tested and its products are currently being evaluated. Motion parallax introduced by moving the point of observation during display is an important depth cue, which, when added to the perspective parallax creates a very realistic appearing display.

Suomi, V. E.↗

Testing of a failure accommodation system on a highly flexible grid

This paper presents a scheme for on-line failure detection in systems with appreciable structural dynamics. The design is suboptimal because of extensive computational requirements of optimal schemes. To accomplish failure detection, the innovations sequence of a finite order Kalman filter is examined. Because of the heavy dependence of the system on the zero-mean character of the innovations sequence of the filter much attention has been given to the design and evaluation of the Kalman filters used. The filter designs are based on modal models of the structural dynamics. Two modal models were considered, one based on an analytic finite element model and the other based on empirically derived frequency and damping. Experiments using a grid structure are presented which illustrate operation and performance of the filter designs based on these models. The general character of the results presented is that appreciable errors exist in the filter design based on the finite element model. Substantial improvement results if the design model is modified to include empirically derived frequency and damping.

Montgomery, R. C.↗

Advanced adaptive computational methods for Navier-Stokes simulations in rotorcraft aerodynamics

A phase 2 research and development effort was conducted in area transonic, compressible, inviscid flows with an ultimate goal of numerically modeling complex flows inherent in advanced helicopter blade designs. The algorithms and methodologies therefore are classified as adaptive methods, which are error estimation techniques for approximating the local numerical error, and automatically refine or unrefine the mesh so as to deliver a given level of accuracy. The result is a scheme which attempts to produce the best possible results with the least number of grid points, degrees of freedom, and operations. These types of schemes automatically locate and resolve shocks, shear layers, and other flow details to an accuracy level specified by the user of the code. The phase 1 work involved a feasibility study of h-adaptive methods for steady viscous flows, with emphasis on accurate simulation of vortex initiation, migration, and interaction. Phase 2 effort focused on extending these algorithms and methodologies to a three-dimensional topology.

Stowers, S. T.↗

Multi-grid domain decomposition approach for solution of Navier-Stokes equations in primitive variable form

The new multigrid (or adaptive) pseudospectral element method was carried out for the solution of incompressible flow in terms of primitive variable formulation. The desired features of the proposed method include the following: (1) the ability to treat complex geometry; (2) high resolution adapted in the interesting areas; (3) requires minimal working space; and (4) effective in a multiprocessing environment. The approach for flow problems, complex geometry or not, is to first divide the computational domain into a number of fine-grid and coarse-grid subdomains with the inter-overlapping area. Next, it is necessary to implement the Schwarz alternating procedure (SAP) to exchange the data among subdomains, where the coarse-grid correction is used to remove the high frequency error that occurs when the data interpolation from the fine-grid subdomain to the coarse-grid subdomain is conducted. The strategy behind the coarse-grid correction is to adopt the operator of the divergence of the velocity field, which intrinsically links the pressure equation, into this process. The solution of each subdomain can be efficiently solved by the direct (or iterative) eigenfunction expansion technique with the least storage requirement, i.e. O(N(exp 3)) in 3-D and O(N(exp 2)) in 2-D. Numerical results of both driven cavity and jet flow will be presented in the paper to account for the versatility of the proposed method.

Ku, Hwar-Ching↗

Development of an Ion Thruster and Power Processor for New Millennium's Deep Space 1 Mission

The NASA Solar Electric Propulsion Technology Applications Readiness Program (NSTAR) will provide a single-string primary propulsion system to NASA's New Millennium Deep Space 1 Mission which will perform comet and asteroid flybys in the years 1999 and 2000. The propulsion system includes a 30-cm diameter ion thruster, a xenon feed system, a power processing unit, and a digital control and interface unit. A total of four engineering model ion thrusters, three breadboard power processors, and a controller have been built, integrated, and tested. An extensive set of development tests has been completed along with thruster design verification tests of 2000 h and 1000 h. An 8000 h Life Demonstration Test is ongoing and has successfully demonstrated more than 6000 h of operation. In situ measurements of accelerator grid wear are consistent with grid lifetimes well in excess of the 12,000 h qualification test requirement. Flight hardware is now being assembled in preparation for integration, functional, and acceptance tests.

Sovey, James S.↗

Flow Applications of the Least Squares Finite Element Method

The main thrust of the effort has been towards the development, analysis and implementation of the least-squares finite element method (LSFEM) for fluid dynamics and electromagnetics applications. In the past year, there were four major accomplishments: 1) special treatments in computational fluid dynamics and computational electromagnetics, such as upwinding, numerical dissipation, staggered grid, non-equal order elements, operator splitting and preconditioning, edge elements, and vector potential are unnecessary; 2) the analysis of the LSFEM for most partial differential equations can be based on the bounded inverse theorem; 3) the finite difference and finite volume algorithms solve only two Maxwell equations and ignore the divergence equations; and 4) the first numerical simulation of three-dimensional Marangoni-Benard convection was performed using the LSFEM.

Jiang, Bo-Nan↗

Distributed Relaxation for Conservative Discretizations

A multigrid method is defined as having textbook multigrid efficiency (TME) if the solutions to the governing system of equations are attained in a computational work that is a small (less than 10) multiple of the operation count in one target-grid residual evaluation. The way to achieve this efficiency is the distributed relaxation approach. TME solvers employing distributed relaxation have already been demonstrated for nonconservative formulations of high-Reynolds-number viscous incompressible and subsonic compressible flow regimes. The purpose of this paper is to provide foundations for applications of distributed relaxation to conservative discretizations. A direct correspondence between the primitive variable interpolations for calculating fluxes in conservative finite-volume discretizations and stencils of the discretized derivatives in the nonconservative formulation has been established. Based on this correspondence, one can arrive at a conservative discretization which is very efficiently solved with a nonconservative relaxation scheme and this is demonstrated for conservative discretization of the quasi one-dimensional Euler equations. Formulations for both staggered and collocated grid arrangements are considered and extensions of the general procedure to multiple dimensions are discussed.

Diskin, Boris↗

Miniature Ion-Mobility Spectrometer

The figure depicts a proposed miniature ion-mobility spectrometer that would be fabricated by micromachining. Unlike prior ion-mobility spectrometers, the proposed instrument would not be based on a time-of-flight principle and, consequently, would not have some of the disadvantageous characteristics of prior time-of-flight ion-mobility spectrometers. For example, one of these characteristics is the need for a bulky carrier-gas-feeding subsystem that includes a shutter gate to provide short pulses of gas in order to generate short pulses of ions. For another example, there is need for a complex device to generate pulses of ions from the pulses of gas and the device is capable of ionizing only a fraction of the incoming gas molecules; these characteristics preclude miniaturization. In contrast, the proposed instrument would not require a carrier-gas-feeding subsystem and would include a simple, highly compact device that would ionize all the molecules passing through it. The ionization device in the proposed instrument would be a 0.1-micron-thick dielectric membrane with metal electrodes on both sides. Small conical holes would be micromachined through the membrane and electrodes. An electric potential of the order of a volt applied between the membrane electrodes would give rise to an electric field of the order of several megavolts per meter in the submicron gap between the electrodes. An electric field of this magnitude would be sufficient to ionize all the molecules that enter the holes. Ionization (but not avalanche arcing) would occur because the distance between the ionizing electrodes would be less than the mean free path of gas molecules at the operating pressure of instrument. An accelerating grid would be located inside the instrument, downstream from the ionizing membrane. The electric potential applied to this grid would be negative relative to the potential on the inside electrode of the ionizing membrane and would be of a magnitude sufficient to generate a moderate electric field. Positive ions leaving the membrane holes would be accelerated in this electric field. The resulting flux of ions away from the ionization membrane would create a partial vacuum that would draw more of the gas medium through the membrane. The figure depicts a filter electrode and detector electrodes located along the sides of a drift tube downstream from the accelerator electrode. These electrodes would apply a transverse AC electric field superimposed on a ramped DC electric field. The AC field would effect differential transverse dispersal of ions. At a given instant of time, the trajectories of most of the ions would be bent toward the electrodes, causing most of the ions to collide with the electrodes and thereby become neutralized. The DC field would partly counteract the dispersive effect of the AC field, straightening the trajectories of a selected species of ions; the selection would vary with the magnitude of the applied DC field. The straightening of the trajectories of the selected ions would enable them to pass into the region between the detector electrodes. Depending on the polarity of the voltage applied to the detector electrodes, the electric field between the detector electrodes would draw the selected ions to one of these electrodes. Hence, the current collected by one of the detector electrodes would be a measure of the abundance of ions of the selected species. The ramping of the filter- electrode DC voltage would sweep the selection of ions through the spectrum of ionic species.

Hartley, Frank T.↗

Calculating Soil Wetness, Evapotranspiration and Carbon Cycle Processes Over Large Grid Areas Using a New Scaling Technique

Soil wetness typically shows great spatial variability over the length scales of general circulation model (GCM) grid areas (approx 100 km ), and the functions relating evapotranspiration and photosynthetic rate to local-scale (approx 1 m) soil wetness are highly non-linear. Soil respiration is also highly dependent on very small-scale variations in soil wetness. We therefore expect significant inaccuracies whenever we insert a single grid area-average soil wetness value into a function to calculate any of these rates for the grid area. For the particular case of evapotranspiration., this method - use of a grid-averaged soil wetness value - can also provoke severe oscillations in the evapotranspiration rate and soil wetness under some conditions. A method is presented whereby the probability distribution timction(pdf) for soil wetness within a grid area is represented by binning. and numerical integration of the binned pdf is performed to provide a spatially-integrated wetness stress term for the whole grid area, which then permits calculation of grid area fluxes in a single operation. The method is very accurate when 10 or more bins are used, can deal realistically with spatially variable precipitation, conserves moisture exactly and allows for precise modification of the soil wetness pdf after every time step. The method could also be applied to other ecological problems where small-scale processes must be area-integrated, or upscaled, to estimate fluxes over large areas, for example in treatments of the terrestrial carbon budget or trace gas generation.

Sellers, Piers↗

Entropy Stable Staggered Grid Spectral Collocation for the Burgers' and Compressible Navier-Stokes Equations

Staggered grid, entropy stable discontinuous spectral collocation operators of any order are developed for Burgers' and the compressible Navier-Stokes equations on unstructured hexahedral elements. This generalization of previous entropy stable spectral collocation work [1, 2], extends the applicable set of points from tensor product, Legendre-Gauss-Lobatto (LGL) to a combination of tensor product Legendre-Gauss (LG) and LGL points. The new semi-discrete operators discretely conserve mass, momentum, energy and satisfy a mathematical entropy inequality for both Burgers' and the compressible Navier-Stokes equations in three spatial dimensions. They are valid for smooth as well as discontinuous flows. The staggered LG and conventional LGL point formulations are compared on several challenging test problems. The staggered LG operators are significantly more accurate, although more costly to implement. The LG and LGL operators exhibit similar robustness, as is demonstrated using test problems known to be problematic for operators that lack a nonlinearly stability proof for the compressible Navier-Stokes equations (e.g., discontinuous Galerkin, spectral difference, or flux reconstruction operators).

Carpenter, Mark H.↗

Determination of Best Low-Frequency Microwave Antenna Approach for Future High Resolution Measurements from Space

Microwave remote sensing measurements at L-band (~1.2-1.6 GHz) of geophysical parameters such as soil moisture will need to be at higher spatial resolution than current systems (SMOS/ SMAP/ Aquarius) in order to meet the requirements of land surface, ocean, and numerical weather prediction models in the near future, which will operate at ~9-15 km global grids and 1-3 km regional grids in the next few years. In order to make progress toward these needed spatial resolutions, advancements in technology are necessary which would lead to improved effective (i.e. equivalent) antenna size. An architecture trade study was conducted to quantitatively define the value and limits of different microwave technology paths, and to select the most appropriate path to achieve the high spatial resolution required by science in the future without sacrificing performance, accuracy, and global coverage.

soil moisture↗

The GISMO-2 Bolometer Camera

We present the concept for the GISMO-2 bolometer camera) which we build for background-limited operation at the IRAM 30 m telescope on Pico Veleta, Spain. GISM0-2 will operate Simultaneously in the 1 mm and 2 mm atmospherical windows. The 1 mm channel uses a 32 x 40 TES-based Backshort Under Grid (BUG) bolometer array, the 2 mm channel operates with a 16 x 16 BUG array. The camera utilizes almost the entire full field of view provided by the telescope. The optical design of GISM0-2 was strongly influenced by our experience with the GISMO 2 mm bolometer camera which is successfully operating at the 30m telescope. GISMO is accessible to the astronomical community through the regular IRAM call for proposals.

Staguhn, Johannes G.↗