Search NASA⌕ Search

SEARCH · Search NASA

Results for “Material Modeling”

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

Attenuation of waves in a viscoelastic peridynamic medium

The effect of spatial nonlocality on the decay of waves in a dissipative material is investigated. The propagation and decay of waves in a one-dimensional, viscoelastic peridynamic medium is analyzed. Both the elastic and damping terms in the material model are nonlocal. Waves produced by a source with constant amplitude applied at one end of a semi-infinite bar decay exponentially with distance from the source. The model predicts a cutoff frequency that is influenced by the nonlocal parameters. A method for computing the attenuation coefficient explicitly as a function of material properties and source frequency is presented. Here, the theoretical results are compared with direct numerical simulations in the time domain. The relationship between the attenuation coefficient and the group velocity is derived. It is shown that in the limit of long waves (or small peridynamic horizon), Stokes’ law of sound attenuation is recovered.

97 MATHEMATICS AND COMPUTING↗

Being Scientific: Saryu Fensin: Sparking innovations in both materials science and shock physics

As an enterprising materials scientist, Saryu Fensin leverages both experiments and simulations to improve understanding of how materials behave when subjected to extreme forces. Her aim is to use the complementary techniques to design and develop innovative materials for use in extreme conditions such as stockpile-focused applications. “To certify our stockpile with confidence, we need to understand the role a material’s microstructure plays in determining its dynamic properties,” Saryu says. “I’m hoping to solve this puzzle and improve our predictive capability to model material response in extreme conditions. Achieving this kind of specific controlled functionality of materials could have a big impact on the Laboratory’s national security science mission.” Her tools include molecular dynamics, Monte Carlo methods, gas gun experiments, metallography and various microscopy techniques. She looks forward to the prospect of new capabilities that could let her "see" inside materials.

36 MATERIALS SCIENCE↗

Foreign body impact event damage formation in composite structures

This report discusses a methodology that can be used to assess the effect of foreign body impacts on composite structural integrity. The described effort focuses on modeling the effect of a central impact on a 5 3/4 inch filament wound test article. The discussion will commence with details of the material modeling that was used to establish the input properties for the analytical model. This discussion is followed by an overview of the impact assessment methodology. The progress on this effort to date is reviewed along with a discussion of tasks that have yet to be completed.

Bucinell, Ronald B.↗

Multiscale Modeling of Structurally-Graded Materials Using Discrete Dislocation Plasticity Models and Continuum Crystal Plasticity Models

A multiscale modeling methodology that combines the predictive capability of discrete dislocation plasticity and the computational efficiency of continuum crystal plasticity is developed. Single crystal configurations of different grain sizes modeled with periodic boundary conditions are analyzed using discrete dislocation plasticity (DD) to obtain grain size-dependent stress-strain predictions. These relationships are mapped into crystal plasticity parameters to develop a multiscale DD/CP model for continuum level simulations. A polycrystal model of a structurally-graded microstructure is developed, analyzed and used as a benchmark for comparison between the multiscale DD/CP model and the DD predictions. The multiscale DD/CP model follows the DD predictions closely up to an initial peak stress and then follows a strain hardening path that is parallel but somewhat offset from the DD predictions. The difference is believed to be from a combination of the strain rate in the DD simulation and the inability of the DD/CP model to represent non-monotonic material response.

Saether, Erik↗

Test and Analysis Correlation for Sandwich Composite Longitudinal Joint Specimens

The NASA Composite Technology for Exploration (CTE) project is tasked with evaluating methods to analyze and manufacture composite joints for potential use in block upgrades to the Space Launch System (SLS) launch-vehicle structures such as the Payload Attach Fitting (PAF). To perform this task, the CTE project has initiated test and analysis correlation studies for composite joints under various loading conditions. Herein, NASA-developed numerical models are correlated with the experimental results from a series of tension tests. Pretest strain results matched the far-field test data well, but did not capture the nonlinear response in the vicinity of the joint. A refined pretest analytical model was modified to represent progressive failure of the specimens at failure locations observed during the experimental tests. The nonlinear strain response from this progressive failure model predicted the delamination failure load within 15% of the test data, but underpredicted the nonlinearity of the strain response. Further study of composite material models that account for the nonlinear shear response of fabric composites is recommended for the composite joint structures considered in this paper.

Mason, Brian H.↗

Screening sensitivity analysis of a PWR fuel assembly FEM structural model

Within the present project, computational methods to describe the static mechanical response of PWR fuel assemblies (FA) in the reactor core are developed. For this purpose a structural model of a PWR fuel assembly has been implemented with the Finite Element Method (FEM) using the ANSYS APDL Mechanical code. The purpose of the model is to investigate the long-term creep deformation of the fuel assemblies in the reactor core, commonly known as fuel assembly bow. FA bow is considered to be a complex process with a large number of influencing mechanisms and many only roughly known boundary conditions. Uncertainty and sensitivity analyses are a common way to investigate the predictive power and reliability of such complex models. Hence, as a first approach to understand the importance of the many model parameters, a qualitative sensitivity analysis of the model was conducted using the screening method proposed by Morris. In the first part of the paper, the FA structural model and its principal parameters are presented. The second part presents first the different simulation experiments used for the screening analyses. Based on the model description, relevant input parameters are then identified and a typical range of values is determined. The sensitivity analysis is performed in two subsequent steps. First, the model is screened for influential stiffness parameters with the objective of stiffening the FA structure. Second, the influence of the uncertain boundary conditions and material models on the FA creep deformation is tested. For the considered FA design it is concluded that the influence of the involved uncertainties is comparable to the influence of the parameters which can increase the stiffness. An integral approach accounting for these uncertainties is hence necessary to predict FA bow reliably in order to take additional measures for a further reduction of FA bow amplitudes. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

The effectiveness of D 2 pellet injection in reducing intra-ELM and inter-ELM tungsten divertor erosion rates in DIII-D during the Metal Rings Campaign

Abstract Edge localized modes (ELMs) in H-mode plasmas erode plasma-facing components (PFCs) and lead to impurities in the core, reducing confinement. This study analyzes D 2 pellet injection on the DIII-D fusion experiment used as an ELM mitigation technique applied during the 2016 tungsten Metal Rings Campaign to reduce W erosion during ELMs. The 400.9 nm photon wavelength line emission intensity of tungsten atoms (WI) filterscope channels and Langmuir probes were used to infer the gross erosion rate of tungsten-coated tiles installed in the divertor of DIII-D. D 2 mass injection rates ranging from 34 to 41 arbitrary units (A.U.) and no D 2 injection resulted in a similar total W erosion rate during ELMs (intra-ELM). On average, results show a 29% increase in the total gross W erosion rate with intermediate mass injection rates (∼13–23 A.U.) compared to the no pellets and the highest injection rate cases. On average, the fast D 2 mass injection rate cases had 15% less erosion in the inter-ELM phase than the case with no pellets. Generally, higher D 2 mass injection rates increased the ELM frequency, and the highest injection rates reduced the average erosion per ELM and fractional carbon impurities at the top of the pedestal by nearly 40% when compared to the no-pellet case. As expected, a higher D 2 pellet injection rate led to a higher plasma density and lower plasma temperature in the divertor. Additionally, an increasing divertor inter-ELM plasma electron density directly correlated to more frequent pellet injection and a decrease in both the average gross intra-ELM W erosion and the total gross intra-ELM W erosion rate. Simulations of intra-ELM erosion using the ‘free-streaming plus recycling model’ (FSRM) underestimate W erosion during pellet injection by about 30% on average. The discrepancies between the experimental measurements and the FSRM intra-ELM W erosion predictions are postulated to be due to C/W material mixing. A simple analytic mixed-material model is presented and results in better agreement with the experimental data. These results highlight the importance of incorporating the effects of a mixed-material layer in the analysis of PFC erosion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fuel Performance Analysis of Fast Flux Test Facility MFF-3 and -5 Fuel Pins Using BISON with Post Irradiation Examination Data

Using the BISON fuel-performance code, simulations were conducted of an automated process to read initial and operating conditions from the Pacific Northwest National Laboratory (PNNL) database and reports, which contain metallic-fuel data from the Fast Flux Test Facility (FFTF) MFF Experiments. This work builds on previous modeling efforts involving 1977 EBR-II metallic fuel pins from experiments. Coupling the FFTF PNNL reports to BISON allowed for all 338 pins from MFF-3 and MFF-5 campaigns to be simulated. Each BISON simulation contains unique power and flux histories, axial power and flux profiles, and coolant-channel flow rates. Fission-gas release (FGR), fuel axial swelling, cladding profilometry, and burnup were all simulated in BISON and compared to available post-irradiation examination (PIE) data. Cladding profilometry, FGR, and fuel axial swelling simulation results for full-length MFF metallic pins were found to be in agreement with PIE measurements using FFTF physics and models used previously for EBR-II simulations. The main two peaks observed within the cladding profilometry were able to be simulated, with fuel-cladding mechanical interaction (FCMI), fuel-cladding chemical interaction (FCCI), and thermal and irradiation-induced creep being the cause. A U-Pu-Zr hot-pressing model was included in this work to allow pore collapse within the fuel matrix. This allowed better agreement between BISON-simulated cladding profilometry and PIE measurements for the peak caused by FCMI. This work shows that metallic fuel models used to accurately represent fuel performance for smaller EBR-II pins may be used for full-length metallic fuel, such as FFTF MFF assemblies and the Versatile Test Reactor (VTR). As new material models and PIE measurements become available, FFTF MFF assessment cases will be reassessed to further BISON model development.

36 MATERIALS SCIENCE↗

Fission Product Transport in TRISO Particles and Pebbles

This document demonstrates completion of the goals described in the technical narrative of the FOA project titled: ”Modeling and Simulation Development Pathways to Accelerating KP-FHR Licensing” regarding fission product transport in the Kairos-proposed fuel pebble by INL and Kairos Power. Showcased in this report are code developments and simulations in BISON that extend the state of the art in computation and understanding of fission product transport in a TRISO fuel particle and pebble. These enhancements lay the foundation for making predictions of fission product transport that can be used as input in the fuel licensing process. This was achieved by installing existing fuel material models originally used in PARFUME, developing a new failure probability method that is efficient and multi-dimensional, employing material homogenization, and expanding verification and validation simulations to demonstrate the efficacy of the work. All this work is leveraged to spotlight the main deliverable; a three-dimensional model and corresponding demonstration simulation of a pebble, which will serve as the starting point for models used to predict fission product release.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

High-temperature deformation behavior of the extruded Ti-22Al-25Nb alloy fabricated by powder metallurgy

Highlights: • The high-temperature deformation behavior of the extruded Ti-22Al-25Nb alloy was analyzed. • A true strain-compensated Arrhenius-type constitutive model was established. • The processing maps of true strains 0.2, 0.4, and 0.6 were constructed to optimize the processing parameters. • Microstructural evolution and deformation mechanism in different conditions were explained. - Abstract: The hot-deformation behavior of the extruded Ti-22Al-25Nb alloy was investigated by compression testing in the 1213–1333 K temperature range, under a 0.001–1.0 s{sup −1} strain-rate range. The activation energy, Q, at a deformation strain of 0.3, was calculated to be 574.80 kJ mol{sup −1}. A strain-compensated Arrhenius-type constitutive model was established for this alloy, with an average absolute relative error (AARE) and correlation coefficient (R) of 7.64% and 0.994, respectively. Based on the dynamic materials model (DMM), processing maps were developed by combining the power dissipation and instability maps. An instability region was identified in the 1213–1278 K temperature range and 0.022–1.0 s{sup −1} strain-rate range. The microstructure of specimens deformed under different conditions were analyzed to characterize the corresponding deformation mechanisms. The main softening mechanisms in the stability regions were dynamic recrystallization (DRX) and dynamic globularization. The grain size of the specimen deformed at 1313 K/0.001 s{sup −1} was approximately 15.2 μm.

36 MATERIALS SCIENCE↗

Mesoscale Numerical Characterization of Kevlar and Carbon–Kevlar Hybrid Plain-Woven Fabric Compression Behavior

A mesoscale finite element model (FEM) of the plain-woven fabric representative elementary volume (REV/unit cell) is developed. The developed FEM is used to predict the nonlinear compaction behavior of monolithic and intraply hybrid fabric made of high-strength carbon and Kevlar fibers. For geometrical modeling of the fabric REV, TexGen, an open-source geometrical modeling software (developed at Nottingham University), is utilized, and for mechanical characterization, FE simulation tool Abaqus{sup ®} is used. Geometrical nonlinearity of the fabric, including periodic boundary conditions and yarn interactions, is incorporated in the REV FE model. Material nonlinearity caused by transversely isotropic behavior of fibrous yarns is implemented through the user material subroutine form in FE simulation along with local orientation assigned to each yarn in the REV. Material properties of fibers are extracted from the literature and supplier data. The compression behavior of monolithic Kevlar and hybrid carbon–Kevlar (carbon warp and Kevlar weft yarns, C–K) fabric is predicted using the proposed FE model. FEM-predicted results were found to have good agreement with the experimentally obtained ones. Results depict that C–K hybrid fabrics have improved compression strength. Hybridization effect is governed by yarns mutual interaction, yarn area at crossovers, and transverse in-plane properties of the interyarn hybridized fabric. The proposed FE model can be used to predict the mechanical behavior of woven fabrics with different weaving patterns, unit cell geometry, material properties, and loading conditions.

36 MATERIALS SCIENCE↗

Fission Product Transport in TRISO Particles and Pebbles

This document demonstrates completion of the goals described in the technical narrative of the FOA project titled: ”Modeling and Simulation Development Pathways to Accelerating KP-FHR Licensing” regarding fission product transport in the Kairos-proposed fuel pebble by INL and Kairos Power. Showcased in this report are code developments and simulations in BISON that extend the state of the art in computation and understanding of fission product transport in a TRISO fuel particle and pebble. These enhancements lay the foundation for making predictions of fission product transport that can be used as input in the fuel licensing process. This was achieved by installing existing fuel material models originally used in PARFUME, developing a new failure probability method that is efficient and multi-dimensional, employing material homogenization, and expanding verification and validation simulations to demonstrate the efficacy of the work. All this work is leveraged to spotlight the main deliverable; a three-dimensional model and corresponding demonstration simulation of a pebble, which will serve as the starting point for models used to predict fission product release.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Dynamic chemical expansion of thin-film non-stoichiometric oxides at extreme temperatures

Actuator operation in increasingly extreme and remote conditions requires materials that reliably sense and actuate at elevated temperatures, and over a range of gas environments. Design of such materials will rely on high-temperature, high-resolution approaches for characterizing material actuation in situ. In this work, we demonstrate a novel type of high-temperature, low-voltage electromechanical oxide actuator based on the model material Pr x Ce 1-x O 2-δ (PCO). Chemical strain and interfacial stress resulted from electrochemically pumping oxygen into or out of PCO films, leading to measurable film volume changes due to chemical expansion. At 650 °C, nanometre-scale displacement and strain of >0.1% were achieved with electrical bias values <0.1 V, low compared to piezoelectrically driven actuators, with strain amplified fivefold by stress-induced structural deflection. Finally, this operando measurement of films ‘breathing’ at second-scale temporal resolution also enabled detailed identification of the controlling kinetics of this response, and can be extended to other electrochemomechanically coupled oxide films at extreme temperatures.

36 MATERIALS SCIENCE↗

A New Method for Optimizing Hot Processing Parameters of Mg-6.0Zn-0.5Mn-0.5Er Alloy Based on Kinetic Model of Dynamic Recrystallization and Processing Map

Hot deformation behavior of the cast-homogenized Mg-6.0Zn-0.5Mn-0.5Er alloy was studied using dynamic recrystallization (DRX) kinetic model and processing map. The compressing tests were conducted in the temperature range of 250-450 °C and strain rates of 0.001-1 s{sup −1}. According to the evolution of microstructures, under lower strain rates, the main DRX mechanism of Mg-6.0Zn-0.5Mn-0.5Er alloy is continuous DRX (CDRX); twinning induced DRX (TDRX) and CDRX both become the main DRX mechanism under higher strain rates. The DRX kinetics model of Mg-6.0Zn-0.5Mn-0.5Er alloy is calculated as X{sub DRX}=1−exp[−1.9463(((ε−ε{sub c})/(ε{sup ∗}))){sup 1.4608}], which is corresponding to the microstructure evolutions of DRX under different deformation conditions. The contour map of DRX was proposed based on the calculation results of DRX kinetics model. The processing maps are constructed to predict processing parameters of the alloy, and the predictability was evaluated combining with the contour map of DRX and dynamic materials model (DMM) processing map. It is deduced from the microstructures evolution and processing map that the optimum processing domain is mainly at 380-450 °C and 0.01-0.001 s{sup −1}, and 410-420 °C and 0.01-1 s{sup −1}.

36 MATERIALS SCIENCE↗

Elastic-plastic deformation of molybdenum single crystals shocked to 12.5 GPa: Crystal anisotropy effects

Here, to understand crystal anisotropy effects on shock-induced elastic-plastic deformation of molybdenum (Mo), results from high-purity single crystals shocked along [110] and [111] orientations to an elastic impact stress of 12.5 GPa were obtained and compared with the [100] results previously reported [A. Mandal and Y. M Gupta, J. Appl. Phys. 121, 045903 (2017)]. Measured wave profiles showed a time-dependent response, and strong anisotropy was observed in the elastic wave attenuation with the propagation distance, elastic limits, shock speeds, and overall structure of the wave profiles. Resolved shear stresses on {110} $\langle$111$\rangle$ and {112} $\langle$111$\rangle$ slip systems provided insight into the observed anisotropy in elastic wave attenuation and elastic limits and showed that shear stresses, and not longitudinal stresses, are a better measure of strength in shocked single crystals. Under shock compression, resolved shear stresses at elastic limits were comparable to the Peierls stress of screw dislocations in Mo. Elastic wave attenuation was rapid when shear stresses were larger than the Peierls stress. Large differences in the elastic limits under shock and quasi-static loading are likely a consequence of the large Peierls stress value for Mo. Numerically simulated wave profiles, obtained using the dislocation-based plasticity model described in the [100] work, showed good agreement with all measured wave profiles but could not differentiate between the {110} $\langle$111$\rangle$ and {112} $\langle$111$\rangle$ slip systems. Overall, experimental results and corresponding numerical simulations for the three crystal orientations have provided a comprehensive insight into shock-induced elastic-plastic deformation of Mo single crystals, including the development of a continuum material model.

36 MATERIALS SCIENCE↗

Associative Flow Rule Used to Include Hydrostatic Stress Effects in Analysis of Strain-Rate-Dependent Deformation of Polymer Matrix Composites

designing reliable composite engine cases that are lighter than the metal cases in current use. The types of polymer matrix composites that are likely to be used in such an application have a deformation response that is nonlinear and that varies with strain rate. The nonlinearity and the strain-rate dependence of the composite response are due primarily to the matrix constituent. Therefore, in developing material models to be used in the design of impact-resistant composite engine cases, the deformation of the polymer matrix must be correctly analyzed. However, unlike in metals, the nonlinear response of polymers depends on the hydrostatic stresses, which must be accounted for within an analytical model. By applying micromechanics techniques along with given fiber properties, one can also determine the effects of the hydrostatic stresses in the polymer on the overall composite deformation response. First efforts to account for the hydrostatic stress effects in the composite deformation applied purely empirical methods that relied on composite-level data. In later efforts, to allow polymer properties to be characterized solely on the basis of polymer data, researchers at the NASA Glenn Research Center developed equations to model the polymers that were based on a non-associative flow rule, and efforts to use these equations to simulate the deformation of representative polymer materials were reasonably successful. However, these equations were found to have difficulty in correctly analyzing the multiaxial stress states found in the polymer matrix constituent of a composite material. To correct these difficulties, and to allow for the accurate simulation of the nonlinear strain-rate-dependent deformation analysis of polymer matrix composites, in the efforts reported here Glenn researchers reformulated the polymer constitutive equations from basic principles using the concept of an associative flow rule. These revised equations were characterized and validated in an experimental program carried out through a university grant with the Ohio State University, wherein tensile and shear deformation data were obtained for a representative polymer for strain rates ranging from quasi-static to high rates of several hundred per second. Tensile deformation data also were obtained over a variety of strain rates and fiber orientation angles for a representative polymer matrix composite composed using the polymer.

Goldberg, Robert K.↗

(Project 14-6770) An Investigation to Establish Multiphysical Property Dataset of Nuclear Materials Based on in-situ Observations and Measurements

In-core nuclear materials including fuel pins and cladding materials fail due to issues including corrosion, mechanical wear, and pellet cladding interaction. In most such scenario microstructure dependent and corrosioninduced chemistry dependent property changes significantly affect performance of cladding, pellet, and housing. Emphasis of this work was on replace conventional pellet-cladding material models with a new straingradient viscoplasticity model that is informed by transmission electron microscopy (TEM) based measurements and by nanomechanical Raman spectroscopy (NMRS) based measurements. The TEM measurements are quantitative in nature and therefore reveal stress-strain relations with simultaneous insights into mechanisms of deformation at nanoscale. The NMRS measurements reveal the similar information at mesoscale along with additional information on relating local microstructural stresses with applied stresses. The resulting information is used to fit constants in the strain gradient viscoplasticity model as well as to validate one. During TEM measurements, a micro-electro-mechanical system based setup was developed with mechanical actuation, sensing, heating, and electrical loading. Contrary to post-mortem analysis or qualitative visualization, this setup combines direct visualization of the mechanisms behind deformation with measurement of stress, strain, thermal and electrical properties. The unique research philosophy of visualizing the microstructure at high resolution while measuring the properties led to fundamental understanding in grain size and temperature effects on measured mechanical properties such as fracture toughness. A key contribution is the role of mechanical loading boundary conditions to deconvolute the insitu TEM based nanoscale and NMRS based mesoscale data to bulk behavior. First the literature based pellet cladding mechanical interaction model based on the work of Retel’s and Williamson’s in literature work to predict tempurature and stress distribution in cladding and pellet at normal operating condition was analyzed. Later the data was fitted to find constants for a viscoplastic strain gradient model. The developed model still needs to be refined and calibrated using various experimental results. That remains the focus of future work. Overall, a major thrust of the work was therefore on active control of the microstructure (grain size, defect density and types) exploiting the multi-physics coupling in materials. In particular, using experiments the synergy of current density, mechanical stress and temperature were studied to annihilate defects and recrystallize grains. The developed model is being examined for implementation in BISON. Multiple invited talks, international journal publications, and conference publications were performed by students supported on this work. Another output is support multiple PhD and masters thesis students who will be an important asset for future basic nuclear research. Future Work Recommendations: A nuclear reactor operates under significant variations of thermal loads due to energy cycling and mechanical loads due to constraint effects. Significant thermal and chemical diffusion takes place at the pallet-cladding level. While the proposed work established new experimental approach and new dataset for Zircaloy-4, the irradiation level was in the range of 1-2 dpa. Samples with higher dpa need to be examined. Therefore, a continual of support of the performed work is essential. Currently, these are the only experiments that can measure the produced data. The work also needs to be extended to different fuel types and cladding types such as SiC and FeCrAl based claddings. A combination of datasets for these materials can then be used to analyze accurately predict behavior of critical pellet cladding systems in accident scenario with high heat flux and high thermal loads. This is a BIG unknown as if now.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A phase field model for dislocations in hexagonal close packed crystals

In this work, an application of a phase field formulation suitable for modeling the motion of individual partial and full dislocations in hexagonal close packed (HCP) crystals is presented. The formulation incorporates periodic potentials for glide on the distinct HCP slip systems, which are informed here by density functional theory (DFT). The model is applied to simulate the dissociation process starting from an unstable perfect dislocation and ending at its final equilibrium structure for different slip planes and in different HCP metals. The structural characteristics that are predicted for these dislocations include the partial Burgers vectors, dissociation distances, core widths of the partials, and any asymmetries in these quantities. Mg is selected as one of the model materials since its dislocations are the most well studied and it is nearly elastically isotropic. For Mg, it is shown that the predictions for dissociation distances agree with those reported previously by atomic-scale calculations, including density functional theory, for dislocations on the basal < a >, prismatic < a >, and pyramidal type II slip systems. Furthermore, phase field model results are also presented for dislocations in Ti and Zr, which we find develop distinctively different equilibrium structures than Mg.

36 MATERIALS SCIENCE↗