Computer Methods in Applied Mechanics and Engineering
Dingreville, Remi P.; Francis, Noah M.; Pourahmadian, Fatemeh; Lebensohn, Ricardo A.
This work presents a spectral micromechanical formulation for obtaining the full-field and homogenized response of elastic micropolar composites. The algorithm relies on a coupled set of convolution integral equations for the micropolar strains, where periodic Green’s operators associated with a linear homogeneous reference medium are convolved with functions of the Cauchy and couple stress fields that encode the material’s heterogeneity, as well as any potential material nonlinearity. Such convolution integral equations take an algebraic form in the reciprocal Fourier space that can be solved iteratively. In this vein, the fast Fourier transform (FFT) algorithm is leveraged to accelerate the numerical solution, resulting in a mesh-free formulation in which the periodic unit cell representing the heterogeneous material can be discretized by a regular grid of pixels in two dimensions (or voxels in three dimensions). For verification, the numerical solutions obtained with the micropolar FFT solver are compared with analytical solutions for a matrix with a dilute circular inclusion subjected to plane strain loading. The developed computational framework is then used to study length-scale effects and effective (micropolar) moduli of composites with various topological configurations.
A design study was conducted at the National Solar Thermal Test Facility (NSTTF) at Sandia National Laboratories in Albuquerque, NM with the objective of identifying the technical readiness level, performance limits, capital and O&M costs, and expected thermal losses of particle handling and conveyance components in a particle-based CSP plant. Key findings indicated that skips and high temperature particle conveyance technology are available for moving particles up to 615° ± 25° C. This limits the use of mechanical conveyance above the heat exchanger and suggests vertical integration of the hot storage bin and heat exchanger to facilitate direct gravity fed handling of particles. Skip rails and support structures add significant cost and must be factored into cost analysis. Chutes can be a low-cost option for particle handling but uncertainties in tower costs make it difficult to know whether they can be cost effective in areas above the receiver.
The Sandia National Laboratories (SNL) National Solar Thermal Test Facility (NSTTF) conducted efficacy testing on a shut-off isolation valve for use with molten ternary chloride salt. A ball valve was tested under controlled N2 ullage gas pressure and connected with flanged fittings that featured a spiral-wound gasket. The valve assembly consisted of boronized nickel coated SS316 components, with design features that greatly reduce the cost of overall valve assembly. Testing results showed that the valve did not leak, and post-test analysis demonstrated that the ball, seat, packing, and body all survived both the heat loads and the relative corrosive environment. Spiral-wound gaskets for flanged connections used in the system also functioned nominally, with no leaks or signs of failures during post-test analysis. However, testing was ultimately forced to rapidly stop after testing between 500-530°C as the actuator used on the valve failed in the heat, preventing the valve from sealing in the closed position. In addition, salt plugs and salt vapor plating also prevented the test from continuing.
The need for reliable, cost-effective, utility scale energy storage that is universally applicable across different regions is becoming evident with the global transition towards non-polluting renewable energy resources. The operations and management of these energy storage technologies introduces a unique challenge that is inherently different from the conventional energy storage in the form of fossil fuel. The investigation into the business model, value proposition and economic viability of a utility scale thermal energy storage was part of a program sponsored by the United States Department of Energy, called Energy I-Corps. During this program, the project team reached out to a series of industry stakeholders to conduct interviews on the topic of thermal energy storage for utility scale power generation. Specific focus was placed on the business model based on the market needs in the context of the power grid in the United States. The utilization and re-use of infrastructure at existing thermo-electric power plants yielded the most viable business model for the implementation of the form of thermal energy storage discussed here.
A third-generation chloride salt tank system was designed for a 1 MWth pilot-scale system to be investigated at the National Solar Thermal Test Facility (NSTTF) in Albuquerque, NM, USA. This prototype Gen 3, concentrating solar power (CSP) system was designed to facilitate a minimum of 6 hrs. of thermal energy storage (TES) with operational nominal temperatures of 500°C and 720°C for a cold and hot tank respectively. For this investigation, the researchers developed steady and transient computational fluid mechanics (CFD) circulation models to assess thermal-fluid behavior within the tanks, and their respective interactions with environmental heat transfer. The models developed for this novel CSP system design included unique chloride molten salt thermodynamic properties and correlations. The results of this investigation suggest thermal gradients for the steady flow model less 1oC with overall circulation velocities as high as approximately 2.1 m/s. Higher steady flow rates of salt passing into and out of the tanks resulted in smaller thermal gradients than the slower flow rates as the molten salt mixes better (an increase of around 120% in the heat transfer coefficient) at the higher velocities associated with the higher flow rate. The port spacing of 3.85 m was found to have a highly uniform temperature distribution. For the unsteady model, nitrogen flow was found to become appreciably steady after approximately 10 minutes, and resultant molten salt flow was found to increase slowly as the overall salt level rose.
This paper summarizes findings from a small, mixed-method research study examining industry perspectives on the potential for new forms of automation to invigorate the concentrating solar power (CSP) industry. In Fall 2021, the Solar Energy Technologies Office (SETO) of the United States Department of Energy (DOE) funded Sandia National Laboratories to elicit industry stakeholder perspectives on the potential role of automated systems in CSP operations. We interviewed eleven CSP professionals from five countries, using a combination of structured and open comment response modes. Respondents indicated a preference for automated systems that support heliostat manufacturing and installation, calibration, and responsiveness to shifting weather conditions. This pilot study demonstrates the importance of engaging industry stakeholders in discussions of technology research and development, to promote adoptable, useful innovation.
Falling particle receivers are a promising receiver design to couple with particle-based concentrating solar power to help meet future levelized cost of electricity targets in next generation systems. The thermal performance of receivers is critical to the economics of the overall system, and accurate models of particle receivers are necessary to predict the performance in all conditions. A model validation study was performed using falling particle receiver data recently collected at the National Solar Thermal Test Facility at Sandia National Laboratories. The particle outlet temperature, the thermal efficiency of the receiver, and the wind speed and direction around the receiver were measured in 26 steady-state experiments and compared to a corresponding receiver model. The results of this study showed improved agreement with the experimental data over past validation efforts but did not fully meet all predefined validation metrics. Future model improvements were identified to continue to strengthen the modeling capabilities.
The United Sates Department of Energy (DOE) Generation 3 Concentrated Solar Power (CSP) program is interested in higher efficiency power systems at lower costs, potentially with systems utilizing chloride molten salts. Ternary chloride molten salts are corrosive and need to be held at high temperatures to achieve higher power system efficiencies. However, materials and cost of manufacturing of such a facility can be very expensive, particularly using exotic materials that are not always readily available. Materials that can withstand the harsh corrosive and thermal-mechanical environments of high-temperature molten salt systems (>700 ℃) are needed. High temperature systems offer greater thermodynamic efficiency but must also make cost efficient use of corrosion-resistant alloys. To ensure reliable high-performance operation for molten salt power plant designs confidence in materials compatibility with CSP Gen 3 halide salts must be established. This paper will present an analysis of Inconel 625 as an alternative to the costly Haynes 230 at 760℃ for 500 hours. Both metals were tested in an unaltered state as well as a homogenous weld. Each sample was weighed pre- and post-test, with a final composition analysis using Scanning Electron Microscopy (SEM) and Energy Dispersive X-Ray Spectroscopy (EDS). Preliminary findings suggest that Haynes 230 outperformed Inconel 625, but more research at longer durations, 1,000 hours will be required for full reliable assessment.
The purpose of this proposal is to design a new integral critical experiment to investigate the effects of beryllium oxide and high assay low-enriched uranium fuels. this proposal considers using several existing resources at Sandia: (1) the Critical Experiments (SCX) facility and water tank, (2) spare UO2-BeO fuel for the Annular Core Research Reactor (ACRR), and 7uPCX fuel rods from previous benchmark experiments.
The effective management of plastic waste streams to prevent plastic land and water pollution is a growing problem that is also one of the most important challenges in polymer science today. Polymer materials that are stable over their lifetime and can also be cheaply recycled or repurposed as desired could more easily be diverted from waste streams. However, this is difficult for most commodity plastics. It is especially difficult to conceive this with intractable, cross-linked polymers such as rubbers. In this work, we explore the utility of microencapsulated Grubbs’ catalysts for the in-situ depolymerization and reprocessing of polybutadiene (PB) rubber. Second-generation Hoveyda-Grubbs catalyst (HG2) contained within glassy thermoplastic microspheres can be dispersed in PB rubber below the microsphere’s glass transition temperature (Tg) without adverse depolymerization, evidenced by rubber with and without these microspheres obtaining similar shear storage moduli of ≈16 and ≈28 kPa, respectively. The thermoplastic’s Tg can be used to tune the depolymerization temperature, via release of HG2 into the rubber matrix. For example, using poly(lactic acid) (PLA) vs polysulfone results in an 85 and 162 °C depolymerization temperature, respectively. Liquefaction of rubber to a mixture of small molecules and oligomers is demonstrated using a 0.01 mol % catalyst loading using PLA as the encapsulant. At that same catalyst loading, depolymerization occurs to a greater extent in comparison to two ex-situ approaches, including a conventional solvent-assisted method, where it occurs at roughly twice the extent at each given catalyst loading. In addition, depolymerization of the microsphere-loaded rubbers was demonstrated for samples stored under nitrogen for 23 days. Lastly, we show that the depolymerized products can be reprocessed back into solid rubber with a shear storage modulus of ≈32 kPa. Thus, we envision that this approach could be used to recycle and reuse cross-linked rubbers at the end of their product lifetime.
In this paper, an approach for 3D plasma structure diagnostics using tomographic optical emission spectroscopy (Tomo-OES) of a nanosecond pulsed atmospheric pressure plasma jet (APPJ) is presented. In contrast to the well-known Abel inversion, Tomo-OES does not require cylindrical symmetry to recover 3D distributions of plasma light emission. Instead, many 2D angular projections are measured with intensified cameras and the multiplicative algebraic reconstruction technique is used to recover the 3D distribution of light emission. This approach solves the line-of-sight integration problem inherent to optical diagnostics, allowing recovery of localized OES information within the plasma that can be used to better infer plasma parameters within complex plasma structures. Here, Tomo-OES was applied to investigate an APPJ operated with helium in ambient air and impinging on planar and structured dielectric surfaces. Surface charging caused the guided streamer from the APPJ to transition to a surface ionization wave (SIW) that propagated along the surface. The SIW experienced variable geometrical and electrical material properties as it propagated, leading to 3D configurations that were non-symmetric and spatially complex. Light emission from He, N 2 + , and N2 were imaged at ten angular projections and the respective time-resolved 3D emission distributions in the plasma were then reconstructed. The spatial resolution of each tomographic reconstruction was 7.4 µm and the temporal resolution was 5 ns, sufficient to observe the guided streamer and the effects of the structured surface on the SIW. Emission from He showed the core of the jet and emission from N 2 + and N2 indicated effects of entrainment of ambient air. Penning ionization of N2 created a ring or outer layer of N 2 + that spatially converged to form the ‘plasma bullet’ or spatially diverged across a surface as part of a SIW. The SIW entered trenches of size 150 µm, leading to decreases in plasma light emission in regions above the trenches. The plasma light emission was higher in some regions with trenches, possibly due to effects of field enhancement.
The methodology described in this article enables a type of holistic fleet optimization that simultaneously considers the composition and activity of a fleet through time as well as the design of individual systems within the fleet. Often, real-world system design optimization and fleet-level acquisition optimization are treated separately due to the prohibitive scale and complexity of each problem. This means that fleet-level schedules are typically limited to the inclusion of predefined system configurations and are blind to a rich spectrum of system design alternatives. Similarly, system design optimization often considers a system in isolation from the fleet and is blind to numerous, complex portfolio-level considerations. In reality, these two problems are highly interconnected. To properly address this system-fleet design interdependence, we present a general method for efficiently incorporating multi-objective system design trade-off information into a mixed-integer linear programming (MILP) fleet-level optimization. This work is motivated by the authors' experience with large-scale DOD acquisition portfolios. However, the methodology is general to any application where the fleet-level problem is a MILP and there exists at least one system having a design trade space in which two or more design objectives are parameters in the fleet-level MILP.
We describe a direct magneto-optical approach to measuring the magnetic field driven by a narrow pulse width (<10 ns), 20 kA electrical current flow in the transmission line of a high energy pulsed power accelerator. The magnetic field and electrical current are among the most important operating parameters in a pulsed power accelerator and are critical to understanding the properties of the radiation output. However, accurately measuring these fields and electrical currents using conventional pulsed power diagnostics is difficult due to the strength of ionizing radiation and electromagnetic interference. Our approach uses a fiber coupled laser beam with a rare earth element sensing crystal sensor that is highly resistant to electromagnetic interference and does not require external calibration. Here, we focus on device theory, operating parameters, results from an experiment on a high energy pulsed power accelerator, and comparison to a conventional electrical current shunt sensor.
Advancements in space exploration and sample return technology present a unique opportunity to leverage sample return capsules (SRCs) towards studying atmospheric entry of meteoroids and asteroids. Specifically engineered for the secure transport of valuable extraterrestrial samples from interplanetary space to Earth, SRCs offer unexpected benefits that reach beyond their intended purpose. As SRCs enter the Earth’s atmosphere at hypervelocity, they are analogous to naturally occurring meteoroids and thus, for all intents and purposes, can be considered artificial meteors. Furthermore, SRCs are capable of generating shockwaves upon reaching the lower transitional flow regime, and thus can be detected by strategically positioned geophysical instrumentation. NASA’s OSIRIS-REx (Origins, Spectral Interpretation, Resource Identification, and Security-Regolith Explorer) SRC is one of only a handful of artificial objects to re-enter the Earth’s atmosphere from interplanetary space since the end of the Apollo era and it will provide an unprecedented observational opportunity. This review summarizes past infrasound and seismic observational studies of SRC re-entries since the end of the Apollo era and presents their utility towards the better characterization of meteoroid flight through the atmosphere.
Tritium exhibits unique environmental behavior because of its potential interactions with water and organic substances. Modeling the environmental consequences of tritium releases can be relatively complex and thus an evaluation of MACCS is needed to understand what updates, if any, are needed in MACCS to account for the behavior of tritium. We examine documented tritium releases and previous benchmarking assessments to perform a model intercomparison between MACCS and state-of-practice tritium-specific codes UFOTRI and ETMOD to quantify the difference between MACCS and state of practice models for assessing tritium consequences. Additionally, information to assist an analyst in judging whether a postulated tritium release is likely to lead to significant doses is provided.
Hydrogen diffusion in metals and alloys plays an important role in the discovery of new materials for fuel cell and energy storage technology. While analytic models use hand-selected features that have clear physical ties to hydrogen diffusion, they often lack accuracy when making quantitative predictions. Machine learning models are capable of making accurate predictions, but their inner workings are obscured, rendering it unclear which physical features are truly important. To develop interpretable machine learning models to predict the activation energies of hydrogen diffusion in metals and random binary alloys, we create a database for physical and chemical properties of the species and use it to fit six machine learning models. Our models achieve root-mean-squared errors between 98-119 meV on the testing data and accurately predict that elemental Ru has a large activation energy, while elemental Cr and Fe have small activation energies. By analyzing the feature importances of these fitted models, we identify relevant physical properties for predicting hydrogen diffusivity. While metrics for measuring the individual feature importances for machine learning models exist, correlations between the features lead to disagreement between models and limit the conclusions that can be drawn. Instead grouped feature importance, formed by combining the features via their correlations, agree across the six models and reveal that the two groups containing the packing factor and electronic specific heat are particularly significant for predicting hydrogen diffusion in metals and random binary alloys. This framework allows us to interpret machine learning models and enables rapid screening of new materials with the desired rates of hydrogen diffusion.
Freeplay is a common type of piecewise-smooth nonlinearity in dynamical systems, and it can cause discontinuity-induced bifurcations and other behaviors that may bring about undesirable and potentially damaging responses. Prior research has focused on piecewise-smooth systems with two or three distinct regions, but less attention is devoted to systems with more regions (i.e., multi-segmented systems). In this work, numerical analysis is performed on a dynamical system with multi-segmented freeplay, in which there are four stiffness transitions and five distinct regions in the phase space. The effects of the multi-segmented parameters are studied through bifurcation diagram evolution along with induced multi-stable behavior and different bifurcations. These phenomena are interrogated through various tools, such as harmonic balance, basins of attraction, phase planes, and Poincaré section analysis. Results show that among the three multi-segmented parameters, the asymmetry has the strongest effect on the response of the system.
Nonlinear modeling and optimization is a valuable tool for aiding decisions by engineering practitioners, but programming an optimization problem based on a complex electrical, mechanical, or chemical process is a time-consuming and error-prone activity. Therefore, there is a need for model analysis and debugging tools that can detect and diagnose modeling errors. One such tool is the Dulmage–Mendelsohn decomposition, which identifies structurally under- and over-determined subsets in systems of equations and variables by partitioning the bipartite graph of the system. This work provides the necessary background to understand the Dulmage–Mendelsohn decomposition and its application to the analysis of nonlinear optimization problems, demonstrates its use in diagnosing a variety of modeling errors, and introduces software implementations for analyzing nonlinear optimization problems in the Pyomo and JuMP algebraic modeling languages.