Search NASA⌕ Search

SEARCH · Search NASA

Results for “discrete surface”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Inputs and Outputs for paper "The DESC stellarator code suite. Part 1. Quick and accurate equilibria computations"

Contains the DESC and VMEC input and output files used in the paper, as well as the plotting scripts used to create the figures in the paper (Current with the arXiv Mar 31 2022 version, https://arxiv.org/abs/2203.17173v1): 3D equilibrium codes are vital for stellarator design and operation, and high-accuracy equilibria are also necessary for stability studies. This paper details comparisons of two 3D equilibrium codes, VMEC, which uses a steepest-descent algorithm to reach a minimum-energy plasma state, and DESC, which minimizes the MHD force error in real space directly. Accuracy as measured by final plasma energy and satisfaction of MHD force balance, as well as other metrics, will be presented for each code, along with the computation time. It is shown that DESC is able to achieve more accurate solutions, especially near-axis. DESC's global Fourier-Zernike basis also yields the solution everywhere in the plasma volume, not just on discrete flux surfaces. Further, DESC can compute the same accuracy solution as VMEC in an order of magnitude less time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

On-Demand Designing of Cathode Internal Surface Architecture for Dramatic Enhancement of SOFC Performance and Durability

This project is aimed to design and modify the internal surfaces of porous composite cathode from currently commercially viable Solid Oxide Fuel Cells (SOFCs), using additive manufacturing process of Atomic Layer Deposition (ALD). The material systems being investigated are commercial composite electrodes complex three-dimensional topographies. In term of the chemistry of the ALD layer applied on the internal surface of the porous cathode, this project has employed commercially relevant electrolyte, electrocatalyst and noble metal materials set. Such materials are fully compatible with the commercial fuel cells, and this project has developed special nanostructure on the surface of the commercial composite cathodes. The formation of the designed nanoarchitecture on the surface of SOFC cathode has been achieved through precise control of ALD parameters and their effect on overall cell performance and resultant electrochemical reaction mechanism of cathodes has been investigated through full cell electrochemical performance testing and nanostructure characterization by transmission electron microscopy (TEM). Under the support of this award, following has been achieved: (1). For cathode materials in solid oxide fuel cells (SOFCs), such as perovskite mixed conductor La 0.6 Sr 0.4 Co 0.2 Fe 0.8 O 3-x (LSCF), cation surface segregation and consequently losing conductivity and active sites for the oxygen reduction reaction (ORR) are problematic. To mitigate the cation segregation and enhance SOFC durability, further decorating the internal backbone surface using the desired electrocatalysts could be a promoting approach. Commonly, the cation segregation such as Sr is very volatile, so the effective surface decoration is ideally conformal. Nevertheless, the conformal surface coating would inevitably alter the ORR pathways that initially take place on the surface of the backbone. To reveal the impact of the conformal coating on both the catalytic activity and the conductivity of the cathode, the unary electrocatalyst of Pt or CoO x , was applied to the LSCF/SDC composite electrode of inherently functional SOFCs, respectively. Both ALD coating layers evolve strong interaction with the LSCF composite cathode. Upon operations, the Pt coating layer remains conformal on LSCF grain surfaces but turns into discrete particles on SDC grain surfaces. Meanwhile, CoO x conformal coating grows to be the discrete nanograins on both the LSCF and SDC grains. ALD coating of the cathode alone reduces the ohmic resistance up to 28 % for the entire cells. The increased conductivity induced by the ALD coating of Pt or CoO x is ascribed to different mechanisms. For the inherent functional SOFCs, the present study presents a novel and feasible approach to apply a conformal, dense coating layer on the surface of a mixed conductor, simultaneously increasing the conductivity and durability of the SOFC cathode. (2). High resistance of the oxygen electrode still significantly hinders the state-of-the-art Solid Oxide Fuel Cells (SOFCs). In particular, for an oxygen electrode consisting of mixed electronic and ionic conductors, such as perovskite lanthanum strontium cobalt ferrite (LSCF), it deteriorates due to its low chemical stability of the grain surface. Such degradation is often associated with the segregation of cations. To prevent the cation surface segregation and its resultant perovskite phase decomposition, we demonstrate a conformal ultra-thin (7-10 nm) surface heterogeneous coating layer consisting of subjacent discrete Pt nanoparticles capped with a superjacent fully dense conformal CoO x layer. The performance studies indicate the ALD coating reduces the cell series resistance by up to 40 %. The conformal CoO x layer consists of randomly orientated but single-layered nanograins, with high-density intergranular and surface grain boundaries serving as the electrochemical reaction sites and facilitating mass transport. The conformal coating layer appears to have successfully suppressed the Sr outward diffusion and confined the Sr enriched layer to a ~ 2 nm interface perovskite phase between the coating layer and the LSCF grain surface. Moreover, this ultra-thin Sr enriched perovskite layer presumably possesses high oxygen vacancy and high ionic conductivity and further imposes tensile strain to the LSCF grain surfaces. With the combination of a conformal CoO x nanoionics, Sr enriched layer, and its strained interface, the ALD coating induced surface layer is estimated to have a conductivity of ~ 1.27x10 4 S/cm, which is over two orders magnitude of that from LSCF at 750 ºC.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2026

CHESS 2025: Spectrometer orthorectified at-sensor radiance from NEON AOP imaging spectroscopy surveys

This dataset provides Level 1 (L1) orthorectified at-sensor radiance derived from measurements collected by the Imaging Spectrometer-1 (NIS-1) onboard the NEON (National Ecological Observatory Network) Airborne Observation Platform (AOP) for the 2025 Colorado Headwaters Ecological Spectroscopy Study (CHESS). NIS-1 captures light reflected from the Earth’s surface in 426 discrete wavelength bands as raw digital numbers (DNs; Level 0). These data are then calibrated to physical units (uW/cm²·sr·nm) following the processing steps described in the NEON Imaging Spectrometer Level 1B Calibrated Radiance Algorithm Theoretical Basis Document (ATBD; Gallery 2022). The data delivered here are the primary inputs for the surface reflectance product in “Custom surface reflectance, shade masks, and equivalent water thickness maps for the Colorado Headwaters Ecological Spectroscopy Study” (Carroll et al. 2026). For intertemporal comparison, the radiance data here are most directly relatable to the v2 radiance data in “NEON AOP Imaging Spectroscopy Survey of Upper East River Colorado Watersheds: Raw-Space Radiance and Observational Variable Dataset” (Goulden et al. 2018), to which the same processing methodology was applied. Together, the radiance and reflectance data enable users to exploit the unique reflection signatures of different surface objects for land cover classification, foliar trait mapping, plant vigor assessment, water content estimation, trace-element identification, and other scientific applications. The data were acquired over three study domains in the Upper Gunnison river basin: the upper East River watershed (CRBU); Almont Triangle and Taylor Canyon (ALMO); and Upper Taylor River watershed (UPTA) between 2025-06-13 and 2025-07-15. Within each domain, data are delivered by flightline as orthorectified and calibrated hyperspectral rasters in Hierarchical Data Format version 5 (HDF5) format, with radiance values provided in uW/cm²·sr·nm on a fixed, uniform Universal Transverse Mercator (UTM) grid at 1 meter spatial resolution. The radiance rasters include all 426 NIS-1 spectral bands, along with associated quality-assurance (QA) and diagnostic and ancillary layers needed for atmospheric correction workflows. Orthorectified radiance is produced from pushbroom spectrometer observations by applying NEON’s radiometric calibration (including bad pixel masking, dark subtract, dark pedestal shift correction, electronic panel ghost correction, grating ghost correction, deblur correction and flat-fielding) and spectral calibration (using spectral response function band centers and full-width at half-maximum intensity), followed by geolocation and regridding to the fixed grid. CHESS Project Description: The Colorado Headwaters Ecological Spectroscopy Study (CHESS) comprised a multi-week airborne remote sensing and field observation campaign in the Upper Gunnison Basin, Colorado, conducted in June and July of 2025. Airborne remote sensing was conducted by the National Ecological Observatory Network Airborne Observation Platform (NEON AOP), concurrent with a field campaign run by the Rocky Mountain Biological Laboratory (RMBL), the Lawrence Berkeley National Laboratory (LBNL) and SLAC National Accelerator Laboratory Watershed Function Science Focus Area (SFA), and NASA-JPL (Jet Propulsion Laboratory) Earth Surface Mineral Dust Source Investigation (EMIT) program. Between June 10 and July 18, 2025, the NEON AOP flight team collected high-resolution aerial imaging spectroscopy and Light Detection and Ranging (LiDAR) data over three domains: the Upper East River (CRBU), Almont Triangle (ALMO), and the Upper Taylor Basin (UPTA). In coordination with the flights, a field campaign acquired ground-truth observations, including observations of vegetation composition, foliar traits, forest demography, and subsurface properties in 18 core sampling areas within the domains. Additional surface water observations were taken at over 380 point locations. All CHESS campaign datasets can be found within the CHESS ESS-DIVE data portal: https://data.ess-dive.lbl.gov/portals/chess. Funding Acknowledgement: Field and remote-sensing data acquisition was performed under a grant from the National Aeronautics and Space Administration (80NSSC24K1005). This work was also supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

2018 NEON and 2025 CHESS Campaigns↗
DOE OSTI2021

Microbubble dynamics and heat transfer in boiling droplets

Dissipating large heat fluxes from a surface is critically important in numerous industrial and natural applications. Boiling based spray cooling and surface texturing are two of the most promising methods being investigated to address this problem. Although our understanding on these topics has significantly improved over past decades, critical gaps remain in the knowledgebase stymieing the realization of their full potential. As an example, while bubble growth in pool boiling have been investigated in detail, comparatively little is known about how the bubbles evolve inside boiling drops. In the present work, we have investigated for the first time, the microbubble dynamics inside water droplets boiling on superhydrophilic textured substrates using high-speed X-ray phase contrast imaging (XRPCI). Our observations show that the transient bubble density variation follows similar characteristics irrespective of the texture spacing at a given surface temperature. For an example microstructure, we found that the number of discrete bubbles on the surface decreases as temperature is increased although their growth rate increases. Here, we observe that bubble growth is highly non-uniform during the lifetime of a drop on the surface. Initially, bubbles grow under diffusion-limited regime, but at later times they grow as ~t 1.45 due to combined effects of coalescence and evaporation. In some conditions, we found that bubbles shrink dramatically af- ter the initial growth spurt presumably due to severe quenching of the surface, and migration of bubbles on the surface. Using the bubble sizes, for the first time we analyzed the heat flux removed by a single bubble and also by all the bubbles at a given time. We find that the highest dissipation through latent- heat component (~600 W/cm 2 ) occurs just in the beginning and thereafter it decreases. We expect that our findings and the analysis would guide further work on the topic and will aid in the overarching goal of engineering surfaces that are more efficient in boiling heat transfer.

42 ENGINEERING↗
DOE OSTI2025

A U.S. Scientific Community Vision for Sustained Earth Observations of Greenhouse Gases to Support Local to Global Action

Managing carbon stocks in the land, ocean, and atmosphere under changing climate requires a globally‐integrated view of carbon cycle processes at local and regional scales. The growing Earth Observation (EO) record is the backbone of this multi‐scale system, providing local information with discrete coverage from surface measurements and regional information at global scale from satellites. Carbon flux information, anchored by inverse estimates from spaceborne Greenhouse Gas (GHG) concentrations, provides an important top‐down view of carbon emissions and sinks, but currently lacks global continuity at assessment and management scales (<100 km). Partial‐column data can help separate signals in the boundary layer from the overlying atmosphere, providing an opportunity to enhance surface sensitivity and bring flux resolution down from that of column‐integrated data (100–500 km). Based on a workshop held in September 2024, the carbon cycle community envisions a carbon observation system leveraging GHG partial columns in the lower and upper troposphere to weave together information across scales from surface and satellite EO data, and integration of top‐down/bottom‐up analyses to link process understanding to global assessment.

Parazoo, Nicholas C. [California Institute of Tech↗
DOE OSTI2022

Development of A Multi-Spectral Pyrometry Sensor for High-Speed Transient Surface-Temperature Measurements in Combustion-Relevant Harsh Environments

Accurate and high-speed transient surface-temperature measurements of combustion devices including internal combustion (IC) engines, gas turbines, etc., provide validation targets and boundary conditions for computational fluid dynamics models, and are broadly relevant to technology advancements such as performance improvement and emissions reduction. Development and demonstration of a multi-infrared-channel pyrometry-based optical instrument for high-speed surface-temperature measurement is described. The measurement principle is based on multi-spectral radiation thermometry (MRT) and uses surface thermal radiation at four discrete spectral regions and a corresponding emissivity model to obtain surface temperature via non-linear least squares (NLLS) optimization. Rules of thumb for specifying the spectral regions and considerations to avoid interference with common combustion products are developed; the impact of these along with linear and non-linear MRT analysis are assessed as a function of temperature and signal-to-noise ratio. A multi-start method to determine the MRT-solution global optimum is described and demonstrated. The resulting multi-channel transient pyrometry instrument is described along with practical considerations including optical-alignment drift, matching intra-channel transient response, and solution-confidence indicators. The instrument demonstrated excellent >97% accuracy and >99% 2-sigma precision over the 400–800 °C range, with ~20 µs (50 kHz, equivalent to 0.2 cad at 2000 RPM IC-engine operation) transient response in the bench validation.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

The Mass-loss History of the Red Hypergiant VY CMa

Imaging and spectroscopy of the knots, clumps, and extended arcs in the complex ejecta of VY CMa confirm a record of high mass-loss events over the past few hundred years. Hubble Space Telescope/Space Telescope Imaging Spectrograph spectroscopy of numerous small knots close to the star allow us to measure their radial velocities from the strong K i emission and determine their separate motions, spatial orientations, and time since ejecta. Their ages concentrate around 70, 120, 200, and 250 yr ago. A K i emission knot only 50 mas from the star ejected as recently as 1985–1995 may coincide with an H{sub 2}O maser. Comparison with VY CMa’s historic light curve from 1800 to the present shows several knots with ejection times that correspond with extended periods of variability and deep minima. The similarity of this correspondence in VY CMa with the remarkable recent dimming of Betelgeuse and an outflow of gas is apparent. The evidence for similar outflows from the surface of a more typical red supergiant suggests that discrete ejections are more common and surface or convective activity is a major source of mass loss for red supergiants.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Numerical methodology for treating static and dynamic dislocation problems near a free surface

Abstract Simulation of dislocation dynamics enables researchers and scientists to explore the plastic behavior of crystalline materials under loading. Analytic solutions for the stress field due to a linear dislocation segment near a free surface are case-specific, e.g. dealing with either a horizontal segment or a vertical segment, and therefore hard to implement in time-dependent dislocation dynamics simulations as different dislocation segments could be differently oriented. This article presents a generalized numerical framework to find the stress field beneath a free surface due to the presence of a dislocation segment. The framework can be expanded to non-flat surfaces. Also, three-dimensional discrete dislocation dynamics simulations are performed here, which clearly show the effect of free surfaces on the flow stress of a material.

Siddique, Abu Bakar↗
DOE OSTI2025

Computational and Experimental Investigation of Discrete “W” Turbulators for sCO2 Turbine Internal Cooling

Direct-fired supercritical CO2 (sCO2 ) cycles are of particular interest for power generation due to the potential for higher thermal efficiencies with inherent CO2 capture. Direct fired sCO2 turbines must withstand high temperatures and the density of the working fluid is much higher than conventional Brayton cycles. Based on conventional Brayton cycle applications, internal serpentine cooling channels with ribbed turbulator strips are a common approach to achieve the required cooling, while minimizing coolant utilization. Most of the previous studies of coolant channel heat transfer enhancements have focused on serpentine channels using air as the cooling media at Reynold’s numbers typical of Brayton cycles. This paper investigates the Nusselt number augmentation of coolant channels with no features, ribbed turbulators, and discrete “W” turbulators with surface roughness effects included and CO2 as the coolant. The numerical results are compared to experimental results using CO2 as the cooling media at Reynold’s numbers ranging from 100,000 to 300,000. The results show the Nusselt number augmentation plateaus with increasing Reynolds number, in contrast to the decay observed in prior publications with air as the coolant media. Additionally, the CFD results indicate the discrete “W” geometry provides a 10% increase in Nusselt number augmentation relative to the angled ribbed turbulators.

computational fluid dynamics (CFD)↗
DOE OSTI2020

Machine learning analysis of perovskite oxides grown by molecular beam epitaxy

Reflection high-energy electron diffraction (RHEED) is a ubiquitous in situ molecular beam epitaxial (MBE) characterization tool. Although RHEED can be a powerful means for crystal surface structure determination, it is often used as a static qualitative surface characterization method at discrete intervals during a growth. A full analysis of RHEED data collected during the entirety of MBE growths is made possible using principle component analysis (PCA) and $\textit{k}$-means clustering to examine significant boundaries that occur in the temporal clusters grouped from RHEED data and identify statistically significant patterns. This process is applied to data from homoepitaxial SrTiO 3 growths, heteroepitaxial SrTiO 3 grown on scandate substrates, BaSnO 3 films grown on SrTiO 3 substrates, and LaNiO 3 films grown on SrTiO 3 substrates. We report this analysis may provide additional insights into the surface evolution and transitions in growth modes at precise times and depths during growth, and that video archival of an entire RHEED image sequence may be able to provide more insight and control overgrowth processes and film quality.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A Robust Numerical Treatment of Solid-Phase Diffusion in Pseudo Two-Dimensional Lithium-Ion Battery Models

Solid-phase diffusion in active materials of lithium-ion batteries significantly affects charging and safety-related behavior of lithium-ion batteries. Therefore, it is essential to develop an efficient and robust numerical algorithm for solving solid-phase diffusion equations in physics-based battery models. In this work, we discuss the origins of numerical instabilities that can occur when solving the solid-phase diffusion equations using iterative methods. Then, in order to resolve such issues, we propose a simple numerical treatment to the surface flux term of discretized solid-phase diffusion equations. To demonstrate its numerical robustness, the proposed method is implemented into a pseudo two-dimensional (P2D) physics-based battery model and simulations are conducted at wide ranges of operating conditions. Even with extremely poor initial guesses for the Li+ concentrations of the active materials, computations using the proposed method do not diverge and the their computational speeds are comparable to those with conventional initial guesses. Comprehensive tests of the proposed method are also performed with a dynamic current profile based on US06 driving profile and a multi-stage charging profile with very high initial C-rate (12C).

battery modeling↗
DOE OSTI2026

A low-latitude, high-elevation 10Be production-rate calibration from the Rwenzori Mountains, Uganda

Cosmogenic nuclide surface-exposure dating has become a key tool in geomorphology. Robust application of the technique relies on determination of local nuclide production rates using independently dated calibration sites. However, there are few existing calibration sites in certain regions, particularly in low-latitude, high-elevation regions where cosmogenic nuclide production is theoretically more sensitive to changes in Earth's magnetic field and atmospheric conditions. We present a new cosmogenic beryllium-10 (10Be) production-rate calibration from the equatorial Rwenzori Mountains in Uganda. The calibration is based on radiocarbon dating of basal lake sediments (∼21.3 kyr BP) from Lake Mahoma (2900 m asl) and measured 10Be concentrations from boulders on the crest of a moraine which dams the lake. Our results indicate a time-integrated 10Be production rate of 15.83 ± 0.81 at/g/yr at Lake Mahoma. We then use these calibration data and two public online calculators to calculate ages for moraines in a separate Rwenzori valley for which there exists a minimum-limiting radiocarbon age constraint (∼12 kyr BP) on surficial glacial sediments. This enables us to assess the performance of different parameters used to scale production rates spatially and temporally (e.g., scaling frameworks, geomagnetic field reconstructions, atmospheric models) for surface-exposure age calculation in the region. Our results highlight the sensitivity of low-latitude cosmogenic nuclide surface-exposure ages to discrete calculation parameters, such as scaling framework or reconstructions of Earth's magnetic field, and the utility of determining production rates at sites where limiting age data exists for other nearby landforms.

Anderson, N [Dartmouth College]↗
DOE OSTI2026

High-Order Mesh r-Adaptivity with Tangential Relaxation and Guaranteed Mesh Validity

High-order meshes are crucial for achieving optimal convergence rates in curvilinear domains, preserving symmetry, and aligning with key flow features in moving mesh simulations [1], but their quality is challenging to control. In prior work, we have developed techniques based on Target-Matrix Optimization Paradigm (TMOP) to adapt a given high-order mesh to the geometry and solution of the partial differential equation (PDE) [2, 3]. Here, we extend this framework to address two key gaps in the literature for highorder mesh 𝑟-adaptivity. First, we introduce tangential relaxation on curved surfaces using solely the discrete mesh representation, eliminating the need for access to underlying geometry (e.g., CAD model). Second, we ensure a continuously positive Jacobian determinant throughout the domain. This determinant positivity is essential for using the high-order mesh resulting from 𝑟-adaptivity with arbitrary quadrature schemes in simulations. The proposed approach is demonstrated to be robust using a variety of numerical experiments.

Mathematics and Computing↗
DOE OSTI2021

Photoluminescence and Raman mapping of β-Ga 2 O 3

Semi-insulating single crystal β-Ga 2 O 3 is becoming increasingly useful as a substrate for device fabrication. Fe doping is a method for producing such substrates. Along with Fe dopants, β-Ga 2 O 3 :Fe also contains Cr 3+ . Photoluminescence (PL) emission peaks at 690 nm (1.80 eV) and 696 nm (1.78 eV), as well as a broad feature around 709 nm (1.75 eV), are observed in β-Ga 2 O 3 :Fe. PL mapping of the 690 nm emission showed high and low intensity bands due to impurity striations introduced during crystal growth. PL mapping also revealed surface defects showing broad emissions around 983 nm (1.26 eV) and 886 nm (1.40 eV) that were spatially localized, occurring at discrete spots on the sample surface. Raman mapping of an 886 nm emission center revealed peaks at 2878 and 2930 cm -1 , consistent with an organometallic or hydrocarbon compound. Raman mapping of the 983 nm center showed a peak at 2892 cm -1 . Bright UV emission centers showed Raman peaks at 2910 and 2968 cm -1 , which are attributed to Si-CH 3 groups that may originate from silica polishing compounds or annealing in a silica ampoule.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Impingement panel for a turbomachine

An integrated combustor nozzle includes a combustion liner that extends radially between an inner liner segment and an outer liner segment. The combustion liner includes a forward end portion, an aft end portion, a first side wall, and a second side wall. The aft end portion of the combustion liner defines a turbine nozzle. The integrated combustor nozzle further includes an impingement panel having an impingement plate disposed along an exterior surface of one of the inner liner segment or the outer liner segment. The impingement plate defines a plurality of impingement holes that direct coolant in discrete jets towards the exterior surface of the inner liner segment or the outer liner segment. The impingement panel is radially spaced from the exterior surface to form a cooling flow gap therebetween. The impingement panel includes a collection duct that extends from the impingement panel and defines a collection passage.

Berry, Jonathan Dwight↗
DOE OSTI2023

Impingement panel support structure and method of manufacture

An integrated combustor nozzle includes a combustion liner that extends between an inner liner segment and an outer liner segment along a radial direction. The combustion liner including a forward end portion, an aft end portion, a first side wall, and a second side wall. An impingement panel having an impingement plate disposed along an exterior surface of one of the inner liner segment or the outer liner segment. The impingement plate defines a plurality of impingement apertures that direct coolant in discrete jets towards the exterior surface of the inner liner segment or the outer liner segment. The impingement panel includes an inlet portion that extends from the impingement plate to a collection duct. The impingement panel further includes a plurality of supports spaced apart from one another. The plurality of supports extend between, and are coupled to, the inlet portion, the collection duct, and the impingement plate.

Berry, Jonathan Dwight↗
DOE OSTI2026

Robust Containment Queries over Collections of Trimmed NURBS Surfaces via Generalized Winding Numbers

Here, we propose a containment query that is robust to the watertightness of regions bound by trimmed NURBS surfaces, as this property is difficult to guarantee for in-the-wild CAD models. Containment is determined through the generalized winding number (GWN), a mathematical construction that is indifferent to the arrangement of surfaces in the shape. Applying contemporary techniques for the 3D GWN to trimmed NURBS surfaces requires some form of geometric discretization, introducing computational inefficiency to the algorithm and even risking containment misclassifications near the surface. In contrast, our proposed method leverages properties of the 3D solid angle to solve the relevant surface integral using a boundary formulation with rapidly converging adaptive quadrature. Batches of queries are further accelerated by memoizing (i.e., caching and reusing) quadrature node positions and tangents as they are evaluated. We demonstrate that our GWN method is robust to complex trimming geometry in a CAD model, and is accurate up to arbitrary precision at arbitrary distances from the surface. The derived containment query is therefore robust to model non-watertightness while respecting all curved features of the input shape.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Bridging molecular-scale interfacial science with continuum-scale models

Solid–water interfaces are crucial for clean water, conventional and renewable energy, and effective nuclear waste management. However, reflecting the complexity of reactive interfaces in continuum-scale models is a challenge, leading to oversimplified representations that often fail to predict real-world behavior. This is because these models use fixed parameters derived by averaging across a wide physicochemical range observed at the molecular scale. Recent studies have revealed the stochastic nature of molecular-level surface sites that define a variety of reaction mechanisms, rates, and products even across a single surface. To bridge the molecular knowledge and predictive continuum-scale models, we propose to represent surface properties with probability distributions rather than with discrete constant values derived by averaging across a heterogeneous surface. This conceptual shift in continuum-scale modeling requires exponentially rising computational power. By incorporating our molecular-scale understanding of solid–water interfaces into continuum-scale models we can pave the way for next generation critical technologies and novel environmental solutions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗