SEARCH · Engineering Papers

Results for “spatial heterogeneity”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 595 records · Page 33

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Physics-informed semantic inpainting: Application to geostatistical modeling

A fundamental problem in geostatistical modeling is to infer the heterogeneous geological field based on limited measurements and some prior spatial statistics. Semantic inpainting, a technique for image processing using deep generative models, has been recently applied for this purpose, demonstrating its effectiveness in dealing with complex spatial patterns. However, the original semantic inpainting framework incorporates only information from direct measurements, while in geostatistics indirect measurements are often plentiful. In this work, to overcome this limitation, we propose a physics-informed semantic inpainting framework, employing the Wasserstein Generative Adversarial Network with Gradient Penalty (WGAN-GP) and jointly incorporating the direct and indirect measurements by exploiting the underlying physical laws. Our simulation results for a high-dimensional problem with 512 dimensions show that in the new method, the physical conservation laws are satisfied and contribute in enhancing the inpainting performance compared to using only the direct measurements.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

IDAES GTEP 0.1 Release

The Institute for the Design of Advanced Energy Systems (IDAES) Generation and Transmission Expansion Planning (GTEP) package provides a Pyomo-based implementation of a modular, flexible, Generalized Disjunctive Programming (GDP) formulation for power infrastructure planning problems. This package is designed with the following goals in mind: - Abstract GTEP modeling away from any particular case study or fixed modeling assumptions (e.g., technologies, temporal resolution, spatial resolution, policy implications, etc.) - Admit flexible decision sets and heterogeneous parameterization - Allow high-level modeling options to be understood easily, chosen modularly, and changed rapidly

AS↗

DOE OSTI2023

Seismic Characterization of the Blue Mountain Geothermal Field

Subsurface characterization is crucial for geothermal energy exploration and production. Yet hydrothermal reservoirs usually reside in highly fractured and faulted zones where accurate characterization is very challenging because of low signal-to-noise ratios of land seismic data and lack of coherent reflection signals. We perform an active-source seismic characterization for the Blue Mountain geothermal field in Nevada using active seismic data to reveal the elastic medium property complexity and fault distribution at this field. We first employ an unsupervised machine learning method to attenuate groundroll and near-surface guided-wave noise and enhance coherent reflection and scattering signals from noisy seismic data. We then build a smooth initial P-wave velocity model based on an existing magnetotellurics survey result, and use 3D first-arrival traveltime tomography to refine the initial velocity model. We then derive a set of elastic wave velocities and anisotropic parameters using elastic full-waveform inversion, and obtain PP and PS images using elastic reverse-time migration. We identify major faults by analyzing the variations of seismic velocities and anisotropy parameters, and reveal mid- to small-scale faults by applying a supervised machine learning method to the seismic migration images. Our characterization reveals complex velocity heterogeneities and anisotropies, as well as faults, with a high spatial resolution. These results can provide valuable information for optimal placement of future injection and production wells to increase geothermal energy production at the Blue Mountain geothermal power plant.

58 GEOSCIENCES↗

NASA NTRS1978

Nitrogen abundances and isotopic compositions in stony meteorites

Nitrogen contents range from a few parts per million in ordinary chondrites and achondrites to several hundred parts per million in enstatite chondrites and carbonaceous chondrites. Four major isotopic groups are recognized: (1) C1 and C2 carbonaceous chondrites, (2) enstatite chondrites, (3) C3 chondrites, and (4) ordinary chondrites. The major variations are primary, representing isotopic abundances established at the time of condensation and accretion. Secondary processes, such as spallation reactions, solar wind implantation and metamorphic loss may cause small but observable isotopic variations in particular cases. The large isotopic difference between enstatite chondrites and carbonaceous chondrites cannot be accounted for by equilibrium condensation from a homogeneous nebular gas, and requires either unusually large kinetic effects, or a temporal or spatial variation of isotopic composition of the nebula. Nitrogen isotopic heterogeneity in the nebula due to nuclear processes has not been firmly established, but may be required to account for the large variations found within the Allende and Leoville meteorites. The unique carbonaceous chondrite, Renazzo, also requires a special source. It is not yet possible, from the meteoritic data, to establish the mode of accretion of nitrogen onto the primitive earth.

Kung, C.-C.↗

NASA NTRS1981

The two-dimensional diagnostic model for tropospheric OH - An uncertainty analysis

Tropospheric OH is believed to play a major role in atmospheric photochemistry, and, via its chemical interactions OH ultimately may affect the climate, stratospheric ozone levels, and the pH of rain. A valuable tool for predicting the global OH abundance as a function of latitude and altitude is the two-dimensional diagnostic model for tropospheric OH. This model can also be applied to budget studies of species which interact with OH, such as CH4, CO, and O3. In addition to laboratory-obtained kinetic data, important variables include the rate constants for several reactions, the rate of heterogeneous removal of soluble species, and the global abundance and spatial variability in H2O, O3, CO, and NOx. The uncertainties in these parameters imply over a factor of 3 variability in the calculated global OH concentration as well as significant uncertainties in the latitudinal and vertical OH distribution.

Chameides, W. L.↗

NASA NTRS2019

Understanding the Impacts of Soil Moisture Initial Conditions on NWP in the Context of Land-Atmosphere Coupling

The role of soil moisture in NWP has gained more attention in recent years, as studies have demonstrated impacts of land surface states on ambient weather from diurnal to seasonal scales. However, soil moisture initialization approaches in coupled models remain quite diverse in terms of their complexity and observational roots, while assessment using bulk forecast statistics can be simplistic and misleading. In this study, a suite of soil moisture initialization approaches is used to generate short-term coupled forecasts over the U.S. Southern Great Plains using NASA’s Land Information System (LIS) and NASA Unified WRF (NU-WRF) modeling systems. This includes a wide range of currently used initialization approaches, including soil moisture derived from ‘off-the-shelf’ products such as atmospheric models and land data assimilation systems, high-resolution land surface model spinups, and satellite-based soil moisture products from SMAP. Results indicate that the spread across initialization approaches can be quite large in terms of soil moisture conditions and spatial resolution, and that SMAP performs well in terms of heterogeneity and temporal dynamics when compared against high resolution land surface model and in-situ soil moisture estimates. Case studies are analyzed using the local land-atmosphere coupling (LoCo) framework that relies on integrated assessment of soil moisture, surface flux, boundary layer, and ambient weather, with results highlighting the critical role of inherent model background biases. In addition, simultaneous assessment of land vs. atmospheric initial conditions in an integrated, process-level fashion can help address the question of whether improvements in traditional NWP verification statistics are achieved for the right reasons.

Southern Great Plains↗

DOE OSTI2023

TRACER-VOC: Speciated Volatile Organic Compounds at the Tracking Aerosol Convection Interactions Experiment Field Campaign Report

Formation of secondary organic aerosol (SOA) is responsible for the growth of particles into cloud condensation nuclei (CCN)-relevant sizes and accounts for a substantial fraction of aerosol mass in the troposphere. Despite its wide abundance and importance, capturing SOA in regional and climate models is still challenging due to complex chemical transformations involved in production of SOA. Oxidation products of biogenic hydrocarbons, such as isoprene, monoterpenes, and sesquiterpenes account for the largest fraction of SOA mass over the continents. As such, biogenic SOA is an important source of continental CCN, and has been recently shown to nucleate ice in the deposition mode. To better account for the processes involved in SOA formation in global and regional climate models, we need better measurements of sources and transformations of volatile organic compounds (VOCs), the precursors to SOA formation in the atmosphere. Understanding aerosol-deep convection interactions is an important science driver of the Tracking Aerosol Convection Interactions Experiment (TRACER), requiring a characterization of the critical controls on aerosol properties and processes, and how those properties and processes vary spatially and temporally in the Houston area. In that environment, the atmospheric photo-oxidative chemistry that drives the evolution and transformation of aerosol properties and processes is strongly impacted by the emission of VOCs, which have both anthropogenic (e.g., industrial, urban, transportation, biomass burning) and biogenic (e.g., marine) sources in the region. Understanding (and disentangling) the complex role of VOCs in such a heterogenous region requires high-time-resolution, speciated VOC measurements that are spatially resolved. To meet these measurement requirements, the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility’s time-of-flight proton transfer reaction-mass spectrometer (ToF-PTR-MS) was deployed during the TRACER intensive operational period (IOP) from June to September, 2022 at the ancillary site (ANC), which was sited to provide a characterization of the rural atmospheric environment in contrast to the main site.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1986

Strain accumulation and surface deformation along the San Andreas, California

Stressing and rupture of a locked zone adjacent to a creeping fault segment was studied with special reference to strength heterogeneity depthwise and along-strike. The resulting precursory temporal and spatial variations of surface strain rate profiles were compared to geodetic measurements on the San Andreas fault in central California. Crustal deformation in great California earthquake cycles was also studied with special reference to the temporal decay of strain rate observed since the 1957 and 1906 great earthquakes, and comtemporary surface strain rate and velocity profiles at several locations along the San Andreas. The effect of viscoelastic response in the deep aseismic shear zone on the surface deformation behavior was examined. Work was begun on a fundamental reformulation of the crustal deformation problem focusing on the crustal deformation process affected by deep aseismic slip as the slip zone progresses toward an instability and as deep seismic slip continues postseismically, the 3-D nature of the problem due to geometry and material heterogeneity, and the time-dependent source coming from the lithosphere/astenospheric coupling process.

Li, Victor C.↗

DOE OSTI2024

Upsampling Monte Carlo Reactor Simulation Tallies in Depleted Sodium-Cooled Fast Reactor Assemblies Using a Convolutional Neural Network

The computational demand of neutron Monte Carlo transport simulations can increase rapidly with the spatial and energy resolution of tallied physical quantities. Convolutional neural networks have been used to increase the resolution of Monte Carlo simulations of light water reactor assemblies while preserving accuracy with negligible additional computational cost. Here, we show that a convolutional neural network can also be used to upsample tally results from Monte Carlo simulations of sodium-cooled fast reactor assemblies, thereby extending the applicability beyond thermal systems. The convolutional neural network model is trained using neutron flux tallies from 300 procedurally generated nuclear reactor assemblies simulated using OpenMC. Validation and test datasets included 16 simulations of procedurally generated assemblies, and a realistic simulation of a European sodium-cooled fast reactor assembly was included in the test dataset. We show the residuals between the high-resolution flux tallies predicted by the neural network and high-resolution Monte Carlo tallies on relative and absolute bases. The network can upsample tallies from simulations of fast reactor assemblies with diverse and heterogeneous materials and geometries by a factor of two in each spatial and energy dimension. The network’s predictions are within the statistical uncertainty of the Monte Carlo tallies in almost all cases. This includes test assemblies for which burnup values and geometric parameters were well outside the ranges of those in assemblies used to train the network.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2020

Uncovering Uranium Isotopic Heterogeneity of Fuel Pellets from the Fifth Collaborative Materials Exercise of the Nuclear Forensics International Technical Working Group

In 2017, the Nuclear Forensics International Technical Working Group (ITWG) organized their fifth 37 Collaborative Materials Exercise (CMX-5). The exercise samples were two uranium dioxide fuel pellets 38 manufactured from the same starting materials by different processes to have similar bulk isotopic 39 composition, but different spatial uranium isotopic distributions. Sets of identical materials were sent to 40 all participating laboratories, who then utilized their existing nuclear forensic capabilities to 41 independently analyse fuel pellets and identify similarities and differences of the materials’ 42 characteristics. The analytical methods used to probe the fuel pellets included ex situ, such as sectioning 43 or breaking up the pellets and analyzing dissolved pieces using inductively coupled plasma mass 44 spectrometry (ICP-MS), analyzing particles collected from intact or fragmented pellets by secondary ion 45 mass spectrometry (SIMS), as well as in situ methods, such as laser ablation coupled with ICP-MS, 46 autoradiography and nanoSIMS. In this paper we present the results of these independent analyses and 47 compare the capabilities of those nuclear forensic analytical methods to uncover details of the isotopic 48 heterogeneity of uranium fuel pellets.

Nuclear Forensic Analysis of Uranium Fuel Pellets,↗

DOE OSTI2021

Emerging contact force heterogeneity in ordered soft granular media

Under external perturbations, inter-particle forces in disordered granular media are well known to form a heterogeneous distribution with filamentary patterns. Better understanding these forces and the distribution is important for predicting the collective behavior of granular media, the media second only to water as the most manipulated material in global industry. However, studies in this regard so far have been largely confined to granular media exhibiting only geometric heterogeneity, leaving the dimension of mechanical heterogeneity a rather uncharted area. In this work, through a FEM contact mechanics model, we show that a heterogeneous inter-particle force distribution can also emerge from the dimension of mechanical heterogeneity alone. Specifically, we numerically study inter-particle forces in packing of mechanically heterogeneous disks arranged over either a square or a hexagonal lattice and under quasi-static isotropic compression. Our results show that, at the system scale, a hexagonal packing exhibit a more heterogeneous inter-particle force distribution than a square packing does; At the particle scale, for both packing lattices, preliminary analysis shows the consistent coexistence of outliers (i.e., softer disks sustaining larger forces while stiffer disks sustaining smaller forces) in comparison to their homogeneous counterparts, which implies the existence of nonlocal effect. Further analysis on the portion of outliers and on spatial contact force correlations suggest that the hexagonal packing shows more pronounced nonlocal effect over the square packing under small mechanical heterogeneity. However, such trend is reversed when assemblies becomes more mechanically heterogeneous. Lastly, we confirm that, in the absence of particle reorganization events, contact friction merely plays the role of packing stabilization while its variation has little effect on inter-particle forces and their distribution.

36 MATERIALS SCIENCE↗

NASA NTRS2016

Understanding Cirrus Ice Crystal Number Variability for Different Heterogeneous Ice Nucleation Spectra

Along with minimizing parameter uncertainty, understanding the cause of temporal and spatial variability of the nucleated ice crystal number, Ni, is key to improving the representation of cirrus clouds in climate models. To this end, sensitivities of Ni to input variables like aerosol number and diameter provide valuable information about nucleation regime and efficiency for a given model formulation. Here we use the adjoint model of the adjoint of a cirrus formation parameterization (Barahona and Nenes, 2009b) to understand Ni variability for various ice-nucleating particle (INP) spectra. Inputs are generated with the Community Atmosphere Model version 5, and simulations are done with a theoretically derived spectrum, an empirical lab-based spectrum and two field-based empirical spectra that differ in the nucleation threshold for black carbon particles and in the active site density for dust. The magnitude and sign of Ni sensitivity to insoluble aerosol number can be directly linked to nucleation regime and efficiency of various INP. The lab-based spectrum calculates much higher INP efficiencies than field-based ones, which reveals a disparity in aerosol surface properties. Ni sensitivity to temperature tends to be low, due to the compensating effects of temperature on INP spectrum parameters; this low temperature sensitivity regime has been experimentally reported before but never deconstructed as done here.

Cirrus↗

DOE OSTI2026

Spatial Correlations of the Poisson Model for Radiation Transport

Characterizing the relationship between bulk physical properties and mixing in randomly heterogeneous media is a central challenge across many areas of science and engineering. A benchmark model for such studies is the Poisson model, a random tessellation of space by a Poisson process of hyperplanes. In radiation transport studies, the lack of exact expressions for the Poisson model’s spatial multipoint functions has led to approximate methods being used, introducing unquantified sources of error. Here, we recently introduced an exact solution for the Poisson model’s multipoint functions and closely related conditional probability functions (CPFs), providing a new opportunity to understand and reduce these sources of error. In this paper, we enable a more rigorous investigation of radiation transport in stochastic media by applying the recently introduced exact solution for the Poisson model’s CPFs. This paper consists of three main contributions. First, we introduce a unified framework for CPFs of the Poisson model, encompassing the recently introduced exact CPFs as well as the previously introduced atomic mix, nearest-neighbor, and combination CPFs. This framework also includes existing pruning techniques for the approximate CPFs, such as angular exclusion, as well as a novel form of angular exclusion suitable for the exact CPFs. Second, we use the exact CPFs to characterize the spatial regions where each approximate three-point CPF is most accurate, thereby explaining the observed hierarchy of accuracy among the approximate models. Finally, we evaluate material transmittance, reflectance, and flux in a three-dimensional test problem using conditional point sampling, demonstrating the relationship between CPF accuracy and transport simulation accuracy.

Poisson model↗

NASA NTRS1992

Some practicable applications of quadtree data structures/representation in astronomy

Development of quadtree as hierarchical data structuring technique for representing spatial data (like points, regions, surfaces, lines, curves, volumes, etc.) has been motivated to a large extent by storage requirements of images, maps, and other multidimensional (spatially structured) data. For many spatial algorithms, time-efficiency of quadtrees in terms of execution may be as important as their space-efficiency concerning storage conditions. Briefly, the quadtree is a class of hierarchical data structures which is based on the recursive partition of a square region into quadrants and sub-quadrants until a predefined limit. Beyond the wide applicability of quadtrees in image processing, spatial information analysis, and building digital databases (processes becoming ordinary for the astronomical community), there may be numerous further applications in astronomy. Some of these practicable applications based on quadtree representation of astronomical data are presented and suggested for further considerations. Examples are shown for use of point as well as region quadtrees. Statistics of different leaf and non-leaf nodes (homogeneous and heterogeneous sub-quadrants respectively) at different levels may provide useful information on spatial structure of astronomical data in question. By altering the principle guiding the decomposition process, different types of spatial data may be focused on. Finally, a sampling method based on quadtree representation of an image is proposed which may prove to be efficient in the elaboration of sampling strategy in a region where observations were carried out previously either with different resolution or/and in different bands.

Pasztor, L.↗

NASA NTRS2017

Scale Dependence of Cirrus Horizontal Heterogeneity Effects on TOA Measurements: MODIS Brightness Temperatures in the Thermal Infrared - Part I

This paper presents a study on the impact of cirrus cloud heterogeneities on MODIS simulated thermal infrared (TIR) brightness temperatures (BTs) at the top of the atmosphere (TOA) as a function of spatial resolution from 50 meters to 10 kilometers. A realistic 3-D (three-dimensional) cirrus field is generated by the 3DCLOUD model (average optical thickness of 1.4, cloudtop and base altitudes at 10 and 12 kilometers, respectively, consisting of aggregate column crystals of D (sub eff) equals 20 microns), and 3-D thermal infrared radiative transfer (RT) is simulated with the 3DMCPOL (3-D Monte Carlo Polarized) code. According to previous studies, differences between 3-D BT computed from a heterogenous pixel and 1-D (one-dimensional) RT computed from a homogeneous pixel are considered dependent at nadir on two effects: (i) the optical thickness horizontal heterogeneity leading to the plane-parallel homogeneous bias (PPHB); and the (ii) horizontal radiative transport (HRT) leading to the independent pixel approximation error (IPAE). A single but realistic cirrus case is simulated and, as expected, the PPHB mainly impacts the low-spatial resolution results (above approximately 250 meters), with averaged values of up to 5-7 K (thousand), while the IPAE mainly impacts the high-spatial resolution results (below approximately 250 meters) with average values of up to 1-2 K (thousand). A sensitivity study has been performed in order to extend these results to various cirrus optical thicknesses and heterogeneities by sampling the cirrus in several ranges of parameters. For four optical thickness classes and four optical heterogeneity classes, we have found that, for nadir observations, the spatial resolution at which the combination of PPHB and HRT effects is the smallest, falls between 100 and 250 meters. These spatial resolutions thus appear to be the best choice to retrieve cirrus optical properties with the smallest cloud heterogeneity-related total bias in the thermal infrared. For off-nadir observations, the average total effect is increased and the minimum is shifted to coarser spatial resolutions.

TOA↗

DOE OSTI2024

Impact of Pressure-Dependent Interfacial Tension and Contact Angle on Capillary Heterogeneity Trapping of CO2 in Storage Aquifers

Summary Carbon dioxide (CO2) capillary trapping increases the total amount of CO2 that can be effectively immobilized in storage aquifers. This trapping, manifesting itself as accumulated CO2 columns at a continuum scale, is because of capillary threshold effects that occur below low-permeability barriers. Considering that capillary pressure is dictated by heterogeneous pore throat size, the trapped CO2 column height and associated CO2 saturation will vary spatially within a storage aquifer. This variation will be influenced by two pressure-dependent interfacial parameters—CO2/brine interfacial tension (IFT) and CO2/brine/rock contact angle. Our objective is to understand how the pressure dependence of these two parameters affects the heterogeneity of capillary trapped CO2 at a continuum scale. Our conceptual model is a 1D two-zone system with the upper zone being a flow barrier (low permeability) and the lower zone being a flow path (high permeability). The inputs to this model include microfacies-dependent capillary pressure vs. saturation curves and permeability values. The input capillary pressure curves were collected in the literature that represents carbonate microfacies (e.g., dolograinstone) in a prevalent formation in the Permian Basin. We then used the Leverett j-function to scale the capillary pressure curve for the two zones that are assigned with the same or different microfacies. During scaling, we considered the influence of pressure on both the IFT and contact angle of CO2/brine/dolomite systems. We varied the zone permeability contrast ratio from 2 to 50. We then assumed capillary gravity equilibriums and calculated the CO2 saturation buildup corresponding to various trapped CO2 column heights. The CO2 saturation buildup is defined as the CO2 saturation in the lower layer minus that in the upper one. We found that the saturation buildup can be doubled when varying pressure in a storage aquifer, after considering pressure-dependent IFT and contact angles. Thus, assuming these two parameters to be constant across such aquifers would cause large errors in the quantification of capillary trapping of CO2. The whole study demonstrates the importance of considering pressure-dependent interfacial properties in predicting the vertical distribution of capillary trapped CO2. It has important implications in developing a better understanding of leakage risks and consequent storage safety.

Engineering↗

NASA NTRS2008

Stochastic Transport Theory for Investigating the Three-Dimensional Canopy Structure from Space Measurements

Radiation reflected from vegetation canopies exhibits high spatial variation. Satellite-borne sensors measure the mean intensities emanating from heterogeneous vegetated pixels. The theory of radiative transfer in stochastic media provides the most logical linkage between satellite observations and the three-dimensional canopy structure through a closed system of simple equations which contains the mean intensity and higher statistical moments directly as its unknowns. Although this theory has been a highly active research field in recent years, its potential for satellite remote sensing of vegetated surfaces has not been fully realized because of the lack of models of a canopy pair-correlation function that the stochastic radiative transfer equations require. The pair correlation function is defined as the probability of finding simultaneously phytoelements at two points. This paper presents analytical and Monte Carlo generated pair correlation functions. Theoretical and numerical analyses show that the spatial correlation between phytoelements is primarily responsible for the effects of the three-dimensional canopy structure on canopy reflective and absorptive properties. The pair correlation function, therefore, is the most natural and physically meaningful measure of the canopy structure over a wide range of scales. The stochastic radiative transfer equations naturally admit this measure and thus provide a powerful means to investigate the three-dimensional canopy structure from space. Canopy reflectances predicted by the stochastic equations are assessed by comparisons with the PARABOLA measurements from coniferous and broadleaf forest stands in the BOREAS Southern Study Areas. The pair correlation functions are derived from data on tree structural parameters collected during field campaigns conducted at these sites. The simulated canopy reflectances compare well with the PARABOLA data.

Huang, Dong↗

DOE OSTI2020

Effects of temperature on helium bubble behaviour in Fe–9Cr alloy

In the present study, Fe–9Cr model alloys were irradiated by 275 keV He + ions at 500 to 800 °C to a peak damage level and implanted helium content of 0.5 dpa and 0.75 at% He, respectively. Effects of temperature on He bubble behaviour were studied by transmission electron microscopy (TEM) to reveal the temperature effects. The bubble spatial distribution changes from homogenous at 500 and 600 °C to heterogeneous at 700 and 800 °C, due to preferential bubble precipitation on dislocations and grain boundaries at high temperature. Faceted bubbles with low energy faces along the {001} planes were present at all four temperatures and are attributed to surface energy minimization. Bubble precipitation on the GBs became more pronounced with increasing temperature. Meanwhile, the cavity swelling increased with increasing temperature, from 0.59% at 500 °C to 8.21% at 800 °C with a sharp increase above 700 °C. By analysing the temperature dependence of bubble size and density, a transition temperature between low and high temperature regimes was obtained in terms of different activation energies which are related to different bubble nucleation and formation mechanisms.

36 MATERIALS SCIENCE↗