Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “BURNUP”

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 325 records · Page 18

Multi-element isotopic analysis of hot particles from Chornobyl

Microscopic fuel fragments, so-called “hot particles”, were released during the 1986 accident at the Chornobyl nuclear powerplant and continue to contaminate the exclusion zone in northern Ukraine. Isotopic analysis can provide vital information about sample origin, history and contamination of the environment, though it has been underutilized due to the destructive nature of most mass spectrometric techniques, and inability to remove isobaric interference. Recent developments have diversified the range of elements that can be investigated through resonance ionization mass spectrometry (RIMS), notably in the fission products. The purpose of this study is to demonstrate the application of multi-element analysis on hot particles as relates to their burnup, particle formation in the accident, and weathering. The particles were analysed with two RIMS instruments: resonant-laser secondary neutral mass spectrometry (rL-SNMS) at the Institute for Radiation Protection and Radioecology (IRS) in Hannover, Germany, and laser ionization of neutrals (LION) at Lawrence Livermore National Laboratory (LLNL) in Livermore, USA. Comparable results across instruments show a range of burnup dependent isotope ratios for U and Pu and Cs, characteristic of RBMK-type reactors. Results for Rb, Ba and Sr show the influence of the environment, retention of Cs in the particles and time passed since fuel discharge.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Mechanistic verification of empirical UO 2 fuel fracture models

Standard UO 2 fuel pellets used in light-water reactors fracture during irradiation due to the large thermal gradient in the radial direction. Over the decades, numerous researchers have explored fuel cracking from experimental and modeling points of view. To date, there have been both empirical and mechanistic approaches to predict the number of fragments that form in UO 2 . The empirical models only consider maximum power and burnup as inputs. Existing mechanistic approaches for normal operation have not accounted for irradiation effects. Here, this work employs a mechanistic fuel cracking model using the extended finite element method to explore radial crack formation while including a sensitivity analysis that accounts for the randomization of tensile strength within the fuel, the strength randomization criteria (uniform or volume-weighted Weibull), power ramping rates, computational mesh density, maximum power level, and irradiation (burnup) effects. The results indicate that the uncertainty in this mechanistic modeling approach envelopes the predicted values from three different empirical correlations in almost all cases. This means that, for computationally intensive analyses involving UO 2 fragmentation, the empirical correlations can be used. However, since the mechanistic calculations bound those of the empirical correlations, there is confidence in the applicability of the mechanistic approach developed in this work to generate a correlation for fuel types where limited data exists (e.g., doped-UO 2 , U 3 Si 2 ).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

An analysis of fluff formation in metallic fuel via data analyzes from EBR-II experiments and BISON fuel code modeling

During the operation of EBR-II, it was found that a highly porous structure (over 40% area fraction) formed at the top of several fuel columns. Previous work has shown that this structure, designated fluff in this paper, contains a significant fraction of fuel elements (e.g., U and Pu) which could potentially impact neutronics. This work aims in analyzing the formation mechanism of this microstructure so its impact can be incorporated into future metallic fuel modeling codes and algorithms. This paper details a preliminary examination into the formation mechanisms of fluff by performing qualitative and statistical analysis of EBR-II experimental data. Additionally, the operating conditions that have the greatest impact on fluff formation were determined based on this data set. Also, BISON fuel code simulations were used to help postulate potential fluff formation mechanisms. From this analysis it was found that the largest contributors to fluff formation were fuel burnup and composition, with fluff formation exhibiting a roughly linear positive correlation with increasing burnup and a negative correlation with increasing Pu content. It was also found that higher pin operating temperatures decreased fluff formation but only for U-10Zr fuels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Transmission electron microscopy characterization of the fuel-cladding chemical interactions in HT9 cladded U-10Zr fuel

Fuel cladding chemical interaction (FCCI) is a key phenomenon needs to be better understood to establish the design basis for U-10Zr metallic fuel performance. Characterizing the microstructure and chemical composition of FCCI at micron and sub-micron scale is critically important toward a more mechanistic understanding of FCCI phenomenon and its potential effects on cladding integrity and metallic fuel performance. This paper, by using transmission electron microscopy, investigated the FCCI region in HT9 cladded U-10Zr fuel irradiated to 5.7% FIMA burnup at a peak inner cladding temperature of 615 °C in Fast Flux Test Facility (FFTF). Four distinct layers are identified in the FCCI region. The migration of Fe into the fuel side leads to the formation of several U-Zr-Fe ternary phases, including ?-Fe0.5Zr0.32U0.18, e-Fe0.3Zr0.4U0.3, and ?-Fe0.06Zr0.23U0.71, mingled with UFe2, U6Fe, and U phase at various Fe penetration depth up to ~ 150 µm. On the cladding side, grain coarsening and significant lanthanides infiltration along grain boundaries are observed. Laves phase, (Fe,Cr)2(Mo,W), which typically does not exist in fresh HT9, is identified in a wide radial range in the cladding. The typical HT9 martensitic lath structure and pre-existing M23C6 precipitates disappear, partially or completely, depending on the radial distance from the fuel-cladding interface. Those microstructural and compositional changes could cause mechanical degradation in the HT9 cladding. The present characterization results will improve the understanding of FCCI phenomenon and facilitate the development of microstructure-informed FCCI modeling for metallic fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Possible impacts of Mo chemical banding and second phase impurities on the irradiation behavior of monolithic U-10Mo fuels

This study investigated the microstructural behavior of both full-size and mini-size monolithic U-10Mo fuel plates irradiated to high burnup with a focus on the evolution of the second phase impurities in monolithic U-Mo using scanning electron microscopy (SEM), energy dispersive spectroscopy (EDS), and wavelength dispersive spectroscopy (WDS). Key indicators of possible mechanical and thermal compromise include cracks, large fission gas porosity, and interconnection of fission gas pores. For the fission densities evaluated in this work (3.5 × 10 21 fissions/cm 3 -5.1 × 10 21 fissions/cm 3 ), fine porosity can develop along the UC phase boundary; however, the size of the fission gas pores is no more than those observed in the U-Mo fuel phase. Further, other inclusions such as Si-rich second-phase impurities found in the as-fabricated microstructure were difficult to resolve post-irradiation because they can become overshadowed by porosity development in the fuel phase. Additionally, the presence of a Fe-rich sublayer formed in the Zr diffusion layer during fabrication remained enriched in the Zr layer in the irradiated U-10Mo microstructures near the U-Mo/Zr interface; however, based on the burnup assessed in this study the identified impurities did not appear to contribute to notable microstructural degradation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Single-size and cluster dynamics modeling of intra-granular fission gas bubbles in UO 2

For this work, we perform simulations of intra-granular fission gas bubble evolution in UO 2 using both a relatively simple, computationally inexpensive single-size model and a detailed cluster dynamics model. Simulations encompass 36 experimental cases from 4 different databases, covering various temperature and burnup levels. We systematically compare results from the two models to each other and to post-irradiation experimental data of bubble average size and number density. Overall, the model-to-model comparisons reveal an excellent agreement across the set of simulations. This outcome indicates that, in spite of the underlying assumptions, the single-size model provides a good approximation of the complex physical behavior that is more rigorously described by the cluster dynamics model. Qualitatively, both models reproduce the trends of the experimental data with temperature and burnup correctly. Quantitatively, calculated results are either in good agreement with the data or within errors that appear consistent with the inherent uncertainties. Moreover, for the single-size model, we demonstrate and assess a multiscale approach whereby values for the fission gas atom diffusion coefficient from separate atomistic calculations are used. Systematic comparisons to experimental data point out a credible accuracy of the multiscale model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

TEM characterization of two variants of fuel cladding chemical interaction in a HT-9 Clad U-10Zr Fuel. Variant 1: FCCI with a Zr Rind

Here, this study investigated the fuel cladding chemical interaction (FCCI), a key factor that limits operational temperature and burnup, in an HT-9 clad U-10Zr nuclear fuel sample irradiated to a high burnup of 13.1 at.% at a time-averaged peak inner cladding temperature (PICT) of 530 °C. Previous results showed this fuel sample exhibited two distinct levels of FCCI at d. This paper analyzed the FCCI at an azimuthal position showing an interdiffusion layer of <10 µm using transmission electron microscopy to examine chemical and crystallographic nature of phases at the fuel-cladding interface at the nanoscale level. A ZrC layer and a Zr 3 Si phase were identified at the interface; these, along with the relatively low local temperature, potentially contributed to limit interdiffusion, behaving as inhibitors for deleterious interactions. Lanthanides (Ln) partially consumed the ZrC layer and interacted with Fe, forming a Zr-Ln compound and a (Zr,Ce)Fe 2+x phase while also infiltrating up to 4 µm into the cladding. Neither U nor Zr were observed in the cladding, whereas Fe diffused up to 3–5 µm in the fuel. Fe infiltration formed a ternary U-Zr-Fe ε-phase and likely promoted the precipitation of a Cr-rich α’ phase on the cladding interface. Additionally, a Cr-rich χ-phase, likely formed by the dissociation of pre-existing M 23 C 6 carbide precipitates, was identified about 2–5 µm from the fuel-cladding interface. Irradiation-induced nano-voids were also observed in the HT-9 bulk. These findings provide critical insights into FCCI mechanisms at representative irradiation conditions, essential for developing models simulating in-pile metallic fuel behaviors for next-generation reactors.

36 - MATERIALS SCIENCE↗

Post-irradiation examination of AGR-3/4 TRISO fuel compacts using three-dimensional X-ray computed tomography

The AGR-3/4 irradiation tests combined the third and fourth planned irradiation experiments in the US Department of Energy’s Advanced Gas Reactor (AGR) testing campaign of tri-structural isotropic (TRISO) fuel compacts. In this article, we present post-irradiation examination (PIE) using X-ray computed tomography (XCT) of two unirradiated and two irradiated compacts from the AGR-3/4 irradiation tests. The irradiated compacts studied (compact 7–1 and compact 12–4) represent the upper and lower limit of burnup within the AGR-3/4 irradiation experiment. This article presents a detailed quantitative analysis on the post-irradiation structure of TRISO fuel compacts. Various quantitative parameters including shape, size, and packing of kernels, and their spatial distribution, were utilized to gain insights into the structural changes caused by irradiation. The equivalent diameter and sphericity were found to increase and decrease, respectively, in irradiated compact 7–1 due to its higher burnup. Nearest neighbor distance between fuel kernels decreased after irradiation, suggesting irradiation-induced shrinkage of graphitic matrix. Furthermore, each compact in AGR-3/4 irradiation tests contained 20 designed-to-fail (DTF) fuel particles that were meant to act as a source of fission product release to the experiment test train. Furthermore, in the present work, all DTF fuel particles in the four compacts studied were identified, and it was found that they exhibited larger kernel swelling in compact 12–4 and smaller kernel swelling in compact 7–1, compared to the driver particles.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

ENDF/B-VIII.1: Updated Nuclear Reaction Data Library for Science and Applications

The ENDF/B-VIII.1 library is the newest recommended evaluated nuclear data file by the Cross Section Evaluation Working Group (CSEWG) for use in nuclear science and technology applications, and incorporates advances made in the six years since the release of ENDF/B-VIII.0. Among key advances made are that the 239 Pu file was reevaluated by a joint international effort and that updated 16,18 O, 19 F, 28–30 Si, 50–54 Cr, 55 Mn, 54,56,57 Fe, 63,65 Cu, 139 La, 233,235,238 U, and 240,241 Pu neutron nuclear data from the IAEA coordinated INDEN collaboration were adopted. Over 60 neutron dosimetry cross sections were adopted from the IAEA's IRDFF-II library. In addition, the new library includes significant changes for 3 He, 6 Li, 9 Be, 51 V, 88 Sr, 103 Rh, 140,142 Ce, Dy, 181 Ta, Pt, 206–208 Pb, and 234,236 U neutron data, and new nuclear data for the photonuclear, charged-particle and atomic sublibraries. Numerous thermal neutron scattering kernels were reevaluated or provided for the very first time. On the covariance side, work was undertaken to introduce better uncertainty quantification standards and testing for nuclear data covariances. The significant effort to reevaluate important nuclides has reduced bias in the simulations of many integral experiments with particular progress noted for fluorine, copper, and stainless steel containing benchmarks. Data issues hindered the successful deployment of the previous ENDF/B-VIII.0 for commercial nuclear power applications in high burnup situations. These issues were addressed by improving the 238 U and 239,240,241 Pu evaluated data in the resonance region. The new library performance as a function of burnup is similar to the reference ENDF/B-VII.1 library.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Conceptual design study of neutron detectors for safeguards measurement of an irradiated pebble

Nuclear material control and accounting (MC&A) of pebble-bed reactors (PBRs) is challenging because a PBR utilizes hundreds of thousands of identical, unmarked pebbles that are continuously recirculated through the core. To develop tools that enable the implementation of international safeguards, especially in the context of MC&A of spent pebbles, we designed and simulated three neutron detection concepts to determine fissile content in individual pebbles: a differential die-away (DDA) detector, a californium interrogation prompt neutron (CIPN) detector, and a passive neutron albedo reactivity (PNAR) detector using Monte Carlo calculations. Burnup calculations were performed on the spent pebbles from the PBMR-400 classic PBR. The varying neutron and gamma source terms, and isotopic compositions in the spent pebbles calculated at various burnup levels were used in the neutron detector models. DDA was found to be sensitive to the number of passes a pebble has had through the core and to the fissile content contained in a spent pebble. Optimization in the DDA design further increased the neutron count rates and thus reduced counting uncertainty. Meanwhile, passive neutron counting using the same detector body could distinguish pebbles with different numbers of passes, but its response was dominated by neutron-emitting actinides and was not sensitive to fissile content. On the other hand, the PNAR technique was not viable for a single pebble but performed reasonably for a 27-pebble array, which suggested potential use for verification measurements of containers filled with 27 or more spent pebbles.

CIPN↗

Spent fuel nondestructive assay integrated characterization from active neutron, passive neutron, and passive gamma

Spent nuclear fuel comprises a wide range of irradiated isotopic material compositions, and characterization through nondestructive measurements is beneficial in verifying declared parameters before the fuel is placed in storage, final disposal, and/or reprocessed. This paper discusses results from three nondestructive assay instruments, including passive gamma, passive neutron, and active neutron methods, that measured fifty spent fuel assemblies at the Clab interim storage facility in Sweden. Here, integrated analysis of the measurements from the three different instruments allowed parametric assessments of cooling time, burnup, neutron multiplication, fissile mass, initial enrichment, and decay heat of each individual fuel assembly. Passive gamma measurements were found to be the most beneficial in predicting cooling time, passive neutron for determining burnup, active neutron in estimating initial enrichment, and both passive and total neutron for multiplication correlations. Fissile mass was best estimated using any combination of any two of the instruments such that corrections for isotopic changes in the fuel could be accounted for with the first set of measurements and multiplication of the assembly in the second. In conclusion, the nondestructive assay technologies demonstrated through this effort enhance the characterization of spent nuclear fuel assemblies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Multi-modal tomographic imaging system for poolside characterization of nuclear test fuels: Design considerations and studies

Testing and qualification of advanced nuclear fuels involves an iterative process of prototyping, in-pile irradiation testing, and in-situ or ex-situ examination. Fuel restructuring and fission product migration during burnup are among the most important aspects of fuel evolution that affect several important performance characteristics such as heat removal, accident tolerance, and fission product retention. Pool-side non-destructive characterization techniques provide fuel developers with tools to understand fuel evolution at different points of burnup. A design for a compact, submersible, and multi-modal gamma-ray tomography system for imaging irradiated nuclear fuel is presented here. Detector selection, collimator geometry and fabrication, mechanical design, imaging protocol and acquisition speed are discussed. Modeling calculations show that sub-millimeter resolution can be achieved in both transmission computed tomography images as well as in emission computed tomography images in a matter of hours. Several design compromises and fabrication challenges are discussed which should be taken into consideration for future submersible gamma-ray tomography instruments.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

BISON fuel performance modeling optimization for experiment X447 and X447A using axial swelling and cladding strain measurements

With the recent need to qualify new reactor designs such as the Versatile Test Reactor (VTR), fuel performance calculations need to be performed to determine safety criteria of the proposed designs. In order to validate the fuel performance results obtained by a fuel performance code, BISON, for new reactor designs, legacy fuel from EBR-II and FFTF MFF with Post -Irradiation Examination (PIE) data need to be used as validation cases to benchmark models. Here in this work, BISON has been paired with the Fuels Irradiation & Physics Database (FIPD) and IFR Materials Information System (IMIS) to supply PIE data for comparison with simulations of EBR-II experiments X447/X447A. X447/X447A were assessed by implementing models for Fuel Cladding Chemical Interaction (FCCI) within BISON and optimizing the friction coefficient between the fuel surface and the cladding, the anisotropic swelling factor, and the HT9 first thermal creep scalar (which scales the first term in the HT9 creep equation) to best match the PIE axial fuel swelling height and cladding profilometry for all pins in X447/X447A. The optimal values were found using a generic algorithm developed to select different values for the three parameters until end criteria was met and error couldn’t be reduced further. The BISON-simulated cladding profilometry was evaluated using Standard Error of the Estimate (SEE) to account for the profile shape of the cladding profilometry. Optimal values for the friction coefficient, anisotropic fuel swelling factor, and HT9 first thermal creep scalar were found to best fit the BISON simulation results to the PIE measurements found in IMIS and FIPD. Improvements to current models are suggested to account for the underprediction of fuel swelling at low burnups and the overprediction of fuel swelling at higher burnups observed for the axial fuel swelling height. Although two pins in EBR-II X447/X447A (DP70 and DP75) were known to fail due to FCCI, none of the pins simulated in BISON reached a cumulative damage fraction (CDF) above 0.008 with FCCI correlations coupled in the BISON simulations. The error estimate generated for all pins in X447/X447A using optimal values was 209 µm, which is deemed acceptable.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Assessing fuel cycle design options for advanced microreactors

Microreactor (MR) technology have recently gained traction due to their reduced construction cost compared to traditional large reactors. Their small physical geometry results in increased neutron leakage and deteriorated neutron economy which limits the achievable burnup of MR fuel. Refueling strategy currently discussed in MR design community is the battery-like refueling which is involve replacing the entire reactor core at the end of their operational cycle. Here, this paper investigates partial core refueling as means to increase MR fuel burnup and improve fuel cycle economy as a result. Two fuel cycle design strategies have been successful. These are the partial core refueling and partial core refueling combined with assembly rotation. Both strategies have been tested to operate a heatpipe-cooled MR referred to as the eVinci-like core for 72 Effective Full Power Months (EFPMs). The partial core refueling showed to reduce the refueling requirements by 18% compared to the battery-like refueling strategy. This corresponds to 10.5% reduction in the Net Present Value (NPV) of fuel cost over the 72 EFPMs of operation. The partial core refueling combined with rotation case reduces the refueling requirements by 29% (corresponding to 14.5% reduction in NPV of fuel cost) compared to the battery-like refueling strategy.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Refinement and demonstration of a coupled BISON-Griffin workflow for designing targeted TRISO transient experiments in TREAT

The U.S. nuclear industry is expected to deploy tristructural isotropic (TRISO) particle fuel technologies for commercial reactors within the next decade. In previous work, we defined a preliminary transient design space for TRISO fuels, identified potential gaps in the available data, and began to develop multiphysics modeling tools that could be applied to design targeted Transient Reactor Test Facility (TREAT) experiments to fill these gaps. Here, this work builds on that foundation by (1) updating BISON fuel performance and Griffin reactor physics models to reflect the current TREAT experiment tube and capsule designs, (2) coupling the codes to improve the accuracy and usability of the transient design analyses, and (3) demonstrating their use over an expanded design space that includes fuel burnup. The simulated mechanical responses of the TRISO particles were complex functions of fission product accumulation, fission gas release, and irradiation-induced dimensional change in the pyrolytic carbon layers. The predicted tangential stresses in the particles' silicon carbide layers were least compressive for preheated tests involving fresh fuels but remained compressive throughout the ranges of temperature, heat rate, and burnup considered in this work. Finally, comparisons between the potential TREAT transients and historical test reactor irradiations showed that the TREAT tests would produce significantly lower average energy deposition rates, yielding less severe transients with greater relevance to near-term commercial applications. The use of these predictive capabilities has the potential to increase the value of each test, improving the overall efficiency and cost-effectiveness of transient testing for TRISO and other advanced fuels.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Estimation of post-burst vibration-induced fuel dispersal thresholds following LOCA-simulant testing

The nuclear industry intends to increase burnup past the regulatory limit of 62 GWd/MTU, and to be successful, the dispersal behavior of fragmented high-burnup fuel during loss-of-coolant accidents must be better understood. Part of this behavior is fuel fragmentation, relocation, and dispersal. This work reports on testing concerning the dispersal of fuel from rods following burst and pressure blowdown through the burst. The Fuel Relocation Induced by Accident Recovery (FRIAR) system was used to examine surrogate material and pre-burst cladding specimens to evaluate the effects of vibrational loading, mixture composition, and burst size on post-burst dispersal. Vibrational loading resulted in maximum dispersal around 15 Hz when the loading was performed with peak-to-peak oscillation amplitudes at 0.5 and 1 mm. Here, under this aggressive vibrational loading, dispersal was found to become significant as burst widths increased to approximately 5 mm and larger. However, below 4 mm widths, little material dispersed out of the opening, even at 1 mm amplitudes and for more dispersive mixtures (smaller average fragments). At and above 7 mm widths, dispersal was rapid for all mixtures. Finally, dispersal in the threshold region (4–5 mm wide bursts) strongly depended on the mixture. For mixtures with large fragments, dispersal was minimized to a few fragments, as the large fragments proved to be effective at minimizing dispersal. Conversely, mixtures with smaller fragments showed increased dispersal in this region.

Burst geometry↗

Design of a first-of-A-kind instrumented advanced test reactor irradiation Capsule experiment for In situ thermal conductivity measurements of metallic fuel

Metallic fuel undergoes dramatic microstructural changes early in life due to fission gas swelling until ~2–3 at% burnup which affects the conductivity of the material, however the evolution of metallic fuel thermal conductivity during this early phase burnup has never been successfully measured in situ. The Irradiated Material Properties Accelerated Characterization Test (IMPACT) experiment will be the first in a series of experiments to irradiate advanced nuclear metallic fuel specimens with novel embedded thermal conductivity probes in ATR. In the current work the IMPACT experiment final design and supporting analysis is reported in detail. Results are evaluated for various reactor operational conditions to meet the functional requirements of the experiment. Finally, the first iteration of this IMPACT experiment will provide data regarding thermal properties evolution in uranium-zirconium (U10Zr) fuel, but this experiment vehicle is envisioned for future advanced fuels and structural materials irradiations in ATR.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Identifying chemically similar multiphase nanoprecipitates in compositionally complex non-equilibrium oxides via machine learning

Abstract Characterizing oxide nuclear fuels is difficult due to complex fission products, which result from time-evolving system chemistry and extreme operating environments. Here, we report a machine learning-enhanced approach that accelerates the characterization of spent nuclear fuels and improves the accuracy of identifying nanophase fission products and bubbles. We apply this approach to commercial, high-burnup, irradiated light-water reactor fuels, demonstrating relationships between fission product precipitates and gases. We also gain understanding of the fission versus decay pathways of precipitates across the radius of a fuel pellet. An algorithm is provided for quantifying the chemical segregation of the fission products with respect to the high-burnup structure, which enhances our ability to process large amounts of microscopy data, including approaching the atomistic-scale. This may provide a faster route for achieving physics-based fuel performance modeling.

36 MATERIALS SCIENCE↗