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.
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.
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 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.
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.
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.
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.
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.
Machine learning models are promising as surrogates in optimization when replacing difficult to solve equations or black-box type models. This work demonstrates the viability of linear model decision trees as piecewise-linear surrogates in decision-making problems. Linear model decision trees can be represented exactly in mixed-integer linear programming (MILP) and mixed-integer quadratic constrained programming (MIQCP) formulations. Furthermore, they can represent discontinuous functions, bringing advantages over neural networks in some cases. We present several formulations using transformations from Generalized Disjunctive Programming (GDP) formulations and modifications of MILP formulations for gradient boosted decision trees (GBDT). We then compare the computational performance of these different MILP and MIQCP representations in an optimization problem and illustrate their use on engineering applications. We observe faster solution times for optimization problems with linear model decision tree surrogates when compared with GBDT surrogates using the Optimization and Machine Learning Toolkit (OMLT).
Approximating differential operators defined on two-dimensional surfaces is an important problem that arises in many areas of science and engineering. Over the past ten years, localized meshfree methods based on generalized moving least squares (GMLS) and radial basis function finite differences (RBF-FD) have been shown to be effective for this task as they can give high orders of accuracy at low computational cost, and they can be applied to surfaces defined only by point clouds. However, there have yet to be any studies that perform a direct comparison of these methods for approximating surface differential operators (SDOs). The first purpose of this work is to fill that gap. For this comparison, we focus on an RBF-FD method based on polyharmonic spline kernels and polynomials (PHS+Poly) since they are most closely related to the GMLS method. Additionally, we use a relatively new technique for approximating SDOs with RBF-FD called the tangent plane method since it is simpler than previous techniques and natural to use with PHS+Poly RBF-FD. The second purpose of this work is to relate the tangent plane formulation of SDOs to the local coordinate formulation used in GMLS and to show that they are equivalent when the tangent space to the surface is known exactly. The final purpose is to use ideas from the GMLS SDO formulation to derive a new RBF-FD method for approximating the tangent space for a point cloud surface when it is unknown. For the numerical comparisons of the methods, we examine their convergence rates for approximating the surface gradient, divergence, and Laplacian as the point clouds are refined for various parameter choices. We also compare their efficiency in terms of accuracy per computational cost, both when including and excluding setup costs.
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.
This is a SAND Report on cross-correlation data collected at the Redmond Salt Mine. It discusses methods, as well as temporal variability and energy characteristics of the cross-correlation data.
Improved power take-off (PTO) controller design for wave energy converters is considered a critical component for reducing the cost of energy production. However, the device and control design process often remains sequential, with the space of possible final designs largely reduced before the controller has been considered. Control co-design, whereby the device and control design are considered concurrently, has resulted in improved designs in many industries, but remains rare in the wave energy community. In this paper we demonstrate the use of a new open-source code, WecOptTool, for control co-design of wave energy converters, with the aim to make the co-design approach more accessible and accelerate its adoption. Additionally, we highlight the importance of designing a wave energy converter to maximize electrical power, rather than mechanical power, and demonstrate the co-design process while modeling the PTO's components (i.e., drive-train and generator, and their dynamics). We also consider the design and optimization of causal fixed-structure controllers. The demonstration presented here considers the PTO design problem and finds the optimal PTO drive-train that maximizes annual electrical power production. The results show a 22% improvement in the optimal controller and drive-train co-design over the optimal controller for the nominal, as built, device design.
Coe, Ryan G.; Lee, Jantzen; Bacelli, Giorgio B.; Spencer, Steven; Dullea, Kevin; Plueddemann, Albert J.; Buffitt, Derek; Reine, John; Peters, Donald; Spinneken, Johannes; Hamilton, Andrew; Sabet, Sahand; Husain, Salman; Jenne, Dale (Scott); Korde, Umesh; Muglia, Mike; Taylor, Trip; Wade, Eric
The “Pioneer WEC” project is targeted at developing a wave energy generator for the Coastal Surface Mooring (CSM) system within the Ocean Observatories Initiative (OOI) Pioneer Array. The CSM utilizes solar photovoltaic and wind generation systems, along with rechargeable batteries, to power multiple sensors on the buoy and along the mooring line. This approach provides continuous power for essential controller functions and a subset of instruments, and meets the full power demand roughly 70% of the time. Sandia has been tasked with designing a wave energy system to provide additional electrical power and bring the CSM up-time for satisfying the full-power demand to 100%. This project is a collaboration between Sandia and Woods Hole Oceanographic Institution (WHOI), along with Evergreen Innovations, Monterey Bay Aquarium Research Institute (MBARI), Eastern Carolina University (ECU), Johns Hopkins University (JHU), and the National Renewable Energy Laboratory (NREL). This report captures Phase I of an expected two phase project and presents project scoping and concept design results. phase project and presents project scoping and concept design results.
This report provides an overview of some code written to streamline aspects of the finite element simulation process using Sierra/SD. Some specific examples of common simulation tasks are provided.
Large rocket motors may violently explode when exposed to accidental fires. Even hot metal fragments from a nearby accident may penetrate the propellant and ultimately cause thermal ignition. A mechanistic understanding of heated propellants leading to thermal runaway is a major unsolved problem. Here we show that thermal ignition in propellants can be predicted using a universal cookoff model coupled to a micromechanics pressurization model. Our model predicts the time to thermal ignition in cookoff experiments with variable headspace volumes. We found that experiments with headspace volumes are more prone to deformation which distorts pores and causes increased permeability when the propellant expands into this headspace. Delayed ignition with larger headspace volume correlates with lower headspace pressures during decomposition. We found that our predictions matched experimental measurements best when the initial propellant was impermeable to gas flow rather than being permeable. Similar behavior is expected with other energetic materials with rubbery binders. Our model is validated using data from a separate laboratory. We also present an uncertainty analysis using Latin Hypercube Sampling (LHS) of thermal ignition caused by a steel fragment embedded in the propellant.
The Sound Fixing and Ranging (SOFAR) channel in the ocean allows for low frequency sound to travel thousands of kilometers, making it particularly useful for detecting underwater nuclear explosions. Suggestions that an elevated SOFAR-like channel should exist in the stratosphere date back over half a century and imply that sources within this region can be reliably sensed at vast distances. However, this theory has not been supported with evidence of direct observations from sound within this channel. Here we show that an infrasound sensor on a solar hot air balloon recorded the first infrasound detection of a ground truth airborne source while within this acoustic channel, which we refer to as the AtmoSOFAR channel. Our results support the existence of the AtmoSOFAR channel, demonstrate that acoustic signals can be recorded within it, and provide insight into the characteristics of recorded signals. Results also show a lack of detections on ground-based stations, highlighting the advantages of using balloon-borne infrasound sensors to detect impulsive sources at altitude.
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.
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.
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.
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.
As large systems of Li-ion batteries are being increasingly deployed, the safety of such systems must be assessed. Due to the high cost of testing large systems, it is important to extract key safety information from any available experiments. Developing validated predictive models that can be exercised at larger scales offers an opportunity to augment experimental data In this work, experiments were conducted on packs of three Li-ion pouch cells with different heating rates and states of charge (SOC) to assess the propagation behavior of a module undergoing thermal runaway. The variable heating rates represent slow or fast heating that a module may experience in a system. As the SOC decreases, propagation slows down and eventually becomes mitigated. It was found that the SOC boundary between propagation and mitigation was higher at a heating rate of 50 °C/min than at 10 °C/min for these cells. However, due to increased pre-heating at the lower heating rate, the propagation speed increased. Simulations were conducted with a new intra-particle diffusion-limited reaction model for a range of anode particle sizes. Propagation speeds and onset times were generally well predicted, and the variability in the propagation/mitigation boundary highlighted the need for greater uncertainty quantification of the predictions.
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.
Accurate targeting of radioisotope classifiers and estimators requires an understanding of the target problem space. In order to facilitate clear communication on expected model behavior and performance between practitioners and stakeholders on their problems, this questionnaire was created. Stakeholder responses form the basis of a trained model as well as the start of usage requirements for the model as it is integrated with analysis processes or detection systems. This questionnaire may also be useful to machine learning practitioners and gamma spectroscopists developing new algorithms as a starting point for characterizing their problem space, especially if they are using PyRIID.
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.
Multifidelity uncertainty quantification (MF UQ) sampling approaches have been shown to significantly reduce the variance of statistical estimators while preserving the bias of the highest-fidelity model, provided that the low-fidelity models are well correlated. However, maintaining a high level of correlation can be challenging, especially when models depend on different input uncertain parameters, which drastically reduces the correlation. Existing MF UQ approaches do not adequately address this issue. In this work, we propose a new sampling strategy that exploits a shared space to improve the correlation among models with dissimilar parameterization. We achieve this by transforming the original coordinates onto an auxiliary manifold using the adaptive basis (AB) method (Tipireddy and Ghanem, 2014). The AB method has two main benefits: (1) it provides an effective tool to identify the low-dimensional manifold on which each model can be represented, and (2) it enables easy transformation of polynomial chaos representations from high- to low-dimensional spaces. This latter feature is used to identify a shared manifold among models without requiring additional evaluations. We present two algorithmic flavors of the new estimator to cover different analysis scenarios, including those with legacy and non-legacy high-fidelity (HF) data. We provide numerical results for analytical examples, a direct field acoustic test, and a finite element model of a nuclear fuel assembly. For all examples, we compare the proposed strategy against both single-fidelity and MF estimators based on the original model parameterization.
As Machine Learning (ML) continues to advance, it is being integrated into more systems. Often, the ML component represents a significant portion of the system that reduces the burden on the end user or significantly improves task performance. However, the ML component represents an unknown complex phenomenon that is learned from collected data without the need to be explicitly programmed. Despite the improvement in task performance, the models are often black boxes. Evaluating the credibility and the vulnerabilities of ML models poses a gap in current test and evaluation practice. For high consequence applications, the lack of testing and evaluation procedures represents a significant source of uncertainty and risk. To help reduce that risk, here we present considerations to evaluate systems embedded with an ML component within a red-teaming inspired methodology. We focus on (1) cyber vulnerabilities to an ML model, (2) evaluating performance gaps, and (3) adversarial ML vulnerabilities.
The benefits of high-performance unidirectional carbon fiber composites are limited in many cost-driven industries due to the high cost relative to alternative reinforcement fibers. Low-cost carbon fibers have been previously proposed, but the longitudinal compressive strength continues to be a limiting factor or studies are based on simplifications that warrant further analysis. A micromechanical model is used to (1) determine if the longitudinal compressive strength of composites can be improved with noncircular carbon fiber shapes and (2) characterize why some shapes are stronger than others in compression. In comparison to circular fibers, the results suggest that the strength can be increased by 10%–13% by using a specific six-lobe fiber shape and by 6%–9% for a three-lobe fiber shape. A slight increase is predicted in the compressive strength of the study two-lobe fiber but has the highest uncertainty and sensitivity to fiber orientation and misalignment direction. The underlying mechanism governing the compressive failure of the composites was linked to the unique stress fields created by the lobes, particularly the pressure stress in the matrix. This work provides mechanics-based evidence of strength improvements from noncircular fiber shapes and insight on how matrix yielding is altered with alternative fiber shapes.
The development of additively-manufactured (AM) 316L stainless steel (SS) using laser powder bed fusion (LPBF) has enabled near net shape components from a corrosion-resistant structural material. In this article, we present a multiscale study on the effects of processing parameters on the corrosion behavior of as-printed surfaces of AM 316L SS formed via LPBF. Laser power and scan speed of the LPBF process were varied across the instrument range known to produce parts with >99 % density, and the macroscale corrosion trends were interpreted via microscale and nanoscale measurements of porosity, roughness, microstructure, and chemistry. Porosity and roughness data showed that porosity φ decreased as volumetric energy density Ev increased due to a shift in the pore formation mechanism and that roughness Sq was due to melt track morphology and partially fused powder features. Cross-sectional and plan-view maps of chemistry and work function ϕs revealed an amorphous Mn-silicate phase enriched with Cr and Al that varied in both thickness and density depending on Ev. Finally, the macroscale potentiodynamic polarization experiments under full immersion in quiescent 0.6 M NaCl showed significant differences in breakdown potential Eb and metastable pitting. In general, samples with smaller φ and Sq values and larger ϕs values and homogeneity in the Mn-silicate exhibited larger Eb. The porosity and roughness effects stemmed from an increase to the overall number of initiation sites for pitting, and the oxide phase contributed to passive film breakdown by acting as a crevice former or creating a galvanic couple with the SS.
Characterizing interface trap states in commercial wide bandgap devices using frequency-based measurements requires unconventionally high probing frequencies to account for both fast and slow traps associated with wide bandgap materials. The C − ψ S technique has been suggested as a viable quasi-static method for determining the interface trap state densities in wide bandgap systems, but the results are shown to be susceptible to errors in the analysis procedure. This work explores the primary sources of errors present in the C − ψ S technique using an analytical model that describes the apparent response for wide bandgap MOS capacitor devices. Measurement noise is shown to greatly impact the linear fitting routine of the 1 / C S ∗ 2 vs ψ S plot to calibrate the additive constant in the surface potential/gate voltage relationship, and an inexact knowledge of the oxide capacitance is also shown to impede interface trap state analysis near the band edge. In addition, a slight nonlinearity that is typically present throughout the 1 / C S ∗ 2 vs ψ S plot hinders the accurate estimation of interface trap densities, which is demonstrated for a fabricated n-SiC MOS capacitor device. Methods are suggested to improve quasi-static analysis, including a novel method to determine an approximate integration constant without relying on a linear fitting routine.