Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “transition zone”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Modeling and Simulation of Xe-100-type Pebble Bed Gas-Cooled Reactor with SCALE

The US Department of Energy (DOE) announced the Advanced Reactor Demonstration Program (ARDP) to accelerate the deployment of advanced reactor concepts. Awardees of ARDP funds are expected to demonstrate the operation of an advanced reactor within 7 years of receiving the award. X-Energy’s advanced reactor concept, the Xe-100, was selected as one of two advanced reactor concepts to receive funding to demonstrate the operation of its high-temperature gas-cooled pebble-bed reactor before the end of this decade. As a result of this push to bring advanced reactors to maturation and commercialization, transition and deployment scenario studies are being performed under the Systems Analysis and Integration (SA&I) campaign within the DOE Office of Nuclear Energy (DOE-NE) to evaluate the transition of the current US commercial fleet of light-water reactors (LWRs) to a future fleet of advanced reactors consisting of a mix of ARDP type reactor concepts and advanced LWRs. To accurately evaluate the front- and back-end resource requirements, it is important to perform reactor physics calculations to determine the discharge burnup and isotopic content, fuel residence time, as well as other parameters. For this purpose, a summer project funded by the SA&I campaign allowed for the setup of SCALE models for full-core Xe-100 type high-temperature gas-cooled pebble-bed reactor and a Xe-100 type slice using publicly available information. The core-averaged equilibrium compositions and zone-wise equilibrium compositions for the slice and 3D models, respectively, were obtained following an iterative depletion method developed by Bostelmann et al. using SCALE’s reactor physics sequence TRITON. The slice model was used with TRITON to generate burnup-dependent cross section libraries at different temperatures which can be used with SCALE’s ORIGAMI code to rapidly determine fuel inventory and therefore to perform quick sensitivity studies on parameters such as the pebble location in the core. The SCALE/TRITON transport and depletion calculation for the Xe-100 type slice model indicates that the isotopic concentrations are in good agreement at 1,300 effective full power days (EFPD) for 235 U. An analysis of 236 U results match 239 Pu results would seem to indicate a typographical error in Mulder and Boyes wherein the reported results of 236 U and 239 Pu are reversed. In addition to SCALE/TRITON calculations, a new capability within SCALE/ORIGAMI for the simulation of pebble-bed reactors was used to study the burnup sensitivity with respect to the pebble pathway through the core. The SCALE/ORIGAMI results show that pebbles that travel closer to the reflector for the entire depletion history have a higher burnup than pebbles that travel through the middle of the core because of the higher thermal to fast flux ratio near the reflector. Consequently, a pebble’s burnup is strongly affected by the pebble’s pathway for each pass. Additional phenomena such as temperature distributions in the core and different travel times of the pebbles in the individual radial zones further affect the burnup distribution. The sensitivity of the discharge vector to the pebble pathways taken during each pass can be evaluated in the future using SCALE/ORIGAMI now that the SCALE inputs have been established.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Molecular dynamics simulations of inelastic x-ray scattering from shocked copper

By taking the spatial and temporal Fourier transforms of the coordinates of the atoms in molecular dynamics simulations conducted using an embedded-atom-method potential, we calculate the inelastic scattering of x rays from copper single crystals shocked along [001] to pressures of up to 70 GPa. Above the Hugoniot elastic limit, we find that the copious stacking faults generated at the shock front introduce strong quasi-elastic scattering (QES) that competes with the inelastic scattering signal, which remains discernible within the first Brillouin zone; for specific directions in reciprocal space outside the first zone, the QES dominates the inelastic signal overwhelmingly. The synthetic scattering spectra we generate from our Fourier transforms suggest that energy resolutions of order 10 meV would be required to distinguish inelastic from quasi-elastic scattering within the first Brillouin zone of shock-loaded copper. We further note that high-resolution inelastic scattering also affords the possibility of directly measuring particle velocities via the Doppler shift. These simulations are of relevance to future planned inelastic scattering experiments at x-ray Free Electron Laser facilities.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Degradation in interfacial shear strength of carbon fiber/ vinyl ester composites due to long-term exposure to seawater using push-out tests

Here, in this study, single fiber (~7-μm diameter) push-out tests are conducted to evaluate hygrothermal effects on the interfacial shear strength (IFSS) of carbon fiber/vinyl ester (CF/VE) composites. Hygrothermal conditioning is achieved by saturating samples in simulated seawater at 40 °C for two years. An investigation has been conducted on the preparation, validity, and interpretation of the push-out test results. First, the authors present a polishing methodology that results in thin films of CF/VE composites in the thickness range of 15–120 μm and produces an average 41.2% drop in IFSS due to long-term hygrothermal exposure. Using scanning electron microscopy (SEM), we show that during the push-out tests, the failure initiates locally at the zone of minimum bond strength at the bottom (away from the indenter), then propagates along the length of the interface. The influence of radial tensile stresses originating due to bending is found to be negligible. Using the SEM imaging of the pushed-out fibers, we validate the failure of the interface to be the primary source of failure. The associated results are found to depend on the thickness of the interface. We then reevaluate the results using the Weibull distribution, knowing that the failure mechanism is analogous to the weakest link theory. The results show a 25.5% drop in the IFSS of the CFVE composite, measured at an infinitesimal scale due to long-term to hygrothermal conditioning at 40 °C. A significant drop in IFSS was observed after reheating above glass transition temperature (T g ) and cooling.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Evolution of storage monitoring – update in response to commercial and regulatory drivers

Carbon Capture and Storage (CCS) is in transition from first-of-a kind projects and research-orientated pilots to commercially-motivated applications. Monitoring results from many newly developed and planned large scale commercial projects are limited; however, it is worthwhile to assess their evolution and consider new strategies as part of an effort to assess and document best practices. Commercial monitoring is targeted to activities that comply with regulatory drivers and de-risk investments. Commercial monitoring also supports accounting that storage has occurred and is tied to project financing. It deals with long time frames and large volumes injected into multiple wells and multiple projects in favorable areas. We see developing trends toward reproducible workflows that systematically reduce risks and clarify expectations for oversight and long-term surveillance. Monitoring techniques showing increasing trends include injection zone pressure as a history-matching and compliance tool. To reduce cost and environmental impact of time-lapse seismic data collection, deploying new approaches and tools, such as use of fibre and installed sources are increasingly applied. Concern over the risk of induced seismicity by regulatory bodies and the general public has increased, which has also resulted in increased monitoring. Some techniques used in the early research phases have been sidelined or used only in restricted applications. For example, geochemical analyses in the injection zone as well as the environment are now being deployed less than it was in research-oriented programs, except in the US where it is required by the permitting process. Expectations of frequent area-wide near surface monitoring have also decreased.

25 ENERGY STORAGE↗
DOE OSTI2021

Chemical Compositions of Red Giant Stars from Habitable Zone Planet Finder Spectroscopy

We have used the Habitable Zone Planet Finder (HPF) to gather high-resolution, high signal-to-noise near-infrared spectra of 13 field red horizontal branch (RHB) stars, one open cluster giant, and one very metal-poor halo red giant. The HPF spectra cover the 0.81–1.28 μm wavelength range of the zyJ bands, partially filling the gap between the optical (0.4–1.0 μm) and infrared (1.5–2.4 μm) spectra already available for the program stars. We derive abundances of 17 species from LTE-based computations involving equivalent widths and spectrum syntheses, and estimate abundance corrections for the species that are most affected by departures from LTE in RHB stars. Generally good agreement is found between HPF-based metallicities and abundance ratios and those from the optical and infrared spectral regions. Light element transitions dominate the HPF spectra of these red giants, and HPF data can be used to derive abundances from species with poor or no representation in optical spectra (e.g., C i, P i, S i, K i). Attention is drawn to the HPF abundances in two field solar-metallicity RHB stars of special interest: one with an extreme carbon isotope ratio, and one with a rare, very large lithium content. The latter star is unique in our sample in exhibiting very strong He i 10830 Å absorption. The abundances of the open cluster giant concur with those derived from other wavelength regions. Detections of C i and S i in HD 122563 are reported, yielding the lowest metallicity determination of [S/Fe] from more than one multiplet.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Spin excitations in the kagome-lattice metallic antiferromagnet Fe 0.89 Co 0.11 Sn

Kagome-lattice materials have attracted tremendous interest due to the broad prospect for seeking superconductivity, quantum spin liquid states, and topological electronic structures. Among them, the transition-metal kagome lattices are high-profile objects for the combination of topological properties, rich magnetism, and multiple-orbital physics. Here we report an inelastic neutron scattering study on the spin dynamics of a kagome-lattice antiferromagnetic metal Fe 0.89 Co 0.11 Sn. Although the magnetic excitations can be observed up to ~250 meV, well-defined spin waves are only identified below ~90 meV and can be modeled using Heisenberg exchange with ferromagnetic in-plane nearest-neighbor coupling J 1 , in-plane next-nearest-neighbor coupling J 2 , and antiferromagnetic (AFM) interlayer coupling J c under linear spin-wave theory. Above ~90 meV, the spin waves enter the itinerant Stoner continuum and become highly damped particle-hole excitations. At the K point of the Brillouin zone, we reveal a possible band crossing of the spin wave, which indicates a potential Dirac magnon. Finally, our results uncover the evolution of the spin excitations from the planar AFM state to the axial AFM state in Fe 0.89 Co 0.11 Sn, solve the magnetic Hamiltonian for both states, and confirm the significant influence of the itinerant magnetism on the spin excitations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

The surface triple-coupling on single crystalline cathode for lithium ion batteries

Single crystalline (SC) cathode materials, which are less susceptible to micro/nano-cracks formation and offer better structure stability compared to the polycrystalline counterpart, have attained great attention. However, the parasitic side reactions at the cathode-electrolyte interface induces the loss of active species, which consequently leads to continual degradation of the electrochemical performances. Herein, a triple coupling of concentration-gradient Na+, F- co-doping and surface NaF coating are exploited for the first time on SC LiNi 0.5C o 0.2 Mn 0.3 O 2 cathode by the hydrolysis of NaPF 6 . This process regulates the external structure of materials by constructing a “sandwich” configuration from surface to bulk: rock salt - mixing zone - layered phase. The detailed interface transformation mechanism is revealed by Neutron powder diffraction (NPD), spherical aberration corrected high-resolution scanning transmission electron microscopy (HR-STEM), electron energy loss spectroscopy (EELS), and Ar+ sputtering assisted X-ray photoelectron spectroscopy (XPS). The synergistic effects endow the SC cathode with outstanding capacity retentions: 91.3% at 25 °C and 85% at 45 °C, after 500 cycles at 5 C between 3.0 and 4.5 V. In addition, a high full-cell reversible capacity of 168.9 mAh g -1 with a capacity retention of 92.4% is achieved after 300 cycles at 1 C. Multiple characterizations further indicate that these superior results are mainly ascribed to the overall structure integrity of SC material, the thin cathode electrolyte interface, high content of lithium fluoride, and the low solubility of transition metal ions. This work opens a new avenue to construct a benign interface towards high-performance lithium ion batteries.

lithium ion batteries↗
DOE OSTI2023

Attaining freshwater and estuarine-water soil saturation in an ecosystem-scale coastal flooding experiment

Abstract Coastal upland forests are facing widespread mortality as sea-level rise accelerates and precipitation and storm regimes change. The loss of coastal forests has significant implications for the coastal carbon cycle; yet, predicting mortality likelihood is difficult due to our limited understanding of disturbance impacts on coastal forests. The manipulative, ecosystem-scale Terrestrial Ecosystem Manipulation to Probe the Effects of Storm Treatments (TEMPEST) experiment addresses the potential for freshwater and estuarine-water disturbance events to alter tree function, species composition, and ecosystem processes in a deciduous coastal forest in MD, USA. The experiment uses a large-unit (2000 m 2 ), un-replicated experimental design, with three 50 m × 40 m plots serving as control, freshwater, and estuarine-water treatments. Transient saturation (5 h) of the entire soil rooting zone (0–30 cm) across a 2000 m 2 coastal forest was attained by delivering 300 m 3 of water through a spatially distributed irrigation network at a rate just above the soil infiltration rate. Our water delivery approach also elevated the water table (typically ~ 2 m belowground) and achieved extensive, low-level inundation (~ 8 cm standing water). A TEMPEST simulation approximated a 15-cm rainfall event and based on historic records, was of comparable intensity to a 10-year storm for the area. This characterization was supported by showing that Hurricane Ida’s (~ 5 cm rainfall) hydrologic impacts were shorter (40% lower duration) and less expansive (80% less coverage) than those generated through experimental manipulation. Future work will apply TEMPEST treatments to evaluate coastal forest resilience to changing hydrologic disturbance regimes and identify conditions that initiate ecosystem state transitions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Hawai‘i Play Fairway (Final Report)

Most of Hawai'i's geothermal resources are blind—their manifestations, such as hot springs and steam vents, do not appear on the ground surface because the heated water flows far below. With the exception of $K\bar{i}lauea East Rift Zone$, in most areas of Hawai'i, high lateral permeability in the first kilometer below ground surface prevents surface thermal features from developing. As a methodology for discovering these blind resources, Play Fairway Analysis (PFA) involves finding potential locations of blind hydrothermal systems and describing potential geothermal sources in rift-zone settings. Using the PFA to find Hawai'i's geothermal resources, the University of Hawai'i (UH) conducted the Hawai‘i Play Fairway Project, Hawai'i's first statewide geothermal resource assessment since 1985. Sponsored by the U.S. Department of Energy, the Hawai'i Play Fairway Project provided an updated resource assessment, a roadmap for additional exploration activities, and the identification of areas for further exploration. Benefitting from UH's core competency in earth sciences and experienced geothermal researchers, the project comprised three phases. During the first phase, the team identified, compiled, and ranked existing geologic, groundwater, and geophysical datasets relevant to subsurface heat, fluid and permeability. Using a Bayesian statistical approach, the team developed a statistical methodology to integrate these data into a resource probability map. The team evaluated the confidence in the probability value and considered development viability of areas with geothermal resources. With these analyses, the team identified 10 locations in the Hawaiian Islands for exploration activities. For the second phase, the team collected new groundwater data in 10 locations across the state and new geophysical data on $L\bar{a}na‘i, Maui$, and central Hawai'i Island and modeled topographically induced stress to better characterize subsurface permeability. Analyzing the subsurface stresses, the team evaluated the potential for fracture-induced permeability. The team inverted the MT and gravity data to produce 3D models of resistivity and density, respectively, on $L\bar{a}na‘i$, across $Haleakal\bar{a}'s$ SW rift (Maui), and surrounding Mauna Kea (Hawai‘i Island). The team developed and applied a new method for incorporating depth information about resistivity, density, and potential for fracture-induced permeability into the statistical method for computing resource probability in these three focus areas. The team incorporated the new groundwater results with the new geophysical results and the calculations of potential for fracture-induced permeability to produce updated maps of resource probability and confidence. Through combining data from the first and second phases, the team determined locations for further exploration during the third phase. For MT and gravity surveys, the team recommended $Kaua'i's$ $L\bar{i}hu'e$ $Basin$, the east rift of $Maui's$ $Haleakal\bar{a}$ volcano, and the southwest rift of Hawai'i Island's Mauna Loa volcano. The MT and gravity surveys aimed to enable improved confidence in the resource potential in these locations. For drilling deep groundwater well(s), the team recommended Southeast Mauna Kea and $L\bar{a}na's$ $P\bar{a}l\bar{a}wai$ $Basin$. During the third phase, further exploration involved drilling a groundwater well in $L\bar{a}na's$ $P\bar{a}l\bar{a}wai$ $Basin$ and performing more geophysical surveys. We deepened an existing water well proximal to our target area on $L\bar{a}na'i$ due to funding constraints that precluded us from spudding a new well that would exceed 1km depth. Drilling was preceded by a number of substantial elements including: writing an Environmental Assessment and the subsequent legal process, performance of deviation logging, lowering a camera down the well, coordinating site preparation with $P\bar{u}lama$ $L\bar{a}na'i$, shipping the UH-owned rig interisland, procuring supplies, and leading 3 community meetings on $L\bar{a}na'i$. Drilling occurred 24/7 the entire month of June 2019 over which time $L\bar{a}na'i$ $Well$ 10 was deepened from 427 m to 1057 m, with continuous core collected. We measured a roughly linear temperature gradient averaging 42°C/km and a maximum bottom hole temperature of 66°C. This gradient is more than twice the background for Hawai'i and within a range of gradients measured in this depth range for some exploration wells within KERZ. We consider these results encouraging for $L\bar{a}na'i's$ resource potential and recommend following with a slim hole within $L\bar{a}na'i's$ caldera (our target zone) to ~ 2 km. Further, the positive implications such results have for the island of O‘ahu are substantial - the shield stage of O'ahu's volcanoes ended 1-2 My earlier. However, O'ahu uses more electricity than the rest of the islands combined, and the utility recently called for 500-700MW of firm, dispatchable renewable electricity on O'ahu by 2033. In Phase 3, we also collected limited new encouraging groundwater data, and updated our thoughts on the probabilities of fluid and permeability at resource depths (Pr F = 1; Pr P = mostly unconstrained). Ultimately, we advocate for using our final probability of heat, and confidence in this probability, to drive the next phase of exploration. We contend further development of geothermal in Hawai‘i will enable the state to achieve its 100% renewable policy objective and Hawai'i to transition off of fossil fuels through geothermal discovery and development. The project not only produced a large amount of data and expanded the existing knowledge of Hawai'i's geothermal resources, but also produced publications, theses, presentations, core photos, datasets, media reports, television interviews, community events, and a blog. Students and new professionals benefitted from the project's hands-on research experiences and educational opportunities and earned awards and recognition.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Trends in electronic structures and s ± -wave pairing for the rare-earth series in bilayer nickelate superconductor R 3 Ni 2 O 7

The recent discovery of pressure-induced superconductivity in the bilayer La 3 Ni 2 O 7 (LNO) has opened a new platform for the study of unconventional superconductors. In this publication, we investigate theoretically the whole family of bilayer 327-type nickelates R 3 Ni 2 O 7 (R = rare-earth elements) under pressure. From La to Lu, the lattice constants and volume decrease, leading to enhanced in-plane and out-of-plane hoppings, resulting in an effectively reduced electronic correlation U / W. Furthermore, the Ni's t 2g states shift away from the e g states, while the crystal-field splitting between d 3z 2 -r 2 and d x 2 -y 2 is almost unchanged. In addition, six candidates were found to become stable in the Fmmm phase, with increasing values of critical pressure as the atomic number increases. Similar to the case of LNO, the s ± -wave pairing tendency dominates in all candidates, due to the nesting between the M = (π , π) and the X = (π , 0) and Y = (0 , π) points in the Brillouin zone. Then, T c is expected to decrease as the radius of rare-earth (RE) ions decreases. In conclusion, our results suggest that LNO is already the “optimal” candidate, with Ce a close competitor, among the whole of the RE bilayer nickelates, and to increase T c we suggest growing on special substrates with larger in-plane lattice spacings.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Effect of seasonal anoxia on geochemical cycling in a stratified pond: Comparison to cooler pond conditions 40 years ago

Seasonal stratification in temperate lakes deeper than a few meters creates favorable conditions for pronounced vertical redox zones, often resulting in anaerobic hypolimnions and significant geochemical changes. Here, this study examined thermocline formation and trace element behavior in a seasonally stratified pond amid rising air temperatures. Over two years, data were collected from Pond B at the US Department of Energy Savannah River Site in Aiken, South Carolina. Pond B, a man-made monomictic reservoir, received cooling water from a nuclear reactor from 1961 to 1964. Strong thermal stratification forms a distinct thermocline in May and progresses downward until November. Compared to the 1980s, this study shows a delayed onset and extended duration of stratification. The prolonged summer stratification reduces deep water oxygen replenishment, extending hypoxic conditions. Trace and major elements sampled in the water column revealed strong correlations between As, Fe, and Mn profiles, with concentrations increasing by 1–2 orders of magnitude in the anaerobic hypolimnion. This period captured the seasonal transition from winter mixing to summer stratification to fall overturn. Under anoxic conditions, Fe(III) reduces to Fe(II) in the sediment, releasing dissolved iron into the water column. The extended anoxic periods likely promoted arsenic release from sediments. Prolonged anoxia may enhance arsenic mobilization and solubility in the lake. This study illustrates how climate-induced changes in seasonal stratification of contaminated waters can convert contaminant sinks into sources, offering insights into the cycling of arsenic and other dissolved ions in stratified lakes and their implications for water quality management.

Anoxic conditions↗
DOE OSTI2021

Electronic Structure and Stacking Arrangement of Tungsten Disulfide at the Gold Contact

There is an intensive effort to control the nature of attractive interactions between ultrathin semiconductors and metals and to understand its impact on the electronic properties at the junction. In this work, we present a photoelectron spectroscopy study on the interface between WS 2 films and gold, with a focus on the occupied electronic states near the Brillouin zone center (i.e., the Γ point). To delineate the spectra of WS 2 supported on crystalline Au from the suspended WS 2 , we employ a microscopy approach and a tailored sample structure, in which the WS 2 /Au junction forms a semi-epitaxial relationship and is adjacent to suspended WS 2 regions. The photoelectron spectra, as a function of WS 2 thickness, display the expected splitting of the highest occupied states at the Γ point. In multilayer WS 2 , we discovered variations in the electronic states that spatially align with the crystalline grains of underlying Au. Corroborated by density functional theory calculations, we attribute the electronic structure variations to stacking variations within the WS 2 films. We propose that strong interactions exerted by Au grains cause slippage of the interfacing WS 2 layer with respect to the rest of the WS 2 film. Our findings illustrate that the electronic properties of transition metal dichalcogenides, and more generally 2D layered materials, are physically altered by the interactions with the interfacing materials, in addition to the electron screening and defects that have been widely considered.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Resonance from antiferromagnetic spin fluctuations for superconductivity in UTe 2

We report superconductivity originates from the formation of bound (Cooper) pairs of electrons that can move through the lattice without resistance below the superconducting transition temperature $T_c$. Electron Cooper pairs in most superconductors form anti-parallel spin singlets with total spin S = 0, although they can also form parallel spin-triplet Cooper pairs with S = 1 and an odd parity wavefunction. Spin-triplet pairing is important because it can host topological states and Majorana fermions relevant for quantum computation. Because spin-triplet pairing is usually mediated by ferromagnetic (FM) spin fluctuations3, uranium-based materials near an FM instability are considered to be ideal candidates for realizing spin-triplet superconductivity. Indeed, UTe 2 , which has a $T_c$ ≈ 1.6 K, has been identified as a candidate for a chiral spin-triplet topological superconductor near an FM instability, although it also has antiferromagnetic (AF) spin fluctuations. Here we use inelastic neutron scattering (INS) to show that superconductivity in UTe 2 is coupled to a sharp magnetic excitation, termed resonance, at the Brillouin zone boundary near AF order. Because the resonance has only been found in spin-singlet unconventional superconductors near an AF instability, its observation in UTe 2 suggests that AF spin fluctuations may also induce spin-triplet pairing or that electron pairing in UTe 2 has a spin-singlet component.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

PDPTW-DB: MILP-Based Offline Route Planning for PDPTW with Driver Breaks

The Pickup and Delivery Problem with Time Windows (PDPTW) involves optimizing routes for vehicles to meet pickup and delivery requests within specific time constraints, a challenge commonly faced in logistics and transportation. Microtransit, a flexible and demand-responsive service using smaller vehicles within defined zones, can be effectively modeled as a PDPTW. Yet, the need for driver breaks—a key human constraint—is frequently overlooked in PDPTW solutions, despite being necessary for regulatory compliance. This study presents a novel mixed-integer linear programming formulation for the Pickup and Delivery Problem with Time Windows and Driver Breaks (PDPTW-DB). To the best of our knowledge this formulation is the first to consider mandatory periodic driver breaks within optimized Microtransit routes. The proposed model incorporates regulatory compliant break scheduling directly within the vehicle routing optimization framework. By considering driver break requirements as an integral component of the optimization process, rather than as a post-processing step, the model enables the generation of routes that respect hours of service regulations while minimizing operational costs. This integrated approach facilitates the generation of schedules that are operationally efficient and prioritize driver welfare through driver breaks. We work with a public transit agency from the southern USA, and highlight the specific nuances of driver break optimization, and present a Pickup and Delivery Problem with Time Windows formulation for optimizing Microtransit operations and scheduling driver breaks. We validate our approach using real-world data from the transit agency. Our results validate our formulation in producing cost-effective, and regulation-compliant solutions.

Applied Computing, Transportation↗
DOE OSTI2021

Flowback Test Analyses at the Utah Frontier Observatory for Research in Geothermal Energy (FORGE) Site

Abstract In 2017 and 2019, injection testing was carried out in three zones in a vertical well in granite at the Frontier Observatory for Research in Geothermal Energy site near Milford, Utah, USA. In several injection cycles, flowback was implemented rather than shut-in. The goal was to explore an alternative to prolonged shut-in periods for inferring closure stress, formation compressibility, and formation permeability (permeability thickness product). The flowback procedures involved a cyclic flowback/shut-in, while pressure decreased. The flowback data are presented, and analyses are shown. The inferred closure stress(es) from flowback analyses are lower than for equivalent injection cycles that were strictly shut-in. Relatively high formation compressibility obtained from the flowback analyses indicates an extensive, fractured system. This study also includes numerical simulation of the flowback events. The numerical model shows that the rebound pressure is not necessarily the lower bound of the minimum principal stress. The signature of stiffness change can be identified as the process when the depletion mainly transitions from hydraulic fracture to natural fractures from numerical analysis. Overall, flowback potentially has advantages over shut-in because of the reduced time to closure.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Integrated parameter and process learning for hydrologic and biogeochemical modules in Earth System Models

Focus area: Primary focal area #2; secondary focal area #3: Learning about parameters and processes of land surface hydrologic and biogeochemical models in Earth System models by integrating machine learning, physics, and big data. Science challenges: How do we maximally leverage big-data observations to improve hydrobiogeochemical process description and parameterization so that such modules more realistically capture hydrologic and vegetation responses and feedbacks under the future climate? For example, how can we leverage physics, limited observations of vegetation and streamflow to better estimate evapotranspiration, and, relatedly, net primary productivity, especially for drought areas? Vegetation plays a critical role in regional and global water cycles; however, existing vegetation models have failed to predict vegetation response to droughts (McDowell & Xu, 2017) , arctic greening (Keenan & Riley, 2018) , and critical transitions between forest and savanna (Hirota et al., 2011) . These studies suggest that when we build process-based models (PBM) parameterized from regional and global plant traits, we tend to poorly describe plant adaptation and local-scale competition processes. The models and their associated parameters assigned for different regions in the world are not capturing essential heterogeneity in vegetation responses at finer spatial scales. Many parameters of the land surface models control hydrology and vegetation dynamics at the same time. The heterogeneity in vegetation response is a function of (i) plant type, (ii) plant size, (iii) competition and succession, (iv) environmental controls, and (v) local variations due to the unique ecological community that are very difficult to describe (e.g., the size of gaps resulting from fire that facilitated the coexistence of pioneering species). In the demographic models, only factors (i) and (iv) were captured, and plant types were generally described only by leaf phenology and climate zones. With current demographic models, we generally consider more traits to define plant types (i) and calibrate these traits to consider factors (ii), (iii) and (iv); however, it is substantially challenging to scale to regional and global simulations due to trait variations across space (Ali et al., 2016). Moreover, it has been noted that hillslope processes, including ridge-to-valley flow and sunny vs. shady slopes are primary organizers of water, energy, and vegetation (Clark et al., 2015; Fan et al., 2019) . Although gradual improvements in the hydrologic model component in earth system models may reduce this error (at a remarkably slow pace), the long-term, gradual impact of hydrology on plant traits are not well captured. Recent work showed that the hydrologic controls exerted by groundwater and lateral flow are primary regulators of rooting depth (Fan et al., 2017) . Such hydrologic controls have seldom been reflected in vegetation model parameterizations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Two-dimensional overdamped fluctuations of the soft perovskite lattice in CsPbBr 3

Lead halide perovskites exhibit structural instabilities and large atomic fluctuations thought to impact their optical and thermal properties, yet detailed structural and temporal correlations of their atomic motions remain poorly understood. Here, these correlations are resolved in CsPbBr 3 crystals using momentum-resolved neutron and X-ray scattering measurements as a function of temperature, complemented with first-principles simulations. We uncover a striking network of diffuse scattering rods, arising from the liquid-like damping of low-energy Br-dominated phonons, reproduced in our simulations of the anharmonic phonon self-energy. Here, these overdamped modes cover a continuum of wave vectors along the edges of the cubic Brillouin zone, corresponding to two-dimensional sheets of correlated rotations in real space, and could represent precursors to proposed two-dimensional polarons. Further, these motions directly impact the electronic gap edge states, linking soft anharmonic lattice dynamics and optoelectronic properties. These results provide insights into the highly unusual atomic dynamics of halide perovskites, relevant to further optimization of their optical and thermal properties.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A block-spectral adaptive H-/$p$-refinement strategy for shock-dominated problems

An adaptive H-/p-refinement strategy using a novel sensor is devised and tested in a block-spectral compressible Euler code equipped with adaptive-mesh refinement (AMR) and high-order flux-reconstruction numerics. At each Gauss quadrature point (or solution point) within each spectral block (or mesh element) the discrete velocity jump ΔU = ∂U/∂y 1 Δy 1 + ∂V/∂y 2 Δy 2 + ∂W/∂y 3 Δy 3 is calculated and normalized by the local speed of sound, a. Here, the grid spacing, Δx i , is calculated in each direction as the distance between auxiliary Gauss-Lobatto points, staggered relative to the solution points. The polynomial order is increased from p = 0 to p = p max in regions of weak compression, (ΔU/a) crit < ΔU/a < 0 and kept at p = p max in regions of flow expansion ΔU/a ≥ 0, while staying at the H = 0 base mesh level. Regions experiencing strong compressions, i.e. ΔU/a < (ΔU/a) crit , are H-refined up to H = H max where H max is applied at the location of maximum compression, ΔU/a = min(ΔU/a) in the domain, while keeping p = 0 to guarantee robustness and monotonicity of the solution in the H refined region. The critical value of (ΔU/a) crit = -0.06 is found to effectively separate smooth and non-smooth solution regions, supported by a 1D detonation initiation test case in ideal gas and a shock-to-detonation transition in high explosives. Using this value, the Sod shock tube, Shu-Osher problem, double Mach reflection and a 2D detonation in a high-explosive are simulated with the proposed adaptive H-/p-refinement. In the Sod shock tube case, p-refinement resolves the (weak) contact discontinuity while H-refinement enhances the grid resolution in the shock exploiting the monotonicity of the p = 0 reconstruction. For the Shu-Osher problem, p-refinement captures the small-scale oscillations trailing the shock that would be otherwise attenuated, while H-refinement triggered by the ΔU-sensor appropriately tracks the shock. In the double Mach reflection problem, H-refinement confines the numerical diffusion around the reflected shock while p-refinement recaptures many physical features trailing the shock. Finally, in the 2D high-explosive detonation case, H-refinement follows the leading shock and resolves the curvature of the detonation wave, while p-refinement adds resolution to the trailing reaction zone. Finally, the proposed methodology is tested in a detonation-wave propagation test case in high-explosives with numerical predictions comparing favorably against experiments.

97 MATHEMATICS AND COMPUTING↗