Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variational analysis”

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

First-Principles Cost Estimation of a Sodium Fast Reactor Nuclear Plant

A multi-tiered cost analysis is performed to estimate full costs of a nuclear power plant (NPP) based on sodium-cooled fast reactor (SFR) technology. To address the lack of fully transparent cost estimations from past undertakings for NPPs, we have developed a detailed and first-principles-based cost estimate for a generalized SFR NPP. Our intent is to achieve a high degree of transparency with our cost assumptions and develop a cost model that is flexible and easily extendable to variations in NPP design and other nuclear reactor types. Furthermore, we strive to achieve a clear organization of costs and complete identification of key cost drivers based on first principles. To this end, the cost results of our analysis as given in Table 24 and Table 25 are organized and categorized into a code of accounts (COA) under development at Idaho National Laboratory (INL). Varying degrees of first-principles methods are employed, such as design for manufacture and assembly® (DFMA® ), to elucidate costs in all process levels of the plant equipment, buildings and site structures, personnel, and other miscellaneous but significant cost elements. These approaches have been successfully applied in past cost analysis projects and are designed to enable rapid and flexible cost estimation. Application of these techniques for evaluating NPP costs is similar in concept to the full, detailed estimation of construction and fabrication costs determined in a later stage of NPP development. Note that our approach tries to avoid use of other past analysis results and data such as those from the legacy Energy Economic Data Base (EEDB) Program, as these resources are based on historical NPP costs and thus may not be indicative of new reactor technologies or construction and fabrication/manufacturing techniques. However, we provide a comparison of our SFR NPP cost results in Table 90 against those included in the EEDB for a representative pressurized water reactor (PWR).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

CdTe Core: Final Technical Report (FTR)

CdTe is presently the cost-leading thin-film PV technology, directly competing with Si at scale, even when domestically manufactured. While an impressive technology, its efficiency remains much below the detailed balance limit with the largest cause due to its low photovoltage and fill factor. To realize gains, the carrier concentration, minority carrier lifetime, and interface recombination all need to be improved simultaneously over historic levels. Using a new defect chemistry (group V doping instead of copper) has been identified as a viable route using single crystal systems. This project focused on implementing this new defect chemistry in scalable, polycrystalline thin-film photovoltaic CdTe devices with tasks focusing improvements to the front interface, absorber, and rear interface as well as capability development & stakeholder engagement. The goal of the project was to establish a strategy using devices, test structures, detailed characterization, and modeling to quantify the sources of losses in state-of-the-art CdTe photovoltaic devices. Using this strategy and advanced synthesis, losses at the front interface, absorber, and rear interface were worked on in parallel. The final objective was to significantly improve the voltage deficit in CdTe devices to enable improvements in photovoltage and efficiency that can be implemented by industry in the near-term. Over the course of the project, the team developed new characterization techniques, analysis, and modeling which were then applied to state-of-the-art materials generated internally and collaboratively. In particular to enable rapid progress, NREL worked closely with First Solar where NREL grew complete devices as well as ones that interleaved process steps where First Solar had completed different steps such as absorber growth or absorber growth and activation using their baseline methods. Using detailed characterization and analysis including photoemission, photoluminescence, and scanning probe techniques enabled understanding of the loss pathways and area for improvements in our own and First Solar s materials. Ultimately, this contributed to the first series of new world record CdTe efficiencies since 2016, culminating in a 23.1% certified cell that was P-doped along with As-doped cells of similar performance. Internally, NREL improved the statistical variation in baseline As-doped devices and improved average photovoltage by over 100 mV. This was done through an improvement in absorber quality, changed front interface, and improved back contact. In addition to materially improving the fabrication processes at NREL, characterization, analysis, and modeling were developed and disseminated. NREL also played a pivotal role in community building over the course of this project working closely with the Cadmium Telluride Accelerator Consortium. NREL worked in a series of collaborations with academic and industry partners, leveraging knowledge and innovations from this project, as well as helped organize a series of workshops to ensure rapid progress in the field. Working closely with the academic community has led to a dissemination of knowledge; working with First Solar as increased US competitiveness First Solar expanded domestic production to ~10 GW and opened new facilities.

14 SOLAR ENERGY↗
DOE OSTI2026

Telomere-to-telomere assemblies of chromosome 10 reveal complex adaptive variation of 3-ketoacyl-CoA-synthases in Populus trichocarpa likely driven by Helitrons

The model woody plant Populus trichocarpa displays an atypical alkene-diverse wax cuticle likely driven by copy number variation (CNV) of 3-ketoacyl-CoA synthases ( KCS ), which has been difficult to confirm with short-read assemblies. Long-read sequencing enables the development of telomere-to-telomere resources to detect cryptic variation, including CNVs, which are currently missed. Integrating this information can improve genomic prediction for breeding and provide insights into the evolutionary basis of important traits. Our analysis of 78 long-read haplotypes from chromosome 10 identified more than twice as many KCS genes as previously reported, and numerous intragenic non-synonymous substitutions. Random Forest predictive models highlighted the importance of Potri.010G079500 in producing very long chain alkenes; however, its absence did not predict previously reported alkene-deficient phenotypes. Instead, alkene levels are best predicted by the combinations of KCS copies. Additionally, amino acid substitutions clustered around ligand and donor binding pockets, suggesting they contribute to differing wax cuticle composition. Finally, each KCS gene and copy was linked to a Helitron transposon. A phylogenetic analysis suggests Helitrons are the evolutionary mechanism for generating KCS tandem arrays. Long-read generated telomere-to-telomere assemblies of P. trichocarpa chromosome 10 revealed large-effect loci critical to genetic studies that are unattainable from short-reads. This new resource produced novel insights into genome structure and function, and a novel mechanism for generating tandem gene duplication. Our results highlight that, given current challenges in annotation and assembly, detailed and focused long-read sequences are key to interpreting complex genomic regions that contain tandem copy number variants.

09 BIOMASS FUELS↗
DOE OSTI2025

S AP F LOWER : an automated tool for sap flow data preprocessing, gap-filling, and analysis using deep learning

Sap flow, a critical process in plant water use and ecosystem water cycles, is often measured using thermal dissipation probes (TDP) due to their ease of installation and continuous data collection. However, sap flow data frequently include noise, outliers, and gaps, creating challenges for analysis and requiring substantial manual processing. We developed S AP F LOWER , a tool that automates data preprocessing, model training, gap-filling, sapwood area scaling and modeling, and water use analysis. It integrates autocleaning, machine learning and deep learning models (e.g. random forest, Gaussian process regression, long short-term memory (LSTM), bidirectional LSTM (BiLSTM)), and efficient workflows to process sap flow data. S AP F LOWER can remove over 90% of noisy data while preserving legitimate variations and achieve high accuracy in gap-filling based on user-determined parameters. Random forest, LSTM, and BiLSTM models reduced root mean square error to 10% or less for long-term gaps. Model training and prediction can be performed efficiently within seconds. S AP F LOWER significantly enhances the efficiency and accessibility of TDP data analysis by automating complex tasks, enabling researchers without programming expertise to employ advanced techniques. Future improvements will focus on species-specific corrections for TDP and support for additional measurement methods. S AP F LOWER is openly available on GitHub (https://github.com/JiaxinWang123/SapFlower) and Zenodo (doi: 10.5281/zenodo.13665919).

ecosystem water balance↗
DOE OSTI2023

An inexact semismooth Newton method with application to adaptive randomized sketching for dynamic optimization

In many applications, one can only access the inexact gradients and inexact hessian times vector products. Thus it is essential to consider algorithms that can handle such inexact quantities with a guaranteed convergence to solution. An inexact adaptive and provably convergent semismooth Newton method is considered to solve constrained optimization problems. In particular, dynamic optimization problems, which are known to be highly expensive, are the focus. A memory efficient semismooth Newton algorithm is introduced for these problems. The source of efficiency and inexactness is the randomized matrix sketching. Further, applications to optimization problems constrained by partial differential equations are also considered.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

The influence of cooling rate on condensation of iron, aluminum, and uranium oxide nanoparticles

Fundamental observations of particle size distributions are needed to develop models that predict the fate and transport of radioactive materials in the atmosphere following a nuclear incident. The extent of material transport is influenced by the time scales of particle formation processes (e.g., condensation, coagulation). In this study, we investigated the influence of cooling time scales on size distributions of uranium, aluminum, and iron oxide particles that are synthesized separately under identical run conditions inside the controlled environment of an argon plasma flow reactor. Two distinct temperature distributions are imposed along the flow reactor by varying the argon flow rate downstream of the plasma torch. The vaporized reactants of uranium, aluminum, and iron are cooled from about 5000K to 1000K before they are collected on silicon wafers for ex situ scanning electron microscope analysis. The microscope images show that the sizes of the largest aluminum and iron oxide particles heavily depend on the cooling time scales, whereas significant size variation with cooling rate is not observed for uranium oxide particles. In addition, the size distribution of aluminum oxide particles exhibits the broadest range among all three metal oxides studied. We performed simulations of particle size distributions using a kinetic model that couples gas phase oxidation chemistry with particle formation processes, including nucleation, condensation, and coagulation. The model results demonstrate the strong sensitivity of particle size distribution to different cooling histories (i.e., temperature vs residence time) along the flow reactor. In conclusion, the kinetic model also helps identify directions for future research to improve the predictions.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A Bayesian inferencing framework for ultrasound wave speed measurements in metal additive manufacturing

Process-related changes during metal additive manufacturing introduce microstructural variability in the material properties of printed parts, directly affecting component reliability. Accurate estimation of these property variations with part performance are essential for quality assurance. Ultrasound testing offers a non-destructive means to estimate mechanical properties and detect defects; however, conventional analysis methods often neglect the influence of microstructural variability, limiting their effectiveness. Here, this research presents a Bayesian inference technique for quantifying wave speed uncertainty from ultrasound measurements of metal additive manufactured parts. By integrating prior ultrasound data with a Bayesian model, the proposed approach generates posterior density estimates of wave speed that systematically account for manufacturing-induced variability and uncertainty. The novelty of this research lies in applying a Bayesian framework to analyze experimental ultrasound measurements within the context of metal additive manufacturing variability. The method enhances the accuracy of wave speed estimation by 64%, defect position by 50% and increases confidence associated with wave speed variance by 30% across different porosity levels, thereby providing a robust foundation for improved decision-making and increased reliability in additively manufactured components.

Additive manufacturing↗
DOE OSTI2020

Radiation shielding analysis of the small-angle X-ray scattering flight tube end station of APS upgrade project

The CSSI beamline of the APSU will have a 20.5 m long, 2.75 m diameter SS end station designed to accept pink beam. In this work, the radiation shielding analysis of this end station is analyzed using STAC8 and FLUKA codes. Effect of mirror properties such as reflectivity, coating, inclination and roughness as well as variation in the steel composition are also studied. 8 mm thick SS is found to be adequate if the source is defined by a pinhole and the white beam is reflected with two mirrors inclined at 3.0 mradian or more. The dose rates from monochromatic beams are found to be below the desired levels when the XOP calculated bandwidth (BW) are used while the use of a flat 0.1% BW requires the mirror inclinations to be above 2.5 mradian or more. STAC8 results are consistently higher than the FLUKA results by about 1.5–2.0 times.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Modeling and simulation of earthquake soil structure interaction excited by inclined seismic waves

Presented is an application of wave potential formulation (WPF) together with domain reduction method (DRM) to modeling earthquake soil structure interaction (ESSI) behavior in horizontally layered ground under inclined incident seismic waves. Wave potential formulation is used to develop a spatially varying, inclined seismic wave field from incident Primary (P) and Secondary (S) waves that propagate through layered ground. Developed seismic wave field is then used to develop effective forces for Domain Reduction Method that are then used for analyzing ESSI response of a soil structure system. Developed methodology, called WPF-DRM, is verified using analytic solution for a free field response of layered ground subjected to inclined incident waves. Developed WPF-DRM methodology is illustrated through analysis of an ESSI response of a deeply embedded structure, a small modular reactor (SMR) subjected to incident S wave polarized in vertical plane (SV) with variation in inclinations and frequencies. Presented example highlights the influences of incident wave inclination and frequency on ESSI response of analyzed SMR.

58 GEOSCIENCES↗
DOE OSTI2022

Regional Variation in Light-Duty Plug-in Electric Vehicle Emissions

Automobile electrification is viewed as one of the key requirements for deep decarbonization of transportation. However, the exact magnitude of decarbonization from switching from a combustion of petroleum-based fuels to grid-derived electricity is a complex issue. Variations in local grid mixes lead to differing carbon intensities of electricity in different regions. This complexity is compounded by geographical differences in vehicle characteristics. This analysis considers these geographic differences to produce a historical assessment of the fuel-cycle carbon emissions of plug-in electric vehicles (PEVs) in the United States from 2011 to 2021. We find that PEVs in the United States decreased in electricity-derived carbon intensity, from 187 grams per mile to 110 grams per mile from 2011 to 2021 due to improvements in the electric grid and vehicle efficiency. Through 2021, PEVs emitted a total of 13.6 million metric tons of greenhouse gas (GHG) emissions in the United States, including both emissions from electricity and gasoline, while driving a total of 100 billion miles. Of this total, 8.39 million metric tons was due to electricity consumption. Over the same distance, comparable gasoline vehicles would have emitted between 29 and 41 million metric tons of GHG, leading to a reduction of 15 to 27 million metric tons. Relative to national-level assessments that do not factor in this regional difference, these electric vehicles resulted in 13% greater reductions in GHG emissions than previously calculated. In the future, announced goals to further decarbonize the electricity sector will continue to reduce the emissions of PEVs, including reducing emissions rates for vehicles already on the road, further improving the benefits relative to the unchanging emissions of contemporaneous gasoline vehicles.

02 PETROLEUM↗
DOE OSTI2022

DISTRIBUTION OF ION-EXCHANGEABLE RARE EARTH ELEMENTS AND CRITICAL MINERALS IN COAL-HOSTED CLAYS OF THE POWDER RIVER BASIN, WY

Increasing demand for a more robust domestic supply of rare earth elements (REE) and critical minerals (CM) has led to significant investigation into unconventional sources. Coal and coal byproducts from the Powder River Basin (PRB) of Wyoming are potential sources of REE but lack complete characterization of mineral hosts including ion-adsorption clays (IACs). IACs, which provide a substantial portion of the global REE supply, form when hydrated REE complexes are physiosorbed to negatively charged sites on clay surfaces. Typical IAC minerals include kaolinite, illite, and smectite. REE complexes are easily leached from IACs via ion exchange, making them economically attractive. However, research regarding the partitioning behaviors of strategic metals into PRB coal-hosted clays is limited. This project aims to quantify the modes of occurrence and distribution of REE and CM in clay minerals derived from PRB coals. Strategic metals investigated in this project have historic production or promising deposits within the PRB and/or high anticipated supply risk, and include: REE (lanthanides and Y), V, Ti, platinum group elements (PGE; Pt, Pd, Rh, Ru, Ir, Os), Al, Co, Ga, Zr, and W. Semi-quantitative compositional data from portable x-ray fluorescence spectroscopy (pXRF) on coal core from the PRB is consistent with the known correlation between REE and Al content, suggesting increased REE abundance will be found in clay-rich units. pXRF data also identifies elevated V, Ti, Co, and Zr correlated with Al. Preliminary x-ray diffraction (XRD) data indicates kaolinite, illite, and dickite are abundant minerals in overburden and clay-rich partings. Quantitative analysis via inductively coupled plasma mass spectrometry (ICP-MS) provides insights into REE and CM fractionation in coals as we compare trace metal occurrence to variations in mineralogy. Additionally, we apply scanning electron microscopy (SEM) and energy dispersive x-ray spectroscopy (EDS) data to correlate clay size and morphology to trace metal occurrence. Our findings supplement a broader effort to couple REE and CM extraction with existing domestic coal production and contribute to the fundamental understanding of adsorption processes that occur in low-temperature basin environments.

Stuart, Sophia↗
DOE OSTI2024

Collaborative Research: Vlasov–Maxwell Simulations to resolve electron heating and dissipation, in quasiperpendicular shocks (Final Technical Report)

The project is to study electron heating at quasi-perpendicular shocks. We have analyzed measurements from the Magnetosphere Multiscale (MMS) spacecraft and particle-in-cell simulations to address the question. One key discovery is that waves and are excited in the shock transition region, which further develop into complicated structures and thin current layers, which contribute to accelerate and heat electrons. We further investigated the excitation and evolution of these waves and structures, and their association with ion and electron dynamics. In addition, the analysis technique of Vlasov equation terms is further applied to study electron dynamics at the dayside magnetopause, to understand the effect of density, velocity and temperature variations on electron distributions and kinetics.

58 GEOSCIENCES↗
DOE OSTI2023

Assessment and Non-Destructive Evaluation of the Influence of Residual Solvent on a Two-Part Epoxy-Based Adhesive Using Ultrasonics

Polymers are increasingly being used in higher demanding applications due to their ability to tailor the properties of structures while allowing for a weight and cost reduction. Solvents play an important role in the manufacture of polymeric structures since they allow for a reduction in the polymer’s viscosity or assist with the dispersion of fillers into the polymer matrix. However, the incorrect removal of the solvent affects both the physical and chemical properties of polymeric materials. The presence of residual solvent can also negatively affect the curing kinetics and the final quality of polymers. Destructive testing is mainly performed to characterize the properties of these materials. However, this type of testing involves using lab-type equipment that cannot be taken in-field to perform in situ testing and requires a specific sample preparation. Here, a method is presented to non-destructively evaluate the curing process and final viscoelastic properties of polymeric materials using ultrasonics. In this study, changes in longitudinal sound speed were detected during the curing of an aerospace epoxy adhesive as a result of variations in polymer chemistry. To simulate the presence of residual solvent, samples containing different weight percentages of isopropyl alcohol were manufactured and tested using ultrasonics. Thermogravimetric analysis was used to show changes in the decomposition of the adhesive due to the presence of IPA within the polymer structure. Adding 2, 4, and 6 wt.% of IPA decreased the adhesive’s lap shear strength by 40, 58, and 71%, respectively. Ultrasonics were used to show how the solvent influenced the curing process and the final sound speed of the adhesive. Young’s modulus and Poisson’s ratio were determined using both the longitudinal and shear sound speeds of the adhesive. Using ultrasonics has the potential to non-invasively characterize the quality of polymers in both an in-field and manufacturing settings, ensuring their reliability during use in demanding applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Investigation into Sizing Photovoltaic with Energy Storage for Off-Grid Transactive Scenarios

In this study, a novel sizing methodology was developed for centralized and interconnected operating strategies of transactive microgrids and several variables were investigated including starting month, initial charge of battery, load variability, unit cost of solar panels and energy storage, number of systems, climate, and required reliability to determine their effect on total cost. The centralized strategy improved cost by seven to ten percent compared to the isolated strategy in every case. The interconnected strategy saved an incremental amount of money consistently compared to the isolated standard. The number of connected systems was not a strong effect. It was thought that increasing the number of systems would increase the benefit of energy sharing. Climate zones studied (“Cold”; “Hot-Dry/Mixed Dry”; “Mixed Humid”; and “Cold but with lower solar irradiation”) showed a large variation on cost with the Hot-Dry/Mixed Dry being the least expensive and Cold, with lower solar irradiation being the most expensive. Cost sensitivity analysis was performed showing that the unit cost of solar has a greater effect on the total cost. Required reliability of power, measured in outage hours, exhibited an inverse relation with cost.

14 SOLAR ENERGY↗
DOE OSTI2023

A variational encoder–decoder approach to precise spectroscopic age estimation for large Galactic surveys

Constraints on the formation and evolution of the Milky Way Galaxy require multidimensional measurements of kinematics, abundances, and ages for a large population of stars. Ages for luminous giants, which can be seen to large distances, are an essential component of studies of the Milky Way, but they are traditionally very difficult to estimate precisely for a large data set and often require careful analysis on a star-by-star basis in asteroseismology. Because spectra are easier to obtain for large samples, being able to determine precise ages from spectra allows for large age samples to be constructed, but spectroscopic ages are often imprecise and contaminated by abundance correlations. Here we present an application of a variational encoder–decoder on cross-domain astronomical data to solve these issues. The model is trained on pairs of observations from APOGEE and Kepler of the same star in order to reduce the dimensionality of the APOGEE spectra in a latent space while removing abundance information. The low dimensional latent representation of these spectra can then be trained to predict age with just ∼1000 precise seismic ages. We demonstrate that this model produces more precise spectroscopic ages (∼ 22 per cent overall, ∼ 11 per cent for red-clump stars) than previous data-driven spectroscopic ages while being less contaminated by abundance information (in particular, our ages do not depend on [α/M]). We create a public age catalogue for the APOGEE DR17 data set and use it to map the age distribution and the age-[Fe/H]-[α/M] distribution across the radial range of the Galactic disc.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Doppler Lidar Measurements of Wind Variability and LLJ Properties in Central Oklahoma during the August 2017 Land–Atmosphere Feedback Experiment

Abstract Low-level jets (LLJs) are an important nocturnal source of wind energy in the U.S. Great Plains. An August 2017 lidar-based field-measurement campaign [the Land–Atmosphere Feedback Experiment (LAFE)] studied LLJs over the central SGP site in Oklahoma and found nearly equal occurrences of the usual southerly jets and postfrontal northeasterly jets—typically rare during this season—for an opportunity to compare the two types of LLJs during this month. Southerly winds were stronger than the northeasterlies by more than 4 m s −1 on average, reflecting a significantly higher frequency of winds stronger than 12 m s −1 . The analysis of this dataset has been expanded to other SGP Doppler-lidar sites to quantify the variability of winds and LLJ properties between sites of different land use. Geographic variations of winds over the study area were noted: on southerly wind nights, the winds blew stronger at the highest, westernmost sites by 2 m s −1 , whereas on the northeasterly flow nights, the easternmost sites had the strongest wind speeds. Lidar measurements at 5 sites during August 2017, contrasted to the 2016–21 summertime data, revealed unusual wind and LLJ conditions. Temporal hodographs using hourly averaged winds at multiple heights revealed unorganized behavior in the turbulent stable boundary layer (SBL) below the jet nose. Above the nose, some nights showed veering qualitatively similar to inertial oscillation (IO) behavior, but at amplitudes much smaller than expected for an IO, whereas other nights showed little veering. Vertical hodographs had a linear shape in the SBL, indicating little directional shear there, and veering above, resulting in a hook-shaped hodograph with height. Significance Statement Doppler-lidar measurements at the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) were used to quantify the variability of winds and low-level jet (LLJ) properties between five sites of different land use and wind regimes across this area. Knowledge of wind and LLJ structure and dynamics is important for many applications, and strong southerly LLJ winds at night are an important resource for wind energy. The analysis of multiyear (2016–21) summertime LLJ parameters provided insight into the LLJ climatology in this part of the Great Plains.

Meteorology & Atmospheric Sciences↗
DOE OSTI2021

Local Multimodal Electro-Chemical-Structural Characterization of Solid-Electrolyte Grain Boundaries

Typical models of polycrystalline ionic materials treat the grain boundary properties as single valued, without consideration of the full range of values that define the macroscopically measured average. Here a unique experimental platform suitable for local multimodal characterization of individual grain boundaries in bicrystal fibers is reported. A variation of three orders of magnitude in the grain boundary conductivity of ceria is observed, as measured across six individual bicrystals by both alternating current impedance spectroscopy and direct current (D.C.) linear sweep voltammetry. Nonlinear behavior of the D.C. measurements is consistent with resistance due to a space charge effect. Time-of-flight secondary ion beam spectroscopy reveals a correlation between grain boundary resistance and the concentration of impurities Si, Al, and Ca segregated at the grain boundaries, although the bulk concentrations of these impurities are negligible. Electron backscatter diffraction analysis of the crystal orientations suggests a correlation between the misorientation across the grain boundaries and grain boundary resistance. These correlations point towards a grain boundary resistance that arises from impurity-generated space charge effects and variations in impurity concentration and hence resistivity driven by the energetics of impurity segregation to grain boundaries of differing surface energies.

36 MATERIALS SCIENCE↗