Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “scale coupling”

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.

At least 523 records · Page 29

Global linear and nonlinear gyrokinetic simulations of tearing modes

Abstract To better understand multi-scale interactions between global tearing modes and microturbulence in the Madison Symmetric Torus (MST) reversed-field pinch (RFP), the global gyrokinetic code Gene is modified to describe global tearing mode instability via a shifted Maxwellian (SM) distribution consistent with experimental equilibria. The implementation of the SM is tested and benchmarked by comparisons with different codes and models. Good agreement is obtained in code-code and code-theory comparisons. Linear stability of tearing modes of a non-reversed MST discharge is studied. A collisionality scan is performed to the lowest order unstable modes ( n = 5, n = 6) and shown to behave consistently with theoretical scaling. The nonlinear evolution is simulated, and saturation is found to arise from mode coupling and transfer of energy from the most unstable tearing mode to small-scale stable modes mediated by the m = 2 tearing mode. The work described herein lays the foundation for nonlinear simulation and analysis of the interaction of tearing modes and gyroradius-scale instabilities in RFP plasmas.

Physics↗

Historical evaluation and future projections of compound heatwave and drought extremes over the conterminous United States in CMIP6 *

Abstract Independently, both droughts and heatwaves can induce severe impacts on human and natural systems. However, when these two climate extremes occur concurrently in a given region, their compound impacts are often more pronounced. With the improvement in both the spatiotemporal resolution and representation of complex climate processes in the global climate models (GCMs), they are increasingly used to study future changes in these extremes and associated regional impacts. However, GCM selection for such impact assessments is generally based on historical performance and/or future mean changes, without considering individual or compound extremes. In contrast, this study evaluates historical performance and projected changes in heatwaves, droughts, and compound heatwave-droughts using an ensemble of GCMs from the latest Phase 6 of Coupled Models Intercomparison Project at a regional scale across the conterminous United States. Additionally, we explore the inter-model differences in the projected changes that are associated with various characteristics of extremes and the choice of drought indices. Our analysis reveals considerable variation among the GCMs, as well as substantial differences in the projected changes based on the choice of drought indices and region of interest. For example, the projected increases in both the frequency and intensity of drought and associated compound extreme days, based on the standardized precipitation evapotranspiration index far exceed those derived from the standard precipitation index. Further, the largest changes in the frequency of compound extremes are projected over the Southwest, South Central, and parts of the Southeast while the smallest changes are projected over the Northeast. Overall, this study provides important insights for the interpretation and selection of GCMs for future assessment studies that are crucial for the development of regional adaptation strategies in the face of climate change.

54 ENVIRONMENTAL SCIENCES↗

Faster Randomized Dynamical Decoupling

We present a randomized dynamical decoupling (DD) protocol that can substantially improve the performance of any given deterministic DD scheme for suppressing coherent noise by using no more than two additional pulses. Our construction is implemented by probabilistically applying sequences of pulses, which, when combined, effectively eliminate the error terms that scale linearly with the system-environment coupling strength. As a result, we show that a randomized protocol using a few pulses can outperform deterministic DD protocols that require considerably more pulses. Furthermore, we prove that the randomized protocol provides an improvement compared to deterministic DD sequences that aim to reduce the error in the system’s Hilbert space, such as Uhrig DD, which had been previously regarded to be optimal. To rigorously evaluate the performance, we introduce new analytical methods suitable for analyzing higher-order DD protocols that might be of independent interest. Here, we also present numerical simulations confirming the significant advantage of using randomized protocols compared to widely used deterministic protocols.

Quantum algorithms & computation↗

Mapping and probing Froggatt-Nielsen solutions to the quark flavor puzzle

The Froggatt-Nielsen (FN) mechanism is an elegant solution to the flavor problem. In its minimal application to the quark sector, the different quark types and generations have different charges under a 𝑈⁢(1)𝑋 flavor symmetry. The SM Yukawa couplings are generated below the flavor breaking scale with hierarchies dictated by the quark charge assignments. Only a handful of charge assignments are generally considered in the literature. We analyze the complete space of possible charge assignments with |𝑋 𝑞𝑖 | ≤ 4 and perform both a set of Bayesian-inspired numerical scans and an analytical spurion analysis to identify those charge assignments that reliably generate SM-like quark mass and mixing hierarchies. The resulting set of top-20 flavor charge assignments significantly enlarges the viable space of FN models but is still compact enough to enable focused phenomenological study. We then apply our numerical methodology to demonstrate that these distinct charge assignments result in the generation of correlated flavor-violating four-quark operators characterized by significantly varied strengths, potentially differing substantially from the possibilities previously explored in the literature. Future precision measurement of Δ⁢𝐹 = 2 observables, along with increasingly accurate SM predictions, may therefore enable us to distinguish among otherwise equally plausible FN charges, thus shedding light on the UV structure of the flavor sector.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Fluid-driven fractures in granular media: Insights from numerical investigations

We investigate the mechanisms of opening-mode fracture initiation in granular media. The study is based on a simulation of grain-scale fluid-grain interactions through a coupled numerical approach in which the discrete element method is used to solve for the mechanics of a solid granular medium, and computational fluid dynamics is used to model fluid flow and drag forces. In this work, we present benchmark problems with analytical solutions and validate this numerical model against experiments on a viscous-drag-driven cavity in the literature. Additional simulation results show fracture initiation mechanisms in a random granular packing subjected to constant boundary stresses and to fluid injection with a localized source. The dimensionless variable $F_s/F_{sk}$ (ratio of seepage force $F_s$ and skeletal force $F_{sk}$) incorporates the impacts of physical properties and injection parameters including fluid viscosity, injection velocity, grain size, and effective stresses, and it has been used as a criterion separating regimes of fluid invasion and drag-driven fracture opening. Our simulation results show that $F_s/F_{sk}$ in combination with $τ_1$ (ratio of diffusion time from hydromechanical coupling and injection time) serves as a prediction of fracture opening within granular packing. We suggest a simple criterion ($F_s/F_{sk}$ > 1 or $τ_1$ > 0.17) that is valid for various types of granular media and injection conditions to determine if fracture opening will occur. Among other applications, this study is useful to predict the initiation and propagation of fractures in natural sediments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Synchronous and Concurrent Multidomain Computing Method for Cloud Computing Platforms

We present a numerical method for synchronous and concurrent solution of transient elastodynamics problem where the computational domain is divided into subdomains that may reside on separate computational platforms. Here, this work employs the variational multiscale discontinuous Galerkin (VMDG) method to develop interdomain transmission conditions for transient problems. The fine-scale modeling concept leads to variationally consistent coupling terms at the common interfaces. The method admits a large class of time discretization schemes, and decoupling of the solution for each subdomain is achieved by selecting any explicit algorithm. Numerical tests with a manufactured solution problem show optimal convergence rates. The energy history in a free vibration problem is in agreement with that of the solution from a monolithic computational domain.

97 MATHEMATICS AND COMPUTING↗

pyTCR: A tropical cyclone rainfall model for python

pyTCR is a climatology software package developed in the Python programming language. It integrates the capabilities of several legacy physical models and increases computational efficiency to allow rapid estimation of tropical cyclone (TC) rainfall consistent with the large-scale environment. Specifically, pyTCR implements a horizontally distributed and vertically integrated model [Zhu et al., 2013] for simulating rainfall driven by TCs. Along storm tracks, rainfall is estimated by computing the cross-boundary-layer, upward water vapor transport caused by different mechanisms including frictional convergence, vortex stretching, large-scale baroclinic effect (i.e., wind shear), topographic forcing, and radiative cooling [Lu et al., 2018]. The package provides essential functionalities for modeling and interpreting spatio-temporal TC rainfall data. pyTCR requires a limited number of model input parameters, making it a convenient and useful tool for analyzing rainfall mechanisms driven by TCs. To sample rare (most intense) rainfall events that are often of great societal interest, pyTCR adapts and leverages outputs from a statistical-dynamical TC downscaling model [Lin et al., 2023] capable of rapidly generating a large number of synthetic TCs given a certain climate. As a result, pyTCR significantly reduces computational effort and improves the efficiency in capturing extreme TC rainfall events at the tail of the distributions from limited datasets. Furthermore, the TC downscaling model is forced entirely by large-scale environmental conditions from reanalysis data or coupled General Circulation Models (GCMs), simplifying the projection of TC-induced rainfall and wind speed under future climate using pyTCR. Finally, pyTCR can be coupled with hydrological and wind models to assess risks associated with independent and compound events (e.g., storm surges and freshwater flooding).

54 ENVIRONMENTAL SCIENCES↗

Roadmap for thermal property measurements of Molten Salt Reactor systems

Thermal properties of molten salts, both with and without fissile and fertile material, are fundamental to modeling and simulation that supports design, construction, operation, and accident progression evaluation of molten salt reactors. This roadmap provides guidance to the research community on the types and amounts of required data for base, or carrier, salt systems, recommends simulated burnup fuel salt compositions to evaluate, and provides a technical foundation for a focused salt property evaluation matrix. The property values generated from both computation and experimental methods are used to develop models that constitute the Molten Salt Thermal Properties Database (MSTDB). While the building blocks of MSTDB are fundamental in nature, its data facilitates constructing engineering models that allow stakeholders to investigate the behavior of specific compositions of interest and supports broader modeling and simulation of MSRs through coupling to multi-physics, multi-scale mass accountancy tools.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Frequency Comb-Based Remote Methane Observation Network (Final Scientific/Technical Report)

ARPA-E investment in long-range, frequency comb-based methane detection has resulted in the first scalable, cost-effective continuous methane emissions system. At the onset of this ARPA-E MONITOR award, the dual-frequency comb spectrometer was a sprawling, expensive, laboratory-confined device which had only recently been demonstrated capable of measuring atmospheric trace gases over open-paths. Bringing the technology to the point of being able to monitor for and characterize methane emissions at oil and natural gas production facilities seemed an almost impossible task. Over the subsequent months and years, our team transitioned the Nobel-prize-winning frequency-comb-laser spectrometer technology to a field-deployed regional methane leak detection system. The spectrometer enables highly sensitive near-infrared absorption measurements of methane along kilometer-scale laser beam paths. The measurements are coupled with an atmospheric modeling and inversion framework to triangulate the location of emission sources and quantify the emission rates.

03 NATURAL GAS↗

Advanced Processing of Rare Earth Elements and Critical Minerals from Acid Mine Drainage Feedstocks

Ongoing work by the project team identified a process for extracting rare earths and critical mineral (REE/CM) concentrates from acid mine drainage from coal mines, a major stream pollutant. The process removes harmful metals and separates REE/CM into economically attractive concentrates. This project explored the processing of those concentrates into individually separated, high purity metals, suitable for domestic manufacturers. The team brought unique experience in working with AMD-based feedstocks since the inception of USDOE’s coal-based REE/CM program in 2016. The project team has pioneered the recovery of rare earth/critical minerals (REE/CM) from acid mine drainage (AMD). This environmentally beneficial strategy will incentivize AMD treatment while recovering materials critical to the Nation’s defense and advanced manufacturing industries. In addition, the project will investigate the co-production of five high purity CMs from AMD, namely Co, Mn, Li, Cr, Ge. This team’s upstream and midstream recovery and processing technology has been demonstrated at the continuous, small scale pilot scale. In conjunction with the West Virginia Department of Environmental Protection, the team is constructing a full-scale, integrated, 500 gpm AMD treatment plant coupled with an REE/CM facility that will produce 21 tons REE/CM per year. Commissioning is scheduled for September 2021. This will provide feedstock for the team’s proprietary downstream processes that met the project objective: production of individually separated, high purity (ISHP) REE/CM oxides and their further processing to elemental metals (REM). To facilitate entry into the U.S. market, we have teamed up with companies that would use our products to strengthen their domestic supply chains and grow REE/CM based manufacturing in the U.S. They include Rockwell Automation, Rivian, Endress+Hauser, Inc., Northrop Grumman and Hela Novel Metals.

01 COAL, LIGNITE, AND PEAT↗

Assessing Impacts of Waves on Hub-Height Winds off the U.S. West Coast Using Lidar Buoys and Coupled Modeling Approaches

Given the importance of offshore wind energy development to the U.S. clean energy targets, it is vital to be able to characterize the wind resource in that environment accurately. Toward that end, two Bureau of Ocean Energy Management buoys equipped with Doppler lidar are being maintained by Pacific Northwest National Laboratory on behalf of the Department of Energy and deployed to regions of potential offshore wind development. In addition to standard meteorological and oceanographic measurements, the buoys document the wind profile between about 40 m and 250 m above the sea surface through Doppler lidar retrievals. After a multiyear deployment of two buoys along the U.S. East Coast, the buoys were redeployed to the U.S. West coast from 2020 – 2022 to locations near the Humboldt and Morro Bay lease areas. The buoys provide nearly continuous, multiyear datasets that can be used to evaluate predictions of hub-height (~100 m) wind speed for standard atmospheric models in the region. In the absence of measurements at the study site, offshore wind developers rely on model-based data to assess site conditions. Potential sources of model error in this environment include under-resolution or misrepresentation of coastal topographically forced flows, marine boundary layer dynamics and the evolution of their associated cloud and turbulence fields, the role of upwelling and other currents on surface heat fluxes into the boundary layer, and the impact of wave fields on surface momentum fluxes and thus the wind speed profile. In particular, most predictive models of wind speed do not predict wave fields at all, relying on parameterizations to represent their effects. In thus study, we focus on evaluating the role of wind / wave interactions on modeled hub-height wind speed and error by using the Coupled Ocean–Atmosphere–Wave–Sediment–Transport Modeling System to capture two-way interactions between an atmospheric model (Weather Research and Forecasting (WRF)) and a wave model (WAVEWATCHIII (WW3)) and compare to both stand-alone WRF and one-way coupled WRF / WW3 configurations. Our approach is similar to that used in Gaudet et al. (2022) to evaluate wind / wave coupling over the U.S. East Coast, but applied to the very different environment of the U.S. West Coast. We show examples for two cases, a cold-season frontal case and a warm-season low-level jet case. We find that wind / wave coupling makes little impact on model error for these cases at the location of the lidar buoys, for which other misrepresentations of model physics seems to be responsible for model-observation discrepancies. However, domain-wide evaluations, which also make use of the National Buoy Data Center network, show that a two-way coupling approach is less prone to systematic errors in the hub-height wind field than the one-way coupled approach. WRF resolution of kilometer-scale or less is needed to properly capture the sharp wind speed gradients that can be found along the coastline, and WW3 simulations driven by the downscaled WRF produce better bulk and spectral wave fields when compared to observations. Implications of the results for wind resource characterization are then discussed.

17 WIND ENERGY↗

Time-dependent phenomena in correlated materials

Understanding time-dependent processes and light-matter interaction in strongly correlated materials, and the interplay between electronic, orbital, vibrational, and spin degrees of freedom, is a cornerstone of condensed matter. These mechanisms can be proven by measuring the response of the systems to time-dependent perturbations. The corresponding time scales are dictated by the way light couples to the different excitations, and how these excitations exchange energy and momentum. Our research advances our understanding of these processes, and the interpretation of different equilibrium and time-resolved spectroscopies. Our project encompasses two main themes: (i) developing and refining computational techniques to study non-equilibrium spectroscopies including non-perturbative effects and (ii) applications to non-equilibrium phenomena. We have developed a new computational approach that works directly in the time domain: by including all the degrees of freedom involved in the scattering process (e.g. electrons, photons, neutrons), we solve the time dependent problem: a faithful numerical simulation of the experiment. By measuring the energy and momentum of the outgoing particles, we can extract information about the energy and momentum absorbed by the system. Prior to our work, people attempting to model and calculate non-equilibrium spectral functions relied on a description of the scattering cross section based on a formulation in the frequency domain, a treatment that is extremely cumbersome and complex. Our technique works in and out of equilibrium and can reproduce spectra by several spectroscopic techniques, such as time-resolved photoemission, neutron scattering, Raman, X-ray spectroscopies (RIXS, Auger, XAS, XMCD), and, by not relying on analytical approximations, yields results that reveal novel overlooked transient mechanisms. These tools provide sorely needed intuition for understanding the phenomenology of strongly correlated materials and will help experimentalists in identifying signatures of relevant excitations in pump-probe experiments, such as those conducted in DOE supported facilities.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Increased Crystal Field Drives Intermediate Coupling and Minimizes Decoherence in Tetravalent Praseodymium Qubits

Crystal field (CF) control of rare-earth (RE) ions has been employed to minimize decoherence in qubits and to enhance the effective barrier of single-molecule magnets. The CF approach has been focused on the effects of symmetry on dynamic magnetic properties. Herein, the magnitude of the CF is increased via control of the RE oxidation state. The enhanced 4f metal-ligand covalency in Pr 4+ gives rise to CF energy scales that compete with the spin-orbit coupling of Pr 4+ and thereby shifts the paradigm from the ionic ζ SOC >> V CF limit, used to describe trivalent RE-ion, to an intermediate coupling (IC) regime. We examine Pr 4+ -doped perovskite oxide lattices (BaSnO 3 and BaZrO 3 ). These systems are defined by IC which quenches orbital angular momentum. Therefore, the single-ion spin-orbit coupled states in Pr 4+ can be chemically tuned. We demonstrate a relatively large hyperfine interaction of A iso = 1800 MHz for Pr 4+ , coherent manipulation of the spin with Q M = 2Ω R $T$ m reaching up to ~400 for 0.1Pr:BSO at $T$ = 5 K, and significant improvement of the temperature at which $T$ m is limited by $T$ 1 ($T$* = 60 K) compared to other RE ion qubits.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Representing DC-Coupled PV+Battery Hybrids in a Capacity Expansion Model

Recent technology cost and performance improvements and the federal investment tax credit (ITC) have driven growing interest in coupling solar photovoltaic (PV) and battery systems. Combining these technologies into co-located or hybridized PV+battery systems has the potential to lower costs and increase energy output relative to multiple independent systems. In this work, we provide an overview of PV+battery systems and demonstrate methods for incorporating them into NREL’s Regional Energy Deployment System (ReEDS) capacity expansion model. Although the methods are applied to a specific model, we anticipate that the approaches used here can be useful for informing PV+battery method development for other capacity expansion models. The implemented method relies heavily on the native representations of PV and battery technologies; therefore, the focus of this work is on capturing and parameterizing the interactions between them for a configuration in which the PV and battery technologies share a single bi-directional inverter. This work also demonstrates the impacts of including PV+battery systems in the ReEDS optimization for the conterminous United States through 2050. In particular, we perform parametric sensitivities for input assumptions that are uncertain and expected to influence PV+battery deployment levels, including (a) the cost of PV+battery systems relative to independent PV and battery systems, (b) the battery component’s qualification for the ITC, and (c) future cost trajectories for PV and battery systems. We find that PV+battery deployment could occur throughout the conterminous United States if there are cost savings associated with DC coupling PV and battery technologies. If even modest (5%) cost savings can be achieved (through a shared inverter and balance-of-system costs), then approximately one-third of utility-scale PV deployment through 2050 adopts the DC-coupled hybrid configuration, resulting in total PV+battery deployment that exceeds the magnitude of PV+battery projects in U.S. interconnection queues in 2020. If greater cost savings can be achieved through DC coupling (e.g., due to a growing amount of shared balance-of-system costs, reduced financial risk, or modularity) or more rapid cost and performance improvements are realized for PV and battery technologies, then total PV+battery deployment and the share of PV and battery deployment that adopts the hybrid configuration grows (to >50%). In all cases, growing PV+battery deployment primarily displaces independent PV and battery technologies, indicating the strong competition between the hybrid and independent configurations comprising technologies with similar performance characteristics.

14 SOLAR ENERGY↗

Massively parallel phase-field simulations targeting exascale

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.

Exascale↗

Multirate partitioned Runge–Kutta methods for coupled Navier–Stokes equations

Earth system models are complex integrated models of atmosphere, ocean, sea ice, and land surface. Coupling the components can be a significant challenge due to the difference in physics, temporal, and spatial scales. Further, this study explores multirate partitioned Runge-Kutta methods for the fluid-fluid interaction problem and demonstrates its parallel performance by using the PETSc library. We consider compressible Navier-Stokes equations with gravity coupled through a rigid-lid interface. Our large-scale numerical experiments reveal that multirate partitioned Runge-Kutta coupling schemes (1) can conserve total mass; (2) have second-order accuracy in time; and (3) provide favorable strong- and weak-scaling performance on modern computing architectures. We also show that the speedup factors of multirate partitioned Runge-Kutta methods match theoretical expectations over their base (single-rate) method.

54 ENVIRONMENTAL SCIENCES↗

Coupled Monte Carlo Transport and Conjugate Heat Transfer for Wire-Wrapped Bundles Within the MOOSE Framework

Cardinal is an open-source application that couples OpenMC Monte Carlo transport and NekRS computational fluid dynamics (CFD) to the Multiphysics Object-Oriented Simulation Environment (MOOSE), closing neutronics and thermal-fluid gaps in conducting high-resolution multiscale and multiphysics analyses of nuclear systems. We first provide a brief introduction to Cardinal's software design, data mapping, and coupling strategy to highlight our approach to overcoming common challenges in high-fidelity multiphysics simulations. Here we then present two Cardinal simulations for hexagonal pin bundles. The first is a validation of Cardinal's conjugate heat transfer coupling of NekRS's Reynolds-Averaged Navier Stokes model with MOOSE's heat conduction physics for a bare seven-pin Freon-12 bundle flow experiment. Predictions for pin surface temperatures under three different heating modes agree reasonably well with experimental data and similar CFD modeling from the literature. The second simulation is a multiphysics coupling of OpenMC, NekRS, and BISON for a reduced-scale, seven-pin wire-wrapped version of an Advanced Burner Reactor bundle. Wire wraps are approximated using a momentum source model, and coupled predictions are provided for velocity, temperature, and power distribution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗