Engineering Papers⌕ Search

SEARCH · Engineering Papers

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.

At least 451 records · Page 25

A Heterogeneous System for Eagle Detection, Deterrent, and Wildlife Collision Detection for Wind Turbines (Final Technical Report)

This report summarizes the design, implementation, and test of an integrated system for automated detection and deterrence of eagles, with included wind turbine blade strike detection and imaging functionality. A machine learning approach was used in conjunction with a 360° camera system for automated detection and classification of golden eagles. This was developed using footage obtained from trained golden eagles and other raptors, in collaboration with wildlife biologists and professional bird handlers. Oregon State University developed a visual deterrent system, which uses inflatable anthropomorphic sculptures with random, kinetic motion to deter eagles, and conducted limited field testing on live eagles; the deterrent can be triggered by the visual detection of eagles using the vision system. Finally, a multi-sensor module was developed that is mounted at the turbine blade root. This module measures vibration and other motions to detect blade strikes, and an integrated on-blade camera captures an image of any impacting objects. Long-term, this blade strike detection system is intended to support an automatic monitoring and certification system for the eagle detection and deterent system. Independent field testing of each system component is described. Testing of the integrated system on an operational wind turbine was conducted across three separate field tests. This includes multi-day fields tests on a General Electric 1.5MW wind turbine at the National Renewable Energy Laboratory (NREL) National Wind Technology Center (NWTC) in Boulder, CO in October 2018 and July 2019; installation procedures, test procedures, and a summary of collected data are presented. A third multi-day on-turbine field test is also presented, which was performed using a General Electric 1.5MW wind turbine at the North American Wind Research and Training Center (NAWRTC) at Mesalands Community College, Tucumcari, NM in April 2019. Across these field tests, the vision system was demonstrated using unmanned aerial vehicles (UAV), and the eagle classification algorithm was not tested; the visual deterrent system was demonstrated, including automatic, remote deployment following surrogate visual detections; and, multi-sensor on-blade data was recorded across multiple wind turbine operational conditions and through more than 100 surrogate blade strikes using soft projectiles, including the successful demonstration of automatic image capture of striking objects. This data set was also used for offline development and validation of enhanced collision detection algorithms. As summarized in this report, the development and field validation of an integrated detection, deterrent, and blade collision detection system represents a critical proof of concept for future technology development of related detection and deterrent technologies, where both deterrent as well as collision detection recording devices are needed for future siting, monitoring, and operation of wind turbine installations, both onshore and offshore.

17 WIND ENERGY↗

Evolution of the electrothermal instability from thick rod z pinches subject to dynamically and statically applied axial magnetic field

LDRD Project 229427 aimed to determine how electrothermal instability (ETI) driven heating on a z-pinch rod pulsed with intense current evolves under mixed magnetic field (azimuthal + axial) conditions, which is pertinent to pulsed-power-driven magnetically-insulated transmission lines and physics targets. Experiments focused on diagnosing ETI-driven heating from deliberately-machined and well-characterized micron-scale surface defects (referred to as engineered defects or ED). Prior to the start of this project, understanding of how unmagnetized (B z =0) ED evolve had been obtained—simulations largely reproduce the experimentally observed high temperature spots which develop at the poles of bare/uncoated ED. Project 229427 extended the Mykonos Facility ED experimental platform to include axial field. In the first class of experiments, axial field was provided “dynamically” via a helical return can (HRC). In this case, B z and B θ rise at the same rate. Generally, the HRC generated magnetic field at a fixed polarization angle Φ B =arctan(B z /B θ )=15° on the rod's surface. In the second class of experiments, axial field was provided “statically” via a slow-rising (millisecond) external Helmholtz coil pair. In this case, B z was effectively constant/static throughout the 100 ns rise of the Mykonos current. For either case, a primary goal was to determine whether ETI provides a helical seed perturbation for the subsequent growth of the helical magneto Rayleigh-Taylor modes observed in MagLIF (static B z ) and dynamic screw pinch (DSP, dynamic B z ) experiments. When dynamic field was applied using an HRC, emissions from individual ED aligned toward Φ B , while emissions from ED within pairs elongated and preferentially merged along Φ B . These data strongly support that for a randomized defect distribution, heating from nearby current-density perturbations will favorably merge about Φ B to generate an extended seed perturbation that aligns toward the surface-field polarization, and this may impact the orientation of subsequent MRT growth on imploding liners. The results from the static field experiments were largely inconclusive, as any ETI heating rotation, if present, was obscured/overwhelmed by local/random heating from ED rim imperfections.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Deflation as a method of variance reduction for estimating the trace of a matrix inverse

Many fields require computing the trace of the inverse of a large, sparse matrix. The typical method used for such computations is the Hutchinson method which is a Monte Carlo (MC) averaging over matrix quadratures. To improve its convergence, several variance reductions techniques have been proposed. In this paper, we study the effects of deflating the near null singular value space. We make two main contributions. First, we analyze the variance of the Hutchinson method as a function of the deflated singular values and vectors. Although this provides good intuition in general, by assuming additionally that the singular vectors are random unitary matrices, we arrive at concise formulas for the deflated variance that include only the variance and mean of the singular values. We make the remarkable observation that deflation may increase variance for Hermitian matrices but not for non-Hermitian ones. This is a rare, if not unique, property where non-Hermitian matrices outperform Hermitian ones. The theory can be used as a model for predicting the benefits of deflation. Second, we use deflation in the context of a large scale application of "disconnected diagrams" in Lattice QCD. On lattices, Hierarchical Probing (HP) has previously provided an order of magnitude of variance reduction over MC by removing "error" from neighboring nodes of increasing distance in the lattice. Although deflation used directly on MC yields a limited improvement of 30% in our problem, when combined with HP they reduce variance by a factor of over 150 compared to MC. For this, we pre-computated 1000 smallest singular values of an ill-conditioned matrix of size 25 million. Furthermore, using PRIMME and a domain-specific Algebraic Multigrid preconditioner, we perform one of the largest eigenvalue computations in Lattice QCD at a fraction of the cost of our trace computation.

97 MATHEMATICS AND COMPUTING↗

Information Leakage Analysis using a Co-design-Based Fault Injection Technique on a RISC-V Microprocessor

The RISC-V instruction set architecture open licensing policy has spawned a hive of development activity, making a range of implementations publicly available. The environments in which RISC-V operates have expanded correspondingly, driving the need for a generalized approach to evaluating the reliability of RISC-V implementations under adverse operating conditions or after normal wear-out periods. Fault injection (FI) refers to the process of changing the state of registers or wires, either permanently or momentarily, and then observing execution behavior. The analysis provides insight into the development of countermeasures that protect against the leakage or corruption of sensitive information which might occur because of unexpected execution behavior. In this paper, we develop a hardware-software co-design architecture that enables fast, configurable fault emulation and utilize it for information leakage and data corruption analysis. Modern System-on-chip FPGAs enable building an evaluation platform where control elements run on a processor(s) (PS) simultaneously with the target design running in the programmable logic (PL). Software components of the FI system introduce faults and report execution behavior. A pair of RISC-V FI-instrumented implementations are created and configured to execute the Advanced Encryption Standard and Twister algorithms. Key and plaintext information leakage and degraded pseudo-random sequences are both observed in the output for a subset of the emulated faults.

42 ENGINEERING↗

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↗

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)↗

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↗

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↗

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↗

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↗

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↗

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↗

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.↗

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↗

Realistic multisite lattice-gas modeling and KMC simulation of catalytic surface reactions: Kinetics and multiscale spatial behavior for CO-oxidation on metal (100) surfaces

A realistic molecular-level description of catalytic reactions on single-crystal metal surfaces can be provided by stochastic multisite lattice-gas (msLG) models. This approach has general applicability, although in this report, we will focus on the example of CO-oxidation on the unreconstructed fcc metal (100) or M(100) surfaces of common catalyst metals M = Pd, Rh, Pt and Ir (i.e., avoiding regimes where Pt and Ir reconstruct). These models can capture the thermodynamics and kinetics of adsorbed layers for the individual reactants species, such as CO/M(100) and O/M(100), as well as the interaction and reaction between different reactant species in mixed adlayers, such as (CO + O)/M(100). The msLG models allow population of any of hollow, bridge, and top sites. This enables a more flexible and realistic description of adsorption and adlayer ordering, as well as of reaction configurations and configuration-dependent barriers. Adspecies adsorption and interaction energies, as well as barriers for various processes, constitute key model input. The choice of these energies is guided by experimental observations, as well as by extensive Density Functional Theory analysis. Model behavior is assessed via Kinetic Monte Carlo (KMC) simulation. We also address the simulation challenges and theoretical ramifications associated with very rapid diffusion and local equilibration of reactant adspecies such as CO. These msLG models are applied to describe adsorption, ordering, and temperature programmed desorption (TPD) for individual CO/M(100) and O/M(100) reactant adlayers. In addition, they are also applied to predict mixed (CO + O)/M(100) adlayer structure on the nanoscale, the complete bifurcation diagram for reactive steady-states under continuous flow conditions, temperature programmed reaction (TPR) spectra, and titration reactions for the CO-oxidation reaction. Extensive and reasonably successful comparison of model predictions is made with experimental data. Furthermore, we discuss the possible transition from traditional mean-field-type bistability and reaction kinetics for lower-pressure to multistability and enhanced fluctuation effects for moderate- or higher-pressure. Behavior in the latter regime reflects a stronger influence of adspecies interactions and also lower diffusivity in the higher-coverage mixed adlayer. We also analyze mesoscale spatiotemporal behavior including the propagation of reaction diffusion fronts between bistable reactive and inactive states, and associated nucleation-mediated transitions between these states. This behavior is controlled by complex surface mass transport processes, specifically chemical diffusion in mixed reactant adlayers for which we provide a precise theoretical formulation. The msLG models together with an appropriate treatment of chemical diffusivity enable equation-free heterogeneous coupled lattice-gas (HCLG) simulations of spatiotemporal behavior. In addition, msLG + HCLG modeling can describe coverage variations across polycrystalline catalysts surfaces, pressure variations across catalyst surfaces in microreactors, and could be incorporated into a multiphysics framework to describe mass and heat transfer limitations for high-pressure catalysis. (C) 2013 Elsevier Ltd. All rights reserved.

36 MATERIALS SCIENCE↗

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↗

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↗

Environmental Impacts of Wind Power Development on the Population Biology of Greater Prairie-Chickens

1. We investigated the impacts of wind power development on the demography, movements, and population genetics of Greater Prairie-Chickens (Tympanuchus cupido) at three sites in northcentral and eastern Kansas for a 7-year period. Only 1 of 3 sites was developed for wind power, the 201MW Meridan Way Wind Power Facility at the Smoky Hills site in northcentral Kansas. Our project report is based on population data for prairie chickens collected during a 2-year preconstruction period (2007-2008), a 3-year postconstruction period (2009-2011) and one final year of lek surveys (2012). Where relevant, we present preconstruction data from our field studies at reference sites in the northern Flint Hills (2007-2009) and southern Flint Hills (2006-2008). 2. We addressed seven potential impacts of wind power development on prairie chickens: lek attendance, mating behavior, use of breeding habitat, fecundity rates, natal dispersal, survival rates, and population numbers. Our analyses of pre- and postconstruction impacts are based on an analysis of covariance design where we modeled population performance as a function of treatment period, distance to eventual or actual site of the nearest wind turbine, and the interaction of these factors. Our demographic and movement data from the 6-year study period at the Smoky Hills site included 23 lek sites, 251 radio-marked females monitored for 287 bird-years, and 264 nesting attempts. Our genetic data were based on genotypes of 1,760 females, males and chicks that were screened with a set of 27 microsatellite markers that were optimized in the lab. 3. In our analyses of lek attendance, the annual probability of lek persistence during the preconstruction period was ~0.9. During the postconstruction period, distance to nearest turbine did not have a significant effect on the probability of lek persistence. However, the probability of lek persistence increased from 0.69 at 0 m to 0.89 at 30 km from turbines, and most abandoned lek sites were located <5 km from turbines. Probability of lek persistence was significantly related to habitat and number of males. Leks had a higher probability of persistence in grasslands than agricultural fields, and increased from ~0.2 for leks of 5 males, to >0.9 for leks of 10 or more males. Large leks in grasslands should be a higher priority for conservation. Overall, wind power development had a weak effect on the annual probability of lek persistence. 3. We used molecular methods to investigate the mating behavior of prairie chickens. The prevailing view for lek-mating grouse is that females mate once to fertilize the clutch and that conspecific nest parasitism is rare. We found evidence that females mate multiple times to fertilize the clutch (8-18% of broods, 4-38% of chicks) and will parasitize nests of other females during egg-laying (~17% of nests). Variable rates of parentage were highest in the fragmented landscapes at the Smoky Hills field site, and were lower at the Flint Hills field site. Comparisons of the pre- and postconstruction periods showed that wind energy development did not affect the mating behaviors of prairie chickens. 4. We examined use of breeding habitats by radio-marked females and conducted separate analyses for nest site selection, and movements of females not attending nests or broods. The landscape was a mix of native prairie and agricultural habitats, and nest site selection was not random because females preferred to nest in grasslands. Nests tended to be closer to turbines during the postconstruction period and there was no evidence of behavioral avoidance of turbines by females during nest site selection. Movements of females not attending nests or broods showed that females crossed the site of the wind power development at higher rates during the preconstruction period (20%) than the postconstruction period (11%), and that movements away from turbines were more frequent during the postconstruction period. Thus, wind power development appears to affect movements in breeding habitats but not nest site selection of female prairie chickens during the breeding season. 5. We tested the effects of wind power development on five components of female fecundity: timing of clutch initiation, clutch size of first nests and renests, nest survival, and hatchability of eggs. Average date of clutch initiation was 26 April, clutch size was 12.7 and 10.6 eggs for first nests and renests, probability of nest survival was low at 0.18, but egg hatchability was high at 0.79. Wind power development had no impact on reproductive effort or nesting success, and all five components of fecundity were not related to treatment period or distance to turbine. Nest survival was the main factor limiting reproductive output of female prairie chickens and most losses were due to predation. Daily nest survival was strongly related to vegetative cover at the nest. Changes to rangeland management practices that would double nesting cover from 2.5 to 5 dm would triple the probability of nest survival from 0.17 to 0.52. Grass and forb cover had weak positive effects on daily nest survival whereas shrub cover, proximity to woodlands, and recent rainfall had negative effects. Reproductive performance of prairie chickens is low in managed rangelands in northcentral Kansas and efforts to improve range conditions and reduce predator activity would aid recovery of prairie chicken populations. 6. We used molecular methods to investigate patterns of natal dispersal in prairie chickens. High rates of nest failure limited the number of young that we could sample. Direct detections of natal dispersal were limited because survival of newly hatched chicks to become adults were low and because we were unable to detect dispersal distances outside of our study area. Direct observations of natal movements were limited and were inadequate to make conclusions about the potential impacts of wind energy development on natal dispersal. Spatial correlograms of genetic distance among males at leks were a more sensitive measure of population structure, and indicated a weak effect of wind energy development on the spatial genetic structure of prairie chickens. 7. We tested the effects of wind power development on female survival with time-to-event models, and on residual body mass of males with analyses of covariance. Distance to turbine and the interaction of distance and treatment period had no effect on female survival. Contrary to predictions of negative impacts of wind power development, the probability of female was lowest during the preconstruction period (0.274) and increased significantly during the postconstruction period (0.543). Inspection of hazard functions indicated that the difference in annual survival could be attributed to a higher risk of mortality during the lekking season in the preconstruction period. We suggest that wind power development may have improved ecological conditions for prairie chickens by disrupting the foraging behavior of diurnal raptors that kill prairie chickens at lek sites. In support of this idea, raptor kills tended to be farther from turbines during the postconstruction period whereas mammalian kills were closer. Analyses of the major causes of mortality did not support our hypothesis because the odds of raptor predation were greater after development was completed. Most mortality losses of radio-marked females were due to predation, and losses to collision mortality or harvest were rare events. Low rates of natural mortality during fall and winter imply that harvest is likely to be additive mortality in prairie chickens. Wind energy development reduced the residual body mass of male Greater Prairie-Chickens at lek sites near turbines. Low values of residual body mass could have a negative impact on individual survival or fecundity rates, or may be related to predation risk and flight performance of males displaying at open lek sites. 8. We tested for impacts of wind power development on population numbers of prairie chickens with monitoring of male numbers at leks, and with genetic measures of population structure. Lek counts indicated that wind power development did not affect the population size of prairie chickens. Peak counts of males at leks were recorded the first year after construction was completed and the highest rates of population change were observed during the interval when the wind power facility was constructed. Population numbers of prairie chickens near and distant from turbines appeared to covary in parallel, probably because bird numbers were linked by dispersal movements. Estimates of population viability based on genetic diversity, effective population size and rates of population exchange did not show annual changes and were unaffected by wind development during our study. Estimates of relatedness among males at the same and different leks suggested that wind power development has either reduced dispersal rates or changed settlement patterns, leading to higher rates of relatedness among males displaying at the same lek site. 9. Greater Prairie-Chickens were not strongly affected by wind power development in Kansas. Negative impacts of wind power development included a trend for reductions in lek persistence near turbines, behavioral avoidance of turbines by females during their breeding season movements, and changes in the genetic structure of males at leks that were consistent with reduced dispersal or recruitment rates. We found no impacts of wind power development on nest site selection, female reproductive effort or nesting success, or population numbers. Positive impacts of wind power development included an increase in female survival rates. We hypothesized that the unexpected increase in female survival was related to changes in trophic interactions and disruption of the foraging behavior of raptors that kill prairie chickens at lek site. 10. Research funding for this project included a grant from the 20% Wind by 2030 Program of the Department of Energy (this final report), and grants from the Kansas Department of Wildlife, Parks and Tourism, the National Fish and Wildlife Federation, and initial funding from the National Wind Coordinating Collaborative. Research products from data collected during the preconstruction period have included five peer-reviewed research articles and two PhD dissertations at Kansas State University. Additional manuscripts are in review for possible publication in 2013-2014.

17 WIND ENERGY↗