The operation of Wave Energy Converter (WEC) devices can pose many challenging problems to the Water Power Community. A key research question is how to significantly improve the performance of these WEC devices through improving the control system design. This report summarizes an effort to analyze and improve the performance of WEC through the design and implementation of control systems. Controllers were selected to span the WEC control design space with the aim of building a more comprehensive understanding of different controller capabilities and requirements. To design and evaluate these control strategies, a model scale test-bed WEC was designed for both numerical and experimental testing (see Section 1.1). Seven control strategies have been developed and applied on a numerical model of the selected WEC. This model is capable of performing at a range of levels, spanning from a fully-linear realization to varying levels of nonlinearity. The details of this model and its ongoing development are described in Section 1.2.
Ruehl, Kelley M.; Chang, G.; Jones, C.A.; Roberts, J.; Chartrand, C.
Modeled nearshore wave propagation was investigated downstream of simulated wave energy converters (WECs) to evaluate overall near- and far-field effects of WEC arrays. Model sensitivity to WEC characteristics and WEC array deployment scenarios was evaluated using a modified version of an industry standard wave model, Simulating WAves Nearshore (SWAN), which allows the incorporation of device-specific WEC characteristics to specify obstacle transmission. The sensitivity study illustrated that WEC device type and subsequently its size directly resulted in wave height variations in the lee of the WEC array. Wave heights decreased up to 30% between modeled scenarios with and without WECs for large arrays (100 devices) of relatively sizable devices (26 m in diameter) with peak power generation near to the modeled incident wave height. Other WEC types resulted in less than 15% differences in modeled wave height with and without WECs, with lesser influence for WECs less than 10 m in diameter. Wave directions and periods were largely insensitive to changes in parameters. However, additional model parameterization and analysis are required to fully explore the model sensitivity of peak wave period and mean wave direction to the varying of the parameters.
In the event of contamination in a water distribution network (WDN), source identification (SI) methods that analyze sensor data can be used to identify the source location(s). Knowledge of the source location and characteristics are important to inform contamination control and cleanup operations. Various SI strategies that have been developed by researchers differ in their underlying assumptions and solution techniques. The following manuscript presents a systematic procedure for testing and evaluating SI methods. The performance of these SI methods is affected by various factors including the size of WDN model, measurement error, modeling error, time and number of contaminant injections, and time and number of measurements. This paper includes test cases that vary these factors and evaluates three SI methods on the basis of accuracy and specificity. The tests are used to review and compare these different SI methods, highlighting their strengths in handling various identification scenarios. These SI methods and a testing framework that includes the test cases and analysis tools presented in this paper have been integrated into EPA's Water Security Toolkit (WST), a suite of software tools to help researchers and others in the water industry evaluate and plan various response strategies in case of a contamination incident. Finally, a set of recommendations are made for users to consider when working with different categories of SI methods.
This Sandia National Laboratories, New Mexico Environmental Restoration Operations (ER) Consolidated Quarterly Report (ER Quarterly Report) fulfills all quarterly reporting requirements set forth in the Resource Conservation and Recovery Act Facility Operating Permit, and the Compliance Order on Consent.
Here, grain boundary-engineered materials are of immense interest for their corrosion resistance, fracture resistance and microstructural stability. This work contributes to a larger goal of understanding both the structure and thermodynamic properties of grain boundaries vicinal (within ±30°) to the Σ3(1 1 1) <1 1¯0> (coherent twin) boundary which is found in grain boundary-engineered materials. The misoriented boundaries vicinal to the twin show structural changes at elevated temperatures. In the case of nickel, this transition temperature is substantially below the melting point and at temperatures commonly reached during processing, making the existence of such boundaries very likely in applications. Thus, the thermodynamic stability of such features is thoroughly investigated in order to predict and fully understand the structure of boundaries vicinal to twins. Low misorientation angle grain boundaries (|θ| ≲ 16°) show distinct ±1/3(1 1 1) disconnections which accommodate misorientation in opposite senses. The two types of disconnection have differing low-temperature structures which show different temperature-dependent behaviours with one type undergoing a structural transition at approximately 600 K. At misorientation angles greater than approximately ±16°, the discrete disconnection nature is lost as the disconnections merge into one another. Free energy calculations demonstrate that these high-angle boundaries, which exhibit a transition from a planar to a faceted structure, are thermodynamically more stable in the faceted configuration.
Formulas for incremental or parallel computation of second order central moments have long been known, and recent extensions of these formulas to univariate and multivariate moments of arbitrary order have been developed. Formulas such as these, are of key importance in scenarios where incremental results are required and in parallel and distributed systems where communication costs are high. We survey these recent results, and improve them with arbitrary-order, numerically stable one-pass formulas which we further extend with weighted and compound variants. We also develop a generalized correction factor for standard two-pass algorithms that enables the maintenance of accuracy over nearly the full representable range of the input, avoiding the need for extended-precision arithmetic. We then empirically examine algorithm correctness for pairwise update formulas up to order four as well as condition number and relative error bounds for eight different central moment formulas, each up to degree six, to address the trade-offs between numerical accuracy and speed of the various algorithms. Finally, we demonstrate the use of the most elaborate among the above mentioned formulas, with the utilization of the compound moments for a practical large-scale scientific application.
First-principles calculations of electrical conductivity (σo) are revisited to determine the atomistic origin of its stochasticity in a distribution generated from sampling 14 ab-initio molecular dynamics configurations from 10 independently quenched models (n = 140) of substoichiometric amorphous Ta2O5, where each structure contains a neutral O monovacancy (VO0). Structural analysis revealed a distinct minimum Ta-Ta separation (dimer/trimer) corresponding to each VO0 location. Bader charge decomposition using a commonality analysis approach based on the σo distribution extremes revealed nanostructural signatures indicating that both the magnitude and distribution of cationic charge on the Ta subnetwork have a profound influence on σo. Furthermore, visualization of local defect structures and their electron densities reinforces these conclusions and suggests σo in the amorphous oxide is best suppressed by a highly charged, compact Ta cation shell that effectively screens and minimizes localized VO0 interaction with the a-Ta2O5 network; conversely, delocalization of VO0 corresponds to metallic character and high σo. The random network of a-Ta2O5 provides countless variations of an ionic configuration scaffold in which small perturbations affect the electronic charge distribution and result in a fixed-stoichiometry distribution of σo; consequently, precisely controlled and highly repeatable oxide fabrication processes are likely paramount for advancement of resistive memory technologies.
There is a rising interest in developing functional electronics using additively manufactured components. Considerations in materials selection and pathways to forming hybrid circuits and devices must demonstrate useful electronic function; must enable integration; and must complement the complex shape, low cost, high volume, and high functionality of structural but generally electronically passive additively manufactured components. This article reviews several emerging technologies being used in industry and research/development to provide integration advantages of fabricating multilayer hybrid circuits or devices. First, we review a maskless, noncontact, direct write (DW) technology that excels in the deposition of metallic colloid inks for electrical interconnects. Second, we review a complementary technology, aerosol deposition (AD), which excels in the deposition of metallic and ceramic powder as consolidated, thick conformal coatings and is additionally patternable through masking. As a result, we show examples of hybrid circuits/devices integrated beyond 2-D planes, using combinations of DW or AD processes and conventional, established processes.
Here, the contribution from loss of Li+ inventory to capacity fade is described for slow rates (C/10) and long-term cycling (up to 80 cycles). It was found through electrochemical testing and ex-situ Raman analysis that at these slow rates, the entirety of capacity loss up to 80 cycles can be explained by loss of Li+ inventory in the cell. The Raman spectrum of LiCoO2 is sensitive to the state of lithiation and can therefore be leveraged to quantify the state of lithiation for individual particles. With these Raman derived estimates, the lithiation state of the cathode in the discharged state is compared to electrochemical data as a function of cycle number. High correlation is found between Raman quantifications of cycleable lithium and the capacity fade. Additionally, the linear relationship between discharge capacity and cell overpotential suggests that the loss of capacity stems from an impedance rise of the electrodes, which based on Li inventory losses, is caused by SEI formation and repair.
Here, the widely recognized property improvements of nanocrystalline (NC) materials have generated significant interest, yet have been difficult to realize in engineering applications due to the propensity for grain growth in these interface-dense systems. While traditional pathways to thermal stabilization can slow the mobility of grain boundaries, recent theories suggest that solute segregation in NC alloy can reduce the grain boundary energy such that thermodynamic stabilization is achieved. Following the predictions of Murdock et al., here we compare for the first time the thermal stability of a predicted NC stable alloy (Fe-10at.% Mg) with a predicted non-NC stable alloy (Fe-10at.% Cu) using the same processing and characterization methodologies. Results indicate improved thermal stability of the Fe-Mg alloy in comparison to the Fe-Cu, and observed microstructures are consistent with those predicted by Monte Carlo simulations.
Fowler, Michael J.; Howard, Marylesa; Luttman, Aaron; Mitchell, Stephen E.; Webb, Timothy J.
One of the primary causes of blur in a high-energy X-ray imaging system is the shape and extent of the radiation source, or ‘spot’. It is important to be able to quantify the size of the spot as it provides a lower bound on the recoverable resolution for a radiograph, and penumbral imaging methods – which involve the analysis of blur caused by a structured aperture – can be used to obtain the spot’s spatial profile. We present a Bayesian approach for estimating the spot shape that, unlike variational methods, is robust to the initial choice of parameters. The posterior is obtained from a normal likelihood, which was constructed from a weighted least squares approximation to a Poisson noise model, and prior assumptions that enforce both smoothness and non-negativity constraints. A Markov chain Monte Carlo algorithm is used to obtain samples from the target posterior, and the reconstruction and uncertainty estimates are the computed mean and variance of the samples, respectively. Synthetic data-sets are used to demonstrate accurate reconstruction, while real data taken with high-energy X-ray imaging systems are used to demonstrate applicability and feasibility.
Maskey, Sabina; Osti, Naresh C.; Grest, Gary S.; Perahia, Dvora
The conformation and dynamics of luminescent polymers collapsed into nanoparticles or polydots were studied using fully atomistic molecular dynamics (MD) simulations, providing a first insight into their internal dynamics. Controlling the conformation and dynamics of confined polymers is essential for realization of the full potential of polydots in nanomedicine and biotechnology. Specifically, the shape and internal dynamics of polydots that consist of highly rigid dialkyl p-phenylene ethynylene (PPE) are probed as a function of temperature. At room temperature, the polydots are spherical without any correlations between the aromatic rings on the PPE backbone. With increasing temperature, they expand and become slightly aspherical; however, the polymers remain confined. The coherent dynamic structure factor reveals that the internal motion of the polymer backbone is arrested, and the side chains dominate the internal dynamics of the polydots. These new soft nanoparticles retain their overall shape and dynamics over an extended temperature range, and their conformation is tunable via their degree of expansion.
Mast, Daniel S.; Lavina, Barbara; Kim, Eunja; Siska, Emily M.; Weck, Philippe F.; Poineau, Frederic; Czerwinski, Kenneth R.; Forster, Paul M.
Here, the ambient temperature equation of state (EoS) of technetium metal has been measured by X-ray diffraction. The metal was compressed using a diamond anvil cell and using a 4:1 methanol-ethanol pressure transmitting medium. The maximum pressure achieved, as determined from the gold pressureEquation of state for technetium from X-ray diffraction and first-principle calculations scale, was 67 GPa. The compression data shows that the HCP phase of technetium is stable up to 67 GPa. The compression curve of technetium was also calculated using first-principles total-energy calculations. Utilizing a number of fitting strategies to compare the experimental and theoretical data it is determined that the Vinet equation of state with an ambient isothermal bulk modulus of B0T=288 GPa and a first pressure derivative of B'=5.9(2) best represent the compression behavior of technetium metal.
While structural damage detection based on flexural vibration shapes, such as mode shapes and steady-state response shapes under harmonic excitation, has been well developed, little attention is paid to that based on longitudinal vibration shapes that also contain damage information. This study originally formulates a slope vibration shape for damage detection in bars using longitudinal vibration shapes. To enhance noise robustness of the method, a slope vibration shape is transformed to a multiscale slope vibration shape in a multiscale domain using wavelet transform, which has explicit physical implication, high damage sensitivity, and noise robustness. These advantages are demonstrated in numerical cases of damaged bars, and results show that multiscale slope vibration shapes can be used for identifying and locating damage in a noisy environment. A three-dimensional (3D) scanning laser vibrometer is used to measure the longitudinal steady-state response shape of an aluminum bar with damage due to reduced cross-sectional dimensions under harmonic excitation, and results show that the method can successfully identify and locate the damage. Slopes of longitudinal vibration shapes are shown to be suitable for damage detection in bars and have potential for applications in noisy environments.
Pre-flown Vestas V-27 wind turbine blades will be used at the new Sandia SWiFT test facility near Lubbock,Texas. Proof load testing the blades in the flapwise direction was needed to verify the blades are still fit for service after 20 years offield operation. Testing was conducted in Building 252(A60)at the NWTC near Boulder,Colorado, between Oct 28(arrival of blades)to Dec 4,2015(completion of testing), under NREL subtasks WE141A50 and WE161A50. Single-point, quasi-static loads were applied to the blade at three independent load stations of 5.92 m, 8.69 m, and 11.62 m, using the facility gantry crane to apply the test load in the positive flapwise direction (low pressure surface in compression). Following the flapwise load cases, a single point load was applied at the 8.69 m station in the negative edgewise direction (leading edge in compression). No damage was observed due to the introduction of the proof loads. Several small pops were heard during testing on each blade, primarily during loading at the 8.69 m station. Compliance (deflection per load) varied less than 7% between all blades for all load cases. Table 3-1 summarizes the applied loads, deflections, and compliance for each blade for the flapwise test cases. Applied test loads are the raw load cell reading minus the total tare load (combined tare load of saddle and tare load for blade self weight). Deflection here is average deflection from both LE and TE string pot data.
Cl-initiated oxidation reactions of three small-chain methyl esters, methyl propanoate (CH3CH2COOCH3; MP), methyl butanoate (CH3CH2CH2COOCH3; MB), and methyl valerate (CH3CH2CH2CH2COOCH3; MV), are studied at 1 or 8 Torr and 550 and 650 K. Products are monitored as a function of mass, time, and photoionization energy using multiplexed photoionization mass spectrometry coupled to tunable synchrotron photoionization radiation. Pulsed photolysis of molecular chlorine is the source of Cl radicals, which remove an H atom from the ester, forming a free radical. In each case, after addition of O2 to the initial radicals, chain-terminating HO2-elimination reactions are observed to be important. Branching ratios among competing HO2-elimination channels are determined via absolute photoionization spectra of the unsaturated methyl ester coproducts. At 550 K, HO2-elimination is observed to be selective, resulting in nearly exclusive production of the conjugated methyl ester coproducts, methyl propenoate, methyl-2-butenoate, and methyl-2-pentenoate, respectively. However, in MV, upon raising the temperature to 650 K, other HO2-elimination pathways are observed that yield methyl-3-pentenoate and methyl-4-pentenoate. In each methyl ester oxidation reaction, a peak is observed at a mass consistent with cyclic ether formation, indicating chain-propagating OH loss/ring formation pathways via QOOH intermediates. Evidence is observed for the participation of resonance-stabilized QOOH in the most prominent cyclic ether pathways. Stationary point energies for HO2-elimination pathways and select cyclic ether formation channels are calculated at the CBS-QB3 level of theory and assist in the assignment of reaction pathways and final products.
Improved models of energy storage systems are needed to enable the electric grid’s adaptation to increasing penetration of renewables. This paper develops a generic empirical model of energy storage system performance agnostic of type, chemistry, design or scale. Parameters for this model are calculated using test procedures adapted from the US DOE Protocol for Uniformly Measuring and Expressing the Performance of Energy Storage. We then assess the accuracy of this model for predicting the performance of the TransPower GridSaver – a 1 MW rated lithium-ion battery system that underwent laboratory experimentation and analysis. The developed model predicts a range of energy storage system performance based on the uncertainty of estimated model parameters. Finally, this model can be used to better understand the integration and coordination of energy storage on the electric grid.
Using the kinetic Monte Carlo simulator, Stochastic Parallel PARticle Kinetic Simulator, from Sandia National Laboratories, a user routine has been developed to simulate mesoscale predictions of a grain structure near a moving heat source. Here, we demonstrate the use of this user routine to produce voxelized, synthetic, three-dimensional microstructures for electron-beam welding by comparing them with experimentally produced microstructures. When simulation input parameters are matched to experimental process parameters, qualitative and quantitative agreement for both grain size and grain morphology are achieved. The method is capable of simulating both single- and multipass welds. As a result, the simulations provide an opportunity for not only accelerated design but also the integration of simulation and experiments in design such that simulations can receive parameter bounds from experiments and, in turn, provide predictions of a resultant microstructure.
We study functional hybrid metasurfaces consisting of metal-dielectric nanoantennas that direct light from an incident plane wave or from localized light sources into a preferential direction. The directionality is obtained by carefully balancing the multipolar contributions to the scattering response from the constituents of the metasurface. The hybrid nanoantennas are composed of a plasmonic gold nanorod acting as a feed element and a silicon nanodisk acting as a director element. In order to experimentally realize this design, we have developed a two-step electron-beam lithography process in combination with a precision alignment step. The optical response of the fabricated sample is measured and reveals distinct signatures of coupling between the plasmonic and the dielectric nanoantenna elements that ultimately leads to unidirectional radiation of light.
Displacement damage caused by ions or neutrons in microelectronic devices can have significant effect on the performance of these devices. Therefore, it is important to predict not only the displacement damage profile, but also its magnitude precisely. Analytical methods and binary collision approximation codes working with amorphous targets use the concept of displacement energy, the energy that a lattice atom has to receive to create a permanent replacement. It was found that this "displacement energy" is direction dependent; it can range from 12 to 32 eV in silicon. Obviously, this model fails in BCA codes that work with crystalline targets, such as Marlowe. Marlowe does not use displacement energy; instead, it uses lattice binding energy only and then pairs the interstitial atoms with vacancies. Then based on the configuration of the Frenkel pairs it classifies them as close, near, or distant pairs, and considers the distant pairs the permanent replacements. Unfortunately, this separation is an ad hoc assumption, and the results do not agree with molecular dynamics calculations. After irradiation, there is a prompt recombination of interstitials and vacancies if they are nearby, within a recombination radius. In order to implement this recombination radius in Marlowe, we used the comparison of MD and Marlowe calculation in a range of ion energies in single crystal silicon target. The calculations showed that a single recombination radius of ∼7.4 Å in Marlowe for a range of ion energies gives an excellent agreement with MD.
A MS Excel program has been written that calculates ion channeling half-angles and minimum yields in cubic bcc, fcc and diamond lattice crystals. All of the tables and graphs in the three Ion Beam Analysis Handbooks that previously had to be manually looked up and read from were programed into Excel in handy lookup tables, or parameterized, for the case of the graphs, using rather simple exponential functions with different power functions of the arguments. The program then offers an extremely convenient way to calculate axial and planar half-angles, minimum yields, effects on half-angles and minimum yields of amorphous overlayers. The program can calculate these half-angles and minimum yields for 〈uvw〉 axes and [hkl] planes up to (5 5 5). The program is open source and available at http://www.sandia.gov/pcnsc/departments/iba/ibatable.html.
We demonstrate intrinsic, linearly polarized lasing from single GaN nanowires using cross-sectional shape control. A two-step top-down fabrication approach was employed to create straight nanowires with controllable rectangular cross-sections. A clear lasing threshold of 444 kW cm-2 and a narrow spectral line width of 0.16 nm were observed under optical pumping at room temperature, indicating the onset of lasing. The polarization was along the short dimension (y-direction) of the nanowire due to the higher transverse confinement factors for y-polarized transverse modes resulting from the rectangular nanowire cross-section. The results show that cross-sectioned shape control can enable inherent control over the polarization of nanowire lasers without additional environment requirements, such as placement onto lossy substrates.
Solar-blind photodetection and photoconductive gain >50 corresponding to a responsivity >8 A/W were observed for β-Ga2O3 Schottky photodiodes. The origin of photoconductive gain was investigated. Current-voltage characteristics of the diodes did not indicate avalanche breakdown, which excludes carrier multiplication by impact ionization as the source for gain. However, photocapacitance measurements indicated a mechanism for hole localization for above-band gap illumination, suggesting self-trapped hole formation. Comparison of photoconductivity and photocapacitance spectra indicated that self-trapped hole formation coincides with the strong photoconductive gain. It is concluded that self-trapped hole formation near the Schottky diode lowers the effective Schottky barrier in reverse bias, producing photoconductive gain. Ascribing photoconductive gain to an inherent property like self-trapping of holes can explain the operation of a variety of β-Ga2O3 photodetectors.
Active self-assembly offers a powerful route for the creation of dynamic multiscale structures that are presently inaccessible with standard microfabrication techniques. One such system uses the translation of microtubule filaments by surface-tethered kinesin to actively assemble nanocomposites with bundle, ring, and spool morphologies. Attempts to observe mechanisms involved in this active assembly system have been hampered by experimental difficulties with performing observation during buffer exchange and photodamage from fluorescent excitation. In the present work, we used a custom microfluidic device to remove these limitations and directly study ring/spool formation, including the earliest events (nucleation) that drive subsequent nanocomposite assembly. Three distinct formation events were observed: pinning, collisions, and induced curvature. Of these three, collisions accounted for the majority of event leading to ring/spool formation, while the rate of pinning was shown to be dependent on the amount of photodamage in the system. We further showed that formation mechanism directly affects the diameter and rotation direction of the resultant rings and spools. Overall, the fundamental understanding described in this work provides a foundation by which the properties of motor-driven, actively assembled nanocomposites may be tailored toward specific applications.
The electrical conductivity of materials under extremes of temperature and pressure is of crucial importance for a wide variety of phenomena, including planetary modeling, inertial confinement fusion, and pulsed power based dynamic materials experiments. There is a dearth of experimental techniques and data for highly compressed materials, even at known states such as along the principal isentrope and Hugoniot, where many pulsed power experiments occur. We present a method for developing, calibrating, and validating material conductivity models as used in magnetohydrodynamic (MHD) simulations. The difficulty in calibrating a conductivity model is in knowing where the model should be modified. Our method isolates those regions that will have an impact. It also quantitatively prioritizes which regions will have the most beneficial impact. Finally, it tracks the quantitative improvements to the conductivity model during each incremental adjustment. In this paper, we use an experiment on Sandia National Laboratories Z-machine to isentropically launch multiple flyer plates and, with the MHD code ALEGRA and the optimization code DAKOTA, calibrated the conductivity such that we matched an experimental figure of merit to +/-1%.
Magnetically applied pressure-shear (MAPS) is a new experimental technique that provides a platform for direct measurement of material strength at extreme pressures. The technique employs an imposed quasi-static magnetic field and a pulsed power generator that produces an intense current on a planar driver panel, which in turn generates high amplitude magnetically induced longitudinal compression and transverse shear waves into a planar sample mounted on the drive panel. In order to apply sufficiently high shear traction to the test sample, a high strength material must be used for the drive panel. Molybdenum is a potential driver material for the MAPS experiment because of its high yield strength and sufficient electrical conductivity. To properly interpret the results and gain useful information from the experiments, it is critical to have a good understanding and a predictive capability of the mechanical response of the driver. In this work, the inelastic behavior of molybdenum under uniaxial compression and biaxial compression-shear ramp loading conditions is experimentally characterized. It is observed that an imposed uniaxial magnetic field ramped to approximately 10 T through a period of approximately 2500 μs and held near the peak for about 250 μs before being tested appears to anneal the molybdenum panel. In order to provide a physical basis for model development, a general theoretical framework that incorporates electromagnetic loading and the coupling between the imposed field and the inelasticity of molybdenum was developed. Based on this framework, a multi-axial continuum model for molybdenum under electromagnetic loading is presented. The model reasonably captures all of the material characteristics displayed by the experimental data obtained from various experimental configurations. In addition, data generated from shear loading provide invaluable information not only for validating but also for guiding the development of the material model for multiaxial loadings.
Recent studies have shown the potential for nanocrystalline metals to possess excellent fatigue resistance compared to their coarse-grained counterparts. Although the mechanical properties of nanocrystalline metals are believed to be particularly susceptible to material defects, a systematic study of the effects of geometric discontinuities on their fatigue performance has not yet been performed. In the present work, nanocrystalline Ni–40 wt%Fe containing both intrinsic and extrinsic defects were tested in tension–tension fatigue. The defects were found to dramatically reduce the fatigue resistance, which was attributed to the relatively high notch sensitivity in the nanocrystalline material. Microstructural analysis within the crack-initiation zones underneath the defects revealed cyclically-induced abnormal grain growth (AGG) as a predominant deformation and crack initiation mechanism during high-cycle fatigue. Furthermore, the onset of AGG and the ensuing fracture is likely accelerated by the stress concentrations, resulting in the reduced fatigue resistance compared to the relatively defect-free counterparts.
Our research group has been developing a fast neutron imaging platform to enhance the capabilities of emergency responders in the localization and characterization of special nuclear material. This mobile imager of neutrons for emergency responders (MINER) is a compact neutron scatter camera optimized to provide omni-directional (4-Pi) imaging with only a ~twofold decrease in sensitivity compared to our much larger neutron scatter cameras. The system performance is tuned for fission energy neutron imaging and spectroscopy, and it also can function as a Compton camera for gamma imaging. Results will be presented relating to detector response as well as several measurement campaigns at external facilities.
Our research group has been developing a fast neutron imaging platform to enhance the capabilities of emergency responders in the localization and characterization of special nuclear material. This mobile imager of neutrons for emergency responders (MINER) is a compact neutron scatter camera optimized to provide omni-directional (4-Pi) imaging with only a ~twofold decrease in sensitivity compared to our much larger neutron scatter cameras. The system performance is tuned for fission energy neutron imaging and spectroscopy, and it also can function as a Compton camera for gamma imaging. Results will be presented relating to detector response as well as several measurement campaigns at external facilities.
Leuty, Gary M.; Tsige, Mesfin; Grest, Gary S.; Rubinstein, Michael
Molecular dynamics simulations of a coarse-grained bead-spring model have been used to study the effects of molecular crowding on the accumulation of tension in the backbone of bottle-brush polymers tethered to a flat substrate. The number of bottle-brushes per unit surface area, Σ, as well as the lengths of the bottle-brush backbones Nbb (50 ≤ Nbb ≤ 200) and side chains Nsc (50 ≤ Nsc ≤ 200) were varied to determine how the dimensions and degree of crowding of bottle-brushes give rise to bond tension amplification along the backbone, especially near the substrate. From these simulations, we have identified three separate regimes of tension. For low Σ, the tension is due solely to intramolecular interactions and is dominated by the side chain repulsion that governs the lateral brush dimensions. With increasing Σ, the interactions between bottle-brush polymers induce compression of the side chains, transmitting increasing tension to the backbone. For large Σ, intermolecular side chain repulsion increases, forcing side chain extension and reorientation in the direction normal to the surface and transmitting considerable tension to the backbone.
Researchers have rigorously verified collision and energy partitioning algorithms to give our users confidence that Aleph can correctly predict these atomic processes for their applications.
This code implements a variety of "deep learning" algorithms from openly published academic journals. These deep learning algorithms allow the user to train neural networks to do various tasks. Such tasks include predicting future values in a time sequence, categorizing images or compressing information. The code is mostly written in the easy to understand Matlab / GNU Octave language, which enables rapid prototyping and understanding for research and educational purposes.
Oxygen reduction at the cathode of fuel cells typically requires a platinum-based material to catalyse the reaction, but lower-cost, more stable catalysts are sought. Here, an intrinsically conductive metal–organic framework based on cheaper elements is shown to be a durable, structurally well-defined catalyst for this reaction.
A minimax estimator has the minimum possible error ("risk") in the worst case. We construct the first minimax estimators for quantum state tomography with relative entropy risk. The minimax risk of nonadaptive tomography scales as O(1/N) - in contrast to that of classical probability estimation, which is O(1/N) - where N is the number of copies of the quantum state used. We trace this deficiency to sampling mismatch: future observations that determine risk may come from a different sample space than the past data that determine the estimate. This makes minimax estimators very biased, and we propose a computationally tractable alternative with similar behavior in the worst case, but superior accuracy on most states.
The role of solute attractive forces on hydrophobic interactions is studied by coordinated development of theory and simulation results for Ar atoms in water. We present a concise derivation of the local molecular field (LMF) theory for the effects of solute attractive forces on hydrophobic interactions, a derivation that clarifies the close relation of LMF theory to the EXP approximation applied to this problem long ago. The simulation results show that change from purely repulsive atomic solute interactions to include realistic attractive interactions diminishes the strength of hydrophobic bonds. For the Ar-Ar rdfs considered pointwise, the numerical results for the effects of solute attractive forces on hydrophobic interactions are opposite in sign and larger in magnitude than predicted by LMF theory. That comparison is discussed from the point of view of quasichemical theory, and it is suggested that the first reason for this difference is the incomplete evaluation within LMF theory of the hydration energy of the Ar pair. With a recent suggestion for the system-size extrapolation of the required correlation function integrals, the Ar-Ar rdfs permit evaluation of osmotic second virial coefficients B2. Those B2's also show that incorporation of attractive interactions leads to more positive (repulsive) values. With attractive interactions in play, B2 can change from positive to negative values with increasing temperatures. This is consistent with the puzzling suggestions of decades ago that B2 ≈ 0 for intermediate cases of temperature or solute size. In all cases here, B2 becomes more attractive with increasing temperature.
Sandia National Laboratories is pursuing a variation of Magneto-Inertial Fusion called Magnetized Liner Inertial Fusion, or MagLIF. The MagLIF approach requires magnetization of the deuterium fuel, which is accomplished by an initial external B-Field and laser-driven pre-heat. Although magnetization is crucial to the concept, it is challenging to couple sufficient energy to the fuel, since laser-plasma instabilities exist, and a compromise between laser spot size, laser entrance window thickness, and fuel density must be found. Ultimately, nonlinear processes in laser plasma interaction, or laser-plasma instabilities (LPI), complicate the deposition of laser energy by enhanced absorption, backscatter, filamentation and beam-spray. We determine and discuss key LPI processes and mitigation methods. Results with and without improvement measures are presented.
The question of how to accurately model turbulent flows is one of the most long-standing open problems in physics. Advances in high performance computing have enabled direct numerical simulations of increasingly complex flows. Nevertheless, for most flows of engineering relevance, the computational cost of these direct simulations is prohibitive, necessitating empirical model closures for the turbulent transport. These empirical models are prone to "model form uncertainty" when their underlying assumptions are violated. Understanding, quantifying, and mitigating this model form uncertainty has become a critical challenge in the turbulence modeling community. This paper will discuss strategies for using machine learning to understand the root causes of the model form error and to develop model corrections to mitigate this error. Rule extraction techniques are used to derive simple rules for when a critical model assumption is violated. The physical intuition gained from these simple rules is then used to construct a linear correction term for the turbulence model which shows improvement over naive linear fits.
Journal of Polymer Science, Part B: Polymer Physics
Wijesinghe, Sidath; Maskey, Sabina; Perahia, Dvora; Grest, Gary S.
The conformation of dinonyl poly para phenylene ethynylenes (PPEs) with carboxylate side chains, equilibrated in solvents of different quality have been studied using molecular dynamics simulations. PPEs are of interest because of their tunable electro-optical properties, chemical diversity, and functionality which are essential in wide range of applications. The polymer conformation determines the conjugation length and their assembly mode and affects electro-optical properties which are critical in current and potential uses. This study investigates the effect of carboxylate fraction on PPEs side chains on the conformation of chains in the dilute limit, in solvents of different quality. The dinonyl PPE chains are modeled atomistically, where the solvents are modeled both implicitly and explicitly. Dinonyl PPEs maintained a stretched out conformation up to a carboxylate fraction f of 0.7 in all solvents studied. The nonyl side chains are extended and oriented away from the PPE backbone in toluene and in implicit good solvent, whereas in water and implicit poor solvent, the nonyl side chains are collapsed toward the PPE backbone. Rotation around the aromatic ring is fast and no long range correlations are seen within the backbone.
For the purposes of making reliable first-principles predictions of defect energies in semiconductors, it is crucial to distinguish between effective-mass-like defects, which cannot be treated accurately with existing supercell methods, and deep defects, for which density functional theory calculations can yield reliable predictions of defect energy levels. The gallium antisite defect GaAs is often associated with the 78/203 meV shallow double acceptor in Ga-rich gallium arsenide. Within a conceptual framework of level patterns, analyses of structure and spin stabilization can be used within a supercell approach to distinguish localized deep defect states from shallow acceptors such as BAs. This systematic approach determines that the gallium antisite supercell results has signatures inconsistent with an effective mass state and cannot be the 78/203 shallow double acceptor. The properties of the Ga antisite in GaAs are described, total energy calculations that explicitly map onto asymptotic discrete localized bulk states predict that the Ga antisite is a deep double acceptor and has at least one deep donor state.
Experiments on the Sandia Z pulsed-power accelerator have demonstrated the ability to produce warm dense matter (WDM) states with unprecedented uniformity, duration, and size, which are ideal for investigations of fundamental WDM properties. For the first time, space-resolved x-ray Thomson scattering (XRTS) spectra from shocked carbon foams were recorded on Z. The large (>20 MA) electrical current produced by Z was used to launch Al flyer plates up to 25 km/s. The impact of the flyer plate on a CH2 foam target produced a shocked state with an estimated pressure of 0.75 Mbar, density of 0.52 g/cm3, and temperature of 4.3 eV. Both unshocked and shocked portions of the foam target were probed with 6.2 keV x-rays produced by focusing the Z-Beamlet laser onto a nearby Mn foil. The data are composed of three spatially distinct spectra that were simultaneously captured with a single spectrometer with high spectral (4.8 eV) and spatial (190 μm) resolutions. Detailed spectral information from three target locations is provided simultaneously: the incident x-ray source, the scattered signal from unshocked foam, and the scattered signal from shocked foam.
The second Sandia Fracture Challenge illustrates that predicting the ductile fracture of Ti-6Al-4V subjected to moderate and elevated rates of loading requires thermomechanical coupling, elasto-thermo-poro-viscoplastic constitutive models with the physics of anisotropy and regularized numerical methods for crack initiation and propagation. We detail our initial approach with an emphasis on iterative calibration and systematically increasing complexity to accommodate anisotropy in the context of an isotropic material model. Blind predictions illustrate strengths and weaknesses of our initial approach. We then revisit our findings to illustrate the importance of including anisotropy in the failure process. Mesh-independent solutions of continuum damage models having both isotropic and anisotropic yields surfaces are obtained through nonlocality and localization elements.
Silicone foams have been used in a variety of applications from gaskets to cushioning pads over a wide range of environments. Particularly, silicone foams are used as a shock mitigation material for shock and vibration applications. Understanding the shock mitigation response, particularly in the frequency domain, is critical for optimal designs to protect internal devices and components more effectively and efficiently. The silicone foams may be subjected to pre-strains during the assembly process which may consequently influence the frequency response with respect to shock mitigation performance. A Kolsky compression bar was modified with pre-compression capabilities to characterize the shock mitigation response of silicone foam in the frequency domain to determine the effect of pre-strain. A silicone sample was also intentionally subjected to repeated pre-strain and dynamic loadings to explore the effect of repeated loading on the frequency response of shock mitigation.
A series of alkali metal yttrium neo-pentoxide ([AY(ONep)4]) compounds were developed as precursors to alkali yttrium oxide (AYO2) nanomaterials. The reaction of yttrium amide ([Y(NR2)3] where R=Si(CH3)3) with four equivalents of H-ONep followed by addition of [A(NR2)] (A=Li, Na, K) or Ao (Ao=Rb, Cs) led to the formation of a complex series of AnY(ONep)3+n species, crystallographically identified as [Y2Li3(μ3-ONep)(μ3-HONep)(μ-ONep)5(ONep)3(HONep)2] (1), [YNa2(μ3-ONep)4(ONep)]2 (2), {[Y2K3(μ3-ONep)3(μ-ONep)4(ONep)2(ηξ-tol)2][Y4K2(μ4-O)(μ3-ONep)8(ONep)4]•ηx-tol]} (3), [Y4K2(μ4-O)(μ3-ONep)8(ONep)4] (3 a), [Y2Rb3(μ4-ONep)3(μ-ONep)6] (4), and [Y2Cs4(μ6-O)(μ3-ONep)6(μ3-HONep)2(ONep)2(ηx-tol)4]•tol (5). Compounds 1–5 were investigated as single source precursors to AYOx nanomaterials following solvothermal routes (pyridine, 185 oC for 24 h). The final products after thermal processing were found by powder X-ray diffraction experiments to be Y2O3 with variable sized particles based on transmission electron diffraction. Energy dispersive X-ray spectroscopy studies indicated that the heavier alkali metal species were present in the isolated nanomaterials.
Nanoparticles (NPs) grafted with organic layers form hybrids able to retain their unique properties through integration into the mesoscopic scale. The organic layer structure and response often determine the functionality of the hybrids on the mesoscopic length scale. Using molecular dynamics (MD) simulations, we probe the conformation of luminescent rigid polymers, dialkyl poly(p-phenylene ethynylene)s (PPE), end-grafted onto a silica nanoparticle in different solvents as the molecular weights and polymer coverages are varied. We find that, in contrast to NP-grafted flexible polymers, the chains are fully extended independent of the solvent. In toluene and decane, which are good solvents, the grafted PPEs chains assume a similar conformation to that observed in dilute solutions. In water, which is a poor solvent for the PPEs, the polymer chains form one large cluster but remain extended. The radial distribution of the chains around the core of the nanoparticle is homogeneous in good solvents, whereas in poor solvents clusters are formed independent of molecular weights and coverages. The clustering is distinctively different from the response of grafted flexible and semiflexible polymers.
In the last few decades, various solutions have been proposed to increase the modulation bandwidth and, consequently, the transmission bit-rate of semiconductor lasers. In this manuscript, we discuss a design procedure for a recently proposed laser cavity realized with the monolithic integration of two distributed Bragg reflector (DBR) lasers allowing one to extend the modulation bandwidth. Such an extension is obtained introducing in the dynamic response a photon-photon resonance (PPR) at a frequency higher than the modulation bandwidth of the corresponding single-section laser. Design guidelines will be proposed, and dynamic small and large signal simulations results, calculated using a finite difference traveling wave (FDTW) numerical simulator, will be discussed to confirm the design results. The effectiveness of the design procedure is verified in a structure with PPR frequency at 35GHz allowing one to obtain an open eye diagram for a non-return-to-zero (NRZ) digital signal up to 80 GHz. Furthermore, the investigation of the rich dynamics of this structure shows that with proper bias conditions, it is possible to obtain also a tunable self-pulsating signal in a frequency range related to the PPR design.
We describe an optical coating design suitable for broad bandwidth high reflection (BBHR) at 45° angle of incidence (AOI), P polarization (Ppol) of femtosecond (fs) laser pulses whose wavelengths range from 800 to 1000 nm. Our design process is guided by quarter-wave HR coating properties. The design must afford low group delay dispersion (GDD) for reflected light over the broad, 200 nm bandwidth in order to minimize temporal broadening of the fs pulses due to dispersive alteration of relative phases between their frequency components. The design should also be favorable to high laser-induced damage threshold (LIDT). We base the coating on TiO 2 /SiO 2 layer pairs produced by means of e-beam evaporation with ion-assisted deposition, and use OptiLayer Thin Film Software to explore designs starting with TiO 2 /SiO 2 layers having thicknesses in a reverse chirped arrangement. This approach led to a design with R > 99% from 800 to 1000 nm and GDD < 20 fs 2 from 843 to 949 nm (45° AOI, Ppol). The design's GDD behaves in a smooth way, suitable for GDD compensation techniques, and its electric field intensities show promise for high LIDTs. Reflectivity and GDD measurements for the initial test coating indicate good performance of the BBHR design. Subsequent coating runs with improved process calibration produced two coatings whose HR bands satisfactorily meet the design goals. For the sake of completeness, we summarize our previously reported transmission spectra and LIDT test results with 800 ps, 8 ps and 675 fs pulses for these two coatings, and present a table of the LIDT results we have for all of our TiO 2 /SiO 2 BBHR coatings, showing the trends with test laser pulse duration from the ns to sub-ps regimes.
Ductile failure of structural metals is relevant to a wide range of engineering scenarios. Computational methods are employed to anticipate the critical conditions of failure, yet they sometimes provide inaccurate and misleading predictions. Challenge scenarios, such as the one presented in the current work, provide an opportunity to assess the blind, quantitative predictive ability of simulation methods against a previously unseen failure problem. Rather than evaluate the predictions of a single simulation approach, the Sandia Fracture Challenge relies on numerous volunteer teams with expertise in computational mechanics to apply a broad range of computational methods, numerical algorithms, and constitutive models to the challenge. This exercise is intended to evaluate the state of health of technologies available for failure prediction. In the first Sandia Fracture Challenge, a wide range of issues were raised in ductile failure modeling, including a lack of consistency in failure models, the importance of shear calibration data, and difficulties in quantifying the uncertainty of prediction [see Boyce et al. (Int J Fract 186:5–68, 2014) for details of these observations]. This second Sandia Fracture Challenge investigated the ductile rupture of a Ti–6Al–4V sheet under both quasi-static and modest-rate dynamic loading (failure in (Formula presented.) 0.1 s). Like the previous challenge, the sheet had an unusual arrangement of notches and holes that added geometric complexity and fostered a competition between tensile- and shear-dominated failure modes. The teams were asked to predict the fracture path and quantitative far-field failure metrics such as the peak force and displacement to cause crack initiation. Fourteen teams contributed blind predictions, and the experimental outcomes were quantified in three independent test labs. Additional shortcomings were revealed in this second challenge such as inconsistency in the application of appropriate boundary conditions, need for a thermomechanical treatment of the heat generation in the dynamic loading condition, and further difficulties in model calibration based on limited real-world engineering data. As with the prior challenge, this work not only documents the ‘state-of-the-art’ in computational failure prediction of ductile tearing scenarios, but also provides a detailed dataset for non-blind assessment of alternative methods.
A well-field within a uranium (U) plume in the groundwater-surface water transition zone was monitored for a 3 year period for water table elevation and dissolved solutes. The plume discharges to the Columbia River, which displays a dramatic spring stage surge resulting from snowmelt. Groundwater exhibits a low hydrologic gradient and chemical differences with river water. River water intrudes the site in spring. Specific aims were to assess the impacts of river intrusion on dissolved uranium (Uaq), specific conductance (SpC), and other solutes, and to discriminate between transport, geochemical, and source term heterogeneity effects. Time series trends for Uaq and SpC were complex and displayed large temporal and well-to-well variability as a result of water table elevation fluctuations, river water intrusion, and changes in groundwater flow directions. The wells were clustered into subsets exhibiting common behaviors resulting from the intrusion dynamics of river water and the location of source terms. Hot-spots in Uaq varied in location with increasing water table elevation through the combined effects of advection and source term location. Heuristic reactive transport modeling with PFLOTRAN demonstrated that mobilized Uaq was transported between wells and source terms in complex trajectories, and was diluted as river water entered and exited the groundwater system. While Uaq time-series concentration trends varied significantly from year-to-year as a result of climate-caused differences in the spring hydrograph, common and partly predictable response patterns were observed that were driven by water table elevation, and the extent and duration of river water intrusion.
We demonstrate a photonic waveguide technology based on a two-material core, in which light is controllably and repeatedly transferred back and forth between sub-micron thickness crystalline layers of Si and LN bonded to one another, where the former is patterned and the latter is not. In this way, the foundry-based wafer-scale fabrication technology for silicon photonics can be leveraged to form lithium-niobate based integrated optical devices. Using two different guided modes and an adiabatic mode transition between them, we demonstrate a set of building blocks such as waveguides, bends, and couplers which can be used to route light underneath an unpatterned slab of LN, as well as outside the LN-bonded region, thus enabling complex and compact lightwave circuits in LN alongside Si photonics with fabrication ease and low cost.
The fast pyrolysis behaviour of pretreated banagrass was examined at four temperatures (between 400 and 600 C) and four residence times (between ∼1.2 and 12 s). The pretreatment used water washing/leaching to reduce the inorganic content of the banagrass. Yields of bio-oil, permanent gases and char were determined at each reaction condition and compared to previously published results from untreated banagrass. Comparing the bio-oil yields from the untreated and pretreated banagrass shows that the yields were greater from the pretreated banagrass by 4 to 11 wt% (absolute) at all reaction conditions. The effect of pretreatment (i.e. reducing the amount of ash, and alkali and alkali earth metals) on pyrolysis products is: 1) to increase the dry bio-oil yield, 2) to decrease the amount of undetected material, 3) to produce a slight increase in CO yield or no change, 4) to slightly decrease CO2 yield or no change, and 5) to produce a more stable bio-oil (less aging). Char yield and total gas yield were unaffected by feedstock pretreatment. Four other tropical biomass species were also pyrolyzed under one condition (450°C and 1.4 s residence time) for comparison to the banagrass results. The samples include two hardwoods: leucaena and eucalyptus, and two grasses: sugarcane bagasse and energy-cane. A sample of pretreated energy-cane was also pyrolyzed. Of the materials tested, the best feedstocks for fast pyrolysis were sugarcane bagasse, pretreated energy cane and eucalyptus based on the yields of 'dry bio-oil', CO and CO2. On the same basis, the least productive feedstocks are untreated banagrass followed by pretreated banagrass and leucaena.
Patrick, Joshua D.; Harvill, Jane L.; Hansen, Clifford W.
We evaluate semiparametric spatio-temporal models for global horizontal irradiance at high spatial and temporal resolution. These models represent the spatial domain as a lattice and are capable of predicting irradiance at lattice points, given data measured at other lattice points. Using data from a 1.2 MW PV plant located in Lanai, Hawaii, we show that a semiparametric model can be more accurate than simple interpolation between sensor locations. We investigate spatio-temporal models with separable and nonseparable covariance structures and find no evidence to support assuming a separable covariance structure. Our results indicate a promising approach for modeling irradiance at high spatial resolution consistent with available ground-based measurements. Such modeling may find application in design, valuation, and operation of fleets of utility-scale photovoltaic power systems.
Geometric partitioning is fast and effective for load-balancing dynamic applications, particularly those requiring geometric locality of data (particle methods, crash simulations). We present, to our knowledge, the first parallel implementation of a multidimensional-jagged geometric partitioner. In contrast to the traditional recursive coordinate bisection algorithm (RCB), which recursively bisects subdomains perpendicular to their longest dimension until the desired number of parts is obtained, our algorithm does recursive multi-section with a given number of parts in each dimension. By computing multiple cut lines concurrently and intelligently deciding when to migrate data while computing the partition, we minimize data movement compared to efficient implementations of recursive bisection. We demonstrate the algorithm's scalability and quality relative to the RCB implementation in Zoltan on both real and synthetic datasets. Our experiments show that the proposed algorithm performs and scales better than RCB in terms of run-time without degrading the load balance. Our implementation partitions 24 billion points into 65,536 parts within a few seconds and exhibits near perfect weak scaling up to 6K cores.
This paper presents an extension of the all-quad meshing algorithm called LayTracks to generate high quality hex-dominant meshes of general solids. LayTracks3D uses the mapping between the Medial Axis (MA) and the boundary of the 3D domain to decompose complex 3D domains into simpler domains called Tracks. Tracks in 3D have no branches and are symmetric, non-intersecting, orthogonal to the boundary, and the shortest path from the MA to the boundary. These properties of tracks result in desired meshes with near cube shape elements at the boundary, structured mesh along the boundary normal with any irregular nodes restricted to the MA, and sharp boundary feature preservation. The algorithm has been tested on a few industrial CAD models and hex-dominant meshes are shown in the Results section. Work is underway to extend LayTracks3D to generate all-hex meshes.
Photo-conductive AFM spectroscopy ('pcAFMs') is proposed as a high-resolution approach for investigating nanostructured photovoltaics, uniquely providing nanoscale maps of photovoltaic (PV) performance parameters such as the short circuit current, open circuit voltage, maximum power, or fill factor. The method is demonstrated with a stack of 21 images acquired during in situ illumination of micropatterned polycrystalline CdTe/CdS, providing more than 42 000 I/V curves spatially separated by ~5 nm. For these CdTe/CdS microcells, the calculated photoconduction ranges from 0 to 700 picoSiemens (pS) upon illumination with ~1.6 suns, depending on location and biasing conditions. Mean short circuit currents of 2 pA, maximum powers of 0.5 pW, and fill factors of 30% are determined. The mean voltage at which the detected photocurrent is zero is determined to be 0.7 V. Significantly, enhancements and reductions in these more commonly macroscopic PV performance metrics are observed to correlate with certain grains and grain boundaries, and are confirmed to be independent of topography. These results demonstrate the benefits of nanoscale resolved PV functional measurements, reiterate the importance of microstructural control down to the nanoscale for 'PV devices, and provide a widely applicable new approach for directly investigating PV materials.
In this paper, we derive a strengthened MILP formulation for certain gas turbine unit commitment problems, in which the ramping rates are no smaller than the minimum generation amounts. This type of gas turbines can usually start-up faster and have a larger ramping rate, as compared to the traditional coal-fired power plants. Recently, the number of this type of gas turbines increases significantly due to affordable gas prices and their scheduling flexibilities to accommodate intermittent renewable energy generation. In this study, several new families of strong valid inequalities are developed to help reduce the computational time to solve these types of problems. Meanwhile, the validity and facet-defining proofs are provided for certain inequalities. Finally, numerical experiments on a modified IEEE 118-bus system and the power system data based on recent studies verify the effectiveness of applying our formulation to model and solve this type of gas turbine unit commitment problems, including reducing the computational time to obtain an optimal solution or obtaining a much smaller optimality gap, as compared to the default CPLEX, when the time limit is reached with no optimal solutions obtained.
The vector is a fundamental data structure, which provides constant-time access to a dynamically-resizable range of elements. Currently, there exist no wait-free vectors. The only non-blocking version supports only a subset of the sequential vector API and exhibits significant synchronization overhead caused by supporting opposing operations. Since many applications operate in phases of execution, wherein each phase only a subset of operations are used, this overhead is unnecessary for the majority of the application. To address the limitations of the non-blocking version, we present a new design that is wait-free, supports more of the operations provided by the sequential vector, and provides alternative implementations of key operations. These alternatives allow the developer to balance the performance and functionality of the vector as requirements change throughout execution. Compared to the known non-blocking version and the concurrent vector found in Intel’s TBB library, our design outperforms or provides comparable performance in the majority of tested scenarios. Over all tested scenarios, the presented design performs an average of 4.97 times more operations per second than the non-blocking vector and 1.54 more than the TBB vector. In a scenario designed to simulate the filling of a vector, performance improvement increases to 13.38 and 1.16 times. This work presents the first ABA-free non-blocking vector. Finally, unlike the other non-blocking approach, all operations are wait-free and bounds-checked and elements are stored contiguously in memory.
Predictions for ductile tearing of a geometrically complex Ti-6Al-4V plate were generated using a Unified Creep Plasticity Damage model in fully coupled thermal stress simulations. Uniaxial tension and butterfly shear tests performed at displacement rates of 0.0254 and 25.4 mm/s were also simulated. Results from these simulations revealed that the material temperature increase due to plastic work can have a dramatic effect on material ductility predictions in materials that exhibit little strain hardening. Furthermore, this occurs because the temperature increase causes the apparent hardening of the material to decrease which leads to the initiation of deformation localization and subsequent ductile tearing earlier in the loading process.
The purpose of this CWL Annual Post-Closure Care Report is to document monitoring, inspection, maintenance, and repair activities conducted during CY 2015 in accordance with Attachment 1 of the PCCP (NMED October 2009 and subsequent revisions). This annual report documents post-closure care activities conducted from January through December 2015 and fulfills the PCCP requirement for annual reporting to the New Mexico Environment Department (NMED).
The Corrective Action Management Unit (CAMU) at Sandia National Laboratories, New Mexico (SNL/NM) consisted of a containment cell, two treatment systems, and four associated waste staging and storage areas, and support areas; all were used for management of remediation wastes between 1997 and 2003. The U.S. Department of Energy and Sandia Corporation (Sandia) operated the CAMU in accordance with the requirements of: (1) The Class III Permit Modification Request for the CAMU (SNL/NM September 1997, as amended); and (2) Module IV of Permit NM5890110518, “Special Conditions Pursuant to the 1984 Hazardous and Solid Waste Amendments (HSWA) to Resource Conservation and Recovery Act (RCRA) for Sandia National Laboratory, U.S. Environmental Protection Agency (EPA) Identification Number NM5890110518,” August 1993, (EPA 1993, as amended) and subsequently administered by the New Mexico Environment Department (NMED). The CAMU underwent closure in 2003 in accordance with the Closure Plan in Appendix D of the Class III Permit Modification (SNL/NM September 1997). The containment cell was closed with wastes in place. Hazardous wastes were removed from all other CAMU systems and areas, and they underwent clean closure. The NMED approved completion of closure in May 2004 (Kieling May 2004). From May 2004 until February 2015, the containment cell underwent post-closure care (i.e., monitoring, inspections, maintenance, and repairs) in accordance with the Closure Plan. On December 19, 2014, NMED renewed the Hazardous Waste Facility Operating Permit (Permit) for Sandia National Laboratories. The renewed Permit was issued January 26, 2015 (NMED January 2015) with an effective date of February 26, 2015. The CAMU is undergoing post-closure care in accordance with the Permit.
Long-term controls were maintained at 21 Solid Waste Management Units (SWMUs) and Areas of Concern (AOCs) in accordance with the requirements of the Long-Term Monitoring and Maintenance Plan for SWMUs and AOCs Granted Corrective Action Complete with Controls in Attachment M of the Resource Conservation and Recovery Act Facility Operating Permit which took effect February 26, 2015. Maintenance and controls at these SWMUs and AOCs are described and documented in this report. The need for maintenance or repair activities were not identified at any of the inspected SWMUs or AOCs. Four repairs identified in 2014 were completed in 2015. Based upon the inspections performed and site conditions observed, the administrative and physical institutional controls in place at the SWMUs and AOCs are effectively providing continued protection of human health and the environment.