Search NASA⌕ Search

SEARCH · Search NASA

Results for “temperature 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.

At least 1,729 records · Page 96

Assessment of Terra MODIS Thermal Emissive Band Calibration Using Cold Targets and Measurements in Lunar Roll Events

Terra MODIS has provided continuous global observations for science research and applications for more than 18 years. The MODIS Thermal emissive bands (TEB) radiometric calibration uses a quadratic function for instrument response. The calibration coefficients are updated using the response of an on-board blackbody (BB) in quarterly warm-up and cool-down (WUCD) events. As instrument degradation and electronic crosstalk of long-wave infrared (LWIR) bands 27 to 30 developed substantial issues, accurate calibration is crucial for a high-quality L1B product. The on-board BB WUCD temperature ranges from 270 K to 315 K and the derived nonlinear response has a relatively large uncertainty for the offset, especially for these LWIR bands, which affects the measurements of low brightness temperature (BT) scenes. In this study, the TEB radiometric calibration impact on the L1B product is assessed using selected cold targets and the measurements during regular lunar rolls. The cold targets include Antarctic Dome Concordia (Dome-C) and deep convective clouds (DCC) for the calibration assessment, focusing on bands 27 to 30. Dome-C area is covered with uniformly-distributed permanent snow, and the atmospheric effect is small and relatively constant. Usually the DCC is treated as an invariant earth target to evaluate the reflective solar band calibration. The DCC can also be treated as a stable target to assess the performance of TEB calibration. During a scheduled lunar observation event with a spacecraft roll maneuver to view the moon through the space view port, the instrument cavity provides a stable reference for calibration assessment. The long-term trending of BT measurements and the relative difference between scan mirror sides and detectors are used for the assessment of the calibration consistency and stability. The comparison of L1B products over the selected targets before and after the calibration coefficients update can be used to assess the impact of a calibration look-up table (LUT) update. This assessment is beneficial for future calibration algorithm and LUT update procedure improvements for enhancing the L1B product quality.

Chang, Tiejun↗

The Algorithm for MODIS Wavelength On-Orbit Calibration using the SRCA

The Spectro-Radiometric Calibration Assembly (SRCA) provides on-orbit spectral calibration of the MODerate resolution Imaging Spectroradiometer (MODIS) reflected solar bands and this paper describes how it is accomplished. The SRCA has two adjacent exit slits: 1) Main slit and 2) Calibration slit. The output from the main slit is measured by a reference silicon photo-diode (SIPD) and then passes through the MODIS. The output from the calibration slit passes through a piece of didymium transmission glass and then it is measured by a calibration SIPD. The centroids of the sharp spectral peaks of a didyrnium glass are utilized as wavelength standards. After normalization using the reference SIPD signal to eliminate the effects of the illuminating source spectra, the calibration SIPD establishes the relationship between the peaks of the didymium spectra and the grating angle; this is accomplished through the grating equation. In the grating equation the monochromator parameters, beta (half angle between the incident and diffractive beams) and theta(sub ff) (offset angle of the grating motor) are determined by matching, in a least square sense, the known centroid wavelengths of the didymium peaks and the calculated centroid grating angles from the calibration SIPD signals for the peaks. A displacement between the calibration SIPD and the reference SIPD complicates the signal processing. The response of the MODIS bands to the SRCA output is also normalized by the reference SIPD signal to eliminate the effect of the source spectrum. That response differs from what is measured by MODIS using the Spectral Measurements Assembly (SpMA), a laboratory double-monochromator, due to the wider slit width of the SRCA. The SRCA slit function, calculated using measurements by the SPMA and the SRCA at the same MODIS temperature, will be used to recover the MODIS spectral response using the SRCA data measured on-orbit.

Montgomery, Harry↗

ARM Research in the Equatorial Western Pacific: A Decade and Counting

The tropical western Pacific (TWP) is an important climatic region. Strong solar heating, warm sea surface temperatures, and the annual progression of the intertropical convergence zone (ITCZ) across this region generate abundant convective systems, which through their effects on the heat and water budgets have a profound impact on global climate and precipitation. In order to accurately evaluate tropical cloud systems in models, measurements of tropical clouds, the environment in which they reside, and their impact on the radiation and water budgets are needed. Because of the remote location, ground-based datasets of cloud, atmosphere, and radiation properties from the TWP region have come primarily from short-term field experiments. While providing extremely useful information on physical processes, these short-term datasets are limited in statistical and climatological information. To provide longterm measurements of the surface radiation budget in the tropics and the atmospheric properties that affect it, the Atmospheric Radiation Measurement program established a measurement site on Manus Island, Papua New Guinea, in 1996 and on the island republic of Nauru in late 1998. These sites provide unique datasets now available for more than 10 years on Manus and Nauru. This article presents examples of the scientific use of these datasets including characterization of cloud properties, analysis of cloud radiative forcing, model studies of tropical clouds and processes, and validation of satellite algorithms. New instrumentation recently installed at the Manus site will provide expanded opportunities for tropical atmospheric science.

tropical regions↗

EHW Approach to Temperature Compensation of Electronics

Efforts are under way to apply the concept of evolvable hardware (EHW) to compensate for variations, with temperature, in the operational characteristics of electronic circuits. To maintain the required functionality of a given circuit at a temperature above or below the nominal operating temperature for which the circuit was originally designed, a new circuit would be evolved; moreover, to obtain the required functionality over a very wide temperature range, there would be evolved a number of circuits, each of which would satisfy the performance requirements over a small part of the total temperature range. The basic concepts and some specific implementations of EHW were described in a number of previous NASA Tech Briefs articles, namely, "Reconfigurable Arrays of Transistors for Evolvable Hardware" (NPO-20078), Vol. 25, No. 2 (February 2001), page 36; Evolutionary Automated Synthesis of Electronic Circuits (NPO- 20535), Vol. 26, No. 7 (July 2002), page 37; "Designing Reconfigurable Antennas Through Hardware Evolution" (NPO-20666), Vol. 26, No. 7 (July 2002), page 38; "Morphing in Evolutionary Synthesis of Electronic Circuits" (NPO-20837), Vol. 26, No. 8 (August 2002), page 31; "Mixtrinsic Evolutionary Synthesis of Electronic Circuits" (NPO-20773) Vol. 26, No. 8 (August 2002), page 32; and "Synthesis of Fuzzy-Logic Circuits in Evolvable Hardware" (NPO-21095) Vol. 26, No. 11 (November 2002), page 38. To recapitulate from the cited prior articles: EHW is characterized as evolutionary in a quasi-genetic sense. The essence of EHW is to construct and test a sequence of populations of circuits that function as incrementally better solutions of a given design problem through the selective, repetitive connection and/or disconnection of capacitors, transistors, amplifiers, inverters, and/or other circuit building blocks. The connection and disconnection can be effected by use of field-programmable transistor arrays (FPTAs). The evolution is guided by a search-andoptimization algorithm (in particular, a genetic algorithm) that operates in the space of possible circuits to find a circuit that exhibits an acceptably close approximation of the desired functionality. The evolved circuits can be tested by mathematical modeling (that is, computational simulation) only, tested in real hardware, or tested in combinations of computational simulation and real hardware.

Stoica, Adrian↗

Large Scale Characteristics and Variability of the Global Sea Ice Cover

More than two decades of satellite passive microwave data are used to study and evaluate the large scale characteristics and the changing state of the sea ice cover in both the Northern and Southern Hemispheres. Satellite data provide day/night almost continuous observation of global sea ice cover thereby enabling quantitative variability studies at various time scales. Despite coarse sensor resolution, spatial detail is provided through the use of sea ice concentrations which are derived using an algorithm that determines the fraction of ice and open water within each satellite footprint. Large seasonal fluctuations in the extent are apparent with those of the Southern Hemisphere having larger amplitudes but less symmetrical seasonal distribution than those of the Northern Hemisphere. The large scale interannual variability of the ice cover has been evaluated globally as well as regionally and in the Northern Hemisphere, the yearly anomaly maps show a predominance of positive values in the 1980s and negative values in the 1990s. Regression analysis show that the ice extent and ice area are on a decline at the rate of -2.0 +/- 0.5% and -3.1 +/- 0.3% per decade, respectively, in the Northern Hemisphere but there are regions like the Bering Sea with positive trends. What is intriguing, however, is that the perennial sea ice cover has been declining at a much faster rate than for the entire hemisphere, i.e., 6.7 +/- 2.4% and 8.3 +/- 2.4 % per decade for ice extent and ice area, respectively. The perennial ice cover consists mainly of thick multiyear ice floes, and its persistent decline would mean a reduction in the average thickness of sea ice and a change in the overall characteristics of the Arctic sea ice cover. Furthermore, the yearly anomaly patterns are coherent with those of surface temperatures derived from 19 years of thermal infrared AVHRR data. The latter also shows that in consolidated ice regions, the average temperature during summer minima has been increasing at about 0.9 +/- 0.6 K per decade.

Comiso, Josefino C.↗

Assimilation of Gridded Terrestrial Water Storage Observations from GRACE into a Land Surface Model

Observations of terrestrial water storage (TWS) from the Gravity Recovery and Climate Experiment (GRACE) satellite mission have a coarse resolution in time (monthly) and space (roughly 150,000 km(sup 2) at midlatitudes) and vertically integrate all water storage components over land, including soil moisture and groundwater. Data assimilation can be used to horizontally downscale and vertically partition GRACE-TWS observations. This work proposes a variant of existing ensemble-based GRACE-TWS data assimilation schemes. The new algorithm differs in how the analysis increments are computed and applied. Existing schemes correlate the uncertainty in the modeled monthly TWS estimates with errors in the soil moisture profile state variables at a single instant in the month and then apply the increment either at the end of the month or gradually throughout the month. The proposed new scheme first computes increments for each day of the month and then applies the average of those increments at the beginning of the month. The new scheme therefore better reflects submonthly variations in TWS errors. The new and existing schemes are investigated here using gridded GRACE-TWS observations. The assimilation results are validated at the monthly time scale, using in situ measurements of groundwater depth and soil moisture across the U.S. The new assimilation scheme yields improved (although not in a statistically significant sense) skill metrics for groundwater compared to the open-loop (no assimilation) simulations and compared to the existing assimilation schemes. A smaller impact is seen for surface and root-zone soil moisture, which have a shorter memory and receive smaller increments from TWS assimilation than groundwater. These results motivate future efforts to combine GRACE-TWS observations with observations that are more sensitive to surface soil moisture, such as L-band brightness temperature observations from Soil Moisture Ocean Salinity (SMOS) or Soil Moisture Active Passive (SMAP). Finally, we demonstrate that the scaling parameters that are applied to the GRACE observations prior to assimilation should be consistent with the land surface model that is used within the assimilation system.

GRACE↗

Progress Toward Generation of a Navier-Stokes Database for a Harrier in Ground Effect

The Harrier YAV-8B aircraft is capable of vertical and short-field take-off and landing (V/STOL) by directing its four exhaust nozzles toward the ground, or conventional flight by rotating its nozzles into a horizontal position. The British Royal Air Force and the United States Marine Corps have used this aircraft for more than 30 years to provide a quick reaction time for troop support, and reduce the need for long runways. The success of this powered-lift (PL) vehicle has also prompted the more recent design of the Joint Strike Fighter (JSF). However there are significant safety issues that must be addressed when operating a PL vehicle in close proximity to the ground. Hot Gas Ingestion (HGI) by the inlets can result in a rapid loss of powered lift; and high-speed jet flows along the ground plane can induce low pressures underneath the vehicle, causing a 'suck-down' effect. Under these conditions, departure from controlled flight may occur. Moreover, unsteady ground vortices and jet fountains can affect the aircraft,s controllability and its proximity to ground troops. The viscous, time-dependent flow fields of PL vehicles are difficult to accurately and efficiently predict using Computational Fluid Dynamics (CFD). A number of researchers have used the time-dependent Reynolds-averaged Navier-Stokes (RANS) equations to compute flows for single and multiple jets in a cross-flow. A few have added some geometric complexity to the problem by computing flows for jet-augmented delta wings near a ground plane. Smith et.al. computed for the first time a single RANS solution about a simplified Harrier. This geometry included a fuselage, wing, leading edge root extension (LERX), inlets, and exhaust nozzles. All of these investigations cite two practical problems with computing these flows: 1) the need for improved solution accuracy; and, 2) the need for faster solution methods. We view the need for faster solution methods as key to improving the solution accuracy and making this class of computation more routine. One can hardly refine grids, explore the use of advanced turbulence models, and generate databases when it takes weeks of dedicated computer time for a single solution. Chaderjian, Ahmad, Pandya, and Murman have focused on reducing the time-to-solution for this very difficult and complex problem through process automation and exploitation of parallel computing. They began with the Harrier geometry reported, and added a deflected wing flap and empennage for greater realism. To date more than 80 solutions have been carried out. This paper will describe this process and progress made in reducing the time required to generate a simple longitudinal force and moment database for a Harrier in ground effect. It shows a typical snap-shot from an unsteady streakline animation, where fluid particles are colored by temperature. The ground vortex and a jet-fountain vortex are highlighted. It also shows a similar streakline image, where HGI occurs due to the vehicle in close proximity to the ground. It is show the mean lift coefficient as a function of angle of attack and height. The angle of attack range was 4 deg less than or = alpha less than or = 10 deg with an increment of 1 degree, and the height range was 10 ft less than or = h less than or = 30ft with an increment of 5 feet. This 35 solution database was extended to over 2500 cases using a monotone cubic-spline interpolation procedure. The suck-down effect (reduction of lift near the ground) is highlighted in the figure. The "cushion effect," the conventional reduction of lift as the vehicle moves out of ground effect, is also indicated. All 35 RANS solutions were obtained using 952 Silicon Graphics Origin 2000 and 3000 processors in dedicated mode for one week. Typically, 112 processors were assigned to each case. Some other cases used fewer processors to utilize all available CPUS. The final paper will report on the automation of the solution process, including: grid generation, job monitoring, solution completion criteria, and post processing. Moreover, improvements in parallel efficiency for a dual time-step algorithm for the RANS equations will also be presented. Results will be discussed in detail using unsteady streakline flow visualization to correlate unsteady flow structures with dominant aerodynamic frequencies. The stability derivatives, CL, and CL, will also be presented.

Chaderjian, Neal M.↗

Validation of POLDER/ADEOS data using a ground-based lidar network: Preliminary results for semi-transparent and cirrus clouds

At mid and tropical latitudes, cirrus clouds are present more than 50% of the time in satellites observations. Due to their large spatial and temporal coverage, and associated low temperatures, cirrus clouds have a major influence on the Earth-Ocean-Atmosphere energy balance through their effects on the incoming solar radiation and outgoing infrared radiation. At present the impact of cirrus clouds on climate is well recognized but remains to be asserted more precisely, for their optical and radiative properties are not very well known. In order to understand the effects of cirrus clouds on climate, their optical and radiative characteristics of these clouds need to be determined accurately at different scales in different locations i.e. latitude. Lidars are well suited to observe cirrus clouds, they can detect very thin and semi-transparent layers, and retrieve the clouds geometrical properties i.e. altitude and multilayers, as well as radiative properties i.e. optical depth, backscattering phase functions of ice crystals. Moreover the linear depolarization ratio can give information on the ice crystal shape. In addition, the data collected with an airborne version of POLDER (POLarization and Directionality of Earth Reflectances) instrument have shown that bidirectional polarized measurements can provide information on cirrus cloud microphysical properties (crystal shapes, preferred orientation in space). The spaceborne version of POLDER-1 has been flown on ADEOS-1 platform during 8 months (October 96 - June 97), and the next POLDER-2 instrument will be launched in 2000 on ADEOS-2. The POLDER-1 cloud inversion algorithms are currently under validation. For cirrus clouds, a validation based on comparisons between cloud properties retrieved from POLDER-1 data and cloud properties inferred from a ground-based lidar network is currently under consideration. We present the first results of the validation.

Chepfer, H.↗

Assessing Global Surface Water Inundation Dynamics Using Combined Satellite Information from SMAP, AMSR2 and Landsat

A method to assess global land surface water (fw) inundation dynamics was developed by exploiting the enhanced fw sensitivity of L-band (1.4 GHz) passive microwave observations from the Soil Moisture Active Passive (SMAP) mission. The L-band fw (fw(sub LBand)) retrievals were derived using SMAP H-polarization brightness temperature (Tb) observations and predefined L-band reference microwave emissivities for water and land endmembers. Potential soil moisture and vegetation contributions to the microwave signal were represented from overlapping higher frequency (Tb) observations from AMSR2. The resulting (fw(sub LBand)) global record has high temporal sampling (1-3 days) and 36-km spatial resolution. The (fw(sub LBand)) annual averages corresponded favourably (R=0.84, p<0.001) with a 250-m resolution static global water map (MOD44W) aggregated at the same spatial scale, while capturing significant inundation variations worldwide. The monthly (fw(sub LBand)) averages also showed seasonal inundation changes consistent with river discharge records within six major US river basins. An uncertainty analysis indicated generally reliable (fw(sub LBand)) performance for major land cover areas and under low to moderate vegetation cover, but with lower accuracy for detecting water bodies covered by dense vegetation. Finer resolution (30-m) (fw(sub LBand)) results were obtained for three sub-regions in North America using an empirical downscaling approach and ancillary global Water Occurrence Dataset (WOD) derived from the historical Landsat record. The resulting 30-m (fw(sub LBand)) retrievals showed favourable spatial accuracy for water (70.71%) and land (98.99%) classifications and seasonal wet and dry periods when compared to independent water maps derived from Landsat-8 imagery. The new (fw(sub LBand)) algorithms and continuing SMAP and AMSR2 operations provide for near real-time, multi-scale monitoring of global surface water inundation dynamics and potential flood risk.

WOD↗

Assessment of GOES-16/ABI middle wave infrared band using references of Himawari-8/AHI and Aqua/MODIS

GOES-16 is the first of the GOES-R series of Geostationary Operational Environmental Satellites (GOES) and was launched on November 19, 2016. The spacecraft was initially in a test position of 89.5° West and reached its operational position (75.2° West) on December 11, 2017. The Himawari-8 spacecraft was launched on October 7, 2014 and is located at 140.7º East. The similar design and similar calibration algorithm between the Advanced Baseline Imager (ABI) on-board GOES-16 and the Advanced Himawari Imager (AHI) on board Himawari-8 makes the importance of inter-comparison. Due to their locations, double difference is an appropriate method for their comparison and Aqua MODIS is one of good references. However, ABI (AHI) midwave-infrared (MWIR) band 7 does not have good matching with Aqua MODIS. In this work, the ABI-AHI comparison and ABI assessment for MWIR band 7 is performed using Aqua bands 20, 22, and 23. The ocean sites under ABI (AHI) at nadir are used for inter-comparison with Aqua MODIS. To enhance the comparison accuracy, a few procedures and corrections have been applied. For MWIR band 7, the ABI-AHI difference is about -0.39K over ocean scene, with the ABI measurement precision being slightly better than that of AHI. This double difference method is also being used for the assessment of ABI consistency before and after re-location on November 30, 2017. Two ocean scenes are selected for the ABI re-location assessments, with measurement precision at nadir providing better measurements than non-zero view angles. The ABI MWIR measurement over the ocean scene at the same view angle before and after re-location shows that the precisions are comparable. The MWIR band brightness temperature (BT) measurement over ocean scene shows a 0.04K difference, while the measurement precision before and after re-location is consistent.

Inter-comparison↗

How Much Global Burned Area Can Be Forecast on Seasonal Time Scales Using Sea Surface Temperatures?

Large-scale sea surface temperature (SST) patterns influence the interannual variability of burned area in many regions by means of climate controls on fuel continuity, amount, and moisture content. Some of the variability in burned area is predictable on seasonal timescales because fuel characteristics respond to the cumulative effects of climate prior to the onset of the fire season. Here we systematically evaluated the degree to which annual burned area from the Global Fire Emissions Database version 4 with small fires (GFED4s) can be predicted using SSTs from 14 different ocean regions. We found that about 48 of global burned area can be forecast with a correlation coefficient that is significant at a p < 0.01 level using a single ocean climate index (OCI) 3 or more months prior to the month of peak burning. Continental regions where burned area had a higher degree of predictability included equatorial Asia, where 92% of the burned area exceeded the correlation threshold, and Central America, where 86% of the burned area exceeded this threshold. Pacific Ocean indices describing the El Nino-Southern Oscillation were more important than indices from other ocean basins, accounting for about 1/3 of the total predictable global burned area. A model that combined two indices from different oceans considerably improved model performance, suggesting that fires in many regions respond to forcing from more than one ocean basin. Using OCI-burned area relationships and a clustering algorithm, we identified 12 hotspot regions in which fires had a consistent response to SST patterns. Annual burned area in these regions can be predicted with moderate confidence levels, suggesting operational forecasts may be possible with the aim of improving ecosystem management.

Chen, Yang↗

On the X-ray spectrum of Kepler's supernova remnant

We have devised a method to do nonequilibrium ionization calculations on the results of two-dimensional hydrodynamical simulations, based on the algorithm of Hughes & Helfand (1985). We have calculated the ionization structure and X-ray emission for a two-dimensional numerical hydrodynamical simulation for the remnant of Kepler's supernova (SN); the hydrodynamical model was presented in a previous paper. In this model, the progenitor of Kepler's SN is assumed to have been a massive runaway star ejected from the Galactic plane. In its red supergiant stage, its dense stellar wind was distorted and compressed into a bow shock by the ram pressure of the tenuous interstellar medium. The subsequent interaction of the supernova ejecta with this asymmetric circumstellar matter produced a strongly asymmetric supernova remnant (SNR). In this paper, we present calculated X-ray spectra for this hydrodynamical model. A comparison with observations implies a moderate overabundance of Fe in Kepler's SNR (only 50% larger than its cosmic value), in contrast to a large (6 to 15) Fe overabundance derived previously. However, we confirm earlier conclusions that Si and S abundances are 2 to 3 times solar. These modest enhancements of Si, S, and Fe may be attributed either to heavy-element enriched SN ejecta or to the initial chemical abundances of the SN progenitor, which originated in the metal-rich inner Galaxy. The comparison of our models with the observed spectra confirm theoretical predictions that moderate electron heating occurs at strong collisionless shock fronts, with the implied electron/mean temperature ratio of approximately 0.5.

Borkowski, Kazimierz J.↗

Exploring and Visualizing A-Train Instrument Data

The succession of US and international satellites that follow each other in close succession, known as the A-Train, affords an opportunity to atmospheric researchers that no single platform could provide: Increasing the number of observations at any given geographic location.. . a more complete "virtual science platform". However, vertically and horizontally, co-registering and regridding datasets from independently developed missions, Aqua, Calipso, Cloudsat, Parasol, and Aura, so that they can be inter-compared can be daunting to some, and may be repeated by many. Scientists will individually spend much of their time and resources acquiring A-Train datasets of interest residing at various locations, developing algorithms to match up and graph datasets along the A-Train track, and search through large amounts of data for areas and/or phenomena of interest. The aggregate amount of effort that can be expended on repeating pre-science tasks could climb into the tens of millions of dollars. The goal of the A-Train Data Depot (ATDD) is to enable free movement of remotely located A-Train data so that they are combined to create a consolidated vertical view of the Earth's Atmosphere along the A-Train tracks. The innovative approach of analyzing and visualizing atmospheric profiles along the platforms track (i.e., time) is accomplished by through the ATDDs Giovanni data analysis and visualization tool. Giovanni brings together data from Aqua (MODIS, AIRS, AMSR-E), Cloudsat (cloud profiling radar) and Calipso (CALIOP, IIR), as well as the Aura (OMI, MLS, HIRDLS, TES) to create a consolidated vertical view of the Earth's Atmosphere along the A-Train tracks. This easy to learn and use exploration tool will allow users to create vertical profiles of any desired A-Train dataset, for any given time of choice. This presentation shows the power of Giovanni by describing and illustrating how this tool facilitates and aids A-Train science and research. A web based display system Giovanni provides users with the capability of creating co-located profile images of temperature and humidity data from the MODIS, MLS and AIRS instruments for a user specified time and spatial area. In addition, Cloud and Aerosol profiles may also be displayed for the Cloudsat and Caliop instruments. The ability to modify horizontal and vertical axis range, data range and dynamic color range is also provided. Two dimensional strip plots of MODIS, AIRS, OM1 and POLDER parameters, co-located along the Cloudsat reference track, can also be plotted along with the Cloudsat cloud profiling data. Center swath pixels for the same parameters can also be shown as line plots overlaying the Cloudsat or Calipso profile images. Images and subsetted data produced in each analysis run may be downloaded. Users truly can explore and discover data specific to their needs prior to ever transferring data to their analysis tools.

Kempler, S.↗

Comparative Study on High-Order Positivity-preserving WENO Schemes

In gas dynamics and magnetohydrodynamics flows, physically, the density and the pressure p should both be positive. In a standard conservative numerical scheme, however, the computed internal energy is obtained by subtracting the kinetic energy from the total energy, resulting in a computed p that may be negative. Examples are problems in which the dominant energy is kinetic. Negative may often emerge in computing blast waves. In such situations the computed eigenvalues of the Jacobian will become imaginary. Consequently, the initial value problem for the linearized system will be ill posed. This explains why failure of preserving positivity of density or pressure may cause blow-ups of the numerical algorithm. The adhoc methods in numerical strategy which modify the computed negative density and/or the computed negative pressure to be positive are neither a conservative cure nor a stable solution. Conservative positivity-preserving schemes are more appropriate for such flow problems. The ideas of Zhang & Shu (2012) and Hu et al. (2012) precisely address the aforementioned issue. Zhang & Shu constructed a new conservative positivity-preserving procedure to preserve positive density and pressure for high-order WENO schemes by the Lax-Friedrichs flux (WENO/LLF). In general, WENO/LLF is too dissipative for flows such as turbulence with strong shocks computed in direct numerical simulations (DNS) and large eddy simulations (LES). The new conservative positivity-preserving procedure proposed in Hu et al. (2012) can be used with any high-order shock-capturing scheme, including high-order WENO schemes using the Roe's flux (WENO/Roe). The goal of this study is to compare the results obtained by non-positivity-preserving methods with the recently developed positivity-preserving schemes for representative test cases. In particular the more difficult 3D Noh and Sedov problems are considered. These test cases are chosen because of the negative pressure/density most often exhibited by standard high-order shock-capturing schemes. The simulation of a hypersonic nonequilibrium viscous shock tube that is related to the NASA Electric Arc Shock Tube (EAST) is also included. EAST is a high-temperature and high Mach number viscous nonequilibrium flow consisting of 13 species. In addition, as most common shock-capturing schemes have been developed for problems without source terms, when applied to problems with nonlinear and/or sti source terms these methods can result in spurious solutions, even when solving a conservative system of equations with a conservative scheme. This kind of behavior can be observed even for a scalar case (LeVeque & Yee 1990) as well as for the case consisting of two species and one reaction (Wang et al. 2012). For further information concerning this issue see (LeVeque & Yee 1990; Griffiths et al. 1992; Lafon & Yee 1996; Yee et al. 2012). This EAST example indicated that standard high-order shock-capturing methods exhibit instability of density/pressure in addition to grid-dependent discontinuity locations with insufficient grid points. The evaluation of these test cases is based on the stability of the numerical schemes together with the accuracy of the obtained solutions.

WENO↗

166 GHz Ice Scattering Signal in Cold-Season Precipitating Cloud Systems

Surface snowfall can be produced by a variety of forcing mechanisms that are accompanied by different in-cloud physical processes and microphysical composition. For instance, shallow cumuliform snowfall forced by marine cold air outbreaks shows distinct geographic and seasonal signatures and plays in important role in annual snowfall production in many regions worldwide. Coincident CloudSat-CALIPSO observations confirm that oceanic regions exhibiting frequent shallow cumuliform snow are typically composed of mixed-phase cloud and precipitation particles and display passive microwave brightness temperature behavior that is dominated by cloud liquid water emission, while combined Global Precipitation Measurement (GPM) radiometer and surface radar studies of intense shallow cumuliform events indicate significant ice particle scattering signals at the highest GPM radiometer frequencies. Passive microwave imager channels in the 37 to 89 GHz range have long been used to detect precipitation-sized ice particle scattering signatures in clouds. Dual polarization plays an important role in distinguishing between scattering signatures and other background variability. The GPM Microwave Imager (GMI) is the first satellite radiometer to possess dual-polarized channels at 166 GHz. These high-frequency channels are more susceptible to emission and attenuation by cloud liquid water and water vapor than lower-frequency channels, but are also more sensitive to ice scattering. This research studies the ice scattering information contained in the 166 GHz channels compared to the historically used 85/89 GHz channels on previous sensors. A 166 GHz scattering index (SI) will be developed to understand how its correlation with non-precipitating ice scattering at higher altitudes also affects the 166 GHz scattering signature in cold-season precipitating cloud systems. The 166 GHz SI will be used to possibly isolate precipitation scattering signature for different snowfall modes (shallow convective vs. deep stratiform). The supercooled liquid water affecting shallow convective systems and drastically reducing the scattering signal could be detected and quantified, providing additional information for convective snow detection and quantification in solid precipitation retrieval algorithms.

Lisa Milani↗

Atmospheric Airborne Pressure Measurements Using the Oxygen A Band for the ASCENDS Mission

We report on airborne atmospheric pressure measurements using new fiber-based laser technology and the oxygen A-band at 765 nm. Remote measurements of atmospheric temperature and pressure are required for a number of NASA Earth science missions and specifically for the Active Sensing of CO2 Emissions Over Nights, Days, and Seasons (ASCENDS) mission. Accurate measurements of tropospheric CO2 on a global scale are very important in order to better understand its sources and sinks and to improve predictions on any future climate change. The ultimate goal of a CO2 remote sensing mission, such as ASCENDS, is to derive the CO2 concentration in the atmosphere in terms of mole fraction in unit of parts-per-million (ppmv) with regard to dry air. Therefore, both CO2 and the dry air number of molecules in the atmosphere are needed in deriving this quantity. O2 is a stable molecule and uniformly mixed in the atmosphere. Measuring the O2 absorption in the atmosphere can thus be used to infer the dry air number of molecules and then used to calculate CO2 concentration. With the knowledge of atmospheric water vapor, we can then estimate the total surface pressure needed for CO2 retrievals. Our work, funded by the ESTO IIP program, uses fiber optic technology and non-linear optics to generate 765 nm laser radiation coincident with the Oxygen A-band. Our pulsed, time gated technique uses several on- and off-line wavelengths tuned to the O2 absorption line. The choice of wavelengths allows us to measure the pressure by using two adjacent O2 absorptions in the Oxygen A-band. Our retrieval algorithm fits the O2 lineshapes and derives the pressure. Our measurements compare favorably with a local weather monitor mounted outside our laboratory and a local weather station.

Riris, Haris↗

Atmospheric Airborne Pressure Measurements Using the Oxygen A Band for the ASCENDS Mission

We report on airborne atmospheric pressure measurements using new fiber-based laser technology and the oxygen A-band at 765 nm. Remote measurements of atmospheric temperature and pressure are required for a number of NASA Earth science missions and specifically for the Active Sensing of CO2 Emissions Over Nights, Days, and Seasons (ASCENDS) mission. Accurate measurements of tropospheric CO2 on a global scale are very important in order to better understand its sources and sinks and to improve predictions on any future climate change. The ultimate goal of a CO2 remote sensing mission, such as ASCENDS, is to derive the CO2 concentration in the atmosphere in terms of mole fraction in unit of parts-per-million (ppmv) with regard to dry air. Therefore, both CO2 and the dry air number of molecules in the atmosphere are needed in deriving this quantity. O2 is a stable molecule and uniformly mixed in the atmosphere. Measuring the O2 absorption in the atmosphere can thus be used to infer the dry air number of molecules and then used to calculate CO2 concentration. With the knowledge of atmospheric water vapor, we can then estimate the total surface pressure needed for CO2 retrievals. Our work, funded by the ESTO IIP program, uses fiber optic technology and non-linear optics to generate 765 nm laser radiation coincident with the Oxygen A-band. Our pulsed, time gated technique uses several on- and off-line wavelengths tuned to the O2 absorption line. The choice of wavelengths allows us to measure the pressure by using two adjacent O2 absorptions in the Oxygen A-band. Our retrieval algorithm fits the O2 lineshapes and derives the pressure. Our measurements compare favorably with a local weather monitor mounted outside our laboratory and a local weather station.

Riris, Haris↗

Validation of CFD/Heat Transfer Software for Turbine Blade Analysis

I am an intern in the Turbine Branch of the Turbomachinery and Propulsion Systems Division. The division is primarily concerned with experimental and computational methods of calculating heat transfer effects of turbine blades during operation in jet engines and land-based power systems. These include modeling flow in internal cooling passages and film cooling, as well as calculating heat flux and peak temperatures to ensure safe and efficient operation. The branch is research-oriented, emphasizing the development of tools that may be used by gas turbine designers in industry. The branch has been developing a computational fluid dynamics (CFD) and heat transfer code called GlennHT to achieve the computational end of this analysis. The code was originally written in FORTRAN 77 and run on Silicon Graphics machines. However the code has been rewritten and compiled in FORTRAN 90 to take advantage of more modem computer memory systems. In addition the branch has made a switch in system architectures from SGI's to Linux PC's. The newly modified code therefore needs to be tested and validated. This is the primary goal of my internship. To validate the GlennHT code, it must be run using benchmark fluid mechanics and heat transfer test cases, for which there are either analytical solutions or widely accepted experimental data. From the solutions generated by the code, comparisons can be made to the correct solutions to establish the accuracy of the code. To design and create these test cases, there are many steps and programs that must be used. Before a test case can be run, pre-processing steps must be accomplished. These include generating a grid to describe the geometry, using a software package called GridPro. Also various files required by the GlennHT code must be created including a boundary condition file, a file for multi-processor computing, and a file to describe problem and algorithm parameters. A good deal of this internship will be to become familiar with these programs and the structure of the GlennHT code. Additional information is included in the original extended abstract.

Kiefer, Walter D.↗