This study demonstrates that containment of municipal and hazardous waste in arid and semiarid environments can be accomplished effectively without traditional, synthetic materials and complex, multi-layer systems. This research demonstrates that closure covers combining layers of natural soil, native plant species, and climatic conditions to form a sustainable, functioning ecosystem will meet the technical equivalency criteria prescribed by the U. S. Environmental Protection Agency. In this study, percolation through a natural analogue and an engineered cover is simulated using the one-dimensional, numerical code UNSAT-H. UNSAT-H is a Richards. equation-based model that simulates soil water infiltration, unsaturated flow, redistribution, evaporation, plant transpiration, and deep percolation. This study incorporates conservative, site-specific soil hydraulic and vegetation parameters. Historical meteorological data are used to simulate percolation through the natural analogue and an engineered cover, with and without vegetation. This study indicates that a 3-foot (ft) cover in arid and semiarid environments is the minimum design thickness necessary to meet the U. S. Environmental Protection Agency-prescribed technical equivalency criteria of 31.5 millimeters/year and 1 x 10{sup -7} centimeters/second for net annual percolation and average flux, respectively. Increasing cover thickness to 4 or 5 ft results in limited additional improvement in cover performance.
Effective, high-performance, networked file systems and storage is needed to solve I/O bottlenecks between large compute platforms. Frequently, parallel techniques such as PFTP, are employed to overcome the adverse effect of TCP's congestion avoidance algorithm in order to achieve reasonable aggregate throughput. These techniques can suffer from end-system bottlenecks due to the protocol processing overhead and memory copies involved in moving large amounts of data during I/O. Moreover, transferring data using PFTP requires manual operation, lacking the transparency to allow for interactive visualization and computational steering of large-scale simulations from distributed locations. This paper evaluates the emerging Internet SCSI (iSCSI) protocol [2] as the file/data transport in order that remote clients can transparently access data through a distributed global file system available to local clients. We started our work characterizing the performance behavior of iSCSI in Local Area Networks (LANs). We then proceeded to study the effect of propagation delay on throughput using remote iSCSI storage and explored optimization techniques to mitigate the adverse effects of long delay in high-bandwidth Wide Area Networks (WANs). Lastly, we evaluated iSCSI in a Storage Area Network (SAN) for a Global Parallel Filesystem. We conducted our benchmark based on typical usage model of large-scale scientific applications at Sandia. We demonstrated the benefit of high-performance parallel VO to scientific applications at the IEEE 2004 Supercomputing Conference, using experiences and knowledge gained from this study.
Techniques for mitigating the adsorption of {sup 137}Cs and {sup 60}Co on metal surfaces (e.g. RAM packages) exposed to contaminated water (e.g. spent-fuel pools) have been developed and experimentally verified. The techniques are also effective in removing some of the {sup 60}Co and {sup 137}Cs that may have been adsorbed on the surfaces after removal from the contaminated water. The principle for the {sup 137}Cs mitigation technique is based upon ion-exchange processes. In contrast, {sup 60}Co contamination primarily resides in minute particles of crud that become lodged on cask surfaces. Crud is an insoluble Fe-Ni-Cr oxide that forms colloidal-sized particles as reactor cooling systems corrode. Because of the similarity between Ni{sup 2+} and Co{sup 2+}, crud is able to scavenge and retain traces of cobalt as it forms. A number of organic compounds have a great specificity for combining with nickel and cobalt. Ongoing research is investigating the effectiveness of chemical complexing agent EDTA with regard to its ability to dissolve the host phase (crud) thereby liberating the entrained {sup 60}Co into a solution where it can be rinsed away.
Using a magnetic pressure drive, an absolute measurement of stress and density along the principal compression isentrope is obtained for solid aluminum to 240 GPa. Reduction of the free-surface velocity data relies on a backward integration technique, with approximate accounting for unknown systematic errors in experimental timing. Maximum experimental uncertainties are {+-}4.7% in stress and {+-}1.4% in density, small enough to distinguish between different equation-of-state (EOS) models. The result agrees well with a tabular EOS that uses an empirical universal zero-temperature isotherm.
Spectral imaging where a complete spectrum is collected from each of a series of spatial locations (1D lines, 2D images or 3D volumes) is now available on a wide range of analytical tools - from electron and x-ray to ion beam instruments. With this capability to collect extremely large spectral images comes the need for automated data analysis tools that can rapidly and without bias reduce a large number of raw spectra to a compact, chemically relevant, and easily interpreted representation. It is clear that manual interrogation of individual spectra is impractical even for very small spectral images (< 5000 spectra). More typical spectral images can contain tens of thousands to millions of spectra, which given the constraint of acquisition time may contain between 5 and 300 counts per 1000-channel spectrum. Conventional manual approaches to spectral image analysis such as summing spectra from regions or constructing x-ray maps are prone to bias and possibly error. One way to comprehensively analyze spectral image data, which has been automated, is to utilize an unsupervised self-modeling multivariate statistical analysis method such as multivariate curve resolution (MCR). This approach has proven capable of solving a wide range of analytical problems based upon the counting of x-rays (SEM/STEM-EDX, XRF, PIXE), electrons (EELS, XPS) and ions (TOF-SIMS). As an example of the MCR approach, a STEM x-ray spectral image from a ZrB2-SiC composite was acquired and analyzed. The data were generated in a FEI Tecnai F30-ST TEM/STEM operated at 300kV, equipped with an EDAX SUTW x-ray detector. The spectral image was acquired with the TIA software on the STEM at 128 by 128 pixels (12nm/pixel) for 100msec dwell per pixel (total acquisition time was 30 minutes) with a probe of approximately the same size as each pixel. Each spectrum in the image had, on average, 500 counts. The calculation took 5 seconds on a PC workstation with dual 2.4GHz PentiumIV Xeon processors and 2Gbytes of RAM and resulted in four chemically relevant components, which are shown in Figure 1. The analysis region was at a triple junction of three ZrB2 grains that contained zirconium oxide, aluminum oxide and a glass phase. The power of unbiased statistical methods, such as MCR as applied here, is that no a priori knowledge of the material's chemistry is required. The algorithms, in this case, effectively reduced over 16,000 2000-channel spectra (64Mbytes) to four images and four spectral shapes (72kbytes), which in this case represent chemical phases. This three order of magnitude compression is achieved rapidly with no loss of chemical information. There is also the potential to correlate multiple analytical techniques like, for example, EELS and EDS in the STEM adding sensitivity to light elements as well as bonding information for EELS to the more comprehensive spectral coverage of EDS.
Reducing agricultural water use in arid regions while maintaining or improving economic productivity of the agriculture sector is a major challenge. Controlled environment agriculture (CEA, or, greenhouse agriculture) affords advantages in direct resource use (less land and water required) and productivity (i.e., much higher product yield and quality per unit of resources used) relative to conventional open-field practices. These advantages come at the price of higher operating complexity and costs per acre. The challenge is to implement and apply CEA such that the productivity and resource use advantages will sufficiently outweigh the higher operating costs to provide for overall benefit and viability. This project undertook an investigation of CEA for livestock forage production as a water-saving alternative to open-field forage production in arid regions. Forage production is a large consumer of fresh water in many arid regions of the world, including the southwestern U.S. and northern Mexico. With increasing competition among uses (agriculture, municipalities, industry, recreation, ecosystems, etc.) for limited fresh water supplies, agricultural practice alternatives that can potentially maintain or enhance productivity while reducing water use warrant consideration. The project established a pilot forage production greenhouse facility in southern New Mexico based on a relatively modest and passive (no active heating or cooling) system design pioneered in Chihuahua, Mexico. Experimental operations were initiated in August 2004 and carried over into early-FY05 to collect data and make initial assessments of operational and technical system performance, assess forage nutrition content and suitability for livestock, identify areas needing improvement, and make initial assessment of overall feasibility. The effort was supported through the joint leveraging of late-start FY04 LDRD funds and bundled CY2004 project funding from the New Mexico Small Business Technical Assistance program at Sandia. Despite lack of optimization with the project system, initial results show the dramatic water savings potential of hydroponic forage production compared with traditional irrigated open field practice. This project produced forage using only about 4.5% of the water required for equivalent open field production. Improved operation could bring water use to 2% or less. The hydroponic forage production system and process used in this project are labor intensive and not optimized for minimum water usage. Freshly harvested hydroponic forage has high moisture content that dilutes its nutritional value by requiring that livestock consume more of it to get the same nutritional content as conventional forage. In most other aspects the nutritional content compares well on a dry weight equivalent basis with other conventional forage. More work is needed to further explore and quantify the opportunities, limitations, and viability of this technique for broader use. Collection of greenhouse environmental data in this project was uniquely facilitated through the implementation and use of a self-organizing, wirelessly networked, multi-modal sensor system array with remote cell phone data link capability. Applications of wirelessly networked sensing with improved modeling/simulation and other Sandia technologies (e.g., advanced sensing and control, embedded reasoning, modeling and simulation, materials, robotics, etc.) can potentially contribute to significant improvement across a broad range of CEA applications.
Video and image data are knowledge-rich sources of information, but their utility for current and future systems is limited without autonomous methods for understanding and characterizing their content. Semantic-based video understanding may benefit systems dedicated to the detection of insiders, alarm patterns, unauthorized activities in material monitoring applications, etc. A direct benefit of this technology is not only intelligent alarm analysis, but the ability to browse and perform query-based searches for useful and interesting information after video data has been acquired and stored. These searches can provide a tremendous benefit for use in intelligence agency, government, military, and DOE site investigations. This report provides an initial investigation into the algorithms and methods needed to characterize and understand video content. Such algorithms include background modeling, detecting dynamic image regions, grouping dynamic pixels into coherent objects, and robust tracking strategies. With solid approaches for addressing these problems, analysis can be performed seeking to recognize distinctive objects and their motions leading to semantic-based video searches.
This SAND report provides the technical progress through October 2004 of the Sandia-led project, %22Carbon Sequestration in Synechococcus Sp.: From Molecular Machines to Hierarchical Modeling,%22 funded by the DOE Office of Science Genomes to Life Program. Understanding, predicting, and perhaps manipulating carbon fixation in the oceans has long been a major focus of biological oceanography and has more recently been of interest to a broader audience of scientists and policy makers. It is clear that the oceanic sinks and sources of CO2 are important terms in the global environmental response to anthropogenic atmospheric inputs of CO2 and that oceanic microorganisms play a key role in this response. However, the relationship between this global phenomenon and the biochemical mechanisms of carbon fixation in these microorganisms is poorly understood. In this project, we will investigate the carbon sequestration behavior of Synechococcus Sp., an abundant marine cyanobacteria known to be important to environmental responses to carbon dioxide levels, through experimental and computational methods. This project is a combined experimental and computational effort with emphasis on developing and applying new computational tools and methods. Our experimental effort will provide the biology and data to drive the computational efforts and include significant investment in developing new experimental methods for uncovering protein partners, characterizing protein complexes, identifying new binding domains. We will also develop and apply new data measurement and statistical methods for analyzing microarray experiments. Computational tools will be essential to our efforts to discover and characterize the function of the molecular machines of Synechococcus. To this end, molecular simulation methods will be coupled with knowledge discovery from diverse biological data sets for high-throughput discovery and characterization of protein-protein complexes. In addition, we will develop a set of novel capabilities for inference of regulatory pathways in microbial genomes across multiple sources of information through the integration of computational and experimental technologies. These capabilities will be applied to Synechococcus regulatory pathways to characterize their interaction map and identify component proteins in these - 4 - pathways. We will also investigate methods for combining experimental and computational results with visualization and natural language tools to accelerate discovery of regulatory pathways. The ultimate goal of this effort is develop and apply new experimental and computational methods needed to generate a new level of understanding of how the Synechococcus genome affects carbon fixation at the global scale. Anticipated experimental and computational methods will provide ever-increasing insight about the individual elements and steps in the carbon fixation process, however relating an organism's genome to its cellular response in the presence of varying environments will require systems biology approaches. Thus a primary goal for this effort is to integrate the genomic data generated from experiments and lower level simulations with data from the existing body of literature into a whole cell model. We plan to accomplish this by developing and applying a set of tools for capturing the carbon fixation behavior of complex of Synechococcus at different levels of resolution. Finally, the explosion of data being produced by high-throughput experiments requires data analysis and models which are more computationally complex, more heterogeneous, and require coupling to ever increasing amounts of experimentally obtained data in varying formats. These challenges are unprecedented in high performance scientific computing and necessitate the development of a companion computational infrastructure to support this effort. More information about this project, including a copy of the original proposal, can be found at www.genomes-to-life.org Acknowledgment We want to gratefully acknowledge the contributions of the GTL Project Team as follows: Grant S. Heffelfinger1*, Anthony Martino2, Andrey Gorin3, Ying Xu10,3, Mark D. Rintoul1, Al Geist3, Matthew Ennis1, Hashimi Al-Hashimi8, Nikita Arnold3, Andrei Borziak3, Bianca Brahamsha6, Andrea Belgrano12, Praveen Chandramohan3, Xin Chen9, Pan Chongle3, Paul Crozier1, PguongAn Dam10, George S. Davidson1, Robert Day3, Jean Loup Faulon2, Damian Gessler12, Arlene Gonzalez2, David Haaland1, William Hart1, Victor Havin3, Tao Jiang9, Howland Jones1, David Jung3, Ramya Krishnamurthy3, Yooli Light2, Shawn Martin1, Rajesh Munavalli3, Vijaya Natarajan3, Victor Olman10, Frank Olken4, Brian Palenik6, Byung Park3, Steven Plimpton1, Diana Roe2, Nagiza Samatova3, Arie Shoshani4, Michael Sinclair1, Alex Slepoy1, Shawn Stevens8, Chris Stork1, Charlie Strauss5, Zhengchang Su10, Edward Thomas1, Jerilyn A. Timlin1, Xiufeng Wan11, HongWei Wu10, Dong Xu11, Gong-Xin Yu3, Grover Yip8, Zhaoduo Zhang2, Erik Zuiderweg8 *Author to whom correspondence should be addressed (gsheffe%40sandia.gov) 1. Sandia National Laboratories, Albuquerque, NM 2. Sandia National Laboratories, Livermore, CA 3. Oak Ridge National Laboratory, Oak Ridge, TN 4. Lawrence Berkeley National Laboratory, Berkeley, CA 5. Los Alamos National Laboratory, Los Alamos, NM 6. University of California, San Diego 7. University of Illinois, Urbana/Champaign 8. University of Michigan, Ann Arbor 9. University of California, Riverside 10. University of Georgia, Athens 11. University of Missouri, Columbia 12. National Center for Genome Resources, Santa Fe, NM Sandia is a multiprogram laboratory operated by Sandia Corporation, a Lockheed Martin Company, for the United States Department of Energy's National Nuclear Security Administration under contract DE-AC04-94AL85000.
Local topological modification is widely used to improve mesh quality after automatic generation of tetrahedral and quadrilateral meshes. These same techniques are also used to support adaptive refinement of these meshes. In contrast, few methods are known for locally modifying the topology of hexahedral meshes. Most efforts to do this have been based on fixed transition templates or global refinement. In contrast, a dual-based 'pillowing' method has been used which, while local, is still quite restricted in its application, and is typically applied in a template-based fashion. In this presentation, I will describe the generalization of a dual-based approach to the local topological modification of hex meshes and its application to clean up hexahedral meshes. A set of three operations for locally modifying hex mesh topology has been shown to reproduce the so-called 'flipping' operations described by Bern et. al as well as other commonly-used refinement templates. I will describe the implementation of these operators and their application to real meshes. Challenging aspects of this work have included visualization of a hex mesh and its dual (especially for poor-quality meshes); the incremental modification of both the primal (i.e. the mesh) and the dual simultaneously; and the interactive steering of these operations with the goal of improving hex meshes which would otherwise have unacceptable quality. These aspects will be discussed in the context of improving hex meshes generated by curve contraction-based whisker weaving. Application of these techniques for improving other hexahedral mesh types, for example those resulting from tetrahedral subdivision, will also be discussed.
Many highly oscillatory circuits have a wide separation of time scales between the underlying oscillation and the behavior of interest. This is particularly true of communication circuits. Multiple-time Partial Differential Equation (MPDE) methods offer substantial speed-up for these circuits by introducing a periodic artificial time variable that represents the highly oscillatory behavior. This leaves just the slowly changing behavior of interest, which can be integrated with much larger steps. One problem of particular interest is the larger initial condition that must be specified for this periodic artificial time variable. One possible solution is to formulate an optimization problem in the hopes of increasing the step sizes taken in the slow time direction. This talk will discuss one possible unconstrained optimization problem for determining this initial condition. Numerical results and comparisons to several other initial condition strategies will be presented in addition to MPDE background research and implementation issues.
Experimental evidence and corresponding theoretical analyses have led to the conclusion that the system composed of Xe hollow atom states, that produce a characteristic Xe(L) spontaneous emission spectrum at 1 {at} 2.9 {angstrom} and arise from the excitation of Xe clusters with an intense pulse of 248 nm radiation propagating in a self-trapped plasma channel, closely represents the ideal situation sought for amplification in the multikilovolt region. The key innovation that is central to all aspects of the proposed work is the controlled compression of power to the level ({approx} 10{sup 20} W/cm{sup 3}) corresponding to the maximum achieved by thermonuclear events. Furthermore, since the x-ray power that is produced appears in a coherent form, an entirely new domain of physical interaction is encountered that involves states of matter that are both highly excited and highly ordered. Moreover, these findings lead to the concept of 'photonstaging', an idea which offers the possibility of advancing the power compression by an additional factor of {approx} 10{sup 9} to {approx} 10{sup 29} W/cm{sup 3}. In this completely unexplored regime, g-ray production ({h_bar}{omega}{sub {gamma}} {approx} 1 MeV) is expected to be a leading process. A new technology for the production of very highly penetrating radiation would then be available. The Xe(L) source at {h_bar}{omega}{sub x} {approx} 4.5 keV can be applied immediately to the experimental study of many aspects of the coupling of intense femtosecond x-ray pulses to materials. In a joint collaboration, the UIC group and Sandia plan to explore the following areas. These are specifically, (1) anomalous electromagnetic coupling to solid state materials, (2) 3D nanoimaging of solid matter and hydrated biological materials (e.g. interchromosomal linkers and actin filaments in muscle), and (3) EMP generation with attosecond x-rays.
Because many solid objects, both stationary and mobile, will be present in an indoor environment, the design of an indoor aerosol cloud finding lidar (light detection and ranging) instrument presents a number of challenges. The cloud finder must be able to discriminate between these solid objects and aerosol clouds as small as 1-meter in depth in order to probe suspect clouds. While a near IR ({approx}1.5-{micro}m) laser is desirable for eye-safety, aerosol scattering cross sections are significantly lower in the near-IR than at visible or W wavelengths. The receiver must deal with a large dynamic range since the backscatter from solid object will be orders of magnitude larger than for aerosol clouds. Fast electronics with significant noise contributions will be required to obtain the necessary temporal resolution. We have developed a laboratory instrument to detect aerosol clouds in the presence of solid objects. In parallel, we have developed a lidar performance model for performing trade studies. Careful attention was paid to component details so that results obtained in this study could be applied towards the development of a practical instrument. The amplitude and temporal shape of the signal return are analyzed for discrimination of aerosol clouds in an indoor environment. We have assessed the feasibility and performance of candidate approaches for a fieldable instrument. With the near-IR PMT and a 1.5-{micro}m laser source providing 20-{micro}J pulses, we estimate a bio-aerosol detection limit of 3000 particles/l.
Solution-based synthesis is a powerful approach for creating nano-structured materials. Although there have been significant recent successes in its application to fabricating nanomaterials, the general principles that control solution synthesis are not well understood. The purpose of this LDRD project was to develop the scientific principles required to design and build unique nanostructures in crystalline oxides and II/VI semiconductors using solution-based molecular self-assembly techniques. The ability to synthesize these materials in a range of different nano-architectures (from controlled morphology nanocrystals to surface templated 3-D structures) has provided the foundation for new opportunities in such areas as interactive interfaces for optics, electronics, and sensors. The homogeneous precipitation of ZnO in aqueous solution was used primarily as the model system for the project. We developed a low temperature, aqueous solution synthesis route for preparation of large arrays of oriented ZnO nanostructures. Through control of heterogeneous nucleation and growth, methods to predicatively alter the ZnO microstructures by tailoring the surface chemistry of the crystals were established. Molecular mechanics simulations, involving single point energy calculations and full geometry optimizations, were developed to assist in selecting appropriate chemical systems and understanding physical adsorption and ultimately growth mechanisms in the design of oxide nanoarrays. The versatility of peptide chemistry in controlling the formation of cadmium sulfide nanoparticles and zinc oxide/cadmium sulfide heterostructures was also demonstrated.
There is a general lack of compact electromagnetic radiation sources between 1 and 10 terahertz (THz). This a challenging spectral region lying between optical devices at high frequencies and electronic devices at low frequencies. While technologically very underdeveloped the THz region has the promise to be of significant technological importance, yet demonstrating its relevance has proven difficult due to the immaturity of the area. While the last decade has seen much experimental work in ultra-short pulsed terahertz sources, many applications will require continuous wave (cw) sources, which are just beginning to demonstrate adequate performance for application use. In this project, we proposed examination of two potential THz sources based on intersubband semiconductor transitions, which were as yet unproven. In particular we wished to explore quantum cascade lasers based sources and electronic based harmonic generators. Shortly after the beginning of the project, we shifted our emphasis to the quantum cascade lasers due to two events; the publication of the first THz quantum cascade laser by another group thereby proving feasibility, and the temporary shut down of the UC Santa Barbara free-electron lasers which were to be used as the pump source for the harmonic generation. The development efforts focused on two separate cascade laser thrusts. The ultimate goal of the first thrust was for a quantum cascade laser to simultaneously emit two mid-infrared frequencies differing by a few THz and to use these to pump a non-linear optical material to generate THz radiation via parametric interactions in a specifically engineered intersubband transition. While the final goal was not realized by the end of the project, many of the completed steps leading to the goal will be described in the report. The second thrust was to develop direct THz QC lasers operating at terahertz frequencies. This is simpler than a mixing approach, and has now been demonstrated by a few groups with wavelengths spanning 65-150 microns. We developed and refined the MBE growth for THz for both internally and externally designed QC lasers. Processing related issues continued to plague many of our demonstration efforts and will also be addressed in this report.
A numerical screening study of the interaction between a penetrator and a geological target with a preformed hole has been carried out to identify the main parameters affecting the penetration event. The planning of the numerical experiment was based on the orthogonal array OA(18,7,3,2), which allows 18 simulation runs with 7 parameters at 3 levels each. The strength of 2 of the array allows also for two-factor interaction studies. The seven parameters chosen for this study are: penetrator offset, hole diameter, hole taper, vertical and horizontal velocity of the penetrator, angle of attack of the penetrator and target material. The analysis of the simulation results has been based on main effects plots and analysis of variance (ANOVA), and it has been performed using three metrics: the maximum values of the penetration depth, penetrator deceleration and plastic strain in the penetrator case. This screening study shows that target material has a major influence on penetration depth and penetrator deceleration, while penetrator offset has the strongest effect on the maximum plastic strain.
The proposed Yucca Mountain repository is anticipated to be the first facility for long-term disposal of commercial spent nuclear fuel and high-level radioactive waste in the United States. The facility, located in the southern Nevada desert, is currently in the planning stages with initial exploratory excavations completed. It is an underground facility mined into the tuffaceous volcanic rocks that sit above the local water table. The focus of the work described in this paper is the development of radionuclide absorbers or 'getter' materials for neptunium (Np), iodine (I), and technetium (Tc) for potential deployment in the repository. 'Getter' materials retard the migration of radionuclides through sorption, reduction, or other chemical and physical processes, thereby slowing or preventing the release and transport of radionuclides. An overview of the objectives and approaches utilized in this work with respect to materials selection and modeling of ion 'getters' is presented. The benefits of the 'getter' development program to the United States Department of Energy (US DOE) are outlined.
We have investigated the potential for intense particle beam surface modification to improve the mechanical properties of materials commonly used in the human body for contact surfaces in, for example, hip and knee implants. The materials studied include Ultra-High Molecular Weight Polyethylene (UHMWPE), Ti-6Al-4Al (titanium alloy), and Co-Cr-Mo alloy. Samples in flat form were exposed to both ion and electron beams (UHMWPE), and to ion beam treatment (metals). Post-analysis indicated a degradation in bulk properties of the UHMWPE, except in the case of the lightest ion fluence tested. A surface-alloyed Hf/Ti layer on the Ti-6Al-4V is found to improve surface wear durability, and have favorable biocompatibility. A promising nanolaminate ceramic coating is applied to the Co-Cr-Mo to improve surface hardness.
Large-scale three dimensional Discrete Element simulations of granular flow in a modified split-bottom Couette cell for packs of up to 180,000 mono-disperse spheres are presented and compared with experiments. We find that the velocity profiles collapse onto a universal curve not only at the surface but also in the bulk of the pack until slip between layers becomes significant. In agreement with experiment, we find similar relations between the cell geometry and parameters involved in rescaling the velocities at the surface and in the bulk. Likewise, a change in the shape of the shear zone is observed as predicted for tall packs once the center of the shear zone is correctly defined; although the transition does not appear to be first order. Finally, the effect of cohesion is considered as a means to test the theoretical predictions.
This report summarizes the results of a five-year effort to understand the mechanisms and develop models that predict the corrosion of refractories in oxygen-fuel glass-melting furnaces. Thermodynamic data for the Si-O-(Na or K) and Al-O-(Na or K) systems are reported, allowing equilibrium calculations to be performed to evaluate corrosion of silica- and alumina-based refractories under typical furnace operating conditions. A detailed analysis of processes contributing to corrosion is also presented. Using this analysis, a model of the corrosion process was developed and used to predict corrosion rates in an actual industrial glass furnace. The rate-limiting process is most likely the transport of NaOH(gas) through the mass-transport boundary layer from the furnace atmosphere to the crown surface. Corrosion rates predicted on this basis are in better agreement with observation than those produced by any other mechanism, although the absolute values are highly sensitive to the crown temperature and the NaOH(gas) concentration at equilibrium and at the edge of the boundary layer. Finally, the project explored the development of excimer laser induced fragmentation (ELIF) fluorescence spectroscopy for the detection of gas-phase alkali hydroxides (e.g., NaOH) that are predicted to be the key species causing accelerated corrosion in these furnaces. The development of ELIF and the construction of field-portable instrumentation for glass furnace applications are reported and the method is shown to be effective in industrial settings.
The large number of government and industry activities supporting the Unit of Action (UA), with attendant documents, reports and briefings, can overwhelm decision-makers with an overabundance of information that hampers the ability to make quick decisions often resulting in a form of gridlock. In particular, the large and rapidly increasing amounts of data and data formats stored on UA Advanced Collaborative Environment (ACE) servers has led to the realization that it has become impractical and even impossible to perform manual analysis leading to timely decisions. UA Program Management (PM UA) has recognized the need to implement a Decision Support System (DSS) on UA ACE. The objective of this document is to research the commercial Knowledge Discovery and Data Mining (KDDM) market and publish the results in a survey. Furthermore, a ranking mechanism based on UA ACE-specific criteria has been developed and applied to a representative set of commercially available KDDM solutions. In addition, an overview of four R&D areas identified as critical to the implementation of DSS on ACE is provided. Finally, a comprehensive database containing detailed information on surveyed KDDM tools has been developed and is available upon customer request.
This paper describes the development of a set of software tools useful for analyzing ultra-wideband (UWB) antennas and structures. These tools are used to perform finite difference time domain (FDTD) simulation of a conical antenna with continuous wave (CW) and UWB pulsed excitations. The antenna is analyzed using spherical coordinate-based FDTD equations that are derived from first principles. The simulation results for CW excitation are compared to simulation and measured results from published sources; the results for UWB excitation are new.
This document contains a description of the verification and validation process used for the RADTRAN 5.5 code. The verification and validation process ensured the proper calculational models and mathematical and numerical methods were used in the RADTRAN 5.5 code for the determination of risk and consequence assessments. The differences between RADTRAN 5 and RADTRAN 5.5 are the addition of tables, an expanded isotope library, and the additional User-Defined meteorological option for accident dispersion. 3
As interest in 3D face recognition increases the importance of the initial alignment problem does as well. In this paper we present a method utilizing the registered 2D color and range image of a face to automatically identify the eyes, nose, and mouth. These features are important to initially align faces in both standard 2D and 3D face recognition algorithms. For our algorithm to run as fast as possible, we focus on the 2D color information. This allows the algorithm to run in approximately 4 seconds on a 640X480 image with registered range data. On a database of 1,500 images the algorithm achieved a facial feature detection rate of 99.6% with 0.4% of the images skipped due to hair obstruction of the face.
A continuum-scale, evolutionary model of bubble nucleation, growth and He release for aging metal tritides is described which accounts for major features of the tritide database. Bubble nucleation, modeled as self-trapping of interstitially diffusing He atoms, occurs during the first few days following tritium introduction into the metal. Bubble growth by dislocation loop punching yields good agreement between He atomic volumes and bubble pressures determined from bulk swelling and 3He NMR data. The bubble spacing distribution determined from NMR is shown to remain fixed with age, justifying the separation of nucleation and growth phases and providing a sensitive test of the growth formulation. Late in life, bubble interactions are proposed to produce cooperative stress effects, which lower the bubble pressure. Helium generated near surfaces and surface-connected porosity accounts for the low-level early helium release. Use of an average ligament stress criterion predicts an onset of inter-bubble fracture in good agreement with the He/Metal ratio observed for rapid He release. From the model, it is concluded that He retention can be controlled through control of bubble nucleation.
The diminished response of thermoluminescent phosphors over time is a well-documented challenge to thermoluminescent dosimetry. Wide ranges in fading rates for various phosphor types have been reported, making it necessary for many external dosimetry programs to perform individual studies on thermoluminescent fade. Sandia National Laboratories currently uses the Thermo 8802 LiF:Mg,Ti thermoluminescent dosimeter (TLD) in its personnel external dosimetry program. Doses received in the field are calculated by applying a fade algorithm published by the manufacturer to TLD readings. Since the algorithm was established by characterizing the diminished response of a TLD similar to the 8802, Sandia chose to model its fade study after the analysis done by Thermo. As a result, the parameters of each experiment were comparable, and data from the two studies were compared to determine whether or not the current algorithm should be modified specifically for use at Sandia. Cards were irradiated using an internal 90Sr/90Y source, and pre- and post-irradiation fading rates were monitored over a period of 18 wk. While significant fading was demonstrated, results closely matched those found in the original Thermo study.
The molten salt Flibe, a combination of lithium and beryllium fluorides studied for molten salt fission reactors, has been proposed as a breeder and coolant for fusion applications. The melting points of 2LiF-BeF2 and LiF-BeF2 are 460°C and 363°C, but LiF-BeF2 is rather viscous and has less lithium for breeding. In the Advanced Power Extraction (APEX) Program, concepts with a free flowing liquid for the first wall and blanket were investigated. Flinabe (a mixture of LiF, BeF2 and NaF) was selected for a molten salt design because a melting temperature below 350°C appeared possible and this provided an attractive operating temperature window for a reactor. To confirm that a ternary salt with a low melting temperature existed, several combinations of the fluoride salts, LiF, NaF and BeF2, were melted in a stainless steel crucible under vacuum. One had an apparent melting temperature of 305°C. The test system, preparation of the mixtures, melting procedures and temperature curves for the melting and cooling are presented along with the apparent melting points. Thermal modeling of the salt pool and crucible is reported in an accompanying paper.
PDC drill bit performance has been greatly improved over the past three decades by innovations in bit design and how these designs are applied. The next leap forward is most likely to come from using high-speed, real-time downhole data to optimize the performance of these sophisticated bits on an application-by-application basis. By effectively managing conditions of lateral, axial and torsional acceleration, damage to cutting structures can be minimized for improved penetration rates. Avoiding these damaging vibrations is essential to increasing bit durability and improving overall drilling economics. This paper describes one set of independent drilling optimization results obtained as part of a series of controlled demonstrations of PDC bits. Sandia National Laboratories on behalf of the U. S. Department of Energy (DOE) managed this work. The effort was organized as a Cooperative Research and Development Agreement (CRADA) established between Sandia and four bit manufacturers with DOE funding and in-kind contributions by the industry partners. The goal of this CRADA was to demonstrate drag bit performance in formations with degrees of hardness typical of those encountered while drilling geothermal wells. The test results indicate that the surface weight-on-bit (WOB), revolutions per minute (RPM) and torque readings traditionally used to guide adjustments in the drilling parameters do not always provide the true picture of what is really taking place at the bit. Instead, a holistic approach combining traditional methods of optimization together with high-speed, real-time data enables far better decisions for improving bit performance and avoiding damaging situations that may have otherwise gone unnoticed.
The Z-Pinch Power Plant uses the results from Sandia National Laboratories' Z accelerator in a power plant application to generate energy pulses using inertial confinement fusion. A collaborative project has been initiated by Sandia to investigate the scientific principles of a power generation system using this technology. Research is under way to develop an integrated concept that describes the operational issues of a 1000 MW electrical power plant. Issues under consideration include: 1-20 gigajoule fusion pulse containment, repetitive mechanical connection of heavy hardware, generation of terawatt pulses every 10 seconds, recycling often thousand tons of steel, and manufacturing of millions of hohlraums and capsules per year. Additionally, waste generation and disposal issues are being examined. This paper describes the current concept for the plant and also the objectives for future research.
An experimental investigation is made into the fluid mechanics and heat transfer of a circular cylinder immersed in a wall-bounded turbulent mixed-convection flow of water. The cylinder is oriented spanwise to the forced channel flow and within the thermal boundary layer of the heated lower wall. The flow channel is capped with a cold, near-adiabatic upper wall producing a fully turbulent gap Rayleigh number of 108. A low-speed crossflow is applied to advect the turbulent thermal plumes over the cylinder surface. We present spatially resolved cylinder-surface heat-flux data alongside 2-D PIV imaging of the streamwise and wall-normal velocity components for two flow conditions in the mixed-convection heat-transfer regime. The measured cylinder-wake flowfield reflects the complex coupling between the separated wake flow, the highly turbulent freestream and the buoyant wall and cylinder boundary layers. A method for measurement of spatially resolved surface heat fluxes based on the measured cylinder-surface temperature distribution and a well-posed two-dimensional solution to the conduction problem in the cylinder wall is presented. The resulting spatially resolved flux measurements show enhanced surface heat transfer, which results from the intense buoyancy generated free-stream turbulence and mixing in the cylinder wake. This work extends the literature on thermal convection with crossflow well into the turbulent regime and is, to our knowledge, the first investigation of surface heat-transfer to an object of engineering importance placed in this type of turbulent mixed-convection flowfield. The data are currently being utilized for validation of mixed-convection turbulence models at Sandia and comparisons between the computational and experimental results are presented.
Forward-to-reverse bias step-recovery measurements were performed on In.07Ga.93N/GaN and Al.36Ga.64N/Al.46Ga.54N quantum-well (QW) light-emitting diodes grown on sapphire. With the QW sampling the minority-carrier hole density at a single position, distinctive two-phase optical decay curves were observed. Using diffusion equation solutions to self-consistently model both the electrical and optical responses, hole transport parameters tp = 758 {+-} 44 ns, Lp = 588 {+-} 45 nm, and up = 0.18 {+-} 0.02 cm2/Vs were obtained for GaN. The mobility was thermally activated with an activation energy of 52 meV, suggesting trap-modulated transport. Optical measurements of sub-bandgap peaks exhibited slow responses approaching the bulk lifetime. For Al.46Ga.54N, a longer lifetime of tp = 3.0 us was observed, and the diffusion length was shorter, Lp = 280 nm. Mobility was an order of magnitude smaller than in GaN, up = 10-2 cm2/Vs, and was insensitive to temperature, suggesting hole transport through a network of defects.
We report micro-Raman studies of self-heating in an AlGaN/GaN heterostructure field-effect transistor using below (visible 488.0 nm) and near (UV 363.8 nm) GaN band-gap excitation. The shallow penetration depth of the UV light allows us to measure temperature rise ({Delta}T) in the two-dimensional electron gas (2DEG) region of the device between drain and source. Visible light gives the average {Delta}T in the GaN layer, and that of the SiC substrate, at the same lateral position. Combined, we depth profile the self-heating. Measured {Delta}T in the 2DEG is consistently over twice the average GaN-layer value. Electrical and thermal transport properties are simulated. We identify a hotspot, located at the gate edge in the 2DEG, as the prevailing factor in the self-heating.
Order-of-accuracy verification is necessary to ensure that software correctly solves a given set of equations. One method to verify the order of accuracy of a code is the method of manufactured solutions. In this study, a manufactured solution has been derived and implemented that allows verification of not only the Euler, Navier-Stokes, and Reynolds-Averaged Navier-Stokes (RANS) equation sets, but also some of their associated boundary conditions (BC's): slip, no-slip (adiabatic and isothermal), and outflow (subsonic, supersonic, and mixed). Order-of-accuracy verification has been performed for the Euler and Navier-Stokes equations and these BC's in a compressible computational fluid dynamics code. All of the results shown are on skewed, non-uniform meshes. RANS results will be presented in a future paper. The observed order of accuracy was lower than the expected order of accuracy in two cases. One of these cases resulted in the identification and correction of a coding mistake in the CHAD gradient correction that was reducing the observed order of accuracy. This mistake would have been undetectable on a Cartesian mesh. During the search for the CHAD gradient correction problem, an unrelated coding mistake was found and corrected. The other case in which the observed order of accuracy was less than expected was a test of the slip BC; although no specific coding or formulation mistakes have yet been identified. After the correction of the identified coding mistakes, all of the aforementioned equation sets and BC's demonstrated the expected (or at least acceptable) order of accuracy except the slip condition.
Rate constants for the thermal dissociation of Si{sub 2}H{sub 6} are predicted with a novel transition state model. The saddle points for dissociation on the Si{sub 2}H{sub 6} potential energy surface are lower in energy than the corresponding separated products, as confirmed by high level ab initio quantum mechanical calculations. Thus, the dissociations of Si{sub 2}H{sub 6} to produce SiH{sub 2} + SiH{sub 4} (R1) and H{sub 3}SiSiH + H{sub 2} (R2) both proceed through tight inner transition states followed by loose outer transition states. The present 'dual' transition state model couples variational phase space theory treatments of the outer transition states with ab initio based fixed harmonic vibrator treatments of the inner transition states to obtain effective numbers of states for the two transition states acting in series. It is found that, at least near room temperature, such a dual transition state model is generally required for the proper description of each of the dissociations. Only at quite high temperatures, i.e., above 2000 K for (R1) and 600 K for (R2), does a single fixed inner transition state provide an adequate description. Similarly, only at quite low temperatures (below 100 and 10 K for (R1) and (R2), respectively) does a single outer transition state provide an adequate description. Pressure dependent rate constants are obtained from solutions to the multichannel master equation. These calculations confirm that dissociation channel (R2) is negligible under conditions relevant to the thermal chemical vapor deposition (CVD) processes. Rate constants for the chemical activation reactions, SiH{sub 2} + SiH{sub 4} {yields} Si{sub 2}H{sub 6} (R-1) and SiH{sub 2} + SiH{sub 4} {yields} H{sub 3}SiSiH + H{sub 2} (R3), are also evaluated within the dual transition state model. It is found that reaction R3 is the dominant channel for low pressures and high temperatures, i.e., below 100 Torr for temperatures above 1100 K.
In this article, we discuss in detail the addition of hydrogen atoms to diacetylene and the reverse dissociation reactions, H + C{sub 4}H{sub 2} {leftrightarrow} i-C{sub 4}H{sub 3} (R1) and H + C{sub 4}H{sub 2} n-C{sub 4}H{sub 3} (R2). The theory utilizes high-level electronic structure methodology to characterize the potential energy surface, Rice-Ramsperger-Kassel-Marcus (RRKM) theory to calculate microcanonical/J-resolved rate coefficients, and a two-dimensional master-equation approach to extract phenomenological (thermal) rate coefficients. Comparison is made with experimental results where they are available. The rate coefficients k{sub 1}(T, p) and k{sub 2}(T, p) are cast in forms that can be used in chemical kinetic modeling. In addition, we predict values of the heats of formation of i-C{sub 4}H{sub 3} and n-C{sub 4}H{sub 3} and discuss their importance in flame chemistry. Our basis-set extrapolated, quadratic-configuration-interaction with single and double excitations (and triple excitations added perturbatively), QCISD(T), predictions of these heats of formation at 298 K are 130.8 kcal/mol for n-C{sub 4}H{sub 3} and 119.3 kcal/mol for the i-isomer; multireference CI calculations with a nine-electron, nine-orbital, complete-active-space (CAS) reference wavefunction give just slightly larger values for these parameters. Our results are in good agreement with the recent focal-point analysis of Wheeler et al. (J. Chem. Phys. 2004, 121, 8800-8813), but they differ substantially for {Delta} H{sub f 298}{sup 0}(n-C{sub 4}H{sub 3}) with the earlier diffusion Monte Carlo predictions of Krokidis et al.
Many accelerators at Sandia National Laboratories utilize the Rimfire gas switch for high-voltage, high-power switching. Future accelerators will have increased performance requirements for switching elements. When designing improved versions of the Rimfire switch, there is a need for quick and accurate simulation of the electrical effects of geometry changes. This paper presents an advanced circuit model of the Rimfire switch that can be used for these simulations. The development of the model is shown along with comparisons to past models and experimental results.
The Sandia Lightning Simulator at Sandia National Laboratories can provide up to 200 kA for a simulated single lightning stroke, 100 kA for a subsequent stroke, and hundreds of Amperes of continuing current. It has recently been recommissioned after a decade of inactivity and the single-stroke capability demonstrated. The simulator capabilities, basic design components, upgrades, and diagnostic capabilities are discussed in this paper.
An experimental investigation is made into the fluid mechanics and heat transfer of a circular cylinder immersed in a wall-bounded turbulent mixed-convection flow of water. The cylinder is oriented spanwise to the forced channel flow and within the thermal boundary layer of the heated lower wall. The flow channel is capped with a cold, near-adiabatic upper wall producing a fully turbulent gap Rayleigh number of 10{sup 8}. A low-speed crossflow is applied to advect the turbulent thermal plumes over the cylinder surface. We present spatially resolved cylinder-surface heat-flux data alongside 2-D PIV imaging of the streamwise and wall-normal velocity components for two flow conditions in the mixed-convection heat-transfer regime. The measured cylinder-wake flowfield reflects the complex coupling between the separated wake flow, the highly turbulent freestream and the buoyant wall and cylinder boundary layers. A method for measurement of spatially resolved surface heat fluxes based on the measured cylinder-surface temperature distribution and a well-posed two-dimensional solution to the conduction problem in the cylinder wall is presented. The resulting spatially resolved flux measurements show enhanced surface heat transfer, which results from the intense buoyancy generated free-stream turbulence and mixing in the cylinder wake. This work extends the literature on thermal convection with crossflow well into the turbulent regime and is, to our knowledge, the first investigation of surface heat-transfer to an object of engineering importance placed in this type of turbulent mixed-convection flowfield. The data are currently being utilized for validation of mixed convection turbulence models at Sandia and comparisons between the computational and experimental results are presented.
The WLF equation is typically used to describe the dependence of polymer mobility on temperature at atmospheric pressure. Tests at different pressures would at least require different WLF parameterization. Completely different tests, for example, probing the temperature dependence of mobility at constant density, would require even greater modifications. By performing molecular dynamics simulations on simple chain molecules equilibrated at different thermodynamic states, we have shown that the mobility depends in a more general sense on the potential energy density of the system. That is, mobilities for any equilibrated state collapse onto one master curve when plotted against the potential energy density. Moreover, this relationship can be fit by either a 'generalized' WLF equation or by a power-law relationship observed in critical phenomena. When this mobility relationship is used within a rheologically simple, thermodynamically consistent, viscoelastic framework, quantitative agreement is seen between experimental data and theoretical predictions on a range of tests covering enthalpy relaxation to mechanical yield to physical aging.
In this study, we describe the extension of the 2-d preliminary design bluff body drag estimation tool developed by De Chant1 to apply for 3-d flows. As with the 2-d method, the 3-d extension uses a combined approximate Green's function/Gram-Charlier series approach to retain the body geometry information. Whereas, the 2-d methodology relied solely upon the use of small disturbance theory for the inviscid flow field associated with the body of interest to estimate the near-field initial conditions, e.g. velocity defect, the 3-d methodology uses both analytical (where available) and numerical inviscid solutions. The defect solution is then used as an initial condition in an approximate 3-d Green's function solution. Finally, the Green's function solution is matched to the 3-d analog of the classical 2-d Gram-Charlier series and then integrated to yield the net form drag on the bluff body. Preliminary results indicate that drag estimates computed are of accuracy equivalent to the 2-d method for flows with large separation, i.e. less than 20% relative error. As was the lower dimensional method, the 3-d concept is intended to be a supplement to turbulent Navier-Stokes and experimental solution for estimating drag coefficients over blunt bodies.
The Method of Nearby Problems is employed to generate exact solutions to equations 'nearby' the steady and unsteady Burgers equation. Burgers equation is chosen because of the existence of exact solutions, and these exact solutions are discussed. Legendre polynomials are used to derive the exact solutions to the nearby problems, and the application of Legendre polynomials for both 1D and 2D problems is also discussed. Results are presented for the steady-state Burgers equation corresponding to a viscous shock wave for Reynolds numbers of 8, 16, and 512. The low Reynolds number cases are well approximated by 10th order Legendre polynomial fits, while the high Reynolds number case is not. The unsteady Burgers equation corresponding to coalescence of two viscous shock waves at a Reynolds number of 8 is also examined. Preliminary results indicate that further investigation is required to accurately capture this 2D solution.
We performed calculations to investigate the classical theories of chain branching and thermal--run--away that lead to the rapid oxidation of fuels. Mathematically, both theories infer the existence of eigenvalues with positive real parts i.e., explosive modes. We found in studies of homogeneous hydrogen--air and the methane--air mixtures that when ignition is initiated by a sufficiently high initial temperature, the transient response of the system exhibits two stages. The first stage is characterized by the existence of explosive modes. The ensuing second stage consists of fast exponential decay modes that bring the system to its equilibrium point. We demonstrated with two examples that the existence of explosive modes is not a necessary condition for the existence of a premixed flame. Homogeneous ignition calculations for mixtures with an initial concentration of radical species suggest that the diffusive transport of radical species is probably responsible for the lack of explosive modes in premixed flames.
Modal analysis of three-dimensional structures frequently involves finite element discretizations with millions of unknowns and requires computing hundreds or thousands of eigenpairs. In this presentation we review methods based on domain decomposition for such eigenspace computations in structural dynamics. We distinguish approaches that solve the eigenproblem algebraically (with minimal connections to the underlying partial differential equation) from approaches that tightly couple the eigensolver with the partial differential equation.
Low-temperature combustion concepts that utilize cooled EGR, early/retarded injection, high swirl ratios, and modest compression ratios have recently received considerable attention. To understand the combustion and, in particular, the soot formation process under these operating conditions, a modeling study was carried out using the KIVA-3V code with an improved phenomenological soot model. This multi-step soot model includes particle inception, surface growth, surface oxidation, and particle coagulation. Additional models include a piston-ring crevice model, the KH/RT spray breakup model, a droplet wall impingement model, a wall heat transfer model, and the RNG k-{var_epsilon} turbulence model. The Shell model was used to simulate the ignition process, and a laminar-and-turbulent characteristic time combustion model was used for the post-ignition combustion process. A low-load (IMEP=3 bar) operating condition was considered and the predicted in-cylinder pressures and heat release rates were compared with measurements. Predicted soot mass, soot particle size, soot number density distributions and other relevant quantities are presented and discussed. The effects of variable EGR rate (0-68%), injection pressure (600-1200 bar), and injection timing were studied. The predictions demonstrate that both EGR and retarded injection are beneficial for reducing NO{sub x} emissions, although the former has a more pronounced effect. Additionally, higher soot emissions are typically predicted for the higher EGR rates. However, when the EGR rate exceeds a critical value (over 65% in this study), the soot emissions decrease. Reduced soot emissions are also predicted when higher injection pressures or retarded injection timings are employed. The reduction in soot with retarded injection is less than what is observed experimentally, however.
Two-color resonant four-wave-mixing spectroscopy (TC-RFWM) is used to investigate ground-state energy transfer of hydroxyl radical in atmospheric-pressure flames. Two amplified distributed-feedback dye lasers produce 50-ps, nearly transform-limited, infrared (IR) and ultraviolet pulses. The infrared pump laser is tuned to individual rovibrational transitions of OH X {sup 2}{pi}{sub 3/2} (v{prime}=1, N{prime}) {l_arrow} X {sup 2}{pi}{sub 3/2} (v{double_prime}=0, N{double_prime}), and the ultraviolet pulse probes either the directly pumped or collisionally populated intermediate levels via A{sup 2}{Sigma}{sup +} (v*=1, N*) {l_arrow} X{sup 2}{pi}{sub 3/2}(v{prime}=1, N{prime}). By time-delaying the probe pulse with respect to the pump pulse, and appropriately constraining the polarizations of each of the four fields taking part in the wave-mixing process, we are able to independently and unambiguously measure the moments of the rotational angular momentum distribution in single rotational levels of the ground state. We present measurements of population, alignment, and orientation decay in X {sup 2}{pi}{sub 3/2} for several flame conditions. These experiments provide data necessary for the development of accurate models for diagnostic techniques using saturating laser pulses.
An alternative theory of solid mechanics, known as the peridynamic theory, formulates problems in terms of integral equations rather than partial differential equations. This theory assumes that particles in a continuum interact with each other across a finite distance, as in molecular dynamics. Damage is incorporated in the theory at the level of these two-particle interactions, so localization and fracture occur as a natural outgrowth of the equation of motion and constitutive models. A numerical method for solving dynamic problems within the peridynamic theory is described. Accuracy and numerical stability are discussed. Examples illustrate the properties of the method for modeling brittle dynamic crack growth.
Three methods that were used to measure the chemical changes associated with oxidative degradation of polymeric materials are presented. The first method is based on the nuclear activation of {sup 18}O in an elastomer that was thermally aged in an {sup 18}O{sub 2} atmosphere. Second, the alcohol groups in a thermally aged elastomer were derivatized with trifluoroacetic anhydride and their concentration measured via {sup 19}F NMR spectroscopy. Finally, a respirometer was used to directly measure the oxidative rates of a polyurethane foam as a function of aging temperature. The measurement of the oxidation rates enabled acceleration factors for oxidative degradation of these materials to be calculated.
An environmentally friendly method and materials study for desalinating inland brackish waters (i.e., coal bed methane produced waters) using a set of ion-exchange materials is presented. This desalination process effectively removes anions and cations in separate steps with minimal caustic waste generation. The anion-exchange material, hydrotalcite (HTC), exhibits an ion-exchange capacity (IEC) of {approx} 3 mequiv g{sup -1}. The cation-exchange material, an amorphous aluminosilicate permutite-like material, (Na{sub x+2y}Al{sub x}Si{sub 1-x}O{sub 2+y}), has an IEC of {approx}2.5 mequiv g{sup -1}. These ion-exchange materials were studied and optimized because of their specific ion-exchange capacity for the ions of interest and their ability to function in the temperature and pH regions necessary for cost and energy effectiveness. Room temperature, minimum pressure column studies (once-pass through) on simulant brackish water (total dissolved solids (TDS) = 2222 ppm) resulted in water containing TDS = 25 ppm. A second once-pass through column study on actual produced water (TDS = {approx}11,000) with a high carbonate concentration used an additional lime softening step and resulted in a decreased TDS of 600 ppm.
The Arsenic Water Technology Partnership program is a multi-year program funded by a congressional appropriation through the Department of Energy. The program is designed to move technologies from benchscale tests to field demonstrations. It will enable water utilities, particularly those serving small, rural communities and Indian tribes, to implement the most cost-effective solutions to their arsenic treatment needs. As part of the Arsenic Water Technology Partnership program, Sandia National Laboratories is carrying out field demonstration testing of innovative technologies that have the potential to substantially reduce the costs associated with arsenic removal from drinking water. The scope for this work includes: (1) Selection of sites and identification of technologies for pilot demonstrations; (2) Laboratory studies to develop rapid small-scale test methods; and (3) Pilot-scale studies at community sites involving side-by-side tests of innovative technologies. The goal of site selection is to identify sites that allow examination of treatment processes and systems under conditions that are relevant to different geochemical settings throughout the country. A number of candidate sites have been identified through reviews of groundwater quality databases, conference proceedings and discussions with state and local officials. These include sites in New Mexico, Arizona, Colorado, Oklahoma, Michigan, and California. Candidate technologies for the pilot tests are being reviewed through vendor forums, proof-of-principle benchscale studies managed by the American Water Works Association Research Foundation (AwwaRF) and the WERC design contest. The review considers as many potential technologies as possible and screens out unsuitable ones by considering data from past performance testing, expected costs, complexity of operation and maturity of the technology. The pilot test configurations will depend on the site-specific conditions such as access, power availability, waste disposal options and availability of permanent structures to house the test. Conducting pilot tests for media comparison at all sites in need of arsenic treatment would be extremely time consuming and costly. Laboratory studies are being conducted using rapid small-scale column tests (RSSCTs) to predict the performance of pilot-scale adsorption columns. RSSCTs are a rapid and inexpensive method of investigating innovative technologies while varying water quality and/or system design. RSSCTs are scaled down columns packed with smaller diameter adsorption media that receive higher hydraulic loading rates to significantly reduce the duration of experiments. Results for RSSCTs can be obtained in a matter of days to a few weeks, whereas pilot tests can take a number of months to over a year. In the pilot tests, the innovative technologies will be evaluated in terms of adsorptive capacity for arsenic; robustness of performance with respect to water quality parameters including pH, TDS, foulants such as Fe, Mn, silica, and organics, and other metals and radionuclides; and potentially deleterious effects on the water system such as pipe corrosion from low pH levels, fluoride removal, and generation of disinfection by-products. The new arsenic MCL will result in modification of many rural water systems that otherwise would not require treatment. Simultaneous improvement of water quality in systems that will require treatment for other contaminants such as uranium, radon and radium would be an added benefit of this program.
The use of rotational echo adiabatic passage double resonance (REAPDOR) solid-state nuclear magnetic resonance (NMR) to determine the site location of an adsorbed polar molecule in a zeolite cage is presented. Nitrogen-15 labeled ammonia is used as a probe molecule to characterize the initial adsorption site in 3A zeolite molecular sieves. The relative position of the ammonia adsorption site in the cage is determined by measuring the internuclear distance between the N on ammonia and both a Na cation site and an Al framework environment using {sup 15}N/{sup 23}Na and {sup 15}N/{sup 27}Al REAPDOR NMR experiments, respectively. The measured internuclear distances are similar to a specific ammonia adsorption site for the zeolite 4A ammonia sorption complex located using X-ray diffraction. Additional details regarding the ammonia hydrogen-bonding environment can be extracted from {sup 1}H/{sup 23}Na and {sup 1}H/{sup 27}Al REAPDOR NMR measurements.
Sandia National Laboratories has previously tested a capability to impose a 7.5 g-rms (30 g peak) radial vibration load up to 2 kHz on a 25 lb object with superimposed 50 g acceleration at its centrifuge facility. This was accomplished by attaching a 3,000 lb Unholtz-Dickie mechanical shaker at the end of the centrifuge arm to create a 'Vibrafuge'. However, the combination of non-radial vibration directions, and linear accelerations higher than 50g's are currently not possible because of the load capabilities of the shaker and the stresses on the internal shaker components due to the combined centrifuge acceleration. Therefore, a new technique using amplified piezo-electric actuators has been developed to surpass the limitations of the mechanical shaker system. They are lightweight, modular and would overcome several limitations presented by the current shaker. They are 'scalable', that is, adding more piezo-electric units in parallel or in series can support larger-weight test articles or displacement/frequency regimes. In addition, the units could be mounted on the centrifuge arm in various configurations to provide a variety of input directions. The design along with test results will be presented to demonstrate the capabilities and limitations of the new piezo-electric Vibrafuge.
The influence of thermal stratification on auto-ignition at constant volume and high pressure is studied by Direct Numerical Simulation (DNS) with complex H{sub 2}/air chemistry with a view to providing better understanding of combustion processes in homogeneous charge compression ignition engines. In particular the dependence of overall ignition progress on initial mixture conditions is determined. The propagation speed of ignition fronts that emanate from 'hot spots' given by a temperature spectrum is monitored by using the displacement velocity of a scalar that tracks the location of maximum heat release. The evolution of the front velocity is compared for different initial temperature distributions and the role of scalar dissipation of heat and mass is identified. It is observed that both deagrative as well as spontaneous ignition front propagation occur depending upon the local temperature gradient. It is found that the ratio of the instantaneous front speed to the deflagrative speed is a good measure of the local mode of propagation. This is verified by examining the energy and species balances. A parametric study in the amplitudes of the initial temperature fluctuation is performed and shows that this parameter has a significant influence on the observed combustion mode. Higher levels of stratification lead to more front-like structures. Predictions of the multi-zone model are presented and explained using the diagnostics developed.
The impacts of small niobium additions to processing, microstructure, and electrical properties in the Zr-rich lead zirconate titanate ceramics (PZT 95/5) were investigated. The influence of niobium content on dielectric responses and the characteristics of ferroelectric behaviors, as well as the relative phase stability and the hydrostatic pressure induced ferroelectric-to- antiferroelectric phase transformation are reported. Results indicate that increasing the niobium concentration in the solid solutions enhances densification, refines the microstructure, decreases dielectric constant and spontaneous polarization, and stabilizes the ferroelectric phase. The stabilization of ferroelectric phase with respect to the antiferroelectric phase near PZT 95/5 composition dramatically increases the pressure required for the ferroelectric-to-antiferroelectric phase transformation. These observations were correlated to the creation of A-site vacancies and a slight modification of the crystal structure. The importance of these composition-property relationships on device application will be presented.
Optical actuation of microelectromechanical systems (MEMS) is advantageous for applications for which electrical isolation is desired. Thirty-two polycrystalline silicon opto-thermal actuators, optically-powered MEMS thermal actuators, were designed, fabricated, and tested. The design of the opto-thermal actuators consists of a target for laser illumination suspended between angled legs that expand when heated, providing the displacement and force output. While the amount of displacement observed for the opto-thermal actuators was fairly uniform for the actuators, the amount of damage resulting from the laser heating ranged from essentially no damage to significant amounts of damage on the target. The likelihood of damage depended on the target design with two of the four target designs being more susceptible to damage. Failure analysis of damaged targets revealed the extent and depth of the damage.
A three-dimensional tungsten photonic crystal is thermally excited and shown to emit light at a narrow band, {lambda} = 3.3-4.25 {micro}m. The emission is experimentally observed to exceed that of the free-space Planck radiation over a wide temperature range, T = 475-850 K. it is proposed that an enhanced density of state associated with the propagating electromagnetic Bloch waves in the photonic crystal is responsible for this experimental finding.
This paper provides a brief overview of the fields of biological micro-electromechanical systems (bioMEMs) and associated nanobiotechnologies, collectively denoted as BioMicroNano. Although they are developing at a very rapid pace and still redefining themselves, several stabilized areas of research and development can be identified. Six major areas are delineated, and specific examples are discussed and illustrated. Various applications of the technologies are noted, and potential market sizes are compared.
Simplified models that are based on macroscopic force balances and droplet-geometry approximations are presented for predicting the onset of instability leading to removal of water droplets at the gas diffusion layer (GDL)/gas flow channel (GFC) interface. Visualization experiments are carried out to observe the formation, growth, and removal or instability of the water droplets at the GDL/GFC interface of a simulated polymer electrolyte fuel cell cathode. Droplet-instability diagrams or windows computed by the simplified models are compared with those measured experimentally, and good agreement is obtained. Two-dimensional flow simulations employing the finite element method coupled with an arbitrary Lagrangian-Eulerian formulation for determining the liquid/gas interface position are also performed to assess the simplified cylindrical-droplet model. Necessary conditions for preventing fully grown droplets from lodging in the flow channel are derived using the simplified models. It is found that droplet removal can be enhanced by increasing flow channel length or mean gas flow velocity, decreasing channel height or contact angle hysteresis, or making the GDL/GFC interface more hydrophobic.
Concerns over the illicit trafficking of radiological and nuclear materials were focused originally on the lack of security and accountability of such material throughout the former Soviet states. This is primarily attributed to the frequency of events that have occurred involving the theft and trafficking of critical material components that could be used to construct a Radiological Dispersal Device (RDD) or even a rudimentary nuclear device. However, with the continued expansion of nuclear technology and the deployment of a global nuclear fuel cycle these materials have become increasingly prevalent, affording a more diverse inventory of dangerous materials and dual-use items. To further complicate the matter, the list of nuclear consumers has grown to include: (1) Nation-states that have gone beyond the IAEA agreed framework and additional protocols concerning multiple nuclear fuel cycles and processes that reuse the fuel through reprocessing to exploit technologies previously confined to the more industrialized world; (2) Terrorist organizations seeking to acquire nuclear and radiological material due to the potential devastation and psychological effect of their use; (3) Organized crime, which has discovered a lucrative market in trafficking of illicit material to international actors and/or countries; and (4) Amateur smugglers trying to feed their families in a post-Soviet era. An initial look at trafficking trends of this type seems scattered and erratic, localized primarily to a select group of countries. This is not necessarily the case. The success with which other contraband has been smuggled throughout the world suggests that nuclear trafficking may be carried out with relative ease along the same routes by the same criminals or criminal organizations. Because of the inordinately high threat posed by terrorist or extremist groups acquiring the ingredients for unconventional weapons, it is necessary that illicit trafficking of these materials be better understood as to prepare for the sustained global development of the nuclear fuel cycle. Conversely, modeling and analyses of this activity must not be limited in their scope to loosely organized criminal smuggling, but address the problem as a commercial, industrial project for the covert development of nuclear technologies and unconventional weapon development.
The nation's electric power sector is highly interdependent with the economic sectors it serves; electric power needs are driven by economic activity while the economy itself depends on reliable and sustainable electric power. To advance higher level understandings of the vulnerabilities that result from these interdependencies and to identify the loss prevention and loss mitigation policies that best serve the nation, the National Infrastructure Simulation and Analysis Center is developing and using N-ABLE{trademark}, an agent-based microeconomic framework and simulation tool that models these interdependencies at the level of collections of individual economic firms. Current projects that capture components of these electric power and economic sector interdependencies illustrate some of the public policy issues that should be addressed for combined power sector reliability and national economic security.
The thermal interdiffusion of AlSb/GaSb multiquantum wells was measured and the intrinsic diffusivities of Al and Ga determined over a temperature range of 823-948 K for 30-9000 s. The 77-K photoluminescence (PL) was used to monitor the extent of interdiffusion through the shifts in the superlattice luminescence peaks. The chemical diffusion coefficient was quantitatively determined by fitting the observed PL peak shifts to the solution of the Schroedinger equation, using a potential derived from the solution of the diffusion equation. The value of the interdiffusion coefficient ranged from 5.2 x 10{sup -4} to 0.06 nm{sup 2}/s over the conditions studied and was characterized by an activation energy of 3.0 {+-} 0.1 eV. The intrinsic diffusion coefficients for Al and Ga were also determined with higher values for Al than for Ga, described by activation energies of 2.8 {+-} 0.4 and 1.1 {+-} 0.1 eV, respectively.
A number of industrial combustion systems are adopting oxygen-enhanced firing to improve heat transfer characteristics and reduce emissions. The exhaust gas from these systems is dominated by H2O and CO2 and therefore has substantially different gas properties from traditional combustion exhaust. In the past, laser-induced breakdown spectroscopy (LIBS) has been successfully used for the evaluation of alkali aerosol concentrations in air-based combustion systems. This paper presents results of LIBS measurements of alkali concentrations in a laboratory calibration setup and in an oxygen/natural gas container glass furnace. It shows how both gas conditions (composition and temperature) and the molecular form of the alkali species affect the LIBS signals. The paper proposes strategies for mitigating these effects in future applications of LIBS in oxygen-enhanced combustion systems.
The dimensionless extinction coefficient (K{sub e}) of soot must be known to quantify laser extinction measurements of soot concentration and to predict optical attenuation through smoke clouds. Previous investigations have measured K{sub e} for post-flame soot emitted from laminar and turbulent diffusion flames and smoking laminar premixed flames. This paper presents the first measurements of soot K{sub e} from within laminar diffusion flames, using a small extractive probe to withdraw the soot from the flame. To measure K{sub e}, two laser sources (635 nm and 1310 nm) were coupled to a transmission cell, followed by gravimetric sampling. Coannular diffusion flames of methane, ethylene and nitrogen-diluted kerosene burning in air were studied, together with slot flames of methane and ethylene. K{sub e} was measured at the radial location of maximum soot volume fraction at several heights for each flame. Results for K{sub e} at both 635 nm and 1310 nm for ethylene and kerosene coannular flames were in the range of 9-10, consistent with the results from previous studies of post-flame soot. The ethylene slot flame and the methane flames have lower K{sub e} values, in some cases as low as 2.0. These lower values of K{sub e} are found to result from the contributions of (a) the condensation of PAH species during the sampling of soot, (b) the wavelength-dependent absorptivity of soot precursor particles, and, in the case of methane, (c) the negligible contribution of soot scattering to the extinction coefficient. RDG calculations of soot scattering, in combination with the measured K{sub e} values, imply that the soot refractive index is in the vicinity of 1.75-1.03i at 635 nm.
A severe fire and explosion occurred at a propane storage yard in Truth or Consequences, N.M., when a truck ran into the pumping and plumbing system beneath a large propane tank. The storage tank emptied when the liquid-phase excess flow valve tore out of the tank. The ensuing fire engulfed several propane delivery trucks, causing one of them to explode. A series of elevated-temperature stress-rupture tears developed along the top of a 9800 L (2600 gal) truck-mounted tank as it was heated by the fire. Unstable fracture then occurred suddenly along the length of the tank and around both end caps, along the girth welds connecting the end caps to the center portion of the tank. The remaining contents of the tank were suddenly released, aerosolized, and combusted, creating a powerful boiling liquid expanding vapor explosion (BLEVE). Based on metallography of the tank pieces, the approximate tank temperature at the onset of the BLEVE was determined. Metallurgical analysis of the ruptured tank also permitted several hypotheses regarding BLEVE mechanisms to be evaluated. Suggestions are made for additional work that could provide improved predictive capabilities regarding BLEVEs and for methods to decrease the susceptibility of propane tanks to BLEVEs.
43rd AIAA Aerospace Sciences Meeting and Exhibit - Meeting Papers
Chen, Guoquan; Collis, S.S.
A new multidomain/multiphysics computational framework for optimal control of aeroacoustic noise has been developed based on a near-field compressible Navier-Stokes solver coupled with a far-field linearized Euler solver both based on a discontinuous Galerkin formulation. In this approach, the coupling of near- and far-field domains is achieved by weakly enforcing continuity of normal fluxes across a coupling surface that encloses all nonlinearities and noise sources. For optimal control, gradient information is obtained by the solution of an appropriate adjoint problem that involves the propagation of adjoint information from the far-field to the near-field. This computational framework has been successfully applied to study optimal boundary-control of blade-vortex interaction, which is a significant noise source for helicopters on approach to landing. In the model-problem presented here, the noise propogated toward the ground is reduced by 12dB.
An electrically programmable surface acoustic wave (SAW) correlator was recently completed from design through small scale production in support of low power space-based communications for NASA. Three different versions of this RF microsystem were built to satisfy design requirements and overcome packaging and system reliability related issues. Flip-chip packaging and conventional thick film hybrid assembly techniques are compared in the fabrication of this microsystem.
The apparent metal-insulator transition is observed in a high-quality two-dimensional electron system (2DES) in the strained Si quantum well of a Si/Si{sub 1-x}Ge{sub x} heterostructure with mobility {mu} = 1.9 x 10{sup 5} cm{sup 2}/V s at density n = 1.45 x 10{sup 11} cm{sup -2}. The critical density, at which the thermal coefficient of low T resistivity changes sign, is -0.32 x 10{sup 11} cm{sup -2}, a very low value obtained in Si-based 2D systems. The in-plane magnetoresistivity {rho}(B{sub ip}) was measured in the density range, 0.35 x 10{sup 11} < n < 1.45 x 10{sup 11} cm{sup -2}, where the 2DES shows the metallic-like behavior. It first increases and then saturates to a finite value {rho}(B{sub c}) for B{sub ip}>B{sub c} , with B{sub c} the full spin-polarization field. Surprisingly, {rho}(B{sub c})/{rho}(0)-1.8 for all the densities, even down to n = 0.35 x 10{sup 11} cm{sup -2}, only 10% higher than n{sub c}. This is different from that in clean Si metal-oxide-semiconductor field-effect transistors, where the enhancement is strongly density dependent and {rho}(B{sub c})/{rho}(0) appears to diverge as n {yields} n{sub c}. Finally, we show that in the fully spin-polarized regime, dependent on the 2DES density, the temperature dependence of {rho}(B{sub ip}) can be either metallic-like or insulating.
This article describes recent improvements in mapping a highly representative set of the world-wide scientific literature. The process described in this article extends existing work in this area in three major ways. First, we argue that a separate structural analysis of current literature vs. reference literature is required for R&D planning. Second, visualization software is used to improve coverage of the literature while maintaining structural integrity. Third, quantitative techniques for measuring the structural integrity of a map are introduced. Maps with high structural integrity, covering far more of the available literature, are presented.
A dynamic hohlraum is created when an annular z-pinch plasma implodes onto a cylindrical 0.014 g/cc 6-mm-diameter CH{sub 2} foam. The impact launches a radiating shock that propagates toward the axis at {approx}350 {micro}m/ns. The radiation trapped by the tungsten z-pinch plasma forms a {approx}200 eV hohlraum that provides X-rays for indirect drive inertial confinement fusion capsule implosion experiments. We are developing the ability to diagnose the hohlraum interior using emission and absorption spectroscopy of Si atoms added as a tracer to the central portion of the foam. Time- and space-resolved Si spectra are recorded with an elliptical crystal spectrometer viewing the cylindrical hohlraum end-on. A rectangular aperture at the end of the hohlraum restricts the field of view so that the 1D spectrometer resolution corresponds approximately to the hohlraum radial direction. This enables distinguishing between spectra from the unshocked radiation-heated foam and from the shocked foam. Typical spectral lines observed include the Si Ly{alpha} with its He-like satellites and the He-like resonance sequence including He{alpha}, He{beta}, and He{gamma}, along with some of their associated Li-like satellites. Work is in progress to infer the hohlraum conditions using collisional-radiative modeling that accounts for the radiation environment and includes both opacity effects and detailed Stark broadening calculations. These 6-mm-scale radiation-heated plasmas might eventually also prove suitable for testing Stark broadening line profile calculations or for opacity measurements.
Fluorescence correlation spectroscopy (FCS) is used to examine mobility of labeled probes at specific sites in supported bilayers consisting of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) lipid domains in 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC). Those sites are mapped beforehand with simultaneous atomic force microscopy and submicron confocal fluorescence imaging, allowing characterization of probe partitioning between gel DPPC and disordered liquid DOPC domains with corresponding topography of domain structure. We thus examine the relative partitioning and mobility in gel and disordered liquid phases for headgroup- and tailgroup-labeled GM1 ganglioside probes and for headgroup- and tailgroup-labeled phospholipid probes. For the GM1 probes, large differences in mobility between fluid and gel domains are observed; whereas unexpected mobility is observed in submicron gel domains for the phospholipid probes. We attribute the latter to domain heterogeneities that could be induced by the probe. Furthermore, fits to the FCS data for the phospholipid probes in the DOPC fluid phase require two components (fast and slow). Although proximity to the glass substrate may be a factor, local distortion of the probe by the fluorophore could also be important. Overall, we observe nonideal aspects of phospholipid probe mobility and partitioning that may not be restricted to supported bilayers.
We consider linear systems arising from the use of the finite element method for solving a certain class of linear elliptic problems. Our main result is that these linear systems, which are symmetric and positive semidefinite, are well approximated by symmetric diagonally dominant matrices. Our framework for defining matrix approximation is support theory. Significant graph theoretic work has already been developed in the support framework for preconditioners in the diagonally dominant case, and in particular it is known that such systems can be solved with iterative methods in nearly linear time. Thus, our approximation result implies that these graph theoretic techniques can also solve a class of finite element problems in nearly linear time. We show that the quality of our approximation, which controls the number of iterations in the preconditioned iterative solver, depends primarily on a mesh quality measure but not on the problem size or shape of the domain.
The authors explore various possible approaches for generating lowest order and higher order bases for modeling surface currents and their divergence for moment method application to integral equations. The bases developed are defined on curved triangular and quadrilateral elements. All the bases are conveniently defined in parent element coordinates, and each expansion function spans one or two patches.
This article examines how well one can predict the importance of a current paper (a paper that is recently published in the literature). We look at three factors--journal importance, reference importance and author reputation. Citation-based measures of importance are used for all variables. We find that journal importance is the best predictor (explaining 22.3% out of a potential 29.1% of the variance in the data), and that this correlation value varies significantly by discipline. Journal importance is a better predictor of citation in Computer Science than in any other discipline. While the finding supports the present policy of using journal impact statistics as a surrogate for the importance of current papers, it calls into question the present policy of equally weighting current documents in text-based analyses. We suggest that future researchers take into account the expected importance of a document when attempting to describe the cognitive structure of a field.
American Rock Mechanics Association - 40th US Rock Mechanics Symposium, ALASKA ROCKS 2005: Rock Mechanics for Energy, Mineral and Infrastructure Development in the Northern Regions
Bauer, Stephen; Munson, D.E.; Hardy, M.P.; Barrix, J.; McGunegle, B.
Stress measurements have been obtained from within the Norton Mine in support of site characterization activities intended to determine the in situ stress field around the mine. These results together with other measurements in the area permit an estimate of the principal stresses at the mine. Based on the most recent measurements, the maximum (σHmax) and minimum (σHmin) stresses acting in the horizontal plane are oriented nearly east-west and north-south, respectively, and their magnitudes are 5330 psi and 4100 psi, respectively. These values are expected to be essentially uniform within a few hundred feet vertically above and below the mine elevation. The stress acting in the vertical direction has a magnitude of 3270 psi at the mine level. This measured vertical stress is related to the overburden weight according to σv=1.26ρgh (where ρ is the overburden density, g acceleration of gravity, and h overburden depth). The measured vertical stress exceeds the stress calculated from overburden weight by a factor of 1.26. These in situ stresses are assumed to be principal stresses and, as a result, the vertical stress is the minimum principal stress. These measurements are generally consistent in magnitude and direction with two other much older sets of measurements taken in the mine and they are consistent with the east-west trend of the regional in situ principal stress direction. The average of all three sets of measurements, recent and old, in the mine give a maximum horizontal stress of 6110 psi, a minimum horizontal stress of 3630, and a vertical stress of 3030 psi. The directions of the mine excavation development, which normally are oriented according to the principal stresses, are also consistent with the current and past measurements.
Proceedings of the ASME/Pacific Rim Technical Conference and Exhibition on Integration and Packaging of MEMS, NEMS, and Electronic Systems: Advances in Electronic Packaging 2005
