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 1,009 records · Page 56

An automatic multigrid method for the solution of sparse linear systems

An automatic version of the multigrid method for the solution of linear systems arising from the discretization of elliptic PDE's is presented. This version is based on the structure of the algebraic system solely, and does not use the original partial differential operator. Numerical experiments show that for the Poisson equation the rate of convergence of our method is equal to that of classical multigrid methods. Moreover, the method is robust in the sense that its high rate of convergence is conserved for other classes of problems: non-symmetric, hyperbolic (even with closed characteristics) and problems on non-uniform grids. No double discretization or special treatment of sub-domains (e.g. boundaries) is needed. When supplemented with a vector extrapolation method, high rates of convergence are achieved also for anisotropic and discontinuous problems and also for indefinite Helmholtz equations. A new double discretization strategy is proposed for finite and spectral element schemes and is found better than known strategies.

Shapira, Yair↗

Atmospheric Correction at AERONET Locations: A New Science and Validation Data Set

This paper describes an AERONET-based Surface Reflectance Validation Network (ASRVN) and its dataset of spectral surface bidirectional reflectance and albedo based on MODIS TERRA and AQUA data. The ASRVN is an operational data collection and processing system. It receives 50x50 square kilometer subsets of MODIS L1B data from MODAPS and AERONET aerosol and water vapor information. Then it performs an accurate atmospheric correction for about 100 AERONET sites based on accurate radiative transfer theory with high quality control of the input data. The ASRVN processing software consists of L1B data gridding algorithm, a new cloud mask algorithm based on a time series analysis, and an atmospheric correction algorithm. The atmospheric correction is achieved by fitting the MODIS top of atmosphere measurements, accumulated for 16-day interval, with theoretical reflectance parameterized in terms of coefficients of the LSRT BRF model. The ASRVN takes several steps to ensure high quality of results: 1) cloud mask algorithm filters opaque clouds; 2) an aerosol filter has been developed to filter residual semi-transparent and sub-pixel clouds, as well as cases with high inhomogeneity of aerosols in the processing area; 3) imposing requirement of consistency of the new solution with previously retrieved BRF and albedo; 4) rapid adjustment of the 16-day retrieval to the surface changes using the last day of measurements; and 5) development of seasonal back-up spectral BRF database to increase data coverage. The ASRVN provides a gapless or near-gapless coverage for the processing area. The gaps, caused by clouds, are filled most naturally with the latest solution for a given pixels. The ASRVN products include three parameters of LSRT model (k(sup L), k(sup G), k(sup V)), surface albedo, NBRF (a normalized BRF computed for a standard viewing geometry, VZA=0 deg., SZA=45 deg.), and IBRF (instantaneous, or one angle, BRF value derived from the last day of MODIS measurement for specific viewing geometry) for MODIS 500m bands 1-7. The results are produced daily at resolution of 1 km in gridded format. We also provide cloud mask, quality flag and a browse bitmap image. The new dataset can be used for a wide range of applications including validation analysis and science research.

Wang, Yujie↗

Applying a Space-Based Security Recovery Scheme for Critical Homeland Security Cyberinfrastructure Utilizing the NASA Tracking and Data Relay (TDRS) Based Space Network

Protection of the national infrastructure is a high priority for cybersecurity of the homeland. Critical infrastructure such as the national power grid, commercial financial networks, and communications networks have been successfully invaded and re-invaded from foreign and domestic attackers. The ability to re-establish authentication and confidentiality of the network participants via secure channels that have not been compromised would be an important countermeasure to compromise of our critical network infrastructure. This paper describes a concept of operations by which the NASA Tracking and Data Relay (TDRS) constellation of spacecraft in conjunction with the White Sands Complex (WSC) Ground Station host a security recovery system for re-establishing secure network communications in the event of a national or regional cyberattack. Users would perform security and network restoral functions via a Broadcast Satellite Service (BSS) from the TDRS constellation. The BSS enrollment only requires that each network location have a receive antenna and satellite receiver. This would be no more complex than setting up a DIRECTTV-like receiver at each network location with separate network connectivity. A GEO BSS would allow a mass re-enrollment of network nodes (up to nationwide) simultaneously depending upon downlink characteristics. This paper details the spectrum requirements, link budget, notional assets and communications requirements for the scheme. It describes the architecture of such a system and the manner in which it leverages off of the existing secure infrastructure which is already in place and managed by the NASAGSFC Space Network Project.

Cybersecurity↗

Double ion production in mercury thrusters

The development of a model which predicts doubly charged ion density is discussed. The accuracy of the model is shown to be good for two different thruster sizes and a total of 11 different cases. The model indicates that in most cases more than 80% of the doubly charged ions are produced from singly charged ions. This result can be used to develop a much simpler model which, along with correlations of the average plasma properties, can be used to determine the doubly charged ion density in ion thrusters with acceptable accuracy. Two different techniques which can be used to reduce the doubly charged ion density while maintaining good thruster operation, are identified as a result of an examination of the simple model. First, the electron density can be reduced and the thruster size then increased to maintain the same propellant utilization. Second, at a fixed thruster size, the plasma density, temperature and energy can be reduced and then to maintain a constant propellant utilization the open area of the grids to neutral propellant loss can be reduced through the use of a small hole accelerator grid.

Peters, R. R.↗

ERS-1 scatterometer calibration and validation activities at ECMWF. B: From radar backscatter characteristics to wind vector solutions

Calibration and validation activities for the ERS-1 scatterometer were carried out at ECMWF (European Center for Medium range Weather Forecast) complementary to the 'Haltenbanken' field campaign off the coast of Norway. At a Numerical Weather Prediction (NWP) center a wealth of verifying data is available both in time and space. This data is used to redefine the wind retrieval procedure given the instrumental characteristics. It was found that a maximum likelihood estimation procedure to obtain the coefficients of a reformulated sigma deg to wind relationship should use radar measurements in logarithmic rather than physical space, and use winds as the wind components rather than wind speed and direction. Doing this, a much more accurate transfer function than the one currently operated by ESA was derived. Sigma deg measurement space shows no signature of a separation in an upwind solution cone and a downwind solution cone. As such signature was anticipated in ESA's wind direction ambiguity removal algorithm, reconsideration of the procedure is necessary. Despite the fact that revisions have to be made in the process of wind retrieval; a grid potential is shown for scatterometry in meteorology and climatology.

Stoffelen, AD↗

Hurricane Directional Wave Spectrum Spatial Variation in the Open Ocean and at Landfall

The sea surface directional wave spectrum was measured for the first time in all quadrants of a hurricane in open water using the NASA airborne scanning radar altimeter (SRA) carried aboard one of the NOAA WP-3D hurricane hunter aircraft at 1.5 km height. The SRA measures the energetic portion of the directional wave spectrum by generating a topographic map of the sea surface. At 8 Hz, the SRA sweeps a radar beam of 1' half-power width (two-way) across the aircraft ground track over a swath equal to 0.8 of the aircraft height, simultaneously measuring the backscattered power at its 36 GHz (8.3 mm) operating frequency and the range to the sea surface at 64 positions. These slant ranges are multiplied by the cosine of the off-nadir angles to determine the vertical distances from the aircraft to the sea surface. Subtracting these distances from the aircraft height produces the sea surface elevation map. The sea surface topography is interpolated to a uniform grid, transformed by a two dimensional FFT, and Doppler corrected. The open-ocean data were acquired on 24 August 1998 when hurricane Bonnie was east of the Bahamas and moving toward 330 deg at about 5 m/s. Individual waves up to 18 m height were observed and the spatial variation of the wave field was dramatic. The dominant waves generally propagated at significant angles to the downwind direction. At some positions there were three different wave fields of comparable energy crossing each other. The NOAA aircraft spent over five hours within 180 km of the eye, and made five eye penetrations. On 26 August 1998, the NOAA aircraft flew at 2.2 km height when hurricane Bonnie was making landfall near Wilmington, NC, documenting the directional wave spectrum in the region between Charleston, SC and Cape Hatteras, NC. The aircraft flight lines included segments near and along the shoreline as well as far offshore. Animations of the directional wave spectrum spatial variation along the aircraft tracks on the two flights will be presented using a 100: 1 time compression.

Walsh, Edward J.↗

Fpack and Funpack Utilities for FITS Image Compression and Uncompression

Fpack is a utility program for optimally compressing images in the FITS (Flexible Image Transport System) data format (see http://fits.gsfc.nasa.gov). The associated funpack program restores the compressed image file back to its original state (as long as a lossless compression algorithm is used). These programs may be run from the host operating system command line and are analogous to the gzip and gunzip utility programs except that they are optimized for FITS format images and offer a wider choice of compression algorithms. Fpack stores the compressed image using the FITS tiled image compression convention (see http://fits.gsfc.nasa.gov/fits_registry.html). Under this convention, the image is first divided into a user-configurable grid of rectangular tiles, and then each tile is individually compressed and stored in a variable-length array column in a FITS binary table. By default, fpack usually adopts a row-by-row tiling pattern. The FITS image header keywords remain uncompressed for fast access by FITS reading and writing software. The tiled image compression convention can in principle support any number of different compression algorithms. The fpack and funpack utilities call on routines in the CFITSIO library (http://hesarc.gsfc.nasa.gov/fitsio) to perform the actual compression and uncompression of the FITS images, which currently supports the GZIP, Rice, H-compress, and PLIO IRAF pixel list compression algorithms.

Pence, W.↗

Optical and Microphysical Retrievals of Marine Stratocumulus Clouds off the Coast of Namibia from Satellite and Aircraft

Though the emphasis of the Southern Africa Regional Science Initiative 2000 (SAFARI-2000) dry season campaign was largely on emission sources and transport, the assemblage of aircraft (including the high altitude NASA ER-2 remote sensing platform and the University of Washington CV-580, UK MRF C-130, and South African Weather Bureau JRA in situ aircrafts) provided a unique opportunity for cloud studies. Therefore, as part of the SAFARI initiative, investigations were undertaken to assess regional aerosol-cloud interactions and cloud remote sensing algorithms. In particular, the latter part of the experiment concentrated on marine boundary layer stratocumulus clouds off the southwest coast of Africa. Associated with cold water upwelling along the Benguela current, the Namibian stratocumulus regime has received limited attention but appears to be unique for several reasons. During the dry season, outflow of continental fires and industrial pollution over this area can be extreme. From below, upwelling provides a rich nutrient source for phytoplankton (a source of atmospheric sulfur through DMS production as well as from decay processes). The impact of these natural and anthropogenic sources on the microphysical and optical properties of the stratocumulus is unknown. Continental and Indian Ocean cloud systems of opportunity were also studied during the campaign. SAFARI 2000 aircraft flights off the coast of Namibia were coordinated with NASA Terra Satellite overpasses for synergy with the Moderate Resolution Imaging Spectroradiometer (MODIS) and other Terra instruments. MODIS was developed by NASA and launched onboard the Terra spacecraft on December 18, 1999 (and Aqua spacecraft on May 4, 2002). Among the remote sensing algorithms developed and applied to this sensor are cloud optical and microphysical properties that include cloud thermodynamic phase, optical thickness, and effective particle radius of both liquid water and ice clouds. The archived products from these algorithms have applications in climate change studies, climate modeling, numerical weather prediction, and fundamental atmospheric research. The archived MODIS Collection 5 cloud products processing stream will be used to analyze low water cloud scenes off the Namibian and Angolan coasts during SAFARI 2000 time period, as well as other years. Pixel-level Terra and Aqua MODIS retrievals (l. km spatial resolution at nadir) and gridded (1' uniform grid) statistics of cloud optical thickness and effective particle radius will be presented, including joint probability distributions between the two quantities. In addition, perspectives from the MODIS Airborne Simulator, which flew on the ER-2 during SAFARI 2000 providing high spatial resolution retrievals (50 m at nadir), will be presented as appropriate. The H-SAF Program requires an experimental operational European-centric Satellite Precipitation Algorithm System (E-SPAS) that produces medium spatial resolution and high temporal resolution surface rainfall and snowfall estimates over the Greater European Region including the Greater Mediterranean Basin. Currently, there are various types of experimental operational algorithm methods of differing spatiotemporal resolutions that generate global precipitation estimates. This address will first assess the current status of these methods and then recommend a methodology for the H-SAF Program that deviates somewhat from the current approach under development but one that takes advantage of existing techniques and existing software developed for the TRMM Project and available through the public domain.

Platnick, Steven E.↗

Evaluation of a High-Resolution Numerical Weather Prediction Model’s Simulated Clouds Using Observations from CloudSat, GOES-13 and in Situ Aircraft

This study aimed to assess tropical cloud properties predicted by Environment and Climate Change Canada’s Global Environmental Multiscale (GEM)modelwhen run with the Milbrandt–Yau double-moment cloud microphysical scheme and one-way nesting that culminated at a (∼300 km)2 inner domain with 0.25 km horizontal grid spacing. The assessment utilized satellite and in situ data collected during the High Ice Water Content (HIWC) and High Altitude Ice Crystals (HAIC) projects for a mesoscale convective system on 16May 2015 over French Guiana.Data fromCloud- Sat’s cloud-profiling radar and GOES-13’s imager were compared to data either simulated directly by GEM or produced by operating on GEM’s cloud data with both the CFMIP (Cloud Feedback Model Intercomparison Project) Observation Simulator Package (COSP) instrument simulator and a three-dimensional Monte Carlo solar radiative transfer model. In situ observations were made from research aircraft – Canada’s National Research Council Convair-580 and the French SAFIRE Falcon-20 – whose flight paths were aligned with CloudSat’s ground-track. Spatial and temporal shifts of clouds simulated by GEM compared well to GOES-13 imagery. There are, however, differences between simulated and observed amounts of high and low cloud. While GEM did well at predicting ranges of ice-water content (IWC) near 11 km altitude (Falcon-20), it produces too much graupel and snow near 7 km (Convair-580). This produced large differences between CloudSat’s and COSP-generated radar reflectivities and two-way attenuations. On the other hand, CloudSat’s inferred values of IWC agree well with in situ samples at both altitudes. Generally, GEM’s visible reflectances exceeded GOES-13’s on account of having produced too much low-level liquid cloud. It is expected that GEM’s disproportioning of cloud hydrometeors will improve once it includes a better representation of secondary ice production.

Qu, Zhipeng↗

Southern Colorado Disasters: Using NASA Observations to Map Aspen Extent and Recovery Due to Wildfire

Quaking aspen (Populus tremuloides) is an important species for wildlife, watershed health, and ecosystem resilience across its range. Heavy ungulate browsing and factors influenced by a changing climate including seasonal temperature changes and moisture deficit have led to reduced post-fire aspen regeneration rates in southern Colorado. This project partnered with Trinchera Ranch and the Colorado State Forest Service to estimate aspen recovery after the Spring Creek Fire, which ignited in June of 2018. The Southern Colorado Disasters team utilized field measurements and satellite imagery from Landsat Operational Land Imager (OLI), Sentinel-2 MultiSpectral Instrument (MSI), and the Shuttle Radar Topography Mission (SRTM) to train and run several random forest models that detect pre- and post-fire aspen extent. Ocular sampling of over 500 points on high-resolution pre-fire and post-fire images identified percentage aspen cover in 30 x 30-meter grid cells. This process provided training data for regression models, which were able to detect aspen across the landscape for both time periods using multiple remote sensing vegetation health indices. In addition, landscape suitability for aspen regeneration was modeled to provide a guide for managers on where to monitor for aspen regeneration post-fire.

DEVELOP Project Summary↗

Southern Colorado Disasters: Using NASA Earth Observations to Map Aspen Extent and Recovery Due to Wildfire

Quaking aspen (Populus tremuloides) is an important species for wildlife, watershed health, and ecosystem resilience across its range. Heavy ungulate browsing and factors influenced by a changing climate including seasonal temperature changes and moisture deficit have led to reduced post-fire aspen regeneration rates in southern Colorado. This project partnered with Trinchera Ranch and the Colorado State Forest Service to estimate aspen recovery after the Spring Creek Fire, which ignited in June of 2018. The Southern Colorado Disasters team utilized field measurements and satellite imagery from Landsat Operational Land Imager (OLI), Sentinel-2 MultiSpectral Instrument (MSI), and the Shuttle Radar Topography Mission (SRTM) to train and run several random forest models that detect pre- and post-fire aspen extent. Ocular sampling of over 500 points on high-resolution pre-fire and post-fire images identified percentage aspen cover in 30 x 30-meter grid cells. This process provided training data for regression models, which were able to detect aspen across the landscape for both time periods using multiple remote sensing vegetation health indices. In addition, landscape suitability for aspen regeneration was modeled to provide a guide for managers on where to monitor for aspen regeneration post-fire.

DEVELOP Technical Paper↗

Newly Released GPCP Version 3.2 Global Precipitation Datasets at NASA GES DISC

The Global Precipitation Climatology Project (GPCP) is the precipitation component of an internationally coordinated set of (mainly) satellite-based global products dealing with the Earth's water and energy cycles, under the auspices of the Global Water and Energy Experiment (GEWEX) Data and Assessment Panel (GDAP) of the World Climate Research Program. As the follow-on to the GPCP Version 2.X products, GPCP Version 3 (GPCP V3.2) seeks to continue the production of long, homogeneous precipitation record using modern input and calibration datasets. The GPCP V3.2 provides globally complete analyses of surface precipitation on a 0.5°x 0.5° latitude/longitude grid at both monthly and daily intervals, respectively covering 1983 to the present and June 2000 to the present. New data fields have been introduced to better characterize the precipitation, particularly including an estimate of the fraction of the precipitation that is liquid (rain) in both the Monthly and Daily, and a Quality Index for the Monthly. Compared to the operational GPCP V2.3 Monthly, the V3.2 Monthly provides a more reasonable climatology in the Southern Ocean, and increases the global average precipitation by about 4.46%, which is in line with recommendations of recent assessments. However, the two versions have comparable global and regional trends for 1983-2020. Compared to the operational One-Degree Daily (Version 1.3) product, the V3.2 Daily better represents the histogram of precipitation rates, particularly at high values. In this presentation, we will present the latest GPCP V3.2 daily and monthly datasets archived and distributed at NASA Goddard Earth Sciences (GES) Data and Information Services Center (DISC) along with examples from NASA Giovanni.

Precipitation↗

Scalar Atomic Defect-Based Solid-State Self-calibrating Magnetometer (3SM) for Space Plasma Analysis

The Earth’s magnetosphere is a system of multiple, co-located particle populations interacting via plasma waves. Things to understand are driving processes, radiation belt and ring current issues, auroral physics, internal plasma processes, and magnetosphere-ionosphere mapping issues. These plasma-physics processes enable the Earth’s magnetosphere to evolve in response to temporal changes in the solar wind and they underlie the phenomena of space weather, which impacts spacecraft systems, astronauts, radio communications, and ground based electric-power grids. The objective of 3SM is to measure magnetic field strength and to calibrate the Vector Magnetometer (VM) device to maintain absolute accuracy during a mission. The 3SM features make the instrument preferably suited not only for the traditional role of scalar magnetometers as absolute references for the calibration of the on-board vector instruments, but also for extended operational capacities, such as higher frequency scalar measurements (of potential interest for magnetosphere studies for the low frequency part of the spectrum) or autonomous scalar / vector operations. Diamond has been the solid-state platform of choice for quantum device technologies for some time, however, it suffers from difficulties such as scalability, integration, and cost. While the diamond platform is very useful for quantum technologies, further development is needed to make it viable. The material platform of choice for NASA Glenn’s Quantum Sensing And Spin Physics (Q-SASP) is silicon carbide (SiC). This is due to the much higher industry development of the SiC material platform for high-power and high-temperature electronics. It leverages both the decades-long SiC development expertise and infrastructure at NASA Glenn and its growing capabilities in quantum metrology. To make SiC devices usable for quantum technologies such as quantum sources, a much deeper understanding of defects is needed. Q-SASP is developing quantum metrology capabilities to evaluate the energy structure, defect formation energy, band structure augmentation, generation/recombination rates, and limits of dipole-dipole coupling in non-metal implanted SiC devices. This can be achieved by analysis of zero-field splitting, low-field resonance, and singlet-triplet mixing through various forms of Electrically Detectable Magnetic Resonance (EDMR) and Near-Zero Field Magnetic Resonance (NZFMR) spectroscopy. This work will discuss recent system developments, device developments, computational modeling, and spectroscopy results and analysis of defects created by non-metal implantations in SiC devices. The defect formation energies of V Si ,V C , V C V Si , N C V Si , N Si , N C in 4H-SiC are previously reported values in other research [2]-[3]. The defect formation energies of PSi and PC were calculated in GPAW [fig 1A]. The basic underlying mechanism of the zero-field phenomenon is the mixing of singlet and triplet states [4]-[5]. In most spin-dependent transport, two electron spins are involved, and thus one must consider each of their interactions with the field. We investigated the electronic and magnetic properties of 4H-SiC and 6H-SiC. The defect formation energy helps us determine what types of defects we are observing in the SiC EDMR experiment. They have very low formation energy (it is negative). The phosphorus substitution in 4H-SiC is a very stable defect. The band diagrams provide us with vital information about how the electronic properties of SiC (such as band gap) change as we add non-metal defects. The zero-field splitting parameters allow us to study the inflection point in the NZFMR [fig 1B]. We clearly observed zero-field splitting. We also noted that the zero-field splitting remained constant with changing bias. We aspect it zero-field splitting to remain constant while the hyperfine and exchange interaction perturbations shift under the influence of an external magnetic field. This is the essence of quantum magnetometry and self-calibration.

space plasma↗

The Geoscience Laser Altimeter System (GLAS) for the ICESAT Mission

The Laser In space Technology Experiment, Shuttle Laser Altimeter and the Mars Observer Laser Altimeter have demonstrated accurate measurements of atmospheric backscatter and Surface heights from space. The recent MOLA measurements of the Mars surface have 40 cm vertical resolution and have reduced the global uncertainty in Mars topography from a few km to about 5 m. The Geoscience Laser Altimeter System (GLAS) is a next generation lidar for Earth orbit being developed as part of NASA's Icesat Mission. The GLAS design combines a 10 cm precision surface lidar with a sensitive dual wavelength cloud and aerosol lidar. GLAS will precisely measure the heights of the Earth's polar ice sheets, establish a grid of accurate height profiles of the Earth's land topography, and profile the vertical backscatter of clouds and aerosols on a global scale. GLAS is being developed to fly on a small dedicated spacecraft in a polar orbit with a 590 630 km altitude at inclination of 94 degrees. GLAS is scheduled to launch in the summer 2001 and to operate continuously for a minimum of 3 years with a goal of 5 years. The primary mission for GLAS is to measure the seasonal and annual changes in the heights of the Greenland and Antarctic ice sheets. GLAS will continuously measure the vertical distance from orbit to the Earth's surface with 1064 nm pulses from a ND:YAG laser at a 40 Hz rate. Each 5 nsec wide laser pulse is used to produce a single range measurement, and the laser spots have 66 m diameter and about 170 m center-center spacings. When over land GLAS will profile the heights of the topography and vegetation. The GLAS receiver uses a 1 m diameter telescope and a Si APD detector. The detector signal is sampled by an all digital receiver which records each surface echo waveform with I nsec resolution and a stored echo record lengths of either 200, 400, or 600 samples. Analysis of the echo waveforms within the instrument permits discrimination between cloud and surface echoes. Ground based echo analysis permits precise ranging, determining the roughness or slopes of the surface as well as the vertical distributions of vegetation illuminated by the laser. Accurate knowledge of the laser beam's pointing angle is needed to prevent height biases when over sloped surfaces. For surfaces with 2 deg. slopes, knowledge of pointing angle of the beam's centroid to about 8 urad is needed to achieve 10 cm height accuracy. GLAS uses a stellar reference system (SRS) to determine the pointing angle of each laser firing relative to inertial space. The SRS uses a high precision star camera oriented toward local zenith and a gyroscope to determine the inertial orientation of the SRS optical bench. The far field pattern of each laser is measured pulse relative to the star camera with a laser reference system (LRS). Optically measuring each laser far field pattern relative to the orientation of the star camera and gyroscope permits the precise pointing angle of each laser pulse to be determined. GLAS will also determine the vertical distributions of clouds and aerosols by measuring the vertical profile of laser energy backscattered by the atmosphere at both 1064 and 532 nm. The 1064 nm measurements use the Si APD detector and profile the height and vertical structure of thicker clouds. The measurements at 532 nm use new highly sensitive photon counting, detectors, and measure the height distributions of very thin Clouds and aerosol layers. With averaging these can be used to determine the height of the planetary boundary layer. The instrument design and expected performance will be discussed.

Abshire, James B.↗

Calibration and Application of Plume Radiation Code: PARRAD

The next generation of Reusable Launch Vehicle (RLV) design requires careful considerations towards Thermal Protection System (TPS) not only in the high-heating regions, also in the base region. The use of hydrocarbon based liquid propulsion engine and the construction, launch and operational cost limitations places a significant burden on the base region TPS design. The radiation from the rocket-plume is expected to dominate base heating during ascent. The objective of this paper is to evaluate the accuracy of a plume radiation modeling with the PARRAD code. The original PARRAD code has been modified to work with single and multiple /overset grids and plume-field solution generated by computational codes, such as GASP. These modifications will be detailed in the proposed paper. To establish accuracy, simulations of the Space Shuttle Main Engine nozzle were performed and the radiative heating predicted by PARRAD are compared with ground based measurements. These results and pertinent discussion of issues related to accuracy will be detailed. In addition, predictions for the X-34 RLV, yet to be performed, will be presented and the significance (or lack) of base heating due to plume radiation will be discussed.

Venkatapathy, Ethiraj↗

Application of a Detached Eddy Simulation Approach with Finite-Rate Chemistry to Mars-Relevant Retropropulsion Operating Environments

Human-scale Mars vehicles will require retropropulsion for descent and landing, replacing heritage supersonic parachute systems with an extended phase of powered flight. Due to the limitations of terrestrial testing in Mars-relevant conditions, design and analysis will increasingly rely on computational modeling and simulation. This paper provides an overview of a computational campaign investigating the aerodynamics of a Mars lander concept along various points on a powered descent trajectory including supersonic, transonic, and subsonic conditions using finite-rate chemistry. Simulations using unstructured grids containing billions of elements are performed at scale using thousands of Graphics Processing Units, enabling run-times of a few days for each simulation presented. At each freestream condition, significant minor species concentrations are observed external to the nozzles in the large mixing region upstream of the vehicle. While the flowfields are highly non-stationary in all cases, the mean integrated forces and moments on the vehicle remain small in comparison to the deceleration provided through retropropulsion.

Ashley Korzun↗

Internal erosion rates of a 10-kW xenon ion thruster

A 30 cm diameter divergent magnetic field ion thruster, developed for mercury operation at 2.7 kW, was modified and operated with xenon propellant at a power level of 10 kW for 567 h to evaluate thruster performance and lifetime. The major differences between this thruster and its baseline configuration were elimination of the three mercury vaporizers, use of a main discharge cathode with a larger orifice, reduction in discharge baffle diameter, and use of an ion accelerating system with larger acceleration grid holes. Grid thickness measurement uncertainties, combined with estimates of the effects of reactive residual facility background gases gave a minimum screen grid lifetime of 7000 h. Discharge cathode orifice erosion rates were measured with three different cathodes with different initial orifice diameters. Three potential problems were identified during the wear test: the upstream side of the discharge baffle eroded at an unacceptable rate; two of the main cathode tubes experienced oxidation, deformation, and failure; and the accelerator grid impingement current was more than an order of magnitude higher than that of the baseline mercury thruster. The charge exchange ion erosion was not quantified in this test. There were no measurable changes in the accelerator grid thickness or the accelerator grid hole diameters.

Rawlin, Vincent K.↗

Internal erosion rates of a 10-kW xenon ion thruster

A 30 cm diameter divergent magnetic field ion thruster, developed for mercury operation at 2.7 kW, was modified and operated with xenon propellant at a power level of 10 kW for 567 h to evaluate thruster performance and lifetime. The major differences between this thruster and its baseline configuration were elimination of the three mercury vaporizers, use of a main discharge cathode with a larger orifice, reduction in discharge baffle diameter, and use of an ion accelerating system with larger acceleration grid holes. Grid thickness measurement uncertainties, combined with estimates of the effects of reactive residual facility background gases gave a minimum screen grid lifetime of 7000 h. Discharge cathode orifice erosion rates were measured with three different cathodes with different initial orifice diameters. Three potential problems were identified during the wear test: the upstream side of the discharge baffle eroded at an unacceptable rate; two of the main cathode tubes experienced oxidation, deformation, and failure; and the accelerator grid impingement current was more than an order of magnitude higher than that of the baseline mercury thruster. The charge exchange ion eroison was not quantified in this test. There were no measurable changes in the accelerator grid thickness or the accelerator grid hole diameters.

Rawlin, Vincent K.↗