Laser propagation experiments using four beams of the National Ignition Facility to deliver up to 35 kJ of laser energy at 351 nm laser wavelength to heat magnetized liner inertial fusion-scale (1 cm-long), hydrocarbon-filled gas pipe targets to ∼keV electron temperatures have demonstrated energy coupling >20 kJ with essentially no backscatter in 15% critical electron density gas fills with 0-19 T applied axial magnetic fields. The energy coupling is also investigated for an electron density of 11.5% critical and for applied field strengths up to 24 T at both densities. This spans a range of Hall parameters 0 < ω c e τ e i ≲2, where a Hall parameter of 0.5 is expected to reduce electron thermal conduction across the field lines by a factor of 4-5 for the conditions of these experiments. At sufficiently high applied field strength (and therefore Hall parameter), the measured laser propagation speed through the targets increases in the measurements, consistent with reduced perpendicular electron thermal transport; this reduces the coupled energy to the target once the laser burns through the gas pipe. The results compare well with a 1D analytic propagation model for inverse Bremsstrahlung absorption.
The cost of photovoltaic (PV) modules has declined by 85% since 2010. To achieve this reduction, manufacturers altered module designs and bill of materials; changes that could affect module durability and reliability. To determine if these changes have affected module durability, we measured the performance degradation of 834 fielded PV modules representing 13 module types from 7 manufacturers in 3 climates over 5 years. Degradation rates (Rd) are highly nonlinear over time, and seasonal variations are present in some module types. Mean and median degradation rate values of −0.62%/year and −0.58%/year, respectively, are consistent with rates measured for older modules. Of the 23 systems studied, 6 have degradation rates that will exceed the warranty limits in the future, whereas 13 systems demonstrate the potential of achieving lifetimes beyond 30 years, assuming Rd trends have stabilized.
Nakagawa, Seiji; Kibikas, William M.; Chang, Chun; Kneafsey, Timothy; Dobson, Patrick; Samuel, Abraham; Bruce, Stephen; Kaargeson-Loe, Nils; Bauer, Stephen J.
This work summarizes the findings of a reduced order model (ROM) study performed using Sierra ROM module Pressio_Aria on Sandia National Laboratories' (SNL) Crash-Burn L2 milestone thermal model with pristine geometry. Comparisons are made to full order model (FOM) results for this same Crash-Burn model using Sierra multiphysics module Aria.
Marine energy generation technologies such as wave and tidal power have great potential in meeting the need for renewable energy in the years ahead. Yet, many challenges remain associated with marine-based systems because of the corrosive environment. Conventional materials like metals are subject to rapid corrosive breakdown, crippling the lifespan of structures in such environments. Fiber-reinforced polymer composites offer an appealing alternative in their strength and corrosion resistance, but can experience degradation of mechanical properties as a result of moisture absorption. An investigation is conducted to test the application of a technique for micromechanical analysis of composites, known as multicontinuum theory and demonstrated in past works, as a mechanism for predicting the effects of prolonged moisture absorption on the performance of fiber-reinforced composites. Experimental tensile tests are performed on composite coupons with and without prolonged exposure to a salt water solution to obtain stiffness and strength properties. Multicontinuum theory is applied in conjunction with micromechanical modeling to deduce the effects of moisture absorption on the behavior of constituent materials within the composites. The results are consistent with experimental observations when guided by known mechanisms and trends from previous studies, indicating multicontinuum theory as a potentially effective tool in predicting the long-term performance of composites in marine environments.
Strong gas-mineral interactions or slow adsorption kinetics require a molecular-level understanding of both adsorption and diffusion for these interactions to be properly described in transport models. In this combined molecular simulation and experimental study, noble gas adsorption and mobility is investigated in two naturally abundant zeolites whose pores are similar in size (clinoptilolite) and greater than (mordenite) the gas diameters. Simulated adsorption isotherms obtained from grand canonical Monte Carlo simulations indicate that both zeolites can accommodate even the largest gas (Rn). However, gas mobility in clinoptilolite is significantly hindered at pore-limiting window sites, as seen from molecular dynamics simulations in both bulk and slab zeolite models. Experimental gas adsorption isotherms for clinoptilolite confirm the presence of a kinetic barrier to Xe uptake, resulting in the unusual property of reverse Kr/Xe selectivity. Finally, a kinetic model is used to fit the simulated gas loading profiles, allowing a comparison of trends in gas diffusivity in the zeolite pores.
Physics-constrained machine learning is emerging as an important topic in the field of machine learning for physics. One of the most significant advantages of incorporating physics constraints into machine learning methods is that the resulting model requires significantly less data to train. By incorporating physical rules into the machine learning formulation itself, the predictions are expected to be physically plausible. Gaussian process (GP) is perhaps one of the most common methods in machine learning for small datasets. In this paper, we investigate the possibility of constraining a GP formulation with monotonicity on three different material datasets, where one experimental and two computational datasets are used. The monotonic GP is compared against the regular GP, where a significant reduction in the posterior variance is observed. The monotonic GP is strictly monotonic in the interpolation regime, but in the extrapolation regime, the monotonic effect starts fading away as one goes beyond the training dataset. Imposing monotonicity on the GP comes at a small accuracy cost, compared to the regular GP. The monotonic GP is perhaps most useful in applications where data are scarce and noisy, and monotonicity is supported by strong physical evidence.
Many important engineering and scientific applications such as cement slurries, foams, crude oil, and granular avalanches involve the concept of yield stress. Therefore, modeling yield stress fluids in different flow configurations, including the accurate prediction of the yield surface, is important. In this paper, we present a computational model based on the finite element method to study the flow of yield stress fluids in a thin mold and compare the results with data from flow visualization experiments. We use the level set method to describe the interface between the filling fluid and air. We use polypropylene glycol as a model Newtonian fluid and Carbopol for the model yield stress fluid, as the Carbopol solution demonstrates yielding without thixotropy. To describe the yielding and shear-thinning behavior, we use a generalized Newtonian constitutive equation with a Bingham–Carreau–Yasuda form. We compare the results obtained from the mold filling experiments with the results from the three-dimensional (3D) model and from a reduced-order Hele-Shaw (HS) model that is two-dimensional, including the effect of shear-thinning along the thin direction only approximately. We show that both the 3D and the HS model can capture the experimental meniscus shape reasonably well for all the fluids considered at three different flow rates. This indicates that the shape evolution is insensitive to the dimensionality of the model. However, the viscosity and yield surfaces predicted by the 3D and HS models are different. The HS model underestimates the high viscosity and unyielded regions compared to the estimation by the 3D model.
Polymer concrete (PC) has been used to replace cement concrete when harsh service conditions exist. Polymers have a high carbon footprint when considering their life cycle analysis, and with increased climate change concerns and the need to reduce greenhouse gas emission, bio-based polymers could be used as a sustainable alternative binder to produce PC. This paper examines the development and characterization of a novel bio-polymer concrete (BPC) using bio-based polyurethane used as the binder in lieu of cement, modified with benzoic acid and carboxyl-functionalized multi-walled carbon nanotubes (MWCNTs). The mechanical performance, durability, microstructure, and chemical properties of BPC are investigated. Moreover, the effect of the addition of benzoic acid and MWCNTs on the properties of BPC is studied. The new BPC shows relatively low density, appreciable compressive strength between 20–30 MPa, good tensile strength of 4 MPa, and excellent durability resistance against aggressive environments. The new BPC has a low carbon footprint, 50% lower than ordinary Portland cement concrete, and can provide a sustainable concrete alternative in infrastructural applications.
In order to better understand the complex pooling and vaporization of a liquid hydrogen spill, Sandia National Laboratories is conducting a highly instrumented, controlled experiment inside their Shock Tube Facility. Simulations were run before the experiment to help with the planning of experimental conditions, including sensor placement and cross wind velocity. This paper describes the modeling used in this planning process and its main conclusions. Sierra Suite’s Fuego, an in-house computational fluid dynamics code, was used to simulate a RANS model of a liquid hydrogen spill with five crosswind velocities: 0.45, 0.89, 1.34, 1.79, and 2.24 m/s. Two pool sizes were considered: a diameter of 0.85 m and a diameter of 1.7. A grid resolution study was completed on the smaller pool size with a 1.34 m/s crosswind. A comparison of the length and height of the plume of flammable hydrogen vaporizing from the pool shows that the plume becomes longer and remains closer to the ground with increasing wind speed. The plume reaches the top of the facility only in the 0.45 m/s case. From these results, we concluded that it will be best for the spacing and location of the concentration sensors to be reconfigured for each wind speed during the experiment.
Metal-oxide composites are commonly used in high temperature environments for their thermal stability and high melting points. Commonly employed with refractory oxides or carbides such as ZrC and HfC, these materials may be improved with the use of a low density, high melting point ceramic such as CeO2. In this work, the consolidation of W-CeO2 metal matrix composites in the high CeO2 concentration regime is explored. The CeO2 concentrations of 50, 33, and 25 wt.%, the CeO2 particle size from nanometer to micrometer, and various hot isostatic pressing temperatures are investigated. Decreasing the CeO2 concentration is observed to increase the composite density and increase the Vickers hardness. The CeO2 oxidation state is observed to be a combination of Ce3+ and Ce4+, which is hypothesized to contribute to the porosity of the composites. The hardness of the metal-oxide composite can be improved more than 2.5 times compared to pure W processed by the same route. This work offers processing guidelines for further consolation of oxide-doped W composites.
The production of biochar from biomass and industrial wastes provides both environmental and economic sustainability. An effective way to ensure the sustainability of biochar is to produce high value-added activated carbon. The desirable characteristic of activated carbon is its high surface area for efficient adsorption of contaminants. Feedstocks can include a number of locally available materials with little or negative value, such as orchard slash and crop residue. In this context, it is necessary to determine and know the conversion effects of the feedstocks to be used in the production of activated carbon. In the study conducted for this purpose; several samples (piñon wood, pecan wood, hardwood, dried grass, Wyoming coal dust, Illinois coal dust, Missouri coal dust, and tire residue) of biomass and industrial waste products were investigated for their conversion into activated carbon. Small samples (approximately 0.02 g) of the feedstocks were pyrolyzed under inert or mildly oxidizing conditions in a thermal analyzer to determine their mass loss as a function of temperature and atmosphere. Once suitable conditions were established, larger quantities (up to 0.6 g) were pyrolyzed in a tube furnace and harvested for characterization of their surface area and porosity via gas sorption analysis. Among the samples used, piñon wood gave the best results, and pyrolysis temperatures between 600 and 650 °C gave the highest yield. Slow pyrolysis or hydrothermal carbonization have come to the fore as recommended production methods for the conversion of biochar, which can be produced from biomass and industrial wastes, into activated carbon.
Evaluation of high temperature (HT) microcontrollers for Geothermal applications. The goal was to evaluate the limited microcontrollers capable of operating above 200 degrees Celsius. The knowledge obtained from this research has been documented for contributing to future development of high temperature geothermal tools. Along with this research, a HT data link for communications was developed and evaluated, also for HT geothermal tools.
Here we report on AlGaN high electron mobility transistor (HEMT)-based logic development, using combined enhancement- and depletion-mode transistors to fabricate inverters with operation from room temperature up to 500°C. Our development approach included: (a) characterizing temperature-dependent carrier transport for different AlGaN HEMT heterostructures, (b) developing a suitable gate metal scheme for use in high temperatures, and (c) over-temperature testing of discrete devices and inverters. Hall mobility data (from 30°C to 500°C) revealed the reference GaN-channel HEMT experienced a 6.9x reduction in mobility, whereas the AlGaN channel HEMTs experienced about a 3.1x reduction. Furthermore, a greater aluminum contrast between the barrier and channel enabled higher carrier densities in the two-dimensional electron gas for all temperatures. The combination of reduced variation in mobility with temperature and high sheet carrier concentration showed that an Al-rich AlGaN-channel HEMT with a high barrier-to-channel aluminum contrast is the best option for an extreme temperature HEMT design. Three gate metal stacks were selected for low resistivity, high melting point, low thermal expansion coefficient, and high expected barrier height. The impact of thermal cycling was examined through electrical characterization of samples measured before and after rapid thermal anneal. The 200-nm tungsten gate metallization was the top performer with minimal reduction in drain current, a slightly positive threshold voltage shift, and about an order of magnitude advantage over the other gates in on-to-off current ratio. After incorporating the tungsten gate metal stack in device fabrication, characterization of transistors and inverters from room temperature up to 500°C was performed. The enhancement-mode (e-mode) devices’ resistance started increasing at about 200°C, resulting in drain current degradation. This phenomenon was not observed in depletion-mode (d-mode) devices but highlights a challenge for inverters in an e-mode driver and d-mode load configuration.
Mcglone, Joe F.; Ghadi, Hemant; Cornuelle, Evan; Armstrong, Andrew A.; Burns, George B.; Feng, Zixuan; Uddin Bhuiyan, A.F.M.A.; Zhao, Hongping; Arehart, Aaron R.; Ringel, Steven A.
The impact of 1.8 MeV proton irradiation on metalorganic chemical vapor deposition grown (010) β-Ga2O3 Schottky diodes is presented. It is found that after a 10.8 × 10 13 cm - 2 proton fluence the Schottky barrier height of (1.40 ± 0.05 eV) and the ideality factor of (1.05 ± 0.05) are unaffected. Capacitance-voltage extracted net ionized doping curves indicate a carrier removal rate of 268 ± 10 cm - 1. The defect states responsible for the observed carrier removal are studied through a combination of deep level transient and optical spectroscopies (DLTS/DLOS) as well as lighted capacitance-voltage (LCV) measurements. The dominating effect on the defect spectrum is due to the EC-2.0 eV defect state observed in DLOS and LCV. This state accounts for ∼ 75% of the total trap introduction rate and is the primary source of carrier removal from proton irradiation. Of the DLTS detected states, the EC-0.72 eV state dominated but had a comparably smaller contribution to the trap introduction. These two traps have previously been correlated with acceptor-like gallium vacancy-related defects. Several other trap states at EC-0.36, EC-0.63, and EC-1.09 eV were newly detected after proton irradiation, and two pre-existing states at EC-1.2 and EC-4.4 eV showed a slight increase in concentration after irradiation, together accounting for the remainder of trap introduction. However, a pre-existing trap at EC-0.40 eV was found to be insensitive to proton irradiation and, therefore, is likely of extrinsic origin. The comprehensive defect characterization of 1.8 MeV proton irradiation damage can aid the modeling and design for a range of radiation tolerant devices.
Mohottalalage, Supun S.; Kosgallana, Chathurika; Meedin, Shalika; Connor, Gary S.'.; Grest, Gary S.; Perahia, Dvora
Ionizable polymers form dynamic networks with domains controlled by two distinct energy scales, ionic interactions and van der Waals forces; both evolve under elongational flows during their processing into viable materials. A molecular level insight of their nonlinear response, paramount to controlling their structure, is attained by fully atomistic molecular dynamics simulations of a model ionizable polymer, polystyrene sulfonate. As a function of increasing elongational flow rate, the systems display an initial elastic response, followed by an ionic fraction-dependent strain hardening, stress overshoot, and eventually strain-thinning. As the sulfonation fraction increases, the chain elongation becomes more heterogeneous. Finally, flow-driven ionic assembly dynamics that continuously break and reform control the response of the system.
Automation of rate-coefficient calculations for gas-phase organic species became possible in recent years and has transformed how we explore these complicated systems computationally. Kinetics workflow tools bring rigor and speed and eliminate a large fraction of manual labor and related error sources. In this paper we give an overview of this quickly evolving field and illustrate, through five detailed examples, the capabilities of our own automated tool, KinBot. We bring examples from combustion and atmospheric chemistry of C-, H-, O-, and N-atom-containing species that are relevant to molecular weight growth and autoxidation processes. The examples shed light on the capabilities of automation and also highlight particular challenges associated with the various chemical systems that need to be addressed in future work.
Bifurcations are commonly encountered during force controlled swept and stepped sine testing of nonlinear structures, which generally leads to the so-called jump-down or jump-up phenomena between stable solutions. There are various experimental closed-loop control algorithms, such as control-based continuation and phase-locked loop, to stabilize dynamical systems through these bifurcations, but they generally rely on specialized control algorithms that are not readily available with many commercial data acquisition software packages. A recent method was developed to experimentally apply sequential continuation using the shaker voltage that can be readily deployed using commercially available software. By utilizing the stabilizing effects of electrodynamic shakers and the force dropout phenomena in fixed frequency voltage control sine tests, this approach has been demonstrated to stabilize the unstable branch of a nonlinear system with three branches, allowing for three multivalued solutions to be identified within a specific frequency bandwidth near resonance. Recent testing on a strongly nonlinear system with vibro-impact nonlinearity has revealed jumping behavior when performing sequential continuation along the voltage parameter, like the jump phenomena seen during more traditional force controlled swept and stepped sine testing. Here, this paper investigates the stabilizing effects of an electrodynamic shaker on strongly nonlinear structures in fixed frequency voltage control tests using both numerical and experimental methods. The harmonic balance method is applied to the coupled shaker-structure system with an electromechanical model to simulate the fixed voltage control tests and predict the stabilization for different parameters of the model. The simulated results are leveraged to inform the design of a set of experiments to demonstrate the stabilization characteristics on a fixture-pylon assembly with a vibro-impact nonlinearity. Through numerical simulation and experimental testing on two different strongly nonlinear systems, the various parameters that influence the stability of the coupled shaker-structure are revealed to better understand the performance of fixed frequency voltage control tests.
Vibrational spectroscopy is a nondestructive technique commonly used in chemical and physical analyses to determine atomic structures and associated properties. However, the evaluation and interpretation of spectroscopic profiles based on human-identifiable peaks can be difficult and convoluted. To address this challenge, we present a reliable protocol based on supervised manifold learning techniques meant to connect vibrational spectra to a variety of complex and diverse atomic structure configurations. As an illustration, we examined a large database of virtual vibrational spectroscopy profiles generated from atomistic simulations for silicon structures subjected to different stress, amorphization, and disordering states. We evaluated representative features in those spectra via various linear and nonlinear dimensionality reduction techniques and used the reduced representation of those features with decision trees to correlate them with structural information unavailable through classical human-identifiable peak analysis. We show that our trained model accurately (over 97% accuracy) and robustly (insensitive to noise) disentangles the contribution from the different material states, hence demonstrating a comprehensive decoding of spectroscopic profiles beyond classical (human-identifiable) peak analysis.
Hydrocarbon polymers are used in a wide variety of practical applications. In the field of dynamic compression at extreme pressures, these polymers are used at several high energy density (HED) experimental facilities. One of the most common polymers is poly(methyl methacrylate) or PMMA, also called Plexiglass® or Lucite®. Here, we present high-fidelity, hundreds of GPa range experimental shock compression data measured on Sandia's Z machine. We extend the principal shock Hugoniot for PMMA to more than threefold compression up to 650 GPa and re-shock Hugoniot states up to 1020 GPa in an off-Hugoniot regime, where experimental data are even sparser. These data can be used to put additional constraints on tabular equation of state (EOS) models. The present results provide clear evidence for the need to re-examine the existing tabular EOS models for PMMA above ∼120 GPa as well as perhaps revisit EOSs of similar hydrocarbon polymers commonly used in HED experiments investigating dynamic compression, hydrodynamics, or inertial confinement fusion.
Accurately modeling the impact force used in the analysis of loosely constrained cantilevered pipes conveying fluid is imperative. If little information is known of the motion-limiting constraints used in experiments, the analysis of the system may yield inaccurate predictions. Here in this work, multiple forcing representations of the impact force are defined and analyzed for a cantilevered pipe that conveys fluid. Depending on the representation of the impact force, the dynamics of the pipe can vary greatly when only the stiffness of the constraints is known from experiments. Three gap sizes of the constraints are analyzed, and the representation of the impact force used to analyze the system is found to significantly affect the response of the pipe at each gap size. An investigation on the effects of the vibro-impact force representation is performed through using basin of attraction analysis and nonlinear characterization of the system’s response.