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,621 records · Page 90
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Slot-die-coating operability windows for polymer electrolyte membrane fuel cell cathode catalyst layers

Roll-to-roll (R2R) slot-die coating of polymer electrolyte membrane fuel cell (PEMFC) catalyst layers represents a scalable deposition method for producing 10–20 m 2 ·min -1 of catalyst-coated gas diffusion layers (GDLs). This high-throughput production technique will help lower the cost of PEMFC catalyst layers. The uniformity of the wet layer applied by slot die deposition is affected by process parameters such as substrate speed, vacuum pressure applied at the upstream meniscus, gap between the slot die lips and substrate, ink rheology, and other ink and substrate properties. The set of conditions for producing a defect-free coating with a dilute ink typically requires little to no upstream vacuum pressure, so suitable operating conditions can be found easily through trial and error and operator intuition. However, the higher viscosity of more concentrated inks dramatically shifts the range of settings that result in a homogeneous coating to higher vacuum levels, which are harder to find through hit or miss. A predictive model showing the range of operable conditions decreases material wastage inherent in experimentally searching for suitable parameters. In this study, the defect-free coating parameter window is explored experimentally and theoretically for two concentrations of PEFC cathode inks. Both a full capillary hydrodynamic model and a computationally cheaper viscocapillary model successfully predict the experimentally determined coating window within the experimental and model uncertainty limits for inks with 5.3 wt% and 12.0 wt% solids ink while maintaining the 0.1 mg Pt ·cm -2 Pt areal loading target. Here, this paper demonstrates a viable pathway for meeting the $30/kW net ultimate cost target of the United States Department of Energy (U.S. DOE) Hydrogen Fuel Cells Technologies Office (HFTO). The concentrated ink lowers the thermal energy and capital expenditure (CapEx) budget of the coating process by decreasing the amount of time, energy, and floorspace required for drying the coating.

25 ENERGY STORAGE↗
DOE OSTI2014

Hydrogen passivation of titanium impurities in silicon: Effect of doping conditions

While the contamination of solar silicon by fast diffusing transition metals can be now limited through gettering, much attention has been drawn to the slow diffusing species, especially the early 3d and 4d elements. To some extent, hydrogen passivation has been successful in healing many deep centers, including transition metals in Si. Recent deep-level transient spectroscopy (DLTS) measurements concerning hydrogen passivation of Ti revealed the existence of at least four electrical levels related to Ti{sub i}H{sub n} in the upper-half of the gap. These findings challenge the existing models regarding both the current level assignment as well as the structure/species involved in the defects. We revisit this problem by means of density functional calculations and find that progressive hydrogenation of interstitial Ti is thermodynamically stable in intrinsic and n-doped Si. Full passivation may not be possible to attain in p-type Si as Ti{sub i}H{sub 3} and Ti{sub i}H{sub 4} are metastable against dissociation and release of bond-centered protons. All DLTS electron traps are assigned, namely, E40′ to Ti{sub i}H(-/0), E170′ to Ti{sub i}H{sub 3}(0/+), E(270) to Ti{sub i}H{sub 2}(0/+), and E170 to Ti{sub i}H(0/+) transitions. Ti{sub i}H{sub 4} is confirmed to be electrically inert.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Hierarchical calibration and validation for modeling bench-scale solvent-based carbon capture. Part 1: Non-reactive physical mass transfer across the wetted wall column

A hierarchical model calibration and validation are proposed here for quantifying the confidence level of mass transfer prediction using a computational fluid dynamics (CFD) model, where the solvent-based carbon dioxide (CO 2 ) capture is simulated and the predicted results are compared to the corresponding bench-scale experimental data. Two unit problems with an increasing level of complexity are proposed to break down the complex physical/chemical processes of solvent-based CO 2 capture into relatively simpler problems to separate the effects of physical transport and chemical reaction. This paper focuses on the calibration and validation of the first unit problem, i.e., the CO 2 mass transfer across a falling ethanolamine (MEA) film in the absence of chemical reaction. This problem is investigated both experimentally and numerically using non-reactive nitrous oxide (N 2 O) as a surrogate for CO 2 . To capture motion of the gas-liquid interface, a volume of fluid method is employed with a one-fluid formulation to compute the mass transfer between the two phases. Parallel bench-scale experiments are designed and conducted to validate and calibrate the CFD models using a general Bayesian calibration approach. Two important transport parameters, Henry's constant and gas diffusivity, are calibrated to produce the posterior distributions, which will be used as input for the second unit problem to address the chemical absorption of CO 2 across the MEA falling film, where both mass transfer and chemical reaction are involved. Finally, mass transfer coefficients predicted by CFD are also compared with those predicted by the traditional/empirical correlations under both steady and wavy falling film conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Holographic Formation of Diffraction Elements for Transformation of Light Beams in Liquid Crystal – Photopolymer Compositions

A theoretical model of holographic formation of diffractive optical elements for transformation of light beam field into Bessel-like fields in liquid crystal – photopolymer (LC–PPM) composite materials with a dyesensitizer is developed. Results of numerical modeling of kinetics ofvariation of the refractive index of a material in the process of formation with different relationships between the photopolymerization rates and diffusion processes are presented. Based on the results of numerical simulation, it is demonstrated that when the photopolarization process dominates, the diffractive element being formed is distorted. This leads to a change in the light field distribution at its output and consequently, to ineffective transformation of the reading beam. Thus, the necessity of optimizing of the recording conditions and of the prepolymeric composition to increase the transformation efficiency of light beam fields is demonstrated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

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↗
DOE OSTI2015

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↗
DOE OSTI2016

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↗
DOE OSTI2023

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↗
DOE OSTI2025

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↗
DOE OSTI2014

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↗
DOE OSTI2017

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↗
NASA NTRS1985

Noble gases and the early history of the Earth: Inappropriate paradigms and assumptions inhibit research and communication

The development of models as tracers of nobel gases through the Earth's evolution is discussed. A new set of paradigms embodying present knowledge was developed. Several important areas for future research are: (1) measurement of the elemental and isotopic compositions of the five noble gases in a large number of terrestrial materials, thus better defining the composition and distribution of terrestrial noble gases; (2) determinations of relative diffusive behavior, chemical behavior, and the distribution between solid and melt of noble gases under mantle conditions are urgently needed; (3) disequilibrium behavior in the nebula needs investigation, and the behavior of plasmas and possible cryotrapping on cold nebular solids are considered.

Huss, G. R.↗
NASA NTRS1991

Nonequilibrium chemical and radiation coupling phenomena in AOTV flowfields

A flowfield model for the nonequilibrium stagnation region of high altitude entry vehicles which includes nonequilibrium chemistry, multitemperature, viscous, conduction, and diffusion effects is presented. It contains coupled nongray nonequilibrium radiative transfer for atoms and molecules and local thermodynamic nonequilibrium phenomena. Comparison with Fire 2 flight data verifies that the model is reasonably accurate; and it has been applied to two AFE trajectory points, a high speed return from Mars, a series of points at 80 km for 12 to 16 km/sec, and three altitudes at 16 km/sec. Based on these results shock slip is significant, radiation cooling/coupling is minor at AFE conditions but important by 14 km/sec and dominant at 16 km/sec, radiation for the AFE is small but important and primarily molecular, above 12 km/sec atomic radiation is a significant or dominant portion of the total heating, and local thermodynamic nonequilibrium is important and should be included in all models.

Carlson, Leland A.↗
NASA NTRS1991

The upper ionospheres of Jupiter and Saturn

The electron-density profiles of Jupiter's and Saturn's ionospheres are modeled with a 1D chemical diffusive model incorporating measured parameters of the neutral atmospheric structure. The Voyager RSS ion-density data are modeled by accounting for ion chemistry, the H3(+) recombination-rate coefficient, and H2O chemistry. Electron-density peaks for both ionospheres are 900-1000 km lower in the model than the measured values, and the role of vibrational excitation of H2 is discussed in converting H(+) to H2(+) and H3(+). Vertical ion flow is also considered which can maintain the plasma peaks under the conditions of electrical fields or horizontal neutral winds. An assumed influx of H2O molecules of a specific quantity is theorized to reproduce the measured values on Saturn when combined with a vertical plasma drift.

Majeed, Tariq↗
DOE OSTI2019

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↗
DOE OSTI2020

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↗
NASA NTRS1996

Soot and Radiation Measurements in Microgravity Jet Diffusion Flames

The subject of soot formation and radiation heat transfer in microgravity jet diffusion flames is important not only for the understanding of fundamental transport processes involved but also for providing findings relevant to spacecraft fire safety and soot emissions and radiant heat loads of combustors used in air-breathing propulsion systems. Our objectives are to measure and model soot volume fraction, temperature, and radiative heat fluxes in microgravity jet diffusion flames. For this four-year project, we have successfully completed three tasks, which have resulted in new research methodologies and original results. First is the implementation of a thermophoretic soot sampling technique for measuring particle size and aggregate morphology in drop-tower and other reduced gravity experiments. In those laminar flames studied, we found that microgravity soot aggregates typically consist of more primary particles and primary particles are larger in size than those under normal gravity. Comparisons based on data obtained from limited samples show that the soot aggregate's fractal dimension varies within +/- 20% of its typical value of 1.75, with no clear trends between normal and reduced gravity conditions. Second is the development and implementation of a new imaging absorption technique. By properly expanding and spatially-filtering the laser beam to image the flame absorption on a CCD camera and applying numerical smoothing procedures, this technique is capable of measuring instantaneous full-field soot volume fractions. Results from this technique have shown the significant differences in local soot volume fraction, smoking point, and flame shape between normal and reduced gravity flames. We observed that some laminar flames become open-tipped and smoking under microgravity. The third task we completed is the development of a computer program which integrates and couples flame structure, soot formation, and flame radiation analyses together. We found good agreements between model predictions and experimental data for laminar and turbulent flames under both normal and reduced gravity. We have also tested in the laboratory the techniques of rapid-insertion fine-wire thermocouples and emission pyrometry for temperature measurements. These techniques as well as laser Doppler velocimetry and spectral radiative intensity measurement have been proposed to provide valuable data and improve the modeling analyses.

Ku, Jerry C.↗
DOE OSTI2017

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↗