SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
The vast majority of nuclear weapon yield is generated in a regime of extreme temperature, pressure, and density. Without explosive nuclear testing, HED experiments on ICF facilities are the only way to approach and probe these extreme conditions that affect nuclear weapon performance.
Tracer gases, whether they are chemical or isotopic in nature, are useful tools in examining the flow and transport of gaseous or volatile species in the underground. One application is using detection of short-lived argon and xenon radionuclides to monitor for underground nuclear explosions. However, even chemically inert species, such as the noble gases, have bene observed to exhibit non-conservative behavior when flowing through porous media containing certain materials, such as zeolites, due to gas adsorption processes. This report details the model developed, implemented, and tested in the open source and massively parallel subsurface flow and transport simulator PFLOTRAN for future use in modeling the transport of adsorbing tracer gases.
Moment tensors (MTs) have long been used in earthquake and explosion source analysis, and there has been a renewed interest in how they can inform us about the seismic source, particularly in the geophysical monitoring community due to its application in event identification and yield analysis. However, parameter uncertainties in seismic MT inversion are rarely available. The inverse procedure often does not quantify MT model errors such as event location, data noise and Earth model that are essential for estimating solution robustness. To address this need, we propose to adopt the Bayesian probabilistic framework to incorporate uncertainties in MT inversions. In this study, we present the theoretical background of a probabilistic Bayesian framework for MT inversion accounting for model and measurements errors and illustrate the implementation of the method using a synthetic example.
B225 is operated by the Weapons and Complex Integration (WCI) Directorate and is classified as a Low Hazard facility for radiological hazards and Light Science Industry (LSI) for chemical, radiation generating devices (RGDs) and industrial hazards. Explosives and biological hazards are not authorized for B225. B225 is located in the southwest quadrant of Lawrence Livermore National Laboratory (LLNL) less than 600 meters from the nearest site boundary. B225 is a single-story concrete tilt-up building with a steel deck roof. The primary mission of this facility is to apply nondestructive evaluation (NDE) methods to materials, components, and assemblies. B225 also contains the polymer additive manufacturing laboratories. The floor plans for B225 are provided in Appendix A.
This Safety Basis Document (SBD) covers operations in the Weapons and Complex Integration (WCI) Building 225 (B225) facility. B225 is located in the southwestern quadrant of Lawrence Livermore National Laboratory (LLNL) and is operated and owned by WCI. This SBD satisfies the requirements of LLNL ES&H Manual Document 3.1, Nonnuclear Safety Basis Program (Ref. 2). B225 is classified as a Low Hazard Facility based on the inventory of radiological materials. Other hazards in B225: chemicals, radiation generating devices, and industrial hazards are managed at the Light Science and Industry (LSI) level. Biological and explosive hazards are not authorized in B225. There are four general Administrative Controls (ACs) identified as the B225 Operational Safety Requirements (OSRs). These administrative controls are deemed sufficient to ensure safe operations with respect to co-located workers and the public. Other safety features and administrative requirements are required by the LLNL ES&H Manual, including those to protect the facility workers and the environment, but these have not been selected as credited controls to protect the co-located workers or the public, which are the target populations of concern for the SBD.
The main purpose of the Carbonara experiment is to provide much-needed data to inform physics modeling of boron carbide (B 4 C) subjected to a sweeping high explosive (HE) detonation. The proposed experiment will include three separate shots, each employing a rectangular slab of HE (4” × 3” × 1”) sandwiched between two wedges of B 4 C, see Fig.1. In each case the top B 4 C wedge tapers from a height of 2.5 cm to 0.5 cm, while the bottom wedge tapers from 1.5 cm to 0.5 cm. The wedges are 4” in the direction of the proton beam, and 3” in the direction transverse to the beam. Three shots are currently planned, each one using a different HE: PBX 9501, PBX 9502, and PBX 9701. The HE will be initiated using a line wave generator (LWG) to preserve planar symmetry, in conjunction with a Comp. B booster to smooth out the profile for the initial condition on the main charge and to aid in initiation of the PBX 9502 shot. In each case the booster spans the entire length and height of the main HE slab and extends 0.5” in toward the LWG.
This project aimed to address the challenge of effectively welding dissimilar materials—high-strength steel and high-strength aluminum for creating lightweight, multi-material automotive components. For automotive companies, reducing weight of a vehicle is critical task regulated by the government to solve the issue of greenhouse gas emissions. Production of lightweight cars and trucks can be achieved by substitution of current all-steel structures with multi-material lightweight structures that include high strength-to-weight-ration materials such as high-strength steels, aluminum alloys, magnesium alloys, titanium alloys, coupled with lightweight designs. This requires dissimilar metal welding, which is challenging for state of the art joining processes such as resistance spot welding. The cycle of melting-cooling-freezing during traditional welding that can easily ruin the designed outstanding properties of the advanced base metals, such as aluminum alloys, making the welded area much weaker than the base metals. To weld two different metals with great difference in melting points, such as aluminum and steel, it’s even more difficult or impossible because of the formation of brittle intermetallic compounds at the welded interface. In this project, a novel welding method, developed at OSU, was selected for validation and development. This novel technology enables welding by impact without melting and proves to be robust to join various dissimilar lightweight metals. Termed as vaporizing foil actuator welding or VFAW, the technology uses a thin aluminum foil that is rapidly vaporized by a high current pulse to produce an explosive-like pressure pulse to drive one metallic piece into another at the high speed required for impact welding. This project entailed development of the early-stage welding technology in terms of (a) the consumables, the welding apparatus and the power sources, (b) coupon scale screening of many material combinations including corrosion studies, (c) computational modeling and design of the welded interface as well as of the multi-material prototype component, and (d) mechanical testing for strength and durability at coupon scale and to a certain extent the prototype scale. The all-steel engine cradle of 2016 Chevrolet Cruze was chosen as the baseline prototype component. The target set for the project was to demonstrate a 20% weight reduction at a cost premium of less than $\$ $5/lb saved without compromising on baseline mechanical properties. At project completion, a 12% lighter prototype component was demonstrated with an estimated cost premium of $\$ $9.8/lb saved. Besides prototype level demonstration of the technology, this project also enabled elevation of the technology’s readiness level to where a hydraulically actuated welding head was developed and made ready for deployment at a research and development facility for Tier 1 automotive supplier.
Understanding the physics of collisionless, magnetized plasmas is key to unlocking a number of important problems in space and astrophysics. Magnetic reconnection is a fundamental mechanism which can explosively release stored magnetic energy, convert it to plasma heat and flows, and accelerate particles. In this project we developed technique and conducted experiments and simulations of magnetic reconnection in laser produces plasmas. Through this project, we also developed experimental techniques to generate magnetized collisionless shocks, and we therefore expanded the scope of the project to develop these experiments. These projects make fundamental progress on understanding the dynamics of magnetized laser-produced plasmas. In additional to plasma astrophysics, these results have broad interest for laboratory high energy density plasmas, since magnetic fields can be used to confine and control plasmas in many contexts and may offer new capabilities toward laboratory ignited plasmas.
This report represents the 1st shot (in a series of 8) conducted on September 15, 2022. One 10 lb C4 charge (along with ~200g of Potassium Bromide (KBr)) was detonated inside 9920 North Pad Boom Box. Noise sampling was performed at several points on Site 9920 to characterize the noise mitigation provided by the block structure. This data will help inform safe locations for Members of the Workforce (MOWs) to be located during future testing with similar net explosive weights. During the test, all MOW/site visitors were bunkered inside Building 9926/Mobile Firing Control Point (MFCP) to prevent personnel exposure to any hazards associated with the testing
Tritium breeding is fundamentally required for a sustainable fusion fuel cycle, yet the technological readiness of blanket technology lags far behind other fusion systems. Breeder concepts are divided into solid and liquid media, where solid breeders typically rely on a sweep gas, such as helium, to carry away tritium from lithium containing ceramic materials, and liquid breeders produce tritium from lithium containing eutectics (e.g., PbLi) or molten salts (FLiBe). In each case, tritium must be harvested from the breeding medium. A promising method for tritium extraction is through a vacuum permeator, in which a concentration gradient from the tritium- containing fluid promotes diffusion through a membrane with high hydrogen permeability to the vacuum. This technology has been demonstrated for hydrogen gas systems using Pd and PdAg permeators, but relatively little work has been done to test tritium extraction from PbLi. A Tritium Extraction eXperimental (TEX) loop is being designed to test tritium extraction in a vacuum permeator configuration. The system design is such that it will allow the testing of tritium extraction from both helium and PbLi. A phased approach is being taken that will allow testing of small specimens for fundamental permeation measurements, to multi-meter component testing at near-prototypic conditions. A molten PbLi loop is challenging due to the toxic and explosive nature of Pb and Li, respectively, radiological concerns by introducing tritium, and high temperatures involved in such a system. In addition, PbLi corrosion is a significant issue at high temperatures (>400C). The TEX system will not employ a neutron source for volumetric production of tritium. Herein we present the design and methods used for 1) pumping PbLi, 2) introducing deuterium or tritium into the PbLi, 3) quantifying the amount of deuterium in the loop, 4) extracting deuterium and tritium, and 5) quantifying the total amount of extracted deuterium or tritium from the permeator. In addition, the safety design for operating such a system will be discussed.
Material identification (ID) using radiography is a problem that has implicitly existed since the discovery of the X-ray. Since 9-11, this problem has received renewed attention for homeland security applications of the generic, “what is in the box?” type1. Examples of items potentially in the box include: roadside bombs, nuclear weapons, contraband shipments and the like. In these cases, the “what” can include: explosives, projectiles, fissile material, or drugs. Here we investigate the ability of 60Co gamma-rays, 14MeV neutrons, and their combination to tease apart information about density, thickness and atomic number of items in the box. That information subset can be combined with other types of independent information (e.g. neutron and gammaspectroscopy, intelligence, photographs etc.) to assert a more complete picture of what, exactly, is in the box? That more complete information set can form the basis for actions including: remote detonation, evacuations, disablement, detention, search, arrest etc.
Solving the Euler equations is a fundamental component of simulating many physical phenomena ranging from high explosives to astrophysics. An equation of state (EOS) is a required piece that relates any two thermodynamic quantities to all other thermodynamic values. The presence of multiple materials within a control volume further complicates the solution requiring additional equations to describe the interaction of materials at a sub-grid level. Equations of state themselves can also come in many forms ranging from simple algebraic relations to more complicated differential equation models that describe material interactions over a broad range of physical conditions. In the latter case, the EOS is often pre-computed at a given set of grid points and provided in a tabular form where additional properties can be derived from the interpolation functions.
Enhanced heat transfer model has been developed due to the turbulence in the cracks and pores in the burning fractured HE. This turbulence then serves in the initialization of the turbulence in the product gas. Enhanced heat transfer is also modeled in the turbulent product gas. Numerical simulations of the onset of HE burning and transition to detonation have been performed with and without the turbulent heat transfer model. Computational results indicate that onset of chemical reactions and transition to detonation is induced by the shockwave formed by the pressure rise caused by the heating. In other words, processes similar to shock-to-detonation transition (SDT) could become the driving mechanism for HE explosion and detonation in the case of heating rather than the temperature alone.
Flash radiography of hydrodynamic experiments driven by high explosives is a well-known diagnostic technique in use at many laboratories. At Los Alamos, the Dual Axis Radiographic Hydrodynamic Test (DARHT) facility two linear induction electron accelerators (LIAs) make the bremsstrahlung radiographic source spots for point projection radiographs from orthogonal views. A new LIA, called Scorpius, is presently under development to advance this technology. To better understand electron-beam physics in these LIAs, numerical simulations are frequently performed with the objective of improving the radiography. At Los Alamos we frequently use the TRAK ray-trace and LSP particle-in-cell (PIC) codes to simulate the injector, and the XTR and LAMDA envelope/centroid codes along with LSP to simulate transport of the accelerated beam through the LIAs. The LIA simulations need the injected beam parameters as initial conditions for calculating beam transport and stability. The determination of these initial conditions is the topic of this note.
This report describes the second edition of the Pseudocap Strength models for porous rocks implemented in GEODYN material library. The first model was developed in 2007 and calibrated for concrete. Then, the model parameters were calibrated based on triaxial tests reported for limestones and sandstones of various porosities. In these models some key parameters were chosen as functions of the reference porosity,Φ. Since then multiple modifications were implemented in the model, therefore, it has been recalibrated for some common porous materials such as limestones, sandstones, alluvium, tuffs and granite. Two types of models are described in this repot. The first type (called Pseudocap Model or PM) is for rocks from a specific location. Parameters were calibrated for several specific geologic materials. The second type (called Generic Pseudocap Model or GPM) is useful for the sites where only basic information (rock type, porosity) is available. Generic models include built-in correlations between porosities and other mechanical properties observed for certain rock types. The models of both types were validated by comparing not only to quasi-static triaxial tests for these materials but also to shock Hugoniot data and spherical explosion data for some materials. All models were derived in the frame of isotropic plasticity. They are designed to be used in explicit finite element/difference codes. Tangent stiffness tensor is not provided but can be calculated numerically for the model to be used in implicit finite element codes. For an isotropic material the stress can be decomposed into volumetric and deviatoric parts. The volumetric part is modeled using an Equation of state (EOS) which calculates the pressure and the bulk sound speed as functions of the internal specific energy and density. Here a simple, Mie-Gruneisen EOS is presented, but tabulated EOS (LEOS) provided by the library can be used as well. The stress is limited by the yield surface which depends on three invariants of the stress tensor and specific internal energy. In addition, to capture the strain-rate dependence a simple multiplier is used for the yield surface which depends on the equivalent plastic strain rate. The failure surface (the ultimate yield, Y f , defined later) is chosen in the Hoek-Brown form, commonly used in rock mechanics. It includes measurable parameters such as Unconfined Compressive Strength (UCS) as well as scale parameters characterizing the quality of the rock such as GSI (Geologic Strength Index). Thus, even though the model is calibrated for small samples it offers a way to extrapolate the strength to the field scale using geological characterization of the rock mass. The model captures effects of brittle-ductile transition in rocks by introducing a cap multiplier to the yield function. The rate of dilatancy (bulking) is proportional to the slope of the yield surface affected by the cap. Therefore, it takes place only at low confinements when the pressure is less than the brittle-ductile transition pressure, P BD . On the contrary, the porous compaction takes place at pressures higher than P BD . The cap moves as the porosity is compacted or new porosity is generated due to dilatancy. The porous compaction is modeled using an evolution equation which includes deviatoric stress so that the onset of compaction corresponds to the cap surface. The model captures effects of shear-enhanced compaction which is an important for porous rocks. Section 2 describes the modeling framework and Section 3 presents the model calibration procedure. Section 4 compares experimental data for various rocks versus model predictions. The model parameters used for this comparison are given in Appendix. The files with material constants are available with the latest GEODYN material library distribution.
During this project, we assessed the percentage of events reported by the International Data Centre (IDC) that are candidates for evaluation using data fusion methods. We calculated the percentage of waveform (seismic, hydroacoustic, infrasound, SHI) events that may be associated with radionuclide (RN) detections and vice versa. For this task, data fusion means identifying RN detections that may be associated with SHI events. The task of producing a fused SHI and RN bulletin is included in the IDC operational manual. As of 2022, the production of a fused event bulletin remains incomplete, largely because methods for producing the bulletin are under development. Nonetheless, simple methods for data fusion have been developed and preliminary software tools are available at the IDC to associate SHI events and RN detections. This report is a first look at the percentage of SHI events that may be candidates for fusion with RN detections or bulletins. One year of IDC bulletin data (2021) is considered in this report. It is important to note that there were no known nuclear explosions during this time period. Hence, we expect that events during this time period would generate SHI and RN signals. If many SHI events are candidates for fusion with RN data, greater priority should be placed on establishing screening procedures before and/or during the fusion process. By contrast, if the number of candidate events is diminishingly small, then data fusion is not a driver for further development of event or detection screening at the IDC. LLNL reviewed the events reported in the Standard Event Bulletin (SEB) and determined which events are candidates for fusion with the Reviewed Radionuclide Reports (RRR). Much of the work involves establishing a metric or metrics for determining whether a Reviewed Event Bulletin (REB) event is a candidate for data fusion.
GnarlyX is a new hydrocode for direct numerical simulations of the microstructural behavior of high explosives at the mesoscale. We summarize the computational framework for multidimensional, Eulerian multi-material hydrodynamics coupled to EOS and hyperelastic, plastic constitutive models. We present 1D verification tests of multi-material only and combined multimaterial and strength capabilities with comparisons to exact solutions of shock states resulting from an incident shock impacting the material interface of PBX9502 and aluminum. We show that GnarlyX performs well in capturing the resulting shock waves in examining numerical convergence with exact solutions. In later work, we will summarize the thermomechanics and multi-dimensional, parallel computing capabilities in GnarlyX with multi-dimensional verification tests.
We present an engineering-oriented summary of statistical methods for estimation of the reliability (or equivalently, failure probability) of one-shot devices such as explosive detonators and other weapon components. Estimates may be given as single points or intervals, based on pass/fail tests, margin analysis, computational models, expert judgment, or a combination of these. We focus on highly reliable devices for which few or no failures are expected to occur in testing.