Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coupled 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,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Reducing warpage in a hybrid large-scale additive manufacturing and compression molding process

In recent years, a hybrid manufacturing process, developed by combining extrusion-based large-scale additive manufacturing (AM) and compression molding (CM) techniques, has shown promising outcomes for producing structurally functional parts. The process can be used with both short fiber-reinforced composites and neat polymers and hence, even multi-material parts can be manufactured easily. This process offers the advantages of structural enhancement by having a desired fiber orientation using a large-scale AM process, as well as rapid manufacturing capability using a CM process. In the large-scale AM process, the alignment of fibers in the deposition direction enables significant improvement in the mechanical properties of the manufactured parts. However, the anisotropy resulting from the directional arrangement of fibers also introduces challenges related to warpage in the produced parts. This study aims to identify the causes of warpage and propose strategies to mitigate it. The research involves the use of preforms manufactured through the large-scale AM, which are then combined with the neat resin for CM manufacturing. A finite element-based numerical simulation model is developed, employing a sequentially coupled thermomechanical approach. Through a parametric study using the simulation models, optimization of printing direction and preform geometry is performed to minimize warpage. This contributes to the advancement and wider adoption of AM/CM hybrid manufacturing to produce structurally functional parts.

Jo, Eonyeon↗
DOE OSTI2020

Risk Importance Ranking of Fire Data Parameters to Enhance Fire PRA Model Realism

Fire is historically and analytically a significant contributor to nuclear power plant risk. The level of fire risk and the methods, tools and data for modeling this risk is highly debated by experts. One area of debate is the input data used in fire modeling and how to deal with this data’s high uncertainty. This report outlines initial work performed for determining the key parameters causing this uncertainty and how it propagates into nuclear power plant models. This research paves the way for the development of methods to reduce fire data uncertainty used in modeling. The Nuclear Regulatory Commission has mandated that nuclear power plants perform fire risk modeling. However, there are several issues with the current risk modeling implementation that affect the results. Approved modeling methods can be overly conservative and often do not match plant experience. Also, the data used in the modeling can have high uncertainties and is influenced by expert judgement. To evaluate input data uncertainty, researchers performed an initial review of several fire experiments done at Sandia National Laboratories. Uncertainties for fire data can come from many sources, such as experiment design constraints, environmental conditions, or other plant-specific aspects. There are many different significant and insignificant parameters driving the uncertainty. Additionally, the uncertainty of the different input data used in the fire modeling could have a significant or insignificant effect on the entire plant risk. A four-step methodology was developed to perform Integrated Probabilistic Risk Assessment Importance Ranking. A demonstration case using these steps was set up and three of the four steps were completed in fiscal year (FY) 2019 and the fourth step done FY 2020. These steps are: 1. The qualitative analysis of potential sources was conducted with the following items identified for the demonstration. • Maximum heat release rate • Time to maximum heat release rate • Duration of max heat release rate • Time to decay • Thermal conductivity of concrete • Specific heat of concrete • Density of concrete • Cable jacket thickness 2. A quantitative characterization of dominant sources of uncertainty was performed. A list of distributions and determined values of the dominant sources is shown in Appendix A. 3. A quantitative screening of the potential sources of uncertainty using Morris Elementary Effects Analysis was completed. An experimental model using the physics-based fire modeling tool Fire Dynamics Simulator was developed and coupled with the Risk Analysis Virtual Environment. The Morris analysis identified at least two parameters that can be eliminated as significant contributors (specific heat of concrete and cable jacket thickness). 4. Global importance measure (Global IM) analysis to generate a comprehensive ranking based on their influence on the plant risk. In this research, a moment-independent Global IM is used since it can address (a) uncertainty in the input parameters of the fire model, (b) uncertainty in the risk outputs, and (c) non-linearity and interactions among input parameters in the fire model, more accurately than the correlation-based and variance-based global methods. The observations from the research showed that, depending on the initial and boundary conditions of the fire scenarios, fire-induced damage could have a very small probability and could be dominated by the tail of the uncertainty distribution; hence, the accuracy of the correlation-based and variance-based methods is questionable. The moment-independent Global IM analysis in this research provides a better understanding of how experimental uncertainty data affects industry’s plant models and where improvements in that data will have the largest benefit for improving fire modeling accuracy in causing core damage. Among the five unscreened parameters obtained from the Morris EE analysis, the Global IM analysis results for the case study indicated that max heat release rate and fire location are the most important parameters. The report also outlines benefits of using a unified computational platform that integrates the underlying simulations (e.g., a fire progression model), quantitative screening (using the Morris EE method), and the Global IM analysis. A unified platform can (i) facilitate the ranking of input parameters considering multiple key fire scenarios simultaneously, rather than considering one scenario at a time, (ii) contribute to more explicit and accurate treatment of dependencies at multiple levels of Fire PRA, (iii) facilitate the sampling-based uncertainty quantification for Fire PRA, and (iv) help generating both “industry-wide” and "plant-specific" ranking of uncertainty sources in Fire PRA. Future research should be done to include additional parameters such as detection/suppression or cable fire spread. Adding a PRA software such as SAPHIRE to the RAVEN platform would help with plant model integration and improve treatment of fire-induced dependency. The I-PRA risk importance ranking methodology offered in this report can provide valuable information for efficiently (a) enhancing the realism of Fire PRA for existing plants and (b) supporting the development of Dynamic Fire PRA for advanced reactors and new plants.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Phase-field modeling of crack growth under coupled creep-fatigue

Crack growth under coupled creep-fatigue strongly influences service life of metallic components operating at high temperatures. However, there lack numerical models that enable direct simulations of crack growth and plasticity development under coupled creep-fatigue. Here, in this study, a unified phase-field model is developed for simulating crack growth under fatigue, creep, and coupled creep-fatigue. This model is able to reproduce the Paris laws for both brittle and ductile materials in pure fatigue mode, which shows that the degradation rate due to cyclic stress in ductile materials is generally larger than that in brittle materials. Fatigue simulations with different yield strengths demonstrate that plasticity facilitates crack growth in typical ductile materials. High-throughput simulations are performed under coupled creep-fatigue conditions, with varying hold-stress levels and hold-time. By analyzing crack growth rate under different hold-time, a creep-fatigue interaction term is obtained. Simulation results demonstrate that the creep-fatigue interaction is caused by smaller stress gradient and enhanced degradation rate ahead of the crack tips under coupled creep-fatigue, compared to the situation of pure fatigue. This work reveals the origin of different Paris-law exponents in brittle and ductile materials and demonstrates how the interplay between creep and fatigue affects crack growth.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Free-spray characteristics and spray-wall interactions of methanol on a gasoline direct injector under flash-boiling and non-flash-boiling conditions

Methanol is considered a promising alternative fuel for internal combustion engines (ICEs) due to its high-octane number, fast laminar flame speed, and elevated latent heat of vaporization, all of which support higher compression ratios and improved thermal efficiency. However, its substantial latent heat of vaporization also poses cold-start challenges, such as misfire and fuel film deposition. This study aims to investigate methanol spray morphology and spray-wall interaction using the Spray M injector from the Engine Combustion Network within a constant-pressure flow vessel. A recently developed unified numerical framework capable of modeling both flash and non-flash boiling sprays is validated against experimental liquid volume fraction data acquired via 3-D computed tomography. Here, the results reveal that flash boiling significantly alters the spray morphology, leading to smaller droplets and spray collapse due to enhanced air-entrainment-induced turbulence. Quantitative agreement between experiments and simulations confirms this behavior. Coupled 0-D equilibrium and 3-D computational fluid dynamics analyses show that flash boiling accelerates evaporation and reduces fuel residence time, while non-flash conditions maintain a persistent liquid core more susceptible to wall wetting. Wall temperature diagnostics reveal that spray collapse alters heat transfer patterns by shifting cooling effects. Mixture fraction analysis indicates that evaporation is primarily governed by shear-layer turbulence, though deviations from adiabatic equilibrium mixing emerge under low-turbulence conditions. Finally, increasing fuel, ambient, and wall temperatures reduces wall wetting and film thickness, mitigating cold-start risks. These findings enhance the understanding of methanol sprays’ behavior and support its adoption as a viable, alternative fuel for ICEs.

Engine Combustion Network↗
DOE OSTI2021

Reduced-order modeling of near-field THMC coupled processes for nuclear waste repositories in shale

Performance assessment (PA) of geologic radioactive waste repositories requires three-dimensional simulation of highly nonlinear, thermo-hydro-mechanical-chemical (THMC), multiphase flow and transport processes across many kilometers and over tens to hundreds of thousands of years. Integrating the effects of a near-field geomechanical process (i.e. buffer swelling) into coupled THC simulations through reduced-order modeling, rather than through fully coupled geomechanics, can reduce the dimensionality of the problem and improve computational efficiency. In this study, PFLOTRAN simulations model a single waste package in a shale host rock repository, where re-saturation of a bentonite buffer causes the buffer to swell and exert stress on a highly fractured disturbed rock zone (DRZ). Three types of stress-dependent permeability functions (exponential, modified cubic, and Two-part Hooke’s law models) are implemented to describe mechanical characteristics of the system. Overall, our modeling study suggests that compressing fractures reduces DRZ permeability, which could influence the rate of radionuclide transport and exchange with corrosive species in host rock groundwater that could accelerate waste package degradation. Less permeable shale host rock delays buffer swelling, consequently retarding DRZ permeability reduction as well as chemical transport within the barrier system.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1999

Large Eddy Simulation of Cirrus Clouds

The Regional Atmospheric Modeling System (RAMS) with mesoscale interactive nested-grids and a Large-Eddy Simulation (LES) version of RAMS, coupled to two-moment microphysics and a new two-stream radiative code were used to investigate the dynamic, microphysical, and radiative aspects of the November 26, 1991 cirrus event. Wu (1998) describes the results of that research in full detail and is enclosed as Appendix 1. The mesoscale nested grid simulation successfully reproduced the large scale circulation as compared to the Mesoscale Analysis and Prediction System's (MAPS) analyses and other observations. Three cloud bands which match nicely to the three cloud lines identified in an observational study (Mace et al., 1995) are predicted on Grid #2 of the nested grids, even though the mesoscale simulation predicts a larger west-east cloud width than what was observed. Large-eddy simulations (LES) were performed to study the dynamical, microphysical, and radiative processes in the 26 November 1991 FIRE 11 cirrus event. The LES model is based on the RAMS version 3b developed at Colorado State University. It includes a new radiation scheme developed by Harrington (1997) and a new subgrid scale model developed by Kosovic (1996). The LES model simulated a single cloud layer for Case 1 and a two-layer cloud structure for Case 2. The simulations demonstrated that latent heat release can play a significant role in the formation and development of cirrus clouds. For the thin cirrus in Case 1, the latent heat release was insufficient for the cirrus clouds to become positively buoyant. However, in some special cases such as Case 2, positively buoyant cells can be embedded within the cirrus layers. These cells were so active that the rising updraft induced its own pressure perturbations that affected the cloud evolution. Vertical profiles of the total radiative and latent heating rates indicated that for well developed, deep, and active cirrus clouds, radiative cooling and latent heating could be comparable in magnitude in the cloudy layer. This implies that latent heating cannot be neglected in the construction of a cirrus cloud model. The probability density function (PDF) of w was analyzed to assist in the parameterization of cloud-scale velocities in large-scale models. For the more radiatively-driven, thin cirrus case, the PDFs are approximately Gaussian. However, in the interior of the deep, convectively unstable case, the PDFs of w are multi-modal and very broad, indicating that parameterizing cloud-scale motions for such clouds can be very challenging. The results of this research are described in detail in a paper submitted to the Journal of Atmospheric Science (Wu and Cotton, 1999), which is enclosed as Appendix 2. Using soundings extracted from a mesoscale simulation of the November 26, 1991 cirrus event, the radiative effects on vapor deposition/sublimation of ice crystals was studied using a two-dimensional cloud-resolving model (CRM) version of RAMS, coupled to an explicit bin-resolving microphysics. The CRM simulations of the November 26, 1991 cirrus event demonstrate that the radiative impact on the diffusional growth (or sublimation) of ice crystals is significant. In this case, the ice particles experienced radiative warming. Model results show that radiative feedbacks in the diffusional growth of ice particles can be very complex. Radiative warming of an ice particle will restrict the particle's diffusional growth. In the case of radiative warming, ice particles larger than a certain size will experience so much radiative warming that surface ice saturation vapor pressures become large enough to cause sublimation of the larger crystals, while smaller crystals are growing by vapor deposition. However, ice mass production can be enhanced in the case of radiative cooling of an ice particle. For the November 26, 1991 cirrus event, radiative feedback results in significant reduction in the total ice mass, especially in the production of large ice crystals, and consequently, both radiative and dynamic properties of the cirrus cloud are significantly affected. A complete description of this research has been submitted as a paper to the Journal of Atmospheric Science (Wu et al., 1999), and included as Appendix 3.

Wu, Ting↗
NASA NTRS2018

The SMAP Level-4 ECO Project: Improving Terrestrial Flux Estimates Through Coupled Hydrology-Vegetation Data Assimilation

Simulations of hydrologic and vegetation states as well as water, energy and carbon fluxes from the land surface to the atmosphere are crucial for a wide range of applications, including agricultural advisories, forecasts of (short-term) atmospheric behavior and seasonal weather predictions including forecasts of extreme events, such as heatwaves or droughts. The NASA Soil Moisture Active Passive (SMAP) mission Level-4 Eco-Hydrology (L4-ECO) project aims to improve modeled estimates of the terrestrial water, energy and carbon fluxes and states by developing a fully-coupled hydrology-vegetation data assimilation system. This system is developed around the NASA Goddard Earth Observing System (GEOS) Catchment-CN land surface model, which combines land hydrology and energy balance components of the GEOS Catchment model with dynamic vegetation components of the Community Land Model version 4. Catchment-CN fully couples the terrestrial water, energy and carbon cycles, allowing feedbacks from the land hydrology to the biosphere and vice versa.Here, we implement a calibration of the Catchment-CN vegetation parameterization against observations of the fraction of absorbed photosynthetically active radiation (FPAR) from the Moderate Resolution Imaging Spectroradiometer (MODIS) to improve the model's standalone skill. Later, the DA algorithm used to produce the SMAP L4 soil moisture product will be adapted to Catchment-CN to assimilate SMAP brightness temperatures and inform the model's land hydrology component. Finally, the DA system will be further extended to assimilate MODIS FPAR observations in order to constrain the model's dynamic vegetation component.In this presentation, we demonstrate that the Catchment-CN parameter calibration leads to more realistic vegetation simulations and reduces the root mean squared error between modeled and observed vegetation states across the model's various plant functional types. We also show that the assimilation of SMAP observations is able to improve the average correlation, bias and unbiased RMSE between the modeled surface and root zone soil moisture estimates, and ground observations from the SMAP core validation sites.

Kolassa, Jana↗
DOE OSTI2021

Multi-Scale Land-Atmosphere Interactions: Modeling Convective Processes from Plants to Planet

Research accomplishments include: 1)Two case studies of the effects of heterogeneous soil moisture and surface energy budgets on organization and propagation of convective precipitation during MC3E. 2) Development and evaluation of a new approach for simulating the effect of heterogeneous soil moisture at ARM-SGP using an innovative modeling approach. 3) Investigation of the effects of spatial coupling scale using multidecade global simulations in CESM with the multiscale modeling framework. 4) Exploration of changes to future precipitation intensity resulting from two different climate change scenarios using the multiscale model. 5) Provision of the new cloud-scale coupled multiscale Earth System Model to the larger community through the CESM process.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Bioinspired mechanically interlocking holey graphene@SiO 2 anode

Mechanically interlocking structures that can enhance adhesion at the interface and regulate the stress distribution have been widely observed in biological systems. Inspired by the biological structures in the wings of beetles, we synthesized a holey graphene@SiO 2 anode with strong mechanical interlocking, characterized it electrochemically, and explained its performance by finite element analysis and density functional calculations. The mechanically interlocking structure enhances lithium-ion (Li + ) storage by transmitting the strain from SiO 2 to the holey graphene and by a mechano-electrochemical coupling effect. The interlocking fit hinders the abscission of SiO 2 and the distinctive structure reduces the stress and strain of SiO 2 during (de)lithiation. The positive mechano-electrochemical coupling effect preserves the amount of electrochemically active phase (Li x Si) during cycles and facilitates Li + diffusion. Therefore, the capacity shows only a slight attenuation after 8000 cycles (cycling stability), and the specific capacity is ~1200 mA h g –1 at 5 A/g (rate-performance). This study furnishes a novel way to design high-performance Li + /Na + /K + /Al 3+ anodes with large volume expansion.

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

An ontology to represent synthetic building occupant characteristics and behavior

Since the introduction of the occupant behavior Drivers-Needs-Actions-Systems (DNAS) framework in 2013, researchers have used the framework or further developed it based on their case studies, which include efforts to collect new data on occupant behaviors. The effort is often costly for the relatively few new data points added. Problems emerge when the already collected data do not meet the modelers' interoperability requirements. Previous studies addressed this issue by developing more sophisticated ontologies that enable integration with other datasets and synthetic data methodologies that would meet unique research applications. This paper presents an extension of the DNAS framework for the representation of synthetic occupant data to support various applications and use cases across the building life cycle. An agent-based modeling application is one of our motivations that requires more elaborate characteristics of an occupant-agent or a group-of-agent. The extension, built upon a review of the literature, introduces new elements to the framework that fall into five categories, including socio-economic, geographical location, activities, subjective values, and individual and collective adaptive actions. On-going research includes identifying occupant datasets and developing data fusion methods to generate synthetic occupants, as well as to demonstrate its applications in agent-based modeling coupled with building performance simulation.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2022

Potassium carbonate decomposition modeling within rotating detonation engines for direct power extraction applications

We report that recently, there has been a significant interest in detonation-based combustion systems, such as rotating detonation engines (RDE), due to potential performance advantages in propulsion and energy applications. This includes improvements in thermodynamic efficiency, inherently high gas velocities, and the possibility of achieving an increase in total pressure (i.e. pressure gain) through the combustion process. Due to approximating a constant volume combustion process, these devices also tend to exhibit extremely high local gas temperatures relative to a comparable constant pressure combustor. These advantages overlap well with desirable performance characteristics of direct power extraction (DPE) technologies, such as a magnetohydrodynamic (MHD) generator. Typically, in DPE systems hot combustion products are seeded with an easily ionizable material such as potassium carbonate (K 2 CO 3 ) in order to boost the electrical conductivity. However, due to the short gas residence times within an RDE, it was unclear whether forming an electrically conductive combustion plasma would be feasible for integration with a downstream MHD generator. A model is presented which describes the heating, decomposition, and ionization of solid K 2 CO 3 particles and aqueous solutions of K 2 CO 3 in water, for a given initial particle diameter. This model was combined with available computational fluid dynamics (CFD) data for an oxygen-methane RDE in a one-way coupled Eulerian-Lagrangian framework to predict particle trajectories and the corresponding heating, decomposition, and ionization histories. Electrical conductivities were computed using a previously developed model, and a method was proposed to determine an equivalent average electrical conductivity. Results show that particle sizes below ~30 µm are able to fully decompose before reaching the exit of the RDE. While the one-way coupled nature of the simulations precluded rigorous evaluation of the effects of seed material on detonability, a substantial temperature reduction is expected at the detonation wave front due to heating and decomposition. A preliminary comparison is presented between the RDE and an equivalent constant pressure adiabatic combustor, showing a potential performance advantage for the RDE.

42 ENGINEERING↗
DOE OSTI2020

Assessment of reverse gun taylor cylinder experimental configuration

Experimental efforts for Taylor-anvil impact tests have often been limited to near room temperature. The ‘Reverse Gun’ method proposed by Gust in 1982 allows for the Taylor impact specimen to be uniformly heated without temperature losses before impact. Through the use of finite element analysis, we explore two topics in this work. First, we examine whether the reverse gun experimental configuration is comparable to the traditional Taylor-anvil setup. Second, we assess the accuracy of several commonly employed flow strength models in terms of their ability to predict the reverse gun experimental results which involve dynamic loading conditions and complex thermo-mechanical coupling. The reverse gun simulations are performed for tantalum targets at initial temperatures in the range 295 K to 1295 K and velocities from 135 m/s to 242 m/s. We show that with suitable care in the modeling of the preheated reverse gun experiments one can make valuable assessments of flow strength models. Given the conditions explored, these observations probe the thermal softening, strain hardening, and strain rate sensitivity of the material.

42 ENGINEERING↗
DOE OSTI2020

Integrated hydrogeophysical modelling and data assimilation for geoelectrical leak detection

Time-lapse electrical resistivity tomography (ERT) measurements provide indirect observations of hydrological processes in the Earth's shallow subsurface at high spatial and temporal resolution. ERT has been used in the past decades to detect leaks and monitor the evolution of associated contaminant plumes. Specifically, inverted resistivity images allow visualization of the dynamic changes in the structure of the plume. However, existing methods do not allow the direct estimation of leak parameters (e.g. leak rate, location, etc.) and their uncertainties. We propose an ensemble-based data assimilation framework that evaluates proposed hydrological models against observed time-lapse ERT measurements without directly inverting for the resistivities. Each proposed hydrological model is run through the parallel coupled hydro-geophysical simulation code PFLOTRAN-E4D to obtain simulated ERT measurements. The ensemble of model proposals is then updated using an iterative ensemble smoother. In this paper, we demonstrate the proposed framework on synthetic and field ERT data from controlled tracer injection experiments. Our results show that the approach allows joint identification of contaminant source location, initial release time, and solute loading from the cross-borehole time-lapse ERT data, alongside with an assessment of uncertainties in these estimates. We demonstrate a reduction in site-wide uncertainty by comparing the prior and posterior plume mass discharges at a selected image plane. This framework is particularly attractive to sites that have previously undergone extensive geological investigation (e.g., nuclear sites). It is well suited to complement ERT imaging and we discuss practical issues in its application to field problems.

58 GEOSCIENCES↗
DOE OSTI2025

A statistical approach to screening isotopic signatures in monitoring for underground nuclear explosions

The ability to differentiate between atmospheric radionuclide signatures from underground nuclear explosions (UNEs) and signals from other sources, such as medical isotope-production facilities and nuclear reactors, can be critical to the detection and monitoring of unannounced, low-yield nuclear events. Signatures having anomalously high amplitudes, compared to background levels, remain the best indicator in screening for a UNE. However, isotopic composition can further validate a suspected UNE signature, but separation from any atmospheric background composition is first necessary. To date, evaluating the challenges of performing this separation has typically involved comparing an observed background with a highly idealized deterministic model of radioxenon signature production by a UNE that does not consider the influence of post-detonation chemical/physical processes in the detonation cavity or the subsequent gas transport mechanisms that can also affect the isotopic composition of the detected gas signature. In addition, purely deterministic models, as previously employed, overlook the uncertainty inherent in estimating critical parameters characterizing the UNE and its detonation environment. In this paper, we create detailed, multi-parameter models of radionuclide evolution using the widely accepted England and Rider post-detonation radionuclide decay-chain network coupled to detailed models simulating physical production and transport processes affecting the gas signature. Because these models are governed by uncertain parameters including barometric fluctuations, realistic ranges of variation for each of the parameters influencing isotopic composition are then defined. A Latin-Hypercube sampling approach is used to obtain a random distribution of isotopic production and gas transport results associated with a given value of each parameter. We apply these results to background histories of two stations, one providing 4-isotope background measurements and the other providing two-isotope measurements associated with the 2013 DPRK announced UNE.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2023

Kinetically versus thermodynamically controlled factors governing elementary pathways of GaP(111) surface oxidation

GaP and related III-V semiconductors have attracted interest as high-efficiency photoelectrochemical electrodes for solar water splitting. However, their efficacy is linked to the presence, identity, and integrity of native surface oxides, which are structurally and chemically complex and evolve during operation. Using ambient pressure X-ray photoelectron spectroscopy (APXPS) coupled with ab initio simulations, we track key chemical motifs expressed during evolution of GaP(111) surface oxides, their associated reaction kinetics, and their correlations with electronic properties. We identify two distinct thermal regimes corresponding to kinetically and thermodynamically controlled oxidation. Below 600 K, exposure to O2 generates kinetically facile Ga-O-Ga configurations, whereas higher temperatures cause activated oxygen to insert into Ga–P bonds as part of a thermodynamically driven transformation into a complex, heterogeneous 3D network of surface POx (1≤x ≤ 4) groups and Ga2O3 species, the latter of which eventually dominates upon depletion of surface phosphorus. Here our study highlights the critical competition between kinetic and thermodynamic factors during GaP oxidation, yielding insights for fabricating stable III-P-based photoelectrodes with precisely engineered surface properties.

08 HYDROGEN↗
DOE OSTI2026

Arctic Deployment of a Fully Integrated Self-Powered Drifting Buoy Harvesting Wave Energy via a Triboelectric Nanogenerator

The Arctic Ocean remains one of the most poorly sampled regions on Earth, where improved in situ environmental monitoring is vital for advancing oceanographic and atmospheric studies. However, data collection efforts are constrained by the short operational lifespans and high costs of conventional systems. Drifting buoys powered by pendulum-driven wave energy harvesters offer a cost-effective alternative, yet earlier designs have neither been optimized for real-world wave conditions nor validated in the Arctic. In this study, we develop a self-powered drifting buoy that integrates a pendulum-driven triboelectric nanogenerator (TENG) system with a mechanical motion rectifier, a high-gear-ratio transmission, and power management circuits. Through coupled buoy–pendulum dynamic simulations and laboratory testing using a motion simulator, we identify an optimal pendulum mass of 1.6 kg (12.7% of total buoy weight) that maximizes energy output while maintaining buoy stability. Laboratory experiments achieved average power outputs of 12.7 mW under Arctic-like wave and temperature conditions. The system was successfully deployed in the Bering Sea, where it generated 11 J of energy in 3.1 m waves, marking the first Arctic deployment of a TENG-based drifting buoy for sea surface temperature monitoring. This work establishes a cost-effective framework for designing self-powered Arctic monitoring platforms and advances the feasibility of long-term environmental observations in real Arctic waters.

marine enerby↗
DOE OSTI2020

Benchmark of interfacial area concentration approaches for the two-fluid model in gas-dispersed condensing flow

Accurate predictions of interfacial area concentration (IAC) are required for closure of the condensation phase change term. The interfacial area transport equation (IATE) can dynamically predict interfacial area concentration and is increasingly being used in the benchmarking of phase change flows over static based correlations which require a flow regime map. For higher flow regimes, the two-group two-fluid IATE model can separately model the transport of spherical bubbles and cap/slug bubbles. Recently, a modified formulation of the one-group and two-group two-fluid model with interfacial area transport equation was proposed for gas-dispersed condensing flows. Three datasets have been used for validation consisting of condensing two-phase flow data in the unheated section of a vertical annulus with conditions spanning a range of pressures, inlet subcoolings and mass fluxes for benchmarking the various IAC approaches. For bubbly flows, the coupled void transport - IATE simulation results in good agreement in the predictions of the void fraction and the interfacial area concentration under nearly all sixty conditions. Considering the proposed two-group two-fluid IATE model formulation, the group-2 Nusselt number correlation is investigated to account for the appropriate heat transfer length scale for group-2 bubbles. Validation of the proposed two-group two-fluid model against the new dataset shows good agreement with experimental data with significant improvement over the one-group IATE model for conditions with substantial group-2 void fraction throughout the simulation domain.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗