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.
Implementation of the Progressive Hedging algorithm written in Julia utilizing JuMP. This is used to solve multistage optimization problems (e.g. Stochastic programs or capacity expansion models). The user uses the StructJuMP extension of JuMP to create the model and build up the scenario tree.
Time-to-solution for structure-based screening of massive chemical databases for COVID-19 drug discovery has been decreased by an order of magnitude, and a virtual laboratory has been deployed at scale on up to 27,612 GPUs on the Summit supercomputer, allowing an average molecular docking of 19,028 compounds per second. Over one billion compounds were docked to two SARS-CoV-2 protein structures with full optimization of ligand position and 20 poses per docking, each in under 24 hours. GPU acceleration and high-throughput optimizations of the docking program produced 350× mean speedup over the CPU version (50× speedup per node). GPU acceleration of both feature calculation for machine-learning based scoring and distributed database queries reduced processing of the 2.4 TB output by orders of magnitude. The resulting 50× speedup for the full pipeline reduces an initial 43 day runtime to 21 hours per protein for providing high-scoring compounds to experimental collaborators for validation assays.
The development of heme oxygenase (HO) inhibitors is critical in dissecting and understanding the HO system and for potential therapeutic applications. We have established a program to design and optimize HO inhibitors using structure-activity relationships in conjunction with X-ray crystallographic analyses. One of our previous complex crystal structures revealed a putative secondary hydrophobic binding pocket which could be exploited for a new design strategy by introducing a functional group that would fit into this potential site. To test this hypothesis and gain further insights into the structural basis of inhibitor binding, we have synthesized and characterized 1-(1H-imidazol-1-yl)-4,4-diphenyl-2-butanone (QC-308). Using a carbon monoxide (CO) formation assay on rat spleen microsomes, the compound was found to be {approx}15 times more potent (IC{sub 50} = 0.27{+-}0.07 {mu}M) than its monophenyl analogue, which is already a potent compound in its own right (QC-65; IC{sub 50} = 4.0{+-}1.8 {mu}M). The crystal structure of hHO-1 with QC-308 revealed that the second phenyl group in the western region of the compound is indeed accommodated by a definitive secondary proximal hydrophobic pocket. Thus, the two phenyl moieties are each stabilized by distinct hydrophobic pockets. This 'double-clamp' binding offers additional inhibitor stabilization and provides a new route for improvement of human heme oxygenase inhibitors.
Many advanced technical tools are available to prevent attacks on national infrastructure. Nevertheless, while traditional analyses of security problems have succeeded in producing good technical solutions, they have often ignored the human factor integral to these problems. Human attackers (who may be individuals or state-level attackers) expend substantial effort to breach security because they have the incentive for doing so. People involved in implementing security follow individual incentives, which need not align with global security concerns; consequently, desired security solutions are often implemented poorly, or not at all. This complex interplay between individual incentives and global (organizational and/or national) goals can be modeled and analyzed using game theoretic techniques. By analyzing not only what is possible, but also what is motivated, a holistic approach to security problems can be developed, informing policy and providing tools to policy makers. We study game theoretic models that unify several current incentive-based approaches to security, and develop simulation-based and mathematical optimization methods for analyzing such models that exploit the high-performance computing capabilities at Sandia. Our first model studies security in interdependent settings, offering a scalable local search heuristic to approximate optimal security decisions in general, and a linear programming approach, coupled with simulations of consequences, to optimally compute security in an important special case. Our second class of models addresses security patrolling problems when an adversary gets to observe the patrol location. We present a general framework, based on stochastic games, for computing optimal security policies in such settings, and present more scalable tools that apply in the important special cases. Our third contribution is a model of security that involves many defenders, but only models non-adaptive attackers (or natural disasters, inadvertent errors, etc). In this model, we demonstrate that the security decisions of many players result in global security configuration that is not very far from optimal, and is much more resilient to environment changes that an optimal solution. This positive effect dissipates, however, when the number of decision makers becomes too large.
LANL has developed a process to generate a progressive family of system models for a fissile solution system. This family includes a dynamic system simulation comprised of coupled nonlinear differential equations describing the time evolution of the system. Neutron kinetics, radiolytic gas generation and transport, and core thermal hydraulics are included in the DSS. Extensions to explicit operation of cooling loops and radiolytic gas handling are embedded in these systems as is a stability model. The DSS may then be converted to an implementation in Visual Studio to provide a design team the ability to rapidly estimate system performance impacts from a variety of design decisions. This provides a method to assist in optimization of the system design. Once design has been generated in some detail the C++ version of the system model may then be implemented in a LabVIEW user interface to evaluate operator controls and instrumentation and operator recognition and response to off-normal events. Taken as a set of system models the DSS, Visual Studio, and LabVIEW progression provides a comprehensive set of design support tools.
From beneath the surface of the earth, we currently obtain about 80-percent of the energy our nation consumes each year. In the future we have the potential to generate billions of watts of electrical power from clean, green, geothermal energy sources. Our planet’s subsurface can also serve as a reservoir for storing energy produced from intermittent sources such as wind and solar, and it could provide safe, long-term storage of excess carbon dioxide, energy waste products and other hazardous materials. However, it is impossible to underestimate the complexities of the subsurface world. These complexities challenge our ability to acquire the scientific knowledge needed for the efficient and safe exploitation of its resources. To more effectively harness subsurface resources while mitigating the impacts of developing and using these resources, the U.S. Department of Energy established SubTER – the Subsurface Technology and Engineering RD&D Crosscut team. This DOE multi-office team engaged scientists and engineers from the national laboratories to assess and make recommendations for improving energy-related subsurface engineering. The SubTER team produced a plan with the overall objective of “adaptive control of subsurface fractures and fluid flow.”This plan revolved around four core technological pillars—Intelligent Wellbore Systems that sustain the integrity of the wellbore environment; Subsurface Stress and Induced Seismicity programs that guide and optimize sustainable energy strategies while reducing the risks associated with subsurface injections; Permeability Manipulation studies that improve methods of enhancing, impeding and eliminating fluid flow; and New Subsurface Signals that transform our ability to see into and characterize subsurface systems. The SubTER team developed an extensive R&D plan for advancing technologies within these four core pillars and also identified several areas where new technologies would require additional basic research. In response, the Office of Science, through its Office of Basic Energy Science (BES), convened a roundtable consisting of 15 national lab, university and industry geoscience experts to brainstorm basic research areas that underpin the SubTER goals but are currently underrepresented in the BES research portfolio. Held in Germantown, Maryland on May 22, 2015, the round-table participants developed a basic research agenda that is detailed in this report. Highlights include the following: -A grand challenge calling for advanced imaging of stress and geological processes to help understand how stresses and chemical substances are distributed in the subsurface—knowledge that is critical to all aspects of subsurface engineering; -A priority research direction aimed at achieving control of fluid flow through fractured media; -A priority research direction aimed at better understanding how mechanical and geochemical perturbations to subsurface rock systems are coupled through fluid and mineral interactions; -A priority research direction aimed at studying the structure, permeability, reactivity and other properties of nanoporous rocks, like shale, which have become critical energy materials and exhibit important hallmarks of mesoscale materials; -A cross-cutting theme that would accelerate development of advanced computational methods to describe heterogeneous time-dependent geologic systems that could, among other potential benefits, provide new and vastly improved models of hydraulic fracturing and its environmental impacts; -A cross-cutting theme that would lead to the creation of “geo-architected materials” with controlled repeatable heterogeneity and structure that can be tested under a variety of thermal, hydraulic, chemical and mechanical conditions relevant to subsurface systems; -A cross-cutting theme calling for new laboratory studies on both natural and geo-architected subsurface materials that deploy advanced high-resolution 3D imaging and chemical analysis methods to determine the ;rates and mechanisms of fluid-rock processes, and to test predictive models of such phenomena. Many of the key energy challenges of the future demand a greater understanding of the subsurface world in all of its complexity. This greater under- standing will improve the ability to control and manipulate the subsurface world in ways that will benefit both the economy and the environment. This report provides specific basic research pathways to address some of the most fundamental issues of energy-related subsurface engineering.
The goals of this project are to examine novel and existing techniques for proton transport; determine portability and performance on new hardware; and develop template GPU code that can be expanded.
To explore the new energy frontier, a new generation of particle accelerators is needed. Muon colliders are a promising alternative, if muon cooling can be made to work. Muons are 200 times heavier than electrons, so they produce less synchrotron radiation, and they behave like point particles. However, they have a short lifetime of 2.2 $\mathrm{\mu s}$ and the beam is more difficult to cool than an electron beam. The Muon Accelerator Program (MAP) was created to develop concepts and technologies required by a muon collider. An important effort has been made in the program to design and optimize a muon beam cooling system. The goal is to achieve the small beam emittance required by a muon collider. This work explores a final ionization cooling system using magnetic quadrupole lattices with a low enough $\beta^{\star} $ region to cool the beam to the required limit with available low Z absorbers.
This document is a notional technology implementation plan (TIP) for the development, testing, and qualification of a prototypic fuel element to support design and construction of a nuclear thermal propulsion (NTP) engine, specifically its pre-flight ground test. This TIP outlines a generic methodology for the progression from non-nuclear out-of-pile (OOP) testing through nuclear in-pile (IP) testing, at operational temperatures, flows, and specific powers, of an NTP fuel element in an existing test reactor. Subsequent post-irradiation examination (PIE) will occur in existing radiological facilities. Further, the methodology is intended to be nonspecific with respect to fuel types and irradiation or examination facilities. The goals of OOP and IP testing are to provide confidence in the operational performance of fuel system concepts and provide data to program leadership for system optimization and fuel down-selection. The test methodology, parameters, collected data, and analytical results from OOP, IP, and PIE will be documented for reference by the NTP operator and the appropriate regulatory and oversight authorities. Final full-scale integrated testing would be performed separately by the reactor operator as part of the preflight ground test.
The purpose of this project was to design and build a set of resonant cavity beam diagnostics including a monopole cavity, two dipole cavities, and a quadrupole mode cavity. An array of such diagnostics placed in an accelerator FODO lattice can be used to reconstruct beam emittance and envelope information on a pulse-by-pulse basis. In this project a more efficient quadrupole cavity shape was discovered which could provide the required sensitivity with just two cells vs. nine cells as originally envisioned. This new design was studied and optimized with RF simulation programs, followed by mechanical specifications and drawings, and fabrication of the diagnostics. The device was tested with a relativistic electron beam at Brookhaven National Laboratory ATF accelerator.
This report summarizes the work performed under the project "Linear Programming in Strongly Polynomial Time." Linear programming (LP) is a classic combinatorial optimization problem heavily used directly and as an enabling subroutine in integer programming (IP). Specifically IP is the same as LP except that some solution variables must take integer values (e.g. to represent yes/no decisions). Together LP and IP have many applications in resource allocation including general logistics, and infrastructure design and vulnerability analysis. The project was motivated by the PI's recent success developing methods to efficiently sample Voronoi vertices (essentially finding nearest neighbors in high-dimensional point sets) in arbitrary dimension. His method seems applicable to exploring the high-dimensional convex feasible space of an LP problem. Although the project did not provably find a strongly-polynomial algorithm, it explored multiple algorithm classes. The new medial simplex algorithms may still lead to solvers with improved provable complexity. We describe medial simplex algorithms and some relevant structural/complexity results. We also designed a novel parallel LP algorithm based on our geometric insights and implemented it in the Spoke-LP code. A major part of the computational step is many independent vector dot products. Our parallel algorithm distributes the problem constraints across processors. Current commercial and high-quality free LP solvers require all problem details to fit onto a single processor or multicore. Our new algorithm might enable the solution of problems too large for any current LP solvers. We describe our new algorithm, give preliminary proof-of-concept experiments, and describe a new generator for arbitrarily large LP instances.
This report is part of the collaboration with the Nuclear Criticality Safety Program and constitutes the “CED-2 report” of the experiments proposed in the SPRF/CX facility of SANDIA National Laboratories to contribute to the validation of Molybdenum in thermal energy spectrum. It supplements the study carried out in the CED-1 report by providing a more precise design of the proposed experiments and presents a complete evaluation of the experimental uncertainties. In the CED-1 report, corresponding to the preliminary design of the experiments, IRSN identified several promising configurations that could be implemented in the SPRF/CX facility of SNL. They use metallic molybdenum either in the form of foils between the UO 2 pellets inside BUCCX rods, or in the form of sleeves surrounding 7uPCX rods. In this CED-2 report, a new design using molybdenum sleeves surrounding BUCCX rods is also analyzed and the experimental uncertainties of all the proposed configurations are estimated. Given the constraints to carry out the experiments and the costs associated with the manufacturing of the foils, the experiments using molybdenum foils are discarded. Similarly, the design using molybdenum sleeves around BUCCX rods leads to significant manufacturing costs which could certainly be partially offset by use of existing grids but which would lead to slightly lower values of the integral sensitivities to the cross sections of molybdenum. Consequently, the design using molybdenum sleeves around 7uPCX rods is favored insofar as it presents higher $k_{eff}$ sensitivities to the molybdenum cross sections and that there is still possibility to reduce the number of sleeves to reduce the cost of experiments. Regarding feedback on nuclear data in support of criticality safety assessment, one of the advantages of the work described here is that the proposed experiments are not correlated with the MIRTE and FP experiments using molybdenum. In addition, they lead to higher sensitivity values in the thermal energy range than the existing experiments and partly “cover” the first resonance of the 95 Mo, which is a real advantage for improving knowledge on molybdenum cross sections.
The impetus for this work is the Proliferation Resistance Optimization (PRO-X) program, which was created by the National Nuclear Security Administration (NNSA) of the US Department of Energy (DOE) to provide a framework for developing reactor designs that minimize the production of special nuclear materials while maximizing performance for peaceful uses. The Argentine company INVAP and Argonne National Laboratory collaborated to share knowledge and expertise regarding the thermal analysis of research reactors. Working together, they designed a series of test problems based on the Open Pool Australian Light Water (OPAL) reactor, whose fuel is in the form of flat plates cooled by water flowing through narrow rectangular channels. ANL used their PLTEMP/ANL code. INVAP used their TERMIC code (and its new multi-plate version, TERMIC-MP). Then, INVAP and ANL compared the two sets of results. 1-D models were used during the collaboration since for research reactor analysis, largely 1-D thermal-hydraulics models are commonly used because they are simple to apply and technically defendable. The collaboration between the two organizations aims for improvements in thermal-hydraulic models. Improvements to methods and models can reduce modeling uncertainties and excessive conservatism, allowing greater reactor performance without reducing predicted safety margins. Therefore, these improvements potentially expand the design space, making additional designs feasible, which, in turn, may achieve better fuel utilization and proliferation resistance in the designs. This report provides justification for the lateral node size (stripe width) chosen for the 1-D computer models during the collaboration. The lateral heat conduction in OPAL fuel plates is evaluated using the two-dimensional conduction capability of PLTEMP/ANL to determine a representative hot-stripe width that can be conservatively used for one-dimensional steady-state thermal-hydraulic analysis. Basically, a representative hot-stripe width needs to be selected so that the actual heat flux averaged over this stripe width is greater than or equal to the maximum heat flux obtained by the two-dimensional analysis. Detailed results in this report show that a hot-stripe width of 8.125 mm can be conservatively used for the one-dimensional thermal-hydraulics analysis of all plates. Therefore, eight equal stripes (8.125 mm each) are recommended to be modeled in the MCNP and PLTEMP/ANL full core models of the OPAL reactor.
In this work, we propose a Generalized Disjunctive Programming (GDP) model that optimizes both long-term capacity planning (such as the number and size of dispatchable/renewable generators, batteries, and transmission lines) and hourly operation (such as on/off schedules of dispatchable generators, power output from each generator, and power flow) to maximize power system reliability while minimizing total cost and CO2 emissions.
This report describes a hypothetical facility for production of medical radioisotopes via activation under the Proliferation Resistance and Optimization (PRO-X) program. The facility uses neutron activation of non-special nuclear material (SNM) to produce the medical isotopes 131 I and 99 Mo at a throughput of 60 Ci/week of 131 I and 5 Ci/week of 99 Mo. The hypothetical design was carried out using a 10 MWt research reactor. The precursors used for the activation process were TeO2 for 131 I and MoO 3 for 99 Mo. The processes are performed in 3 hot cells used for target receipt, extraction, purification low specific activity (LSA) generator introduction, and packaging. A fourth hotcell is used for waste processing. The hot cell processing area takes up a footprint of 15.4 m 2 with the total footprint of the facility, including space for administrative offices, non-rad labs, quality assurance, and radiation buffer areas set at 763 m 2 . Waste is produced at a weekly rate of 257.8 g low activity solid waste and 8032.7 mL of low activity liquid waste, 8032 mL of which is water. This baseline hypothetical facility for production of medical isotopes via activation was then compared and contrasted to the hypothetical facility for production of medical isotopes via fission products to show the differences in approach for the two production modes. The two production modes had several highlighted differences including the overall facility and hot cell layout, the type and amount of waste produced by the respective facilities, and economic factors impacting production mode. Finally, a decision tree for which production mode might be more beneficial for an entrant into medical isotope production was developed based on the differences examined and the desired output of medical isotopes desired by the entrant.
Market power exercise in the electricity markets distorts market prices and diminishes social welfare. Many markets have implemented market power mitigation processes to eliminate the impact of such behavior. The design of mitigation mechanisms has a direct influence on investors' profitability and thus mid-/long-term resource adequacy. In order to evaluate the effectiveness of the existing market power mitigation mechanisms, this paper proposes a mitigation-aware strategic bidding model and studies the bidding strategies of the market participants under current practice. The proposed bidding model has a bilevel structure with strategic participant's profit maximization problem in the upper level and the dispatch problem for market operators in the lower level. In particular, the consideration of potential offer mitigation is incorporated as upper-level constraints based on the conduct and impact tests. This bilevel problem is reduced to a single-level mixed-integer linear program using the KKT optimality conditions, duality theory, and linearization. Numerical results illustrate how a strategic player can exercise market power to achieve a higher profit even under the current market power mitigation process and we analyze the social impact that the market power exercise results.
The Cherenkov Telescope Array (CTA) is a ground-based observatory for very high-energy (10 GeV to 100 TeV) gamma rays, planned for operation starting in 2018. It will be an array of dozens of optical telescopes, known as Atmospheric Cherenkov Telescopes (ACTs), of 8 m to 24 m diameter, deployed over an area of more than 1 square km, to detect flashes of Cherenkov light from showers initiated in the Earth's atmosphere by gamma rays. CTA will have improved angular resolution, a wider energy range, larger fields of view and an order of magnitude improvement in sensitivity over current ACT arrays such as H.E.S.S., MAGIC and VERITAS. Several institutions have proposed a research and development program to eventually contribute 36 medium-sized telescopes (9 m to 12 m diameter) to CTA to enhance and optimize its science performance. The program aims to construct a prototype of an innovative, Schwarzschild-Couder telescope (SCT) design that will allow much smaller and less expensive cameras and much larger fields of view than conventional Davies-Cotton designs, and will also include design and testing of camera electronics for the necessary advances in performance, reliability and cost. We report on the progress of the mid-sized SCT development program.
This briefing describes recent work undertaken on the Sparrow Project and results of this work. It describes experiments comparing the use of Genie with 2 classes with 3 classes for the problem of ship delineation. It also describes some preliminary work in the area of the optimization of segmentation techniques.