Search NASA⌕ Search

SEARCH · Search NASA

Results for “FRACTURES”

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

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

At least 1,387 records · Page 77

Imaging and Endovascular Treatment of Bleeding Pelvic Fractures: Review Article

Pelvic fractures are potentially life-threatening injuries with high mortality rates, mainly due to intractable pelvic arterial bleeding. However, concomitant injuries are frequent and may also be the cause of significant blood loss. As treatment varies depending on location and type of hemorrhage, timely imaging is of critical importance. Contrast-enhanced CT offers fast and detailed information on location and type of bleeding. Angiography with embolization for pelvic fracture hemorrhage, particularly when performed early, has shown high success rates as well as low complication rates and is currently accepted as the first method of bleeding control in pelvic fracture-related arterial hemorrhage. In the current review imaging workup, patient selection, technique, results and complications of pelvic embolization are described.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Energy Budget of Brittle Fracturing in Granite Under Stress Relaxation and Creep

Creep and relaxation are the two major time-dependent fracturing processes in rocks. While a considerable amount of research has been done in understanding these two mechanisms, critical gaps remain regarding how different energy components evolve during time-dependent fracturing processes in rocks. In this study, a series of relaxation and creep experiments were conducted on prismatic Barre granite specimens in the laboratory to estimate the energy budget of brittle fracturing in granite. For the input energy, the work done by the machine (W) is calculated and for the output energy the radiated seismic energy (E R ), released in the form of acoustic emissions (AEs), is calculated as the only measurable output energy component in the conducted experiments. The low-frequency plateau (Ω 0 ) and corner frequency (f 0 ) for each AE waveform was estimated by fitting the observed AE spectra with the theoretical spectra using the Omega model. These parameters were used to estimate the seismic moments (M 0 ) based on the radiation pattern for the double couple (shear) and non-double-couple (non-shear) events. The range of f 0 and M 0 varied from 150 to 750 kHz and 10 –4 to 10 –1 N m, respectively. Moment magnitude (M w ) varied in a wider range from –9 to –6 in creep and –8.5 to –7 in relaxation. Stress drops (Δ σ ) and source radius (r) were estimated for the AEs using Brune’s model. Here, the results report on three primary observations: (1) the effects of different source mechanisms on the estimated source parameters showed that M 0 and Δ σ were higher for DC events as compared to NDC in both relaxation and creep. (2) The radiation efficiency in the case of creep is 70% higher as compared to relaxation and, (3) the stress drop estimated in relaxation and creep demonstrated a breakdown in scaling with the seismic moment.

58 GEOSCIENCES↗

Preface to the special volume on the third Sandia Fracture Challenge

The mounting reliance on computational simulations to predict all aspects of the lifecycle of a mechanical system, from fabrication to failure, has prompted the mechanics community to selfassess its abilities to perform those predictions. Benchmark problems in mechanics that compare simulations that use different computational approaches with experiments have sprung up lately, including the NIST AM-Bench looking at additively manufactured (AM) materials (https://www.nist.gov/ambench),the Contact-Mechanics Challenge (Miiser, 2017) considering adhesion between two nominally flat surfaces, Numisheet providing semiannual benchmarking activities in sheet metal forming (http://numisheet2018.org),and the Sandia Fracture Challenge (SFC) (Boyce, 2014 and Boyce, 2016) investigating ductile failure. The previous SFCs have shown that progress has been made in computations of ductile failure, but improvements still can be made, hence the third Sandia Fracture Challenge (SFC3), the subject of this Special Volume. The most recent installment of SFC is building on previous successes and tackling the difficult problem of fracture in an AM 316L stainless steel structure.

42 ENGINEERING↗

Implications of Grain-Scale Mineralogical Heterogeneity for Radionuclide Transport in Fractured Media

The geological disposal of nuclear waste is based on the multi-barrier concept, comprising various engineered and natural barriers, to confine the radioactive waste and isolate it from the biosphere. Some of the planned repositories for high-level nuclear waste will be hosted in fractured crystalline rock formations. The potential of these formations to act as natural transport barriers is related to two coupled processes: diffusion into the rock matrix and sorption onto the mineral surfaces available in the rock matrix. Different in situ and laboratory experiments have pointed out the ubiquitous heterogeneous nature of the rock matrix: mineral surfaces and pore space are distributed in complex microstructures and their distribution is far from being homogeneous (as typically assumed by Darcy-scale coarse reactive transport models). We use a synthetically generated fracture–matrix system to assess the implications of grain-scale physical and mineralogical heterogeneity on cesium transport and retention. The resulting grain-scale reactive transport model is solved using high-performance computing technologies, and the results are compared with those derived from two alternative models, denoted as upscaled models, where mineral abundance is averaged over the matrix volume. In the grain-scale model, the penetration of cesium into the matrix is faster and the penetration front is uneven and finger-shaped. The analysis of the cesium breakthrough curves computed at two different points in the fracture shows that the upscaled models provide later first-arrival time estimates compared to the grain-scale model. The breakthrough curves computed with the three models converge at late times. These results suggest that spatially averaged upscaled parameters of sorption site distribution can be used to predict the late-time behavior of breakthrough curves but could be inadequate to simulate the early behavior.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Effects of the Electron Beam Welding Process on the Microstructure, Tensile, Fatigue and Fracture Properties of Nickel Alloy Nimonic 80A

The purpose of this study was to evaluate rotary bending high-cycle fatigue properties and crack growth of Nimonic 80A-based metal and electron beam-welded joints. All the tests were performed at room temperature. Fracture surfaces under high-cycle fatigue and fatigue crack growth were observed by scanning electron microscopy. Microstructure, hardness and tensile properties were also evaluated in order to understand the effects on the fatigue results obtained. It was found that the tensile properties, hardness and high-cycle fatigue properties of the welded joint are lower than the base metal. The fracture surface of the high-cycle fatigue shows that fatigue crack initiated from the surface under the high stress amplitude and from the subsurface under the low stress amplitude. The effect of the welding process on the statistical fatigue data was studied with a special focus on probabilistic life prediction and probabilistic lifetime limits. The fatigue crack growth rate versus stress intensity factor range data were obtained from the fatigue crack growth tests. From the results, it was evident that the fatigue crack growth rates of the welded are higher than the base metal. The mechanisms and fracture modes of fatigue crack growth of welded specimens were found to be related to the stress intensity factor range ΔK. In addition, the effective fatigue crack propagation thresholds and mismatch of welded joints were described and discussed.

36 MATERIALS SCIENCE↗

The combined effects of hydrogen and aging condition on the deformation and fracture behavior of a precipitation-hardened nickel-base superalloy

The effect of hydrogen (H) on the deformation behavior of Monel K-500 in various isothermal heat treatment conditions (non-aged, under-aged, peak-aged, and over-aged) was assessed via uniaxial mechanical testing. H-charged and non-charged specimens were strained to failure to facilitate a comparison of ductility, fracture surface morphology, strength, and work hardening behavior. For all examined heat treatment conditions, H charging leads to a significant reduction in ductility, which is accompanied by a consistent change in fracture surface morphology from ductile microvoid coalescence to brittle intergranular fracture. While H charging led to a systematic enhancement in the yield strength of all heat treatments, the three age-hardened conditions exhibited a more than 2-fold increase relative to the non-aged heat treatment. This suggests that H modifies the dislocation–precipitate interactions, which also manifest themselves through changes in work hardening metrics related to the dislocation storage and recovery rates. In particular, the H-charged peak-aged specimen exhibited a significant increase in initial hardening (dislocation storage) rate relative to the H-charged under-aged specimen. Transmission electron microscopy of these samples revealed the onset of widespread dislocation looping in the H-charged peak-aged sample, in addition to the planar slip bands characteristic of the non-charged condition. This result suggests that hydrogen induces the particle shearing-to-looping transition at smaller particle sizes. Possible mechanistic explanations for this observed behavior are presented.

36 MATERIALS SCIENCE↗

On the strength and fracture toughness of an additive manufactured CrCoNi medium-entropy alloy

An additively manufactured (nominally equiatomic) CrCoNi alloy was processed by laser powder bed fusion (LPBF). At ambient temperatures (298 K), this medium-entropy alloy displayed a yield strength, σ y of ~691 ± 9 MPa, and an ultimate tensile strength, σ u of ~926 ± 15.2 MPa; at cryogenic temperatures (77 K), yield and tensile strengths increased respectively to σ y ~ 944 ± 6 MPa and σ u ~ 1382 ± 11 MPa. These strength levels are 57 and 44% higher than that of the wrought alloy, due to strengthening from the solidification cellular structures intertwined with dislocations in the LPBF CrCoNi. The crack-initiation fracture toughness, K JIc was measured to be ~183.7 ± 28 MPa√m at 298 K; this value marginally decreased by ~4% to ~176 ± 11 MPa√m at 77 K. These K JIc values of the LPBF CrCoNi were 11% and 35% lower than the wrought CrCoNi alloy at 298 K and 77 K, respectively. Here, the resistance to crack growth of the LPBF CrCoNi from its hierarchical micro- and meso-structures was evaluated using nonlinear-elastic fracture mechanics by measuring R-curve behavior in the form of the J-integral as a function of crack extension. The specific features of the hierarchical microstructures at different length-scales provide a basis for the strengthening and toughening properties of this additively manufactured medium-entropy alloy. This correlation between the deformation and the hierarchical microstructures at different length-scales may provide future guidance for improving the fracture toughness properties of medium-entropy alloys.

36 MATERIALS SCIENCE↗

Atomistic modeling of tensile deformation and fracture of carbon fibers: Nanoscale stress redistribution, effect of local structural characteristics and nanovoids

The mechanisms of tensile deformation and fracture of carbon fibers are investigated in large-scale molecular dynamics simulations. The simulations performed for computational samples with different structural characteristics enable analysis of the structural sensitivity of the mechanical properties. It is found that the carbon fiber nanostructure can be represented by a fine (nanoscale) mixture of stiff and soft regions. Despite the nanoscale structural heterogeneity of the computational fibers, their fracture strain and tensile strength exceed the experimental values by almost an order of magnitude. Here, this observation is explained by the ability of the carbon fiber nanostructure, largely consisting of curved and folded layers of turbostratic carbon, to exhibit continuous nanoscale load redistribution during the deformation. The transient, fluid nature of a percolating cluster of stress-holding structural elements prevents stress localization at the early stage of the deformation and enables the fiber nanostructure to adapt to the increasing load. Overall, the results of the simulations suggest that the design of advanced carbon fiber manufacturing and post-processing methods focused on the elimination of critical mesoscopic structural defects (internal voids, foreign inclusions, and surface flaws) can yield a substantial increase in the strength and fracture strain of carbon fibers.

36 MATERIALS SCIENCE↗

A variational phase-field framework for thermal softening and dynamic ductile fracture

A variational phase field model for dynamic ductile fracture is presented. Here, the model is designed for elasto-viscoplastic materials subjected to rapid deformations in which the effects of heat generation and material softening are dominant. The variational framework allows for the consistent inclusion of plastic dissipation in the heat equation as well as thermal softening. It employs a coalescence function to degrade fracture energy during regimes of high plastic flow. A variationally consistent form of the Johnson–Cook model is developed for use with the framework. Results from various benchmark problems in dynamic ductile fracture are presented to demonstrate capabilities. In particular, the ability of the model to regularize shear band formation and subsequent damage evolution in two- and three-dimensional problems is demonstrated. Importantly, these phenomena are naturally captured through the underlying physics without the need for phenomenological criteria such as stability thresholds for the onset of shear band formation.

42 ENGINEERING↗

In-situ imaging of flexure-induced fracture in tape-laminate composites using high-resolution X-ray computed tomography

Here, this study presents development of an in-situ flexure test for imaging progressive inter- and intralaminar fracture in tape-laminate composites using X-ray computed tomography (CT). The intent of this test is to provide detailed experimental observations of ply-level damage that can be used to validate existing, and develop new, progressive damage analysis (PDA) tools. The test consists of a vertically mounted specimen which is flexed using two eccentric compressive loads using an in-situ uniaxial load stage. The flexure specimen contains a starter notch mid-span which promotes initiation of composite failure within the X-ray field of view. Specimens with two different laminate stacking sequences were tested, imaged, analyzed. For a quasi-isotropic laminate with large angle changes between adjacent plies, there was near simultaneous growth of transverse cracks and delaminations below the midplane of the laminate. For a laminate with small angles between adjacent plies, there was extensive formation of transverse crack networks which penetrated the laminate thickness without delamination growth. In addition to imaging fracture, the X-ray CT data from both specimen types were used to quantify the variability in thickness and analyze the local orientation of individual plies. Overall, the proposed test and the image-data analysis methodology provided an important insight into the fracture processes in tape laminates and highlights the inherent ply-level geometrical variabilities that should be accounted for in PDA simulations.

42 ENGINEERING↗

Poromechanical cohesive interface element with combined Mode I-II cohesive zone elastoplasticity for simulating fracture in fluid-saturated porous media

A combined Mode I-II cohesive zone (CZ) elasto-plastic constitutive model, and a two-dimensional (2D) cohesive interface element (CIE) are formulated and implemented at small strain within an ABAQUS User Element (UEL) for simulating 2D crack nucleation and propagation in fluid-saturated porous media. Here, the CZ model mitigates problems of convergence for the global Newton-Raphson solver within ABAQUS, which when combined with a viscous stabilization procedure allows for simulation of post-peak response under load control for coupled poromechanical finite element analysis, such as concrete gravity dam stability analysis. Verification examples are presented, along with a more complex ambient limestone-concrete wedge fracture experiment, water-pressurized concrete wedge experiment, and concrete gravity dam stability analyses. A calibration procedure for estimating the CZ parameters is demonstrated with the limestone-concrete wedge fracture process. For the water-pressurized concrete wedge fracture experiment it is shown that the inherent time-dependence of the poromechanical CIE analysis provides a good match with experimental force versus displacement results at various crack mouth opening rates, yet misses the pore water pressure evolution ahead of the crack tip propagation. This is likely a result of the concrete being partially-saturated in the experiment, whereas the finite element analysis assumes fully water saturated concrete. For the concrete gravity dam analysis, it is shown that base crack opening and associated water uplift pressure leads to a reduced Factor of Safety, which is confirmed by separate analytical calculations.

97 MATHEMATICS AND COMPUTING↗

Screening improved recovery methods in tight-oil formations by injecting and producing through fractures

Tight-oil reservoirs exhibit two characteristic behaviors that limit potential for improved recovery: (i) limited fluid movement from the unfractured matrix limits the potential to sweep additional oil towards the production wells with injection of a displacing fluid, (ii) the wettability of these reservoirs tend to be oil-wet which holds oil in relatively smaller pores, as compared to gas or water, making it difficult to mobilize that trapped oil. Primary oil production from these formations can decline to half of initial rates in the first year due to low permeability in the unfractured reservoir matrix, and a large amount of unrecovered oil remains in smaller pores that is not able to move out. To address these challenges, this study investigates the potential to improve oil recovery from tight-oil reservoirs by a method that includes injecting a fluid into the fractured reservoir to produce hydrocarbons from adjacent fractures intersecting the same wellbore. The injection and production through fractures has a potential economic advantage over huff-n-puff scheme such that there is no lag-period between injection and production. We compare the estimates of incremental oil recovery obtained by injecting water, immiscible CO 2 , and surfactant over the oil recovered by primary depletion (with no injection). The results suggest that injection of fluids may not always improve recovery from tight-oil reservoirs, especially in a heterogeneous reservoir (having different rock types) whose average matrix permeability is lower than 0.01 mD. In cases where recovery is improved over primary depletion, water flood and surfactant flood perform almost equally well with no noticeable difference between the two. Although, surfactant results in favorable fluid properties to mobilize oil, the reason surfactant did not perform well in tight-oil formation studied here is because the unswept oil is inaccessible due to small permeability, and not because of high capillary pressure.

42 ENGINEERING↗

Tensile mechanical properties and fracture behaviors of nickel-based superalloy 718 in the presence of hydrogen

Hydrogen embrittlement of a nickel-based superalloy IN718 was investigated using slow strain rate tensile tests. Post-mortem observation of fractured samples was performed to explore hydrogen-assisted failure mechanisms of the alloy. The results reveal that hydrogen charging reduces yield strength, tensile strength, fracture strain and work hardening rate. With increasing current density, yield strength and tensile strength reduce linearly and fracture strain decreases exponentially. Furthermore, the crack initiation and propagation in hydrogen-charged region depends on the distribution of δ phase in the alloy. For needle-shaped δ phase within the grains, the nucleation of voids takes place at the intersections between dislocation slip bands and δ phase due to the hydrogen-enhanced localized plasticity (HELP)-assisted shear localization and possible hydrogen agglomeration. For δ phase along the grain boundaries, the impingement of slip bands and local hydrogen accumulation at γ-matrix/δ phase interfaces as well as hydrogen-enhanced decohesion (HEDE)-assisted decohesion lead to the void nucleation at the interfaces. Because of the decoration of δ phase at the grain boundaries, hydrogen-assisted cracking preferentially propagates along the grain boundaries. It is therefore suggested that the synergistic interplay of HELP mechanism and HEDE mechanism can be used to explain the embrittlement of the alloy.

08 HYDROGEN↗

Implementation and experimental validation of nonlocal damage in a large-strain elasto-viscoplastic FFT-based framework for predicting ductile fracture in 3D polycrystalline materials

Ductile materials, such as metal alloys, can undergo substantial deformation before failure. Additionally, these materials are usually of polycrystalline composition and exhibit strongly anisotropic behavior at small length scales. Previously developed fast Fourier transform (FFT)-based models can model ductile fracture of isotropic materials or the elastic–plastic behavior of anisotropic polycrystalline materials; however, there remains a need to couple both capabilities. This work extends a large-strain FFT-based crystal plasticity model to simulate ductile fracture of polycrystalline materials. In this work, a triaxiality-based continuum damage mechanics (CDM) formulation is incorporated into a large-strain elasto-viscoplastic FFT (LS-EVPFFT) framework. The CDM formulation is augmented with an integral-based nonlocal regularization approach that correctly handles gas-phase material necessary to model unconstrained surfaces. To validate the damage-enabled LS-EVPFFT framework, mesoscale copper tensile coupons were machined using microwire electrical discharge machining and experimentally characterized using electron backscatter diffraction. In-situ optical digital image correlation was performed during uniaxial testing to provide a side-by-side comparison of the experimental and computational strain fields and stress–strain responses. The damage-enabled LS-EVPFFT framework can simulate the complete macroscopic stress–strain response of ductile polycrystals to failure. The model reproduces necking behavior that qualitatively agrees with experimental observations. By leveraging the relatively low computational cost of the damage-enabled LS-EVPFFT framework, the framework presented here allows the ductile fracture response of 3D polycrystalline materials to be tractably predicted.

36 MATERIALS SCIENCE↗

Weight function procedure for reduced order fracture analysis of arbitrary flaws in cylindrical pressure vessels

Fracture mechanics calculations using the finite element method can be computationally expensive, which makes them challenging to use in engineering evaluations of flaws in pressure vessels. The weight function (WF) technique is a reduced order fracture modeling method that is widely used to greatly reduce these computational costs. Its computational efficiency is essential for use in probabilistic fracture mechanics evaluations of embrittled nuclear reactor pressure vessels (RPVs), due to the large number of sampled flaws that must be evaluated. Although the WF technique is general, it is typically used only for axis-aligned flaws. Recent discoveries of off-axis flaws in operating nuclear reactors have necessitated detailed simulations of such flaws and the interactions between neighboring flaws. This paper presents a generalized WF approach applicable for analyzing arbitrary flaw geometries in thick-walled cylindrical vessels, including surface-breaking and subsurface flaws, which can either axis-aligned or off-axis, and can account for interactions with other flaws. Herein, this approach is demonstrated on representative simulations of multiple flaw geometries in a RPV subjected to transient loading conditions. In all cases, the WF technique gives good comparison with benchmark results from direct simulation with greatly reduced computational effort.

42 ENGINEERING↗

Risk Factors for Pelvic Insufficiency Fractures in Locally Advanced Cervical Cancer Following Intensity Modulated Radiation Therapy

Purpose: To investigate the incidence of and risk factors for pelvic insufficiency fracture (PIF) after definitive chemoradiation therapy for locally advanced cervical cancer (LACC). Methods and Materials: We analyzed 101 patients with LACC treated from 2008-2014. Patients received weekly cisplatin and underwent external beam radiation therapy with 45 Gy in 25 fractions (node-negative patients) or 50 Gy in 25 fractions with a simultaneous integrated boost of 60 Gy in 30 fractions (node-positive patients). Pulsed dose rate magnetic resonance imaging guided adaptive brachytherapy was given in addition. Follow-up magnetic resonance imaging was performed routinely at 3 and 12 months after the end of treatment or based on clinical indication. PIF was defined as a fracture line with or without sclerotic changes in the pelvic bones. D{sub 50%} and V{sub 55Gy} were calculated for the os sacrum and jointly for the os ileum and pubis. Patient- and treatment-related factors including dose were analyzed for correlation with PIF. Results: The median follow-up period was 25 months. The median age was 50 years. In 20 patients (20%), a median of 2 PIFs (range, 1-3 PIFs) were diagnosed; half were asymptomatic. The majority of the fractures were located in the sacrum (77%). Age was a significant risk factor (P<.001), and the incidence of PIF was 4% and 37% in patients aged ≤50 years and patients aged >50 years, respectively. Sacrum D{sub 50%} was a significant risk factor in patients aged >50 years (P=.04), whereas V{sub 55Gy} of the sacrum and V{sub 55Gy} of the pelvic bones were insignificant (P=.33 and P=.18, respectively). A dose-effect curve for sacrum D{sub 50%} in patients aged >50 years showed that reduction of sacrum D{sub 50%} from 40 Gy{sub EQD2} to 35 Gy{sub EQD2} reduces PIF risk from 45% to 22%. Conclusions: PIF is common after treatment of LACC and is mainly seen in patients aged >50 years. Our data indicate that PIFs are not related to lymph node boosts but rather to dose and volume associated with irradiation of the elective pelvic target. Reducing the prescribed elective dose from 50 to 45 Gy may reduce the risk of PIF considerably.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Coupling geophysical, geological, geochemical and mineralogical assessments to examine preferential contaminant transport pathways in interbedded fractured bedrock

This study shows that a multi-faceted approach, combining borehole geophysical logging and surface seismic P-wave first-arrival tomography with confirmatory coring, well installation, and chemical and mineralogical analysis, is effective for identifying difficult-to-locate preferential contaminant transport pathways in deeper fractured bedrock. Seismic tomography detected porous 10–20 m wide elongated fractured conduits that allow acidic groundwater contaminated with uranium (U) and nitrate (NO 3 - ) to migrate within interbedded shale-limestone bedrock over 1000 m from a former disposal facility (S-3 Ponds site) located at the DOE Y-12 National Security Complex in Tennessee (USA). Conventional drilling techniques would easily miss these conduits because they are oriented parallel with fractured bedding planes. Synchrotron analysis of aquifer solids revealed that > 95 % of the U is hexavalent (U VI ). This uranyl (UO 2 2+ ) species is coordinated with carbonate, iron oxide, silicate and phosphate minerals within cemented white to yellow precipitates, which contain U concentrations as high as ∼21.6 % by weight fraction. Identifying the presence of these mineral phases, enables a further understanding of the potential effectiveness of remediation actions. The combination of methodologies presented here can also be applied to other explorations, such as the detection of water supply.

Fate and transport↗

Effects of the addition of boron nitride nanoplate on the fracture toughness, flexural strength, and Weibull Distribution of hydroxyapatite composites prepared by spark plasma sintering

Hydroxyapatite (HA) has inherently low fracture toughness and low flexural strength, thus limiting it from wide scale application as an implant material in the biomedical field. To increase the fracture toughness and flexural strength, HA composites were fabricated by adding boron nitride nanoplatelets (BNNP) as reinforcement. Spark plasma sintering was utilized to achieve fine grain structure. The addition of BNNP facilitated grain size refinement. The BNNP reinforced HA composites exhibited increased fracture toughness (2.3 MPa m 1/2 ) and flexural strength (79.79 MPa) of HA over previous published values (1.0 MPa m 1/2 ). Despite that the Weibull Distribution indicated a sacrifice in mechanical reliability, all the composites fabricated in this study showed a low probability of failure and a factor of safety (~ 5.6) that is consistent with that of human bones (~ 6). In addition, the current study provides an approach to statistically design sintering parameters and mechanical loading for fabrication of ceramics.

36 MATERIALS SCIENCE↗