Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Localized basis sets”

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 163 records · Page 9

Hydrogen in disordered titania: connecting local chemistry, structure, and stoichiometry through accelerated exploration

Hydrogen incorporation in native surface oxides of metal alloys often controls the onset of metal hydriding, with implications for materials corrosion and hydrogen storage. A key representative example is titania, which forms as a passivating layer on a variety of titanium alloys for structural and functional applications. These oxides tend to be structurally diverse, featuring polymorphic phases, grain boundaries, and amorphous regions that generate a disparate set of unique local environments for hydrogen. Here, we introduce a workflow that can efficiently and accurately navigate this complexity. First, a machine learning force field, trained on ab initio molecular dynamics simulations, was used to generate amorphous configurations. Density functional theory calculations were then performed on these structures to identify local oxygen environments, which were compared against experimental observations. Second, to classify subtle differences across the disordered configuration space, we employ a graph-based sampling procedure. Finally, local hydrogen binding energies and hopping kinetics are computed using exhaustive density functional theory calculations on representative configurations. Here, we leverage this methodology to show that hydrogen binding energetics are described by local oxygen coordination, which in turn is affected by stoichiometry, and form the basis of hopping kinetics and diffusion. Together these results imply that hydrogen incorporation and transport in TiO x can be tailored through compositional engineering, with implications for improving the performance and durability of titanium-derived alloys in hydrogen environments.

36 MATERIALS SCIENCE↗

Feature Detection and Curve Fitting Using Fast Walsh Transforms for Shock Tracking: Applications

Walsh functions form an orthonormal basis set consisting of square waves. Square waves make the system well suited for detecting and representing functions with discontinuities. Given a uniform distribution of 2p cells on a one-dimensional element, it has been proven that the inner product of the Walsh Root function for group p with every polynomial of degree < or = (p - 1) across the element is identically zero. It has also been proven that the magnitude and location of a discontinuous jump, as represented by a Heaviside function, are explicitly identified by its Fast Walsh Transform (FWT) coefficients. These two proofs enable an algorithm that quickly provides a Weighted Least Squares fit to distributions across the element that include a discontinuity. The detection of a discontinuity enables analytic relations to locally describe its evolution and provide increased accuracy. Time accurate examples are provided for advection, Burgers equation, and Riemann problems (diaphragm burst) in closed tubes and de Laval nozzles. New algorithms to detect up to two C0 and/or C1 discontinuities within a single element are developed for application to the Riemann problem, in which a contact discontinuity and shock wave form after the diaphragm bursts.

Gnoffo, Peter A.↗

Impacts of a Prescribed Fire on Air Quality in Central New Mexico

A short-duration but high-impact air quality event occurred on 28 November 2018 along the Rio Grande Valley of New Mexico. This fire occurred outside the typical wildfire season, and greatly impacted the air quality in Socorro, NM, and the surroundings. Measurements were taken during the event using an aerosol light scattering technique (integrating nephelometer) and a particulate mass concentration monitor (DustTrak PM optical monitor). The instruments sampled the ambient air during the event on the campus of the New Mexico Institute of Mining and Technology in Socorro, New Mexico. The peak values on a 5-min basis of light scattering and the PM mass concentration reached 470 Mm -1 and 270 µg/m 3 , respectively. We examined the meteorological context of the event using local meteorological data and back trajectories using the NOAA HYSPLIT model to determine atmospheric transport and possible sources. Several fires, both prescribed and wildfires, occurred in the region including a prescribed burn at Bosque del Apache National Wildlife Refuge (17 km south-southeast of the receptor site). The data suggest that the prescribed burn at Bosque del Apache was the dominant contributor due to transport evidence and the event’s narrow spatiotemporal extent. The increasing importance of restoring ecosystem function using prescribed fire in wildland fire management will likely lead to more frequent air quality impacts and sets up policy tradeoffs that require a balance between these public goals. This study examines the evidence of the effects of a prescribed fire in a protected wildland area impacting the air quality in a nearby populated area.

54 ENVIRONMENTAL SCIENCES↗

3D Source Reconstruction Using Coded Aperture Gamma-Ray Imaging

Recent measurements with coded-aperture imagers demonstrate material mass determination in a holdup setting to an accuracy within a few percent. This capability is of particular interest to the Surplus Plutonium Disposition (SPD) project, which aims to dilute and dispose of surplus plutonium oxide. Gamma-ray imagers can be used to determine holdup without interrupting normal operations. In this work, we examine techniques for 3D source localization and mass determination using gamma-ray imagers. Coded-aperture imagers provide excellent source localization within the 2D image plane; however, multiple imagers operating in tandem are necessary to identify source location in 3D space. A Maximum Likelihood Expectation-Maximization (MLEM) method for fitting detector mappings is a powerful tool for accomplishing this task. MLEM allows 3D source localization to be simultaneously constrained using multiple gamma-ray imagers by constructing the basis for the MLEM fit using detector mappings from different detector locations stitched together. Each of these basis points represents a singular response from a source in 3D space and is generated using Monte Carlo simulations of sources placed individually at different locations throughout the imager’s field of view. Additionally, implementing knowledge of the physical equipment in the simulations of the glovebox used for the SPD project incorporates attenuation effects that are needed to calculate material holdup.

Laminack, Alex↗

A Pulse Generation Framework with Augmented Program-aware Basis Gates and Criticality Analysis

Near-term intermediate scale quantum (NISQ) de- vices are subject to considerable noise and short coherence time. Consequently, it is critical to minimize circuit execution latency. Traditionally, each basis gate of a transpiled circuit is decoded into a fixed episode of the device control pulses. Recently, people started to investigate merged pulse generation for customized gates through quantum optimal control (QOC). However, existing QOC approaches face the challenges of (i) restricted search space due to prohibitive compilation overhead; (ii) suboptimal end-to-end performance due to aggressive local optimization and falsely introduced dependency among the customized gates; (iii) inadequate adaptivity towards system calibration, which is critical for NISQ devices. In this work, we propose PAQOC, a novel QOC framework that can (i) automatically detect frequently encountered gate patterns in the logical circuit by modeling the problem as a subgraph mining process and reuse these patterns to enable much larger search space exploration (i.e., program aware); (ii) systemically construct customized gate-set based on the impact to the overall program latency (i.e., criticality-aware); and (iii) quickly adapt to system re-calibration thanks to the small-scale pattern-based gate generation (i.e., adaptivity-aware). PAQOC achieves a good tradeoff between circuit performance and compilation time, allowing fully automatic, single stop, ad- hoc customized pulse generation for more efficient execution of user programs on NISQ devices. Evaluations using fifteen applications show that PAQOC can achieve on average 1.95× speedup of the circuit latency and achieve on average 36.7% reduction in compilation overhead. With PAQOC, circuits can run faster with reduced noise, allowing deeper circuits to be tested within the coherence time of present NISQ platforms.

Chen, Yanhao↗

Improved Compression of Wavelet-Transformed Images

A recently developed data-compression method is an adaptive technique for coding quantized wavelet-transformed data, nominally as part of a complete image-data compressor. Unlike some other approaches, this method admits a simple implementation and does not rely on the use of large code tables. A common data compression approach, particularly for images, is to perform a wavelet transform on the input data, and then losslessly compress a quantized version of the wavelet-transformed data. Under this compression approach, it is common for the quantized data to include long sequences, or runs, of zeros. The new coding method uses prefixfree codes for the nonnegative integers as part of an adaptive algorithm for compressing the quantized wavelet-transformed data by run-length coding. In the form of run-length coding used here, the data sequence to be encoded is parsed into strings consisting of some number (possibly 0) of zeros, followed by a nonzero value. The nonzero value and the length of the run of zeros are encoded. For a data stream that contains a sufficiently high frequency of zeros, this method is known to be more effective than using a single variable length code to encode each symbol. The specific prefix-free codes used are from two classes of variable-length codes: a class known as Golomb codes, and a class known as exponential-Golomb codes. The codes within each class are indexed by a single integer parameter. The present method uses exponential-Golomb codes for the lengths of the runs of zeros, and Golomb codes for the nonzero values. The code parameters within each code class are determined adaptively on the fly as compression proceeds, on the basis of statistics from previously encoded values. In particular, a simple adaptive method has been devised to select the parameter identifying the particular exponential-Golomb code to use. The method tracks the average number of bits used to encode recent runlengths, and takes the difference between this average length and the code parameter. When this difference falls outside a fixed range, the code parameter is updated (increased or decreased). The Golomb code parameter is selected based on the average magnitude of recently encoded nonzero samples. The coding method requires no floating- point operations, and more readily adapts to local statistics than other methods. The method can also accommodate arbitrarily large input values and arbitrarily long runs of zeros. In practice, this means that changes in the dynamic range or size of the input data set would not require a change to the compressor. The algorithm has been tested in computational experiments on test images. A comparison with a previously developed algorithm that uses large code tables (generated via Huffman coding on training data) suggests that the data-compression effectiveness of the present algorithm is comparable to the best performance achievable by the previously developed algorithm.

Kiely, Aaron↗

Evaluation of the spatial and temporal measurement requirements of remote sensors for monitoring regional air pollution episodes

The spatial and temporal measurement requirements of satellite sensors for monitoring regional air pollution episodes were evaluated. Use was made of two sets of data from the Sulfate Regional Experiment (SURE), which provided the first ground-based aerosol measurements from a regional-scale station network. The sulfate data were analyzed for two air pollution episode cases. The results of the analysis indicate that the key considerations required for episode mapping from satellite sensors are the following: (1) detection of sulfate levels exceeding 20 micron-g/cu m; (2) capability to view a broad area (of the order of 1500 km swath) because of regional extent of pollution episodes; (3) spatial resolution sufficient to detect variations in sulfate levels of greater than 10 micron-g/cu m over distances of the order of 50 to 75 km; (4) repeat coverage at least on a daily basis; and (5) satellite observations during the mid to late morning local time, when the sulfate levels have begun to increase after the early morning minimum levels, and convective-type cloud cover has not yet increased to the amount reached later in the afternoon. Analysis of the satellite imagery shows that convective clouds can obscure haze patterns. Additional parameters based on spectral analysis include wavelength and bandwidth requirements.

Burke, H. H. K.↗

Quantitative Characterization of Hyper-Local Atmospheric Greenhouse Gas Sources

Atmospheric greenhouse gas (GHG) emissions are often characterized using stationary, tower-based sensors. Ground based sensors reside in the turbulent boundary layer and are subject to intense concentration impulses from hyper-local (<100m) point sources of emissions. These high frequency spikes are often filtered out in broader emission flux studies, losing valuable information about how hyper-local sources influence receptors. In this study, we investigated how empirical atmospheric data can be used to locate and quantify a concurrently measured hyper-local point source in a dense urban setting. An eddy covariance style tower and a low-cost sensor tower were deployed in various locations around an urban, hyper-local CO2/CH4 emissions source (a continuously measured restaurant exhaust vent). A model using different processing and statistical techniques was built to examine the most effective procedures for source isolation, directional location, and emission quantification. Using excess concentrations above a minimum baseline, we identify the source using bivariate polar plots and quantify the relationship between source size, receptor distance, and statistical proxies. Furthermore, we find that varying statistical thresholds allows for identification of less influential sources which are drowned out by larger or closer sources. Finally, we show that large sources can be effectively characterized using low-cost sensors, a valuable outcome informing how networks for monitoring larger areas could be implemented. This work may provide a basis for source identification and monitoring protocols for networks that feature sensors influenced by hyper-local point sources, subject to site-specific assumptions.

54 ENVIRONMENTAL SCIENCES↗

An evaluation of global organic aerosol schemes using airborne observations

Chemical transport models have historically struggled to accurately simulate the magnitude and variability of observed organic aerosol (OA), with previous studies demonstrating that models significantly underestimate observed concentrations in the troposphere. In this study, we explore two different model OA schemes within the standard GEOS-Chem chemical transport model and evaluate the simulations against a suite of 15 globally distributed airborne campaigns from 2008 to 2017, primarily in the spring and summer seasons. These include the ATom, KORUS-AQ, GoAmazon,FRAPPE, SEAC4RS, SENEX, DC3, CalNex, OP3, EUCAARI, ARCTAS and ARCPAC campaigns and provide broad coverage over a diverse set of atmospheric composition regimes – anthropogenic, biogenic, pyrogenic and remote. The schemes include significant differences in their treatment of the primary and secondary components of OA – a “simple scheme” that models primary OA (POA) as non-volatile and takes a fixed-yield approach to secondary OA (SOA) formation and a “complex scheme” that simulates POA as semi-volatile and uses a more sophisticated volatility basis set approach for non-isoprene SOA, with an explicit aqueous uptake mechanism to model isoprene SOA. Despite these substantial differences, both the simple and complex schemes perform comparably across the aggregate dataset in their ability to capture the observed variability (with an R2 of 0.41 and 0.44, respectively). The simple scheme displays greater skill in minimizing the overall model bias (with a normalized mean bias of 0.04 compared to 0.30 for the complex scheme). Across both schemes, the model skill in reproducing observed OA is superior to previous model evaluations and approaches the fidelity of the sulfate simulation within the GEOS-Chem model. However,there are significant differences in model performance across different chemical source regimes, classified here into seven categories.Higher-resolution nested regional simulations indicate that model resolution is an important factor in capturing variability in highly localized campaigns, while also demonstrating the importance of well-constrained emissions inventories and local meteorology, particularly over Asia. Our analysis suggests that a semi-volatile treatment of POA is superior to anon-volatile treatment. It is also likely that the complex scheme parameterization overestimates biogenic SOA at the global scale. While this study identifies factors within the SOA schemes that likely contribute to OA model bias (such as a strong dependency of the bias in the complex scheme on relative humidity and sulfate concentrations), comparisons with the skill of the sulfate aerosol scheme in GEOS-Chem indicate the importance of other drivers of bias, such as emissions, transport and deposition, that are exogenous to the OA chemical scheme.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Robust integration schemes for generalized viscoplasticity with internal-state variables. Part 1: Theoretical developments and applications

This two-part report is concerned with the development of a general framework for the implicit time-stepping integrators for the flow and evolution equations in generalized viscoplastic models. The primary goal is to present a complete theoretical formulation, and to address in detail the algorithmic and numerical analysis aspects involved in its finite element implementation, as well as to critically assess the numerical performance of the developed schemes in a comprehensive set of test cases. On the theoretical side, the general framework is developed on the basis of the unconditionally-stable, backward-Euler difference scheme as a starting point. Its mathematical structure is of sufficient generality to allow a unified treatment of different classes of viscoplastic models with internal variables. In particular, two specific models of this type, which are representative of the present start-of-art in metal viscoplasticity, are considered in applications reported here; i.e., fully associative (GVIPS) and non-associative (NAV) models. The matrix forms developed for both these models are directly applicable for both initially isotropic and anisotropic materials, in general (three-dimensional) situations as well as subspace applications (i.e., plane stress/strain, axisymmetric, generalized plane stress in shells). On the computational side, issues related to efficiency and robustness are emphasized in developing the (local) interative algorithm. In particular, closed-form expressions for residual vectors and (consistent) material tangent stiffness arrays are given explicitly for both GVIPS and NAV models, with their maximum sizes 'optimized' to depend only on the number of independent stress components (but independent of the number of viscoplastic internal state parameters). Significant robustness of the local iterative solution is provided by complementing the basic Newton-Raphson scheme with a line-search strategy for convergence. In the present first part of the report, we focus on the theoretical developments, and discussions of the results of numerical-performance studies using the integration schemes for GVIPS and NAV models.

Saleeb, Atef F.↗

Robust integration schemes for generalized viscoplasticity with internal-state variables. Part 2: Algorithmic developments and implementation

This two-part report is concerned with the development of a general framework for the implicit time-stepping integrators for the flow and evolution equations in generalized viscoplastic models. The primary goal is to present a complete theoretical formulation, and to address in detail the algorithmic and numerical analysis aspects involved in its finite element implementation, as well as to critically assess the numerical performance of the developed schemes in a comprehensive set of test cases. On the theoretical side, the general framework is developed on the basis of the unconditionally-stable, backward-Euler difference scheme as a starting point. Its mathematical structure is of sufficient generality to allow a unified treatment of different classes of viscoplastic models with internal variables. In particular, two specific models of this type, which are representative of the present start-of-art in metal viscoplasticity, are considered in applications reported here; i.e., fully associative (GVIPS) and non-associative (NAV) models. The matrix forms developed for both these models are directly applicable for both initially isotropic and anisotropic materials, in general (three-dimensional) situations as well as subspace applications (i.e., plane stress/strain, axisymmetric, generalized plane stress in shells). On the computational side, issues related to efficiency and robustness are emphasized in developing the (local) interative algorithm. In particular, closed-form expressions for residual vectors and (consistent) material tangent stiffness arrays are given explicitly for both GVIPS and NAV models, with their maximum sizes 'optimized' to depend only on the number of independent stress components (but independent of the number of viscoplastic internal state parameters). Significant robustness of the local iterative solution is provided by complementing the basic Newton-Raphson scheme with a line-search strategy for convergence. In the present second part of the report, we focus on the specific details of the numerical schemes, and associated computer algorithms, for the finite-element implementation of GVIPS and NAV models.

Li, Wei↗

Test of model of equivalence of tree height growth and transpiration rates in percolation-based phenomenology for root-soil interaction

We recently introduced a phenomenological model to represent root-soil interactions that accounts for variability in tree growth on local scales and within species (Hunt et al., 2020). Here, using publicly available datasets, we examine the developed model to provide additional empirical evidence, which also gives novel insight into two particular tree species, Eucalyptus regnans and Eucalyptus sieberri. In particular, we discovered two novel aspects of the model: a power-law time dependence with a different power-law exponent than usually heretofore applied, and apparently a direct equivalence of results for transpiration and growth rates. Also, although our original expectation was that growth and transpiration were merely proportional, comparison of the entire BAAD data set (Falster et al., 2015) for tree height with our model (Hunt et al., 2020; Fig. 1) already suggested approximate equality of these disparate measurements (volume/area, height). In this present manuscript we find further evidence for equivalence of the two measurements and provide additional theoretical basis for a model assumption regarding constancy of sap flow rates.

54 ENVIRONMENTAL SCIENCES↗

A new convective adjustment scheme

A marked divergence in philosophy is presented. The primary objective of the proposed parameterization scheme is to ensure that the local vertical temperature and moisture structures, which in nature are strongly constrained by convection, be realistic in the large scale model. The concept of a quasi-equilibrium between the cloud field and the large-scale forcing has been well established. This means that convective regions have characteristic temperature and moisture structures that can be documented observationally and used as the basis of a convective adjustment procedure. In the scheme proposed here, the temperature and moisture structures will be adjusted simultaneously toward observed quasi-equilibrium structures. This ensures that on the grid scale a global model always maintains a realistic vertical temperature and moisture structure in the presence of convection. The saturation point formulation of moist thermodynamics will be used to introduce the observational and theoretical basis of the proposed convective adjustment. The scheme is then applied to a series of data sets from GATE, BOMEX, ATEX and an Arctic air-mass transformation to show the sensitivity of the scheme to different parameters and develop a parameter set suitable for both shallow and deep convection in a global model.

Betts, A. K.↗

Limits to the lunar atmosphere

Apollo UV spectrometer experiment set limits on the density of oxygen of less than 500/cu cm, and the Apollo Lunar Atmospheric Composition Experiment data imply a value less than 50/cu cm above the subsolar point. These limits are surprisingly small relative to the measured value for sodium. A simple consideration of sources and sinks predicts significantly greater densities of oxygen. It is possible but doubtful that the Apollo measurements occurred during an epoch in which source rates were small. A preferential loss process for oxygen on the darkside of the moon is considered in which ionization by electron capture in surface collisions leads to escape through acceleration in the local electric field. Cold trapping in permanently shadowed regions as a net sink is considered and discounted, but the episodic nature of cometary insertion may allow formation of ice layers which act as a stabilized source of OH. On the basis of an assumed meteoroid impact source, a possible emission brightness of 50 R in the OH(A - X)(0,0) band above the lunar bright limb is predicted.

Morgan, T. H.↗

Global Stratigraphy of Venus: Analysis of a Random Sample of Thirty-Six Test Areas

The age relations between 36 impact craters with dark paraboloids and other geologic units and structures at these localities have been studied through photogeologic analysis of Magellan SAR images of the surface of Venus. Geologic settings in all 36 sites, about 1000 x 1000 km each, could be characterized using only 10 different terrain units and six types of structures. These units and structures form a major stratigraphic and geologic sequence (from oldest to youngest): (1) tessera terrain; (2) densely fractured terrains associated with coronae and in the form of remnants among plains; (3) fractured and ridged plains and ridge belts; (4) plains with wrinkle ridges; (5) ridges associated with coronae annulae and ridges of arachnoid annulae which are contemporary with wrinkle ridges of the ridged plains; (6) smooth and lobate plains; (7) fractures of coronae annulae, and fractures not related to coronae annulae, which disrupt ridged and smooth plains; (8) rift-associated fractures; and (9) craters with associated dark paraboloids, which represent the youngest 1O% of the Venus impact crater population (Campbell et al.), and are on top of all volcanic and tectonic units except the youngest episodes of rift-associated fracturing and volcanism; surficial streaks and patches are approximately contemporary with dark-paraboloid craters. Mapping of such units and structures in 36 randomly distributed large regions (each approximately 10(exp 6) sq km) shows evidence for a distinctive regional and global stratigraphic and geologic sequence. On the basis of this sequence we have developed a model that illustrates several major themes in the history of Venus. Most of the history of Venus (that of its first 80% or so) is not preserved in the surface geomorphological record. The major deformation associated with tessera formation in the period sometime between 0.5-1.0 b.y. ago (Ivanov and Basilevsky) is the earliest event detected. In the terminal stages of tessera fon-nation, extensive parallel linear graben swarms representing a change in the style of deformation from shortening to extension were formed on the tessera and on some volcanic plains that were emplaced just after, and perhaps also during the latter stages of the major compressional phase of tessera emplacement. Our stratigraphic analyses suggest that following tessera formation, extensive volcanic flooding resurfaced at least 85% of the planet in the form of the presently-ridged and fractured plains. Several lines of evidence favor a high flux in the post-tessera period but we have no independent evidence for the absolute duration of ridged plains emplacement. During this time, the net state of stress in the lithosphere apparently changed from extensional to compressional, first in the form of extensive ridge belt development, followed by the formation of extensive wrinkle ridges on the flow units. Subsequently, there occurred local emplacement of smooth and lobate plains units which are presently essentially undefortned. The major events in the latest 10% of the presently preserved history of Venus (less than 50 m.y. ago) are continued rifting and some associated volcanism, and the redistribution of eolian material largely derived from impact crater deposits.

Basilevsky, Alexander T.↗

The force on the flex: Global parallelism and portability

A parallel programming methodology, called the force, supports the construction of programs to be executed in parallel by an unspecified, but potentially large, number of processes. The methodology was originally developed on a pipelined, shared memory multiprocessor, the Denelcor HEP, and embodies the primitive operations of the force in a set of macros which expand into multiprocessor Fortran code. A small set of primitives is sufficient to write large parallel programs, and the system has been used to produce 10,000 line programs in computational fluid dynamics. The level of complexity of the force primitives is intermediate. It is high enough to mask detailed architectural differences between multiprocessors but low enough to give the user control over performance. The system is being ported to a medium scale multiprocessor, the Flex/32, which is a 20 processor system with a mixture of shared and local memory. Memory organization and the type of processor synchronization supported by the hardware on the two machines lead to some differences in efficient implementations of the force primitives, but the user interface remains the same. An initial implementation was done by retargeting the macros to Flexible Computer Corporation's ConCurrent C language. Subsequently, the macros were caused to directly produce the system calls which form the basis for ConCurrent C. The implementation of the Fortran based system is in step with Flexible Computer Corporations's implementation of a Fortran system in the parallel environment.

Jordan, H. F.↗

Assessment of the Quality of the Version 1.07 Temperature-Versus-Pressure Profiles of the Middle Atmosphere from TIMED/SABER

The quality of the retrieved temperature-versus-pressure (or T(p)) profiles is described for the middle atmosphere for the publicly available Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) Version 1.07 (V1.07) data set. The primary sources of systematic error for the SABER results below about 70 km are (1) errors in the measured radiances, (2) biases in the forward model, and (3) uncertainties in the corrections for ozone and in the determination of the reference pressure for the retrieved profiles. Comparisons with other correlative data sets indicate that SABER T(p) is too high by 1-3 K in the lower stratosphere but then too low by 1 K near the stratopause and by 2 K in the middle mesosphere. There is little difference between the local thermodynamic equilibrium (LTE) algorithm results below about 70 km from V1.07 and V1.06, but there are substantial improvements/differences for the non-LTE results of V1.07 for the upper mesosphere and lower thermosphere (UMLT) region. In particular, the V1.07 algorithm uses monthly, diurnally averaged CO2 profiles versus latitude from the Whole Atmosphere Community Climate Model. This change has improved the consistency of the character of the tides in its kinetic temperature (T(sub k)). The T(sub k) profiles agree with UMLT values obtained from ground-based measurements of column-averaged OH and O2 emissions and of the Na lidar returns, at least within their mutual uncertainties. SABER T(sub k) values obtained near the mesopause with its daytime algorithm also agree well with the falling sphere climatology at high northern latitudes in summer. It is concluded that the SABER data set can be the basis for improved, diurnal-to-interannual-scale temperatures for the middle atmosphere and especially for its UMLT region.

Remsberg, E. E.↗

Analysis of TIMS performance subjected to simulated wind blast

The results of the performance of the Thermal Infrared Multispectral Scanner (TIMS) when it is subjected to various wind conditions in the laboratory are described. Various wind conditions were simulated using a 24 inch fan or combinations of air jet streams blowing toward either or both of the blackbody surfaces. The fan was used to simulate a large volume of air flow at moderate speeds (up to 30 mph). The small diameter air jets were used to probe TIMS system response in reaction to localized wind perturbations. The maximum nozzle speed of the air jet was 60 mph. A range of wind directions and speeds were set up in the laboratory during the test. The majority of the wind tests were conducted under ambient conditions with the room temperature fluctuating no more than 2 C. The temperature of the high speed air jet was determined to be within 1 C of the room temperature. TIMS response was recorded on analog tape. Additional thermistor readouts of the blackbody temperatures and thermocouple readout of the ambient temperature were recorded manually to be compared with the housekeeping data recorded on the tape. Additional tests were conducted under conditions of elevated and cooled room temperatures. The room temperature was varied between 19.5 to 25.5 C in these tests. The calibration parameters needed for quantitative analysis of TIMS data were first plotted on a scanline-by-scanline basis. These parameters are the low and high blackbody temperature readings as recorded by the TIMS and their corresponding digitized count values. Using these values, the system transfer equations were calculated. This equation allows us to compute the flux for any video count by computing the slope and intercept of the straight line that relates the flux to the digital count. The actual video of the target (the lab floor in this case) was then compared with a simulated target. This simulated target was assumed to be a blackbody at emissivity of .95 degrees and the temperature was assumed to be at ambient temperature as recorded by the TIMS for each scanline. Using the slope and the intercept the flux corresponding to this target was converted into digital counts. The counts were observed to have a strong correlation with the actual video as recorded by the TIMS. The attached graphs describe the performance of the TIMS when compressed air is blown at each one of the blackbodies at different speeds. The effect of blowing a fan and changing the room temperature is also being analyzed. Results indicate that the TIMS system responds to variation in wind speed in real time and maintains the capability to produce accurate temperatures on a scan line basis.

Jaggi, S.↗