Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Conditional diffusion models”

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,027 records · Page 57

Chromium segregation-induced oxide evolution in Ni-10Cr alloys during high-temperature oxidation

The oxidation behavior of a Ni-10(wt%)Cr alloy under high-temperature O 2 conditions is investigated using transmission electron microscopy and first-principles calculations. Results reveal that chromium segregation plays a central role in driving the evolution of complex oxide phase structures during oxidation. At low Cr concentrations, Cr preferentially segregates to NiO grain boundaries or internal pores, substituting for Ni atoms and forming Ni(Cr)O solid solutions. As Cr content increases, enhanced diffusion promotes Cr penetration into the NiO lattice, leading to the formation of multiphase oxide structures. First-principles modeling corroborates these findings: at low Cr concentrations, Cr atoms favor surface and grain-boundary segregation, while higher concentrations lead to Cr aggregation within the NiO bulk. Furthermore, the integrated experimental-theoretical approach provides atomistic insights into Cr-mediated mass transport mechanisms during alloy oxidation and offers valuable guidance for controlling oxide growth kinetics and phase stability in Ni-Cr alloys, with implications for improving oxidation resistance in high-temperature structural applications.

36 MATERIALS SCIENCE↗

Self-consistent description of electrokinetic phenomena in particle-based simulations

A new computational method is presented for study suspensions of charged particles undergoing fluctuating hydrodynamic and electrostatic interactions. The proposed model is appropriate for polymers, proteins, and porous particles embedded in a continuum electrolyte. A self-consistent Langevin description of the particles is adopted in which hydrodynamic and electrostatic interactions are included through a Green's function formalism. An Ewald-like split is adopted in order to satisfy arbitrary boundary conditions for the Stokeslet and Poisson Green functions, thereby providing a formalism that is applicable to any geometry and that can be extended to deformable objects. The convection-diffusion equation for the continuum ions is solved simultaneously considering Nernst-Planck diffusion. In conclusion, the method can be applied to systems at equilibrium and far from equilibrium. Its applicability is demonstrated in the context of electrokinetic motion, where it is shown that the ionic clouds associated with individual particles can be severely altered by the flow and concentration, leading to intriguing cooperative effects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spectrally resolved modulated infrared radiometry of photothermal, photocarrier, and photoluminescence response of CdSe crystals: Determination of optical, thermal, and electronic transport parameters

Spectrally resolved modulated infrared radiometry (SR-MIRR) with super-band gap photoexcitation is introduced as a self-consistent method for semiconductor characterization (CdSe crystals grown under different conditions). Starting from a theoretical model combining the contributions of the photothermal (PT) and photocarrier (PC) signal components, an expression is derived for the thermal-to-plasma wave transition frequency f{sub tc} which is found to be wavelength-independent. The deviation of the PC component from the model at high frequency is quantitatively explained by a quasi-continuous distribution of carrier recombination lifetimes. The integral, broad frequency band (0.1 Hz–1 MHz) MIRR measurements simultaneously yielded the thermal diffusivity a, the effective IR optical absorption coefficient β{sub eff}, and the bulk carrier lifetime τ{sub c}. Spectrally resolved frequency scans were conducted with interchangeable IR bandpass filters (2.2–11.3 μm) in front of the detector. The perfect spectral match of the PT and PC components is the direct experimental evidence of the key assumption in MIRR that de-exciting carriers are equivalent to blackbody (Planck) radiators. The exploitation of the β spectrum measured by MIRR allowed determining the background (equilibrium) free carrier concentration n{sub 0}. At the shortest wavelength (3.3 μm), the photoluminescence (PL) component supersedes the PC one and has distinct features. The average sample temperature influences the PC component but not the PT one.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Electron kinetics in a positive column of AC discharges in a dynamic regime

Abstract We have performed hybrid kinetic-fluid simulations of a positive column in alternating current (AC) argon discharges over a range of driving frequencies f and gas pressure p for the conditions when the spatial nonlocality of the electron energy distribution function (EEDF) is substantial. Our simulations confirmed that the most efficient conditions of plasma maintenance are observed in the dynamic regime when time modulations of mean electron energy (temperature) are substantial. The minimal values of the root mean square electric field and the electron temperature have been observed at f/p values of about 3 kHz Torr −1 in a tube of radius R = 1 cm. The ionization rate and plasma density reached maximal values under these conditions. The numerical solution of a kinetic equation allowed accounting for the kinetic effects associated with spatial and temporal nonlocality of the EEDF. Using the kinetic energy of electrons as an independent variable, we solved an anisotropic tensor diffusion equation in phase space. We clarified the role of different flux components during electron diffusion in phase space over surfaces of constant total energy. We have shown that the kinetic theory uncovers a more exciting and rich physics than the classical ambipolar diffusion (Schottky) model. Non-monotonic radial distributions of excitation rates, metastable densities, and plasma density have been observed in our simulations at pR > 6 Torr cm. The predicted off-axis plasma density peak in the dynamic regime has never been observed in experiments so far. We hope our results stimulate further experimental studies of the AC positive column. The kinetic analysis could help uncover new physics even for such a well-known plasma object as a positive column in noble gases.

Physics↗

Amorphization and siliconization of silicon carbide as a first wall material

The understanding and prediction of silicon carbide (SiC) material evolution exposed to SOL plasma conditions is of prime interest because SiC represents a promising main chamber wall plasma-facing material for next-step fusion devices (low hydrogenic diffusion, good mechanical and thermal properties under neutron irradiation). Gross and net Si erosion rates from SiC surfaces in contact with a well-diagnosed L-mode plasma in the DIII-D tokamak have been simulated and the surface concentrations of impurities have been tracked as a function of time. Coupled simulation of surface model and impurity transport demonstrates amorphization of crystalline SiC exposed to L-mode plasma due to the accumulation of displacement damages under ion irradiation. This affects the lifetime of SiC plasma facing components. Surface evolution is tightly coupled to impurity transport in the plasma and therefore needs to be integrated with impurity transport simulations to effectively predict Si erosion rates and sub-surface concentrations as a function of time. The simulation workflow couples a semi-analytical surface model to the impurity transport code GITR. The surface model is a homogeneous mixed-material model that tracks physical & chemical sputtering and reflection of impurities. Gross erosion is primarily influenced by the background plasma parameters and redeposition patterns are mainly influenced by the prompt redeposition due to the gyro-orbits of impurity ions. Although crystalline form of SiC is preferable for fusion wall applications because of resistance to neutron irradiation, this work indicates that crystalline SiC will undergo amorphization under D plasma contact with implications of higher sputtering and fuel retention. These results direct us to explore the effects of amorphization on crystalline SiC and further the physics basis of SiC usage as first wall material for fusion environments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Comparison of fission product release predictions using PARFUME with results from the AGR-1 safety tests

Safety tests were conducted on fourteen fuel compacts from AGR-1, the first irradiation experiment of the Advanced Gas Reactor (AGR) Fuel Development and Qualification program, at temperatures ranging from 1600 to 1800°C to determine fission product release at temperatures that bound reactor accident conditions. The PARFUME (PARticle FUel ModEl) code was used to predict the release of fission products silver, cesium, strontium, and krypton from fuel compacts containing tristructural isotropic (TRISO) coated particles during the safety tests, and the predicted values were compared with experimental results. Preliminary comparisons between PARFUME predictions and post-irradiation examination (PIE) results of the safety tests show different trends in the prediction of the fractional release depending on the species, and it leads to different conclusions regarding the diffusivities used in the modeling of fission product transport in TRISO-coated particles: • For silver, the diffusivity in silicon carbide (SiC) might be over-estimated by a factor of at least 102 to 103 at 1600°C and 1700°C, and at least 10 to 102 at 1800°C. The diffusivity of silver in uranium oxy-carbide (UCO) might also be over-estimated, but the available data are insufficient to allow definitive conclusions to be drawn. • For cesium, the diffusivity in UCO might be over-estimated by a factor of at least 102 to 103 at 1600°C, 105 at 1700°C, and 103 at 1800°C. The diffusivity of cesium in SiC might also over-estimated, by a factor of 10 at 1600°C and 103 at 1700°C, based upon the comparisons between calculated and measured release fractions from intact particles. There is no available estimate at 1800°C since all the compacts heated up at 1800°C contain particles with failed SiC layers whose release dominates the release from intact particles. • For strontium, the diffusivity in SiC might be over-estimated by a factor of 10 to 102 at 1600 and 1700°C, and 102 to 103 at 1800°C. These values might be somewhat over-estimated because the strontium retention during irradiation cannot be assessed a priori, which affects the magnitude of the calculated release during safety testing. The diffusivity of strontium in UCO cannot be derived from these heating tests, but it is assumed to be modeled correctly using the IAEA recommended value for kernel diffusivity. • For krypton, there is no reliable release data for compacts heated up at 1600°C, which includes all the compacts containing only intact particles. At 1700 and 1800°C, comparisons show an over-prediction of the release from compacts containing particles with failed SiC by 1 to 1.5 orders of magnitude. The available data from these heating tests do not allow to determine which of the TRISO-coating’s layers diffusivities are under or over-estimated.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

An Infrared Study of the Dust Properties and Geometry of the Arched Filaments H ii Region with SOFIA/FORCAST

Massive stellar clusters provide radiation (∼10{sup 7}--10{sup 8} L{sub ⊙}) and winds (∼1000 km s{sup −1}) that act to heat dust and shape their surrounding environment. In this paper, the Arched Filaments in the Galactic center were studied to better understand the influence of the Arches cluster on its nearby interstellar medium (ISM). The Arched Filaments were observed with the Faint Object InfraRed CAMera for the SOFIA Telescope at 19.7, 25.2, 31.5, and 37.1 μm. Color–temperature maps of the region created with the 25.2 and 37.1 μm data reveal relatively uniform dust temperatures (70–100 K) over the extent of the filaments (∼25 pc). Distances between the cluster and the filaments were calculated assuming equilibrium heating of standard-size ISM dust grains (∼0.1 μm). The distances inferred by this method are in conflict with the projected distance between the filaments and the cluster, although this inconsistency can be explained if the characteristic grain size in the filaments is smaller (∼0.01 μm) than typical values. DustEM models of selected locations within the filaments show evidence of depleted abundances of polycyclic aromatic hydrocarbons (PAHs) by factors of ∼1.6–10 by mass compared to the diffuse ISM. The evidence for both PAH depletion and a smaller characteristic grain size points to processing of the ISM within the filaments. We argue that the eroding of dust grains within the filaments is not likely attributable to the radiation or winds from the Arches cluster, but may be related to the physical conditions in the Galactic center.

79 ASTRONOMY AND ASTROPHYSICS↗

Multiscale Modeling of CO 2 Migration and Trapping in Fractured Reservoirs with Validation by Model Comparison and Real-Site Applications (Final Technical Report)

This report documents the accomplishments achieved during the project titled “Multiscale Modeling of CO 2 Migration and Trapping in Fractured Reservoirs with Validation by Model Comparison and Real-Site Applications” funded by the US Department of Energy, Office of Fossil Energy. The objectives of the project were to develop modeling capabilities for predicting CO 2 and brine migration in fractured reservoirs during geologic carbon storage and to investigate the feasibility of carbon storage in fractured reservoirs using the newly developed modeling capabilities. To achieve these objectives new mass transfer functions were developed to adapt existing dual-continuum approaches – a commonly used approach in hydrocarbon reservoir modeling – to the CO 2 -brine system. In the dual-continuum approach the undisturbed rock matrix and the fractures are modeled as two separate continua, which are coupled through mass exchange of CO 2 and brine between the two continua. Mass transfer functions represent the mass exchange between fractures and the rock matrix, without requiring highly resolved numerical grids of individual fractures, thus allowing large modeling domains relevant to geologic carbon storage questions (tens of meters to hundreds of kilometers). Mass transfer functions for the exchange of both free-phase CO 2 and brine between fractures and rock matrix were developed for both the injection phase and the post-injection phase. During the injection phase the fractures are filled with CO 2 , while the rock matrix is initially filled with brine. Therefore, a mass transfer function was developed that is based on vertical displacement of brine by CO 2 in the rock matrix due to gravitational forces. Once CO 2 injection ceases and CO 2 migrates away from the direct vicinity of the injection site, the fractures will begin to re-fill with brine, while the rock matrix will have high CO 2 saturations. For these conditions a mass transfer functions was developed that is based on spontaneous imbibition, as the pore spaces in the rock matrix are much smaller than in the fractures. For formations with low permeability or high capillary entry pressure of the rock matrix, free-phase CO 2 may not be able to enter the rock matrix. However, free-phase CO 2 at the interface between fractures and rock matrix can dissolve into brine resident in the rock matrix. Therefore, additional mass transfer functions were developed based on aqueous-phase diffusion of CO 2 through the rock matrix. These mass transfer functions are based on existing solutions for diffusion and were expanded to be applicable to a variety of block shapes relevant to dual-continuum modeling. Implementations of the newly developed mass transfer functions showed that they give accurate results when compared to highly resolved models. In addition to developing new transfer functions, a new vertically-integrated dual-continuum modeling framework was developed. This development takes advantage of fast vertical segregation of CO 2 and brine in the fractures due to the density difference between CO 2 and brine and the high permeability of the fractures. Therefore, a vertically-integrated vertical equilibrium model for the fractures is coupled with a vertically-distributed model for the rock matrix. The newly developed approach significantly reduces the computational effort, but is only accurate for fracture permeabilities of 100 mD or higher. Finally, the newly developed models were applied to real and hypothetical geologic carbon storage sites to investigate the feasibility of using fractured reservoirs for carbon storage. The general results were, that storage of CO 2 – both as free-phase and dissolved – is feasible, but that storage capacities are likely lower than in unfractured systems. While the high permeability of the fractures reduces the injection pressure (or permits higher injection rates), fast migration of CO 2 may limit the overall storage volume, because CO 2 may escape laterally with not enough of the rock matrix being swept for effective storage.

20 FOSSIL-FUELED POWER PLANTS↗

Predicting secondary organic aerosol phase state and viscosity and its effect on multiphase chemistry in a regional-scale air quality model

Atmospheric aerosols are a significant public health hazard and have substantial impacts on the climate. Secondary organic aerosols (SOAs) have been shown to phase separate into a highly viscous organic outer layer surrounding an aqueous core. This phase separation can decrease the partitioning of semi-volatile and low-volatile species to the organic phase and alter the extent of acid-catalyzed reactions in the aqueous core. A new algorithm that can determine SOA phase separation based on their glass transition temperature (T g ), oxygen to carbon (O:C) ratio and organic mass to sulfate ratio, and meteorological conditions was implemented into the Community Multiscale Air Quality Modeling (CMAQ) system version 5.2.1 and was used to simulate the conditions in the continental United States for the summer of 2013. SOA formed at the ground/surface level was predicted to be phase separated with core–shell morphology, i.e., aqueous inorganic core surrounded by organic coating 65.4 % of the time during the 2013 Southern Oxidant and Aerosol Study (SOAS) on average in the isoprene-rich southeastern United States. Our estimate is in proximity to the previously reported ~70 % in literature. The phase states of organic coatings switched between semi-solid and liquid states, depending on the environmental conditions. The semi-solid shell occurring with lower aerosol liquid water content (western United States and at higher altitudes) has a viscosity that was predicted to be 10 2 –10 12 Pa s, which resulted in organic mass being decreased due to diffusion limitation. Organic aerosol was primarily liquid where aerosol liquid water was dominant (eastern United States and at the surface), with a viscosity <10 2 Pa s. Phase separation while in a liquid phase state, i.e., liquid–liquid phase separation (LLPS), also reduces reactive uptake rates relative to homogeneous internally mixed liquid morphology but was lower than aerosols with a thick viscous organic shell. The sensitivity cases performed with different phase-separation parameterization and dissolution rate of isoprene epoxydiol (IEPOX) into the particle phase in CMAQ can have varying impact on fine particulate matter (PM 2.5 ) organic mass, in terms of bias and error compared to field data collected during the 2013 SOAS. This highlights the need to better constrain the parameters that govern phase state and morphology of SOA, as well as expand mechanistic representation of multiphase chemistry for non-IEPOX SOA formation in models aided by novel experimental insights.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Scaling of thermal stratification in outlet plena of SFRs with gallium as a surrogate fluid

Due to its high thermal diffusivity, liquid sodium can develop thermally stratified regions in the hot and cold pools of SFRs under certain off normal operating conditions. The temperature gradients in the fluid give rise to mixed convection flow in which flow velocities result both from pressure gradients and from volume forces induced by density gradients. High thermal diffusivity and low momentum diffusivity leads to reduced thermal mixing and enhanced stratification. Thermal stratification, if not looked carefully in the design, can lead to severe consequences in normal operation with forced coolant circulation and during off-normal natural circulation. The resulting hydraulic and thermal loads can cause the failure modes in structures and components critical for reactor safety. Thermal stratification in the SFR pools can lead to thermal fatigue, which can lead to structural and material failure, or pose serious concerns for the reactor safety. In addition to thermo-structural issues, the reactivity changes are significantly impacted due to density or temperature driven movement of liquid sodium. Temperature changes or gradients in the Upper Instrumentation Structures (UIS) of SFRs can lead to thermal contraction, a case that can lead to positive reactivity insertion after core shut-down. The motivation of this work is to design a scaled down experimental facility which can generate high fidelity data for thermal stratification physics to validate 3-dimensional (3D) Computational Fluid Dynamics (CFD) models. These validated CFD codes can be used to improve the 1-dimensional (1D) system level codes such as SAS4A/SASSYS-1 which are currently not capable of simulating thermal stratification accurately. These improvements in 1D system level SAS4A/SASSYS-1 codes can significantly improve the safety analyses and design of SFRs. The GaTE facility outlet plenum model is designed based on physics based similarities. Thermal stratification physics is governed primarily by the Richardson numbers (Ri). The expected low flow rates during off-normal conditions can lead to a high degree of thermal stratification in the outlet plenum with a corresponding Ri = 530. Therefore, the inlet flow velocities and inlet hydraulic diameter for the gallium coolant entering the test vessel are required to achieve same Ri. The design of test vessel in GaTE facility is finalized to a 1/20. scale outlet plenum height of the ABTR. This design of the test facility will allow the decoupling of the thermal stratification physics at low flow rates with the transient natural circulation. In other words, it will enable effective control to study thermal stratification at very low flow rates. Numerical studies are conducted with a 3D CFD code to examine the thermal behavior of the selected test design. Simulation results show that thermal gradients are present only in the vertical or axial direction which are similar to the behavior reported in previous studies on full scale SFR models. Future work will involve the experimental studies using the GaTE facility with the high fidelity instrumentation.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Dry Reforming of Ethane and Butane with CO 2 over PtNi/CeO 2 Bimetallic Catalysts

Dry reforming is a potential process to convert CO 2 and light alkanes into syngas (H 2 and CO), which can be subsequently transformed to chemicals and fuels. Here in this work, PtNi bimetallic catalysts have been investigated for dry reforming of ethane and butane using both model surfaces and supported powder catalysts. The PtNi bimetallic catalyst shows an improvement in both activity and stability as compared to the corresponding monometallic catalysts. The formation of PtNi alloy and the partial reduction of Ce 4+ to Ce 3+ under reaction conditions are demonstrated by in-situ Ambient Pressure X-ray Photoemission Spectroscopy (AP-XPS), X-ray Diffraction (XRD) and X-ray Absorption Fine Structure (XAFS) measurements. A Pt-rich bimetallic surface is revealed by Diffuse Reflectance Infrared Fourier Transform Spectroscopy (DRIFTS) following CO adsorption. Combined in-situ experimental results and Density Functional Theory (DFT) calculations suggest that the Pt-rich PtNi bimetallic surface structure would weaken the binding of surface oxygenates/carbon species and reduce the activation energy for C-C bond scission, leading to an enhanced dry reforming activity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

FunDiff: diffusion models over function spaces for physics-informed generative modeling

Recent advances in generative modeling-particularly diffusion models and flow matching-have been widely used for synthesizing discrete data such as images and videos. However, adapting these models to physical applications remains challenging, as the quantities of interest are continuous functions governed by complex physical laws. To address this, we introduce FunDiff, an efficient and robust framework for generative modeling in function spaces. FunDiff combines a latent diffusion process with a function autoencoder architecture to handle input functions with varying discretizations, generates continuous functions that can be evaluated at arbitrary locations, and seamlessly incorporate physical priors. These priors are enforced through architectural constraints or physics-informed loss functions, ensuring that generated samples satisfy fundamental physical laws. We theoretically establish minimax optimality guarantees for density estimation in function spaces, demonstrating that diffusion-based estimators achieve optimal convergence rates under suitable regularity conditions. We further demonstrate the practical effectiveness of FunDiff across diverse applications in fluid dynamics and solid mechanics. Empirical results indicate that our method can generate physically consistent samples with high fidelity to the target distribution, and exhibit robustness to noisy and low-resolution data.

Wang, Sifan [Yale University, New Haven, CT (Unite↗

The physics of long- and intermediate-wavelength asymmetries of the hot spot: Compression hydrodynamics and energetics

To achieve ignition with inertial confinement fusion (ICF), it is important to under- stand the effect of asymmetries on the hydrodynamics and energetics of the compres- sion. This paper describes a theoretical model for the compression of distorted hot spots, and quantitative estimates using hydrodynamic simulations. The asymmetries are categorized into low (Ι < 6) and intermediate (Ι < A < 40) modes by comparison of the wavelength with the thermal-diffusion scale length. Long-wavelength modes introduce substantial nonradial motion, whereas intermediate-wavelength modes in- volve more cooling by thermal ablation. We discover that for distorted hot spots, the measured neutron-averaged properties can be very different from the real hydro- dynamic conditions. This is because mass ablation driven my thermal conduction introduces flows in the Rayleigh–Taylor bubbles, this results in pressure variation, in addition to temperature variation between the bubbles and the neutron-producing region (~1 keV for intermediate modes). The differences are less pronounced for long-wavelength asymmetries since the bubbles are relatively hot and sustain fusion reactions. The yield degradation$-$ with respect to the symmetric$-$ results primarily from a reduction in the hot-spot pressure for low modes and from a reduction in burn volume for intermediate modes. It is shown that the degradation in internal energy of the hot-spot is equivalent for both categories, and is equal to the total residual energy in the shell including the bubbles. This quantity is correlated with the shell residual kinetic energy for low-modes, and includes the kinetic energy in the bubbles for mid-modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Mechanical behavior at high temperatures of highly oxygen- or hydrogen-enriched α and prior-β phases of Zirconium alloys - 2016-0063

During a hypothetical loss-of-coolant accident, zirconium alloy fuel claddings can be loaded by internal pressure and exposed to steam at high temperatures (HTs) up to 1,200 deg. C and then cooled and quenched in water. A significant fraction of the oxygen that reacts with the cladding during HT oxidation diffuses beneath the oxide through the metallic substrate. This induces a progressive transformation of the metallic β{sub Zr}-phase layer into an intermediate layer of α{sub Zr}(O) phase containing up to 7 weight % (wt.%) oxygen. Furthermore, in some specific conditions, the cladding may rapidly absorb a significant amount of hydrogen during steam exposition at HTs. As a β{sub Zr} stabilizer, hydrogen would mainly diffuse and concentrate up to several thousands of parts per million by weight (wppm) into the inner β{sub Zr}-phase layer. This study intends to provide new and more comprehensive data on the HT mechanical behavior of the α{sub Zr}(O) and (prior-)β{sub Zr} phases containing high contents of oxygen and hydrogen, respectively. Model samples produced from M5{sup R} and Zircaloy-4 cladding tubes homogeneously charged in oxygen (less than or equal to 6 wt.%) and hydrogen (less than or equal to 3,000 wppm), respectively, were prepared. Their mechanical behavior was determined under vacuum between 800 and 1,100 deg. C for the oxygen-enriched α{sub Zr} phase and in air between 700 and 20 deg. C after cooling from the β{sub Zr} temperature domain for the hydrogen-enriched (prior-)β{sub Zr} phase. The α{sub Zr} phase was substantially strengthened and embrittled by oxygen. The contribution of the α{sub Zr}(O) layer to the HT creep behavior of an oxidized fuel cladding tube subjected to internal pressure is evaluated by finite element analysis. Mechanical strength and ductility of the model (prior-)β{sub Zr} phase appear to be affected by hydrogen contents of 2000-3000 wppm in ways that depend on temperature. (authors)

36 MATERIALS SCIENCE↗

Physical and Thermodynamic Properties of Secondary Organic Materials for Modeling (Final Technical Report)

Secondary organic material constitutes more than 50% of the submicron component of atmospheric aerosol particles. This organic material is produced by oxidation of biogenic and anthropogenic volatile organic compounds. The physicochemical properties of these aerosol particles and their variability with atmospheric conditions influence many climate-relevant endpoints, such as rates of chemical reactions, production of secondary organic material, number-diameter distributions, direct radiative forcing, and cloud condensation nuclei. The results of the project research highlight a nonlabile-to-labile transition for material representative of anthropogenic secondary organic particulate matter at a threshold humidity range. This behavior differs markedly from materials representing biogenic aerosol sources, which lack such a transition. More specifically, the mass lability of atmospheric organic particulate matter, meaning its tendency to evaporate, is an important property in the mechanisms governing the climate effects. This property is affected by both the volatility and the molecular diffusivity of the constituent organic species. Via a series of laboratory experiments, the project directly determined the mass labilities of films of secondary organic material representative of anthropogenic and biogenic organic particulate matter (PM). The evaporation rates for anthropogenic organic PM increased above a threshold relative humidity (RH) between 20% and 30%, indicating rapid partitioning above a transition RH. Below the threshold, the characteristic time for equilibration was estimated at up to 1 week for a typically sized particle. In contrast, for films representing biogenic PM, no RH threshold was observed, suggesting equilibrium partitioning is rapidly obtained for all RHs. The effective molecular diffusivity for the biogenic case was at least 1000 times greater than that of the anthropogenic case. The results allow for an improved representation of anthropogenic and biogenic secondary organic aerosol formation in climate models. The results of the project also provide new measurements on how organic aerosol particles take up water and serve as cloud condensation nuclei, thus explaining the observed complex thermodynamic behaviors using a simple thermodynamic model. The project results reconciled the discrepancy between sub-saturation hygroscopicity and super-saturation cloud condensation nuclei activity for organic aerosol particles. The key to the explanation is the phase separation between hydrophobic and hydrophilic organics and the induced surface tension lowering. The results are essential for the assessment of the aerosol indirect effect on climate. In more detail, hygroscopic growth and cloud condensation nuclei activation are key processes for accurately modeling the climate impacts of organic particulate matter. Nevertheless, the microphysical mechanisms of these processes remain unresolved. The project results focus on complex thermodynamic behaviors, including humidity-dependent hygroscopicity, diameter-dependent cloud condensation nuclei activity, and liquid-liquid phase separation in the laboratory for biogenically derived secondary organic material representative of similar atmospheric organic particulate matter. These behaviors can be explained by the non-ideal mixing of water with hydrophobic and hydrophilic organic components. The non-ideality-driven liquid-liquid phase separation further enhances water uptake and induces lowered surface tension at high relative humidity, which results in a lower barrier to cloud condensation nuclei activation. By comparison, secondary organic material representing anthropogenic sources does not exhibit complex thermodynamic behavior. The combined results highlight the importance of detailed thermodynamic representations of the hygroscopicity and cloud condensation nuclei activity in models of the Earth’s climate system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Design of anti-ring back reflectors for thin-film solar cells based on three-dimensional optical and electrical modeling

The optical and electrical properties of a photonic-plasmonic nanostructure on the back contact of thin-film solar cells were investigated numerically through the three-dimensional (3D) finite-difference time-domain method and the 3D Poisson and drift-diffusion solver. The focusing effect and the Fabry-Perot resonances are identified as the main mechanisms for the enhancement of the optical generation rate as well as the short circuit current density. However, the surface topography of certain nanopattern structures is found to reduce the internal electrostatic field of the device, thus limiting charge collection. The optimized conditions for both optics and electronics have been analyzed in this paper.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Heterogeneous Two-Phase Pillars in Epitaxial NiFe 2 O 4 -LaFeO 3 Nanocomposites

Self-assembled epitaxial oxide nanocomposites have been explored for a wide range of applications, including multiferroic and magnetoelectric properties, plasmonics, and catalysis. These so-called “vertically aligned nanocomposites” form spontaneously during the deposition process when segregation into two phases is energetically favorable as compared to a solid solution. However, there has been surprisingly little work understanding the driving forces that govern the synthesis of these materials, which can include point defect energetics, surface diffusion, and interfacial energies. To explore these factors, La-Ni-Fe-O films have been synthesized by molecular beam epitaxy and it is shown that these phase segregate into spinel-perovskite nanocomposites. Using complementary scanning transmission electron microscopy and atom-probe tomography, the elemental composition of each phase is examined and found that Ni ions are exclusively found in the spinel phase. From correlative analysis, a model for the relative favorability of the Ni2+ and Ni3+ valences under the growth conditions is developed. It is shown that multidimensional characterization techniques provide previously unobserved insight into the growth process and complex driving forces for phase segregation.

36 MATERIALS SCIENCE↗

Colloidal State Machines as Smart Tracers for Chemical Reactor Analysis

A widely utilized tool in reactor analysis is passive tracers that report the residence time distribution, allowing estimation of the conversion and other properties of the system. Recently, advances in microrobotics have introduced powered and functional entities with sizes comparable to some traditional tracers. This has motivated the concept of Smart Tracers that could record the local chemical concentrations, temperature, or other conditions as they progress through reactors. Herein, the design constraints and advantages of Smart Tracers by simulating their operation in a laminar flow reactor model conducting chemical reactions of various orders are analyzed. It is noted that far fewer particles are necessary to completely map even the most complex concentration gradients compared with their conventional counterparts. Design criteria explored herein include sampling frequency, memory storage capacity, and ensemble number necessary to achieve the required accuracy to inform a reactor model. Cases of severe particle diffusion and sensor noise appear to bind the functional upper limit of such probes and require consideration for future design. The results of the study provide a starting framework for applying the new technology of microrobotics to the broad and impactful set of problems classified as chemical reactor analysis.

97 MATHEMATICS AND COMPUTING↗