Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model based development”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,369 records · Page 76

Dust Aerosol Impact on North Africa Climate: A GCM Investigation of Aerosol-Cloud-Radiation Interactions Using A-Train Satellite Data

The climatic effects of dust aerosols in North Africa have been investigated using the atmospheric general circulation model (AGCM) developed at the University of California, Los Angeles (UCLA). The model includes an efficient and physically based radiation parameterization scheme developed specifically for application to clouds and aerosols. Parameterization of the effective ice particle size in association with the aerosol indirect effect based on cloud and aerosol data retrieved from A-Train satellite observations have been employed in the climate model simulations. Offline simulations reveal that the direct solar, IR, and net forcings by dust aerosols generally increase with increasing aerosol optical depth (AOD). When the dust semi-direct effect is included with the presence of ice clouds, positive IR radiative forcing is enhanced, since ice clouds trap substantial IR radiation, while the positive solar forcing with dust aerosols alone has been changed to negative values due to the strong reflection of solar radiation by clouds, indicating that cloud forcing could exceed aerosol forcing. With the aerosol indirect effect, the net cloud forcing is generally reduced for ice water path (IWP) larger than 20 g m-2. The magnitude of the reduction increases with IWP. AGCM simulations show that the reduced ice crystal mean effective size due to the aerosol first indirect effect result in less OLR and net solar flux at the top of the atmosphere over the cloudy area of the North Africa region because ice clouds with smaller size trap more IR radiation and reflect more solar radiation. The precipitation in the same area, however, increases due to the aerosol indirect effect on ice clouds, corresponding to the enhanced convection as indicated by reduced OLR. The increased precipitation seems to be associated with enhanced ice water contents in this region. The 200 mb radiative heating rate shows more cooling with the aerosol indirect effect since greater cooling is produced at the cloud top with smaller ice crystal size. The 500 mb omega indicates strong upward motion, which, together with the increased cooling effect, results in the increased ice water contents. Adding the aerosol direct effect into the model simulation reduces the precipitation in the normal rainfall band over North Africa, where precipitation is shifted to the south and the northeast produced by the absorption of sunlight and the subsequent heating of the air column by dust particles. As a result, rainfall is drawn further inland to the northeast. This study represents the first attempt to quantify the climate impact of aerosol indirect effect using a GCM in connection with A-train satellite data. The parameterization for the aerosol first indirect effect developed in this study can be readily incorporated for application to any other GCMs.

54 ENVIRONMENTAL SCIENCES↗

Agro-waste shaped porous Al2O3/Ni composites: Corrosion resistance performance and artificial neural network modelling

Highlights: • Porous Al{sub 2}O{sub 3}/Ni composites were developed with agro-waste pore-forming agents. • The RH and SCB graded composites showed Ni{sub 3}Al{sub 2}SiO{sub 8} and NiAl{sub 2}O{sub 4} spinels respectively. • Superior chemical stability was demonstrated by the Al{sub 2}O{sub 3}-2Ni-RH composite. • The ANN model demonstrated a good performance with R{sup 2} value of approx. 0.99. - Abstract: In the present study, an analysis on the combined effect of nickel (Ni) reinforcement and pore former type in characterizing the corrosion behavior of composite porous alumina ceramics was performed. In order to showcase the potential of the new porous ceramics, pore-forming agents (PFAs) from rice husk (RH) and sugarcane bagasse (SCB) were used in shaping the plain and composite porous alumina samples having sample formulation of Al{sub 2}O{sub 3}-xNi-PFA; x = 0, 2, 4, 6 and 8 wt%. Results showed that the emergence of a highly stable Ni{sub 3}Al{sub 2}SiO{sub 8} spinelloid phase in the RH-graded composites enhanced their chemical stability in the corrosive mediums (10 wt% NaOH and 20 wt% H{sub 2}SO{sub 4}) relative to the plain and the corresponding SCB-graded counterparts. An artificial neural network (ANN) model has been developed for predicting the corrosion behavior of the plain and composite porous alumina ceramics based on the experimental data. The developed ANN model satisfactorily predicted the percent mass losses of the porous ceramics in strong alkali and strong acid solutions with coefficient of determination (R{sup 2}) of approximately 0.99.

36 MATERIALS SCIENCE↗

Physics-informed machine learning of the Lagrangian dynamics of velocity gradient tensor

Reduced models describing the Lagrangian dynamics of the velocity gradient tensor (VGT) in homogeneous isotropic turbulence (HIT) are developed under the physics-informed machine learning (PIML) framework. We consider the VGT at both Kolmogorov scale and coarse-grained scale within the inertial range of HIT. Building reduced models requires resolving the pressure Hessian and subfilter contributions, which is accomplished by constructing them using the integrity bases and invariants of the VGT. The developed models can be expressed using the extended tensor basis neural network (TBNN) introduced by Ling et al. [J. Fluid Mech. 807, 155 (2016)]. Physical constraints, such as Galilean invariance, rotational invariance, and incompressibility condition, are thus embedded in the models explicitly. Our PIML models are trained on the Lagrangian data from a high-Reynolds number direct numerical simulation (DNS). To validate the results, we perform a comprehensive out-of-sample test. We observe that the PIML model provides an improved representation for the magnitude and orientation of the small-scale pressure Hessian contributions. Statistics of the flow, as indicated by the joint PDF of second and third invariants of the VGT, show good agreement with the “ground-truth” DNS data. A number of other important features describing the structure of HIT are reproduced by the model successfully. We have also identified challenges in modeling inertial range dynamics, which indicates that a richer modeling strategy is required. This helps us identify important directions for future research, in particular towards including inertial range geometry into the TBNN.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Plasma flow healing of magnetic islands in stellarators

Recent experiments from the large helical device (LHD) demonstrate a correlation between the 'healing' of vacuum magnetic islands in stellarators and changes in the plasma flow. A model explaining this phenomenon is developed based on self-consistent torque balance and island evolution equations. In conventional stellarators, neoclassical flow damping physics plays an important role in establishing the flow profiles. The balance of neoclassical damping and cross-field viscosity produces a radial boundary layer for the plasma rotation profile outside the separatrix of a locked magnetic island. The width of this boundary layer decreases as the plasma becomes less collisional. Associated with these flow effects are plasma currents flowing in the island region that attempt to suppress island formation. These currents are enhanced as the collisionality drops making magnetic island healing occur more readily in high temperature conventional stellarators. The analytic theory produces a critical {beta} for healing that scales monotonically with collisionality and is in qualitative agreement with LHD observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Human in vitro metabolism of an environmental mixture of polycyclic aromatic hydrocarbons (PAH) found at the Portland Harbor Superfund Site

Polycyclic aromatic hydrocarbons (PAHs) are widespread environmental contaminants that pose health risks to humans. Toxicity testing approaches of PAHs have evolved from traditional rodent models to New Approach Methodologies (NAMs), such as high-throughput screening in zebrafish, enabling rapid evaluation of chemical hazards. However, translating toxicity findings from laboratory systems to humans remains difficult due to complexity and species-specific differences. Chemical dosimetry modeling offers a quantitative framework to bridge this gap, but its accuracy depends on robust knowledge of PAH metabolism. The objective of this study was to measure human metabolism rates of Supermix-10, the ten most abundant PAHs found at the Portland Harbor Superfund Site, to support development of human pharmacokinetic models. We incubated individual PAHs from Supermix-10 in pooled human liver microsomes and quantified parent PAH disappearance using high-performance liquid chromatography (HPLC) with UV and florescent detection. To assess the potential of mixture interactions, we also measured metabolism of all 10 compounds in an equimolar mixture and compared rates of parent disappearance to those observed for individual PAHs. All Supermix-10 PAHs demonstrated rapid parent compound disappearance in human hepatic microsomes. PAHs grouped into three metabolism patterns: high metabolism rates and capacity (2-methylnaphthalene, acenaphthylene, fluorene, naphthalene), high affinity metabolism that rapidly achieves low-level saturation (benzo[a]anthracene, chrysene), and moderate metabolism rates and capacity (fluoranthene, pyrene, retene, phenanthrene). Smaller PAHs exhibited faster metabolism, and higher metabolism rates correlated inversely with molecular weight. When incubated in an equimolar mixture, Supermix-10 demonstrated significantly slower metabolism (47–89 %) compared to metabolism of individual PAHs at the same concentration. These findings enhance our understanding of PAH metabolism in humans and demonstrate significant mixture interactions under the conditions tested. Furthermore, our findings offer insights into the metabolic behavior of Supermix-10 and provide critical metabolism rate data to support the development of physiological based pharmacokinetic (PBPK) models. Dosimetry models can translate PAH chemical dosimetry from high-throughput testing platforms, like zebrafish and cellular system assays, to human exposures enhancing the accuracy and reliability of PAH risk assessments.

2-methylnaphthalene↗

Calibration verification for stochastic agent-based disease spread models

Accurate disease spread modeling is crucial for identifying the severity of outbreaks and planning effective mitigation efforts. To be reliable when applied to new outbreaks, model calibration techniques must be robust. However, current methods frequently forgo calibration verification (a stand-alone process evaluating the calibration procedure) and instead use overall model validation (a process comparing calibrated model results to data) to check calibration processes, which may conceal errors in calibration. In this work, we develop a stochastic agent-based disease spread model to act as a testing environment as we test two calibration methods using simulation-based calibration, which is a synthetic data calibration verification method. The first calibration method is a Bayesian inference approach using an empirically-constructed likelihood and Markov chain Monte Carlo (MCMC) sampling, while the second method is a likelihood-free approach using approximate Bayesian computation (ABC). Simulation-based calibration suggests that there are challenges with the empirical likelihood calculation used in the first calibration method in this context. These issues are alleviated in the ABC approach. Despite these challenges, we note that the first calibration method performs well in a synthetic data model validation test similar to those common in disease spread modeling literature. We conclude that stand-alone calibration verification using synthetic data may benefit epidemiological researchers in identifying model calibration challenges that may be difficult to identify with other commonly used model validation techniques.

60 APPLIED LIFE SCIENCES↗

Technical Report on Waveform Fit Metrics for Global Models

The new WAVEFORMS Initiative in the Ground-based Nuclear Detonation Detection (GNDD) program includes an increased emphasis on the development of Earth models and methods to predict entire seismic and acoustic waveforms more accurately. In general, this increased emphasis is predicated on the need to better characterize seismic events and provide improved model-based discrimination between event types including earthquakes and explosions. More specifically, while current moment tensor inversion methods tend to work well for larger events (M>~4) using tuned 1-D Earth models, the development of state-of-the-art 3-D models and methods is required for the prediction of shorter period waves over large areas for discrimination of smaller events. There is no standard metric for model-based waveform prediction accuracy used in the waveform modeling/inversion community. However, there are several popular waveform misfit definitions; and minimizing the corresponding objective functions is the goal of waveform inversion. Some example misfit definitions employed for adjoint waveform tomography include measures of simple travel time differences (e.g. Tape et al., 2010), cross-correlation travel time differences (e.g. Luo and Schuster, 1991), multi-taper frequency dependent methods (e.g. Lei et al., 2020), time-frequency phase misfit functions (e.g. Fichtner 2010; Rodgers et al., 2022), normalized cross-correlation methods (e.g. Tao et al., 2018), and others. In some cases, these misfit definitions also involve complicated weighting schemes and summations over multiple frequency bands making it difficult to duplicate the misfit measurement with alternative models and datasets. Although each of the misfit definitions mentioned above are useful for developing waveform models, the actual misfit values are not usually meaningful outside of a given project, model, and/or dataset. Therefore, it is difficult to understand and communicate model performance for predicting waveforms and comparing to other models and/or new model iterations with a different dataset. Therefore, there is a need for a generalized method for evaluating overall model performance that is independent from the specific misfit chosen to develop the waveform models that is also intuitive and meaningful. In this report, we describe a new metric we refer to as ‘Percent of Correlated Signal’. The following sections describe and demonstrate the metric with a case study event and a more rigorous test using a random selection of globally distributed events. While the focus here is on global tomography models, the metric is meant to applicable to regional ‘wiggle-for-wiggle’ waveform models/studies as well.

58 GEOSCIENCES↗

Produced Water and Waste Heat-aided Blowdown Water Treatment: Using Chemical and Energy Synergisms for Value Creation

The project objective was to develop a cooling blowdown water (BDW) treatment process utilizing produced water (PW) and low-grade heat to maximize water reuse and saleable by-product generation while reducing chemical and energy footprints of the treatment. The proposed treatment process consists of mixing, softening, organics and suspended solids removal, reverse osmosis (RO), thermal desalination, and brine electrolysis. BDW samples collected from a local coal-fired power plant and PW samples from two shale gas production wells were used in this study. Each treatment unit was first designed and tested to quantify its treatment efficiency, and its chemical and energy requirements. In addition, a process model was developed and model simulations were conducted based on the experimental results and literature data to optimize the treatment process. A techno-economic analysis was conducted to quantify chemical and energy savings as well as production of 10-lb brine as a saleable product. With the field-collected BDW and PW samples, mixing experiments determined a volumetric mixing ratio 10:1 (BDW:PW) resulted in the best performance of multivalent ions removal and largest chemical savings for softening. Softening of the BDW/PW mixtures using alkaline chemicals (Na 2 CO 3 and NaOH) achieved 95%-100% removal of scaling-forming cations (Ca, Mg, Fe, Ba, Sr) and 60% of silicon, and 10% of total organic carbon (TOC). The mixing and softening treatments yielded an effluent with total dissolved solids (TDS) concentration of 23 g/L. Activated carbon (AC) filtration removed TOC to a low level (< 3 mg/L) and further removed remaining scale-forming divalent metals and silica from the softened water. The AC filtration resulted in a slight reduction of TDS from 23 g/L to 20 g/L, leaving behind only mostly monovalent ions (i.e., sodium and chloride) in the filtered water. These pretreatments yielded a feed water that met the criteria of the downstream reverse osmosis (RO) to prevent membrane fouling. A cross-flow RO system was used to further concentrate the TDS of the AC effluent. Various factors including TDS, pH, and applied pressure were examined and optimal conditions were determined for the co-treatment process. An integrated process consisting of mixing, softening, AC filtration and RO was used to treat a continuous flow (0.25 – 1.2 L/min, or 0.07 – 0.32 gpm) and successfully generated RO permeate as product water (TDS < 0.5 g/L) for reuse in cooling operation, and a concentrate (TDS ~ 45 g/L) to be further treated in a thermal desalination unit. These flow rates meet the FOA’s criterion of 0.01 – 1 gpm. Overall, the co-treatment of BDW/PW allowed shorter ramp-up time compared to treatment of BDW alone. It resulted in 40% and 55% savings of Na 2 CO 3(s) and NaOH, respectively, compared to treating the BDW and PW individually for the same level of softening. The co-treatment also resulted in a 29% energy saving compared to treatment of BDW only for the level of TDS concentration. A thermal desalination system was designed using CFD simulations and manufactured in the WVU Innovation Hub for further treatment of the RO concentrate to generate 10-lb brine. The system has a design flow rate of 2 gpm and has been successfully tested. A bench-scale brine electrolysis system was developed for on-site generation of chlorine/hypochlorite (Cl 2 /OCl - ) and caustic soda (NaOH) as useful chemicals for the co-treatment process. Using salt solutions (0.5 M and 1 M), the system achieved faradaic efficiencies of 93%-97% and 70%-77% for caustic soda and chlorine/hypochlorite generation, respectively. An economic analysis showed that the electricity costs for on-site generation of these chemicals were significantly lower than the chemical prices offered by suppliers. An industrial-scale process model consisting of mixing, softening, AC filtration, RO, thermal desalination, and brine electrolysis was developed using the Aspen Plus V9 in conjunction with Aspen Custom Modeler V9. The model serves as a solvable Aspen Plus model and as basis to form the costing infrastructure. In addition, techno-economic analysis considering capital, operating, and transportation costs was conducted. An optimization solution showed that produced water for mixing is still advantageous in low quantities. The optimum solution approaches a leveled cost of water (LCW) of 2 $/m 3 which becomes cost competitive with nominal water treatment prices.

20 FOSSIL-FUELED POWER PLANTS↗

Utilizing Adjoint-Based Error Estimates for Surrogate Models to Accurately Predict Probabilities of Events

In thist study, we develop a procedure to utilize error estimates for samples of a surrogate model to compute robust upper and lower bounds on estimates of probabilities of events. We show that these error estimates can also be used in an adaptive algorithm to simultaneously reduce the computational cost and increase the accuracy in estimating probabilities of events using computationally expensive high-fidelity models. Specifically, we introduce the notion of reliability of a sample of a surrogate model, and we prove that utilizing the surrogate model for the reliable samples and the high-fidelity model for the unreliable samples gives precisely the same estimate of the probability of the output event as would be obtained by evaluation of the original model for each sample. The adaptive algorithm uses the additional evaluations of the high-fidelity model for the unreliable samples to locally improve the surrogate model near the limit state, which significantly reduces the number of high-fidelity model evaluations as the limit state is resolved. Numerical results based on a recently developed adjoint-based approach for estimating the error in samples of a surrogate are provided to demonstrate (1) the robustness of the bounds on the probability of an event, and (2) that the adaptive enhancement algorithm provides a more accurate estimate of the probability of the QoI event than standard response surface approximation methods at a lower computational cost.

97 MATHEMATICS AND COMPUTING↗

NWChem: scalable parallel computational chemistry

NWChem is a general purpose computational chemistry code specifically designed to run on distributed memory parallel computers. The core functionality of the code focuses on molecular dynamics, Hartree-Fock and density functional theory methods for both plane-wave basis sets as well as Gaussian basis sets, tensor contraction engine based coupled cluster capabilities and combined quantum mechanics/molecular mechanics descriptions. It was realized from the beginning that scalable implementations of these methods required a programming paradigm inherently different from what message passing approaches could offer. In response a global address space library, the Global Array Toolkit, was developed. The programming model it offers is based on using predominantly one-sided communication. This model underpins most of the functionality in NWChem and the power of it is exemplified by the fact that the code scales to tens of thousands of processors. In this paper the core capabilities of NWChem are described as well as their implementation to achieve an efficient computational chemistry code with high parallel scalability. NWChem is a modern, open source, computational chemistry code1 specifically designed for large scale parallel applications2. To meet the challenges of developing efficient, scalable and portable programs of this nature a particular code design was adopted. This code design involved two main features. First of all, the code is build up in a modular fashion so that a large variety of functionality can be integrated easily. Secondly, to facilitate writing complex parallel algorithms the Global Array toolkit was developed. This toolkit allows one to write parallel applications in a shared memory like approach, but offers additional mechanisms to exploit data locality to lower communication overheads. This framework has proven to be very successful in computational chemistry but is applicable to any engineering domain. Within the context created by the features above NWChem has grown into a general purpose computational chemistry code that supports a wide variety of energy expressions and capabilities to calculate properties based there upon. The main energy expressions are classical mechanics force fields, Hartree-Fock and DFT both for finite systems and condensed phase systems, coupled cluster, as well as QM/MM. For most energy expressions single point calculations, geometry optimizations, excited states, and other properties are available. Below we briefly discuss each of the main energy expressions and the critical points involved in scalable implementations thereof.

99 GENERAL AND MISCELLANEOUS↗

Digital Architecture Planning Model

As part of the U.S. Department of Energy’s Light Water Reactor Sustainability Program, the Digital Architecture (DA) Project focuses on providing a model that nuclear utilities can refer to when planning deployment of advanced technologies. The digital architecture planning model (DAPM) is the methodology for mapping power plant operational and support activities into a DA that unifies all data sources needed by the utilities to operate their plants. The DA is defined as a collection of information technology capabilities needed to support and integrate a wide spectrum of real-time digital capabilities for performance improvements of nuclear power plants. DA can be thought of as integration of the separate instrumentation and control and information systems already in place in nuclear power plants, which are brought together for the purpose of creating new levels of automation in plant work activities. A major objective in DAPM development was to survey all key areas that needed to be reviewed in order for a utility to make knowledgeable decisions regarding needs and plans to implement a DA at the plant. The development was done in two steps. First, researchers surveyed the nuclear industry in order to learn their near-term plans for adopting new advanced capabilities and implementing a network (i.e., wireless and wire) infrastructure throughout the plant, including the power block. Secondly, a literature review covering regulatory documents, industry standards, and technical research reports and articles was conducted. The objective of the review was to identify key areas to be covered by the DAPM, which included the following: 1. The need for a DA and its benefits to the plant 2. Resources required to implement the DA 3. Challenges that need to be addressed and resolved to implement the DA 4. Roles and responsibilities of the DA implementation plan. The DAPM was developed based on results from the survey and the literature review. Model development, including the survey results and conclusions made about the key areas during the literature review, are described in this report.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Development of Versatile Compressor Modeling using Approximation Techniques for Alternative Refrigerants Evaluation

Refrigerants are the life-blood of vapor compression systems that are widely used in Heating, Ventilation, Air-Conditioning, and Refrigeration (HVAC&R) applications. The HVAC&R community is currently transitioning from main-stream refrigerants that have high Global Warming Potential (GWP) to alternative lower-GWP refrigerants. During this transition, it is important to account for the life cycle climate performance of alternative refrigerants since their performance will be different than that of higher-GWP refrigerants. This requires the evaluation of the system performance with the new refrigerants. Unfortunately, it is extremely difficult to predict the realistic performance of new alternative refrigerants without experimental validation. One of the main challenges in this regard is modeling the compressor performance with high fidelity due to the complex interaction of operating parameters, geometry, boundary conditions, and fluid properties. High fidelity compressor models are computationally expensive and require significant pre-processing to evaluate the performance of alternative refrigerants. This paper presents a new approach to modeling compressor performance when alternative refrigerants are used. The new modeling concept relies on using existing compressor performance to create an approximate model that captures the dependence of compressor performance on key operating parameters and fluid properties. The model can be built using a myriad of approximation techniques. This paper focuses on Kriging-based techniques to develop higher fidelity approximate compressor models. Baseline and at least one alternative refrigerant performance data are used to build the model. The model accuracy was evaluated by comparing the model results with compressor performance data using other refrigerants. Preliminary results show that the approximate model can predict the compressor mass flow rate and power consumption within 5%.

Approximation↗

Development of an Energy Services Interface for the EGoT

The Energy Services Interface (ESI) is a set of rules that ensure private, secure, and trustworthy information exchange between Grid Service Providers (GSPs) and utility customers. Large-scale adoption of Distributed Energy Resources (DERs) will be necessary for GSPs to dispatch effective grid services, and to stimulate technological innovations in DER and DERMS (DER Management Systems) technologies, as well as novel grid service programs. The ESI promotes these objectives by advancing a set of rules and interoperability requirements that define bi-directional, service-oriented, logical interfaces between GSPs and customers’ DERs, with expectations for privacy, security, and trust. The ESI rules and interoperability requirements establish boundaries between customers and GSPs that delineate the functions and responsibilities that must be implemented by the developers of Energy Grid of Things (EGoT) ecosystem products. These rules impose constraints on the implementation of a DERMS. The ESI interoperability requirements are based on the Interoperability Maturity Model (IMM), developed by the Grid Modernization Laboratory Consortium. By emphasizing private, secure, and trustworthy information exchange, and by mandating a service-oriented and interoperable architecture, the ESI promotes the development of an EGoT ecosystem that motivates customer participation and technological innovation. Interoperability will encourage innovation by reducing barriers to entry and increasing confidence of stakeholders. Customers will be willing to participate in DER service programs that establish trust and emphasize customer choice. Large-scale customer participation ensures GSPs have ample DER resources to provide grid services that have significant impact on grid reliability and reduce electricity cost for consumers. This in turn signals economic opportunities that encourage innovation, resulting in the development of a robust EGoT ecosystem.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Crack formation in strained β-(AlxGa 1−x ) 2 O 3 films grown on (010) β-Ga 2 O 3 substrates

The cracking and local strain relaxation in (010) (Al x Ga 1−x ) 2 O 3 films grown on Ga 2 O 3 substrates are assessed in terms of film composition and thickness. We utilize x-ray diffraction and electron microscopy techniques combined with simulation and modeling to investigate that film cracking on curvature (flatness) has a directly proportional relationship with film thickness and/or aluminum content. Cross section transmission electron microscopy reveals that cracks along both the (001) and (100) cleavage planes penetrate into the substrate. The diffuse scattered intensity observed in reciprocal space maps (RSMs) is directly correlated with the tilt that is introduced due to the change in the deformation conditions near the cracks. While asymmetric RSMs show that the layers are fully strained, the diffuse scattering distribution in reciprocal space can be interpreted to show that the cracking relaxes and locally tilts the lattice ∼350 nm from the crack edges, which is consistent with the larger radius of curvature associated with the films with higher crack densities. For example, a 200 nm (Al 0.13 Ga 0.87 ) 2 O 3 thick film has an average inter-crack spacing of 3.3 µm, so most of the epitaxial layer is fully strained except near the cracks where it deforms elastically and is consistent with the gallium oxide Poisson ratio. A reciprocal space model was developed, which imports the strain and tilt distributions (based on finite element modeling) to match the features observed in the experimental maps. We also note that previous studies involving (Al x Ga 1−x ) 2 O 3 films may show evidence of cracking as observed in their symmetric and asymmetric RSMs.

36 MATERIALS SCIENCE↗

Prediction of Radial Vibration in Switched Reluctance Machines

Origins of vibration in switched reluctance machines (SRMs) are investigated. Accordingly, an input-output model based on the mechanical impulse response of the SRMis developed. The proposed model is derived using an experimental approach. Using the proposed approach, vibration of the stator frame is captured and experimentally verified.

Lin, CJ↗

Modeling and Analysis for Spent Nuclear Fuel Seismic Testing

Note: This is the final public release version of PNNL-31671 Draft, which was previously released to the sponsor for review. The intent is the public release version will be PNNL-31671. There are no significant changes from the earlier draft version. This unlimited distribution report is the deliverable for M3SF-21PN010202014: Modeling and Analysis for Spent Nuclear Fuel Seismic Testing. This report summarizes the modeling, analysis, and test plan support completed by Pacific Northwest National Laboratory for the spent nuclear fuel dry storage system seismic test plan through May of 2021. Test plan preparation is planned to continue until the seismic test is completed in July of 2022. This report covers preliminary structural dynamic model development and computer-aided design of test hardware. Based on preliminary modeling, the strongest earthquakes under consideration in this test program provide mechanical loading on the spent nuclear fuel that is comparable to the 30 cm cask drop scenario, although the loads do not appear to be strong enough to cause significant permanent deformation of the fuel assembly spacer grids. The potential for grid deformation during the test will be assessed when the proposed shake table motion becomes available. The weakest earthquakes considered in this test are expected to be comparable to the mechanical loads witnessed in the multimodal transportation test of 2017. The horizontal canister system is predicted to provide nearly-uniform loading condition on the fuel assemblies it contains, while the vertical cask system is predicted to cause non-uniform dynamic loads on the fuel assemblies. In the horizontal case, gravity keeps the fuel assemblies in contact with one basket wall surface unless the loads are strong enough to cause a separation. In the vertical case, the fuel assemblies are relatively long and slender, and seismic motion in the anticipated test range is predicted to cause the assemblies to lean, tilt, and impact the basket walls throughout the seismic event. These gap closures are a nonlinear force transmission condition, so the vertical system is expected to have more variation and variability than the horizontal system. While the horizontal system is expected to have a relatively more linear response than the vertical system, there is still the potential for nonlinear behavior in the horizontal system because the fuel assemblies are free to slide, bounce, and impact the basket walls if the seismic loads are strong enough. One major conclusion of this study is that the use of mixed fuel assemblies in the canister will be acceptable. There are differences in overall system response when the mass, center of gravity, or gaps are changed, but the changes in response are within the bounds of a system that contains completely homogeneous fuel assemblies. One important observation is that the loading for each individual fuel assembly within the vertical canister is expected to be different from that of the others because of nonlinearities in the system. The horizontal canister system is expected to have a more uniform and more predictable response than the vertical canister, making variations in fuel assembly characteristics easier to account for. One important recommendation that comes from this modeling work is to repeat some of the strongest shake tests at different angles, particularly for the vertical cask system. Modeling predicts that strong seismic motion will cause the fuel assemblies to close initial gaps and impact the fuel basket walls when the canister is in the vertical cask configuration. Each shake test will impose a pre-defined three-dimensional time history on the cask system.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Reduction of a collisional-radiative mechanism for argon plasma based on principal component analysis

This article considers the development of reduced chemistry models for argon plasmas using Principal Component Analysis (PCA) based methods. Starting from an electronic specific Collisional-Radiative model, a reduction of the variable set (i.e., mass fractions and temperatures) is proposed by projecting the full set on a reduced basis made up of its principal components. Thus, the flow governing equations are only solved for the principal components. The proposed approach originates from the combustion community, where Manifold Generated Principal Component Analysis (MG-PCA) has been developed as a successful reduction technique. Applications consider ionizing shock waves in argon. The results obtained show that the use of the MG-PCA technique enables for a substantial reduction of the computational time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Bayesian modeling of inconsistent plastic response due to material variability

The advent of fabrication techniques such as additive manufacturing has focused attention on the considerable variability of material response due to defects and other microstructural aspects. This variability motivates the development of an enhanced design methodology that incorporates inherent material variability to provide robust predictions of performance. In this work, we develop plasticity models capable of representing the distribution of mechanical responses observed in experiments using traditional plasticity models of the mean response and recently developed uncertainty quantification (UQ) techniques. To account for material response variability through variations in physical parameters, we adapt a recent Bayesian embedded modeling error calibration technique. We use Bayesian model selection to determine the most plausible of a variety of plasticity models and the optimal embedding of parameter variability. To expedite model selection, we develop an adaptive importance-sampling-based numerical integration scheme to compute the Bayesian model evidence. In conclusion, we demonstrate that the new framework provides predictive realizations that are superior to more traditional ones, and how these UQ techniques can be used in model selection and assessing the quality of calibrated physical parameters.

36 MATERIALS SCIENCE↗