Programming and Performance Visualization Tools
Explore the source record for details and available documents.
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.
Explore the source record for details and available documents.
Smart cities depend on flexible and secure energy systems to ensure resilient power for critical infrastructure; however, recent weather-related events and cyberattacks have highlighted weaknesses in our energy systems, with the potential for widespread economic and security impacts. As stated by the Executive Office of the President, "the resilience of the US electric grid is a key part of the nation's defense against severe weather." To address the energy delivery security challenge, microgrids are rising as a viable solution that enhances the flexibility and resilience of the distribution grid and boosts the reliability of the local supply for the end-user. Traditionally, high capital investment has been a barrier to large-scale adoption of microgrid technology. Understanding the flexibility and resilience benefits of microgrids and accounting for the associated value streams can make the microgrid's proposition economically viable. In this chapter, microgrids' utility and their potential to serve as a flexible and resilient resource for the utility grid by providing capabilities such as peak shaving, demand response, and frequency regulation is presented. Moreover, other value streams, such as (1) their ability to island during a disaster and sustain critical loads which makes them a robust resilience solution for end-users, in the event of the utility grid outage and (2) microgrids also provide a flexible platform for integrating distributed energy resources in conjunction with storage and conventional generation technologies, strengthen microgrid's role in reducing the over-arching goal of emission reduction. Given the myriad of benefits associated with microgrids, we present strategies which can be employed for making microgrid itself resilient against physical and cyberthreats by employing hardware, software, and personnel training solutions to operate the microgrid before, during, and after a potential disaster. This chapter, thus, provides a holistic study of the microgrid as a resilience resource, for the utility grid, and a self-contained end-user for the end-user.
The use of fossil fuels is the primary source of greenhouse gas emissions but there are alternatives to these especially in the form of biofuels, fuels derived from bioenergy crops. This paper aims to determine farmers’ potential adoption rates of newly introduced bioenergy crops with a specific example of carinata in the state of Georgia. The determination is done using an agent-based modeling technique with two principal assumptions—farmers are profit maximizer and they are influenced by neighboring farmers. Two diffusion parameters (traditional and expansion) are followed along with two willingness (high and low) scenarios to switch at varying production economics to carinata and other prominent traditional field crops (cotton, peanuts, corn) in the study region. We find that a contract prices around $9, $8 and $7 can be a viable option for encouraging farmers to adopt carinata in low, average, and high profit conditions, respectively. Expansion diffusion (that diffuses all over the geographical area), rather than centered to the few places like traditional diffusion at the early stage of adoption in conjunction with higher willingness conditions influences higher adoption rates in the short-term. As such, the model can be used to understand the behavioral economics of carinata in Georgia and beyond, as well as offering a potential tool to study similar bioenergy crops.
This book examines what role fuels have in the transformation of end-use in buildings and industrial processes. Energy efficient technologies are necessary to lower the carbon footprint for a transition towards clean energy in a sustainable manner. Efficient utilization of primary energy resources including renewables to support the current and future energy needs while targeting grid resiliency, energy and environmental security, at an affordable cost is of significant value. Analysis of configurations consisting of heat pumps, fuel driven thermal providers and power systems is presented. Sensitivity of electrical grid's carbon intensity towards carbon footprint in comparison with fuel driven technologies is necessary to recognize the true value proposition of currently available energy solutions for different end consumers. Similarly, the role of low carbon, zero carbon, and carbon negative fuels such as power to gas, power to liquid, hydrogen, biogas, etc. in conjunction with polygeneration technologies are discussed. Transformation of the primary energy resources from conventional fossil fuels to renewable fuels or electricity will have a significant impact on the overall carbon footprint of various end use sectors, including buildings. Hence, this book also examines two different scenarios focused on sensitivity of the pace of decarbonization of electrical grid and fuel supply on operational energy related carbon emissions.
In this work, a woven carbon/epoxy composite was subjected to fatigue crack growth under mixed Mode-I/Mode-II loading to obtain crack growth behavior at different cyclic strain energies. Owing to the woven structure of the material, pure Mode-II fracture is usually a difficult proposition because of friction, interference, and interlock of woven tows in adjacent plies at an interlaminar crack. These limitations were overcome by the use of a novel form of mixed Mode-I/Mode-II specimen, which imposes sufficient crack surface opening (Mode-I) to alleviate ply-ply interactions, but not so much as to obscure the sliding (Mode-II) response. Finally, comparison with pure Mode-I fatigue crack growth data, in conjunction with a fracture interaction criterion, provided a means to extract the Mode-II behavior.
In this paper, we investigate the OpenMP parallelization and optimization of two novel data classification algorithms. The new algorithms are based on graph and PDE solution techniques and provide significant accuracy and performance advantages over traditional data classification algorithms in serial mode. The methods leverage the Nystrom extension to calculate eigenvalue/eigenvectors of the graph Laplacian and this is a self-contained module that can be used in conjunction with other graph-Laplacian based methods such as spectral clustering. We use performance tools to collect the hotspots and memory access of the serial codes and use OpenMP as the parallelization language to parallelize the most time-consuming parts. Where possible, we also use library routines. We then optimize the OpenMP implementations and detail the performance on traditional supercomputer nodes (in our case a Cray XC30), and test the optimization steps on emerging testbed systems based on Intel’s Knights Corner and Landing processors. We show both performance improvement and strong scaling behavior. Finally, a large number of optimization techniques and analyses are necessary before the algorithm reaches almost ideal scaling.
Cast duplex stainless steel piping in light water nuclear reactors expe- rience thermal aging embrittlement during operational service. Interest in extending the operational life to 80 years requires an increased understanding of the microstructural evolution and corresponding changes in mechanical behavior. We analyze the evolution of the microstructure during thermal aging of cast CF-3 and CF-8 stainless steels using electron microscopy and atom probe tomography. The evolution of the mechanical properties is measured concurrently by mechanical methods such as tensile tests, Charpy V-notch tests, and instrumented nanoinden- tation. A microstructure-based finite element method model is developed and uti- lized in conjunction with the characterization results in order to correlate the local stress-strain effects in the microstructure with the bulk measurements. This work is supported by the DOE Nuclear Energy University Programs (NEUP), contract number DE-NE0000724.
In an automotive power electronics package, bonded interface materials play a critical role in providing a smooth heat dissipation path from the device to the coolant. However, these materials are susceptible to failure when the package gets exposed to repeated changes in its thermal environment and internal temperature variations. Under these thermal conditions, the coefficient of thermal expansion mismatch between the different component layers results in thermally induced stresses and consequently the delamination of the material from its adjoining surfaces or crack initiation and propagation within the material. Thermomechanical modeling offers a path to study the complex nature of these bonded materials, guides their design and selection, and, in conjunction with experimental results, plays an important role in predicting their lifetime. In this chapter, two modeling strategies - similar in approach and implementation but different in the underlying theory - are discussed in the context of high-temperature bonded materials. In addition, a short review of sintered silver is provided, as it is widely seen as a potentially promising and reliable bonded material for high-temperature applications.
HPC systems now exploit GPUs within their compute nodes to accelerate program performance. As a result, high-end application development has become extremely complex at the node level. In addition to restructuring the node code to exploit the cores and specialized devices, the programmer may need to choose a programming model such as OpenMP or CPU threads in conjunction with an accelerator programming model to share and manage the different node resources. This comes at a time when programmer productivity and the ability to produce portable code has been recognized as a major concern. In order to offset the high development cost of creating CUDA or OpenCL kernels, directives have been proposed for programming accelerator devices, but their implications are not well known. In this paper, we evaluate the state of the art accelerator directives to program several applications kernels, explore transformations to achieve good performance, and examine the expressivity and performance penalty of using high-level directives versus CUDA. We also compare our results to OpenMP implementations to understand the benefits of running the kernels in the accelerator versus CPU cores.
Although the fragment molecular orbital (FMO) method enables electronic structure calculations with near-linear scaling behavior with respect to system size, the computational cost of ab initio methods typically employed in conjunction with FMO is still prohibitive for routine calculations of very large systems or long timescale molecular dynamics simulations. We, therefore, combined the FMO and density-functional tight-binding (DFTB) method, which is one of the emerging semi-empirical quantum chemical methods, and have demonstrated that FMO-DFTB is capable of performing geometry optimizations for systems containing up to one million atoms using limited computational resources. In this chapter, we will review the basics of the DFTB method first before introducing FMO-DFTB, focusing on the relationship with density functional theory and other FMO methodologies. We also demonstrate the latest scalings of FMO-DFTB2 and DFTB3 using three-dimensional water clusters, showing that the most favorable scaling is O(N1.16). Applications of FMO-DFTB to various systems are briefly summarized, and an outlook to future applications is provided.
In ordinary gravitational theories, any local bulk operator in an entanglement wedge is accompanied by a long-range gravitational dressing that extends to the asymptotic part of the wedge. Islands are the only known examples of entanglement wedges that are disconnected from the asymptotic region of spacetime. In this paper, we show that the lack of an asymptotic region in islands creates a potential puzzle that involves the gravitational Gauss law, independently of whether or not there is a non-gravitational bath. In a theory with long-range gravity, the energy of an excitation localized to the island can be detected from outside the island, in contradiction with the principle that operators in an entanglement wedge should commute with operators from its complement. In several known examples, we show that this tension is resolved because islands appear in conjunction with a massive graviton. We also derive some additional consistency conditions that must be obeyed by islands in decoupled systems. Our arguments suggest that islands might not constitute consistent entanglement wedges in standard theories of massless gravity where the Gauss law applies.
Observables which distinguish boosted topologies from QCD jets are playing an increasingly important role at the Large Hadron Collider (LHC). These observables are often used in conjunction with jet grooming algorithms, which reduce contamination from both theoretical and experimental sources. In this paper we derive factorization formulae for groomed multi-prong substructure observables, focusing in particular on the groomed D 2 observable, which is used to identify boosted hadronic decays of electroweak bosons at the LHC. Our factorization formulae allow systematically improvable calculations of the perturbative D 2 distribution and the resummation of logarithmically enhanced terms in all regions of phase space using renormalization group evolution. They include a novel factorization for the production of a soft subjet in the presence of a grooming algorithm, in which clustering effects enter directly into the hard matching. We use these factorization formulae to draw robust conclusions of experimental relevance regarding the universality of the D 2 distribution in both e+e- and pp collisions. In particular, we show that the only process dependence is carried by the relative quark vs. gluon jet fraction in the sample, no non-global logarithms from event-wide correlations are present in the distribution, hadronization corrections are controlled by the perturbative mass of the jet, and all global color correlations are completely removed by grooming, making groomed D 2 a theoretically clean QCD observable even in the LHC environment. Lastly, we compute all ingredients to one-loop accuracy, and present numerical results at next-to-leading logarithmic accuracy for e + e - collisions, comparing with parton shower Monte Carlo simulations. Results for pp collisions, as relevant for phenomenology at the LHC, are presented in a companion paper.
A search for compositeness of electrons and muons is presented using a data sample of proton-proton collisions at a center-of-mass energy of $ \sqrt{s} = 8 $ TeV collected with the CMS detector at the LHC and corresponding to an integrated luminosity of 19.7 fb$^{−1}$. Excited leptons (ℓ$^{*}$) produced via contact interactions in conjunction with a standard model lepton are considered, and a search is made for their gauge decay modes. The decays considered are ℓ$^{*}$ → ℓγ and ℓ$^{*}$ → ℓZ, which give final states of two leptons and a photon or, depending on the Z-boson decay mode, four leptons or two leptons and two jets. The number of events observed in data is consistent with the standard model prediction. Exclusion limits are set on the excited lepton mass, and the compositeness scale Λ. For the case M$_{ℓ*}$ = Λ the existence of excited electrons (muons) is excluded up to masses of 2.45 (2.47) TeV at 95% confidence level. Neutral current decays of excited leptons are considered for the first time, and limits are extended to include the possibility that the weight factors f and f′, which determine the couplings between standard model leptons and excited leptons via gauge mediated interactions, have opposite sign.
Minimal gauged U 1 L μ - L τ models can provide for an additional source for the muon anomalous magnetic moment however it is difficult to accommodate the discrepancy in the electron magnetic moment in tandem. Owing to the relative sign between the discrepancies in these quantities, it seems unlikely that they arise from the same source. We show that a supersymmetric (SUSY) gauged U 1 L μ - L τ model can accommodate both the muon and electron anomalous magnetic moments in a very simple and intuitive scenario, without utilizing lepton flavor violation. The currently allowed parameter space in this kind of a scenario is constrained from the latest LHC and various low energy experimental data,e.g., recent COHERENT data, CCFR, Borexino, BaBaR, supernova etc. These constraints, in conjunction with the requirement to explain both lepton magnetic moments, lead to an upper bound on the first generation slepton mass, a lower bound on the second generation slepton mass and constricts the allowed range for the new gauge boson mass and coupling. The scheme can be probed at the ongoing COHERENT and Coherent CAPTAIN-Mills experiments and at future experiments, e.g., DUNE, BELLE-II etc.
The computation of the parton distribution functions (PDF) or distribution amplitudes (DA) of hadrons from first principles lattice QCD constitutes a central open problem in high energy nuclear physics. In this study, we present and evaluate the efficiency of several numerical methods, well established in the study of inverse problems, to reconstruct the full x-dependence of PDFs. Our starting point are the so called Ioffe time PDFs, which are accessible from Euclidean time simulations in conjunction with a matching procedure. Using realistic mock data tests, we find that the ill-posed incomplete Fourier transform underlying the reconstruction requires careful regularization, for which both the Bayesian approach as well as neural networks are efficient and flexible choices.
We study the discrimination of quark-initiated jets from gluon-initiated jets in monojet searches for dark matter using the technique of averaged jet energy profiles. We demonstrate our results in the context of effective field theories of dark matter interactions with quarks and gluons, but our methods apply more generally to a wide class of models. Different effective theories of dark matter and the standard model backgrounds each have a characteristic quark/gluon fraction for the leading jet. When used in conjunction with the traditional cut-and-count monojet search, the jet energy profile can be used to set stronger bounds on contact interactions of dark matter. In the event of a discovery of a monojet excess at the 14 TeV LHC, contact interactions between dark matter with quarks or with gluons can be differentiated at the 95% confidence level. For a given rate at the LHC, signal predictions at direct detection experiments for different dark matter interactions can span five orders of magnitude. Lastly, the ability to identify these interactions allows us to make a tighter connection between LHC searches and direct detection experiments.
A natural extension of the Pasterski-Shao-Strominger (PSS) prescription is described, enabling the map of Minkowski space amplitudes with massive spinning external legs to the celestial sphere to be performed. An integral representation for the conformal primary wave function (CPW) of massive spinning bosons on the celestial sphere is derived explicitly for spin-one and -two. By analogy with the spin-zero case, the spinning bulk-to-boundary propagator on Euclidean AdS is employed to extend the massive CPW integral representation to arbitrary integer spin, and to describe the appropriate inverse transform of massive spinning CPWs back to the plane wave basis in Minkowski space. Subsequently, a massive spin-s momentum operator representation on the celestial sphere is determined, and used in conjunction with known Lorentz generators to derive Poincaré symmetry constraints on generic massive spinning two-, three- and four-point celestial amplitude structures. Finally, as a consistency check, three-point Minkowski space amplitudes of two massless scalars and a spin-one or -two massive boson are explicitly mapped to the celestial sphere, and the resulting three-point function coefficients are confirmed to be in exact agreement with the results obtained from Poincaré symmetry constraints.
Results are reported from a search for supersymmetric particles in proton-proton collisions in the final state with a single, high transverse momentum lepton, multiple jets, including at least one b-tagged jet, and large missing transverse momentum. The data sample corresponds to an integrated luminosity of 2.3 fb$^{−1}$ at $ \sqrt{s}=13 $ TeV, recorded by the CMS experiment at the LHC. The search focuses on processes leading to high jet multiplicities, such as gluino pair production with $ \tilde{\mathrm{g}}\to \mathrm{t}\overline{\mathrm{t}}{\tilde{\chi}}_1^0 $ . The quantity M$_{J}$ , defined as the sum of the masses of the large-radius jets in the event, is used in conjunction with other kinematic variables to provide discrimination between signal and background and as a key part of the background estimation method. The observed event yields in the signal regions in data are consistent with those expected for standard model backgrounds, estimated from control regions in data. Exclusion limits are obtained for a simplified model corresponding to gluino pair production with three-body decays into top quarks and neutralinos. Gluinos with a mass below 1600 GeV are excluded at a 95% confidence level for scenarios with low $ {\tilde{\chi}}_1^0 $ mass, and neutralinos with a mass below 800 GeV are excluded for a gluino mass of about 1300 GeV. For models with two-body gluino decays producing on-shell top squarks, the excluded region is only weakly sensitive to the top squark mass.