10th International Conference on Probabilistic Safety Assessment and Management 2010, PSAM 2010
Dang, Vinh N.; Forester, John A.; Mosleh, Ali
The Office of Nuclear Regulatory Research (RES) of the U.S. Nuclear Regulatory Commission is sponsoring work in response to a Staff Requirements Memorandum (SRM) directing an effort to establish a single human reliability analysis (HRA) method for the agency or guidance for the use of multiple methods. One motivation is the variability in Human Failure Event (HFE) probabilities estimated by different analysts and methods. This work considers that a reduction of the variability in the HRA quantification outputs must address three sources: differences in the scope and implementation of qualitative analysis, the qualitative output-quantitative input interface, and the diversity of algorithms for estimating failure probabilities from these inputs. Two companion papers (Mosleh et al. and Hendrickson et al.) describe a proposed qualitative analysis approach The development of the corresponding quantification approach considers a number of alternatives including a module-based hybrid method and a data-driven quantification scheme. This paper presents on-going work and the views of the contributors.
Cyber security analysis tools are necessary to evaluate the security, reliability, and resilience of networked information systems against cyber attack. It is common practice in modern cyber security analysis to separately utilize real systems computers, routers, switches, firewalls, computer emulations (e.g., virtual machines) and simulation models to analyze the interplay between cyber threats and safeguards. In contrast, Sandia National Laboratories has developed new methods to combine these evaluation platforms into a cyber Live, Virtual, and Constructive (LVC) testbed. The combination of real, emulated, and simulated components enables the analysis of security features and components of a networked information system. When performing cyber security analysis on a target system, it is critical to represent realistically the subject security components in high fidelity. In some experiments, the security component may be the actual hardware and software with all the surrounding components represented in simulation or with surrogate devices. Sandia National Laboratories has developed a cyber LVC testbed that combines modeling and simulation capabilities with virtual machines and real devices to represent, in varying fidelity, secure networked information system architectures and devices. Using this capability, secure networked information system architectures can be represented in our testbed on a single computing platform. This provides an "experiment-in-a-box" capability. The result is rapidly produced, large scale, relatively low-cost, multi-fidelity representations of networked information systems. These representations enable analysts to quickly investigate cyber threats and test protection approaches and configurations.
The launch of nuclear materials requires special care to minimize the risk of adverse effects to human health and the environment. This paper describes the special sources of risk that are inherent to the launch of radioactive materials and provides insights into the analysis and control of these risks that have been gained through the experience of previous US launches. Historically, launch safety has been achieved by eliminating, to the greatest degree possible, the potential for energetic insults to affect the radioactive material. For those insults that cannot be precluded, designers minimize the likelihood, magnitude and duration of their interaction with the material. Finally, when a radioactive release cannot be precluded, designers limit the magnitude and spatial extent of its dispersal.
Phillips, Stan D.; Moen, Kurt A.; Najafizadeh, Laleh; Diestelhorst, Ryan M.; Sutton, Akil K.; Cressler, John D.; Vizkelethy, Gyorgy; Dodd, Paul E.; Marshall, Paul W.
Human Reliability Analysis (HRA) methods have been developed primarily to provide information for use in probabilistic risk assessments analyzing nuclear power plant (NPP) operations. Despite this historical focus on the control room, there has been growing interest in applying HRA methods to other NPP activities such as dry cask storage operations (DCSOs) in which spent fuel is transferred into dry cask storage systems. This paper describes a successful application of aspects of the "A Technique for Human Event Analysis" (ATHEANA) HRA approach [1, 2] in performing qualitative HRA activities that generated insights on the potential for dropping a spent fuel cask during DCSOs. This paper provides a description of the process followed during the analysis, a description of the human failure event (HFE) scenario groupings, discussion of inferred human performance vulnerabilities, a detailed examination of one HFE scenario and illustrative approaches for avoiding or mitigating human performance vulnerabilities that may contribute to dropping a spent fuel cask.
39th ASES National Solar Conference 2010, SOLAR 2010
Gupta, Vipin P.; Boudra, Will; Kuszmaul, Scott S.; Rosenthal, Andrew; Cisneros, Gaby; Merrigan, Tim; Miller, Ryan; Dominick, Jeff
In May 2007, Forest City Military Communities won a US Department of Energy Solar America Showcase Award. As part of this award, executives and staff from Forest City Military Communities worked side-by-side with a DOE technical assistance team to overcome technical obstacles encountered by this large-scale real estate developer and manager. This paper describes the solar technical assistance that was provided and the key solar experiences acquired by Forest City Military Communities over an 18 month period.
The Department of Energy's 2008 Yucca Mountain Performance Assessment represents the culmination of more than two decades of analyses of post-closure repository performance in support of programmatic decision making for the proposed Yucca Mountain repository. The 2008 performance assessment summarizes the estimated long-term risks to the health and safety of the public resulting from disposal of spent nuclear fuel and high-level radioactive waste in the proposed Yucca Mountain repository. The standards at 10 CFR Part 63 request several numerical estimates quantifying performance of the repository over time. This paper summarizes the key quantitative results from the performance assessment and presents uncertainty and sensitivity analyses for these results.
We present an RDFS closure algorithm, specifically designed and implemented on the Cray XMT supercomputer, that obtains inference rates of 13 million inferences per second on the largest system configuration we used. The Cray XMT, with its large global memory (4TB for our experiments), permits the construction of a conceptually straightforward algorithm, fundamentally a series of operations on a shared hash table. Each thread is given a partition of triple data to process, a dedicated copy of the ontology to apply to the data, and a reference to the hash table into which it inserts inferred triples. The global nature of the hash table allows the algorithm to avoid a common obstacle for distributed memory machines: the creation of duplicate triples. On LUBM data sets ranging between 1.3 billion and 5.3 billion triples, we obtain nearly linear speedup except for two portions: file I/O, which can be ameliorated with the additional service nodes, and data structure initialization, which requires nearly constant time for runs involving 32 processors or more.
Innovative energy system optimization models are deployed to evaluate novel fuel cell system (FCS) operating strategies, not typically pursued by commercial industry. Most FCS today are installed according to a 'business-as-usual' approach: (1) stand-alone (unconnected to district heating networks and low-voltage electricity distribution lines), (2) not load following (not producing output equivalent to the instantaneous electrical or thermal demand of surrounding buildings), (3) employing a fairly fixed heat-to-power ratio (producing heat and electricity in a relatively constant ratio to each other), and (4) producing only electricity and no recoverable heat. By contrast, models discussed here consider novel approaches as well. Novel approaches include (1) networking (connecting FCSs to electrical and/or thermal networks), (2) load following (having FCSs produce only the instantaneous electricity or heat demanded by surrounding buildings), (3) employing a variable heat-to-power ratio (such that FCS can vary the ratio of heat and electricity they produce), (4) co-generation (combining the production of electricity and recoverable heat), (5) permutations of these together, and (6) permutations of these combined with more 'business-as-usual' approaches. The detailed assumptions and methods behind these models are described in Part I of this article pair.
The U.S. Strategic Petroleum Reserve stores crude oil in 62 solution-mined caverns in salt domes located in Texas and Louisiana. Historically, three-dimensional geomechanical simulations of the behavior of the caverns have been performed using a power law creep model. Using this method, and calibrating the creep coefficient to field data such as cavern closure and surface subsidence, has produced varying degrees of agreement with observed phenomena. However, as new salt dome locations are considered for oil storage facilities, pre-construction geomechanical analyses are required that need site-specific parameters developed from laboratory data obtained from core samples. The multi-mechanism deformation (M-D) model is a rigorous mathematical description of both transient and steady-state creep phenomena. Recent enhancements to the numerical integration algorithm within the model have created a more numerically stable implementation of the M-D model. This report presents computational analyses to compare the results of predictions of the geomechanical behavior at the West Hackberry SPR site using both models. The recently-published results using the power law creep model produced excellent agreement with an extensive set of field data. The M-D model results show similar agreement using parameters developed directly from laboratory data. It is also used to predict the behavior for the construction and operation of oil storage caverns at a new site, to identify potential problems before a final cavern layout is designed. Copyright 2010 ARMA, American Rock Mechanics Association.
Here we demonstrate the suitability of robust nucleic acid affinity reagents in an integrated point-of-care diagnostic platform for monitoring proteomic biomarkers indicative of astronaut health in spaceflight applications. A model thioaptamer[1] targeting nuclear factor-kappa B (NF-KB) is evaluated in an on-chip electrophoretic gel-shift assay for human serum. Key steps of i) mixing sample with the aptamer, ii) buffer exchange, and iii) preconcentration of sample were successfully integrated upstream of fluorescence-based detection. Challenges due to i) nonspecific interactions with serum, and ii) preconcentration at a nanoporous membrane are discussed and successfully resolved to yield a robust, rapid, and fully-integrated diagnostic system.
We present a method for counting white blood cells that is uniquely compatible with centrifugation based microfluidics. Blood is deposited on top of one or more layers of density media within a microfluidic disk. Spinning the disk causes the cell populations within whole blood to settle through the media, reaching an equilibrium based on the density of each cell type. Separation and fluorescence measurement of cell types stained with a DNA dye is demonstrated using this technique. The integrated signal from bands of fluorescent microspheres is shown to be proportional to their initial concentration in suspension.
We present a platform that combines patterned photopolymerized polymer monoliths with living radical polymerization (LRP) to develop a low cost microfluidic based immunoassay capable of sensitive (low to sub pM) and rapid (<30 minute) detection of protein in 100 μL sample. The introduction of LRP functionality to the porous monolith allows one step grafting of functionalized affinity probes from the monolith surface while the composition of the hydrophilic graft chain reduces non-specific interactions and helps to significantly improve the limit of detection.
We have designed, fabricated, and characterized a digital microfluidic (DMF) platform to function as a central hub for interfacing multiple lab-on-a-chip sample processing modules towards automating the preparation of clinically-derived DNA samples for ultrahigh throughput sequencing (UHTS). The platform enables plug-and-play installation of a two-plate DMF device with consistent spacing, offers flexible connectivity for transferring samples between modules, and uses an intuitive programmable interface to control droplet/electrode actuations. Additionally, the hub platform uses transparent indium-tin oxide (ITO) electrodes to allow complete top and bottom optical access to the droplets on the DMF array, providing additional flexibility for various detection schemes.
We report on advancements of our microscale isoelectric fractionation (μIEFr) methodology for fast on-chip separation and concentration of proteins based on their isoelectric points (pI). We establish that proteins can be fractionated depending on posttranslational modifications into different pH specific bins, from where they can be efficiently transferred to downstream membranes for additional processing and analysis. This technology can enable on-chip multidimensional glycoproteomics analysis, as a new approach to expedite biomarker identification and verification.
We report on a scalable electrostatic process to transfer epitaxial graphene to arbitrary glass substrates, including Pyrex and Zerodur. This transfer process could enable wafer-level integration of graphene with structured and electronically-active substrates such as MEMS and CMOS. We will describe the electrostatic transfer method and will compare the properties of the transferred graphene with nominally-equivalent 'as-grown' epitaxial graphene on SiC. The electronic properties of the graphene will be measured using magnetoresistive, four-probe, and graphene field effect transistor geometries [1]. To begin, high-quality epitaxial graphene (mobility 14,000 cm2/Vs and domains >100 {micro}m2) is grown on SiC in an argon-mediated environment [2,3]. The electrostatic transfer then takes place through the application of a large electric field between the donor graphene sample (anode) and the heated acceptor glass substrate (cathode). Using this electrostatic technique, both patterned few-layer graphene from SiC(000-1) and chip-scale monolayer graphene from SiC(0001) are transferred to Pyrex and Zerodur substrates. Subsequent examination of the transferred graphene by Raman spectroscopy confirms that the graphene can be transferred without inducing defects. Furthermore, the strain inherent in epitaxial graphene on SiC(0001) is found to be partially relaxed after the transfer to the glass substrates.
We report on the novel room temperature method of synthesizing advanced nuclear fuels; a method that virtually eliminates any volatility of components. This process uses radiolysis to form stable nanoparticle (NP) nuclear transuranic (TRU) fuel surrogates and in-situ heated stage TEM to sinter the NPs. The radiolysis is performed at Sandia's Gamma Irradiation Facility (GIF) 60Co source (3 x 10{sup 6} rad/hr). Using this method, sufficient quantities of fuels for research purposes can be produced for accelerated advanced nuclear fuel development. We are focused on both metallic and oxide alloy nanoparticles of varying compositions, in particular d-U, d-U/La alloys and d-UO2 NPs. We present detailed descriptions of the synthesis procedures, the characterization of the NPs, the sintering of the NPs, and their stability with temperature. We have employed UV-vis, HRTEM, HAADF-STEM imaging, single particle EDX and EFTEM mapping characterization techniques to confirm the composition and alloying of these NPs.
This presentation on wind energy discusses: (1) current industry status; (2) turbine technologies; (3) assessment and siting; and (4) grid integration. There are no fundamental technical barriers to the integration of 20% wind energy into the nation's electrical system, but there needs to be a continuing evolution of transmission planning and system operation policy and market development for this to be most economically achieved.
Magnetic Liner Inertial Fusion (MagLIF) [S. A. Slutz, et al., Phys. Plasmas 17 056303 (2010)] is a promising new concept for achieving >100 kJ of fusion yield on Z. The greatest threat to this concept is the Magneto-Rayleigh-Taylor (MRT) instability. Thus an experimental campaign has been initiated to study MRT growth in fast-imploding (<100 ns) cylindrical liners. The first sets of experiments studied aluminum liner implosions with prescribed sinusoidal perturbations (see talk by D. Sinars). By contrast, this poster presents results from the latest sets of experiments that used unperturbed beryllium (Be) liners. The purpose for using Be is that we are able to radiograph 'through' the liner using the 6-keV photons produced by the Z-Beamlet backlighting system. This has enabled us to obtain time-resolved measurements of the imploding liner's density as a function of both axial and radial location throughout the field of view. This data is allowing us to evaluate the integrity of the inside (fuel-confining) surface of the imploding liner as it approaches stagnation.
We present the preliminary design of a Z experiment intended to observe the growth of several hydrodynamic instabilities (RT, RM, and KH) in a high-energy-density plasma. These experiments rely on the Z-machine's unique ability to launch cm-sized slabs of cold material (known as flyer plates) to velocities of several times 10 km/s. During the proposed experiment, the flyer plate will impact a cm-sized target with an embedded interface that has a prescribed sinusoidal perturbation. The flyer plate will generate a strong shock that propagates into the target and later initiates unstable growth of the perturbation. The goal of the experiment is to observe the perturbation at various stages of its evolution as it transitions from linear to non-linear growth, and finally to a fully turbulent state.
Laser-accelerated proton beams can be used in a variety of applications, e.g. ultrafast radiography of dense objects or strong electromagnetic fields. Therefore high energies of tens of MeV are required. We report on proton-acceleration experiments with a 150 TW laser system using mm-sized thin foils and mass-reduced targets of various thicknesses. Thin- foil targets yielded maximum energies of 50 MeV. A further reduction of the target dimensions from mm-size to 250 x 250 x 25 microns increased the maximum proton energy to >65 MeV, which is comparable to proton energies measured only at higher-energy, Petawatt-class laser systems. The dependence of the maximum energy on target dimensions was investigated, and differences between mm-sized thin foils and mass-reduced targets will be reported.
The objectives of this presentation are: (1) To determine if healthcare settings serve as intensive transmission environments for influenza epidemics, increasing effects on communities; (2) To determine which mitigation strategies are best for use in healthcare settings and in communities to limit influenza epidemic effects; and (3) To determine which mitigation strategies are best to prevent illness in healthcare workers.
Cigarette smoking presented the most significant public health challenge in the United States in the 20th Century and remains the single most preventable cause of morbidity and mortality in this country. A number of System Dynamics models exist that inform tobacco control policies. We reviewed them and discuss their contributions. We developed a theory of the societal lifecycle of smoking, using a parsimonious set of feedback loops to capture historical trends and explore future scenarios. Previous work did not explain the long-term historical patterns of smoking behaviors. Much of it used stock-and-flow to represent the decline in prevalence in the recent past. With noted exceptions, information feedbacks were not embedded in these models. We present and discuss our feedback-rich conceptual model and illustrate the results of a series of simulations. A formal analysis shows phenomena composed of different phases of behavior with specific dominant feedbacks associated with each phase. We discuss the implications of our society's current phase, and conclude with simulations of what-if scenarios. Because System Dynamics models must contain information feedback to be able to anticipate tipping points and to help identify policies that exploit leverage in a complex system, we expanded this body of work to provide an endogenous representation of the century-long societal lifecycle of smoking.
Negative bias temperature instability is an issue of critical importance as the space electronics industry evolves because it may dominate the reliability lifetime. Understanding its physical origin is therefore essential in determining how best to search for methods of mitigation. It has been suggested that the magnitude of the effect is strongly dependent on circuit operation conditions (static or dynamic modes). In the present work, we examine the time constants related to the charging and recovery of trapped charged induced by NBTI in HfSiON gate dielectric devices. In previous work, we avoided the issue of charge relaxation during acquisition of the I{sub ds}(V{sub gs}) curve by invoking a continuous stressing technique whereby {Delta}V{sub th} was extracted from a series of single point I{sub ds} measurements. This method relied heavily on determination of the initial value of the source-drain current (I{sub ds}{sup o}) prior to application of gate-source stress. In the present work we have used a new pulsed measurement system (Keithley SCS 4200-PIV) which not only removes this uncertainty but also permits dynamic measurements in which devices are AC stressed (Fig. 1a) or subjected to cycles of continued DC stresses followed by relaxation (Fig. 1b). We can now examine the charging and recovery characteristics of NBTI with higher precision than previously possible. We have performed NBTI stress experiments at room temperature on p-channel MOSFETs made with HfSiON gate dielectrics. In all cases the devices were stressed in the linear regime with V{sub ds}=-0.1V. We have defined two separate waveforms/pulse trains as illustrated in Fig 1. These were applied to the gate of the MOSFET. Firstly we examined the charging characteristics by applying an AC stress at 2.5MHz or 10Hz for different times. For a 50% duty cycle this corresponded to V{sub gs} = - 2V pulses for 200ns or 500ms followed by V{sub gs} = 0V pulses for 200ns or 500ms recovery respectively. In between 'bursts' of AC stress cycles, the I{sub ds}(V{sub gs}) characteristic in the range (-0.6V, -1.3V) was measured in 10.2 {micro}s. V{sub th} was extracted directly from this curve, or from a single I{sub ds} point normalized to the initial I{sub ds}{sup o} using our previous method. The resulting I{sub ds}/I{sub ds}{sup o} curves are compared; in Fig 2, the continuous stress results are included. In the second method, we examined the recovery dynamic by holding V{sub gs} = 0V for a finite amount of time (range 100 ns to 100 ms) following stress at V{sub gs} = - 2V for various times. In Fig 3 we compare |{Delta}V{sub th}(t)| results for recovery times of 100ms, 1ms, 100{micro}s, 50{micro}s, 25{micro}s, 10{micro}s, 100ns, and DC (i.e. no recovery) The data in Fig 2 shows that with a high frequency stress (2.5MHz) devices undergo significantly less (but finite) current degradation than devices stressed at 10Hz. This appears to be limited by charging and not by recovery. Fig 3 supports this hypothesis since for 100ns recovery periods, only a small percentage of the trapped charge relaxes. Detailed explanation of these experiments will be presented at the conference.
Loki-Infect 3 is a desktop application intended for use by community-level decision makers. It allows rapid construction of small-scale studies of emerging or hypothetical infectious diseases in their communities and evaluation of the potential effectiveness of various containment strategies. It was designed with an emphasis on modularity, portability, and ease of use. Our goal is to make this program freely available to community workers across the world.
Injection of CO2 into formations containing brine is proposed as a long-term sequestration solution. A significant obstacle to sequestration performance is the presence of existing wells providing a transport pathway out of the sequestration formation. To understand how heterogeneity impacts the leakage rate, we employ two dimensional models of the CO2 injection process into a sandstone aquifer with shale inclusions to examine the parameters controlling release through an existing well. This scenario is modeled as a constant-rate injection of super-critical CO2 into the existing formation where buoyancy effects, relative permeabilities, and capillary pressures are employed. Three geologic controls are considered: stratigraphic dip angle, shale inclusion size and shale fraction. In this study, we examine the impact of heterogeneity on the amount and timing of CO2 released through a leaky well. Sensitivity analysis is performed to classify how various geologic controls influence CO2 loss. A 'Design of Experiments' approach is used to identify the most important parameters and combinations of parameters to control CO2 migration while making efficient use of a limited number of computations. Results are used to construct a low-dimensional description of the transport scenario. The goal of this exploration is to develop a small set of parametric descriptors that can be generalized to similar scenarios. Results of this work will allow for estimation of the amount of CO2 that will be lost for a given scenario prior to commencing injection. Additionally, two-dimensional and three-dimensional simulations are compared to quantify the influence that surrounding geologic media has on the CO2 leakage rate.
Ethanol and ethanol/gasoline blends are being widely considered as alternative fuels for light-duty automotive applications. At the same time, HCCI combustion has the potential to provide high efficiency and ultra-low exhaust emissions. However, the application of HCCI is typically limited to low and moderate loads because of unacceptably high heat-release rates (HRR) at higher fueling rates. This work investigates the potential of lowering the HCCI HRR at high loads by using partial fuel stratification to increase the in-cylinder thermal stratification. This strategy is based on ethanol's high heat of vaporization combined with its true single-stage ignition characteristics. Using partial fuel stratification, the strong fuel-vaporization cooling produces thermal stratification due to variations in the amount of fuel vaporization in different parts of the combustion chamber. The low sensitivity of the autoignition reactions to variations of the local fuel concentration allows the temperature variations to govern the combustion event. This results in a sequential autoignition event from leaner and hotter zones to richer and colder zones, lowering the overall combustion rate compared to operation with a uniform fuel/air mixture. The amount of partial fuel stratification was varied by adjusting the fraction of fuel injected late to produce stratification, and also by changing the timing of the late injection. The experiments show that a combination of 60-70% premixed charge and injection of 30-40 % of the fuel at 80{sup o}CA before TDC is effective for smoothing the HRR. With CA50 held fixed, this increases the burn duration by 55% and reduces the maximum pressure-rise rate by 40%. Combustion stability remains high but engine-out NO{sub x} has to be monitored carefully. For operation with strong reduction of the peak HRR, ISNO{sub x} rises to around 0.20 g/kWh for an IMEP{sub g} of 440 kPa. The single-cylinder HCCI research engine was operated naturally aspirated without EGR at 1200 rpm, and had low residual level using a CR = 14 piston.
A silicon photonics based integrated optical phase locked loop is utilized to synchronize a 10.2 GHz voltage controlled oscillator with a 509 MHz mode locked laser, achieving 32 fs integrated jitter over 300 kHz bandwidth.
Mudstone mechanical testing is often limited by poor core recovery and sample size, preservation and preparation issues, which can lead to sampling bias, damage, and time-dependent effects. A micropillar compression technique, originally developed by Uchic et al. 2004, here is applied to elasto-plastic deformation of small volumes of mudstone, in the range of cubic microns. This study examines behavior of the Gothic shale, the basal unit of the Ismay zone of the Pennsylvanian Paradox Formation and potential shale gas play in southeastern Utah, USA. Precision manufacture of micropillars 5 microns in diameter and 10 microns in length are prepared using an ion-milling method. Characterization of samples is carried out using: dual focused ion - scanning electron beam imaging of nano-scaled pores and distribution of matrix clay and quartz, as well as pore-filling organics; laser scanning confocal (LSCM) 3D imaging of natural fractures; and gas permeability, among other techniques. Compression testing of micropillars under load control is performed using two different nanoindenter techniques. Deformation of 0.5 cm in diameter by 1 cm in length cores is carried out and visualized by a microscope loading stage and laser scanning confocal microscopy. Axisymmetric multistage compression testing and multi-stress path testing is carried out using 2.54 cm plugs. Discussion of results addresses size of representative elementary volumes applicable to continuum-scale mudstone deformation, anisotropy, and size-scale plasticity effects. Other issues include fabrication-induced damage, alignment, and influence of substrate.
Arctic sea ice is an important component of the global climate system and, due to feedback effects, the Arctic ice cover is changing rapidly. Predictive mathematical models are of paramount importance for accurate estimates of the future ice trajectory. However, the sea ice components of Global Climate Models (GCMs) vary significantly in their prediction of the future state of Arctic sea ice and have generally underestimated the rate of decline in minimum sea ice extent seen over the past thirty years. One of the contributing factors to this variability is the sensitivity of the sea ice state to internal model parameters. A new sea ice model that holds some promise for improving sea ice predictions incorporates an anisotropic elastic-decohesive rheology and dynamics solved using the material-point method (MPM), which combines Lagrangian particles for advection with a background grid for gradient computations. We evaluate the variability of this MPM sea ice code and compare it with the Los Alamos National Laboratory CICE code for a single year simulation of the Arctic basin using consistent ocean and atmospheric forcing. Sensitivities of ice volume, ice area, ice extent, root mean square (RMS) ice speed, central Arctic ice thickness,and central Arctic ice speed with respect to ten different dynamic and thermodynamic parameters are evaluated both individually and in combination using the Design Analysis Kit for Optimization and Terascale Applications (DAKOTA). We find similar responses for the two codes and some interesting seasonal variability in the strength of the parameters on the solution.
Measurements of the electrical and thermal transport properties of one-dimensional nanostructures (e.g., nanotubes and nanowires) typically are obtained without detailed knowledge of the specimen's atomic-scale structure or defects. To address this deficiency we have developed a microfabricated, chip-based characterization platform that enables both transmission electron microscopy (TEM) of atomic structure and defects as well as measurement of the thermal transport properties of individual nanostructures. The platform features a suspended heater line that contacts the center of a suspended nanostructure/nanowire that was placed using in-situ scanning electron microscope nanomanipulators. One key advantage of this platform is that it is possible to measure the thermal conductivity of both halves of the nanostructure (on each side of the central heater), and this feature permits identification of possible changes in thermal conductance along the wire and measurement of the thermal contact resistance. Suspension of the nanostructure across a through-hole enables TEM characterization of the atomic and defect structure (dislocations, stacking faults, etc.) of the test sample. As a model study, we report the use of this platform to measure the thermal conductivity and defect structure of GaN nanowires. The utilization of this platform for the measurements of other nanostructures will also be discussed.
Sandia National Laboratories is collaborating with the National Research Council (NRC) Canada and the National Renewable Energy Laboratory (NREL) to develop a decision-support model that will evaluate the tradeoffs associated with high-latitude algae biofuel production co-located with wastewater, CO2, and waste heat. This project helps Canada meet its goal of diversifying fuel sources with algae-based biofuels. The biofuel production will provide a wide range of benefits including wastewater treatment, CO2 reuse and reduction of demand for fossil-based fuels. The higher energy density in algae-based fuels gives them an advantage over crop-based biofuels as the 'production' footprint required is much less, resulting in less water consumed and little, if any conversion of agricultural land from food to fuel production. Besides being a potential source for liquid fuel, algae have the potential to be used to generate electricity through the burning of dried biomass, or anaerobically digested to generate methane for electricity production. Co-locating algae production with waste streams may be crucial for making algae an economically valuable fuel source, and will certainly improve its overall ecological sustainability. The modeling process will address these questions, and others that are important to the use of water for energy production: What are the locations where all resources are co-located, and what volumes of algal biomass and oil can be produced there? In locations where co-location does not occur, what resources should be transported, and how far, while maintaining economic viability? This work is being funded through the U.S. Department of Energy (DOE) Biomass Program Office of Energy Efficiency and Renewable Energy, and is part of a larger collaborative effort that includes sampling, strain isolation, strain characterization and cultivation being performed by the NREL and Canada's NRC. Results from the NREL / NRC collaboration including specific productivities of selected algal strains will eventually be incorporated into this model.
Concerns over rising concentrations of greenhouse gases in the atmosphere have resulted in serious consideration of policies aimed at reduction of anthropogenic carbon dioxide (CO2) emissions. If large scale abatement efforts are undertaken, one critical tool will be geologic sequestration of CO2 captured from large point sources, specifically coal and natural gas fired power plants. Current CO2 capture technologies exact a substantial energy penalty on the source power plant, which must be offset with make-up power. Water demands increase at the source plant due to added cooling loads. In addition, new water demand is created by water requirements associated with generation of the make-up power. At the sequestration site however, saline water may be extracted to manage CO2 plum migration and pressure build up in the geologic formation. Thus, while CO2 capture creates new water demands, CO2 sequestration has the potential to create new supplies. Some or all of the added demand may be offset by treatment and use of the saline waters extracted from geologic formations during CO2 sequestration. Sandia National Laboratories, with guidance and support from the National Energy Technology Laboratory, is creating a model to evaluate the potential for a combined approach to saline formations, as a sink for CO2 and a source for saline waters that can be treated and beneficially reused to serve power plant water demands. This presentation will focus on the magnitude of added U.S. power plant water demand under different CO2 emissions reduction scenarios, and the portion of added demand that might be offset by saline waters extracted during the CO2 sequestration process.
Because the potential effects of climate change are more severe than had previously been thought, increasing focus on uncertainty quantification is required for risk assessment needed by policy makers. Current scientific efforts focus almost exclusively on establishing best estimates of future climate change. However, the greatest consequences occur in the extreme tail of the probability density functions for climate sensitivity (the 'high-sensitivity tail'). To this end, we are exploring the impacts of newly postulated, highly uncertain, but high-consequence physical mechanisms to better establish the climate change risk. We define consequence in terms of dramatic change in physical conditions and in the resulting socioeconomic impact (hence, risk) on populations. Although we are developing generally applicable risk assessment methods, we have focused our initial efforts on uncertainty and risk analyses for the Arctic region. Instead of focusing on best estimates, requiring many years of model parameterization development and evaluation, we are focusing on robust emergent phenomena (those that are not necessarily intuitive and are insensitive to assumptions, subgrid-parameterizations, and tunings). For many physical systems, under-resolved models fail to generate such phenomena, which only develop when model resolution is sufficiently high. Our ultimate goal is to discover the patterns of emergent climate precursors (those that cannot be predicted with lower-resolution models) that can be used as a 'sensitivity fingerprint' and make recommendations for a climate early warning system that would use satellites and sensor arrays to look for the various predicted high-sensitivity signatures. Our initial simulations are focused on the Arctic region, where underpredicted phenomena such as rapid loss of sea ice are already emerging, and because of major geopolitical implications associated with increasing Arctic accessibility to natural resources, shipping routes, and strategic locations. We anticipate that regional climate will be strongly influenced by feedbacks associated with a seasonally ice-free Arctic, but with unknown emergent phenomena.
Uncertainty quantification in complex climate models is challenged by the sparsity of available climate model predictions due to the high computational cost of model runs. Another feature that prevents classical uncertainty analysis from being readily applicable is bifurcative behavior in climate model response with respect to certain input parameters. A typical example is the Atlantic Meridional Overturning Circulation. The predicted maximum overturning stream function exhibits discontinuity across a curve in the space of two uncertain parameters, namely climate sensitivity and CO2 forcing. We outline a methodology for uncertainty quantification given discontinuous model response and a limited number of model runs. Our approach is two-fold. First we detect the discontinuity with Bayesian inference, thus obtaining a probabilistic representation of the discontinuity curve shape and location for arbitrarily distributed input parameter values. Then, we construct spectral representations of uncertainty, using Polynomial Chaos (PC) expansions on either side of the discontinuity curve, leading to an averaged-PC representation of the forward model that allows efficient uncertainty quantification. The approach is enabled by a Rosenblatt transformation that maps each side of the discontinuity to regular domains where desirable orthogonality properties for the spectral bases hold. We obtain PC modes by either orthogonal projection or Bayesian inference, and argue for a hybrid approach that targets a balance between the accuracy provided by the orthogonal projection and the flexibility provided by the Bayesian inference - where the latter allows obtaining reasonable expansions without extra forward model runs. The model output, and its associated uncertainty at specific design points, are then computed by taking an ensemble average over PC expansions corresponding to possible realizations of the discontinuity curve. The methodology is tested on synthetic examples of discontinuous model data with adjustable sharpness and structure.
CACTUS (Code for Axial and Cross-flow TUrbine Simulation) is a turbine performance simulation code, based on a free wake vortex method, under development at Sandia National Laboratories (SNL) as part of a Department of Energy program to study marine hydrokinetic (MHK) devices. The current effort builds upon work previously done at SNL in the area of vertical axis wind turbine simulation, and aims to add models to handle generic device geometry and physical models specific to the marine environment. An overview of the current state of the project and validation effort is provided.
We present results from a recently developed multiscale inversion technique for binary media, with emphasis on the effect of subgrid model errors on the inversion. Binary media are a useful fine-scale representation of heterogeneous porous media. Averaged properties of the binary field representations can be used to characterize flow through the porous medium at the macroscale. Both direct measurements of the averaged properties and upscaling are complicated and may not provide accurate results. However, it may be possible to infer upscaled properties of the binary medium from indirect measurements at the coarse scale. Multiscale inversion, performed with a subgrid model to connect disparate scales together, can also yield information on the fine-scale properties. We model the binary medium using truncated Gaussian fields, and develop a subgrid model for the upscaled permeability based on excursion sets of those fields. The subgrid model requires an estimate of the proportion of inclusions at the block scale as well as some geometrical parameters of the inclusions as inputs, and predicts the effective permeability. The inclusion proportion is assumed to be spatially varying, modeled using Gaussian processes and represented using a truncated Karhunen-Louve (KL) expansion. This expansion is used, along with the subgrid model, to pose as a Bayesian inverse problem for the KL weights and the geometrical parameters of the inclusions. The model error is represented in two different ways: (1) as a homoscedastic error and (2) as a heteroscedastic error, dependent on inclusion proportionality and geometry. The error models impact the form of the likelihood function in the expression for the posterior density of the objects of inference. The problem is solved using an adaptive Markov Chain Monte Carlo method, and joint posterior distributions are developed for the KL weights and inclusion geometry. Effective permeabilities and tracer breakthrough times at a few 'sensor' locations (obtained by simulating a pump test) form the observables used in the inversion. The inferred quantities can be used to generate an ensemble of permeability fields, both upscaled and fine-scale, which are consistent with the observations. We compare the inferences developed using the two error models, in terms of the KL weights and fine-scale realizations that could be supported by the coarse-scale inferences. Permeability differences are observed mainly in regions where the inclusions proportion is near the percolation threshold, and the subgrid model incurs its largest approximation. These differences also reflected in the tracer breakthrough times and the geometry of flow streamlines, as obtained from a permeameter simulation. The uncertainty due to subgrid model error is also compared to the uncertainty in the inversion due to incomplete data.
Subsurface containment of CO2 is predicated on effective caprock sealing. Many previous studies have relied on macroscopic measurements of capillary breakthrough pressure and other petrophysical properties without direct examination of solid phases that line pore networks and directly contact fluids. However, pore-lining phases strongly contribute to sealing behavior through interfacial interactions among CO2, brine, and the mineral or non-mineral phases. Our high resolution (i.e., sub-micron) examination of the composition of pore-lining phases of several continental and marine mudstones indicates that sealing efficiency (i.e., breakthrough pressure) is governed by pore shapes and pore-lining phases that are not identifiable except through direct characterization of pores. Bulk X-ray diffraction data does not indicate which phases line the pores and may be especially lacking for mudstones with organic material. Organics can line pores and may represent once-mobile phases that modify the wettability of an originally clay-lined pore network. For shallow formations (i.e., < {approx}800 m depth), interfacial tension and contact angles result in breakthrough pressures that may be as high as those needed to fracture the rock - thus, in the absence of fractures, capillary sealing efficiency is indicated. Deeper seals have poorer capillary sealing if mica-like wetting dominates the wettability. We thank the U.S. Department of Energy's National Energy Technology Laboratory and the Office of Basic Energy Sciences, and the Southeast and Southwest Carbon Sequestration Partnerships for supporting this work.
Injection of CO2 into underground rock formations can reduce atmospheric CO2 emissions. Caprocks present above potential storage formations are the main structural trap inhibiting CO2 from leaking into overlying aquifers or back to the Earth's surface. Dissolution and precipitation of caprock minerals resulting from reaction with CO2 may alter the pore network where many pores are of the micrometer to nanometer scale, thus altering the structural trapping potential of the caprock. However, the distribution, geometry and volume of pores at these scales are poorly characterized. In order to evaluate the overall risk of leakage of CO2 from storage formations, a first critical step is understanding the distribution and shape of pores in a variety of different caprocks. As the caprock is often comprised of mudstones, we analyzed samples from several mudstone formations with small angle neutron scattering (SANS) and high-resolution transmission electron microscopy (TEM) imaging to compare the pore networks. Mudstones were chosen from current or potential sites for carbon sequestration projects including the Marine Tuscaloosa Group, the Lower Tuscaloosa Group, the upper and lower shale members of the Kirtland Formation, and the Pennsylvanian Gothic shale. Expandable clay contents ranged from 10% to approximately 40% in the Gothic shale and Kirtland Formation, respectively. During SANS, neutrons effectively scatter from interfaces between materials with differing scattering length density (i.e., minerals and pores). The intensity of scattered neutrons, I(Q), where Q is the scattering vector, gives information about the volume and arrangement of pores in the sample. The slope of the scattering data when plotted as log I(Q) vs. log Q provides information about the fractality or geometry of the pore network. On such plots slopes from -2 to -3 represent mass fractals while slopes from -3 to -4 represent surface fractals. Scattering data showed surface fractal dimensions for the Kirtland formation and one sample from the Tuscaloosa formation close to 3, indicating very rough surfaces. In contrast, scattering data for the Gothic shale formation exhibited mass fractal behavior. In one sample of the Tuscaloosa formation the data are described by a surface fractal at low Q (larger pores) and a mass fractal at high Q (smaller pores), indicating two pore populations contributing to the scattering behavior. These small angle neutron scattering results, combined with high-resolution TEM imaging, provided a means for both qualitative and quantitative analysis of the differences in pore networks between these various mudstones.
We consider the problem of placing a limited number of sensors in a municipal water distribution network to minimize the impact over a given suite of contamination incidents. In its simplest form, the sensor placement problem is a p-median problem that has structure extremely amenable to exact and heuristic solution methods. We describe the solution of real-world instances using integer programming or local search or a Lagrangian method. The Lagrangian method is necessary for solution of large problems on small PCs. We summarize a number of other heuristic methods for effectively addressing issues such as sensor failures, tuning sensors based on local water quality variability, and problem size/approximation quality tradeoffs. These algorithms are incorporated into the TEVA-SPOT toolkit, a software suite that the US Environmental Protection Agency has used and is using to design contamination warning systems for US municipal water systems.
The acoustic field generated during a Direct Field Acoustic Test (DFAT) has been analytically modeled in two space dimensions using a properly phased distribution of propagating plane waves. Both the pure-tone and broadband acoustic field were qualitatively and quantitatively compared to a diffuse acoustic field. The modeling indicates significant non-uniformity of sound pressure level for an empty (no test article) DFAT, specifically a center peak and concentric maxima/minima rings. This spatial variation is due to the equivalent phase among all propagating plane waves at each frequency. The excitation of a simply supported slender beam immersed within the acoustic fields was also analytically modeled. Results indicate that mid-span response is dependent upon location and orientation of the beam relative to the center of the DFAT acoustic field. For a diffuse acoustic field, due to its spatial uniformity, mid-span response sensitivity to location and orientation is nonexistent.
Mercury intrusion porosimetry (MIP) is an often-applied technique for determining pore throat distributions and seal analysis of fine-grained rocks. Due to closure effects, potential pore collapse, and complex pore network topologies, MIP data interpretation can be ambiguous, and often biased toward smaller pores in the distribution. We apply 3D imaging techniques and lattice-Boltzmann modeling in interpreting MIP data for samples of the Cretaceous Selma Group Chalk. In the Mississippi Interior Salt Basin, the Selma Chalk is the apparent seal for oil and gas fields in the underlying Eutaw Fm., and, where unfractured, the Selma Chalk is one of the regional-scale seals identified by the Southeast Regional Carbon Sequestration Partnership for CO2 injection sites. Dual focused ion - scanning electron beam and laser scanning confocal microscopy methods are used for 3D imaging of nanometer-to-micron scale microcrack and pore distributions in the Selma Chalk. A combination of image analysis software is used to obtain geometric pore body and throat distributions and other topological properties, which are compared to MIP results. 3D data sets of pore-microfracture networks are used in Lattice Boltzmann simulations of drainage (wetting fluid displaced by non-wetting fluid via the Shan-Chen algorithm), which in turn are used to model MIP procedures. Results are used in interpreting MIP results, understanding microfracture-matrix interaction during multiphase flow, and seal analysis for underground CO2 storage.
Photovoltaic systems are often priced in $/W{sub p}, where Wp refers to the DC power rating of the modules at Standard Test Conditions (1000 W/m{sup 2}, 25 C cell temperature) and $ refers to the installed cost of the system. However, the true value of the system is in the energy it will produce in kWhs, not the power rating. System energy production is a function of the system design and location, the mounting configuration, the power conversion system, and the module technology, as well as the solar resource. Even if all other variables are held constant, the annual energy yield (kWh/kW{sup p}) will vary among module technologies because of differences in response to low-light levels and temperature. Understanding energy yield is a key part of understanding system value. System performance models are used during project development to estimate the expected output of PV systems for a given design and location. Performance modeling is normally done by the system designer/system integrator. Often, an independent engineer will also model system output during a due diligence review of a project. A variety of system performance models are available. The most commonly used modeling tool for project development and due diligence in the United States is probably PVsyst, while those seeking a quick answer to expected energy production may use PVWatts. In this paper, we examine the variation in predicted energy output among modeling tools and users and compare that to measured output.
To-date, the application of high-performance computing resources to Semantic Web data has largely focused on commodity hardware and distributed memory platforms. In this paper we make the case that more specialized hardware can offer superior scaling and close to an order of magnitude improvement in performance. In particular we examine the Cray XMT. Its key characteristics, a large, global shared-memory, and processors with a memory-latency tolerant design, offer an environment conducive to programming for the Semantic Web and have engendered results that far surpass current state of the art. We examine three fundamental pieces requisite for a fully functioning semantic database: dictionary encoding, RDFS inference, and query processing. We show scaling up to 512 processors (the largest configuration we had available), and the ability to process 20 billion triples completely in-memory.
Dissolved CO2 from geological CO2 sequestration may react with dissolved minerals in fractured rocks or confined aquifers and cause mineral precipitation. The overall rate of reaction can be limited by diffusive or dispersive mixing, and mineral precipitation can block pores and further hinder these processes. Mixing-induced calcite precipitation experiments were performed by injecting solutions containing CaCl2 and Na2CO3 through two separate inlets of a micromodel (1-cm x 2-cm x 40-microns); transverse dispersion caused the two solutions to mix along the center of the micromodel, resulting in calcite precipitation. The amount of calcite precipitation initially increased to a maximum and then decreased to a steady state value. Fluorescent microscopy and imaging techniques were used to visualize calcite precipitation, and the corresponding effects on the flow field. Experimental micromodel results were evaluated with pore-scale simulations using a 2-D Lattice-Boltzmann code for water flow and a finite volume code for reactive transport. The reactive transport model included the impact of pH upon carbonate speciation and calcite dissolution. We found that proper estimation of the effective diffusion coefficient and the reaction surface area is necessary to adequately simulate precipitation and dissolution rates. The effective diffusion coefficient was decreased in grid cells where calcite precipitated, and keeping track of reactive surface over time played a significant role in predicting reaction patterns. Our results may improve understanding of the fundamental physicochemical processes during CO2 sequestration in geologic formations.
Multi-scale binary permeability field estimation from static and dynamic data is completed using Markov Chain Monte Carlo (MCMC) sampling. The binary permeability field is defined as high permeability inclusions within a lower permeability matrix. Static data are obtained as measurements of permeability with support consistent to the coarse scale discretization. Dynamic data are advective travel times along streamlines calculated through a fine-scale field and averaged for each observation point at the coarse scale. Parameters estimated at the coarse scale (30 x 20 grid) are the spatially varying proportion of the high permeability phase and the inclusion length and aspect ratio of the high permeability inclusions. From the non-parametric, posterior distributions estimated for these parameters, a recently developed sub-grid algorithm is employed to create an ensemble of realizations representing the fine-scale (3000 x 2000), binary permeability field. Each fine-scale ensemble member is instantiated by convolution of an uncorrelated multiGaussian random field with a Gaussian kernel defined by the estimated inclusion length and aspect ratio. Since the multiGaussian random field is itself a realization of a stochastic process, the procedure for generating fine-scale binary permeability field realizations is also stochastic. Two different methods are hypothesized to perform posterior predictive tests. Different mechanisms for combining multi Gaussian random fields with kernels defined from the MCMC sampling are examined. Posterior predictive accuracy of the estimated parameters is assessed against a simulated ground truth for predictions at both the coarse scale (effective permeabilities) and at the fine scale (advective travel time distributions). The two techniques for conducting posterior predictive tests are compared by their ability to recover the static and dynamic data. The skill of the inference and the method for generating fine-scale binary permeability fields are evaluated through flow calculations on the resulting fields using fine-scale realizations and comparing them against results obtained with the ground truth fine-scale and coarse-scale permeability fields.
We present advances with a 32 element scalable, segmented dual mode imager. Scaling up the number of cells results in a 1.4 increase in efficiency over a system we deployed last year. Variable plane separation has been incorporated which further improves the efficiency of the detector. By using 20 cm diameter cells we demonstrate that we could increase sensitivity by a factor of 6. We further demonstrate gamma ray imaging in from Compton scattering. This feature allows for powerful dual mode imaging. Selected results are presented that demonstrate these new capabilities.
The neutron scatter camera was originally developed for a range of SNM detection applications. We are now exploring the feasibility of applications in treaty verification and warhead monitoring using experimentation, maximum likelihood estimation method (MLEM), detector optimization, and MCNP-PoliMi simulations.
A rational approach was used to design polymeric materials for thin-film electronics applications, whereby theoretical modeling was used to determine synthetic targets. Time-dependent density functional theory calculations were used as a tool to predict the electrical properties of conjugated polymer systems. From these results, polymers with desirable energy levels and band-gaps were designed and synthesized. Measurements of optoelectronic properties were performed on the synthesized polymers and the results were compared to those of the theoretical model. From this work, the efficacy of the model was evaluated and new target polymers were identified.
Decisions made to address climate change must start with an understanding of the risk of an uncertain future to human systems, which in turn means understanding both the consequence as well as the probability of a climate induced impact occurring. In other words, addressing climate change is an exercise in risk-informed policy making, which implies that there is no single correct answer or even a way to be certain about a single answer; the uncertainty in future climate conditions will always be present and must be taken as a working-condition for decision making. In order to better understand the implications of uncertainty on risk and to provide a near-term rationale for policy interventions, this study estimates the impacts from responses to climate change on U.S. state- and national-level economic activity by employing a risk-assessment methodology for evaluating uncertain future climatic conditions. Using the results from the Intergovernmental Panel on Climate Change's (IPCC) Fourth Assessment Report (AR4) as a proxy for climate uncertainty, changes in hydrology over the next 40 years were mapped and then modeled to determine the physical consequences on economic activity and to perform a detailed 70-industry analysis of the economic impacts among the interacting lower-48 states. The analysis determines industry-level effects, employment impacts at the state level, interstate population migration, consequences to personal income, and ramifications for the U.S. trade balance. The conclusions show that the average risk of damage to the U.S. economy from climate change is on the order of $1 trillion over the next 40 years, with losses in employment equivalent to nearly 7 million full-time jobs. Further analysis shows that an increase in uncertainty raises this risk. This paper will present the methodology behind the approach, a summary of the underlying models, as well as the path forward for improving the approach.
Motivated by social network data mining problems such as link prediction and collaborative filtering, significant research effort has been devoted to computing topological measures including the Katz score and the commute time. Existing approaches typically approximate all pairwise relationships simultaneously. In this paper, we are interested in computing: the score for a single pair of nodes, and the top-k nodes with the best scores from a given source node. For the pairwise problem, we apply an iterative algorithm that computes upper and lower bounds for the measures we seek. This algorithm exploits a relationship between the Lanczos process and a quadrature rule. For the top-k problem, we propose an algorithm that only accesses a small portion of the graph and is related to techniques used in personalized PageRank computing. To test the scalability and accuracy of our algorithms we experiment with three real-world networks and find that these algorithms run in milliseconds to seconds without any preprocessing.
We describe an integrated information management system for an independent spent fuel dry-storage installation (ISFSI) that can provide for (1) secure and authenticated data collection, (2) data analysis, (3) dissemination of information to appropriate stakeholders via a secure network, and (4) increased public confidence and support of the facility licensing and operation through increased transparency. This information management system is part of a collaborative project between Sandia National Laboratories, Taiwan Power Co., and the Fuel Cycle Materials Administration of Taiwan's Atomic Energy Council, which is investigating how to implement this concept.
A novel multiphase shock tube has been constructed to test the interaction of a planar shock wave with a dense gas-solid field of particles. The particle field is generated by a gravity-fed method that results in a spanwise curtain of 100-micron particles producing a volume fraction of about 15%. Interactions with incident shock Mach numbers of 1.67 and 1.95 are reported. High-speed schlieren imaging is used to reveal the complex wave structure associated with the interaction. After the impingement of the incident shock, transmitted and reflected shocks are observed, which lead to differences in flow properties across the streamwise dimension of the curtain. Tens of microseconds after the onset of the interaction, the particle field begins to propagate downstream, and disperse. The spread of the particle field, as a function of its position, is seen to be nearly identical for both Mach numbers. Immediately downstream of the curtain, the peak pressures associated with the Mach 1.67 and 1.95 interactions are about 35% and 45% greater than tests without particles, respectively. For both Mach numbers tested, the energy and momentum fluxes in the induced flow far downstream are reduced by about 30-40% by the presence of the particle field.
The very high temperature reactor (VHTR) concept is being developed by the US Department of Energy (DOE) and other groups around the world for the future generation of electricity at high thermal efficiency (> 48%) and co-generation of hydrogen and process heat. This Generation-IV reactor would operate at elevated exit temperatures of 1,000-1,273 K, and the fueled core would be cooled by forced convection helium gas. For the prismatic-core VHTR, which is the focus of this analysis, the velocity of the hot helium flow exiting the core into the lower plenum (LP) could be 35-70 m/s. The impingement of the resulting gas jets onto the adiabatic plate at the bottom of the LP could develop hot spots and thermal stratification and inadequate mixing of the gas exiting the vessel to the turbo-machinery for energy conversion. The complex flow field in the LP is further complicated by the presence of large cylindrical graphite posts that support the massive core and inner and outer graphite reflectors. Because there are approximately 276 channels in the VHTR core from which helium exits into the LP and a total of 155 support posts, the flow field in the LP includes cross flow, multiple jet flow interaction, flow stagnation zones, vortex interaction, vortex shedding, entrainment, large variation in Reynolds number (Re), recirculation, and mixing enhancement and suppression regions. For such a complex flow field, experimental results at operating conditions are not currently available. Instead, the objective of this paper is to numerically simulate the flow field in the LP of a prismatic core VHTR using the Sandia National Laboratories Fuego, which is a 3D, massively parallel generalized computational fluid dynamics (CFD) code with numerous turbulence and buoyancy models and simulation capabilities for complex gas flow fields, with and without thermal effects. The code predictions for simpler flow fields of single and swirling gas jets, with and without a cross flow, are validated using reported experimental data and theory. The key processes in the LP are identified using phenomena identification and ranking table (PIRT). It may be argued that a CFD code that accurately simulates simplified, single-effect flow fields with increasing complexity is likely to adequately model the complex flow field in the VHTR LP, subject to a future experimental validation. The PIRT process and spatial and temporal discretizations implemented in the present analysis using Fuego established confidence in the validation and verification (V and V) calculations and in the conclusions reached based on the simulation results. The performed calculations included the helicoid vortex swirl model, the dynamic Smagorinsky large eddy simulation (LES) turbulence model, participating media radiation (PMR), and 1D conjugate heat transfer (CHT). The full-scale, half-symmetry LP mesh used in the LP simulation included unstructured hexahedral elements and accounted for the graphite posts, the helium jets, the exterior walls, and the bottom plate with an adiabatic outer surface. Results indicated significant enhancements in heat transfer, flow mixing, and entrainment in the VHTR LP when using swirling inserts at the exit of the helium flow channels into the LP. The impact of using various swirl angles on the flow mixing and heat transfer in the LP is qualified, including the formation of the central recirculation zone (CRZ), and the effect of LP height. Results also showed that in addition to the enhanced mixing, the swirling inserts result in negligible additional pressure losses and are likely to eliminate the formation of hot spots.
The development of turbulent spots in a hypersonic boundary layer was studied on the nozzle wall of the Boeing/AFOSR Mach-6 Quiet Tunnel. Under quiet flow conditions, the nozzle wall boundary layer remains laminar and grows very thick over the long nozzle length. This allows the development of large turbulent spots that can be readily measured with pressure transducers. Measurements of naturally occurring wave packets and developing turbulent spots were made. The peak frequencies of these natural wave packets were in agreement with second-mode computations. For a controlled study, the breakdown of disturbances created by spark and glow perturbations were studied at similar freestream conditions. The spark perturbations were the most effective at creating large wave packets that broke down into turbulent spots. The flow disturbances created by the controlled perturbations were analyzed to obtain amplitude criteria for nonlinearity and breakdown as well as the convection velocities of the turbulent spots. Disturbances first grew into linear instability waves and then quickly became nonlinear. Throughout the nonlinear growth of the wave packets, large harmonics are visible in the power spectra. As breakdown begins, the peak amplitudes of the instability waves and harmonics decrease into the rising broad-band frequencies. Instability waves are still visible on either side of the growing turbulent spots during this breakdown process.
In recent years there has been an unstable supply of the critical diagnostic medical isotope 99Tc. Several concepts and designs have been proposed to produce 99Mo the parent nuclide of 99Tc, at a commercial scale sufficient to stabilize the world supply. This work lays out a testing and experiment plan for a proposed 2 MW open pool reactor fueled by Low Enriched Uranium (LEU) 99Mo targets. The experiments and tests necessary to support licensing of the reactor design are described and how these experiments and tests will help establish the safe operating envelop for a medical isotope production reactor is discussed. The experiments and tests will facilitate a focused and efficient licensing process in order to bring on line a needed production reactor dedicated to supplying medical isotopes. The Target Fuel Isotope Reactor (TFIR) design calls for an active core region that is approximately 40 cm in diameter and 40 cm in fuel height. It contains up to 150 cylindrical, 1-cm diameter, LEU oxide fuel pins clad with Zircaloy (zirconium alloy), in an annular hexagonal array on a {approx}2.0 cm pitch surrounded, radially, by a graphite or a Be reflector. The reactor is similar to U.S. university reactors in power, hardware, and safety/control systems. Fuel/target pin fabrication is based on existing light water reactor fuel fabrication processes. However, as part of licensing process, experiments must be conducted to confirm analytical predictions of steady-state power and accident conditions. The experiment and test plan will be conducted in phases and will utilize existing facilities at the U.S. Department of Energy's Sandia National Laboratories. The first phase is to validate the predicted reactor core neutronics at delayed critical, zero power and very low power. This will be accomplished by using the Sandia Critical Experiment (CX) platform. A full scale TFIR core will be built in the CX and delayed critical measurements will be taken. For low power experiments, fuel pins can be removed after the experiment and using Sandia's metrology lab, relative power profiles (radially and axially) can be determined. In addition to validating neutronic analyses, confirming heat transfer properties of the target/fuel pins and core will be conducted. Fuel/target pin power limits can be verified with out-of-pile (electrical heating) thermal-hydraulic experiments. This will yield data on the heat flux across the Zircaloy clad and establish safety margin and operating limits. Using Sandia's Annular Core Research Reactor (ACRR) a 4 MW TRIGA type research reactor, target/fuel pins can be driven to desired fission power levels for long durations. Post experiment inspection of the pins can be conducted in the Auxiliary Hot Cell Facility to observe changes in the mechanical properties of the LEU matrix and burn-up effects. Transient tests can also be conducted at the ACRR to observe target/fuel pin performance during accident conditions. Target/fuel pins will be placed in double experiment containment and driven by pulsing the ACRR until target/fuel failure is observed. This will allow for extrapolation of analytical work to confirm safety margins.
Although planar heterostructures dominate current solid-state lighting architectures (SSL), 1D nanowires have distinct and advantageous properties that may eventually enable higher efficiency, longer wavelength, and cheaper devices. However, in order to fully realize the potential of nanowire-based SSL, several challenges exist in the areas of controlled nanowire synthesis, nanowire device integration, and understanding and controlling the nanowire electrical, optical, and thermal properties. Here recent results are reported regarding the aligned growth of GaN and III-nitride core-shell nanowires, along with extensive results providing insights into the nanowire properties obtained using cutting-edge structural, electrical, thermal, and optical nanocharacterization techniques. A new top-down fabrication method for fabricating periodic arrays of GaN nanorods and subsequent nanorod LED fabrication is also presented.
Cielo, a Cray XE6, is the Department of Energy NNSA Advanced Simulation and Computing (ASC) campaign's newest capability machine. Rated at 1.37 PFLOPS, it consists of 8,944 dual-socket oct-core AMD Magny-Cours compute nodes, linked using Cray's Gemini interconnect. Its primary mission objective is to enable a suite of the ASC applications implemented using MPI to scale to tens of thousands of cores. Cielo is an evolutionary improvement to a successful architecture previously available to many of our codes, thus enabling a basis for understanding the capabilities of this new architecture. Using three codes strategically important to the ASC campaign, and supplemented with some micro-benchmarks that expose the fundamental capabilities of the XE6, we report on the performance characteristics and capabilities of Cielo.
A conceptual structure for performance assessments (PAs) for radioactive waste disposal facilities and other complex engineered facilities based on the following three basic conceptual entities is described: EN1, a probability space that characterizes aleatory uncertainty; EN2, a function that predicts consequences for individual elements of the sample space for aleatory uncertainty; and EN3, a probability space that characterizes epistemic uncertainty. The implementation of this structure is illustrated with results from PAs for the Waste Isolation Pilot Plant and the proposed Yucca Mountain repository for high-level radioactive waste.
In this work, we describe a novel design for a H2SO 4decomposer. The decomposition of H2SO4 to produce SO2is a common processing operation in the sulfur-based thermochemical cycles for hydrogen production where acid decomposition takes place at 850°C in the presence of a catalyst. The combination of high temperature and sulfuric acid creates a very corrosive environment that presents significant design challenges. The new decomposer design is based on a bayonet-type heat exchanger tube with the annular space packed with a catalyst. The unit is constructed of silicon carbide and other highly corrosion resistant materials. The new design integrates acid boiling, superheating, decomposition and heat recuperation into a single process and eliminates problems of corrosion and failure of high temperature seals encountered in previous testing using metallic construction materials. The unit was tested by varying the acid feed rate and decomposition temperature and pressure.
This report presents the results of an aging experiment that was established in FY09 and completed in FY10 for the Sandia MEMS Passive Shock Sensor. A total of 37 packages were aged at different temperatures and times, and were then tested after aging to determine functionality. Aging temperatures were selected at 100 C and 150 C, with times ranging from as short as 100 hours to as long as 1 year to simulate a predicted aging of up to 20 years. In all of the tests and controls, 100% of the devices continued to function normally.
Mandell, John F.; Samborsky, Daniel D.; Agastra, Pancasatya; Sears, Aaron T.; Wilson, Timothy J.; Ashwill, Thomas D.; Laird, Daniel L.
This report presents an analysis of trends in fatigue results from the Montana State University program on the fatigue of composite materials for wind turbine blades for the period 2005-2009. Test data can be found in the SNL/MSU/DOE Fatigue of Composite Materials Database which is updated annually. This is the fifth report in this series, which summarizes progress of the overall program since its inception in 1989. The primary thrust of this program has been research and testing of a broad range of structural laminate materials of interest to blade structures. The report is focused on current types of infused and prepreg blade materials, either processed in-house or by industry partners. Trends in static and fatigue performance are analyzed for a range of materials, geometries and loading conditions. Materials include: sixteen resins of three general types, five epoxy based paste adhesives, fifteen reinforcing fabrics including three fiber types, three prepregs, many laminate lay-ups and process variations. Significant differences in static and fatigue performance and delamination resistance are quantified for particular materials and process conditions. When blades do fail, the likely cause is fatigue in the structural detail areas or at major flaws. The program is focused strongly on these issues in addition to standard laminates. Structural detail tests allow evaluation of various blade materials options in the context of more realistic representations of blade structure than do the standard test methods. Types of structural details addressed in this report include ply drops used in thickness tapering, and adhesive joints, each tested over a range of fatigue loading conditions. Ply drop studies were in two areas: (1) a combined experimental and finite element study of basic ply drop delamination parameters for glass and carbon prepreg laminates, and (2) the development of a complex structured resin-infused coupon including ply drops, for comparison studies of various resins, fabrics and pry drop thicknesses. Adhesive joint tests using typical blade adhesives included both generic testing of materials parameters using a notched-lap-shear test geometry developed in this study, and also a series of simulated blade web joint geometries fabricated by an industry partner.
The purpose of this report is to describe the methods commonly used to measure heat flux in fire applications at Sandia National Laboratories in both hydrocarbon (JP-8 jet fuel, diesel fuel, etc.) and propellant fires. Because these environments are very severe, many commercially available heat flux gauges do not survive the test, so alternative methods had to be developed. Specially built sensors include 'calorimeters' that use a temperature measurement to infer heat flux by use of a model (heat balance on the sensing surface) or by using an inverse heat conduction method. These specialty-built sensors are made rugged so they will survive the environment, so are not optimally designed for ease of use or accuracy. Other methods include radiometers, co-axial thermocouples, directional flame thermometers (DFTs), Sandia 'heat flux gauges', transpiration radiometers, and transverse Seebeck coefficient heat flux gauges. Typical applications are described and pros and cons of each method are listed.
For the U.S. Nuclear Regulatory Commission (NRC) Extremely Low Probability of Rupture (xLPR) pilot study, Sandia National Laboratories (SNL) was tasked to develop and evaluate a probabilistic framework using a commercial software package for Version 1.0 of the xLPR Code. Version 1.0 of the xLPR code is focused assessing the probability of rupture due to primary water stress corrosion cracking in dissimilar metal welds in pressurizer surge nozzles. Future versions of this framework will expand the capabilities to other cracking mechanisms, and other piping systems for both pressurized water reactors and boiling water reactors. The goal of the pilot study project is to plan the xLPR framework transition from Version 1.0 to Version 2.0; hence the initial Version 1.0 framework and code development will be used to define the requirements for Version 2.0. The software documented in this report has been developed and tested solely for this purpose. This framework and demonstration problem will be used to evaluate the commercial software's capabilities and applicability for use in creating the final version of the xLPR framework. This report details the design, system requirements, and the steps necessary to use the commercial-code based xLPR framework developed by SNL.
Sandia National Laboratories (SNL) conducts pioneering research and development in Micro-Electro-Mechanical Systems (MEMS) and solar cell research. This dissertation project combines these two areas to create ultra-thin small-form-factor crystalline silicon (c-Si) solar cells. These miniature solar cells create a new class of photovoltaics with potentially novel applications and benefits such as dramatic reductions in cost, weight and material usage. At the beginning of the project, unusually low efficiencies were obtained in the research group. The intention of this research was thus to investigate the main causes of the low efficiencies through simulation, design, fabrication, and characterization. Commercial simulation tools were used to find the main causes of low efficiency. Once the causes were identified, the results were used to create improved designs and build new devices. In the simulations, parameters were varied to see the effect on the performance. The researched parameters were: resistance, wafer lifetime, contact separation, implant characteristics (size, dosage, energy, ratio between the species), contact size, substrate thickness, surface recombination, and light concentration. Out of these parameters, it was revealed that a high quality surface passivation was the most important for obtaining higher performing cells. Therefore, several approaches for enhancing the passivation were tried, characterized, and tested on cells. In addition, a methodology to contact and test the performance of all the cells presented in the dissertation under calibrated light was created. Also, next generation cells that could incorporate all the optimized layers including the passivation was designed, built, and tested. In conclusion, through this investigation, solar cells that incorporate optimized designs and passivation schemes for ultrathin solar cells were created for the first time. Through the application of the methods discussed in this document, the efficiency of the solar cells increased from below 1% to 15% in Microsystems Enabled Photovoltaic (MEPV) devices.
We have developed a high sensitivity (<5 fTesla/{radical}Hz), fiber-optically coupled magnetometer to detect magnetic fields produced by the human brain. This is the first demonstration of a noncryogenic sensor that could replace cryogenic superconducting quantum interference device (SQUID) magnetometers in magnetoencephalography (MEG) and is an important advance in realizing cost-effective MEG. Within the sensor, a rubidium vapor is optically pumped with 795 laser light while field-induced optical rotations are measured with 780 nm laser light. Both beams share a single optical axis to maximize simplicity and compactness. In collaboration with neuroscientists at The Mind Research Network in Albuquerque, NM, the evoked responses resulting from median nerve and auditory stimulation were recorded with the atomic magnetometer and a commercial SQUID-based MEG system with signals comparing favorably. Multi-sensor operation has been demonstrated with two AMs placed on opposite sides of the head. Straightforward miniaturization would enable high-density sensor arrays for whole-head magnetoencephalography.
