Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Parallel projection algorithm”

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.

113 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

Performance of an Optimized Eta Model Code on the Cray T3E and a Network of PCs

In the year 2001, NASA will launch the satellite TRIANA that will be the first Earth observing mission to provide a continuous, full disk view of the sunlit Earth. As a part of the HPCC Program at NASA GSFC, we have started a project whose objectives are to develop and implement a 3D cloud data assimilation system, by combining TRIANA measurements with model simulation, and to produce accurate statistics of global cloud coverage as an important element of the Earth's climate. For simulation of the atmosphere within this project we are using the NCEP/NOAA operational Eta model. In order to compare TRIANA and the Eta model data on approximately the same grid without significant downscaling, the Eta model will be integrated at a resolution of about 15 km. The integration domain (from -70 to +70 deg in latitude and 150 deg in longitude) will cover most of the sunlit Earth disc and will continuously rotate around the globe following TRIANA. The cloud data assimilation is supposed to run and produce 3D clouds on a near real-time basis. Such a numerical setup and integration design is very ambitious and computationally demanding. Thus, though the Eta model code has been very carefully developed and its computational efficiency has been systematically polished during the years of operational implementation at NCEP, the current MPI version may still have problems with memory and efficiency for the TRIANA simulations. Within this work, we optimize a parallel version of the Eta model code on a Cray T3E and a network of PCs (theHIVE) in order to improve its overall efficiency. Our optimization procedure consists of introducing dynamically allocated arrays to reduce the size of static memory, and optimizing on a single processor by splitting loops to limit the number of streams. All the presented results are derived using an integration domain centered at the equator, with a size of 60 x 60 deg, and with horizontal resolutions of 1/2 and 1/3 deg, respectively. In accompanying charts we report the elapsed time, the speedup and the Mflops as a function of the number of processors for the non-optimized version of the code on the T3E and theHIVE. The large amount of communication required for model integration explains its poor performance on theHIVE. Our initial implementation of the dynamic memory allocation has contributed to about 12% reduction of memory but has introduced a 3% overhead in computing time. This overhead was removed by performing loop splitting in some of the high demanding subroutines. When the Eta code is fully optimized in order to meet the memory requirement for TRIANA simulations, a non-negligeable overhead may appear that may seriously affect the efficiency of the code. To alleviate this problem, we are considering implementation of a new algorithm for the horizontal advection that is computationally less expensive, and also a new approach for marching in time.

Kouatchou, Jules↗
NASA NTRS2010

The OMPS Limb Profiler Instrument: Two-Dimensional Retrieval Algorithm

The upcoming Ozone Mapper and Profiler Suite (OMPS), which will be launched on the NPOESS Preparatory Project (NPP) platform in early 2011, will continue monitoring the global distribution of the Earth's middle atmosphere ozone and aerosol. OMPS is composed of three instruments, namely the Total Column Mapper (heritage: TOMS, OMI), the Nadir Profiler (heritage: SBUV) and the Limb Profiler (heritage: SOLSE/LORE, OSIRIS, SCIAMACHY, SAGE III). The ultimate goal of the mission is to better understand and quantify the rate of stratospheric ozone recovery. The focus of the paper will be on the Limb Profiler (LP) instrument. The LP instrument will measure the Earth's limb radiance (which is due to the scattering of solar photons by air molecules, aerosol and Earth surface) in the ultra-violet (UV), visible and near infrared, from 285 to 1000 nm. The LP simultaneously images the whole vertical extent of the Earth's limb through three vertical slits, each covering a vertical tangent height range of 100 km and each horizontally spaced by 250 km in the cross-track direction. Measurements are made every 19 seconds along the orbit track, which corresponds to a distance of about 150km. Several data analysis tools are presently being constructed and tested to retrieve ozone and aerosol vertical distribution from limb radiance measurements. The primary NASA algorithm is based on earlier algorithms developed for the SOLSE/LORE and SAGE III limb scatter missions. All the existing retrieval algorithms rely on a spherical symmetry assumption for the atmosphere structure. While this assumption is reasonable in most of the stratosphere, it is no longer valid in regions of prime scientific interest, such as polar vortex and UTLS regions. The paper will describe a two-dimensional retrieval algorithm whereby the ozone distribution is simultaneously retrieved vertically and horizontally for a whole orbit. The retrieval code relies on (1) a forward 2D Radiative Transfer code (to model limb radiances within a non-uniform atmosphere and evaluate 2D analytical partial derivatives) and (2) an optimal estimator inversion routine. The algorithm uses the typically sparse nature of the kernel matrices as well as fast matrix inversion techniques to allow for fast inversion of limb data with efficient memory management (as was done for MIPAS data processing). While the method has so far only been developed in the context of Single Scatter, the paper will show how the CPU intensive Multiple Scatter modeling can be implemented using parallel CPU processing. Initial results will be presented in terms of retrieved ozone profiles and code performance.

Rault, Didier F.↗
NASA NTRS2000

High Order Discontinuous Gelerkin Methods for Convection Dominated Problems with Application to Aeroacoustics

This project is about the investigation of the development of the discontinuous Galerkin finite element methods, for general geometry and triangulations, for solving convection dominated problems, with applications to aeroacoustics. On the analysis side, we have studied the efficient and stable discontinuous Galerkin framework for small second derivative terms, for example in Navier-Stokes equations, and also for related equations such as the Hamilton-Jacobi equations. This is a truly local discontinuous formulation where derivatives are considered as new variables. On the applied side, we have implemented and tested the efficiency of different approaches numerically. Related issues in high order ENO and WENO finite difference methods and spectral methods have also been investigated. Jointly with Hu, we have presented a discontinuous Galerkin finite element method for solving the nonlinear Hamilton-Jacobi equations. This method is based on the RungeKutta discontinuous Galerkin finite element method for solving conservation laws. The method has the flexibility of treating complicated geometry by using arbitrary triangulation, can achieve high order accuracy with a local, compact stencil, and are suited for efficient parallel implementation. One and two dimensional numerical examples are given to illustrate the capability of the method. Jointly with Hu, we have constructed third and fourth order WENO schemes on two dimensional unstructured meshes (triangles) in the finite volume formulation. The third order schemes are based on a combination of linear polynomials with nonlinear weights, and the fourth order schemes are based on combination of quadratic polynomials with nonlinear weights. We have addressed several difficult issues associated with high order WENO schemes on unstructured mesh, including the choice of linear and nonlinear weights, what to do with negative weights, etc. Numerical examples are shown to demonstrate the accuracies and robustness of the methods for shock calculations. Jointly with P. Montarnal, we have used a recently developed energy relaxation theory by Coquel and Perthame and high order weighted essentially non-oscillatory (WENO) schemes to simulate the Euler equations of real gas. The main idea is an energy decomposition under the form epsilon = epsilon(sub 1) + epsilon(sub 2), where epsilon(sub 1) is associated with a simpler pressure law (gamma)-law in this paper) and the nonlinear deviation epsilon(sub 2) is convected with the flow. A relaxation process is performed for each time step to ensure that the original pressure law is satisfied. The necessary characteristic decomposition for the high order WENO schemes is performed on the characteristic fields based on the epsilon(sub l) gamma-law. The algorithm only calls for the original pressure law once per grid point per time step, without the need to compute its derivatives or any Riemann solvers. Both one and two dimensional numerical examples are shown to illustrate the effectiveness of this approach.

Shu, Chi-Wang↗
NASA NTRS1999

Experimental Investigation of Pool Boiling Heat Transfer Enhancement in Microgravity in the Presence of Electric Fields

In boiling high heat fluxes are possible driven by relatively small temperature differences, which make its use increasingly attractive in aerospace applications. The objective of the research is to develop ways to overcome specific problems associated with boiling in the low gravity environment by substituting the buoyancy force with the electric force to enhance bubble removal from the heated surface. Previous studies indicate that in terrestrial applications nucleate boiling heat transfer can be increased by a factor of 50, as compared to values obtained for the same system without electric fields. The goal of our research is to experimentally explore the mechanisms responsible for EHD heat transfer enhancement in boiling in low gravity conditions, by visualizing the temperature distributions in the vicinity of the heated surface and around the bubble during boiling using real-time holographic interferometry (HI) combined with high-speed cinematography. In the first phase of the project the influence of the electric field on a single bubble is investigated. Pool boiling is simulated by injecting a single bubble through a nozzle into the subcooled liquid or into the thermal boundary layer developed along the flat heater surface. Since the exact location of bubble formation is known, the optical equipment can be aligned and focused accurately, which is an essential requirement for precision measurements of bubble shape, size and deformation, as well as the visualization of temperature fields by HI. The size of the bubble and the frequency of bubble departure can be controlled by suitable selection of nozzle diameter and mass flow rate of vapor. In this approach effects due to the presence of the electric field can be separated from effects caused by the temperature gradients in the thermal boundary layer. The influence of the thermal boundary layer can be investigated after activating the heater at a later stage of the research. For the visualization experiments a test cell was developed. All four vertical walls of the test cell are transparent, and they allow transillumination with laser light for visualization experiments by HI. The bottom electrode is a copper cylinder, which is electrically grounded. The copper block is heated with a resistive heater and it is equipped with 6 thermocouples that provide reference temperatures for the measurements with HI. The top electrode is a mesh electrode. Bubbles are injected with a syringe into the test cell through the bottom electrode. The working fluids presently used in the interferometric visualization experiments, water and PF 5052, satisfy requirements regarding thermophysical, optical and electrical properties. A 30kV power supply equipped with a voltmeter allows to apply the electric field to the electrodes during the experiments. The magnitude of the applied voltage can be adjusted either manually or through the LabVIEW data acquisition and control system connected to a PC. Temperatures of the heated block are recorded using type-T thermocouples, whose output is read by a data acquisition system. Images of the bubbles are recorded with 35mm photographic and 16mm high-speed cameras, scanned and analyzed using various software packages. Visualized temperature fields HI allows the visualization of temperature fields in the vicinity of bubbles during boiling in the form of fringes. Typical visualized temperature distributions around the air bubbles injected into the thermal boundary layer in PF5052 are shown. The temperature of the heated surface is 35 C. The temperature difference for a pair of fringes is approximately 0.05 C. The heat flux applied to the bottom surface is moderate, and the fringe patterns are regular. In the image a bubble penetrating the thermal boundary layer is visible. Because of the axial symmetry of the problem, simplified reconstruction techniques can be applied to recover the temperature field. The thermal plume developing above the heated surface for more intensive heating is shown. The temperature distribution in the liquid is clearly 3D, and tomographic techniques have to be applied to recover the temperature distribution in such a physical situation. A sequence of interferometric images showing the temperature distribution around the rising bubble, recorded with a high-speed camera is shown. Again, the temperature distribution is 3D, and a more complex approach to the evaluation, the tomographic reconstruction has to be taken. Measurement of the temperature distribution from the fringe pattern temperature distributions that yield important information regarding heat transfer are determined. Two algorithms that allow the quantitative evaluation of interferometric fringe patterns and the reconstruction of temperature fields during boiling have been developed at the Heat Transfer Laboratory of the Johns Hopkins University. In the first algorithm the bubble is assumed to be axially symmetrical, which significantly reduces the computational effort for quantifying the temperature distribution around the bubble. For this purpose the thermal boundary layer around the bubble is divided into equidistant concentric shells, and the refractive index is assumed to be constant in each of the shells. Since large temperature gradients are expected in the vicinity of the bubble during boiling, the deflection of the light beam cannot be neglected in boiling experiments. Since the exit angle of the light beam is known, this allows to account for the deflections and phase shifts outside the boundary layer (in the bulk fluid and in the windows of the test cell). Three dimensional temperature distributions in the vicinity of the bubble are reconstructed using tomographic techniques. In tomography, the measurement volume is sliced into 2D planes. In the present study these planes are parallel to the heated surface. The objective is to determine the values of the field parameter of interest in form of the field function in these 2D planes. The field parameter is the change of the refractive index of the liquid in the measurement volume caused by temperature changes. By superimposing data for many 2D planes recorded at the same time instant, the 3D temperature distribution in the measurement volume is recovered.

Herman, Cila↗
NASA NTRS1998

The Athena Mars Rover Science Payload

The Mars Surveyor missions that will be launched in April of 2001 will include a highly capable rover that is a successor to the Mars Pathfinder mission's Sojourner rover. The design goals for this rover are a total traverse distance of at least 10 km and a total lifetime of at least one Earth year. The rover's job will be to explore a site in Mars' ancient terrain, searching for materials likely to preserve a record of ancient martian water, climate, and possibly biology. The rover will collect rock and soil samples, and will store them for return to Earth by a subsequent Mars Surveyor mission in 2005. The Athena Mars rover science payload is the suite of scientific instruments and sample collection tools that will be used to perform this job. The specific science objectives that NASA has identified for the '01 rover payload are to: (1) Provide color stereo imaging of martian surface environments, and remotely-sensed point discrimination of mineralogical composition. (2) Determine the elemental and mineralogical composition of martian surface materials. (3) Determine the fine-scale textural properties of these materials. (4) Collect and store samples. The Athena payload has been designed to meet these objectives. The focus of the design is on field operations: making sure the rover can locate, characterize, and collect scientifically important samples in a dusty, dirty, real-world environment. The topography, morphology, and mineralogy of the scene around the rover will be revealed by Pancam/Mini-TES, an integrated imager and IR spectrometer. Pancam views the surface around the rover in stereo and color. It uses two high-resolution cameras that are identical in most respects to the rover's navigation cameras. The detectors are low-power, low-mass active pixel sensors with on-chip 12-bit analog-to-digital conversion. Filters provide 8-12 color spectral bandpasses over the spectral region from 0.4 to 1.1 micron Narrow-angle optics provide an angular resolution of 0.28 mrad/pixel, nearly a factor of four higher than that of the Mars Pathfinder and Mars Surveyor '98 cameras. Image compression will be performed using a wavelet compression algorithm. The Mini-Thermal Emission Spectrometer (Mini-TES) is a point spectrometer operating in -the thermal IR. It produces high spectral resolution (5 /cm) image cubes with a wavelength range of 5-40 gm, a nominal signal/noise ratio of 500:1, and a maximum angular resolution of 7 mrad (7 cm at a distance of 10 in). The wavelength region over which it operates samples the diagnostic fundamental absorption features of rockforming minerals, and also provides some capability to see through dust coatings that could tend to obscure spectral features. The mineralogical information that Mini-TES provides will be used to select from a distance the rocks and soils that will be investigated in more detail and ultimately sampled. Mini-TES is derived from the MO/MGS TES instrument, but is significantly smaller and simpler. The instrument uses an 8-cm Cassegrain telescope, a Michelson interferometer, and uncooled pyroelectric detectors. Along with its mineralogical capabilities, Mini-TES can provide information on the thermophysical properties of rocks and soils. Viewing upward, it can also provide temperature profiles through the martian atmospheric boundary layer. Elemental and Mineralogical Composition: Once promising samples have been identified from a distance using Pancam/Mini-TES, they will be studied in detail using up to three compositional sensors that can be placed directly against them by an Instrument Arm. The two compositional sensors, presently on the payload are an Alpha-Proton-X-Ray Spectrometer (APXS), and a Mossbauer Spectrometer. The APXS is derived closely from the instrument that flew on Mars Pathfinder. Radioactive alpha sources and three detection modes (alpha, proton, and x-ray) provide elemental abundances of rocks and soils to complement and constrain mineralogical data. The Athena APXS will have a revised mechanical design that will cut down significantly on backscattering of alpha particles from martian atmospheric carbon. It will also include a target of known elemental composition that will be used for calibration purposes. The Athena Mossbauer Spectrometer is a diagnostic instrument for the mineralogy and oxidation state of Fe-bearing phases, which are particularly important on Mars. The instrument measures the resonant absorption of gamma rays produced by a Co-57 source to determine splitting of nuclear energy levels in Fe atoms that is related to the electronic environment surrounding them. It has been under development for space flight for many years at the Technical University of Darmstadt. The Mossbauer Spectrometer (and the other arm instruments) will be able to view a small permanent magnet array that will attract magnetic particles in the martian soil. The payload may also include a Raman Spectrometer. If included, the Raman Spectrometer will provide precise identification of major and minor mineral phases. It requires no sample preparation, and is also sensitive to organics. Fine-Scale Texture: The Instrument Arm a also carries a Microscopic Imager that will obtain high-resolution monochromatic images of the same materials for which compositional data will be obtained. Its spatial resolution is 20 micron/pixel over a 1 cm depth of field, and 40 micron/pixel over a 1-cm depth of field. Like Pancam, it uses the same active pixel sensor detectors and electronics as the rover's navigation cameras. The Instrument Arm is a three degree-of-freedom arm that uses designs and components from the Mars Pathfinder and Mars Surveyor '98 projects. Its primary function is instrument positioning. Along with the instruments noted above, it also carries a brush that can be used to remove dust and other loose coatings from rocks. Sample Collection and Storage: Martian rock and soil samples will be collected using a low-power rotary coring drill called the Mini-Corer. An important characteristic of this device is that it can obtain intact samples of rock from up to 5 cm within strong boulders and bedrock, Nominal core dimensions are 8xl7 mm. The Mini-Corer drills a core to the commanded depth in a rock, shears it off, retains it, and extracts it. It can also acquire samples of loose soil, using soil sample cups that are pressed downward into loose material. The Mini-Corer can drill at angles from vertical to 45' off vertical. It has six interchangeable bits for long life. Mechanical damage to the sample during drilling is minimal, and heating is negligible. After acquisition, the sample may be viewed by the arm instruments, and/or placed in one of 104 compartments in the Sample Container. A subset of the acquired samples may be replaced with other samples obtained later if desired. The Sample Container has no moving parts, and is mounted external to the rover for easy removal by the Mars Surveyor 2005 flight system. Operation of the rover will make extensive use of automated onboard navigation and hazard avoidance capabilities. Otherwise, use of onboard autonomy is minimal. Data downlink capability is about 40 Mbit/sol, and the use of the Mars Surveyor '01 orbiter for data relay imposes a limit of at most two command cycles per sol. Because of the significant amount of time available between command cycles, all payload elements will be operated sequentially, rather than in parallel.; this approach also significantly simplifies operations and minimizes peak power usage. The landing site for the '01 rover has not been selected yet. Site selection will make as full use as possible of Mars Global Surveyor data, and will involve substantial input from the broad Mars science community. Summary: The following table describes the mass, power, providers, and key scientific objectives of all the major elements of the Athena payload. Additional Athena payload information may be found at: http://astrosun.tn.cornell.edu/athena/index.html. Additional information contained in the original.

Squyes, S. W.↗