Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “collective variables”

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 181 records · Page 10

Algorithm for Detection of Ground and Canopy Cover in Micropulse Photon-Counting Lidar Altimeter Data in Preparation for the ICESat-2 Mission

NASA's Ice, Cloud and Land Elevation Satellite-II (ICESat-2) mission is a decadal survey mission (2016 launch). The mission objectives are to measure land ice elevation, sea ice freeboard, and changes in these variables, as well as to collect measurements over vegetation to facilitate canopy height determination. Two innovative components will characterize the ICESat-2 lidar: 1) collection of elevation data by a multibeam system and 2) application of micropulse lidar (photon-counting) technology. A photon-counting altimeter yields clouds of discrete points, resulting from returns of individual photons, and hence new data analysis techniques are required for elevation determination and association of the returned points to reflectors of interest. The objective of this paper is to derive an algorithm that allows detection of ground under dense canopy and identification of ground and canopy levels in simulated ICESat-2 data, based on airborne observations with a Sigma Space micropulse lidar. The mathematical algorithm uses spatial statistical and discrete mathematical concepts, including radial basis functions, density measures, geometrical anisotropy, eigenvectors, and geostatistical classification parameters and hyperparameters. Validation shows that ground and canopy elevation, and hence canopy height, can be expected to be observable with high accuracy by ICESat-2 for all expected beam energies considered for instrument design (93.01%-99.57% correctly selected points for a beam with expected return of 0.93 mean signals per shot (msp), and 72.85%-98.68% for 0.48 msp). The algorithm derived here is generally applicable for elevation determination from photoncounting lidar altimeter data collected over forested areas, land ice, sea ice, and land surfaces, as well as for cloud detection.

Algorithms↗

The Little Ice Age and Human-Environmental Interactions in the Central Balkans: Insights from a New Serbian Paleorecord

This paper presents a 600-year well-dated, high-resolution Central Balkan paleo-record including the Little Ice Age (LIA; 1450-850 CE). Utilizing pollen-based REVEALS modelling estimates, geochemical indicators, rarefaction analyses and the AMS 14C-based Bacon age model, this first-hand record from the Sava Basin reveals the transformation of the Central Balkan landscape involving linkages between changing climatic and socio-political regimes. The pre-LIA interval (1370-1418 CE) in the Sava Region reveals a wooded steppe and increased cultivation under warmer/stable climatic and socio-political conditions. In contrast, the LIA interval in the region is expressed through continuous transitions between forest and grassland, extensive land erosion and stressed agriculture, potentially as a collective artifact of the climatic variability and human impact associated with socio- political stressors of the time. The post-LIA/Industrial Era interval (1850-2012 CE) in the Sava Region shows an overall increase in woodland as well as agriculture following the exit of the Ottomans. However, increased population pressures and the subsequent onset of the Industrial Revolution as well as increased trade led to intense deforestation throughout the 20th Century, which continued during the Socialist period.

Societal impacts of climate change↗

Use of an Alternative Conceptual Model of Vadose Zone Heterogeneity to Evaluate Past Tank Leaks and Other Unplanned Releases within a Tank Farm at the Hanford Site - 20061

Washington State Department of Ecology (Ecology) requested that U.S. Department of Energy (DOE)-Office of River Protection (ORP) consider an evaluation of effects of fine-grained thin sediment layers on transport with a separate alternative conceptual model in its evaluation of the potential impact of tank leaks and other unplanned releases within the tank farm at the Hanford Site. For this alternative model, Ecology recommended that the model be developed based on the general framework of fine-grained units identified by a stakeholder group in their interpretations of variability in moisture content data collected in the vicinity of Waste Management Area (WMA) C. In discussions with DOE-ORP, Ecology acknowledged that the underlying data and interpretations of the occurrence and lateral continuity of the fine-grained thin layers identified by a stakeholder group are uncertain. However, DOE-ORP agreed to Ecology's recommendation and has provided support for the requested evaluation that involved development of an alternative model based on the general framework of the unpublished report by the stakeholder. DOE-ORP considered this evaluation to be a hypothetical evaluation of vadose zone heterogeneities at WMA C. General observations from the range of simulation cases examined in the evaluation of the effects of hypothetical vadose zone heterogeneities at WMA C are as follows. - The movement of the center of mass of the simulated plumes was generally vertically downward below the source for all simulations, including those that incorporated the hypothetical fine-grained units. - All simulations that incorporated hypothetical heterogeneity produced additional plume spreading over what was produced in simulations using Equivalent Homogeneous Media (EHM) model(s)a. The spreading resulted in a broadening of the fringes of the plume, resulting in a wider region of low concentration, but lower peak concentrations associated with the center of mass of the plume. - Simulations that used the silty-sand hydraulic properties, suggested by Ecology for the hypothetical fine-grained units, generally produced similar spreading and slightly lower peak mass flux at the water table when compared to the EHM modeling results. a The EHM-based models used at WMA C do not explicitly include small-scale fine-grained heterogeneities used in the model advocated by Ecology. - Simulations that used the silty-sand hydraulic properties, suggested by Ecology for the hypothetical fine-grained units, generally produced less spreading and an earlier arrival of mass flux at the water table when compared to the use of another set of hydraulic properties from a silty-sand sample collected at a nearby disposal facility. - The EHM representation of the vadose zone generally produced higher peak mass flux and an earlier occurrence of peak fluxes at the water table compared to all analyses incorporating additional hypothetical heterogeneity. Results of this alternative model evaluation of past leaks provided some insight into the transport effects of vadose zone heterogeneities at WMA C that helped resolve Ecology's comments and issues related to the effects of vadose zone heterogeneities on past leaks and losses from the WMA C tank farms area. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Estimation of soil classes and their relationship to grapevine vigor in a Bordeaux vineyard: advancing the practical joint use of electromagnetic induction (EMI) and NDVI datasets for precision viticulture

Working within a vineyard in the Pessac Léognan Appellation of Bordeaux, France, this study documents the potential of using simple statistical methods with spatially-resolved and increasingly available electromagnetic induction (EMI) geophysical and normalized difference vegetation index (NDVI) datasets to accurately estimate Bordeaux vineyard soil classes and to quantitatively explore the relationship between vineyard soil types and grapevine vigor. First, co-located electrical tomographic tomography (ERT) and EMI datasets were compared to gain confidence about how the EMI method averaged soil properties over the grapevine rooting depth. Then, EMI data were used with core soil texture and soil-pit based interpretations of Bordeaux soil types (Brunisol, Redoxisol, Colluviosol and Calcosol) to estimate the spatial distribution of geophysically-identified Bordeaux soil classes. A strong relationship (r = 0.75, p < 0.01) was revealed between the geophysically-identified Bordeaux soil classes and NDVI (both 2 m resolution), showing that the highest grapevine vigor was associated with the Bordeaux soil classes having the largest clay fraction. The results suggest that within-block variability of grapevine vigor was largely controlled by variability in soil classes, and that carefully collected EMI and NDVI datasets can be exceedingly helpful for providing quantitative estimates of vineyard soil and vigor variability, as well as their covariation. The method is expected to be transferable to other viticultural regions, providing an approach to use easy-to-acquire, high resolution datasets to guide viticultural practices, including routine management and replanting.

54 ENVIRONMENTAL SCIENCES↗

Temporal relationships between spectral response and agronomic variables of a corn canopy

Attention is given to an experiment in which spectral radiance data collected in three spectral regions are related to corn canopy variables. The study extends the work of Tucker et al. (1979) in that more detailed measurements of corn canopy variables were made using quantitative techniques. Wet and dry green leaf biomass is considered along with the green leaf area index, chlorotic leaf biomass, chlorotic leaf area, and leaf water content. In addition, spectral data were collected with a hand-held radiometer having Landsat-D Thematic Mapper (TM) bands TM3 (0.63-0.69 micrometers), TM4 (0.76-0.90 micrometers), and TM5 (1.55-1.75 micrometers). TM3, TM4, and TM5 seem to be well situated spectrally for making remotely sensed measurements related to chlorophyll concentration, leaf density, and leaf water content.

Kimes, D. S.↗

Characterization of Tank 11H Samples from Tank Closure Cesium Removal 1A (TCCR 1A) Batch 1

Savannah River Mission Completion (SRMC) is currently operating the Tank Closure Cesium Removal 1A (TCCR 1A) process to remove 137 Cs from tank waste supernate using an ion exchange process. The TCCR 1A unit processes dissolved salt fed from Tank 10H through a series of ion exchange columns containing crystalline silicotitanate (CST, IONSIV TM R9120-B, 30x60) and the effluent is then discharged to Tank 11H. In support of the TCCR 1A program, SRNL analyzed samples taken from Tank 11H (without tank mixing) before, during, and at the completion of TCCR 1A Batch 1 processing. Tank 11H serves as the receipt tank for the filtered and cesium removed product from the TCCR 1A system. Processing of Batch 1 commenced on January 13, 2022 and completed on February 17, 2022, after processing approximately 70,100 gallons. A pre-production sample was collected from the heel remaining in Tank 11H just before Batch 1 processing began. In addition, five interim surface samples were collected from Tank 11H during processing, and both a surface and a variable depth sample (VDS, ~7” from tank bottom) were collected just after processing completed. Analysis of all samples included density and gamma spectroscopy, in addition to a more comprehensive suite of analytes for the pre- and post-production samples. The density of the pre-production sample was the highest of all samples measured and was then observed to decrease for the first two interim samples, followed by becoming fairly consistent for the remainder of the samples (~1.3 g/mL). The density of 1.3 g/mL is similar to the density measured for one of the three Tank 10H qualification samples (HTF-10-21-126). The 137 Cs activity was found to decrease as additional decontaminated effluent from the TCCR columns was added to Tank 11H during processing; however, an increase in activity was observed during periods of no processing which can be attributed to leaching of 137 Cs from the known solids in Tank 11H. The Cs isotope concentrations in the Tank 11H postproduction surface sample were determined to be 99.7-99.8% lower than the concentrations measured in the Tank 10H feed as measured by mass spectrometry.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

JPL field measurements at the Finney County, Kansas, test site, October 1976: Meteorological variables, surface reflectivity, surface and subsurface temperatures

Data collected at the Finney County, Kansas test site as part of the Joint Soil Moisture Experiment (JSME) are presented here, prior to analysis, to provide all JSME investigators with an immediate source of primary information. The ground-truth measurements were taken to verify and complement soil moisture data taken by microwave and infrared sensors during aircraft overflights. Measurements were made of meteorological variables (air speed, temperature, relative humidity, and rainfall), surface reflectivity, and temperatures at and below the surface.

Kahle, A. B.↗

Computational Imaging for Intelligence in Highly Scattering Aerosols (Final Report)

Natural and man-made degraded visual environments pose major threats to national security. The random scattering and absorption of light by tiny particles suspended in the air reduces situational awareness and causes unacceptable down-time for critical systems and operations. To improve the situation, we have developed several approaches to interpret the information contained within scattered light to enhance sensing and imaging in scattering media. These approaches were tested at the Sandia National Laboratory Fog Chamber facility and with tabletop fog chambers. Computationally efficient light transport models were developed and leveraged for computational sensing. The models are based on a weak angular dependence approximation to the Boltzmann or radiative transfer equation that appears to be applicable in both the moderate and highly scattering regimes. After the new model was experimentally validated, statistical approaches for detection, localization, and imaging of objects hidden in fog were developed and demonstrated. A binary hypothesis test and the Neyman-Pearson lemma provided the highest theoretically possible probability of detection for a specified false alarm rate and signal-to-noise ratio. Maximum likelihood estimation allowed estimation of the fog optical properties as well as the position, size, and reflection coefficient of an object in fog. A computational dehazing approach was implemented to reduce the effects of scatter on images, making object features more readily discernible. We have developed, characterized, and deployed a new Tabletop Fog Chamber capable of repeatably generating multiple unique fog-analogues for optical testing in degraded visual environments. We characterized this chamber using both optical and microphysical techniques. In doing so we have explored the ability of droplet nucleation theory to describe the aerosols generated within the chamber, as well as Mie scattering theory to describe the attenuation of light by said aerosols, and correlated the aerosol microphysics to optical properties such as transmission and meteorological optical range (MOR). This chamber has proved highly valuable and has supported multiple efforts inclusive to and exclusive of this LDRD project to test optics in degraded visual environments. Circularly polarized light has been found to maintain its polarization state better than linearly polarized light when propagating through fog. This was demonstrated experimentally in both the visible and short-wave infrared (SWIR) by imaging targets made of different commercially available retroreflective films. It was found that active circularly polarized imaging can increase contrast and range compared to linearly polarized imaging. We have completed an initial investigation of the capability for machine learning methods to reduce the effects of light scattering when imaging through fog. Previously acquired experimental long-wave images were used to train an autoencoder denoising architecture. Overfitting was found to be a problem because of lack of variability in the object type in this data set. The lessons learned were used to collect a well labeled dataset with much more variability using the Tabletop Fog Chamber that will be available for future studies. We have developed several new sensing methods using speckle intensity correlations. First, the ability to image moving objects in fog was shown, establishing that our unique speckle imaging method can be implemented in dynamic scattering media. Second, the speckle decorrelation over time was found to be sensitive to fog composition, implying extensions to fog characterization. Third, the ability to distinguish macroscopically identical objects on a far-subwavelength scale was demonstrated, suggesting numerous applications ranging from nanoscale defect detection to security. Fourth, we have shown the capability to simultaneously image and localize hidden objects, allowing the speckle imaging method to be effective without prior object positional information. Finally, an interferometric effect was presented that illustrates a new approach for analyzing speckle intensity correlations that may lead to more effective ways to localize and image moving objects. All of these results represent significant developments that challenge the limits of the application of speckle imaging and open important application spaces. A theory was developed and simulations were performed to assess the potential transverse resolution benefit of relative motion in structured illumination for radar systems. Results for a simplified radar system model indicate that significant resolution benefits are possible using data from scanning a structured beam over the target, with the use of appropriate signal processing.

58 GEOSCIENCES↗

XANES spectroscopy of sulfides stable under reducing conditions

X-ray absorption near-edge structure (XANES) spectroscopy is a powerful technique to quantitatively investigate sulfur speciation in geologically complex materials such as minerals, glasses, soils, organic compounds, industrial slags, and extraterrestrial materials. This technique allows non-destructive investigation of the coordination chemistry and oxidation state of sulfur species ranging from sulfide (2–oxidation state) to sulfate (6+ oxidation state). Each sulfur species has a unique spectral shape with a characteristic K-edge representing the s → p and d hybridization photoelectron transitions. As such, sulfur speciation is used to measure the oxidation state of samples by comparing the overall XANES spectra to that of reference compounds. Although many S XANES spectral standards exist for terrestrial applications under oxidized conditions, new sulfide standards are needed to investigate reduced (oxygen fugacity, f O 2 , below IW) silicate systems relevant for studies of extraterrestrial materials and systems. Sulfides found in certain meteorites (e.g., enstatite chondrites and aubrites) and predicted to exist on Mercury, such as CaS (oldhamite), MgS (niningerite), and FeCr 2 S 4 (daubréelite), are stable at f O 2 below IW-3 but rapidly oxidize to sulfate and/or produce sulfurous gases under terrestrial surface conditions. XANES spectra of these compounds collected to date have been of variable quality, possibly due to the unstable nature of certain sulfides under typical (e.g., oxidizing) laboratory conditions. A new set of compounds was prepared for this study and their XANES spectra are analyzed for comparison with potential extraterrestrial analogs. S K -edge XANES spectra were collected at Argonne National Lab for FeS (troilite), MnS (alabandite), CaS (oldhamite), MgS (niningerite), Ni 1–x S, NiS 2 , CaSO 4 (anhydrite), MgSO 4 , FeSO 4 , Fe 2 (SO 4 ) 3 , FeCr 2 S 4 (daubréelite), Na 2 S, Al 2 S 3 , Ni 7 S 6 , and Ni 3 S 2 ; the latter five were analyzed for the first time using XANES. These standards expand upon the existing S XANES end-member libraries at a higher spectral resolution (0.25 eV steps) near the S K -edge. Processed spectra, those that have been normalized and “flattened,” are compared to quantify uncertainties due to data processing methods. Future investigations that require well-characterized sulfide standards, such as the ones presented here, may have important implications for understanding sulfur speciation in reduced silicate glasses and minerals with applications for the early Earth, Moon, Mercury, and enstatite chondrites.

Anzures, Brendan A.↗

A new class of variable-radii diffraction optics for high-resolution x-ray spectroscopy at the National Ignition Facility (invited)

A new class of crystal shapes has been developed for x-ray spectroscopy of point-like or small (a few mm) emission sources. These optics allow for dramatic improvement in both achievable energy resolution and total throughput of the spectrometer as compared with traditional designs. This class of crystal shapes, collectively referred to as the Variable-Radii Spiral (VR-Spiral), utilize crystal shapes in which both the major and minor radii are variable. A crystal using this novel VR-Spiral shape has now been fabricated for high-resolution Extended X-ray Absorption Fine Structure (EXAFS) experiments targeting the Pb-L 3 (13.0 keV) absorption edge at the National Ignition Facility. The performance of this crystal has been characterized in the laboratory using a microfocus x-ray source, showing that high-resolution high-throughput EXAFS spectra can be acquired using this geometry. Importantly, these successful tests show that the complex three-dimensional crystal shape is manufacturable with the required precision needed to realize the expected performance of better than 5 eV energy resolution while using a 30 mm high crystal. An improved generalized mathematical form for VR-Spiral shapes is also presented allowing improved optimization as compared to the first sinusoidal-spiral based design. Finally, this new formulation allows VR-Spiral spectrometers to be designed at any magnification with optimized energy resolution at all energies within the spectrometer bandwidth.

47 OTHER INSTRUMENTATION↗

Pressure Transducer Measurement Variability in Deep Wells Screened Across the Water Table

Abstract Automated water level measurements collected using vented pressure transducers in deep wells screened across the water table may exhibit a greater response to barometric pressure changes than the true water level. The cause was hypothesized to be disequilibrium in barometric pressure between the wellbores and land surface due to air exchange with the deep vadose zone. In this study, vented and nonvented pressure transducers were installed and operated simultaneously in two deep wells screened across the water table. A vent tube open to the atmosphere at land surface allowed for barometric compensation of the vented transducers. Two nonvented transducers were installed in each well, one submerged in the water and one above the water surface. The difference in readings allowed for barometric compensation. Manual measurements were also collected. It was confirmed that measurements from the vented transducers exhibited greater variability in response to barometric pressure changes than the nonvented transducers and manual measurements. Comparison of the downhole barometric pressure measurements to values from a nearby meteorology station showed the response in the wells to changes in barometric pressure was time‐lagged and attenuated. Thus, the reference pressure from land surface supplied to the vented transducers was not representative of the air pressure within the wells. This caused fluctuations of the transducer readings in response to barometric pressure changes to be greater than the true water level change. This issue can be resolved by the use of nonvented pressure transducers.

McDonald, John P.↗

Understanding the structural and morphological effects of synthesis route on NpO 2

The availability of actinide standard materials for use in nuclear safeguard applications is critical, as is thorough characterization thereof. Although accurate trace element compositions and isotopic considerations are paramount for deployment of reference standards, structural characterization is also essential towards accurately describing the chemical form and potential matrix effects in candidate materials. Here, to this end, samples of NpO 2 were synthesized via a direct denitration (DD) method and probed with powder X-ray diffraction (PXRD), Raman spectroscopy, and scanning electron microscopy (SEM) for structural and morphological characterization and comparison with NpO 2 materials produced via modified direct denitration (MDD). PXRD confirmed the bulk identity of NpO 2 , and no additional phases were identified using this method. Analysis of Raman data collected using a 532 nm excitation wavelength indicates that samples are mostly phase pure; however, some variability in spectral features is observed. Analysis of additional spectroscopic data collected with a 785 nm excitation wavelength revealed variability in the relative intensity of spectral features. Raman spectroscopy indicates that the sample is primarily NpO 2 ; however, additional signals indicate possible structural disorder, oxidized species, or potential contributions from other Np phases. To further investigate the possibility of additional phase contributions within the sample of NpO 2 , Raman spectroscopic mapping was employed to examine the homogeneity of the sample produced via DD. From this analysis, we determined that despite variability in the intensity of Raman-active vibrational modes, consistent spectra are obtained throughout the area of the sample investigated. SEM images show aggregates with variable sizes and shapes, with rounded, primary particles possessing an average diameter of approximately 100 nm. Comparison of the results of these multimodal analyses to the literature indicates that the crystal chemical, spectroscopic, and microstructural properties of NpO 2 vary based on synthesis method, even if X-ray diffraction data indicate that the bulk phase is NpO 2 .

Direct denitration↗

Experimental Droplet Impingement on Several Two-Dimensional Airfoils with Thickness Ratios of 6 to 16 Percent

The rate and area of cloud droplet impingement on several two-dimensional swept and unswept airfoils were obtained experimentally in the NACA Lewis icing tunnel with a dye-tracer technique. Airfoil thickness ratios of 6 to 16 percent; angles of attack from 0 deg to 12 deg, and chord sizes from 13 to 96 inches were included in the study. The data were obtained at 152 knots and are extended to other conditions by dimensionless impingement parameters. In general, the data show that the total and local collection efficiencies and impingement limits are primary functions of the modified inertia parameter (in which airspeed, droplet size, and body size are the most significant variables) and the airfoil thickness ratio. Local collection efficiencies and impingement limits also depend on angle of attack. Secondary factors affecting impingement characteristics are airfoil shape, camber, and sweep angle. The impingement characteristics obtained experimentally for the airfoils were within +/-10 percent on the average of the characteristics calculated from theoretical trajectories. Over the range of conditions studied, the experimental data demonstrate that a specific method can be used to predict the impingement characteristics of swept airfoils with large aspect ratios from the data for unswept airfoils of the same series.

Gelder, Thomas F.↗

SULI Technical Report: WIRE (Draft)

Local climate undoubtedly affects facility maintenance requirements; however, knowledge of the link between weather and facility condition over time is grounded in the 1990s and needs to be updated. WIRE aims to quantify the relationship between climate and infrastructure by calculating climate-based service life estimates to be used in BUILDER, a program the Department of Energy (DOE) uses for facility management. Phase One of this project is a case study on window-or-wall-mounted AC units. The independent variables influencing service life are load, maintenance quality, and effective age, with load calculated from three weather variables: temperature, humidity, and solar radiation. Given sufficient data, service life can be calculated as a function of these variables via regression. WIRE is currently collecting data. Preliminary estimated service life projections as a function of load derived from Lawrence Livermore National Laboratory (LLNL) data support the need for updated analysis on the relationship between climate and facility condition. With the groundwork laid, WIRE can generate deliverables immediately upon arrival of data. The methods used for this case study constitute a framework for deriving the relationship between weather variables and service life for other components.

42 ENGINEERING↗

Single-Use Destructive Assay for Uranium Hexafluoride Sampling

Sampling uranium hexafluoride (UF6) for the determination of enrichments by destructive analysis (DA) is a critical component in the International Atomic Energy Agency’s layered safeguards approach for uranium processing facilities. Typically, gram-quantity UF6 samples are collected during inspections and stored under tag-and-seal until transportation to an off-site analytical laboratory. The shipping times can be long, and evolving restrictions on radioactive/corrosive materials shipments may increasingly limit the IAEA’s ability to transport UF6 samples easily. Pacific Northwest National Laboratory has developed a low-cost UF6 sampling technology called Single-Use Destructive Assay (SUDA) that addresses these challenges, as well as provides DA sample geometries that can be tailored for different analytical methods, including potential on-site analyses. The SUDA samplers, along with a unique holder, are designed for direct attachment to existing taps at uranium processing facilities, allowing gaseous UF6 to come into direct contact with a zeolite film. The SUDA technology features the ability to capture uranium in a more easily shipped and handled form as the solid, more stable, and relatively less hazardous hydrated uranyl fluoride (UO2F2•nH2O), which is formed through the controlled hydrolysis of UF6. We have recently simulated uranium collection under enrichment plant sampling conditions to further improve our understanding of SUDA sampling. Presented here is our recent work on measuring the relationship between sampling conditions and uranium collection, which includes control of the uranium-mass-to-zeolite ratio and assessing variable UF6 gas and sampling parameters that can affect collection using the SUDA sampler.

Pope, Timothy R.↗

BOREAS TF-10 NSA-YJP Tower Flux, Meteorological, and Porometry Data

The BOREAS TF-10 team collected tower flux and meteorological data at two sites, a fen and a young jack pine forest, near Thompson, Manitoba, Canada, as part of BOREAS. A preliminary data set was assembled in August 1993 while field testing the instrument packages, and at both sites data were collected from 15-Aug to 31-Aug. The main experimental period was in 1994, when continuous data were collected from the young jack pine site from 23-May to 20-Sep. Upon examination of the 1994 data set, it became clear that the behavior of the heat, water, and carbon dioxide fluxes throughout the whole growing season was an important scientific question, and that the 1994 data record was not sufficiently long to capture the character of the seasonal behavior of the fluxes. Thus, the young jack pine site was operated from 08-May to 07-Nov in 1996 in order to collect data from spring melt to autumn freeze-up. All variables are presented as 30-minute averages. Supporting data were also collected to describe the surface#s state and to provide the information, in association with the flux data, to build SVAT models. For the young jack pine site, these supporting data included stomatal conductance measurements. The data are stored in tabular ASCII files.

McCaughey, J. Harry↗

Solving Ordinary Differential Equations

Initial-value ordinary differential equation solution via variable order Adams method (SIVA/DIVA) package is collection of subroutines for solution of nonstiff ordinary differential equations. There are versions for single-precision and double-precision arithmetic. Requires fewer evaluations of derivatives than other variable-order Adams predictor/ corrector methods. Option for direct integration of second-order equations makes integration of trajectory problems significantly more efficient. Written in FORTRAN 77.

Krogh, F. T.↗

Two-Phase Thermal Switching System for a Small, Extended Duration Lunar Surface Science Platform

This paper describes a novel thermal control system for the Warm Electronics Box (WEB) on board a small lunar surface lander intended to support science activities anywhere on the lunar surface for an extended duration of up to 6 years. Virtually all lander electronics, which collectively dissipate about 60 W in the reference mission, are contained within the WEB. These devices must be maintained below 323 K (with a goal of 303 K) during the nearly 15-earth-day lunar day, when surface temperatures can reach 390K, and above 263 K during the nearly 15-earth-day lunar night, when surface temperatures can reach 100K. Because of the large temperature swing from lunar day-to-night, a novel thermal switching system was required that would be able to provide high conductance from WEB to radiator(s) during the hot lunar day and low (or negligible) conductance during the cold lunar night. The concept that was developed consists of ammonia variable conductance heat pipes (VCHPs) to collect heat from WEB components and a polymer wick propylene loop heat pipe (LHP) to transport the collected heat to the radiator(s). The VCHPs autonomously maximize transport when the WEB is warm and autonomously shut down when the WEB gets cold. The LHP autonomously shuts down when the VCHPs shut down. When the environment transitions from lunar night to day, the VCHPs and LHP autonomously turn back on. Out of 26 analyzed systems, this novel arrangement was able to best achieve the combined goals of zero control power, autonomous operation, long life, low complexity, low T, and landed tilt tolerance.

Bugby, David C.↗