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 Xyce™ Parallel Electronic Simulator has been written to support the simulation needs of Sandia National Laboratories’ electrical designers. Xyce™ is a SPICE-compatible simulator with the ability to solve extremely large circuit problems on large-scale parallel computing platforms, but also includes support for most popular parallel and serial computers.
Wind loading and the fluctuating pressure loads it creates on PV panel surfaces are associated with multiple degradation mechanisms and failures. Modest wind speeds create reversing loads that can initiate cell cracks and weather cracked cells. Stronger wind speeds and extreme weather events can lead to larger scale forces and the aerodynamic instability known as torsional galloping. All these effects are dependent on the complex coupling between wind speed, panel orientation, and a myriad of other hardware and site-specific factors. In this work, we present the latest developments from our work to build an open-source, high-performance computing (HPC) fluid dynamics solver to predict and mitigate these effects. This simulation package allows users to easily specify different array layouts, solar-tracking angles, panel geometries, and weather conditions before automatically generating a refined computational mesh and solving for the unsteady loading on each panel surface. Small domains (e.g., a single panel row in isolation) can be solved on a modern laptop, while larger domains or very high-fidelity studies can be solved on distributed or HPC resources with minimal modifications to the underlying problem specification. We present preliminary case studies obtained using this simulation package and highlight how increased wind speeds combined with sub-optimal tracking angles can exacerbate degradation drivers.
Direct Numerical Simulations of realistic combustion devices are extremely challenging due to the wide separation of scales in the simulation, for example an internal combustion (IC) engine chamber, and the flame thickness of a high-pressure flame. The PeleLMeX solver uses adaptive mesh refinement (AMR) to evolve multi-species reacting flows in the low Mach number limit at the Exascale and relies on an embedded boundary (EB) approach to represent complex geometries. In that framework, the EB geometries often give rise to very small cut-cells along the boundary, which translate into extreme ill-conditioning of the pressure-projection, with eigenvalues that span 15-16 orders of magnitude. In this talk, we focus on the case of a typical IC piston bowl geometry for which we present on a novel approach towards solving these nearly singular linear systems with ILU-based, C-AMG smoothers on massively parallel architectures. In particular, we use scaling and equilibration algorithms to handle the non-normality of the upper triangular factors. This enables us to approximate the highly sequential triangular solve algorithm, embedded in the AMG smoothing-solve phase, with Jacobi iterations. This approximation can be written as a convergent Neumann series whose terms are composed of highly parallel sparse matrix vector multiplications. The result is an algorithm that substantially decreases setup and solve time, compared to state-of-the-art, for these challenging linear systems.
Focal Area(s): We will use unsupervised machine learning methods to identify and quantify the influence of large scale natural modes of climate variability to gain insight into the observed and simulated seasonal average and extreme precipitation changes. Science Challenge: A recent paper, led by co-PI Mark Risser, finds that although much of the variability in seasonal average and extreme precipitation over CONUS is unforced, the effect of large-scale modes of circulation variability (such as ENSO, AMO, PNA, etc.) can be detected and attributed. However, it is unclear whether or not unsupervised learning methods can (a) replicate this finding or (b) yield insight into possible nonlinear behavior that was not captured in the initial statistical analysis. Further work would entail extending this framework to other global land areas.
An electric power grid is one of the complex infrastructures, and because of its complex nature, there is simply no way that outages can be completely avoided. Thus, a modern society that depends on reliable electric supply requires a resilient electric system that can recover from disruptions while integrating emerging smart grid technologies. Here, this article presents a novel approach for service restoration in a modern power distribution system with controllable switches and distributed generation (DG) resources for any kind of outage. The proposed framework supports both the traditional service restoration using feeder reconfiguration and the grid-forming DG-assisted intentional islands that are dynamically sized using algorithms based on the fault scenario, available resources, and priority of loads. The problem is formulated as a mixed-integer linear program that incorporates critical system connectivity and operating constraints. Simulations are performed to demonstrate the effectiveness of the proposed approach using a large-scale four-feeder 1069-bus three-phase unbalanced distribution test system. It is demonstrated that the framework is effective in utilizing all available resources in quickly restoring the power supply to improve resiliency during extreme events and is scalable for a large-scale unbalanced power distribution system.
Abstract A ∼50 km resolution atmospheric general circulation model (GCM) is used to investigate the impact of radiative interactions on spatial organization of convection, the model's mean state, and extreme precipitation events in the presence of realistic boundary conditions. Mechanism‐denial experiments are performed in which synoptic‐scale feedbacks between radiation and dynamics are suppressed by overwriting the model‐generated atmospheric radiative cooling rates with its monthly varying climatological values. When synoptic‐scale radiative interactions are disabled, the annual mean circulation and precipitation remain almost unchanged, however tropical convection becomes less aggregated, with an increase in cloud fraction and relative humidity in the free troposphere but a decrease in both variables in the boundary layer. Changes in cloud fraction and relative humidity in the boundary layer exhibit more sensitivity to the presence of radiative interactions than variations in the degree of aggregation. The less aggregated state is associated with a decrease in the frequency of extreme precipitation events, coincident with a decrease in the dynamical contribution to the magnitude of extreme precipitation. At regional scales, the spatial contrast in radiative cooling between dry and moist regions diminishes when radiative interactions are suppressed, reducing the upgradient transport of energy, degree of aggregation, and frequency of extreme precipitation events. However, the mean width of the tropical rain belt remains almost unaffected when radiative interactions are disabled. These results offer insights into how radiation‐circulation coupling affects the spatial organization of convection, distributions of clouds and humidity, and weather extremes.
The Clausius–Clapeyron (CC) scaling, which indicates a roughly 7% increase in saturated water vapor per 1 °C increase in temperature, can serve as a strong constraint linking the intensity of precipitation extremes and local temperature. However, the relationship between precipitation extreme and local temperature (referred to as the PE-T relationship) does not always follow the CC scaling and is highly dependent on climate regimes. In this study, we investigated the impacts of different moisture processes and precipitation types on the PE-T relationship over eastern China during the summertime based on convection-permitting model simulations. Consistent with observations, the simulated intensity of precipitation extremes increases with temperature at a rate close to CC (double-CC) scaling below (above) 20 °C. When the temperature exceeds 25 °C, precipitation intensity starts to drop. Precipitation extremes are mainly contributed by the stratiform, MCS (i.e., mesoscale convective system) convective, and non-MCS convective precipitation at low (< 20 °C), medium (20–25 °C), and high (> 25 °C) temperatures, respectively, suggesting that the double-CC scaling occurs when convective types become dominant, while the negative scaling at high temperatures is attributed to the reduced horizontal scale of convection. Corresponding to the reduced intensity of precipitation at high temperatures, there are stronger divergence and subsidence in the low-level atmosphere, which is probably caused by the net cooling associated with the enhanced melting and evaporation of falling hydrometeors due to the lower relative humidity in the low-level atmosphere. Overall, our findings contribute to a deeper understanding of the temperature dependence of precipitation extremes in eastern China.
Large-scale meteorological pattern (LSMP)–based analysis is used in a novel way to understand meteorological conditions before and during short-duration dry spells over the northeastern United States. These LSMPs are useful to assess models and select better-performing models for future projections. Dry-spell events are identified from histograms of consecutive dry days below a daily precipitation threshold. Events lasting 12 days or longer, which correspond to ~10% of dry-spell events, are examined. The 500-hPa stream-function anomaly fields for the first 12 days of each event are time averaged, and k -means clustering is applied to isolate the dry-spell-related LSMPs. The first cluster has a strong low pressure anomaly over the Atlantic Ocean, southeast of the region, and is more common in winter and spring. The second cluster has strong high pressure over east-central North America and is most common during autumn. Over the region, both clusters have negative specific humidity anomalies, negative integrated vapor transport from the north, and subsidence associated with a midlatitude jet stream dipole structure that reinforces upper-level convergence. Subsidence is supported by cold-air advection in the first cluster and the location on the east side of the lower-level high pressure in the second cluster. Extratropical cyclone storm tracks are generally shifted southward of the region during the dry spells. Individual events lie on a continuum between two distinct clusters. These clusters have similar local, but different remote, properties. Although dry spells occur with greater frequency during drought months, most dry spells occur during nondrought months. Significance Statement: This study examines the large-scale weather patterns and meteorological conditions associated with dry-spell events lasting at least 2 weeks while affecting the northeastern United States. A statistical approach groups events together on the basis of similar atmospheric features. We find two distinct sets of patterns that we call large-scale meteorological patterns. These patterns reduce moisture, foster localized sinking, and shift the storm track southward along the Atlantic seaboard, all of which reduce precipitation. Besides greater understanding, knowing the meteorological patterns during short-term dryness in the region provides an important tool to assess how well atmospheric models reproduce these specific patterns. More dry spells occur in non-drought months than in drought months, which means that dry spells can occur without preexisting drought conditions.
Coupling multiphase flow with energy transport due to high temperature heat sources introduces significant new challenges since boiling and condensation processes can lead to dry-out conditions with subsequent re-wetting. The transition between two-phase and single-phase behavior can require changes to the primary dependent variables adding discontinuities as well as extending constitutive nonlinear relations to extreme physical conditions. Practical simulations of large-scale engineered domains lead to Jacobian systems with a very large number of unknowns that must be solved efficiently using iterative methods in parallel on high-performance computers. Performance assessment of potential nuclear repositories, carbon sequestration sites and geothermal reservoirs can require numerous Monte-Carlo simulations to explore uncertainty in material properties, boundary conditions, and failure scenarios. Due to the numerical challenges, standard NR iteration may not converge over the range of required simulations and require more sophisticated optimization method like trust-region. In this study, we use the open-source simulator PFLOTRAN for the important practical problem of the safety assessment of future nuclear waste repositories in the U.S. DOE geologic disposal safety assessment Framework. The simulator applies the PETSc parallel framework and a backward Euler, finite volume discretization. We demonstrate failure of the conventional NR method and the success of trust-region modifications to Newton’s method for a series of test problems of increasing complexity. Trust-region methods essentially modify the Newton step size and direction under some circumstances where the standard NR iteration can cause the solution to diverge or oscillate. Furthermore, we show how the Newton Trust-Region method can be adapted for Primary Variable Switching (PVS) when the multiphase state changes due to boiling or condensation. The simulations with high-temperature heat sources which led to extreme nonlinear processes with many state changes in the domain did not converge with NR, but they do complete successfully with the trust-region methods modified for PVS. This implementation effectively decreased weeks of simulation time needing manual adjustments to complete a simulation down to a day. Finally, we show the strong scalability of the methods on a single node and multiple nodes in an HPC cluster.
The direct use of solar thermal energy provides opportunities for low-cost heating sources for a variety of applications. Ultra-high temperatures around 1000°C are high value and highly useful for energy-demanding industries. Many materials cannot withstand these conditions. In the area of Sustainable Chemicals, further limitations on material stability exist. Combining state-of-the-art materials with new designs provides a promising pathway for harvesting solar thermal energy and performing high temperature chemical processes. However, conventional manufacturing limits the potential for design flexibility. In this project, Additive Manufacturing was combined with advanced materials and new chemical reactor designs. In addition, 24/7 energy is necessary for chemical processing, and designs for ultra-high temperature thermal storage were devised. Specifically, preliminary design of a novel solar thermal receiver was developed in this project and designed to work with thermocatalytic reactors for producing sustainable chemicals and fuels. An ultra-high temperature particle storage system and heat exchangers were proposed to transport ultra-hot air as thermal fluid for the system. On a broader scale, this system could be used to tap solar thermal energy for a centralized facility with capability of transferring that heat to various segments at a full range of temperatures to 1000°C. The project pushed the temperature boundaries past those in current use, and Additive Manufacturing was envisaged for fabricating the receiver to meet requirements of extreme environments. An extensive analysis of silicon carbide additive manufacturing was performed to compare the thermal and mechanical properties of complex geometries compared to conventional material and those manufactured via other methods. The Additive Manufacturing via Binder-Jet printing was optimized and characterized to provide high quality and reproducible components capable of withstanding the proposed extreme environments. The designs for the concentrating solar thermal cavity with ultra-hot air thermal fluid showed high performance in simulations, attributable to the complex optimized geometries of the 3D printed systems. The bright future of Additive Manufacturing with advanced materials developments should provide more options and even higher quality as the technology further develops. Current costs for Additive Manufacturing of advanced ceramics is relatively low, however post-processing of the materials for extreme environments is currently high. There is little industrial-scale infrastructure for these, but it is growing as niche applications become more mainstream. The results of the project can be translated into similar extreme environments for concentrating solar thermal energy as well as its integration with ultra-hot air thermal fluids. A number of industries that require ultra-high temperatures need to electrify or otherwise decarbonize for climate goals, and this project showed that theoretically there is a pathway to do so with direct concentrated solar thermal power.
Atmospheric rivers (ARs), filaments of intense atmospheric moisture transport, play a significant role in delivering moisture poleward into the Arctic and triggering weather extremes. Although previous studies have focused on large‐scale circulations driving these events, this study investigates ARs through attributing their moisture sources using the Community Atmosphere Model version 5 (CAM5) with moisture‐tagging capability. Examining ARs in the Atlantic and Pacific sectors of the Arctic separately revealed distinct contributions from remote versus local, and ocean versus land, moisture sources. Unlike non‐AR events, Arctic ARs primarily draw moisture from their respective ocean basins in lower‐latitude regions during the cold season, and shift to land sources in the warm season. Cold‐season ARs in the Atlantic and Pacific sectors source 73.2% and 85.3% of their moisture from their respective ocean basins at lower latitudes, however warm‐season contributions decrease to 41.3% and 29.4%. In contrast, mid‐ and high‐latitude continents combined contribute 40.1% and 51.3% in the warm season. Trajectory clustering analysis shows that AR moisture sources depend on both their genesis regions and transport pathways. In recent decades, the Arctic has experienced a moistening trend in both the cold and warm seasons. Our results suggest that local and remote sources equally drive the observed cold‐season moistening, whereas remote sources predominantly drive the warm‐season moistening. A better understanding of Arctic AR moisture sources under current climate conditions provides valuable insights into their potential future changes amid the projected heterogeneous northern hemispheric warming.
Abstract There is a need for high-throughput, scale-relevant, and direct electrochemical analysis to understand the corrosion behavior and sensitivity of nuclear materials that are exposed to extreme (high pressure, temperature, and radiation exposure) environments. We demonstrate the multi-scale, multi-modal application of scanning electrochemical cell microscopy (SECCM) to electrochemically profile corrosion alterations in nuclear alloys in a microstructurally resolved manner. Particularly, we identify that both mechanically deformed and irradiated microstructures show reduced charge-transfer resistance that leads to accelerated oxidation. We highlight that the effects of mechanical deformation and irradiation are synergistic, and may in fact, superimpose each other, with implications including general-, galvanic-, and/or irradiation-activated stress-corrosion cracking. Taken together, we highlight the ability of non-destructive, electrochemical interrogations to ascertain how microstructural alterations result in changes in the corrosion tendency of a nuclear alloy: knowledge which has implications to rank, qualify and examine alloys for use in nuclear construction applications.
Photovoltaics (PV), or solar electricity generation, has become the cheapest form of energy in many locations worldwide and, combined with energy storage, has the potential to satisfy most of our electricity needs. PV has grown at an annual compounded growth rate of approximately 30% in the last three decades. Solar energy systems will continue their impressive growth in distributed energy, microgrids, and utility scale, as efforts are made for dependable electricity in an age with increasing extreme weather. However, within this remarkable success lies a new challenge. The growth curve, combined with rapid product innovation and scale up, means that the majority of PV systems are new, without the three years of performance data that have been required in the past to estimate product lifetime. PV reliability has to address this challenge. In this review we present a brief summary of PV reliability starting with brief historical synopsis, detailing some of the technological challenges and present a framework required for long lifetime.
The microstructural length scale of metals changes by orders of magnitude under extreme processing conditions producing a concurrent wide range of mechanical strength and plasticity behaviors. A unified stress-strain σε σ(ε) model is formulated that’s based on superposing the components of asymptotic-curvilinear work hardening Θσ Θ(σ) to qualify and quantify these mechanical behaviors. This approach accounts for the rapid strengthening of metals beyond the initial yield point, through extended steady-state deformation, to the structural instability. The relationship between the softening coefficients cbi c bi of the work hardening formulation Θσ Θ(σ) and strength are found to reveal the microstructural scale in the material. Specifically, the rapid decrease in the slope of the Θσ Θ(σ) curve provides a measure for microstructural size consistent with a functional Hall-Petch relationship of strength. A successful application is shown for the tensile behavior of pure aluminum microstructures that result from extreme plastic deformation by equal-channel angle pressing.
At continuous phase transitions, quantum many-body systems exhibit complex, emergent behavior. Most notably, at a quantum critical point, correlations decay as a power law, with exponents determined by a set of universal scaling dimensions. Experimentally probing such power law correlations is extremely challenging, owing to the interplay between decoherence, the vanishing energy gap, and boundary effects. In this work, we used a Rydberg quantum simulator to adiabatically prepare critical ground states of both a one-dimensional ring and a two-dimensional square lattice. By accounting for and tuning the openness of our quantum system, which is well-captured by a single phenomenological length scale, we directly observed power law correlations and extracted the corresponding scaling dimensions. Our work complements recent studies of quantum criticality that use the Kibble-Zurek mechanism and digital quantum circuits.
The interface thickness in the phase-field (PF) method limits its simulation scales. Consequently, large-scale PF simulations become prohibitively expensive for resolving the extremely fine microstructures that typically form during rapid solidification processing. This challenge is significant in predicting microstructure evolution in metal additive manufacturing and has been identified by the United States Department of Energy’s Exascale Computing Project. Here, to address this, we develop a multi-GPU and MPI-based massively parallel simulation code, utilizing state-of-the-art algorithms, software, and libraries, for large-scale three-dimensional (3D) PF simulations. We report the first GPU-parallel PF simulations on Frontier (currently the second TOP500 exascale cluster) and Summit machines, taking dendritic growth as an example problem. We evaluate the parallel performance of our implementation using scaling studies with more than 24 000 GPUs (among the largest known computations to date) and the acceleration performance using large-scale simulations of dendritic growth in 3D. Finally, massively parallel GPUs in these supercomputers enabled the first coupled multiscale simulations of laser melting and subsequent dendritic solidification on the scale of a full melt-pool, demonstrating the feasibility of performing PF simulations with a point total over 2 billion grid points within an acceptable time.
ABSTRACT V2487 Oph is a recurrent nova with detected eruptions in 1900 and 1998. Startlingly, V2487 Oph shows flares, called ‘Superflares’, with up to 1.10 mag amplitude, fast rises of under one-minute, always with an initial impulsive spike followed by a roughly exponential tail, typically one-hour durations, and with random event times averaging once-per-day. The typical flare energy E is over 1038 erg, while the yearly energy budget is 1041 erg. V2487 Oph Superflares obey three relations; the number distribution of flare energies scales as E−2.34 ± 0.35, the waiting time from one flare to the next is proportional to E of the first event, and flare durations scale as E0.44 ± 0.03. Scenarios involving gravitational energy and nuclear energy fail to satisfy the three relations. The magnetic energy scenario, however, can explain all three relations. This scenario has magnetic field lines above the disc being twisted and amplified by the motions of their footprints, with magnetic reconnection releasing energy that comes out as Superflare light. This exact mechanism is already well known to occur in white light solar flares, in ordinary M-type flare stars, and in the many Superflare stars observed all across the H-R diagram. Superflares on Superflare stars have rise times, light-curve shapes, and durations that are very similar to those on V2487 Oph. So we conclude that the V2487 Oph Superflares are caused by large-scale magnetic reconnection. V2487 Oph is now the most extreme Superflare star, exhibiting the largest known flare energy (1.6 × 1039 erg) and the fastest occurrence rate.
A search is reported for excited $τ$-leptons and leptoquarks in events with two hadronically decaying $τ$-leptons and two or more jets. The search uses proton-proton ($pp$) collision data at $\sqrt{s}$ = 13 TeV recorded by the ATLAS experiment during the Run 2 of the Large Hadron Collider in 2015–2018. The total integrated luminosity is 139 fb -1 . The excited $τ$-lepton is assumed to be produced and to decay via a four-fermion contact interaction into an ordinary $τ$-lepton and a quark-antiquark pair. The leptoquarks are assumed to be produced in pairs via the strong interaction, and each leptoquark is assumed to couple to a charm or lighter quark and a $τ$-lepton. No excess over the background prediction is observed. Excited $τ$-leptons with masses below 2.8 TeV are excluded at 95% CL in scenarios with the contact interaction scale Λ set to 10 TeV. At the extreme limit of model validity where Λ is set equal to the excited $τ$-lepton mass, excited $τ$-leptons with masses below 4.6 TeV are excluded. Leptoquarks with masses below 1.3 TeV are excluded at 95% CL if their branching ratio to a charm quark and a $τ$-lepton equals 1. The analysis does not exploit flavour-tagging in the signal region.