Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Observing System Simulation”

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,441 records · Page 80
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Validation of Thermohydraulic Simulations using RELAP for Critical Dual Purpose Canisters

Simulations using the RELAP5-3D systems analysis code and STAR-CCM+, a high-fidelity CFD code, were performed on a model of a real-scale MPC-32 dual purpose canister undergoing a criticality event. These simulations were performed to provide code-to-code verification of the natural convection capabilities in RELAP5. The STAR-CCM+ simulations revealed various natural convection flow patterns that aided in the cooling of fission heat produced in the fuel rods in the canister. Negligible spatial temperature variations were observed across the canister in the STAR-CCM+ results, supporting the use of RELAP5 in predicting the sub-cooled regime. The surface-averaged rod temperature predicted by RELAP5 compared well with STAR-CCM+, verifying the predictive capability of RELAP5 for sub-cooled conditions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Human Factors and Technologies Design to Improve User Acceptance of Pooled Rideshare for Increasing Transportation System Energy Efficiency

This multi-year project delivered a comprehensive, human-factors-driven framework to understand, model, and improve pooled rideshare (PR) adoption in the United States. Through three large-scale national survey studies involving more than 16,000 participants across multiple cities and demographic groups, the research established one of the most extensive datasets to date on user perceptions, behavioral barriers, and service expectations related to pooled rideshare. These data revealed key human factors barriers of user acceptance of PR and suggested potential actionable experience optimizations that could lead to increased PR usage. This foundational knowledge guided the development of novel human-factors models and behavioral choice models that quantify how psychological, demographic, and trip-level factors influence willingness to pool. Building on these empirical insights, the project developed advanced behavioral modeling tools, including mixed logit and integrated choice and latent variable models, to capture both observable and latent influences on PR adoption. These models significantly improved the ability to predict riders’ acceptance of pooled trips, explaining choice heterogeneity through latent constructs such as safety, service experience, privacy concerns, time sensitivity, and environmental attitudes. Together, these models provide a robust analytical foundation for designing PR systems that more effectively meet user needs. The project translated human-factors insights and behavioral models into actionable technology innovations by extending POLARIS—an agent-based, activity-based travel simulation platform—into a fully functional pooled rideshare simulation environment. New PR modules, acceptance models, and regional scenarios were implemented for Greenville, SC and Austin, TX, enabling high-fidelity validation of algorithmic strategies under realistic demand and traffic conditions. The simulation platform supported the development and evaluation of adaptive discount-based assignment algorithms, enhanced willingness-to-pay formulations, demographic-aware incentive mechanisms, and a proactive joint assignment and repositioning strategy. Simulation results demonstrated substantial gains in pooling uptake, average vehicle occupancy, energy efficiency, and fleet profitability. In Greenville, pooling adoption more than doubled, while reductions in vehicle-miles traveled and energy consumption were significant. In Austin, pooling improvements were achieved with minimal service-quality trade-offs, and profitability increased across all fleet sizes. Through this research, we developed a comprehensive understanding of the human factors barriers that limit user acceptance of pooled rideshare services. These insights enabled the design of human-factors-aware pooled rideshare technologies that more effectively address user concerns and improve adoption rates. By integrating these models into an advanced agent-based simulation framework, we demonstrated that higher adoption of pooled rideshare can lead to measurable improvements in energy efficiency and system performance. Together, these contributions establish a validated pathway from human-centered analysis to technology development and energy-saving outcomes, supporting national goals for more sustainable and efficient mobility systems.

Jia, Yunyi↗
DOE OSTI2024

Dual-pump Coherent anti-Stokes Raman scattering measurements in H 2 /CH 4 Counterflow Flames

Temperature and major species concentration measurements are reported in a Williams Seshadri counterflow burner using dual-pump coherent anti-Stokes Raman scattering (DPCARS). Hencken burner H2/Air flames with equivalence ratio of 0.8-1.4 were used to characterize the accuracy of the system. The mean errors for temperature, XH2/XN2 mole fraction ratio, and XO2/XN2 mole fraction ratio are 1.0%, 14.7%, and 5.5%, respectively. The counterflow burner uses a CH4, H2, and N2 mixture as fuel and synthetic air as oxidizer. Measurements were performed along the counterflow burner central axis for two different fuel compositions: 4% H2/86% CH4/10% N2 and 87%H2/3% CH4/10% N2. The measured temperature and species mole fractions agree well with the 1-D opposed-flow diffusion flame simulation. One observation is that the fuel inlet temperature is higher than 300 K, which suggests the possibility of the stagnation plane residing closer to the fuel side. Lastly, these measurements provide valuable data for testing and validation of advanced 2-D flame models.

Murdock, Benjamin K.↗
DOE OSTI2024

Simulation, Structural, Thermal and Mechanical Properties of the FeTiTaVW High Entropy Alloy

Developing new materials to be applied in extreme environments is an opportunity and a challenge for the future. High entropy alloys are new materials that seem promising approaches to work in nuclear fusion reactors. In this work, FeTaTiVW high entropy alloys were developed and characterized with Molecular Dynamic and Hybrid Molecular Dynamic Monte Carlo simulations. The simulation results show that phase separation originates a lower potential energy per atom and a high level of segregation compared to those of a uniform solid solution. Moreover, the experimental diffractogram of the milled powder shows the formation of a body-centred cubic-type structure and the presence of TiO2. In addition, the microstructure of the consolidated material evidenced three phases: W-rich, Ti-rich, and a phase with all the elements. This phase separation observed in the microstructure agrees with the Hybrid Molecular Dynamic Monte Carlo simulation. Moreover, the consolidated material’s thermal conductivity and specific heat are almost constant from 25 °C to 1000 °C, and linear expansion increases with increasing temperature. On the other hand, specific heat and thermal expansion values are in between CuCrZr and W values (materials chosen for the reactor walls). The FeTaTiVW high entropy alloy evidences a ductile behaviour at 1000 °C. Therefore, the promising thermal properties of this system can be attributed to the multiple phases and systems with different compositions of the same elements, which is exciting for future developments.

Martins, Ricardo (ORCID:0000000177685704)↗
DOE OSTI2020

Post-Calibration Uncertainty Analysis for Travel Times at a Naval Weapons Industrial Reserve Plant

The Naval Weapons Industrial Reserve Plant (NWIRP) in McGregor, Texas began manufacturing explosives in 1980 and several hazardous chemicals were discovered in lakes and streams surrounding the site in 1998. Contaminants traveled to local lakes and streams much faster than initially predicted. This research estimated contaminant travel times and identified locations where monitoring wells should be installed to yield the greatest reductions in uncertainties in travel-time predictions. To this end, groundwater and particle-tracking models for NWIRP site were built to predict hydraulic heads and contaminant travel times. Next, parameter (hydraulic conductivities) uncertainties, parameter identifiabilities, observation (hydraulic heads) worth, and predictive (contaminant travel times) uncertainties were quantified. Parameter uncertainties were reduced by up to 92%; a total of 19 of 158 parameters were at least moderately identifiable; travel-time uncertainties were reduced up to 92%. Additionally, travel-time predictions and post-calibration parameter distributions were generated using the null-space Monte Carlo (NSMC) technique. NSMC predicted that conservative tracers exited the flow system within a year, which matches with field data. Finally, an observations-worth analysis found that additional 11 more measurements would reduce travel-time uncertainties by factors from 1.04 to 4.3 over existing data if monitoring wells were installed at the suggested locations.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2026

Multi-model analysis of the impact of water vapor on the radiative forcing of volcanic aerosols after the 2022 Hunga Eruption

On 15 January 2022, the Hunga volcano eruption released unprecedented amounts of water vapor into the atmosphere alongside a modest amount of SO 2 . In this work we analyse results from multiple Earth system models as part of the Hunga Tonga-Hunga Ha'apai Volcano Impact Model Observation Comparison Project. Our results show a good model agreement over the climatic outcomes of the eruption, overall indicating a significant negative radiative forcing from the Hunga eruption. The multi-model mean of global instantaneous radiative forcing averaged over 2022–2023 is estimated at −0.19 ± 0.06 W m −2 at the top-of-atmosphere (TOA), and −0.16 ± 0.06 W m −2 at the surface. Simulations with free-running meteorology and climatological sea surface temperatures and sea ice yield a global mean TOA forcing of −0.14 ± 0.10 W m −2 across two models for the first 2 years, decreasing to −0.09 ± 0.10 W m −2 on average between 2022 and 2027. However, these global values may be underestimated by about 50 %, considering that recent SO 2 injection retrievals suggest nearly twice the amount than the 0.5 Tg-SO 2 used in the protocol. We also find that the contribution from added stratospheric water vapor is minimal and that the injected SO 2 and the resulting formation of stratospheric sulfate dominate the radiative forcing. However, water vapor played a key role in the initial aerosol growth, leading to a stronger negative radiative forcing during the first 6 months after the eruption compared to simulations without water vapor co-injection.

Geosciences↗
DOE OSTI2020

Sensitivity Study for Forecasting Variables of WRF-Solar Using a Tangent Linear Approach

Integrating solar generation in recent years has highlighted the need for improved accuracy in predicting solar power. Confidence in solar power forecasting can be achieved by designing an ensemble that provides reliable probabilistic information for solar radiation with reduced uncertainty and error. Ideally, ensemble members are created through the optimized perturbation of the initial conditions in numerical weather prediction (NWP) models. Tangent linear models are capable of efficiently investigating the sensitivity of solar radiation to model input parameters because they do not require individual perturbation of each variable. This sensitivity study using tangent linear models provide us the capability to identify the right variables to perturb in an ensemble prediction system. In this study, we developed tangent linear models for WRF-Solar modules that directly impact the computation of solar radiation and the simulation of cloud formation and dissipation including the Fast All-sky Model for Solar Applications (FARMS), the Noah land surface model (LSM), the Thompson microphysics, the Mello-Yamada-Nakanishi-Niino (MYNN) boundary layer parameterization, and the Deng scheme for a shallow-convection parameterization. A sensitivity analysis was conducted under various scenarios based on satellite observations and model simulations from the National Solar Radiation Data Base (NSRDB) and WRF-Solar, respectively. Critical forecasting variables that are highly sensitive to the forecasting of global horizontal irradiance (GHI), direct normal irradiance (DNI), cloud mixing ratio, cloud tendency, cloud fraction, and sensible and latent heat fluxes were determined using the relevant WRF-Solar module. This study will be used as a guidance on future research leading to high-quality probabilistic solar forecasting. In this presentation, we discuss the validation of tangent linear approach for WRF-Solar modules and illustrate how the sensitivity results are valuable in the improvement of probabilistic solar prediction.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Peering down the barrel with DESI DR2: 10 000+ inflows at $z$ < 0.6 reveal how galaxies accrete cold gas

Direct observational constraints on how galaxies acquire their gas remain remarkably limited, hindering our understanding of the baryon cycle. We present a search for down-the-barrel NaI D absorption towards 15.6 million galaxies at $z < 0.6$ in DESI Data Release 2. We use Bayesian evidence ratios to assess whether the absorption requires additional components tracing interstellar gas distinct from the systemic component of the galaxy. We construct a catalogue of 50 088 (27 420) galaxies with moderate (strong) evidence for down-the-barrel absorption. The inferred absorption components are broadly distributed in velocity, with approximately 50% at $v_{\rm flow} < -50$ km/s, 30% within 50 km/s of the systemic velocity and the remaining 20% at $v_{\rm flow} > 50$ km/s. We find strong evidence for a large population of low-velocity, infalling absorbers with velocities $\sim$20 km/s in edge-on galaxies, consistent with radial inflows predicted in simulations. The stronger correlation in early-type galaxies between inflow velocity and stellar velocity dispersion, compared to that with stellar mass, suggests that a portion of these inflows may be associated with accreting satellites. These results reveal the multiple pathways in which galaxies accrete gas at redshift $z < 0.6$ for the first time in a statistically significant sample.

Weng, S. [Marseille, Lab. Astrophys.]↗
DOE OSTI2023

Experimental and Theoretical Studies of the Surface Oxidation Process of Rare–Earth Tritellurides

Recent studies have established Van der Waals (vdW) layered and 2D rare-earth tritellurides (RTe 3 ) as superconductors and near room-temperature charge density wave (CDW) materials. Their environmental stability raises natural concern owing to aging/stability effects observed in other tellurium-based layered crystals. Here, the results establish the stability and environmental aging characteristics of these RTe 3 systems involving a variety of metals such as La, Nd, Sm, Gd, Dy, and Ho. The atomic force microscopy (AFM) and scanning electron microscopy (SEM) results show that all the RTe 3 sheets oxidize to form thin TeO x layers that are primarily confined to the surface, edges, and grain boundaries. Time-resolved in situ Raman spectroscopy measurements are used to understand the kinetics of the oxidization process for different lanthanide metal cations and establish their relative stability/resilience to oxidization. Overall results indicate that the vdW layers show higher air stability as the 4f electron number decreases going from Ho to La, resulting in the most stable LaTe 3 compared to the least stable HoTe 3 . Comprehensive quantum mechanical simulations reveal that environmental degradation originates from a strong oxidizing reaction with O 2 molecules, while humidity (H 2 O) plays a negligible role unless Te vacancies are present. Moreover, the simulations explain the effects of 4f electrons on the work function and Te vacancies formation, which directly impact the aging characteristics of RTe 3 layers. Interestingly, optical and electrical measurements show that the CDW response is still observed in aged RTe 3 layers owing to the presence of underlying pristine/nonoxidized RTe 3 layers, except CDW transition temperatures increase due to the thickness effect. Overall results offer the first in-depth environmental aging studies on these materials, which can be applied to engineer and design their chemical stability, surface properties, and overall CDW characteristics.

36 MATERIALS SCIENCE↗
DOE OSTI2022

LES of HCCI combustion of iso-octane/air in a flat-piston rapid compression machine

Homogeneous Charge Compression Ignition (HCCI) engines promise better efficiency and cleaner emissions than conventional piston engines, but can be challenging to control. Rapid compression machines (RCM) provide a simplified configuration for investigating HCCI combustion behavior, which is necessary for effective control of engine ignition timing and peak pressures. Here, in this study, we assess the utility of large eddy simulations (LES) for predicting HCCI combustion in a 3-D configuration. To this end, LES with finite-rate chemistry employing a 99-species iso-octane/air mechanism of two RCM operating conditions are performed. The RCM configuration under consideration was designed by Strozzi et al. with a flat piston to introduce large amounts of thermal stratification representative of realistic HCCI engine conditions, through the generation of corner vortices. It is shown that the simulation provides reasonable agreement with temperature fluctuations (7% difference), as well as ignition delay in the short ignition case (1 ms difference), while the long ignition case (35 ms difference) highlights more substantial deficiencies that are still within the expected uncertainty from the employed chemical mechanism. Flame propagation modes predicted by LES agree with experimental observations: spontaneous ignition is seen in the short ignition case, while deflagration is more prominent in the long ignition case. Analysis of global and local quantities classify the short ignition case in a mixed ignition regime, and the long ignition case in the mild ignition regime. These results demonstrate the utility of FRC-LES for investigations of multimode combustion regimes of HCCI combustion in a 3-D configuration.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

Contaminant Transport Modeling for Technology Evaluation and Long-Term Monitoring in the Tims Branch Testbed, SC - 20343

Studies conducted at U.S. DOE sites have shown the presence of heavy metals, radionuclides, and volatile organic compounds in surface water, groundwater, and soil as a result of nuclear activity in the Cold War era. Since the 1990's, innovative cleanup methods have been implemented in the Tims Branch watershed at Savannah River Site (SRS) to limit the contaminant flux to the stream that have reduced the contaminant concentrations to acceptable regulatory levels in the dissolved phase. A tin-based treatment which effectively eliminated all local anthropogenic mercury inputs to this ecosystem resulted in a known step function addition of inert tin oxide particles which now serve as a potential tracer for sedimentation and particle transport processes in the stream. The long-term effectiveness of this and other remediation techniques and the potential for remobilization of adsorbed contaminant in sediment during extreme hydrologic conditions however remains unclear. It is therefore important to understand not only the fate and transport of dissolved contaminants, but also the movement of sediment and the relevant interactions with dissolved contaminant. To narrow this knowledge gap, a study is being conducted using the Tims Branch watershed as a stream-scale ecosystem test-bed to identify the primary transport processes of major contaminants of concern (such as mercury, nickel and uranium) with an emphasis on interactions with sediment transport. This involves the development of a fully distributed hydrologic watershed model of the Tims Branch watershed to predict streamflow under extreme weather conditions, as well as the development of a comprehensive contaminant transport model that can properly account for coupled contaminant and sediment transport. Review of relevant research reports and peer-reviewed journals revealed that aside from advection-dispersion transport of dissolved contaminants, adsorption and desorption with suspended solids and bed sediment also play an important role in the transport of those contaminants of concern. To develop the fully distributed hydrologic watershed model, the MIKE SHE 2-dimensional (2D) land surface/3D groundwater model that simulates surface/subsurface hydrologic processes (such as overland flow, evapotranspiration, and infiltration) was coupled with a 1D streamflow model that accounts for stream water hydraulics (such as hydraulic structures, cross-sections, and network). To model contaminant transport, the MIKE 11 streamflow component was coupled with the MIKE 11 AD module that simulates solute transport through advection and dispersion, and MIKE ECO Lab module that accounts for both sediment transport and interactions with dissolved contaminant. At this stage, the development and optimization of the fully distributed hydrologic model has been completed, achieving satisfactory statistical results between observed and predicted discharge as indicated by a root mean square error (RMSE) of 0.039 cms and a Nash-Sutcliffe efficiency coefficient (NSE) of 0.764. The ongoing development of the contaminant transport model has also yielded realistic results from preliminary tests. Results from this study are a key to evaluating the effectiveness of tin (II)-based mercury treatment of wetlands at the SRS site, and are also relevant to evaluating the potential of using this type of novel remediation technology in other mercury-contaminated stream systems. Knowledge acquired from this research will also support interpretation of historical data on the trends of contaminant concentration distribution in Tims Branch, particularly considering the effect of extreme hydrological events on the stream flow and pollutant transport. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Connecting In Situ Stress and Wellbore Deviation to Near-Well Fracture Complexity Using Phase-Field Simulations

The interactions among in situ stress, rock fabric, wellbore geometry, natural fractures, and other natural or man-made defects create highly complex fracture trajectories in the near-wellbore region, far more intricate than those in the far-field. These near-wellbore complexities are critical for the Utah FORGE project and Enhanced Geothermal Systems (EGS) in general. Frictional pressure loss in the near-wellbore region during stimulation can significantly influence the growth of far-field fractures, while pressure losses during circulation serve as a major source of energy dissipation. Near-wellbore fracture complexities are often observable through image logs, offering valuable insights into in situ stress characteristics. However, leveraging this information requires a high-fidelity model capable of capturing the interplay among the diverse factors influencing fracture behavior.

58 GEOSCIENCES↗
DOE OSTI2023

Influence of Pore Surface Structure and Contents on Shock-Induced Collapse and Energy Localization

The majority of computational efforts in the shockwave driven collapse of porosity focus on idealized pore geometries and structures. These give significant insight into the physical mechanisms in play, but omit some potentially important effects of real energetic materials. Here, to address this, we simulate pore collapse in 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) for four types of pore surfaces and fillings: an idealized cylindrical pore, a pore surface coated with polymer, a melted/amorphous pore surface, and a pore filled with gas molecules. We find that these perturbations to the pore structure result in only minor changes to the highly complex collapse mechanisms. However, significant differences in hotspot temperature are observed, with peak temperatures ranging over several hundred Kelvin. The polymer-coated system has the most effect, lowering the hotspot temperature by 300 K, and also limiting the manifestation of intramolecular strain energy that is known to induce mechanochemistry.

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

Interface microstructure effects on dynamic failure behavior of layered Cu/Ta microstructures

Abstract Structural metallic materials with interfaces of immiscible materials provide opportunities to design and tailor the microstructures for desired mechanical behavior. Metallic microstructures with plasticity contributors of the FCC and BCC phases show significant promise for damage-tolerant applications due to their enhanced strengths and thermal stability. A fundamental understanding of the dynamic failure behavior is needed to design and tailor these microstructures with desired mechanical responses under extreme environments. This study uses molecular dynamics (MD) simulations to characterize plasticity contributors for various interface microstructures and the damage evolution behavior of FCC/BCC laminate microstructures. This study uses six model Cu/Ta interface systems with different orientation relationships that are as- created, and pre-deformed to understand the modifications in the plasticity contributions and the void nucleation/evolution behavior. The results suggest that pre-existing misfit dislocations and loading orientations (perpendicular to and parallel to the interface) affect the activation of primary and secondary slip systems. The dynamic strengths are observed to correlate with the energy of the interfaces, with the strengths being highest for low-energy interfaces and lowest for high-energy interfaces. However, the presence of pre-deformation of these interface microstructures affects not only the dynamic strength of the microstructures but also the correlation with interface energy.

42 ENGINEERING↗
DOE OSTI2023

Projecting Future Energy Production from Operating Wind Farms in North America. Part I: Dynamical Downscaling

Abstract New simulations at 12-km grid spacing with the Weather and Research Forecasting (WRF) Model nested in the MPI Earth System Model (ESM) are used to quantify possible changes in wind power generation potential as a result of global warming. Annual capacity factors (CF; measures of electrical power production) computed by applying a power curve to hourly wind speeds at wind turbine hub height from this simulation are also used to illustrate the pitfalls in seeking to infer changes in wind power generation directly from low-spatial-resolution and time-averaged ESM output. WRF-derived CF are evaluated using observed daily CF from operating wind farms. The spatial correlation coefficient between modeled and observed mean CF is 0.65, and the root-mean-square error is 5.4 percentage points. Output from the MPI-WRF Model chain also captures some of the seasonal variability and the probability distribution of daily CF at operating wind farms. Projections of mean annual CF (CF A ) indicate no change to 2050 in the southern Great Plains and Northeast. Interannual variability of CF A increases in the Midwest, and CF A declines by up to 2 percentage points in the northern Great Plains. The probability of wind droughts (extended periods with anomalously low production) and wind bonus periods (high production) remains unchanged over most of the eastern United States. The probability of wind bonus periods exhibits some evidence of higher values over the Midwest in the 2040s, whereas the converse is true over the northern Great Plains. Significance Statement Wind energy is playing an increasingly important role in low-carbon-emission electricity generation. It is a “weather dependent” renewable energy source, and thus changes in the global atmosphere may cause changes in regional wind power production (PP) potential. We use PP data from operating wind farms to demonstrate that regional simulations exhibit skill in capturing actual power production. Projections to the middle of this century indicate that over most of North America east of the Rocky Mountains annual expected PP is largely unchanged, as is the probability of extended periods of anomalously high or low production. Any small declines in annual PP are of much smaller magnitude than changes due to technological innovation over the last two decades.

Meteorology & Atmospheric Sciences↗
DOE OSTI2023

Biogeophysical Effects of Land-Use and Land-Cover Change Not Detectable in Warmest Month

Land-use and land-cover changes (hereafter simply “land use”) alter climates biogeophysically by affecting surface fluxes of energy and water. Yet, near-surface temperature responses to land use across observational versus model-based studies and spatial-temporal scales can be inconsistent. Here we assess the prevalence of the historical land use signal of daily maximum temperatures averaged over the warmest month of the year (t LU ) using regularized optimal fingerprinting for detection and attribution. We use observations from the Climatic Research Unit and Berkeley Earth alongside historical simulations with and without land use from phase 6 of the Coupled Model Intercomparison Project to reconstruct an experiment representing the effects of land use on climate. To assess the signal of land use at spatially resolved continental and global scales, we aggregate all input data across reference regions and continents, respectively. At both scales, land use does not comprise a significantly detectable set of forcings for two of four Earth system models and their multimodel mean. Furthermore, using a principal component analysis, we find that t LU is mostly composed of the nonlocal effects of land use rather than its local effects. These findings show that, at scales relevant for climate attribution, uncertainties in Earth system model representations of land use are too high relative to the effects of internal variability to confidently assess land use.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Zinc oxide (ZnO) hybrid metasurfaces exhibiting broadly tunable topological properties

Extreme light confinement observed in periodic photonic structures, such as the vortex singularities in momentum (k) space, has been associated with their topological nature. Consequently, by exploiting and tuning their topological properties, optical metasurfaces have been demonstrated as an attractive platform for active photonics. However, given the fact that most active media under external excitations can only provide limited refractive index change, the potential advancements offered by the topological character of active metasurfaces have remained mostly unexplored. Zinc oxide (ZnO), which has recently exhibited optically-induced extraordinarily large permittivity modulations at visible and near-infrared frequencies, is an excellent active material for dynamic metasurfaces exhibiting strong tuning. This work demonstrates that a hybrid metasurface consisting of an array of ZnO nanodisks on a silver backplane displays broadly tunable topological properties. In particular, by performing k-space scattering simulations using measured pump-fluence-dependent material properties of ZnO, we study in detail the light reflection from the hybrid metasurface. Our results validate that the large k-space topology tuning of the metasurface can result in enormously strong polarization manipulation of near-infrared light in the vicinity of the topological features. The observed polarization switching effect is highly sensitive to the polarization and wavelength of an incident wave, owing to the symmetry and dispersion characteristics of the proposed system. Our study indicates that leveraging a combination of the extraordinary material properties and the k-space topology, hybrid metasurfaces based on ZnO may open new avenues for creating all-optical switchable metadevices.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Radiation Belt Remediation: A Complex Engineered System (RBR-ACES) (Institutional Computing Report)

We have performed the simulation campaigns in support of LDRD-DR project Radiation Belt Remediation: A Complex Engineered System (RBR-ACES). This project involves modeling wave generation and wave propagation in plasma environment and studying the effect of waves on the distribution of high energy electrons with the focus on electrons appeared after a high-altitude nuclear explosion (HANE). One of the major goal of this research is strong validation of the numerical models with laboratory experiments and space-based experiments and observations.

42 ENGINEERING↗