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.
Inspired by Slater's rules, we focus on developing computationally efficient and physically motivated contractions to reproduce the core-excited atomic orbitals with improved numerical accuracy of calculated core-electron binding energies (CEBEs).
FastCaloSim is a parameterized simulation of the particle energy response and of the energy distribution in the ATLAS calorimeter. It is a relatively small and self-contained package with massive inherent parallelism and captures the essence of GPU offloading via important operations like data transfer, memory initialization, floating point operations, and reduction. It was identified by the High Energy Physics Center for Computational Excellence project as a good testbed for evaluating the performance and ease of portability of programming models. In this paper, we will discuss the results of our evaluation of the porting process to Kokkos, SYCL, Alpaka, OpenMP and std::par (nvc++), and compare performance on NVIDIA, AMD and Intel GPUs, as well as multicore CPUs.
3D numerical simulations of the interaction of a powerful CO 2 laser with hydrogen jets demonstrating the role of ionization on the characteristics of induced wakes are presented. Simulations using SPACE, a parallel relativistic Particle-in-Cell code, are performed in support of the plasma wakefield accelerator experiments being conducted at the Brookhaven National Laboratory (BNL) Accelerator Test Facility (ATF). A novelty of the SPACE code is its set of efficient atomic physics algorithms that compute ionization and recombination rates on the grid and transfer them to particles. The influence of ionization on the spectrum of pump laser has been studied for a range of gas densities.Simulations reproduce both Stokes and anti-Stokes shifts in the spectrum of pump laser, similar to those observed in experiments in the spectrum of probe. A good agreement has been achieved with the experiments on the effect of variation in gas density on Stokes/anti-Stokes intensity. In addition, self-injection and trapping of electrons into the self-modulated wakes has been observed and analyzed. Here, the experimentally validated code SPACE will be used for predictive simulations to guide future experiments at BNL ATF.
Few of those who read the 1963 research note by Graeme A. Bird in Physics of Fluids could have imagined that 50 years later the proposed new numerical technique would have become the dominant numerical technique in molecular gas dynamics. The introduction of the Direct Simulation Monte Carlo (DSMC) method not only has altered the field of molecular gas dynamics but also has influenced fields such as physical chemistry, mathematics, computer science, and aerothermodynamics. Further, the DSMC method has been used to probe into previously uninvestigated theoretical aspects of the Boltzmann equation and has also served as a platform for the development of nonequilibrium chemistry models. The DSMC method’s most noteworthy achievement is that molecular gas dynamics became a practical tool in the hands of aerospace engineers in many situations of spacecraft design and mission analysis.
Sampling from probability distributions of the form 𝝈 ∝ e −𝜷V , where V is a continuous potential, is a fundamental task across physics, chemistry, biology, computer science, and statistics. However, when V is nonconvex, the resulting distribution becomes non-logconcave, and classical methods such as Langevin dynamics often exhibit poor performance. We introduce a quantum algorithm that provably accelerates a broad class of continuous-time sampling dynamics. For Langevin dynamics, our method encodes the target Gibbs measure into the amplitudes of aquantum state, identified as the kernel of a block matrix derived from a factorization of the Witten Laplacian operator. This connection enables Gibbs sampling via singular value thresholding and yields up to a quartic quantum speedup over best-knownclassical Langevin-based methods in the non-logconcave setting. Building on this framework, we further develop the first quantum algorithm that accelerates replica exchange Langevin diffusion, a widely used method for sampling from complex, rugged energy landscapes.
We revisit the exquisite archival radio data for the Type Ic supernova SN 1994I and present a revised model for the supernova (SN) radio emission and a pilot study that aims to constrain the rate of C-band radio transients within the face-on host galaxy, M51 (NGC 5194). We find that the temporal and spectral evolution of the SN 1994I radio emission are well fit by a synchrotron self-absorption model and use this to estimate physical parameters. We compute a pre-explosion mass loss rate of M-dot =3.0×10{sup −5} M{sub ⊙} yr{sup −1} for the progenitor, consistent with those observed from galactic Wolf–Rayet stars. Our model makes different assumptions for the dynamical model for the shockwave interaction than the model previously published by Weiler et al., but our M-dot is consistent with theirs to within errors and assumptions. Drawing from a subset of the archival radio observations from the Very Large Array collected for the monitoring of SN 1994I, we conduct a pilot study to search for previously unidentified transients. Data were primarily taken at a frequency of 4.9 GHz and are logarithmic in cadence, enabling sensitivity to transients with variability timescales ranging from days to months. We find no new transient detections in 31 epochs of data, allowing us to place a 2σ upper limit of 17 deg{sup −2} for the source density of radio transients above 0.5 mJy (L ≳ 4×10{sup 25} erg s{sup −1} Hz{sup −1} at the distance of M51). This study highlights the feasibility of utilizing archival high-cadence radio studies of SN host galaxies to place constraints on the radio transient rate as a function of luminosity in the local universe.
Increasing wind turbine size can have substantial economic benefits. However, increasing wind turbine turbine size requires the use of models that accurately capture the aerodynamic response to large deflections and nonstraight blade geometries. Additionally, the installation of offshore floating turbines requires turbine models that can handle large rotor motion. Free vortex wake (FVW) methods are capable of modeling this type of complex physics while remaining computationally tractable to perform the many simulations necessary for turbine design. In this work, a FVW model is added to the National Renewable Energy Laboratory engineering tool OpenFAST to allow for the aerodynamic modeling of highly flexible turbines along with the aero-hydro-servo-elastic response capabilities of OpenFAST. Code implementation is not complete and thus no results are presented. However, areas of potential modeling improvement are discussed, such as modifications to the circulation calculation and the need for different FVW input model specifications when high-thrust cases are being modeled.
Modeling physical phenomena through computational simulation increasingly relies on generating a collection of related runs, known as an ensemble. In this paper, we explore the challenges we face in developing analysis and visualization systems for large and complex ensemble data sets, which we seek to understand without having to view the results of every simulation run. Implementing approaches and ideas developed in response to this goal, we demonstrate the analysis of a 15K run material fracturing study using Slycat, our ensemble analysis system.
We formulate a particle and force level, activated dynamics-based statistical mechanical theory for the continuous startup nonlinear shear rheology of ultradense glass-forming hard sphere fluids and colloidal suspensions in the context of the elastically collective nonlinear Langevin equation approach and a generalized Maxwell model constitutive equation. Activated structural relaxation is described as a coupled local-nonlocal event involving caging and longer range collective elasticity which controls the characteristic stress relaxation time. Theoretical predictions for the deformation-induced enhancement of mobility, the onset of relaxation acceleration at remarkably low values of stress, strain, or shear rate, apparent power law thinning of the steady-state structural relaxation time and viscosity, a nonvanishing activation barrier in the shear thinning regime, an apparent Herschel–Buckley form of the shear rate dependence of the steady-state shear stress, exponential growth of different measures of a yield or flow stress with packing fraction, and reduced fragility and dynamic heterogeneity under deformation were previously shown to be in good agreement with experiments. The central new question we address here is the defining feature of the transient response—the stress overshoot. In contrast to the steady-state flow regime, understanding the transient response requires an explicit treatment of the coupled nonequilibrium evolution of structure, elastic modulus, and stress relaxation time. We formulate a new quantitative model for this aspect in a physically motivated and computationally tractable manner. Theoretical predictions for the stress overshoot are shown to be in good agreement with experimental observations in the metastable ultradense regime of hard sphere colloidal suspensions as a function of shear rate and packing fraction, and accounting for deformation-assisted activated motion appears to be crucial for both the transient and steady-state responses.
Complex physics and long computation time hinder the adoption of computer aided engineering models in the design of large-format battery cells and systems. A modular, efficient battery simulation model -- the multiscale multidomain (MSMD) model -- was previously introduced to aid the scale-up of Li-ion material and electrode designs to complete cell and pack designs, capturing electrochemical interplay with 3-D electronic current pathways and thermal response. Here, this paper enhances the computational efficiency of the MSMD model using a separation of time-scales principle to decompose model field variables. The decomposition provides a quasi-explicit linkage between the multiple length-scale domains and thus reduces time-consuming nested iteration when solving model equations across multiple domains. In addition to particle-, electrode- and cell-length scales treated in the previous work, the present formulation extends to bus bar- and multi-cell module-length scales. We provide example simulations for several variants of GH electrode-domain models.
First release of version 2 of the Energy Exascale Earth System Model. The atmosphere component remains EAM. Major changes since version 1 include: all column-physics parameterizations are computed on a separate grid that has approximately half the number of points of the dynamics grid, a new nonhydrostatic dynamical core (running in hydrostatic mode) with semi-Lagrangian tracer transport, CLUBB updated from v1 to v2, a new convective trigger (dCAPE/ULL) based on the dynamic Convective Available Potential Energy (CAPE) (dCAPE) and the Unrestricted Launch Level (ULL) concepts is used in ZM. minimum cloud droplet number changed, gravity wave drag energy conservation fixed and new tunings used, dust emission size distribution changed to emit more coarse dust particles The land component is still ELM. Major changes since version 1 include: using SNICAR-AD for radiation in snow to match the sea-ice model and fixing bugs in snow compaction and water state calculation. The ocean component remains MPAS-ocean. Major change since version 1 include: Redi isopycnal mixing has been updated, tested, and tuned in combination with the Gent-McWilliams parameterization, a sign error was fixed in the 3rd-order flux routines, the mesh used in low-resolution coupled cases was modified and the time steps adjusted, new regionally refined meshes were created, one focused on North America and another on the Southern Ocean, including ice shelf cavities. The sea-ice component remains MPAS-seaice. Major changes since version 1 include: A new heat- and freshwater-conserving coupling of frazil ice, turning off of SSH filtering, addition of SNICAR-AD and snow grain aging. The land-ice component remains MPAS-Albany-landIce (MALI) and is a static ice sheet. There is more out-of-the-box support for cryosphere configurations including the new regionally refined configuration around Antarctica with ice shelf cavities. The river model is MOSART. The half-degree river mesh was redone so it no longer treats Black and Caspian seas as ocean. The coupler remains cpl7/MCT. Major changes since version 1 include handling of ice shelf melt fluxes (heat / freshwater exchange with the ocean), and data icebergs. All components allow regional refinement of their meshes and two separate example of refinement, one in and around North America and one around Antarctica and the Southern Ocean, are provided.
This study presents a diagnosis of a multi-resolution approach using the Model for Prediction Across Scales - Atmosphere (MPAS-A) for simulating regional climate. Four AMIP experiments are conducted for 1999-2009. In the first two experiments, MPAS-A is configured using global quasi-uniform grids at 120 km and 30 km grid spacing. In the other two experiments, MPAS-A is configured using variable-resolution (VR) mesh with local refinement at 30 km over North America and South America embedded inside a quasi-uniform domain at 120 km elsewhere. Precipitation and related fields in the four simulations are examined to determine how well the VR simulations reproduce the features simulated by the globally high-resolution model in the refined domain. In previous analyses of idealized aqua-planet simulations, the characteristics of the global high-resolution simulation in moist processes only developed near the boundary of the refined region. In contrast, the AMIP simulations with VR grids are able to reproduce the high-resolution characteristics across the refined domain, particularly in South America. This indicates the importance of finely resolved lower-boundary forcing such as topography and surface heterogeneity for the regional climate, and demonstrates the ability of the MPAS-A VR to replicate the large-scale moisture transport as simulated in the quasi-uniform high-resolution model. Outside of the refined domain, some upscale effects are detected through large-scale circulation but the overall climatic signals are not significant at regional scales. Our results provide support for the multi-resolution approach as a computationally efficient and physically consistent method for modeling regional climate.
Verification of physics models and computer simulations are heavily reliant upon the accuracy of experimental measurements. Calibration of instrument responses becomes an important step to achieve this goal. This paper presents systematic studies of bent potassium acid phthalate (KAP) crystals using Lawrence Berkeley National Laboratories, Advanced Light Source, beamline 9.3.1 in the energy range of 2.3 to 7.5 keV. A set of KAP crystals, gradually bent from flat up to a 50.8 mm cylindrical curvature. The measured integrated reflectivity for this set of KAP crystals shows good agreement with the X-ray Oriented Program (XOP) calculations when adjusting the Debye–Waller temperature factor and using the multilamellar model in the calculations. Significant differences in rocking curve profiles were observed between experimental measurements and theory. A forward convolution model and software code were developed to include experimental parameters, allowing the investigation of the difference between measurements and calculations. After considering the experimental parameters, good agreements were obtained for the rocking curve profiles for all bending radii with a unique set of parameters. Our results show that XOP can be a useful and reliable tool to predict performance of cylindrically bent KAP crystals in this energy range.
Nuclear Safety is a review journal that covers significant developments in the field of nuclear safety. Its scope includes the analysis and control of hazards associated with nuclear energy, operations involving fissionable materials, and the products of nuclear fission and their effects on the environment. Primary emphasis is on safety in reactor design, construction, and operation; however, the safety aspects of the entire fuel cycle, including fuel fabrication, spent-fuel processing, nuclear waste disposal, handling of radioisotopes, and environmental effects of these operations, are also treated. Table of Contents for this issue follows. GENERAL SAFETY CONSIDERATIONS: 1 An Integrated Regional Approach To Risk Management of Industrial Systems, S. Chakraborty and R. Stratton; 9 Annual Technical Meeting of the NRC Cooperative Severe Accident Research Program, E. G. Silver; ACCIDENT ANALYSIS: 13 Deterministic Severe Accident Criteria as Severe Accident Design Criteria and Policy for the New Production Reactor-Heavy Water Reactor, P. T. Rhoads; 20 Proposed Deterministic Severe Accident Criteria for the Heavy Water Reactor-New Production Reactor Containment, K. D. Bergeron, S. E. Slezak, and C. E. Leach; 33 Quasi-Static Core Liquid Level Depression and Long-Term Core Uncovery During a PWR LOCA, Y. Kukita, R. R. Schultz, H. Nakamura, and J. Katayama; 49 LIRA: An Advanced Containment System to Minimize the Accidental Radioactivity Releases, A. Turricchia; 63 Errata to “A Review of Hydrogen Production During Melt/Water Interaction in LWRs,” Vol. 33, No. 4; CONTROL AND INSTRUMENTATION: 64 The Nuplex 80+™ Advanced Control Complex from ABB Combustion Engineering, F. Ridolfo, D. Harmon, and K. Scarola; DESIGN FEATURES: 76 A Scheme for Passive Isolation of the Containment of a Reactor, A. K. Ghosh, V. V. Raj, and K. Kakodkar; 84 Utility Requirements for Safety in the Passive Advanced Light-Water Reactor, T. U. Marston, W. H. Layman, and G. Bockhold, Jr.; WASTE AND SPENT FUEL MANAGEMENT: 91 Activities Related to Waste and Spent Fuel Management, Compiled by M. D. Muhlheim and E. G. Silver; OPERATING EXPERIENCES: 103 Testing Deficiencies in Auxiliary Systems Feedwater Systems, J. D. Kueck; 110 Reactor Shutdown Experience, Compiled by J. W. Cletcher; 113 Selected Safety-Related Events, Compiled by G. A. Murphy; 115 Operating U.S. Power Reactors, Compiled by M. D. Muhlheim and E. G. Silver; RECENT DEVELOPMENTS: 133 General Administrative Activities, Compiled by M. D. Muhlheim and E. G. Silver; 145 Reports, Standards, and Safety Guides, D. S. Queener; 150 Proposed Rule Changes as of Sept. 30, 1992; ANNOUNCEMENTS: 75 Harvard Short Course on Risk: Science, Assessment, and Management; 102 International Conference on Reactor Physics and Reactor Computations; 102 The Fourth International Conference on Nuclear Waste Reprocessing and Waste Management [RECOD ’94]; 132 International Symposium on the Radiological Exposure of the Population of the European Community from Radioactivity in the Mediterranean Sea [MARINA-MED]; 162 International High-Level Radioactive Waste Conference ‘Technology Integration—Above and Below”; 154 The Authors; 158 Indexes to Nuclear Safety, Volume 33.
In this final scientific/technical report covering the period of 3.5 years started on July 1, 2011, we report the accomplishments on the study of thermo-mechanical properties of Mo-based intermetallic compounds under NETL support. These include computational method development, physical properties investigation of Mo-based compounds and alloys. The main focus is on the mechanical and thermo mechanical properties at high temperature since these are the most crucial properties for their potential applications. In particular, recent development of applying ab initio molecular dynamic (AIMD) simulations to the T1 (Mo 5 Si 3 ) and T2 (Mo 5 SiB 2 ) phases are highlighted for alloy design in further improving their properties.
As part of technology foraging for the Responder Technology Alliance, established by the Department of Homeland Science and Technologies First Responders Group, this report summarizes technologies that are relevant in the area of “wearables,” with the potential for use by first responders. The content was collected over the previous month(s) and reproduced from a general Internet search using the term wearables. Additional information is available at the websites provided. This report is not meant to be an exhaustive list nor an endorsement of any technology described herein. Rather, it is meant to provide useful information about current developments in the areas wearable technology. The content is organized by technology function including: Sensors – Devices that detect physiological, particle, and chemical activity; Displays – Heads-up and body-worn visual displays; Power – Wearable power systems including chargers, batteries, self-powering or harvesting technologies, and power supplies; Integrated Communications – Voice and data communications systems utilizing Bluetooth, wireless, hands-free, ergonomically optimized systems, noise-filtering digital speakers or microphones, etc.; Exoskeletons – Whole or partial body suit that enhances mobility and physical performance; Wearable Computers – Body-worn data processing devices; General – Miscellaneous technologies as well as emerging trends or recent advances in the field of wearables.
The High-Energy Physics Forum for Computational Excellence (HEP-FCE) was formed by the Department of Energy as a follow-up to a recent report from the Topical Panel on Computing and the associated P5 recommendation. It is a pilot project distributed across the DOE Labs. During this initial incubation period the Forum is to develop a plan for a robust, long-term organization structure and a functioning web presence for forum activities and outreach, and a study of hardware and software needs across the HEP program. In the following sections we give this working group’s “vision” for aspects and qualities we wish to see in a future HEP-FCE. We then give a prioritized list of technical activities with suggested scoping and deliverables that can be expected to provide cross-experiment benefits. The remaining bulk of the report gives a technical survey of some specific “areas of opportunity” for cross-experiment benefit in the realm of software libs/tools. This survey serves as support for the vision and prioritized list. For each area we describe the ways that cross-experiment benefit is achieved today, as well as describe known failings or pitfalls where such benefit has failed to be achieved and which should be avoided in the future. For both cases, we try to give concrete examples. Each area then ends with an examination of what opportunities exist for improvements in that particular area.
Compute power is continually increasing, but this increased performance is largely found in sophisticated computing devices and supercomputer resources that are difficult to use, resulting in under-utilization. We developed a unified set of programming tools that will allow users to take full advantage of the new technology by allowing them to work at a level abstracted away from the platform specifics, encouraging the use of modern computing systems, including government-funded supercomputer facilities.