Search NASA⌕ Search

SEARCH · Search NASA

Results for “Coefficient of variation”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Universal orbital and magnetic structures in infinite-layer nickelates

We conducted a comparative study of the rare-earth infinite-layer nickelate films, RNiO 2 (R = La, Pr, and Nd), using resonant inelastic x-ray scattering (RIXS). We found that the gross features of the orbital configurations are essentially the same, with minor variations in the detailed hybridization. For low-energy excitations, we unambiguously confirm the presence of damped magnetic excitations in all three compounds. By fitting to a linear spin-wave theory, comparable spin exchange coupling strengths and damping coefficients are extracted, indicating a universal magnetic structure in the infinite-layer nickelates. Interestingly, while signatures of a charge order are observed in LaNiO 2 in the quasielastic region of the RIXS spectrum, it is absent in NdNiO 2 and PrNiO 2 . Furthermore, this prompts further investigation into the universality and the origins of charge order within the infinite-layer nickelates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Seasonal Variation of Dust Aerosol Vertical Distribution in Arctic Based on Polarized Micropulse Lidar Measurement

This study investigates the seasonal variation of dust aerosol vertical distribution using polarized Micropulse lidar (MPL) measurements at the Atmospheric Radiation Measurement (ARM) North Slope of Alaska (NSA) observatory from January 2013 to September 2017. For the first time, multi-year aerosol backscatter coefficients are retrieved at the ARM NSA site from MPL measurements and are consistent with co-located high spectral resolution lidar (HSRL) measurements. The high-quality aerosol backscatter coefficient retrievals are used to derive the particle depolarization ratio (PDR) at the wavelength of 532 nm, which is used to identify the presence of dust aerosols. The annual cycles of the vertical distributions of dust backscatter coefficient and PDR and dust aerosol optical depth (DAOD) show that aerosol loading has a maximum in late winter and early spring but a minimum in late summer and early autumn. Vertically, dust aerosol occurs in the entire troposphere in spring and winter and in the low and middle troposphere in summer and autumn. Because dust aerosols are effective ice nuclei, the seasonality of dust aerosol vertical distribution has important implications for the Arctic climate through aerosol–cloud–radiation interactions, primarily through impacting mixed-phase cloud processes.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Modeling of streamflow in a 30 km long reach spanning 5 years using OpenFOAM 5.x

Abstract. Developing accurate and efficient modeling techniques for streamflow at the tens-of-kilometers spatial scale and multi-year temporal scale is critical for evaluating and predicting the impact of climate- and human-induced discharge variations on river hydrodynamics. However, achieving such a goal is challenging because of limited surveys of streambed hydraulic roughness, uncertain boundary condition specifications, and high computational costs. We demonstrate that accurate and efficient three-dimensional (3-D) hydrodynamic modeling of natural rivers at 30 km and 5-year scales is feasible using the following three techniques within OpenFOAM, an open-source computational fluid dynamics platform: (1) generating a distributed hydraulic roughness field for the streambed by integrating water-stage observation data, a rough wall theory, and a local roughness optimization and adjustment strategy; (2) prescribing the boundary condition for the inflow and outflow by integrating precomputed results of a one-dimensional (1-D) hydraulic model with the 3-D model; and (3) reducing computational time using multiple parallel runs constrained by 1-D inflow and outflow boundary conditions. Streamflow modeling for a 30 km long reach in the Columbia River (CR) over 58 months can be achieved in less than 6 d using 1.1 million CPU hours. The mean error between the modeled and the observed water stages for our simulated CR reach ranges from −16 to 9 cm (equivalent to approximately ±7 % relative to the average water depth) at seven locations during most of the years between 2011 and 2019. We can reproduce the velocity distribution measured by the acoustic Doppler current profiler (ADCP). The correlation coefficients of the depth-averaged velocity between the model and ADCP measurements are in the range between 0.71 and 0.83 at 75 % of the survey cross sections. With the validated model, we further show that the relative importance of dynamic pressure versus hydrostatic pressure varies with discharge variations and topography heterogeneity. Given the model's high accuracy and computational efficiency, the model framework provides a generic approach to evaluate and predict the impacts of climate- and human-induced discharge variations on river hydrodynamics at tens-of-kilometers and decadal scales.

58 GEOSCIENCES↗
DOE OSTI2024

Improving Low‐Cloud Fraction Prediction Through Machine Learning

Abstract In this study, we evaluated the performance of machine learning (ML) models (XGBoost) in predicting low‐cloud fraction (LCF), compared to two generations of the community atmospheric model (CAM5 and CAM6) and ERA5 reanalysis data, each having a different cloud scheme. ML models show a substantial enhancement in predicting LCF regarding root mean squared errors and correlation coefficients. The good performance is consistent across the full spectrums of atmospheric stability and large‐scale vertical velocity. Employing an explainable ML approach, we revealed the importance of including the amount of available moisture in ML models for representing spatiotemporal variations in LCF in the midlatitudes. Also, ML models demonstrated marked improvement in capturing the LCF variations during the stratocumulus‐to‐cumulus transition (SCT). This study suggests ML models' great potential to address the longstanding issues of “too few” low clouds and “too rapid” SCT in global climate models.

Geology↗
DOE OSTI2024

Liquid–Liquid Equilibrium Prediction in Fast Pyrolysis Bio-Oil Systems: A Framework for Incorporating Bio-Oil Complexity

The study of mixtures of bio-oil, water and organic solvents in different proportions can serve as a cost-effective analysis of its content due to the formation of immiscible phases. This manuscript attempts to replicate experimentally determined partition coefficients (K OW ) of relevant species present in fast pyrolysis bio-oil (FPBO). A commercial flowsheeting simulator with surrogate bio-oil model representation is used. Concurrently, pyrolytic lignins in FPBO (‘pyrolignin’) do not have an agreed-upon structural representation, and the literature is ripe with wide variations of said representations. Thus, during the description of FPBO, this pyrolignin fraction was modeled using 20 possible structures (phenolic dimers to tetramers), with the goal of determining the structures for which the experimental data are best described. Two cases were considered: Case 1 normalized the reported experimental mass balance, while Case 2 included the unreported fraction in the mass balance to the total pyroligin. Please, add here a comment on the prediction of the Water oil equilibrium. The best KOW predictions for levoglucosan (LVG) were obtained when the system was modeled with no pyrolignin, presenting an MRE under 10% for both systems WO and BO. Among the possible structures, D2 (dimer), F1 (trimer) and I1, and I3 (tetramers) presented MRE ≤ 13% for both cases.

09 BIOMASS FUELS↗
DOE OSTI2021

A physiological signal derived from sun-induced chlorophyll fluorescence quantifies crop physiological response to environmental stresses in the U.S. Corn Belt

Sun-induced chlorophyll fluorescence (SIF) measurements have shown unique potential for quantifying plant physiological stress. However, recent investigations found canopy structure and radiation largely control SIF, and physiological relevance of SIF remains yet to be fully understood. This study aims to evaluate whether the SIF-derived physiological signal improves quantification of crop responses to environmental stresses, by analyzing data at three different spatial scales within the U.S. Corn Belt, i.e. experiment plot, field, and regional scales, where ground-based portable, stationary and space-borne hyperspectral sensing systems are used, respectively. We found that, when controlling for variations in incoming radiation and canopy structure, crop SIF signals can be decomposed into non-physiological (i.e. canopy structure and radiation, 60% ~ 82%) and physiological information (i.e. physiological SIF yield, Φ F , 17% ~ 31%), which confirms the contribution of physiological variation to SIF. We further evaluated whether Φ F indicated plant responses under high-temperature and high vapor pressure deficit (VPD) stresses. The plot-scale data showed that Φ F responded to the proxy for physiological stress (partial correlation coefficient, r p = 0.40, p < 0.001) while non-physiological signals of SIF did not respond ( p > 0.1). The field-scale Φ F data showed water deficit stress from the comparison between irrigated and rainfed fields, and Φ F was positively correlated with canopy-scale stomatal conductance, a reliable indicator of plant physiological condition (correlation coefficient r = 0.60 and 0.56 for an irrigated and rainfed sites, respectively). The regional-scale data showed Φ F was more strongly correlated spatially with air temperature and VPD ( r = 0.23 and 0.39) than SIF ( r = 0.11 and 0.34) for the U.S. Corn Belt. The lines of evidence suggested that Φ F reflects crop physiological responses to environmental stresses with greater sensitivity to stress factors than SIF, and the stress quantification capability of Φ F is spatially scalable. Utilizing Φ F for physiological investigations will contribute to improve our understanding of vegetation responses to high-temperature and high-VPD stresses.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Chemical Mixing Induced by Internal Gravity Waves in Intermediate-mass Stars

Internal gravity waves can cause mixing in the radiative interiors of stars. We study this mixing by introducing tracer particles into 2D hydrodynamic simulations. Following the work of Rogers & McElwaine, we extend our study to different masses (3, 7, and 20 M ⊙ ) and ages (ZAMS, midMS, and TAMS). The diffusion profiles of these models are influenced by various parameters such as the Brunt–Väisälä frequency, density, thermal damping, the geometric effect, and the frequencies of waves contributing to these mixing profiles. We find that the mixing profile changes dramatically across age. In younger stars, we noted that the diffusion coefficient increases toward the surface, whereas in older stars the initial increase in the diffusion profile is followed by a decreasing trend. We also find that mixing is stronger in more massive stars. Hence, future stellar evolution models should include this variation. In order to aid the inclusion of this mixing in 1D stellar evolution models, we determine the dominant waves contributing to these mixing profiles and present a prescription that can be included in 1D models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Detection of surface water temperature variations of Mongolian lakes benefiting from the spatially and temporally gap-filled MODIS data

Lakes provide critical water resources for human activities and ecosystems, particularly in the Mongolian Plateau (MP), which is characterized by a dry climate and a harsh environment. As a region that is sensitive to anthropogenic warming, tracking lake surface water temperature (LSWT) changes in Mongolian lakes is crucial for understanding the consequences of a warming climate on lake ecosystems. However, the long-term monitoring of LSWT is restricted by the spatiotemporal gaps in the raw imagery of remote sensing-based land surface temperature (LST), e.g., the commonly used Moderate Resolution Imaging Spectroradiometer (MODIS) LST products. This study applied an improved gap-filling method by utilizing the discrete cosine transform-based penalized least squares (DCT-PLS) strategy in the spatial domain combined with the linear interpolation (LI) algorithm in the temporal domain. The method was applied to fill gaps in the LSWT imagery of 12 representative lakes across MP. The randomly sampled high-quality MODIS LSWT values in the spatial and temporal domains were excavated as false data gaps and considered “virtual true” validation datasets. The spatial validation results showed that the estimated LSWT for all the lake cases were comparable with the “virtual true” LSWT values, with the average values of the coefficient of determination, mean absolute error, mean square error, and root mean square error being 0.98, 0.38 °C, 0.45 °C, and 0.59 °C, respectively. Meanwhile, the error of nighttime LSWT results was relatively lower than that of daytime LSWT. For temporal interpolation validation, the LI algorithm exhibited relatively better performance and could more objectively indicate the variation in LSWT. Benefiting from the spatially and temporally well-constrained data, we analyzed the interannual and intra-annual change characteristics of the LSWTs of the 12 lakes. The long-term variations of annual and seasonal mean LSWTs in the 12 selected lakes exhibited no evident trends in 2000–2020, while presented apparent interannual fluctuations. The slight changes in the average LSWTs of the 12 selected lakes were in excellent synchronization with the surrounding LST derived from the reanalysis datasets, confirming the widely reported phenomenon of “global warming hiatus” that occurred in the early 21st century. This study improves the understanding of the LSWT variations in Mongolian lakes in response to global climate change. It has the potential to provide an effective approach for monitoring LSWT changes in other large-scale studies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Characterizing Seasonal Variation of the Atmospheric Mixing Layer Height Using Machine Learning Approaches

As machine learning becomes more integrated into atmospheric science, XGBoost has gained popularity for its ability to assess the relative contributions of influencing factors in the atmospheric boundary layer height. To examine how these factors vary across seasons, a seasonal analysis is necessary. However, dividing data by season reduces the sample size, which can affect result reliability and complicate factor comparisons. To address these challenges, this study replaces default parameters with grid search optimization and incorporates cross-validation to mitigate dataset limitations. Using XGBoost with four years of data from the atmospheric radiation measurement (ARM) (Southern Great Plains (SGP) C1 site, cross-validation stabilizes correlation coefficient fluctuations from 0.3 to within 0.1. With optimized parameters, the R value can reach 0.81. Analysis of the C1 site reveals that the relative importance of different factors changes across seasons. Lower tropospheric stability (LTS, ~0.53) is the dominant factor at C1 throughout the year. However, during DJF, latent heat flux (LHF, 0.44) surpasses LTS (0.22). In SON, LTS (0.58) becomes more influential than LHF (0.18). Further comparisons among the four long-term SGP sites (C1, E32, E37, and E39) show seasonal variations in relative importance. Notably, during JJA, the differences in the relative importance of the three factors across all sites are lower than in other seasons. This suggests that boundary layer development in the summer is not dominated by a single factor, reflecting a more intricate process likely influenced by seasonal conditions such as enhanced convective activity, higher temperatures, and humidity, which collectively contribute to a balanced distribution of parameter impacts. Furthermore, the relative importance of LTS gradually increases from morning to noon, indicating that LTS becomes more significant as the boundary layer approaches its maximum height. Consequently, the LTS in the early morning in autumn exhibits greater relative importance compared to other seasons. This reflects a faster development of the mixing layer height (MLH) in autumn, suggesting that it is easier to retrieve the MLH from the previous day during this period. The findings enhance understanding of boundary layer evolution and contribute to improved boundary layer parameterization.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Experimental study of the partitioning of some platinum group elements (Pd and Ir) between orthopyroxene and silicate melt

Past experiments and observations on natural samples have largely focused on the roles of olivine and chromite in controlling the behaviour of the platinum-group elements (PGE) during melting and solidification, whereas other phases, such as pyroxene, have gone largely uncharacterized. Here, to address this, experiments have been done to measure the partitioning of Pd (with a subset of results for Ir), between orthopyroxene and silicate melt at 1340 °C, 0.1 MPa and log fO 2 of FMQ - 1 to FMQ + 6 (FMQ = Fayalite-Magnetite-Quartz). The X-ray Absorption Near-Edge Structure (XANES) was measured in a subset of experiment glasses. Glass concentrations of Pd (corrected to unit Pd activity) increase from ~6 to ~650 ug/g over the fO 2 range of experiments. The slope of the solubility-fO 2 relation is consistent with Pd 1+ as the dominant oxidation state, with evidence for Pd 0 and Pd 2+ at the lowest and highest experiment fO 2 , respectively. Consistent with this result, the XANES reveal spectral features similar to Pd 0 and Pd 2+ spectral reference materials (specRM) at the most reduced and oxidized synthesis conditions, respectively. Other lines of evidence require the presence of a third melt species, here interpreted to be Pd 1+ . Values of orthopyroxene/melt partition coefficients for Pd (D Pd Opx/melt ) are 0.0051 (+/-0.006) at log fO 2 < ΔFMQ + 3, increasing with fO 2 to a maximum of 0.013 at ~FMQ + 6. Sodium partition coefficients, expected to be similar to Pd, range from 0.0061 (+/-0.00061) at FMQ + 3, increase to 0.007–0.009 at higher fO 2 , but with no clear systematic trend. A value for D Ir Opx-melt of ~0.6 was measured at ~FMQ + 4, indicating significantly more compatible behaviour for Ir relative to Pd. Partitioning results are interpreted in the context of the Blundy-Wood elastic strain model in which the variation in partitioning is related to ionic radius mismatch to an optimal crystallographic site size. Based on the trend in ionic radius with oxidation state, the estimated ionic radius of Pd 1+ in octahedral coordination is similar to Na 1+ , and comparison to previous orthopyroxene-melt partitioning experiments suggests D Pd1+ opx/melt and DNaopx/melt should be nearly identical, consistent with the results of this study. The ionic radius of VI-fold Pd 2+ is close to the optimal M2 site size, so an increased proportion of this species with fO 2 accounts for the larger values of D Pd opx/melt at the highest fO 2 investigated. The much larger partition coefficient for Ir is consistent with the presence of Ir 2+ , whose estimated ionic radius is close to Fe 2+ and Mg 2+ , as well as predictions for the optimal M1 site size. With the assumption that D Pd opx/melt = D Na opx/melt , combined with a revised value for the Pd content of the primitive mantle, a melting model is presented that better reproduces the Pd concentration of high degree melts from sulfide-free mantle sources.

58 GEOSCIENCES↗
DOE OSTI2022

Real-time measurement of parametric influences on the refractive index and length changes in silica optical fibers

In this paper, we present a simple cascaded Fabry-Perot interferometer (FPI) that can be used to measure in real-time the refractive index (RI) and length variation in silica optical fibers caused due to external physical parameters, such as temperature, strain, and radiation. As a proof-of-concept, we experimentally demonstrate real-time monitoring of temperature effects on the RI and length and measure the thermo-optic coefficient (TOC) and thermal expansion coefficient (TEC) by using the cascaded FPI within a temperature range of 21–486°C. The experimental results provide a TEC of 5.53 × 10 −7 /°C and TOC of 4.28 × 10 −6 /°C within the specified temperature range. Such a simple cascaded FPI structure will enable the design of optical sensors to correct for measurement errors by understanding the change in RI and length of optical fiber caused by environment parameters.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Measuring thermomechanical response of large-format printed polymer composite structures via digital image correlation

Large-format additive manufacturing (LFAM) is a branch of additive manufacturing (AM) research with the ability to create large structures typically measuring several meters in scale. LFAM is advantageous for tooling applications, not only because it offers the ability to create complex geometries not easily made using subtractive manufacturing processes, but the cost savings of pelletized feedstock used by these systems result in larger parts printed at faster speeds than traditional AM systems. Fiber reinforced polymer (FRP) is a commonly used feedstock material in LFAM structures because it reduces the distortion experienced during printing. However, FRP introduces highly anisotropic thermomechanical properties and contributes to a nonhomogeneous microstructure that can result in critical distortion of dimensions during tooling. Measuring the global thermomechanical response of LFAM structures requires a more representative method that accounts for not only anisotropic properties but also the nonhomogeneous nature of the final part. This is where traditional techniques to measure thermomechanical response, such as thermomechanical analysis (TMA), fall short as they assume homogeneity. This study evaluated the coefficient of thermal expansion (CTE) of LFAM structures as measured by TMA as compared to a novel digital image correlation oven (DIC Oven) system. The LFAM structures were made from 20 % by weight carbon fiber reinforced acrylonitrile butadiene styrene (CF-ABS). TMA measurements showed significant variations in CTE across a single LFAM bead, confirming the need for a global technique that captures overall thermomechanical response. The CTE values measured using the DIC Oven compared well to average TMA values obtained from localized measurements across the sample. The DIC Oven was also used to quantify the effects of different layer orientations on thermomechanical properties, which cannot be easily captured using TMA. A predictive model was also developed by using localized TMA values across an LFAM bead to predict the overall thermomechanical response of an LFAM structure.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Pozzolanic reactivity of natural pozzolans, ground glasses and coal bottom ashes and implication of their incorporation on the chloride permeability of concrete

It is well understood that the partial replacement of portland cement with pozzolans in sufficient proportions results in improved long-term performance of concrete. However, dwindling availability of the most widely used pozzolan in the construction industry, fly ash, is forecast due to the termination of coal-fired electricity power stations around the world, necessitating the exploration of alternative pozzolans. Experimental investigations are conducted to study the performance of traditional pozzolans such as fly ash, silica fume and natural pozzolans (metakaolin, pumice, perlite and lassenite) together with novel pozzolans, including ground glasses (high-alkali and low-alkali), and ground bottom ash. Twenty-one concrete mixtures are examined for compressive strength, electrical conductivity, electrical resistivity, chloride permeability, and chloride migration coefficient. The reactivity of pozzolans studied is also determined in terms of compressive strength in lime-pozzolan mortar using a modified lime-reactivity test method. Many of the materials tested demonstrate the potential to be used in place of fly ash. However, a wide variation in the performance of these materials is evident which highlights the need for a reliable test to determine the level of reactivity of a pozzolan. The results indicate that the 7-day strength in the modified lime-reactivity test provides a good indication of the pozzolanic reactivity of the material and of how the material can be expected to contribute to the strength and permeability of concrete. The use of electrical resistivity as an indicator of the performance of pozzolans in terms of strength development and chloride penetration resistance when used to partially replace portland cement in concrete is also discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Study of the film boiling heat transfer and two-phase flow interface behavior using image processing

Here, in this study, based on a small-scale quench test facility, the two-phase flow interface behavior during quench transients is visualized and analyzed utilizing an image processing framework. The high-fidelity experimental results obtained for two-phase flow in the current framework can support various studies both in the time domain and in the frequency domain. In particular, visualization of the data obtained from different heating surfaces under different test conditions are used to perform a full-scale transient 2-D vapor film reconstruction. The liquid-vapor interface variations in various heat transfer regimes as well as at the initial film breakup point can be directly obtained through the processed data. Moreover, the temporal variation of the interfacial wave frequency approaching quench is investigated in detail. Based on the high-resolution data obtained for the liquid-vapor interface, the detailed phase velocity and temperature profiles are obtained through theoretic analysis, based on which the film boiling heat transfer coefficient (HTC) can be determined. In addition, an improved film boiling HTC model is developed considering the effects of wall superheat, liquid subcooling temperature, vapor film thickness as well as fluid properties. The model is found to predict film boiling HTC well within 15% error.

42 ENGINEERING↗
DOE OSTI2026

Hanford Tank Waste Matrix Impact on Ion Exchange Performance Using Crystalline Silicotitanate

The removal of radiocesium from Hanford tank waste supernate is a critical step in preparing feed for low-activity waste immobilization. This study evaluated cesium ion exchange performance using crystalline silicotitanate (CST) media in a series of tests designed to evaluate the influence of waste matrix variability on capacity and kinetics. Tank waste supernate subsampled from five Hanford double-shell tanks encompassed a range of sodium, hydroxide, nitrate, and nitrite concentrations in order to assess the impact of feed variability on the performance of the ion exchange system. Both equilibrium and dynamic ion exchange tests were conducted to quantify cesium distribution coefficients and breakthrough behavior under prototypic operating conditions. Results indicated that effective cesium capacity varied by up to a factor of five across the matrices tested, with higher sodium concentrations significantly reducing uptake. Kinetic behavior was similarly matrix-dependent, with solution viscosity contributing to a twofold variation in mass-transfer rates. These results demonstrate the strong dependence of CST ion exchange performance on waste composition and must be incorporated into predictive models for future treatment system design and optimization.

Westesen, Amy M.↗
DOE OSTI2021

Variational Self-Consistent Theory for Beam-Loaded Cavities

A new variational theory is presented for beam loading in microwave cavities. The beam{field interaction is formulated as a dynamical interaction whose stationarity according to Hamilton’s principle will naturally lead to steady-state solutions that indicate how a cavity’s resonant frequency, Q and optimal coupling coefficient will detune as a result of the beam loading. A driven cavity Lagrangian is derived from first principles, including the effects of cavity wall losses, input power and beam interaction. The general formulation is applied to a typical klystron input cavity to predict the appropriate detuning parameters required to maximize the gain (or modulation depth) in the average Lorentz factor boost, (Δ γ). Numerical examples are presented, showing agreement with the general detuning trends previously observed in the literature. The developed formulation carries several advantages for beam-loaded cavity structures. It provides a self-consistent model for the dynamical (nonlinear) beam{field interaction, a procedure for maximizing gain under beam-loading conditions, and a useful set of parameters to guide cavity-shape optimization during the design of beam-loaded systems. Enhanced clarity of the physical picture underlying the problem seems to be gained using this approach, allowing straightforward inclusion or exclusion of different field configurations in the calculation and expressing the final results in terms of measurable quantities. Two field configurations are discussed for the klystron input cavity, using finite magnetic confinement or no confinement at all. Formulating the problem in a language that is directly accessible to the powerful techniques found in Hamiltonian dynamics and canonical transformations may potentially carry an additional advantage in terms of analytical computational gains, under suitable conditions

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

A Parametric, Data-Driven, Non-Intrusive Reduced-Order Model Framework for Crystal Plasticity Simulations of Voids

The influence of the internal structure at micrometer length scales on the deformation of polycrystalline materials can be effectively captured using crystal plasticity finite element methods (CPFEM). However, the complexity and nonlinearity of the deformation equations CPFEM solves demand significant computational power and resources to achieve accurate predictions, limiting its broader application. To address this challenge, we have identified a reduced-order representation of the complex data in order to establish a computationally efficient reduced-order models (ROM) and drastically reduce the computational expense of CPFEM. Specifically, in this work, we developed a parametric, data-driven, and non-intrusive ROM framework for CPFEM using proper orthogonal decomposition (POD) and sparse variational Gaussian process (SVGP) regression for single-crystal microstructures under tensile loading conditions. The developed protocol enables one to compress field into a latent/low-dimensional space described by principal component analysis (PCA) via the singular value decomposition (SVD) algorithm. As a result, the high-dimensional data are reduced to a significantly smaller amount of dimensions with POD bases and POD coefficients. Furthermore, we deployed an ensemble of SVGPs—extended from the classical Gaussian process (GP) regression for scalability and handling big data—in a massively parallel manner to train and predict latent POD coefficients using known POD bases from a set of previously obtained simulations results. Lastly, using the predicted POD coefficients, we reconstructed the full-field results and showed reasonable agreement compared with the true values obtained from running CPFEM. The developed framework is validated with a set of CPFEM simulations of a single embedded void in single-crystal aluminum alloy. While the framework is broadly applicable, this work specifically focuses on single-crystal microstructures, a single load case (e.g., tensile), and a specific void geometry (spherical).

Anisotropy↗
DOE OSTI2023

A combined experimental and crystal plasticity study of grain size effects in magnesium alloys

This work presents a method to incorporate the micro Hall-Petch equation into the crystal plasticity finite element (CPFE) framework accounting for the microstructural features to understand the coupling between grain size, texture, and loading direction in magnesium alloys. The effect of grain size and texture is accounted for by modifying the slip resistances of individual basal and prismatic systems based on the micro Hall-Petch equation. The modification based on the micro Hall-Petch equation endows every slip system at each microstructural point with a slip system-level grain size and maximum compatibility factor, which are in turn used to modify the slip resistance. While the slip-system level grain size is a measure of the grain size, the maximum compatibility factor encodes the effect of the grain boundary on the slip system resistance modification and is computed based on the Luster-Morris factor. The model is calibrated using experimental stress-strain curves of Mg-4Al samples with three different grain sizes from which the Hall-Petch coefficients are extracted and compared with Hall-Petch coefficients predicted using original parameters from previous work. The predictability of the model is then evaluated for a Mg-4Al sample with different texture and three grain sizes subjected to loading in different directions. The calibrated parameters are then used for some parametric studies to investigate the variation of Hall-Petch slope for different degrees of simulated spread in basal texture, variation of Hall-Petch slope with loading direction relative to basal poles for a microstructure with strong basal texture, and variation of yield strength with change in grain morphology. The proposed approach to incorporate the micro Hall-Petch equation into the CPFE framework provides a foundation to quantitatively model more complicated scenarios of coupling between grain size, texture and loading direction in the plasticity of Mg alloys.

Crystal plasticity↗