Search NASA⌕ Search

SEARCH · Search NASA

Results for “conditional random field”

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.

432 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Behavior and mechanisms of Doppler wind lidar error in varying stability regimes

Abstract. Wind lidars are widespread and important tools in atmospheric observations. An intrinsic part of lidar measurement error is due to atmospheric variability in the remote-sensing scan volume. This study describes and quantifies the distribution of measurement error due to turbulence in varying atmospheric stability. While the lidar error model is general, we demonstrate the approach using large ensembles of virtual WindCube V2 lidar performing a profiling Doppler-beam-swinging scan in quasi-stationary large-eddy simulations (LESs) of convective and stable boundary layers. Error trends vary with the stability regime, time averaging of results, and observation height. A systematic analysis of the observation error explains dominant mechanisms and supports the findings of the empirical results. Treating the error under a random variable framework allows for informed predictions about the effect of different configurations or conditions on lidar performance. Convective conditions are most prone to large errors (up to 1.5 m s−1 in 1 Hz wind speed in strong convection), driven by the large vertical velocity variances in convective conditions and the high elevation angle of the scanning beams (62∘). Range-gate weighting induces a negative bias into the horizontal wind speeds near the surface shear layer (−0.2 m s−1 in the stable test case). Errors in the horizontal wind speed and direction computed from the wind components are sensitive to the background wind speed but have negligible dependence on the relative orientation of the instrument. Especially during low winds and in the presence of large errors in the horizontal velocity estimates, the reported wind speed is subject to a systematic positive bias (up to 0.4 m s−1 in 1 Hz measurements in strong convection). Vector time-averaged measurements can improve the behavior of the error distributions (reducing the 10 min wind speed error standard deviation to <0.3 m s−1 and the bias to <0.1 m s−1 in strong convection) with a predictable effectiveness related to the number of decorrelated samples in the time window. Hybrid schemes weighting the 10 min scalar- and vector-averaged lidar measurements are shown to be effective at reducing the wind speed biases compared to cup measurements in most of the simulated conditions, with time averages longer than 10 min recommended for best use in some unstable conditions. The approach in decomposing the error mechanisms with the help of the LES flow field could be extended to more complex measurement scenarios and scans.

17 WIND ENERGY↗
DOE OSTI2021

Quantitative Evaluation of Top Coal Caving Methods at the Working Face of Extra‐Thick Coal Seams Based on the Random Medium Theory

Adopting an effective top coal caving method is the key to enhancing coal recovery and reducing gangue content for the fully mechanized top coal caving working face with extra‐thick coal seams. In this study, the movement of coal particles generated during top coal caving is considered to follow a normal distribution. Then, the caving body and coal‐rock settlement along the working face during the caving process are studied based on both the random media theory and probability theory. Accordingly, the optimal caving interval and caving sequences are determined, and a novel interval symmetrical coal caving method is proposed. The proposed method is systematically verified with results from physical similarity tests, and different caving methods are assessed by field tests. The results show the following: (1) The coal‐rock settlement and the caving body demonstrate clear axial symmetrical features along the working face; the size of the caving body increases as the caving height grows and its shape turns progressively from semicircular to semielliptical with a lower foot of the coal‐rock settlement. (2) The caving interval is derived using the sum of the radii of the coal‐rock settlement curves formed by the two largest caving bodies. (3) The symmetrical caving approach provides a symmetrical space for the subsequent movement of the broken top coal, which enables a uniform development of the caving body. (4) Compared with the traditional sequential coal caving method with the same number of supports, the interval symmetrical caving method results in a 21.7% of coal production increase, 17% caving rate promotion, and a shortened caving time by 23.4%. (5) The interval symmetrical caving method is found to improve the controllability of the caving process at the fully mechanized top coal caving working face. In general, this work presents a theoretical approach to select the optimal caving methods for the fully mechanized caving working face in extra‐thick coal seams for an improved production efficiency of the work face. The results of this study can also provide theoretical significance and referencing value for quantitative analyses of the coal caving methods for work faces with similar geological conditions.

Jiulin, Shi (ORCID:0000000155819062)↗
NASA NTRSDate not supplied

The Effects of Chronic Sleep Restriction on Calorie and Macronutrient Intake

Introduction: Chronic sleep restriction (CSR) has been associated with increased calorie intake and increased consumption of fats and carbohydrates, with inconsistent changes in protein. However, the majority of studies have either been observational field studies with no sleep intervention, or laboratory-based studies where food availability may not have reflected participants’ real-world choices. We hypothesized that calorie, fat, and carbohydrate intake would increase during a week of imposed CSR compared to a week of sleep satiation (SS) among individuals living in their home environment. Methods: Twelve healthy participants (6 females) kept a fixed sleep-wake schedule, with a constant waketime, at home for four weeks (actigraphy confirmed compliance). During weeks one and three, participants maintained 9 hours in bed. During weeks two and four, participants were randomly assigned to experimental weeks of 5 and 9 hours of time-in-bed in a crossover design. Participants documented their food consumption during both experimental weeks using a picture-based meal logging application (MealLogger). Intake of calories and macronutrients were classified by two blinded evaluators. Descriptive statistics were calculated in SAS (Cary, NC). Results: Participants averaged 4.43 ± 0.33 (SD) hours of sleep per night during CSR compared to 7.42 ± 0.42 hours during SS. Participants consumed a daily average of 1812 ± 672 kilocalories, 71 ± 31 grams of total fat, 217 ± 69 grams of carbohydrates, and 84 ± 40 grams of protein during CSR, compared to 1682 ± 514 kilocalories, 68 ± 23 grams of total fat, 198 ± 61 grams of carbohydrates, and 77 ± 32 grams of protein during SS. Conclusion: Preliminary descriptive findings suggest that, on average, participants consumed more calories, from an increase in consumption of each macronutrient group, during a week of sleep restriction compared to a week of sleep satiation. Further analysis is needed to determine whether these differences are statistically different and to identify when calories were consumed in each of the experimental conditions.

chronic sleep restriction↗
DOE OSTI2020

Inferred inflow forecast horizons guiding reservoir release decisions across the United States

Abstract. Medium- to long-range forecasts often guide reservoir release decisions to support water management objectives, including mitigating flood and drought risks. While there is a burgeoning field of science targeted at improving forecast products and associated decision support models, data describing how and when forecasts are applied in practice remain undeveloped. This lack of knowledge may prevent hydrological modelers from developing accurate reservoir release schemes for large-scale, distributed hydrology models that are increasingly used to assess the vulnerabilities of large regions to hydrological stress. We address this issue by estimating seasonally varying, regulated inflow forecast horizons used in the operations of more than 300 dams throughout the conterminous United States (CONUS). For each dam, we take actual forward observed inflows (perfect foresight) as a proxy for forecasted flows available to the operator and then identify for each week of the year the forward horizon that best explains the release decisions taken. Resulting “horizon curves” specify for each dam the inferred inflow forecast horizon as a function of the week of the water year. These curves are analyzed for strength of evidence for contribution of medium- to long-range forecasts in decision making. We use random forest classification to estimate that approximately 80 % of large dams and reservoirs in the US (1553±50 out of 1927 dams with at least 10 Mm3 storage capacity) adopt medium- to long-range inflow forecasts to inform release decisions during at least part of the water year. Long-range forecast horizons (more than 6 weeks ahead) are detected in the operations of reservoirs located in high-elevation regions of the western US, where snowpack information likely guides the release. A simulation exercise conducted on four key western US reservoirs indicates that forecast-informed models of reservoir operations may outperform models that neglect the horizon curve – including during flood and drought conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Electrical Analysis of Pulsed Laser Annealed Poly-Si:Ga/SiOx Passivating Contacts

High-efficiency single crystalline silicon (c-Si) solar cells require precise control of dopant diffusion profiles and highly active doping concentrations through thermal annealing. Conventional furnace annealing has been successfully employed for dopant diffusion in aluminum back-surface field (Al-BSF), passivated emitter and rear contact (PERC), and Topcon cells. However, furnace annealing limits some next-generation polycrystalline silicon (poly-Si) on SiOx passivating contacts. This limitation largely affects p-type passivating contacts by not being able to: 1) control the dopant diffusion accurately to prevent B segregation at the SiOx and c-Si interface; and 2) provide highly activated doping concentrations for low contact resistivity. In this contribution, we use a nanosecond excimer laser to examine the passivation quality and electrical performance of both B- and Ga-doped poly-Si/SiOx passivating contacts. The core of this method is to take advantage of the nonequilibrium nature of the anneal through rapid melting and recrystallizing the poly-Si in a short timescale and to achieve doping concentrations above the solid solubility limit. Simulations were performed on polished surfaces to visualize the doping diffusion profiles under different laser conditions, and secondary ion mass spectrometry (SIMS) was used to verify the experimental diffusion profiles post pulsed laser melting (PLM). The results show the dopant diffusion profiles can be tuned precisely through PLM. The electrical analysis using VdP-Hall measurements reveals that the active doping concentration reached 10^21 cm-3 for B and ~2×10^21 cm-3 for Ga in poly-Si, far exceeding their solid solubility limit in Si, with a nearly 100% dopant activation achieved for B, and 20% for Ga. This highly activated dopant profile results in a low contact resistivity of <15 mO·cm^2 for B-doped contacts and <40 mO·cm^2 for Ga-doped contacts. We demonstrate that the passivation quality of the laser annealed contacts shows a linear dependency with the laser energy density and the number of pulses. PLM was also performed on random pyramid textured samples, and we found that preferential melting exists for the tips of pyramids, which can be pursued further for selective pinhole opening. Additionally, we perform thermal stability tests of these hyperdoped passivating contacts and show that under certain laser conditions, the B-doped samples exhibit good thermal stability with near full activation for temperatures up to 800°C. Lastly, we transfer this technique to a larger area by overlapping laser spots and show uniform passivation quality and electrical performance. We are currently implementing this PLM annealed p-type passivating contacts into double-side back-junction poly-Si/SiOx passivating contact devices. Our results show that PLM can not only enable the next generation of high-efficiency c-Si photovoltaics technology but can also benefit the c-Si integrated circuit industry through hyperdoping other semiconductors with atoms that may otherwise have a low solubility.

doping diffusion profile↗
DOE OSTI2015

The shape of LITTLE THINGS dwarf galaxies DDO 46 and DDO 168: Understanding the stellar and gas kinematics

Determining the shape of dwarf irregular (dIrr) galaxies is controversial because if one assumes that these objects are disks and if these disks are randomly distributed over the sky, then their projected minor-to-major axis ratios should follow a particular statistical distribution, which is not observed. Thus, different studies have led to different conclusions. Some believe that the observed distributions can be explained by assuming the dIrrs are thick disks while others have concluded that dIrrs are triaxial. Fortunately, the central stellar velocity dispersion, σ{sub z,0}, combined with maximum rotation speed, V{sub max}, provides a kinematic measure, V{sub max}/σ{sub z,0}, which gives the three-dimensional shape of a system. In this work, we present the stellar and gas kinematics of DDO 46 and DDO 168 from the Local Irregulars That Trace Luminosity Extremes; The H i Nearby Galaxy Survey (LITTLE THINGS) and determine their respective V{sub max}/σ{sub z,0} values. We used the Kitt Peak National Observatory's Mayall 4 m telescope with the Echelle spectrograph as a long-slit spectrograph, which provided a two-dimensional, 3′-long slit. We acquired spectra of DDO 168 along four position angles (PAs) by placing the slit over the morphological major and minor axes and two intermediate PAs. However, due to poor weather conditions during our observing run for DDO 46, we were able to extract only one useful data point from the morphological major axis. We determined a central stellar velocity dispersion perpendicular to the disk, σ{sub z,0}, of 13.5 ± 8 km s{sup −1} for DDO 46 and 〈σ{sub z,0}〉 of 10.7 ± 2.9 km s{sup −1} for DDO 168. We then derived the maximum rotation speed in both galaxies using the LITTLE THINGS H i data. We separated bulk motions from non-circular motions using a double Gaussian decomposition technique and applied a tilted-ring model to the bulk velocity field. We corrected the observed H i rotation speeds for asymmetric drift and found a maximum velocity, V{sub max}, of 77.4 ± 3.7 and 67.4 ± 4.0 for DDO 46 and DDO 168, respectively. Thus, we derived a kinematic measure, V{sub max}/σ{sub z,0}, of 5.7 ± 0.6 for DDO 46 and 6.3 ± 0.3 for DDO 168. Comparing these values to ones determined for spiral galaxies, we find that DDO 46 and DDO 168 have V{sub max}/σ{sub z,0} values indicative of thin disks, which is in contrast to minor-to-major axis ratio studies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Droplet Growth or Evaporation Does Not Buffer the Variability in Supersaturation in Clean Clouds

Abstract Water vapor supersaturation in clouds is a random variable that drives activation and growth of cloud droplets. The Pi Convection–Cloud Chamber generates a turbulent cloud with a microphysical steady state that can be varied from clean to polluted by adjusting the aerosol injection rate. The supersaturation distribution and its moments, e.g., mean and variance, are investigated for varying cloud microphysical conditions. High-speed and collocated Eulerian measurements of temperature and water vapor concentration are combined to obtain the temporally resolved supersaturation distribution. This allows quantification of the contributions of variances and covariances between water vapor and temperature. Results are consistent with expectations for a convection chamber, with strong correlation between water vapor and temperature; departures from ideal behavior can be explained as resulting from dry regions on the warm boundary, analogous to entrainment. The saturation ratio distribution is measured under conditions that show monotonic increase of liquid water content and decrease of mean droplet diameter with increasing aerosol injection rate. The change in liquid water content is proportional to the change in water vapor concentration between no-cloud and cloudy conditions. Variability in the supersaturation remains even after cloud droplets are formed, and no significant buffering is observed. Results are interpreted in terms of a cloud microphysical Damköhler number (Da), under conditions corresponding to , i.e., the slow-microphysics regime. This implies that clouds with very clean regions, such that is satisfied, will experience supersaturation fluctuations without them being buffered by cloud droplet growth. Significance Statement The saturation ratio (humidity) in clouds controls the growth rate and formation of cloud droplets. When air in a turbulent cloud mixes, the humidity varies in space and time throughout the cloud. This is important because it means cloud droplets experience different growth histories, thereby resulting in broader size distributions. It is often assumed that growth and evaporation of cloud droplets buffers out some of the humidity variations. Measuring these variations has been difficult, especially in the field. The purpose of this study is to measure the saturation ratio distribution in clouds with a range of conditions. We measure the in-cloud saturation ratio using a convection cloud chamber with clean to polluted cloud properties. We found in clouds with low concentrations of droplets that the variations in the saturation ratio are not suppressed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Particle-number distribution in large fluctuations at the tip of branching random walks

Here, we investigate properties of the particle distribution near the tip of one-dimensional branching random walks at large times t , focusing on unusual realizations in which the rightmost lead particle is very far ahead of its expected position, but still within a distance smaller than the diffusion radius ~$\sqrt{t}$. Our approach consists in a study of the generating function $G_{Δx}(λ) = Σ_n$ ${λ^n}p_n(Δx)$ for the probabilities $p_n(Δx)$ of observing $\textit{n}$ particles in an interval of given size $Δ\textit{x}$ from the lead particle to its left, fixing the position of the latter. This generating function can be expressed with the help of functions solving the Fisher-Kolmogorov-Petrovsky-Piscounov (FKPP) equation with suitable initial conditions. In the infinite-time and large-$Δ\textit{x}$ limits, we find that the mean number of particles in the interval grows exponentially with $Δ\textit{x}$, and that the generating function obeys a nontrivial scaling law, depending on $Δ\textit{x}$ and λ through the combined variable $[Δx — f(λ)]^3 / Δx^2$, where $\textit{f}$(λ) ≡ – ln(1 – λ) – ln [– ln(1 – λ)]. From this property, one may conjecture that the growth of the typical particle number with the size of the interval is slower than exponential, but, surprisingly enough, only by a subleading factor at large Δ$\textit{x}$. The scaling we argue is consistent with results from a numerical integration of the FKPP equation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1995

Investigation of computational and spectral analysis methods for aeroacoustic wave propagation

Most computational fluid dynamics (CFD) schemes are not adequately accurate for solving aeroacoustics problems, which have wave amplitudes several orders of magnitude smaller yet with frequencies larger than the flow field variations generating the sound. Hence, a computational aeroacoustics (CAA) algorithm should have minimal dispersion and dissipation features. A dispersion relation preserving (DRP) scheme is, therefore, applied to solve the linearized Euler equations in order to simulate the propagation of three types of waves, namely: acoustic, vorticity, and entropy waves. The scheme is derived using an optimization procedure to ensure that the numerical derivatives preserve the wave number and angular frequency of the partial differential equations being discretized. Consequently, simulated waves propagate with the correct wave speeds and exhibit their appropriate properties. A set of radiation and outflow boundary conditions, compatible with the DRP scheme and derived from the asymptotic solutions of the governing equations, are also implemented. Numerical simulations are performed to test the effectiveness of the DRP scheme and its boundary conditions. The computed solutions are shown to agree favorably with the exact solutions. The major restriction appears to be that the dispersion relations can be preserved only for waves with wave lengths longer than four or five spacings. The boundary conditions are found to be transparent to the outgoing disturbances. However, when the disturbance source is placed closer to a boundary, small acoustic reflections start appearing. CAA generates enormous amounts of temporal data which needs to be reduced to understand the physical problem being simulated. Spectral analysis is one approach that helps us in extracting information which often can not be easily interpreted in the time domain. Thus, three different methods for the spectral analysis of numerically generated aeroacoustic data are studied. First, the capabilities of two traditional methods for spectral analysis, namely, the Blackman-Tukey method and periodogram method, are compared in estimating the spectra of a simple-periodic process. The periodogram is then applied to analyze transitory-deterministic processes. Finally, these two methods are compared with a more recent method, referred as the Weighted-Overlapped-Segment-Averaging (WOSA) method, in estimating the spectra of a chaotic (random-like) process. From the demonstrative case for the spectral analyses of data generated by simple-periodic process, the periodogram method is found to give a better estimate of the steep-sloped spectra than the Blackman-Tukey method. Also, for this problem, the Hanning window is found to perform better with the periodogram method than with the Blackman-Tukey method. Finally, for the spectral analysis of data generated by the chaotic process, the periodogram method does not perform well, whereas, the WOSA and Blackman-Tukey methods give equivalently good results.

Vanel, Florence O.↗
DOE OSTI2024

Molecular mechanisms and environmental adaptations of flagellar loss and biofilm growth of Rhodanobacter under environmental stress

Biofilms aid bacterial adhesion to surfaces via direct and indirect mechanisms, and formation of biofilms is considered as an important strategy for adaptation and survival in suboptimal environmental conditions. However, the molecular underpinnings of biofilm formation in subsurface sediment/groundwater ecosystems where microorganisms often experience fluctuations in nutrient input, pH, and nitrate or metal concentrations are underexplored. Here, we examined biofilm formation under different nutrient, pH, metal, and nitrate regimens of 16 Rhodanobacter strains isolated from subsurface groundwater wells spanning diverse levels of pH (3.5 to 5) and nitrates (13.7 to 146 mM). Eight Rhodanobacter strains demonstrated significant biofilm growth under low pH, suggesting adaptations for survival and growth at low pH. Biofilms were intensified under aluminum stress, particularly in strains possessing fewer genetic traits associated with biofilm formation, findings warranting further investigation. Through random barcode transposon-site sequencing (RB-TnSeq), proteomics, use of specific mutants, and transmission electron microscopy analysis, we discovered flagellar loss under aluminum stress, indicating a potential relationship between motility, metal tolerance, and biofilm growth. Comparative genomic analyses revealed the absence of flagella and chemotaxis genes and the presence of a putative type VI secretion system in the highly biofilm-forming strain FW021-MT20. In this study we identified genetic determinants associated with biofilm growth under metal stress in a predominant environmental genus, Rhodanobacter, and identified traits aiding survival and adaptation to contaminated subsurface environments.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

Final Technical Report for Award SC0008613

Discovering how aerosol particles, present in the atmosphere in sizes of a few nanometers to hundred micrometers, initiate ice crystal formation represents a great challenge. Atmospheric ice nucleation is important because ice crystals alter the radiative properties of clouds and thus climate, and impact precipitation and thus the hydrological cycle. The difficulty in predicting atmospheric ice formation is attributable at least in part, to the diversity of ice nucleation pathways, the physical and chemical complexity of the ice nucleating particles (INPs), and the relatively small numbers of INPs (compared with all other aerosol particles), sometimes less than one in 100000. These factors in turn makes constraining ice nucleation parameterizations for modeling applications challenging. The majority of airborne particles are known to be organic in nature or contain organic biogenic material. The presence of organic material adds to the complexity of the particles and therefore the predictability of ice nucleation events since the organic species can display different phase states, e.g. liquid or solid, in response to temperature and humidity. The award DE-SC0008613 to PI Prof. Daniel Knopf at Stony Brook University, “Relating the Chemical and Physical Properties of Aerosols to the Water Uptake and Ice Nucleation Potential of Particles Collected During the Carbonaceous Aerosols and Radiative Effects Study (CARES)”, allowed examination of laboratory generated aerosol particles and field-collected particles for their propensity to nucleate ice under typical tropospheric conditions and relate ice nucleation to the physicochemical properties of the particles including their morphology and chemical composition. This in turn allowed for development of ice nucleation parameterizations for implementation in cloud models. The award resulted in 10 peer-reviewed publications and more than 20 seminar and conference presentations. We demonstrated that the rate of immersion freezing, when a particle immersed in an aqueous droplet acts as an ice nucleus (IN), can be predicted by knowledge of the IN particle type present and the droplet’s water activity which is equal to ambient relative humidity. Our water activity based immersion freezing model is successful in predicting freezing data including INPs such as mineral dusts, marine biological material, organic species, and surfactant molecules. Its mathematical simplicity makes it an ideal candidate for implementation in cloud and climate models. Furthermore, we could show that this model can reproduce many past laboratory measurements which were generated using a variety of instruments to study immersion freezing. Lastly, we have demonstrated that this model can also be applied to field collected particles. Currently, this novel physical parameterization of immersion freezing is being implemented in a cloud model. We examined the physicochemical properties and the ice nucleation potential of particles collected during CARES applying a novel experimental method that allows identification of the individual INPs within a large population of particles sampled from an ambient environment. Taking advantage of a variety of micro-spectroscopic techniques, we characterized the composition and morphology of IN and non-IN particles present in the airborne population. We developed a new parameterization for quantifying the mixing state of the entire aerosol populations by introducing a mixing state index. We found that the identified INPs belong to the most common particle-type classes observed in the CARES field samples and as such are not special or rare particles. In other words, the INPs can be shown not to be unique in contrast to the common paradigm of being rare and exceptional. Either there are differences between particles acting as IN and particles not acting as IN which are beyond our current detection limit or nucleation occurs randomly on the surface of any one of these compositionally equivalent particles. These results suggest that total particle surface area of the different particle types present in the aerosol population is also a crucial factor when predicting ice nucleation in an air mass. We also observed that ambient organic aerosol particles can initiate ice nucleation and corroborated these findings using laboratory generated organic particles. These experiments demonstrated that information of the organic phase state is important for predicting the ice nucleation pathway and the ability of an organic particle to participate in atmospheric ice crystal formation. This award resulted in substantial new insights in the processes governing immersion freezing, the role of organic aerosol particles in ice cloud formation, and the importance of the ambient aerosol population for prediction of ice nucleation in an air parcel. These findings have significant implications for modeling and field measurement strategies of atmospheric ice nucleation.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2013

What Reliability Engineers Should Know about Space Radiation Effects

Space radiation in space systems present unique failure modes and considerations for reliability engineers. Radiation effects is not a one size fits all field. Threat conditions that must be addressed for a given mission depend on the mission orbital profile, the technologies of parts used in critical functions and on application considerations, such as supply voltages, temperature, duty cycle, and redundancy. In general, the threats that must be addressed are of two types-the cumulative degradation mechanisms of total ionizing dose (TID) and displacement damage (DD). and the prompt responses of components to ionizing particles (protons and heavy ions) falling under the heading of single-event effects. Generally degradation mechanisms behave like wear-out mechanisms on any active components in a system: Total Ionizing Dose (TID) and Displacement Damage: (1) TID affects all active devices over time. Devices can fail either because of parametric shifts that prevent the device from fulfilling its application or due to device failures where the device stops functioning altogether. Since this failure mode varies from part to part and lot to lot, lot qualification testing with sufficient statistics is vital. Displacement damage failures are caused by the displacement of semiconductor atoms from their lattice positions. As with TID, failures can be either parametric or catastrophic, although parametric degradation is more common for displacement damage. Lot testing is critical not just to assure proper device fi.mctionality throughout the mission. It can also suggest remediation strategies when a device fails. This paper will look at these effects on a variety of devices in a variety of applications. This paper will look at these effects on a variety of devices in a variety of applications. (2) On the NEAR mission a functional failure was traced to a PIN diode failure caused by TID induced high leakage currents. NEAR was able to recover from the failure by reversing the current of a nearby Thermal Electric Cooler (turning the TEC into a heater). The elevated temperature caused the PIN diode to anneal and the device to recover. It was by lot qualification testing that NEAR knew the diode would recover when annealed. This paper will look at these effects on a variety of devices in a variety of applications. Single Event Effects (SEE): (1) In contrast to TID and displacement damage, Single Event Effects (SEE) resemble random failures. SEE modes can range from changes in device logic (single-event upset, or SEU). temporary disturbances (single-event transient) to catastrophic effects such as the destructive SEE modes, single-event latchup (SEL). single-event gate rupture (SEGR) and single-event burnout (SEB) (2) The consequences of nondestructive SEE modes such as SEU and SET depend critically on their application--and may range from trivial nuisance errors to catastrophic loss of mission. It is critical not just to ensure that potentially susceptible devices are well characterized for their susceptibility, but also to work with design engineers to understand the implications of each error mode. -For destructive SEE, the predominant risk mitigation strategy is to avoid susceptible parts, or if that is not possible. to avoid conditions under which the part may be susceptible. Destructive SEE mechanisms are often not well understood, and testing is slow and expensive, making rate prediction very challenging. (3) Because the consequences of radiation failure and degradation modes depend so critically on the application as well as the component technology, it is essential that radiation, component. design and system engineers work togetherpreferably starting early in the program to ensure critical applications are addressed in time to optimize the probability of mission success.

DiBari, Rebecca↗
NASA NTRS1997

Foreshock Langmuir Waves for Unusually Constant Solar Wind Conditions: Data and Implications for Foreshock Structure

Plasma wave data are compared with ISEE 1's position in the electron foreshock for an interval with unusually constant (but otherwise typical) solar wind magnetic field and plasma characteristics. For this period, temporal variations in the wave characteristics can be confidently separated from sweeping of the spatially varying foreshock back and forth across the spacecraft. The spacecraft's location, particularly the coordinate D(sub f) downstream from the foreshock boundary (often termed DIFF), is calculated by using three shock models and the observed solar wind magnetometer and plasma data. Scatterplots of the wave field versus D(sub f) are used to constrain viable shock models, to investigate the observed scatter in the wave fields at constant D(sub f), and to test the theoretical predictions of linear instability theory. The scatterplots confirm the abrupt onset of the foreshock waves near the upstream boundary, the narrow width in D(sub f) of the region with high fields, and the relatively slow falloff of the fields at large D(sub f), as seen in earlier studies, but with much smaller statistical scatter. The plots also show an offset of the high-field region from the foreshock boundary. It is shown that an adaptive, time-varying shock model with no free parameters, determined by the observed solar wind data and published shock crossings, is viable but that two alternative models are not. Foreshock wave studies can therefore remotely constrain the bow shock's location. The observed scatter in wave field at constant D(sub f) is shown to be real and to correspond to real temporal variations, not to unresolved changes in D(sub f). By comparing the wave data with a linear instability theory based on a published model for the electron beam it is found that the theory can account qualitatively and semiquantitatively for the abrupt onset of the waves near D(sub f) = 0, for the narrow width and offset of the high-field region, and for the decrease in wave intensity with increasing D(sub f). Quantitative differences between observations and theory remain, including large overprediction of the wave fields and the slower than predicted falloff at large D(sub f) of the wave fields. These differences, as well as the unresolved issue of the electron beam speed in the high-field region of the foreshock, are discussed. The intrinsic temporal variability of the wave fields, as well as their overprediction based on homogeneous plasma theory, are indicative of stochastic growth physics, which causes wave growth to be random and varying in sign, rather than secular.

Cairns, Iver H.↗
NASA NTRS2000

Stereo Imaging Velocimetry of Mixing Driven by Buoyancy Induced Flow Fields

Mixing of two fluids generated by steady and particularly g-jitter acceleration is fundamental towards the understanding of transport phenomena in a microgravity environment. We propose to carry out flight and ground-based experiments to quantify flow fields due to g-jitter type of accelerations using Stereo Imaging Velocimetry (SIV), and measure the concentration field using laser fluorescence. The understanding of the effects of g-jitter on transport phenomena is of great practical interest to the microgravity community and impacts the design of experiments for the Space Shuttle as well as the International Space Station. The aim of our proposed research is to provide quantitative data to the community on the effects of g-jitter on flow fields due to mixing induced by buoyancy forces. The fundamental phenomenon of mixing occurs in a broad range of materials processing encompassing the growth of opto-electronic materials and semiconductors, (by directional freezing and physical vapor transport), to solution and protein crystal growth. In materials processing of these systems, crystal homogeneity, which is affected by the solutal field distribution, is one of the major issues. The understanding of fluid mixing driven by buoyancy forces, besides its importance as a topic in fundamental science, can contribute towards the understanding of how solutal fields behave under various body forces. The body forces of interest are steady acceleration and g-jitter acceleration as in a Space Shuttle environment or the International Space Station. Since control of the body force is important, the flight experiment will be carried out on a tunable microgravity vibration isolation mount, which will permit us to precisely input the desired forcing function to simulate a range of body forces. To that end, we propose to design a flight experiment that can only be carried out under microgravity conditions to fully exploit the effects of various body forces on fluid mixing. Recent flight experiments, by the P.I. through collaboration with the Canadian Space Agency (STS-85, August 1997), aimed at determining the stability of the interface between two miscible liquids inside an enclosure show that a long liquid column (5 cm) under microgravity isolation conditions can be stable, i.e. the interface remains sharp and vertical over a short time scale; thus transport occurs by molecular mass diffusion. On the other hand, when the two liquids were excited from a controlled vibration source (Microgravity Vibration Isolation Mount) two to four mode large amplitude quasi-stationary waves were observed. The data was limited to CCD recording of the dynamics of the interface between the two fluids. We propose to carry out flight experiments to quantify the dynamics of the flow field using Stereo Imaging Velocimetry and measure the concentration field using laser fluorescence. The results will serve as a basis to understand effects of g-jitter on transport phenomena, in this case mass diffusion. As the measurement of the kinematics of the flow field will shed light on the instability mechanism. The research will allow measurement of the flow field in microgravity environment to prove two hypotheses: (1) Maxwell's hypothesis: finite convection always exists in diffusing systems, and (2) Quasi-stationary waves inside a bounded enclosure in a microgravity environment is generated by Kelvin-Helmholtz instability; resonance of the interface which produces incipient mixing is due to Rayleigh-Taylor instability. The first hypothesis can be used as a benchmark experiment to illustrate diffusive mixing. The second hypothesis will lead to the understanding of g-jitter effects on buoyancy driven flow fields which occur in many situations involving materials processing, and other basic fluid physics phenomena. In addition, the second hypothesis will also provide insight in how Rayleigh-Taylor and Kelvin-Helmholtz instabilities propagate concentration fronts during mixing. Measurement of the flow field using SIV is important because it is the flow field which causes instability at the interface between the two fluids. Mixing driven by buoyancy induced flow fields will be addressed both experimentally and computationally. The experimental effort will address the kinematics of mixing: stretching, transport and chaos. Quantification of the mechanisms of mixing will consists of measuring the flow field using the SIV system at Glenn and capturing the dynamics of the interface, to measure mass transport, using a CCD camera. These experiments will be carried out within the framework of Earth's gravity and g-jitter microgravity acceleration as in a Space Shuttle environment or the International Space Station. The g-jitter will be induced and controlled using a tunable vibration isolation platform to isolate against vibration as well as input periodic and random vibration to the system. The parametric range of the microgravity experiment will be extended from the experiments on STS-85 to investigate higher mode quasi-stationary waves (8 to 12), as well as resonance regions which leads to chaos and turbulence. Ground-based experiments will focus on effects of vibration on stably stratified fluid layers in order to scale for possible scenarios in a microgravity environment. These vibrations will be subjected perpendicular to the concentration field on the ground since the parallel case can only be carried out in a microgravity environment. The concept of dynamical similarity will be applied to tune the experiments as closely as possible to a Space Shuttle environment or the International Space Station. The computational effort will take advantage of the Computational Laboratory at Glenn to corroborate the experimental findings with predictions of the dynamics of the flow field using the codes FLUENT (finite difference based) and FIDAP (finite element based). We will investigate two important cases, single-fluid model to address dilute systems with negligible jump in viscosity and the more general two-fluid model which accounts for finite jump in viscosity. Apart from its microgravity relevance, this experiment is well suited to study dynamics in nonlinear systems.

Duval, W. M. B.↗
DOE OSTI2020

Dynamic Behavior of Natural Seep Vents: Analysis of Field and Laboratory Observations and Modeling (Final Scientific/Technical Report)

In this project, we have analyzed data collected by the U.S. Department of Energy (DOE), National Energy Technology Laboratory (NETL) in a high pressure water tunnel (HPWT) and data from two research cruises to natural seeps in the Gulf of Mexico to adapt and validate a numerical model to predict the dynamics of natural seeps in the deep oceans. The HPWT data include video observations of the shrinkage rate of individual methane and natural gas bubbles under simulated deep-water conditions. Field data were collected during two cruises by the Gulf Integrated Spill Research (GISR) Consortium led by Texas A&M University and funded by the Gulf of Mexico Research Initiative (GoMRI). These data included in situ observations from a remotely operated vehicle (ROV) of gas bubbles at two natural seep sites in the Gulf and acoustic observations of the natural seep bubble flares in the ocean water column. The acoustic data were from multibeam echosounders, one mounted in a forward-looking orientation on the ROV and another mounted down-looking in the haul of the ship. All of these laboratory and field data were focused on the dynamics of natural gas bubbles at temperatures and pressures favorable for clathrate hydrate formation between the gas and water. Our analyses of this data focused on understanding the mechanisms responsible for gas bubble dissolution within the hydrate stability zone (HSZ) of the oceans. Ice-like hydrate shells may form on the bubble-water interface under these conditions, and it was unknown how this might affect the mass transfer of gas into the ocean. We were able to extract bubble shrinkage rates from the HPWT datasets. Using this data we determined that mass transfer coefficients with and without a hydrate shell match empirical values for bubbles in contaminated systems (so-called dirty bubbles contaminated by naturally occurring surfactants). We also showed that free gas, and not gas hydrate, is the dominant dissolving phase when the hydrate sub-cooling is below 11 degree Celsius (temperature difference between hydrate the hydrate formation temperature and ambient temperator) or the pressure is reducing as bubbles rise through the ocean water column. Using this mass transfer model, our numerical model of bubble dissolution matched the over 200 HPWT experiments with an average error of 10% for predicting the bubble size at the end of an experiment. From field data in the literature, we also observed that gas bubbles dissolve faster when they are initially released, following mass transfer coefficients for so-called clean-bubbles (those not yet contaminated by surfactants). Shortly after release within the HSZ, a hydrate shell forms on the bubble-water interface, and the mass transfer reduces to rates matching those of dirty bubbles. We correlated this transition time from clean to dirty bubble behavior with the initial bubble surface area and the hydrate sub-cooling. With this model for hydrate formation time and using the mass transfer coefficients deduced from the HPWT data, we validated our numerical model for predicting the rise heights of natural seep flares in the oceans. Flare heights are commonly observed in haul-mounted acoustic multibeam data. The numerical model predicts bubbles to rise high in the ocean water column owing to the slower mass transfer rates for dirty bubbles that accompany the majority of their rise time. We found that the numerical model predictions matched the observed flare heights within 5% to 10% accuracy when we compared the rise heights of the largest bubbles released from the seafloor with the bubbles acoustically visible in the multibeam data. Bubbles become acoustically transparent as they shrink to sizes of order 1 mm in diameter for the multibeam frequencies used in the field. The forward-looking multibeam on the ROV also provided data on the lateral spreading of bubbles in natural seep flares. Our analysis of this data showed that spreading follows a diffusion process, with the effective diffusivity correlating with the wobbling length scale of these ellipsoidal bubbles. When we apply this diffusivity in a random displacement model of bubble spreading, our numerical simulations match closely the lateral spread observed by the M3 in the ocean water column. Finally, we compared the seep model predictions for the acoustic properties of these natural seep plumes with that observed by the acoustic instruments in the field. The M3 and EM 302 observations were converted to relative values of target strength using a calibration we obtained in the laboratory for the M3 and using an algorithm from the manufacturer for the EM 302. Comparing the numerical seep model to these data, we obtain good agreement over the whole height of rise of these bubble flares. This further validates the numerical model. Overall, our validated seep model captures the key dynamics of gas bubbles released from natural seeps in the oceans and helps to predict the fate of methane in the water column.

03 NATURAL GAS↗
DOE OSTI2016

Comfort in High-Performance Homes in a Hot-Humid Climate

TO5 2.3.2: Large-Scale Temperature Study Although space-conditioning systems are required to meet standards such as those in the Air Conditioning Contractors of America (ACCA) Manual RS (Rutkowski 1997) and the ASHRAE Standard 55 (ASHRAE 2013) comfort criteria, little is known about how space-conditioning systems are actually operating in the field and if these systems are meeting those standards. Unconventional space-conditioning systems may be necessary for low-load homes in the future. It is important to know how current systems are functioning in existing homes and to what temperatures occupants typically set the thermostats in their homes. The U.S. Department of Energy's Building America research team IBACOS set out to find some answers. By using low-cost data loggers and occupant surveys to monitor 37 single-family homes that have standard heating, ventilating, and air-conditioning (HVAC) equipment, the team gained insight into how sensitive occupants are to comfort in their homes relative to the ACCA and ASHRAE standards. This study focused on one climate zone (hot-humid); however, future studies may find differences among occupants- ideas of thermal comfort in other climate zones. Data were collected from installed sensors during a 2-month period from late August 2014 through late October 2014. The final aspect of this project was to review the business metrics associated with builders involved in a comfort and performance guarantee program. The results of this study-taken as an aggregate-show that the homes provided room-to-room temperature differences less than 6 degrees F 95% of the time. Temperature differences were less than 4 degrees F 80% of the time. Some homes showed better or worse performance. On average, two-story homes had a 3.3 degrees F temperature difference between rooms; single-story homes had an average of a 2.2 degrees F temperature difference between rooms. Occupants who did not use programmable thermostats had an average set point of 75 degrees F. Homes in which the thermostats were programmed with a setback schedule showed a median baseline temperature value of 74 degrees F and had a median afternoon setback of 75 degrees F. System runtime was analyzed and showed that the room-to-room temperature uniformity in some homes worsened during an on cycle, whereas the uniformity improved in other homes. That is, when the system was running, some rooms were receiving too much or too little air relative to the rest of the house, and the temperatures in those rooms were moving away from the temperatures of the other rooms in the house. Finally, analysis of the thermostat data shows a wide degree of variability among homes in the study. Interpretation of the measured data suggests that 21 occupants made no regular adjustments to the thermostats in their homes, seven had regular setbacks, seven showed varying setbacks, and 10 appeared to follow no pattern (with random adjustments). Ultimately, the data show that current systems are maintaining expected levels of comfort. As expected, it is more challenging to maintain a uniform temperature in all rooms of two-story homes, and future systems should address this issue. Furthermore, the data show that the majority of homes maintained relative humidity levels below 60% without any supplemental dehumidification

ACCA Manual RS↗
DOE OSTI2026

TCAD-Machine Learning Enabled TID Compact Model Development for Commercial SiC MOSFET

We propose a TCAD (Technology Computer Aided Design)-machine learning coupled approach that combines a TCAD tool (Charon), optimization/uncertainty quantification tool (Dakota), surrogate models, and Bayesian learning capabilities. The coupling approach is used for accurate modeling and calibration of total ionizing dose (TID) induced threshold voltage (V th ) shifts in Commercial-Off-The-Shelf (COTS) semiconductor devices and to develop physics-informed TID compact models. This versatile approach is applied to model the TID effect in an exemplar COTS 3.3 kV SiC power MOSFET (Metal-Oxide-Semiconductor Field-Effect Transistor). With the Charon-Dakota coupling, we can determine key device geometry and doping values based on device physics, which are difficult to obtain or not available for COTS devices but important for TCAD simulation; additionally, we can efficiently generate thousands of simulation results in a large parameter space, which makes it possible to develop data-driven surrogate models and perform Bayesian calibration. Utilizing the full tool-coupling approach, we achieve calibrated TCAD simulation models that accurately capture the average TID-induced V th shifts behavior with total doses and V th shifts saturation at high doses as observed in experimental data. More importantly, the calibrated TCAD simulations are obtained with determined TID model parameters (e.g., hole trap density and capture cross section) values that contain well quantified uncertainties. Furthermore, we can isolate and quantify the noises that are not captured by the TCAD models but exist in the measured data due to measurements and devices variabilities. Lastly, the calibrated surrogate models are used to develop physics-informed TID compact models. The method is generalizable to other devices and/or radiation conditions with few modifications and can provide well-determined uncertainties.

COTS↗
DOE OSTI2024

Effect of Heat Treatment on Microstructure and Mechanical Property of 316L Stainless Steel Produced by Laser Powder Bed Fusion

The advanced non-light water reactor designs (Gen IV reactors), including molten salt/ very high temperature/ sodium-cooled and lead-cooled fast reactors, typically operate at higher temperatures and more extreme radiation conditions than light water reactors. An intrinsic part of the deployment and progress of Gen IV reactor designs is selecting the most suitable structural material for a specific application. Additive manufacturing (AM), a fairly new process of making physical, three-dimensional objects from a computer design file, is going to completely change the way of design, build and certify nuclear systems. It offers a range of opportunities to produce complex geometries from existing materials, offers new routes for processing of previously difficult to process materials, allows for design of new high-performance materials, and finally facilitates hybridization of dissimilar materials. This emerging technology has successfully produced cars, wind turbine blade molds and even live cells. It could also open up big opportunities for the nuclear industry to quickly deploy technologies at a fraction of the cost. So far, AM techniques have been preliminarily applied in the field of nuclear reactors, including the classical parts such as the pressure vessel of a small reactor with 508-III steel, the bottom nozzle of a fuel assembly with 304L steel, the fuel cladding with zirconium alloy and the integrated impeller of a pump and the multi-channel valve body with 316L steel [6,7]. The AM applications for operating nuclear reactors started in auxiliary plant components and have slowly migrated to metallic reactors and core components, but many of these are not safety critical components. Although many parts used for nuclear reactors have been fabricated by AM techniques, practical applications in engineering are still a long way off due to the uncertainty factors focused on the processing, material properties, analysis methods and application standards, which feeds the safety and life-cycle of the nuclear reactor. Due to rapid, repeated heating and cooling during production, a high dislocation density was present in the AM material. This microstructure feature is unstable at elevated temperature while high temperature is one of the typical operation environments for nuclear reactors. Thus, it is important to understand the thermal effect on the microstructure of AM material. The objectives of this study are to investigate the effect of heat treatment on the microstructure and mechanical properties of 316L stainless steel produced by laser powder bed fusion additive manufacturing, and to determine an appropriate heat treatment practice that will be applied to the lightweight AM lattice-structured material with the same chemistry. The heat treatment study consisted of annealing the samples at a temperature range of 800 to 1200 oC with a 50 oC increment for different times (1-24 hours), followed by vacuum or air cooling. Microstructural characterization was carried out by Scanning Electron Microscope (SEM). Grain size and crystallographic orientation were investigated by Electron Backscatter Diffraction (EBSD). Vickers hardness tests with a 0.5 kg load were employed to determine the hardness of samples after different heat treatments. After heat treatment, the random crystallographic orientation was preserved, and the volume fraction of high-angle grain boundaries (grain boundary misorientation =15 oC) remained the same. The dislocation density decreased with annealing temperature due to recovery. The fine subgrain structures in the as-printed specimen were quite stable up to 1200 oC. Minimal recrystallization was observed up to 1200 oC. Recrystallization initiated only after 8.5 hours at 1200 oC. The SEM images did not show obvious dependence of microstructure on cooling rate. The hardness of the specimens decreased with increasing annealing temperature as a result of the decrease in dislocation density. It is interesting to note that the AM material showed very similar hardness to the wrought material when annealing at similar temperature, although the microstructures are very different. Annealing at 1050 oC for 1 hour followed by air cooling was selected as the heat treatment procedure for the lattice designed lightweight AM 316L material.

36 MATERIALS SCIENCE↗