Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reference models”

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 1,351 records · Page 75

Preface for the special topic collection honoring Dr. Scott Chambers’ 70th birthday and his leadership in the science and technology of oxide thin films

It is an honor to dedicate this special issue to Dr. Scott A. Chambers, who has had a rewarding and impactful career in surface science, spectroscopy, and thin film synthesis. His research career, spanning from his graduate work in the 1970’s to the present day, was built upon pioneering early work in precision thin film synthesis and spectroscopic characterization that occurred beginning in the 1960’s. Notably, this includes the contributions of both Art Gossard to precision film synthesis by molecular beam epitaxy (MBE) and Chuck Fadley to photoelectron spectroscopy; both Art and Chuck were recently honored with JVSTA commemorative issues of their own. Yet Scott is no mere copycat; he extended and expanded their contributions to further advance the field of surface science, and he applied the same scientific rigor to the emerging field of precision epitaxial oxide synthesis. This rigor was perhaps not always appreciated by the more “enthusiastic” members of the community who tended to draw exciting conclusions from limited data. He was once referred to, fondly, by a collaborator as a “spoilsport” for his penchant for using careful, defensible synthesis and characterization to prove that popular models and assumptions of the day did not stand up to scrutiny.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Prospects for macroscopic dark matter detection at space-based and suborbital experiments

Here we compare two different formalisms for modeling the energy deposition of macroscopically sized/massive quark nuggets (also known as macros) in the Earth's atmosphere. We show that for a reference mass of 1 g, there is a discrepancy in the macro luminosity of about 14 orders of magnitude between the predictions of the two formalisms. Armed with our finding we estimate the sensitivity for macro detection at space-based (Mini-EUSO and POEMMA) and suborbital (EUSO-SPB2) experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Coupled Neutronic and Thermal Hydraulic Analysis of a Natural Circulation Based Small Modular Reactor (SMR) Using VERA-CS

As part of the work supported by a US Department of Energy (DOE) Office of Nuclear Energy Gateway for Accelerated Innovation in Nuclear FY 2020 Voucher, Holtec International subsidiary SMR LLC (hereinafter referred to as Holtec) and Oak Ridge National Laboratory entered into a cooperative research and development agreement (CRADA) to develop coupled multiphysics core models of the Holtec developed small modular reactor (SMR), SMR-160. The scope of the 1 year project was to use the DOE–developed tool set, VERA, to analyze several aspects of a representative SMR-160 core design. The goals of the project were to perform the code-to-code benchmarking and to provide Holtec with a confirmatory calculations to compare against the results with the codes used in the SMR-160 design, identify and resolve areas of discrepancy between the two models to give Holtec a basis for selecting certain tool and model options, and develop a workflow so that the VERA model could be adapted easily to any future changes in the design of the SMR-160 core. The work proceeded in four phases: 1. Developing the initial VERA model based on Holtec-provided core design information; 2. Comparing VERA results with Holtec-provided calculations for the first five cycles; 3. Developing a VERA transient model based on the prior model; and 4. Training Holtec personnel on VERA and hand-off the developed models. This report summarizes at a high level the completed tasks. Details of the tasks and results are reserved for the full report, which will not be publicly released for 5 years per the terms of the CRADA. The report is titled “Coupled Neutronic and Thermal Hydraulic Analysis of a Natural Circulation Based Small Modular Reactor (SMR) using VERA-CS” and has the ORNL technical report number ORNL/TM 2021/2303 (NFE-20-08305)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Considerations for Improving the Technoeconomic Viability of Heat Pipe Microreactors by Employing a Functional Containment Approach

The economic cost of a microreactor is being investigated by leveraging MARVEL cost information in conjunction with some simplified design and analysis activities to present a hypothetical, yet feasible, option for investigation. This bottom-up cost estimate for a hypothetical commercial microreactor will be used to provide guidance to the industry and research communities regarding where priorities should be set to improve the viability of broad microreactor deployment. To support the design and analysis activities required to present a feasible option, some understanding of maximum accidental radionuclide releases and consequences is needed to provide bounding estimates on containment or confinement systems for radionuclide retention. To estimate fission product barrier and retention performance, a gaussian plume dose and dispersion model is employed, with isotopic inventories being propagated between barriers assuming conservative fractional release rates. Demonstrably conservative meteorological assumptions are also assumed. Standard reference release rates are assumed where safety requirements could be easily and economically employed. Results are provided in units of total effective dose equivalent as a function of distance from the release (i.e., the reactor). Initial results indicate that the hypothetical microreactor will likely need to credit the fuel cladding and at least one other barrier, or combination of barriers, to ensure public health and safety as required during a postulated accident. This approach is referred to as “functional containment,” since no single barrier will be relied upon for satisfying regulatory requirements and fundamental safety functions for radionuclide retention. Finally, a set of barrier options are provided so that future economic analyses can select the most viable design path.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Synergistic Retrieval of Temperature and Humidity Profiles from Space-Based and Ground-Based Infrared Sounders Using an Optimal Estimation Method

The atmospheric temperature and humidity profiles of the troposphere are generally measured by radiosondes and satellites, which are essential for analyzing and predicting weather. Nevertheless, the insufficient observation frequencies and low detection accuracy of the boundary layer restricts the description of atmospheric state changes by the temperature and humidity profiles. Therefore, this work focus on retrieving the temperature and humidity profiles using observations of the FengYun-4 (FY-4) Geostationary Interferometric Infrared Sounder (GIIRS) combined with ground-based infrared spectral observations from the Atmospheric Emitted Radiance Interferometer (AERI), which are more accurate than space-based individual retrieval results and have a wider effective retrieval range than ground-based individual retrieval results. Based on the synergistic observations, which are made by matching the space-based and ground-based data with those of different spatial and temporal resolutions, a synergistic retrieval process is proposed to obtain the temperature and humidity profiles at a high frequency under clear-sky conditions based on the optimal estimation method. In this research, using the line-by-line radiative transfer model (LBLRTM) as the forward model for observing simulations, a retrieval experiment was carried out in Qingdao, China, where an AERI is situated. Taking radiosonde data as a reference for comparing the retrieval results of the temperature and humidity profiles of the troposphere, the root-mean-square error (RMSE) of the synergistic retrieval algorithm below 400 hPa is within 2 K for temperature and within 12% for relative humidity. Compared with the GIIRS individual retrieval, the RMSE of temperature and relative humidity for the synergistic method is reduced by 0.13–1.5 K and 2.7–4.4% at 500 hPa, and 0.13–2.1 K and 2.5–7.2% at 900 hPa. Moreover, the forecast index (FI) calculated from the retrieval results shows reasonable consistency with the FIs calculated from the ERA5 reanalysis and from radiosonde data. The synergistic retrieval results have higher temporal resolution than space-based retrieval results and can reflect the changes in the atmospheric state more accurately. Overall, the results demonstrated the promising potential of the synergistic retrieval of temperature and humidity profiles at high accuracy and high temporal resolution under clear-sky conditions from FY-4/GIIRS and AERI.

54 ENVIRONMENTAL SCIENCES↗

Multifidelity Active Learning for Failure Estimation of TRISO Nuclear Fuel

The Tristructural isotropic (TRISO)-coated particle fuel is a robust nuclear fuel proposed to be used for multiple modern nuclear technologies. Therefore, characterizing its safety is vital for the reliable operation of nuclear technologies. However, the TRISO fuel failure probabilities are small and the computational model is time consuming to evaluate them using traditional Monte Carlo-type approaches. In the paper, we present a multifidelity active learning approach to efficiently estimate small failure probabilities given an expensive computational model. Active learning suggests the next best training set for optimal subsequent predictive performance and multifidelity modeling uses cheaper low-fidelity models to approximate the high-fidelity model output. After presenting the multifidelity active learning approach, we apply it to efficiently predict TRISO failure probability and make comparisons to the reference results.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Toroidal modeling of 3D perturbations generated by current filaments in scrape-off layer in tokamak with biased divertor targets

The divertor biasing technique offers a promising alternative to control the edge localized mode (ELM) as well as the divertor heat load in tokamaks, as compared to the resonant magnetic perturbation (RMP) generated by magnetic coils. The linear resistive MHD code MARS-F (Liu et al 2000 Phys. Plasmas 7 3681) is employed to study the plasma response to the n = 2 (n is the toroidal mode number) field perturbation in HL-2A, produced by the divertor biasing current filaments flowing in the scrape-off layer region. The response field due to biasing currents is compared with the RMP field utilized for ELM control in HL-2A. The strength of the radial resonant field perturbation, produced by the biasing currents of 100 A level, is found to be comparable with RMP generated by several kA of ELM control coil currents for the reference plasma. The plasma normal displacement near the X-point and the associated neoclassical toroidal viscosity torque are also computed to be similar between these two techniques. The modeling results thus strongly suggest that the biasing technique can be applied to control ELMs. Moreover, the biasing currents produce field perturbations, including the plasma response, that are localized more near the plasma edge (compared to the RMP counterpart), thus reducing the chance of mode locking associated with core perturbations. Particle orbit tracing also reveals that the biasing current produced magnetic perturbation tends to widen the heat deposition region and induce the strike point splitting of the ion saturation flow on the outer divertor surface, consistent with experimental observations in HL-2A. These toroidal modeling results confirm the possibility of ELM control and plasma exhaust solution by the divertor biasing technique.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Testing Hidden Assumptions of Representativeness in Reach-Scale Studies of Hyporheic Exchange

Field studies of hyporheic exchange in mountain systems are often conducted using short study reaches and a limited number of observations. It is common practice to assume these study reaches represent hyporheic exchange at larger scales or different sites and to infer general relationships among potential causal mechanisms from the limited number of observations. However, these assumptions of representativeness are rarely tested. In this study, we develop numerical models from four segments of mountain streams in different geomorphologic settings and extract shorter reaches to test how representative exchange metrics are in shorter reaches compared to their reference segments. We also map the locations of the representative reaches to determine if a pattern exists based on location. Finally, we compare variance of these shorter within-site reaches to 29 additional reaches across the same basin to understand the impacts of inferring causal mechanisms, for example, the expectation that wide and narrow valley bottoms will yield different hyporheic exchange patterns. Our results show that the location and length strategy of the study reach must be considered before assuming an exchange metric to be representative of anything other than the exact segment studied. Further, it is necessary to quantify within and between site variations before making causal inferences based on observable characteristics, such as valley width or stream morphology. Our findings have implications for future field practices and how those practices are translated into models.

54 ENVIRONMENTAL SCIENCES↗

Assessment of the Griffin Reactor Multiphysics Application Using the Empire Micro Reactor Design Concept

In late 2019, INL and ANL agreed to jointly develop the reactor physics code named Griffin based on the integration of the two code suites, MAMMOTH/Rattlesnake (INL) and MC2 - 3/PROTEUS (ANL). Griffin is being developed based on the MOOSE framework and MOOSE quality assurance procedures. This decision was made to be able to allow DOE-NE to efficiently invest funding to this area and to provide effective and timely support for existing and potential users; the latter includes industry and government organizations who are developing various types of advanced reactors in the near and long term. Since MAMMOTH/Rattlesnake has been developed based on the MOOSE framework, the INL/ANL Griffin development team agreed to build Griffin beginning with a merger of MAMMOTH and Rattlesnake into a single code and moving forward by implementing capabilities from the PROTEUS suite into Griffin. Moving forward, both ANL and INL efforts are equally invested in the Griffin project, with management support, to provide an advanced reactor multiphysics tool to assist in reactor design, optimization, and safety analysis. Much work remains in moving Griffin forward to migrate PROTEUS capabilities and to optimize performance to meet user needs. The main objective of this work is to assess the current status of Griffin capabilities in terms of performance and accuracy, to determine priorities for PROTEUS migration, and to identify capabilities and features to improve for supporting the code integration effort. For this assessment, the Empire micro reactor problem that was developed in the ARPA-E MEITNER program was selected as an advance reactor concept of interest to the technical community. The Empire reactor problem was expanded from its original incomplete specification to be a small heat-pipe-cooled micro reactor core with ~113 cm radius and 70 cm in height, composed of 18 fuel assemblies, 12 control drums, and beryllium radial and axial reflectors. In the current model, using 5 cm axial reflectors specified in the original Empire assembly model, more than 10% of neutrons leak axially and through the empty center safety hole, as well as through heat pipe channels in fuel assembly elements that extend through the top reflector region. Several calculation models of the core were defined for systematic assessment, including 2-D and 3-D fuel assemblies and whole cores with cylindrical boundaries. Cross sections were generated using Serpent 2, and meshes were produced using the Argonne mesh tool or the INL neutronics meshing tools combined with CUBIT. Cross sections and meshes were converted to the ISOXML and Exodus formats, respectively, so that Griffin and PROTEUS could use consistent data for solving the reactor problems. With the prepared cross sections and meshes, PROTEUS was run first to ensure that all input data were correctly generated and input options in terms of angle, mesh, and energy group were accurately determined. Comparisons against Serpent 2 solutions were made in terms of eigenvalue and pin power. The same calculations and comparisons were then conducted using Griffin. For the fuel assembly and whole core problems, the PROTEUS eigenvalues agreed well with reference Serpent 2 solutions within 100 and 30 pcm, respectively, and pin power differences relative to Serpent 2 were overall less than 2.2% and RMS 0.8% for the whole core models. This indicated that all input data were properly prepared. Using the same data, Griffin was run selecting the SAAF-CFEM SN solver with Legendre-Gaussian quadrature and NDA and DSA for acceleration. It was found that the SAAF-CFEM solver of Griffin required finer meshes to achieve eigenvalue and pin power solutions in good agreement with Serpent 2, consequently requiring more memory requirement and longer computation time. On the other hand, the SPH-Diffusion 2-D core calculations performed using Griffin were able to recover the exact eigenvalue from the reference Serpent 2 solutions, resulting in a pin-power distribution with an RMS of 0.6% and maximum absolute difference of less than 1.4%. The runtimes for SPH-Diffusion for the 2-D core were less than 3 minutes on 40 cores. During this evolution of this evaluation, many updates were made in Griffin by the Griffin development team of INL (focusing on software updates) and ANL (reviewing and supporting software updates) to complete this assessment. Observations from the code assessment are presented in the conclusion section of this report, followed by a discussion of recommendations for future work.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

MSD CoP Webinar: Representing Climate Impacts in Scenaros

Context: This webinar was hosted by the MultiSector Dynamics Community of Practice (MSD CoP; https://multisectordynamics.org). Abstract: Scenarios of future emissions and land use have been most commonly produced without accounting for the effects that climate impacts linked to those emissions and land use changes may have. The questions of how important such impacts could be, at the regional or global scale, has become an increasingly cogent one. Using the Global Change Analysis Model (GCAM), we implement climate impacts on water supply, agricultural productivity, and energy demand driven by climatic impact drivers whose evolution over the 21st century is representative of a climate consistent with GCAM's reference scenario and compare the output to the reference at both regional and global scales. We show that the impacts from the version of GCAM that uses GDP as an exogenous input are not sufficient to bend emission pathways at the global scale. However, we see effects emerge at the regional scale in terms of emissions, among other metrics linked to water, land and energy. We expect the impacts to become more significant at an aggregate, global scale in a forthcoming version of the model with endogenous GDP. Meanwhile, we document the mechanisms that drive the difference at regional scales which could still be significant in their impacts for individual economies and populations. Presenters: Dr. Claudia Tebaldi (Joint Global Change Research Institute, Pacific Northwest National Laboratory) Moderator(s): Jennifer Morris (MSD CoP SSG member), Patrick Reed (MSD CoP Facilitation Team Member, Moderator and Organizer) This webinar was held on: November 5, 2024 from 1 PM - 2:15 PM ET

Climate Scenarios↗

Parallel transport modeling of linear divertor simulators with fundamental ion cyclotron heating *

Abstract The Material Plasma Exposure eXperiment (MPEX) is a steady state linear device with the goal to perform plasma material interaction studies at future fusion reactor relevant conditions. A prototype of MPEX referred as ‘Proto-MPEX’ is designed to carry out research and development related to source, heating and transport concepts on the planned full MPEX device. The auxiliary heating schemes in MPEX are based on cyclotron resonance heating with radio frequency (RF) waves. Ion cyclotron heating (ICH) and electron cyclotron heating in MPEX are used to independently heat the ions and electrons and provide fusion divertor conditions ranging from sheath-limited to fully detached divertor regimes at a material target. A hybrid particle-in-cell code- PICOS++ is developed and applied to understand the plasma parallel transport during ICH in MPEX/Proto-MPEX to the target. With this tool, evolution of the distribution function of MPEX/Proto-MPEX ions is modeled in the presence of (a) Coulomb collisions, (b) volumetric particle sources and (c) quasi-linear RF-based ICH. The code is benchmarked against experimental data from Proto-MPEX and simulation data from B2.5 EIRENE. The experimental observation of ‘density-drop’ near the target in Proto-MPEX and MPEX during ICH is demonstrated and explained via physics-based arguments using PICOS++ modeling. In fact, the density drops at the target during ICH in Proto-MPEX/MPEX to conserve the flux and to compensate for the increased flow during ICH. Furthermore, sensitivity scans of various plasma parameters with respect to ICH power are performed for MPEX to investigate its role on plasma transport and particle and energy fluxes at the target.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The mechanisms and seasonal differences of the impact of aerosols on daytime surface urban heat island effect

The urban heat island intensity (UHII) is the temperature difference between urban areas and their rural surroundings. It is commonly attributed to changes in the underlying surface structure caused by urbanization. Air pollution caused by aerosol particles can affect the UHII through changing (1) the surface energy balance by the aerosol radiative effect (ARE) and (2) planetary-boundary-layer (PBL) stability and airflow intensity by modifying thermodynamic structure, which is referred to as the aerosol dynamic effect (ADE). By analyzing satellite data and ground-based observations collected from 2001 to 2010 at 35 cities in China and using the WRF-Chem model, we find that the impact of aerosols on UHII differs considerably: reducing the UHII in summer but increasing the UHII in winter. This seasonal contrast is proposed to be caused by the different strengths of the ARE and ADE between summer and winter. In summer, the ARE on UHII is dominant over the ADE, cooling down surface temperature more strongly in urban areas than in rural areas because of much higher aerosol loading, and offsets the urban heating, therefore weakening UHII. In winter, however, the ADE is more dominant, because aerosols stabilize the PBL more in the polluted condition, weakening the near-surface heat transport over urban areas in both vertical and horizontal directions. This means that the heat accumulated in urban areas is dispersed less effectively, and thus the UHII is enhanced. These findings shed new light on the impact of the interaction between urbanization-induced surface changes and air pollution on urban climate.urban areas and their rural surroundings. It is commonly attributed to changes in the underlying surface structure caused by urbanization. Air pollution caused by aerosol particles can affect the UHII through changing (1) the surface energy balance by the aerosol radiative effect (ARE) and (2) planetary-boundary-layer (PBL) stability and airflow intensity by modifying thermodynamic structure, which is referred to as the aerosol dynamic effect (ADE). By analyzing satellite data and ground-based observations collected from 2001 to 2010 at 35 cities in China and using the WRF-Chem model, we find that the impact of aerosols on UHII differs considerably: reducing the UHII in summer but increasing the UHII in winter. This seasonal contrast is proposed to be caused by the different strengths of the ARE and ADE between summer and winter. In summer, the ARE on UHII is dominant over the ADE, cooling down surface temperature more strongly in urban areas than in rural areas because of much higher aerosol loading, and offsets the urban heating, therefore weakening UHII. In winter, however, the ADE is more dominant, because aerosols stabilize the PBL more in the polluted condition, weakening the near-surface heat transport over urban areas in both vertical and horizontal directions. This means that the heat accumulated in urban areas is dispersed less effectively, and thus the UHII is enhanced. These findings shed new light on the impact of the interaction between urbanization-induced surface changes and air pollution on urban climate.

54 ENVIRONMENTAL SCIENCES↗

Measurement of electrons from beauty-hadron decays in $pp$ and Pb-Pb collisions at $\sqrt{s_{NN}}$ = $5.02$ TeV

The production of electrons from beauty-hadron decays was measured at midrapidity in proton-proton ($pp$) and central Pb-Pb collisions at center-of-mass energy per nucleon-nucleon pair $\sqrt{s_{NN}}$ = 5.02 TeV, using the ALICE detector at the LHC. The cross section measured in p p collisions in the transverse momentum interval $2$ < $p_T$ < $26$ $\mathrm{GeV}$ / $c$ was compared with models based on perturbative quantum chromodynamics calculations. The yield in the 10% most central Pb-Pb collisions, measured in the interval $2$ < $p_T$ < $26$ $\mathrm{GeV}$ / $c$, was used to compute the nuclear modification factor $R_{AA}$, extrapolating the $pp$ reference cross section to $p_T$ larger than 8 $\mathrm{GeV}$ / $c$. The measured $R_{AA}$ shows significant suppression of the yield of electrons from beauty-hadron decays at high $p_T$ and does not show a significant dependence above 8 $\mathrm{GeV}$ / $c$ within uncertainties. The results are described by several theoretical models based on different implementations of the interaction of heavy quarks with a quark-gluon plasma, which predict a smaller energy loss for beauty quarks compared to light and charm quarks.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Void formation in operator growth, entanglement, and unitarity

The structure of the Heisenberg evolution of operators plays a key role in explaining diverse processes in quantum many-body systems. In this paper, we discuss a new universal feature of operator evolution: an operator can develop a void during its evolution, where its nontrivial parts become separated by a region of identity operators. Such processes are present in both integrable and chaotic systems, and are required by unitarity. We show that void formation has important implications for unitarity of entanglement growth and generation of mutual information and multipartite entanglement. We study explicitly the probability distributions of void formation in a number of unitary circuit models, and conjecture that in a quantum chaotic system the distribution is given by the one we find in random unitary circuits, which we refer to as the random void distribution. We also show that random unitary circuits lead to the same pattern of entanglement growth for multiple intervals as in (1 + 1)-dimensional holographic CFTs after a global quench, which can be used to argue that the random void distribution leads to maximal entanglement growth.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Search for charged-lepton flavor violation in top quark production and decay in pp collisions at $ \sqrt{s} $ = 13 TeV

Results are presented from a search for charged-lepton flavor violating (CLFV) interactions in top quark production and decay in pp collisions at a center-of-mass energy of 13 TeV. The events are required to contain one oppositely charged electron-muon pair in the final state, along with at least one jet identified as originating from a bottom quark. The data correspond to an integrated luminosity of 138 fb$^{-1}$, collected by the CMS experiment at the LHC. This analysis includes both the production (q → eμt) and decay (t → eμq) modes of the top quark through CLFV interactions, with q referring to a u or c quark. These interactions are parametrized using an effective field theory approach. With no significant excess over the standard model expectation, the results are interpreted in terms of vector-, scalar-, and tensor-like CLFV four-fermion effective interactions. Finally, observed exclusion limits are set at 95% confidence levels on the respective branching fractions of a top quark to an eμ pair and an up (charm) quark of 0.13 × 10$^{-6}$ (1.31 × 10$^{-6}$), 0.07 × 10$^{-6}$ (0.89 × 10$^{-6}$), and 0.25 × 10$^{-6}$ (2.59 × 10$^{-6}$) for vector, scalar, and tensor CLFV interactions, respectively.[graphic not available: see fulltext]

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

New Technologies for Power System Operation and Analysis, 1st Edition

New Technologies for Power System Operation and Analysis considers the very latest developments in renewable energy integration and system operation, including electricity markets and wide-area monitoring systems and forecasting. Helping readers quickly grasp the essential information needed to address renewable energy integration challenges, this new book looks at basic power system mathematical models, advanced renewable integration and system optimizations from transmission and distribution system sides. Sections cover wind, solar, gas and petroleum, making this a useful reference for all engineers interested in power system operation.

electricity markets↗

Decarbonization scenarios of the U.S. Electricity system and their costs

Decarbonizing the electricity system to zero-carbon emission is crucial for climate change mitigation. Previous studies have shown that such a transition in the United States (U.S.) may lead to higher system cost compared to a business-as-usual case, but it is not well-known how the cost of electricity generation varies at sub-regional level under the transition, and studies have rarely evaluated the trade-off between the cost and avoided climate damages, as well as the potential roles of negative emission technologies (NETs) in the electricity decarbonization. In this work, we present a regionally resolved national model to quantify the cost of decarbonizing the U.S. electricity system under a set of possible scenarios. The results show that, compared to the reference scenario without a decarbonization policy, reaching zero CO 2 emission by 2050 would incur, depending on the scenarios, 335-494 billion USD additional cost to the U.S. electric power sector during 2020-2050. The regional costs of electricity generation ranges from 2.4 to 4.7 cent/kWh, largely due to the generation profiles and renewable resources availability of those regions. The additional costs can be translated to an average CO 2 abatement cost of 29-59 USD/metric ton CO 2 (with 2%-7% discount rates), which are comparable to the social cost of carbon in the literature at around 4% discount rate. The results also show that the cost of mitigating the last few percent CO 2 emission from the U.S. electricity system may exceed the costs of NETs, indicating an opportunity for NETs to contribute to electricity decarbonization.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Solution of the Schrödinger equation for quasi-one-dimensional materials using helical waves

We formulate and implement a spectral method for solving the Schrödinger equation, as it applies to quasi-one-dimensional materials and structures. This allows for computation of the electronic structure of important technological materials such as nanotubes (of arbitrary chirality), nanowires, nanoribbons, chiral nanoassemblies, nanosprings and nanocoils, in an accurate, efficient and systematic manner. Our work is motivated by the observation that one of the most successful methods for carrying out electronic structure calculations of bulk/crystalline systems — the plane-wave method — is a spectral method based on eigenfunction expansion. Our scheme avoids computationally onerous approximations involving periodic supercells often employed in conventional plane-wave calculations of quasi-one-dimensional materials, and also overcomes several limitations of other discretization strategies, e.g., those based on finite differences and atomic orbitals. The basis functions in our method — called helical waves (or twisted waves) — are eigenfunctions of the Laplacian with symmetry adapted boundary conditions, and are expressible in terms of plane waves and Bessel functions in helical coordinates. We describe the setup of fast transforms to carry out discretization of the governing equations using our basis set, and the use of matrix-free iterative diagonalization to obtain the electronic eigenstates. Miscellaneous computational details, including the choice of eigensolvers, use of a preconditioning scheme, evaluation of oscillatory radial integrals and the imposition of a kinetic energy cutoff are discussed. We have implemented these strategies into a computational package called HelicES (Helical Electronic Structure). We demonstrate the utility of our method in carrying out systematic electronic structure calculations of various quasi-one-dimensional materials through numerous examples involving nanotubes, nanoribbons and nanowires. We also explore the convergence properties of our method, and assess its accuracy and computational efficiency by comparison against reference finite difference, transfer matrix method and plane-wave results. We anticipate that our method will find applications in computational nanomechanics and multiscale modeling, for carrying out transport calculations of interest to the field of semiconductor devices, and for the discovery of novel chiral phases of matter that are of relevance to the burgeoning quantum hardware industry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗