In the field of scientific computing there are many specialized programs designed for specific applications in areas such as biology, chemistry, and physics. These applications are often very powerful and extraordinarily useful in their respective domains. However, some suffer from a common problem: a non-intuitive, poorly-designed user interface. The purpose of Optika is to address this problem and provide a simple, viable solution. Using only a list of parameters passed to it, Optika can dynamically generate a GUI. This allows the user to specify parameters values in a fashion that is much more intuitive than the traditional 'input decks' used by some parameterized scientific applications. By leveraging the power of Optika, these scientific applications will become more accessible and thus allow their designers to reach a much wider audience while requiring minimal extra development effort.
Hydrocarbon foams are commonly used in HEDP experiments, and are subject to shock compression from tens to hundreds of GPa. Modeling foams is challenging due to the heterogeneous character of the foam. A quantitative understanding of foams under strong dynamic compression is sought. We use Sandia's ALEGRA-MHD code to simulate 3D mesoscale models of pure poly(4-methyl-1-petene) (PMP) foams. We employ two models of the initial polymer-void structure of the foam and analyze the statistical properties of the initial and shocked states. We compare the simulations to multi-Mbar shock experiments at various initial foam densities and flyer impact velocities. Scatter in the experimental data may be a consequence of the initial foam inhomogeneity. We compare the statistical properties the simulations with the scatter in the experimental data.
Article 690.11 in the 2011 National Electrical Code{reg_sign} (NEC{reg_sign}) requires new photovoltaic (PV) systems on or penetrating a building to include a listed arc fault protection device. Currently there is little experimental or empirical research into the behavior of the arcing frequencies through PV components despite the potential for modules and other PV components to filter or attenuate arcing signatures that could render the arc detector ineffective. To model AC arcing signal propagation along PV strings, the well-studied DC diode models were found to inadequately capture the behavior of high frequency arcing signals. Instead dynamic equivalent circuit models of PV modules were required to describe the impedance for alternating currents in modules. The nonlinearities present in PV cells resulting from irradiance, temperature, frequency, and bias voltage variations make modeling these systems challenging. Linearized dynamic equivalent circuits were created for multiple PV module manufacturers and module technologies. The equivalent resistances and capacitances for the modules were determined using impedance spectroscopy with no bias voltage and no irradiance. The equivalent circuit model was employed to evaluate modules having irradiance conditions that could not be measured directly with the instrumentation. Although there was a wide range of circuit component values, the complex impedance model does not predict filtering of arc fault frequencies in PV strings for any irradiance level. Experimental results with no irradiance agree with the model and show nearly no attenuation for 1 Hz to 100 kHz input frequencies.
A new experimental technique to measure material shear strength at high pressures has been developed for use on magnetohydrodynamic (MHD) drive pulsed power platforms. The technique is referred to as Magnetically-Applied Pressure-Shear (MAPS). By applying an external static magnetic field to the sample region, the MHD drive directly induces a shear stress wave in addition to the usual longitudinal stress wave. Strength is probed by passing this shear wave through a sample material where the transmissible shear stress is limited to the sample strength. The magnitude of the transmitted shear wave is measured via a transverse velocity interferometry system (VISAR) from which the sample strength is determined. The strength of materials is defined as the ability of a material to sustain deviatoric (shear) stresses. Strength is an important aspect of the response of materials subjected to compression to high pressure. Beyond the elastic response, material strength will govern at what pressure and to what extent a material will plastically deform. The MAPS technique cleverly exploits the property that, for a von Mises yield criterion at a given longitudinal stress, the maximum amplitude shear wave that can be transmitted is limited by the strength at that stress level. Successful fielding of MAPS experiments to measure shear stresses relies upon correct numerical simulation of the experiment. Complex wave interactions among forward and reflected longitudinal and shear waves, as well as the advancing magnetic diffusion front of the MHD drive, can make the design of the experiment complicated. Careful consideration must be given to driver, sample, and anvil materials; to the thicknesses of the driver, sample and anvil layers; as well as to the timing of the interacting waves. This paper will present and analyze the 2D MHD simulations used to design the MAPS experiments. The MAPS experiments are modeled using Sandia's ALEGRA-MHD simulation code. ALEGRA-MHD is an operator-split, multi-physics, multi-material, arbitrary lagrangian-eulerian code developed to model magnetic implosion, ceramic fracture, and electromagnetic launch. We will detail the numerical investigations into MHD shear generation, longitudinal and shear stress coupling, timing of wave interactions, and transmission of shear at material interfaces.
Radiation portals normally incorporate a dedicated neutron counter and a gamma-ray detector with at least some spectroscopic capability. This paper describes the design and presents characterization data for a detection system called PVT-NG, which uses large polyvinyl toluene (PVT) detectors to monitor both types of radiation. The detector material is surrounded by polyvinyl chloride (PVC), which emits high-energy gamma rays following neutron capture reactions. Assessments based on high-energy gamma rays are well suited for the detection of neutron sources, particularly in border security applications, because few isotopes in the normal stream of commerce have significant gamma ray yields above 3 MeV. Therefore, an increased count rate for high-energy gamma rays is a strong indicator for the presence of a neutron source. The sensitivity of the PVT-NG sensor to bare {sup 252}Cf is 1.9 counts per second per nanogram (cps/ng) and the sensitivity for {sup 252}Cf surrounded by 2.5 cm of polyethylene is 2.3 cps/ng. The PVT-NG sensor is a proof-of-principal sensor that was not fully optimized. The neutron detector sensitivity could be improved, for instance, by using additional moderator. The PVT-NG detectors and associated electronics are designed to provide improved resolution, gain stability, and performance at high-count rates relative to PVT detectors in typical radiation portals. As well as addressing the needs for neutron detection, these characteristics are also desirable for analysis of the gamma-ray spectra. Accurate isotope identification results were obtained despite the common impression that the absence of photopeaks makes data collected by PVT detectors unsuitable for spectroscopic analysis. The PVT detectors in the PVT-NG unit are used for both gamma-ray and neutron detection, so the sensitive volume exceeds the volume of the detection elements in portals that use dedicated components to detect each type of radiation.
At the end of his life, Stephen Jay Kline, longtime professor of mechanical engineering at Stanford University, completed a book on how to address complex systems. The title of the book is 'Conceptual Foundations of Multi-Disciplinary Thinking' (1995), but the topic of the book is systems. Kline first establishes certain limits that are characteristic of our conscious minds. Kline then establishes a complexity measure for systems and uses that complexity measure to develop a hierarchy of systems. Kline then argues that our minds, due to their characteristic limitations, are unable to model the complex systems in that hierarchy. Computers are of no help to us here. Our attempts at modeling these complex systems are based on the way we successfully model some simple systems, in particular, 'inert, naturally-occurring' objects and processes, such as what is the focus of physics. But complex systems overwhelm such attempts. As a result, the best we can do in working with these complex systems is to use a heuristic, what Kline calls the 'Guideline for Complex Systems.' Kline documents the problems that have developed due to 'oversimple' system models and from the inappropriate application of a system model from one domain to another. One prominent such problem is the Procrustean attempt to make the disciplines that deal with complex systems be 'physics-like.' Physics deals with simple systems, not complex ones, using Kline's complexity measure. The models that physics has developed are inappropriate for complex systems. Kline documents a number of the wasteful and dangerous fallacies of this type.
The following report presents an assessment of existing capabilities in Sierra/Fuego applied to modeling several aspects of grid-to-rod-fretting (GTRF) including: fluid dynamics, heat transfer, and fluid-structure interaction. We compare the results of a number of Fuego simulations with relevant sources in the literature to evaluate the accuracy, efficiency, and robustness of using Fuego to model the aforementioned aspects. Comparisons between flow domains that include the full fuel rod length vs. a subsection of the domain near the spacer show that tremendous efficiency gains can be obtained by truncating the domain without loss of accuracy. Thermal analysis reveals the extent to which heat transfer from the fuel rods to the coolant is improved by the swirling flow created by the mixing vanes. Lastly, coupled fluid-structure interaction analysis shows that the vibrational modes of the fuel rods filter out high frequency turbulent pressure fluctuations. In general, these results allude to interesting phenomena for which further investigation could be quite fruitful.
The intense magnetic field generated by the Z accelerator at Sandia National Laboratories is used as a pressure source for material science studies. A current of ∼20 MA can be delivered to the loads used in experiments on a time scale of ∼100-600 ns. Magnetic fields (pressures) exceeding 1200 T (600 GPa) have been produced in planar configurations. In one application we have developed, the magnetic pressure launches a flyer plate to ultra-high velocity in a plate impact experiment; equation of state data is obtained on the Hugoniot of a material that is shock compressed to multi-megabar pressure. This capability has been enhanced by the recent development of a planar stripline configuration that increases the magnetic pressure for a given current. Furthermore, the cross sectional area of a stripline flyer plate is larger than in previous coaxial loads; this improves the planarity of the flyer thereby reducing measurement uncertainty. Results of experiments and multi-dimensional magneto hydrodynamic (MHD) simulation are presented for ultra-high velocity aluminum and copper flyer plates. Aluminum flyer plates with dimensions ∼25 mm by ∼13 mm by ∼1 mm have been launched to velocities up to ∼45 km/s; for copper the peak velocity is ∼22 km/s. The significance of these results is that part of the flyer material remains solid at impact with the target; an accomplishment that is made possible by shaping the dynamic pressure (current) ramp so that the flyer compresses quasi-isentropically (i.e., shocklessly) during acceleration.
The SCREAMER simulation code is widely used at Sandia National Laboratories for designing and simulating pulsed power accelerator experiments on super power accelerators. A preliminary parameter study of Z with a magnetic switching retrofit illustrates the utility of the automating script for optimizing pulsed power designs. SCREAMER is a circuit based code commonly used in pulsed-power design and requires numerous iterations to find optimal configurations. System optimization using simulations like SCREAMER is by nature inefficient and incomplete when done manually. This is especially the case when the system has many interactive elements whose emergent effects may be unforeseeable and complicated. For increased completeness, efficiency and robustness, investigators should probe a suitably confined parameter space using deterministic, genetic, cultural, ant-colony algorithms or other computational intelligence methods. I have developed SAE2 - a user-friendly, deterministic script that automates the search for optima of pulsed-power designs with SCREAMER. This manual demonstrates how to make input decks for SAE2 and optimize any pulsed-power design that can be modeled using SCREAMER. Application of SAE2 to magnetic switching on model of a potential Z refurbishment illustrates the power of SAE2. With respect to the manual optimization, the automated optimization resulted in 5% greater peak current (10% greater energy) and a 25% increase in safety factor for the most highly stressed element.
This document is a reference guide to the Xyce Parallel Electronic Simulator, and is a companion document to the Xyce Users Guide. The focus of this document is (to the extent possible) exhaustively list device parameters, solver options, parser options, and other usage details of Xyce. This document is not intended to be a tutorial. Users who are new to circuit simulation are better served by the Xyce Users Guide. The Xyce Parallel Electronic Simulator has been written to support, in a rigorous manner, the simulation needs of the Sandia National Laboratories electrical designers. It is targeted specifically to run on large-scale parallel computing platforms but also runs well on a variety of architectures including single processor workstations. It also aims to support a variety of devices and models specific to Sandia needs. This document is intended to complement the Xyce Users Guide. It contains comprehensive, detailed information about a number of topics pertinent to the usage of Xyce. Included in this document is a netlist reference for the input-file commands and elements supported within Xyce; a command line reference, which describes the available command line arguments for Xyce; and quick-references for users of other circuit codes, such as Orcad's PSpice and Sandia's ChileSPICE.
This manual describes the use of the Xyce Parallel Electronic Simulator. Xyce has been designed as a SPICE-compatible, high-performance analog circuit simulator, and has been written to support the simulation needs of the Sandia National Laboratories electrical designers. This development has focused on improving capability over the current state-of-the-art in the following areas: (1) Capability to solve extremely large circuit problems by supporting large-scale parallel computing platforms (up to thousands of processors). Note that this includes support for most popular parallel and serial computers; (2) Improved performance for all numerical kernels (e.g., time integrator, nonlinear and linear solvers) through state-of-the-art algorithms and novel techniques. (3) Device models which are specifically tailored to meet Sandia's needs, including some radiation-aware devices (for Sandia users only); and (4) Object-oriented code design and implementation using modern coding practices that ensure that the Xyce Parallel Electronic Simulator will be maintainable and extensible far into the future. Xyce is a parallel code in the most general sense of the phrase - a message passing parallel implementation - which allows it to run efficiently on the widest possible number of computing platforms. These include serial, shared-memory and distributed-memory parallel as well as heterogeneous platforms. Careful attention has been paid to the specific nature of circuit-simulation problems to ensure that optimal parallel efficiency is achieved as the number of processors grows. The development of Xyce provides a platform for computational research and development aimed specifically at the needs of the Laboratory. With Xyce, Sandia has an 'in-house' capability with which both new electrical (e.g., device model development) and algorithmic (e.g., faster time-integration methods, parallel solver algorithms) research and development can be performed. As a result, Xyce is a unique electrical simulation capability, designed to meet the unique needs of the laboratory.
This report describes the supercritical carbon dioxide (S-CO{sub 2}) direct cycle gas fast reactor (SC-GFR) concept. The SC-GFR reactor concept was developed to determine the feasibility of a right size reactor (RSR) type concept using S-CO{sub 2} as the working fluid in a direct cycle fast reactor. Scoping analyses were performed for a 200 to 400 MWth reactor and an S-CO{sub 2} Brayton cycle. Although a significant amount of work is still required, this type of reactor concept maintains some potentially significant advantages over ideal gas-cooled systems and liquid metal-cooled systems. The analyses presented in this report show that a relatively small long-life reactor core could be developed that maintains decay heat removal by natural circulation. The concept is based largely on the Advanced Gas Reactor (AGR) commercial power plants operated in the United Kingdom and other GFR concepts.
The behavior of carbon fiber aircraft composites was studied in adverse thermal environments. The effects of resin composition and fiber orientation were measured in two test configurations: 102 by 127 millimeter (mm) test coupons were irradiated at approximately 22.5 kW/m{sup 2} to measure thermal response, and 102 by 254 mm test coupons were irradiated at approximately 30.7 kW/m{sup 2} to characterize piloted flame spread in the vertically upward direction. Carbon-fiber composite materials with epoxy and bismaleimide resins, and uni-directional and woven fiber orientations, were tested. Bismaleimide samples produced less smoke, and were more resistant to flame spread, as expected for high temperature thermoset resins with characteristically lower heat release rates. All materials lost approximately 20-25% of their mass regardless of resin type, fiber orientation, or test configuration. Woven fiber composites displayed localized smoke jetting whereas uni-directional composites developed cracks parallel to the fibers from which smoke and flames emanated. Swelling and delamination were observed with volumetric expansion on the order of 100% to 200%. The purpose of this work was to provide validation data for SNL's foundational thermal and combustion modeling capabilities.
This report compiles 3-D finite element analyses performed to evaluate the stability of Strategic Petroleum Reserve (SPR) caverns over multiple leach cycles. When oil is withdrawn from a cavern in salt using freshwater, the cavern enlarges. As a result, the pillar separating caverns in the SPR fields is reduced over time due to usage of the reserve. The enlarged cavern diameters and smaller pillars reduce underground stability. Advances in geomechanics modeling enable the allowable pillar to diameter ratio (P/D) to be defined. Prior to such modeling capabilities, the allowable P/D was established as 1.78 based on some very limited experience in other cavern fields. While appropriate for 1980, the ratio conservatively limits the allowable number of oil drawdowns and hence limits the overall utility and life of the SPR cavern field. Analyses from all four cavern fields are evaluated along with operating experience gained over the past 30 years to define a new P/D for the reserve. A new ratio of 1.0 is recommended. This ratio is applicable only to existing SPR caverns.