Search NASA⌕ Search

SEARCH · Search NASA

Results for “BOUNDARY”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Elucidating grain boundary energy minimization mechanisms in textured Ca-doped alumina with inclination-dependent Monte Carlo Potts simulations

The grain growth behavior of textured Ca-doped alumina is compared to Monte Carlo Potts (MCP) simulations to investigate the effect of anisotropic grain boundary (GB) properties on local boundary migration. Experimentally, the growth of textured Ca-doped alumina results in highly elongated grains. The relative GB energy distribution is measured using the thermal groove method before and after heat treating at 1600°C, finding that high energy GBs are eliminated during grain growth. No significant difference in the GB energy distributions is found between the long and short axes of the elongated grains, suggesting that anisotropic mobility may be responsible for the grain shape. However, MCP simulations with anisotropic mobility as a function of plane inclination do not result in grains with distinct morphologies, regardless of the degree of anisotropy introduced. The final grain shape after grain growth of textured Ca-doped alumina resembles that of the MCP simulations using an anisotropic GB energy as a cosine function of plane inclination. Several energy functions are tested and only those that mathematically impose a torque (second derivative of energy with respect to the plane inclination angle) result in elongated grains. Although area reduction is the dominant energy minimization mechanism, these results suggest that local GB migration is affected by anisotropic GB energy and torque and alternative mechanisms like GB replacement and reorientation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Quantifying and visualizing the microscopic degrees of freedom of grain boundaries in the Wigner–Seitz cell of the displacement-shift-complete lattice

Here we introduce a grain boundary (GB) translation vector, t ws , to describe and quantify the domain of the microscopic degrees of freedom of GBs. It has long been recognized that for fixed macroscopic degrees of freedom of a GB there exists a large multiplicity of states characterized by different relative grain translations. More recently another degree of freedom, [n], the number of GB atoms, has emerged and is now recognized as an equally important component of GB structural multiplicity. In this work, we show that all GB microstates can be uniquely characterized by their value of t ws , which is located within the Wigner–Seitz (WS) cell of the Displacement-Shift-Complete lattice (DSCL) of the GB. The GB translation vector captures information about both the translation state and the number of GB atoms. We show that the density of GB microstates inside the WS cell of the DSCL is not uniform and can form clusters that correspond to different GB phases. Further, the vectors connecting the centers of the clusters correspond to the Burgers vectors of GB phase junctions, which can be predicted without building the junctions. Using t ws , we quantify GB excess shear and argue that it is defined up to a DSCL vector, which has implications for thermodynamic equilibrium conditions. Additionally, this work generalizes the definition of the number of GB atoms [n] to asymmetric boundaries.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Reduction of vanadium diffusivity within copper grain boundaries due to enhanced binding

We investigate permeation of vanadium (V) through copper (Cu) layers in V/Cu/niobium (Nb) and molybdenum (Mo)/Cu/V tri-layer composites and compare it with the permeation behavior of Nb. Despite their low solubility in single-crystal Cu, both V and Nb permeate readily due to short-circuit transport through Cu grain boundaries (GBs). However, Nb permeates preferentially to V. Atomistic simulations suggest that this difference is due to lower V diffusivity at general GBs in Cu, compared to Nb. Here, the lower V diffusivity, in turn, arises from stronger V binding to the boundary, which manifests as a greater equilibrium GB concentration as well as a higher activation energy for GB diffusion.

36 MATERIALS SCIENCE↗
DOE OSTI2026

GBOpt: Grain boundary structure optimization using Monte Carlo and evolutionary algorithms

Polycrystalline materials are made of many small crystals separated by grain boundaries (GBs), whose atomic structure strongly influences material properties. Because the structure of a GB determines its properties, the optimal structure must be known in order to determine those impacts. There are many ways of placing atoms in the GB region, but the optimal structure is defined as the one that gives the lowest value of a target property (typically energy). GB structure optimization has been successfully demonstrated using stochastic and evolutionary methods, but no reusable, community-maintained open-source workflow has been developed. GBOpt (Grain Boundary Optimization) is an open-source Python package that creates that workflow, where we have presently implemented two approaches: Markov Chain Monte Carlo, and genetic algorithm based on elite selection. We demonstrate this capability by successfully reproducing the known optimal structures of a specific GB in two materials, and point interested readers to the GitHub repository for additional examples, including optimization for different properties. Both of the implemented approaches recovered the known structures, with the genetic algorithm approach finding the optimal structure faster on average.

99 - GENERAL AND MISCELLANEOUS↗
DOE OSTI2025

On nonlocal problems with Neumann boundary conditions: scaling and convergence for nonlocal operators and solutions

Formulations of Neumann-type boundary conditions for boundary value problems in the nonlocal framework are beset with difficulties, some related to the choice of a proper scaling. Here we identify a space-dependent scaling for a nonlocal Neumann operator, for which we prove linear in δ (δ being the radius for the support for the kernel) convergence of the Neumann operator and $\mathcal{O}$(δ 2 ) convergence of solutions to their classical counterparts. The pointwise-like convergence of the nonlocal normal operator is cast as a new type of two-scale operator-point convergence, which we call condensated convergence . The results hold for general integrable kernels, a setting which is favored in numerical simulations. We support this analysis with numerical convergence studies using a piecewise linear discontinuous Galerkin discretization and show an $\mathcal{O}$(δ 2 ) rate of convergence of solutions, also exhibiting an $\mathcal{O}$(h 2 ) convergence, where h is the mesh size.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Outflow Boundary Conditions for Turbine-Integrated Rotating Detonation Combustors

This study examines outflow boundary conditions (BCs) in computational fluid dynamics (CFD) simulations of a transition duct with and without guide vanes that converts supersonic flow exiting a rotating detonation combustor (RDC) to subsonic flow to drive a turbine. Since the flow exiting the transition duct has swirling shock waves with significant spatial and temporal variations in pressure, temperature, and Mach number, imposing proper BCs poses a challenge. To ensure all swirling shock waves exit the transition duct without creating non-physical reflected waves at its outlet, this study examined three outflow BCs: (1) the average pressure imposed at the duct’s outlet, (2) a nonreflecting BC (NRBC) with a specified average pressure imposed at the duct’s outlet, (3) the average pressure imposed at the outlet of an extension duct made up of a buffer layer and a sponge layer. This study is based on the three-dimensional, unsteady density-weighted-ensemble-averaged continuity, Navier–Stokes, and energy equations for a thermally perfect gas closed by the realizable k–ε model and “enhanced” wall functions. The results obtained show that imposing an average pressure at the transition duct’s outlet produces spurious waves that degrade the physical meaningfulness of the solution. When the NRBC was applied, swirling shock waves exited the duct’s outlet without creating spurious waves. However, its usage requires the gas to be thermally, as well as calorically, perfect, which this study shows could be a concern. By imposing the average pressure at the outlet of an extension duct, the gas does not need to be calorically perfect. The results obtained show the effects of the sponge layer’s length and coarsening ratio on damping nonuniformities in non-physical reflected waves to ensure the flow exiting the transition duct’s outlet can do so as if there are no boundaries present and has the desired average pressure—even though the BC is applied at the extension duct’s outlet.

gas turbines↗
DOE OSTI2025

Bulk‐Boundary Correspondence of Semimetal Ru 3 Sn 7 and Topological Surface States on Chemically Realistic Terminations

Ru 3 Sn 7 is experimentally demonstrated as an efficient catalyst, with potential utilization of topological surface states for hydrogen evolution reaction. Despite its promising catalytic performance, the topological nature of Ru 3 Sn 7 remains uncertain. Particularly, the bulk-boundary correspondence has not yet been established, hence hindering a rigorous justification of its topologically-protected surface states. In this work, the bulk topology of Ru 3 Sn 7 is detailed using first-principles calculations and the topological quantum chemistry formalism. Ru3Sn7 turns out to be an enforced semimetal possessing symmetry-protected crossings within a set of bands near the Fermi level, which are enforced and prescribed by the violations of symmetry-prescribed compatibility relations. Moreover, the surface states and the associated origin from the same set of entangled bands are identified, thereby establishing the bulk-boundary correspondence. To evaluate the effects of chemical modifications, the response of topological surface states to various surface terminations, stoichiometry, and oxidation is examined. The surface structures are globally optimized, and the phase diagrams for various experimental conditions are built. It is shown that, due to changes in the local chemical environment, the original surface states are significantly altered. Modified surface bands can be observed near the Fermi level on surface terminations that preserve the C 4v symmetry.

36 MATERIALS SCIENCE↗
DOE OSTI2026

The Spatial and Temporal Variability of the Clear Convective Boundary Layer at the ARM SGP Supersite

The convective boundary layer (CBL), also known as the mixing layer, constitutes the critical lower segment of the atmosphere that significantly influences daily human activities. The growing demand for precise weather forecasts is driven by the requirements of agriculture, transportation, and routine societal functions. Here, to enhance understanding of the CBL, this study investigates the spatiotemporal variability in the CBL and its controlling factors using four-year Doppler lidar, surface flux, and profiling measurements at five ARM Southern Great Plains sites within a 100 km radius. This investigation utilizes data collected exclusively under clear-sky conditions or scattered low-cloud conditions. Results reveal significant spatial differences in CBL evolutions. Daily mixing layer heights (MLHs) vary up to 1 km (30% of the mean) in late afternoon. There is a clear east–west contrast: western sites (C1, E32, E37) exhibit higher summer MLH (1.9–2.1 km) and vertical velocity variances (1.0–1.2 m 2 s −2 ) than eastern sites (1.6–1.8 km), reversing in winter. Temporally, the MLH peaks at 70% of the sunrise–sunset interval, the lagging heat flux (HF) peaks at 50%; and the seasonal MLH maxima lag the HF by approximately one month, influenced by nighttime PBL (planetary boundary layer) properties. The HF and lower tropospheric stability are the main factors of influence for the CBL, but site-specific dependencies highlight the critical roles of local factors, underscoring the need for including them in CBL modeling.

ARM SGP site↗
DOE OSTI2025

Morphologies of dealloying corrosion attack at grain boundaries

Dealloying corrosion at grain boundaries severely compromises the performance of polycrystalline materials across a wide variety of technological applications. The impact of this phenomenon depends upon the morphology and rate of intergranular dealloying, which can range from planar to wormhole-like patterns that rapidly advance into the alloy. Using 2D and 3D multi-phase field simulations, we reveal how diverse microstructures result from a fundamental interplay between alloy composition and a grain boundary migration mechanism that alters diffusional pathways of dealloying. Inside intergranular dealloying channels, corrosion product buildup can spawn new channels that branch into grain interiors, such that alloys can be degraded from the inside out. These processes further lead to an atypical coarsening mode assisted by diffusion in the dealloying agent. We summarize a unifying explanation for dealloying morphology selection in polycrystalline alloys, which provides an important step towards their optimization for dealloying corrosion environments.

Corrosion↗
DOE OSTI2024

A Modeling framework for flocculated cohesive sediment transport in the current bottom boundary layer

Cohesive sediment transport, where its settling velocity is controlled by the flocculation process, is a crucial component in determining biochemical cycles, fate of pollutants, and morphodynamics in many aquatic ecosystems. In this study, a modeling framework is presented to investigate how flocculation influences cohesive sediment transport in the current bottom boundary layer in dilute conditions, consistent with the calibration range of the flocculation model. From a local analysis of floc dynamics in homogenous turbulence, we identify that the floc size distribution is mainly controlled by floc cohesion and yield strength. The uncertainty in fractal dimension plays a minor role for the floc size but it influences the resulting floc density and settling velocity. The transport analysis in the current boundary layer shows that the flocculation process alters the vertical distribution of the settling velocity and hence the sediment concentration with a strong dependence on cohesion, floc yield strength, and floc structure. When the flocs are more susceptible to breaking, a well-mixed concentration profile is obtained. In contrast, for flocs with higher cohesion or yield strength, higher concentration with a sharp gradient is observed close to the bed. Overall, the settling velocity exhibits a low vertical variability within 20 % of the depth-averaged value except near the bed. Further, this suggests that using a depth-averaged settling velocity yields acceptable predictions of the sediment concentration profiles, especially for flocs with lower cohesion.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Flexible dynamic boundary microgrid operation considering network and load unbalances

Flexible microgrids with dynamic boundaries have recently been introduced in the literature. With the ability to reconfigure the topology of the microgrids dynamically through remotely controlled switches, flexible microgrids with dynamic boundaries can further improve the resiliency and energy efficiency of microgrids with distributed energy resources (DERs). This paper focuses on the optimal operation considering one of the predominant characteristics of microgrids and distribution systems – unbalanced networks and loads. In existing literature, balanced modeling of microgrids is more common due to its attractive simplicity. The three-phase power unbalance has not been considered as a constraint on the generation units in a microgrid. Further, negative sequence constraints have also been neglected. In this article, we propose a set of constraints that is specifically related to the capabilities of inverter interfaced resources to supply unbalanced current/power when the microgrid is islanded from the main distribution grid. We incorporate the new set of constraints into two optimization formulations leveraging two convex relaxations of the three-phase power flow equations: mixed-integer linear programming (MILP) and mixed-integer semidefinite programming (MISDP) that optimize the dispatch of controllable switches and DERs in the microgrid. The algorithms are then extended to networked microgrids with grid-forming sources. We test the algorithms on a realistic community microgrid model in Puerto Rico as well as standardized IEEE distribution test feeders. The testing results demonstrate the performance of the proposed algorithms. The MILP is fast and scalable, and the MISDP enforces the negative sequence voltage constraints.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Finite-element boundary-integral simulation of thin wires and inhomogeneous penetrable bodies in subsurface multilayered anisotropic media

With the prevailing presence of drilling wells near the subsurface in mature oil and gas fields, the application of electromagnetic methods can be particularly challenging where the electromagnetic field is affected by the steel casing. In the past decades, borehole-to-surface and crosswell electromagnetic methods have been utilized for monitoring of reservoir and underground CO 2 storage. This paper presents a unified finite-element boundary-integral (FEBI) method capable of simultaneously modeling the complex electromagnetic interactions between thin metallic wires (representing steel casings) with 3D trajectory and arbitrary 3D inhomogeneous penetrable bodies (such as CO 2 plumes or hydrocarbon reservoirs) within anisotropic multilayered subsurface environments. Unlike existing approaches that treat these components separately or require dense discretization, or are limited to vertical wells, our unified formulation preserves flexible electromagnetic coupling while delivering improved computational efficiency. Assuming the background formation is multilayered anisotropic media, the surface integral equation method is applied to model the thin wires and boundaries of the inhomogeneous bodies. Meanwhile, the finite element method is applied to model the volume of inhomogeneous bodies. Here, the performance of the proposed FEBI method is assessed through comparison with reference numerical results and its practical significance is demonstrated through CO 2 plume monitoring scenarios.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Heterogeneous catalysis: Optimal performance at a phase boundary?

Most of the industrially used heterogeneous catalysts have been discovered by trial and error, and despite decades of experience, the discovery of new catalysts continues to be extremely challenging. The drive to uncover guiding principles in catalyst design is more present than ever. We share a series of observations indicating that optimal catalysts typically function at characteristic phase boundaries (e.g., abrupt changes in adsorbate coverage, catalyst structure, etc.) accessed in the reaction conditions. The catalyst exploits the associated instability—the desire to exist in multiple states simultaneously—as a driving force for chemical transformations. In other words, phase boundaries are good places to start the catalyst search, and indeed, we should focus on at least two phases at once rather than just one. Here, we substantiate this claim with several studies that combine statistical operando modeling and experiments. Transpiring from these observations is a hitherto unrecognized vector in catalyst discovery.

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

Aitken Mode Aerosols Buffer Decoupled Mid‐Latitude Boundary Layer Clouds Against Precipitation Depletion

Aerosol-cloud-precipitation interactions are a leading source of uncertainty in estimating climate sensitivity. Remote marine boundary layers where accumulation mode (~100–400 nm diameter) aerosol concentrations are relatively low are very susceptible to aerosol changes. These regions also experience heightened Aitken mode aerosol (~10–100 nm) concentrations associated with ocean biology. Aitken aerosols may significantly influence cloud properties and evolution by replenishing cloud condensation nuclei and droplet number lost through precipitation (i.e., Aitken buffering). We use a large-eddy simulation with an Aitken-mode enabled microphysics scheme to examine the role of Aitken buffering in a mid-latitude decoupled boundary layer cloud regime observed on 15 July 2017 during the Aerosol and Cloud Experiments in the Eastern North Atlantic flight campaign: cumulus rising into stratocumulus under elevated Aitken concentrations (~100–200 mg -1 ). In situ measurements are used to constrain and evaluate this case study. Our simulation accurately captures observed aerosol-cloud-precipitation interactions and reveals time-evolving processes driving regime development and evolution. Aitken activation into the accumulation mode in the cumulus layer provides a reservoir for turbulence and convection to carry accumulation aerosols into the drizzling stratocumulus layer above. Further Aitken activation occurs aloft in the stratocumulus layer. Together, these activation events buffer this cloud regime against precipitation removal, reducing cloud break-up and associated increases in heterogeneity. We examine cloud evolution sensitivity to initial aerosol conditions. With halved accumulation number, Aitken aerosols restore accumulation concentrations, maintain droplet number similar to original values, and prevent cloud break-up. Without Aitken aerosols, precipitation-driven cloud break-up occurs rapidly. In this regime, Aitken buffering sustains brighter, more homogeneous clouds for longer.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

The vertical-velocity skewness in the atmospheric boundary layer without buoyancy and Coriolis effects

One of the main features of near-neutral atmospheric boundary layer (ABL) turbulence is the positive vertical velocity skewness $Sk_w$ above the roughness sublayer or the buffer region in smooth-walls. The $Sk_w$ variations are receiving renewed interest in many climate-related parameterizations of the ABL given their significance to cloud formation and to testing sub-grid schemes for Large Eddy Simulations (LES). The vertical variations of $Sk_w$ are explored here using wind tunnel and flume experiments collected above smooth, rough, and permeable-walls in the absence of buoyancy and Coriolis effects. These laboratory experiments form a necessary starting point to probe the canonical structure of $Sk_w$ as they deal with a key limiting case (i.e., near-neutral conditions). Diagnostic models based on cumulant expansions, realizability constraints, and constant mass flux approach routinely employed in the convective boundary layer as well as prognostic models based on third-order budgets are used to explain variations in $Sk_w$ for the idealized laboratory conditions. The failure of flux-gradient relations to model $Sk_w$ from the gradients of the vertical velocity variance σ$_w^2$ are explained and corrections based on models of energy transport offered. Novel links between the diagnostic and prognostic models are also featured, especially for the inertial term in the third-order budget of the vertical velocity fluctuation. The co-spectral properties of w′/σ w vs w′ 2 /σ$_w^2$ are also presented for the first time to assess the dominant scales governing $Sk_w$ in the inner and outer layers, where w′ is the fluctuating vertical velocity and σ w is the vertical velocity standard deviation.>

Boundary layer flow↗
DOE OSTI2024

Recommendations on setup in simulating atmospheric gravity waves under conventionally neutral boundary layer conditions

Wind farm-induced atmospheric gravity waves have been the subject of recent research as they can impact wind farm performance. Pressure variations associated with gravity waves can contribute to the global blockage effect and wind farm wake recovery. Therefore, accurate numerical simulation of flow fields, where wind-farm-induced gravity waves may be produced, is important. Three main considerations in such simulations are the overall domain size, the use of Rayleigh damping near domain boundaries to dampen gravity waves, and advection damping at the inlet to prevent spurious oscillations. Often these considerations are treated ad hoc rather than systematically. This work aims to test and extend the systematic modelling of internal gravity waves proposed in a preliminary investigation to modelling of both internal and trapped gravity waves. The preliminary study identifies the length scales to set the domain and damping layer sizes and the time scale to configure the Rayleigh damping coefficient but under linearly stratified conditions. Large eddy simulations of flow through a wind farm canopy are performed under conventionally neutral boundary layer (CNBL) conditions to test the validity of proposed setups for CNBL conditions. Background atmospheric parameters, such as Froude number (Fr), inversion height (H i ), and inversion layer Froude number (Fr i ) control most of the atmospheric gravity wave characteristics. We validated for CBNL conditions that the effective wavelengths of the internal gravity waves are the correct length scale to configure the domain size and damping layer thickness. Likewise, the optimum damping coefficient to dampen the internal gravity waves relates to the free atmosphere's buoyancy frequency or buoyant perturbations' time scale. We infer that the damping coefficient in the inversion layer may relate to the inversion buoyancy frequency to effectively dampen the trapped gravity waves. Moreover, the advection damping length is linked to the horizontal wavelength of the trapped gravity waves in the inversion layer to prevent spurious waves at the inlet by retaining wave energy accumulation.

17 WIND ENERGY↗
DOE OSTI2025

Size-dependent attraction of Cu solutes to clusters formed at Ag grain boundaries

We report a size-dependent solute clustering mechanism at grain boundaries in a sputtered ultrafine-grained Ag-Cu alloy, where large Cu clusters form despite weak individual solute-solute interactions. X-ray diffraction confirms limited Cu solubility in the Ag matrix, while scanning transmission electron microscopy reveals Cu clustering at both ordinary GBs and GB junctions. Density functional theory calculations show that 12-atom Cu clusters are energetically preferred, while smaller three-atom clusters are significantly less stable. Additional calculations demonstrate a marked increase in solute-cluster attraction energy with cluster size. As a result, these findings point to a previously unrecognized pathway for grain-boundary solute clustering in immiscible systems, driven by collective solute-cluster interactions, with implications for segregation behavior and stability in nanocrystalline and ultrafine-grained alloys.

Alloys↗
DOE OSTI2024

Statistically Resolved Planetary Boundary Layer Height Diurnal Variability Using Spaceborne Lidar Data

The Planetary Boundary Layer Height (PBLH) significantly impacts weather, climate, and air quality. Understanding the global diurnal variation of the PBLH is particularly challenging due to the necessity of extensive observations and suitable retrieval algorithms that can adapt to diverse thermodynamic and dynamic conditions. This study utilized data from the Cloud-Aerosol Transport System (CATS) to analyze the diurnal variation of PBLH in both continental and marine regions. By leveraging CATS data and a modified version of the Different Thermo-Dynamics Stability (DTDS) algorithm, along with machine learning denoising, the study determined the diurnal variation of the PBLH in continental mid-latitude and marine regions. The CATS DTDS-PBLH closely matches ground-based lidar and radiosonde measurements at the continental sites, with correlation coefficients above 0.6 and well-aligned diurnal variability, although slightly overestimated at nighttime. In contrast, PBLH at the marine site was consistently overestimated due to the viewing geometry of CATS and complex cloud structures. The study emphasizes the importance of integrating meteorological data with lidar signals for accurate and robust PBLH estimations, which are essential for effective boundary layer assessment from satellite observations.

54 ENVIRONMENTAL SCIENCES↗