Search NASA⌕ Search

SEARCH · Search NASA

Results for “Percolation”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Dynamics on modular networks with heterogeneous correlations

We develop a new ensemble of modular random graphs in which degree-degree correlations can be different in each module, and the inter-module connections are defined by the joint degree-degree distribution of nodes for each pair of modules. We present an analytical approach that allows one to analyze several types of binary dynamics operating on such networks, and we illustrate our approach using bond percolation, site percolation, and the Watts threshold model. The new network ensemble generalizes existing models (e.g., the well-known configuration model and Lancichinetti-Fortunato-Radicchi networks) by allowing a heterogeneous distribution of degree-degree correlations across modules, which is important for the consideration of nonidentical interacting networks.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Metal-insulator transition in Na{sub x}WO{sub 3}: Photoemission spectromicroscopy study

We have investigated the validity of percolation model, which is quite often invoked to explain the metal-insulator transition in sodium tungsten bronzes, Na{sub x}WO{sub 3} by photoelectron spectromicroscopy. The spatially resolved direct spectromicroscopic probing on both the insulating and metallic phases of high quality single crystals of Na{sub x}WO{sub 3} reveals the absence of any microscopic inhomogeneities embedded in the system within the experimental limit. Neither any metallic domains in the insulating host nor any insulating domains in the metallic host have been found to support the validity of percolation model to explain the metal-insulator transition in Na{sub x}WO{sub 3}.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Growth-induced non-planar magnetic anisotropy in FeCoZr-CaF{sub 2} nanogranular films: Structural and magnetic characterization

The relation between nanoscale structure, local atomic order and magnetic properties of (FeCoZr){sub x}(CaF{sub 2}){sub 100−x} (29 ≤ x ≤ 73 at. %) granular films is studied as a function of metal/insulator fraction ratio. The films of a thickness of 1–6 μm were deposited on Al-foils and glass-ceramic substrates, by ion sputtering of targets of different metal/insulator contents. Structural characterization with X-ray and electron diffraction as well as transmission electron microscopy revealed that the films are composed of isolated nanocrystalline bcc α-FeCo(Zr) alloy and insulating fcc CaF{sub 2} matrix. They grow in a columnar structure, where elongated metallic nanograins are arranged on top of each other within the columns almost normal to the substrate surface. Mössbauer spectroscopy and magnetometry results indicate that their easy magnetization axes are oriented at an angle of 65°–74° to the surface in films with x between 46 and 74, above the electrical percolation threshold, which is attributed to the growth-induced shape anisotropy. Interatomic distances characteristic for metallic state of α-FeCo(Zr) nanograins were revealed by X-ray Absorption Spectroscopy. The results show a lack of surface oxidation of the alloy nanograins, so the growth-induced orientation of nanograins in the films cannot be attributed to this effect. The study is among the first to report a growth-induced non-planar magnetic anisotropy in metal/insulator granular films above the percolation threshold and to reveal the origin of it.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Magnetic interactions in γ-Fe{sub 2}O{sub 3}@SiO{sub 2} nanocomposites

Structural and magnetic properties of γ-Fe{sub 2}O{sub 3} nanoparticles of 4 nm diameter, dispersed into silica SiO{sub 2} matrix with a wide range value of volume fraction (0.05 to 1), were investigated. Produced γ-Fe{sub 2}O{sub 3}@SiO{sub 2} nanocomposites consist of an assembly of very small single domain magnetic object (<10 nm) with a random distribution of both inter-particle distance and direction of particle magnetic moment. We focused on the determination of a magnetic percolation threshold, defined as the magnetic particles concentration value above which the magnetic properties of isolated particles vanished in favor of a magnetic collective behavior induced by magnetic interactions. A percolation threshold value of 0.63 was obtained by a local probe technique such as {sup 57}Fe Mössbauer spectrometry and confirmed by global magnetic measurements through zero-field cooled, field-cooled, and ac susceptibility data. Below this threshold, dynamic ac magnetization measurements show a thermally activated Arrhenius dependence of the blocking temperature of superparamagnetic nanoparticles and above this critical value, one observes a slowing down of their dynamic properties, which lead toward the establishment of a spin-glass like state.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Form fluctuations of polymer loaded spherical microemulsions studied by neutron scattering and dielectric spectroscopy

We investigate the structure and shell dynamics of the droplet phase in water/AOT/octane microemulsions with polyethyleneglycol (M{sub W} = 1500) molecules loaded in the droplets. Size and polydispersity of the droplets is determined with small angle X-ray scattering and small angle neutron scattering experiments. Shell fluctuations are measured with neutron spin echo spectroscopy and related to the dynamic percolation seen in dielectric spectroscopy. Shell fluctuations are found to be well described by the bending modulus of the shell and the viscosities inside and outside the droplets. Addition of the polymer decreases the modulus for small droplets. For large droplets the opposite is found as percolation temperature shifts to higher values.

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

The effect of molecular mobility on electronic transport in carbon nanotube-polymer composites and networks

A multiscale modeling approach to the prediction of electrical conductivity in carbon nanotube (CNT)–polymer composite materials is developed, which takes into account thermally activated molecular mobility of the matrix and the CNTs. On molecular level, a tight-binding density functional theory and non-equilibrium Green's function method are used to calculate the static electron transmission function in the contact between two metallic carbon nanotubes that corresponds to electron transport at 0 K. For higher temperatures, the statistical distribution of effective contact resistances is considered that originates from thermal fluctuations of intermolecular distances caused by molecular mobility of carbon nanotube and the polymer matrix. Based on this distribution and using effective medium theory, the temperature dependence of macroscopic electrical resistivity for CNT-polymer composites and CNT mats is calculated. The predicted data indicate that the electrical conductivity of the CNT-polymer composites increases linearly with temperature above 50 K, which is in a quantitative agreement with the experiments. Our model predicts a slight nonlinearity in temperature dependence of electric conductivity at low temperatures for percolated composites with small CNT loading. The model also explains the effect of glass transition and other molecular relaxation processes in the polymer matrix on the composite electrical conductivity. The developed multiscale approach integrates the atomistic charge transport mechanisms in percolated CNT-polymer composites with the macroscopic response and thus enables direct comparison of the prediction with the measurements of macroscopic material properties.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Effect of epitaxial strain on small-polaron hopping conduction in Pr{sub 0.7}(Ca{sub 0.6}Sr{sub 0.4}){sub 0.3}MnO{sub 3} thin films

We investigated the epitaxial strain effect on the small-polaron hopping conduction properties in Pr{sub 0.7}(Ca,Sr){sub 0.3}MnO{sub 3} (PCSMO) films. An increase in the carrier localization, as evidenced by the enhancement of the small-polaron activation energy E{sub A} in the high temperature region, was obtained by increasing the epitaxial lattice strain in either the tensile or compressive strained film. Furthermore, it was found that the magnitude of E{sub A}, and thus the carrier localization, strongly depends on the sign of the lattice strain, which explains the diverse percolative transport behaviors in PCSMO films with different types of strains. Meanwhile, similar dependencies on the strain of the films were also obtained for the electron-phonon interaction, characterized by the calculated small-polaron coupling constant. Our results reveal that the type of lattice strain plays a crucial role in determining the degree of localization of charge carriers and the electron-phonon coupling strength, which is important for understanding the lattice strain-induced various novel properties related to phase separation and percolation behaviors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

The role of interfacial metal silicates on the magnetism in FeCo/SiO{sub 2} and Fe{sub 49%}Co{sub 49%}V{sub 2%}/SiO{sub 2} core/shell nanoparticles

We have investigated the role of spontaneously formed interfacial metal silicates on the magnetism of FeCo/SiO{sub 2} and Fe{sub 49%}Co{sub 49%}V{sub 2%}/SiO{sub 2} core/shell nanoparticles. Element specific x-ray absorption and photoelectron spectroscopy experiments have identified the characteristic spectral features of metallic iron and cobalt from within the nanoparticle core. In addition, metal silicates of iron, cobalt, and vanadium were found to have formed spontaneously at the interface between the nanoparticle core and silica shell. X-ray magnetic circular dichroism experiments indicated that the elemental magnetism was a result of metallic iron and cobalt with small components from the iron, cobalt, and vanadium silicates. Magnetometry experiments have shown that there was no exchange bias loop shift in the FeCo nanoparticles; however, exchange bias from antiferromagnetic vanadium oxide was measured in the V-doped nanoparticles. These results showed clearly that the interfacial metal silicates played a significant role in the magnetism of these core/shell nanoparticles, and that the vanadium percolated from the FeCo-cores into the SiO{sub 2}-based interfacial shell.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

On the thermal stability of physical vapor deposited oxide-hardened nanocrystalline gold thin films

We describe a correlation between electrical resistivity and grain size for PVD synthesized polycrystalline oxide-hardened metal-matrix thin films in oxide-dilute (<5 vol. % oxide phase) compositions. The correlation is based on the Mayadas-Shatzkes (M-S) electron scattering model, predictive of grain size evolution as a function of composition in the oxide-dilute regime for 2 μm thick Au-ZnO films. We describe a technique to investigate grain boundary (GB) mobility and the thermal stability of GBs based on in situ electrical resistivity measurements during annealing experiments, interpreted using a combination of the M-S model and the Michels et al. model describing solute drag stabilized grain growth kinetics. Using this technique, activation energy and pre-exponential Arrhenius parameter values of E{sub a} = 21.6 kJ/mol and A{sub o} = 2.3 × 10{sup −17} m{sup 2}/s for Au-1 vol. % ZnO and E{sub a} = 12.7 kJ/mol and A{sub o} = 3.1 × 10{sup −18} m{sup 2}/s for Au-2 vol. % ZnO were determined. In the oxide-dilute regime, the grain size reduction of the Au matrix yielded a maximum hardness of 2.6 GPa for 5 vol. % ZnO. A combined model including percolation behavior and grain refinement is presented that accurately describes the composition dependent change in electrical resistivity throughout the entire composition range for Au-ZnO thin films. The proposed correlations are supported by microstructural characterization using transmission electron microscopy and electron diffraction mapping for grain size determination.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Optimizing the fabrication process and interplay of device components of polymer solar cells using a field-based multiscale solar-cell algorithm

Both the device composition and fabrication process are well-known to crucially affect the power conversion efficiency of polymer solar cells. Major advances have recently been achieved through the development of novel device materials and inkjet printing technologies, which permit to improve their durability and performance considerably. In this work, we demonstrate the usefulness of a recently developed field-based multiscale solar-cell algorithm to investigate the influence of the material characteristics, like, e.g., electrode surfaces, polymer architectures, and impurities in the active layer, as well as post-production treatments, like, e.g., electric field alignment, on the photovoltaic performance of block-copolymer solar-cell devices. Our study reveals that a short exposition time of the polymer bulk heterojunction to the action of an external electric field can lead to a low photovoltaic performance due to an incomplete alignment process, leading to undulated or disrupted nanophases. With increasing exposition time, the nanophases align in direction to the electric field lines, resulting in an increase of the number of continuous percolation paths and, ultimately, in a reduction of the number of exciton and charge-carrier losses. Moreover, we conclude by modifying the interaction strengths between the electrode surfaces and active layer components that a too low or too high affinity of an electrode surface to one of the components can lead to defective contacts, causing a deterioration of the device performance. Finally, we infer from the study of block-copolymer nanoparticle systems that particle impurities can significantly affect the nanostructure of the polymer matrix and reduce the photovoltaic performance of the active layer. For a critical volume fraction and size of the nanoparticles, we observe a complete phase transformation of the polymer nanomorphology, leading to a drop of the internal quantum efficiency. For other particle-numbers and -sizes, we observe only a local perturbation of the nanostructure, diminishing the number of continuous percolation paths to the electrodes and, therefore, reducing the device performance. From these investigations, we conclude that our multiscale solar-cell algorithm is an effective approach to investigate the impact of device materials and post-production treatments on the photovoltaic performance of polymer solar cells.

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

The influence of hydrogen on the chemical, mechanical, optical/electronic, and electrical transport properties of amorphous hydrogenated boron carbide

Because of its high electrical resistivity, low dielectric constant (κ), high thermal neutron capture cross section, and robust chemical, thermal, and mechanical properties, amorphous hydrogenated boron carbide (a-B{sub x}C:H{sub y}) has garnered interest as a material for low-κ dielectric and solid-state neutron detection applications. Herein, we investigate the relationships between chemical structure (atomic concentration B, C, H, and O), physical/mechanical properties (density, porosity, hardness, and Young's modulus), electronic structure [band gap, Urbach energy (E{sub U}), and Tauc parameter (B{sup 1/2})], optical/dielectric properties (frequency-dependent dielectric constant), and electrical transport properties (resistivity and leakage current) through the analysis of a large series of a-B{sub x}C:H{sub y} thin films grown by plasma-enhanced chemical vapor deposition from ortho-carborane. The resulting films exhibit a wide range of properties including H concentration from 10% to 45%, density from 0.9 to 2.3 g/cm{sup 3}, Young's modulus from 10 to 340 GPa, band gap from 1.7 to 3.8 eV, Urbach energy from 0.1 to 0.7 eV, dielectric constant from 3.1 to 7.6, and electrical resistivity from 10{sup 10} to 10{sup 15} Ω cm. Hydrogen concentration is found to correlate directly with thin-film density, and both are used to map and explain the other material properties. Hardness and Young's modulus exhibit a direct power law relationship with density above ∼1.3 g/cm{sup 3} (or below ∼35% H), below which they plateau, providing evidence for a rigidity percolation threshold. An increase in band gap and decrease in dielectric constant with increasing H concentration are explained by a decrease in network connectivity as well as mass/electron density. An increase in disorder, as measured by the parameters E{sub U} and B{sup 1/2}, with increasing H concentration is explained by the release of strain in the network and associated decrease in structural disorder. All of these correlations in a-B{sub x}C:H{sub y} are found to be very similar to those observed in amorphous hydrogenated silicon (a-Si:H), which suggests parallels between the influence of hydrogenation on their material properties and possible avenues for optimization. Finally, an increase in electrical resistivity with increasing H at <35 at. % H concentration is explained, not by disorder as in a-Si:H, but rather by a lower rate of hopping associated with a lower density of sites, assuming a variable range hopping mechanism interpreted in the framework of percolation theory.

36 MATERIALS SCIENCE↗
DOE OSTI2015

The effect of filler aspect ratio on the electromagnetic properties of carbon-nanofibers reinforced composites

The effect of filler aspect ratio on the electromagnetic properties of epoxy-amine resin reinforced with carbon nanofibers is here investigated. A heat treatment at 2500 °C of carbon nanofibers seems to increase their aspect ratio with respect to as-received ones most likely due to a lowering of structural defects and the improvement of the graphene layers within the dixie cup conformation. These morphological differences revealed by Raman's spectroscopy and scanning electron microscopy analyses may be responsible for the different electrical properties of the resulting composites. The DC characterization of the nanofilled material highlights an higher electrical conductivity and a lower electrical percolation threshold for the heat-treated carbon nanofibers based composites. In fact, the electrical conductivity is about 0.107 S/m and 1.36 × 10{sup −3} S/m for the nanocomposites reinforced with heat-treated and as received fibers, respectively, at 1 wt. % of nanofiller loading, while the electrical percolation threshold falls in the range [0.05–0.32]wt. % for the first nanocomposites and above 0.64 wt. % for the latter. Moreover, also a different frequency response is observed since the critical frequency, which is indicative of the transition from a resistive to a capacitive-type behaviour, shifts forward of about one decade at the same filler loading. The experimental results are supported by theoretical and simulation studies focused on the role of the filler aspect ratio on the electrical properties of the nanocomposites.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Coarse-grained model of water diffusion and proton conductivity in hydrated polyelectrolyte membrane

Using dissipative particle dynamics (DPD), we simulate nanoscale segregation, water diffusion, and proton conductivity in hydrated sulfonated polystyrene (sPS). We employ a novel model [Lee et al. J. Chem. Theory Comput. 11(9), 4395-4403 (2015)] that incorporates protonation/deprotonation equilibria into DPD simulations. The polymer and water are modeled by coarse-grained beads interacting via short-range soft repulsion and smeared charge electrostatic potentials. The proton is introduced as a separate charged bead that forms dissociable Morse bonds with the base beads representing water and sulfonate anions. Morse bond formation and breakup artificially mimics the Grotthuss mechanism of proton hopping between the bases. The DPD model is parameterized by matching the proton mobility in bulk water, dissociation constant of benzenesulfonic acid, and liquid-liquid equilibrium of water-ethylbenzene solutions. The DPD simulations semi-quantitatively predict nanoscale segregation in the hydrated sPS into hydrophobic and hydrophilic subphases, water self-diffusion, and proton mobility. As the hydration level increases, the hydrophilic subphase exhibits a percolation transition from isolated water clusters to a 3D network. The analysis of hydrophilic subphase connectivity and water diffusion demonstrates the importance of the dynamic percolation effect of formation and breakup of temporary junctions between water clusters. The proposed DPD model qualitatively predicts the ratio of proton to water self-diffusion and its dependence on the hydration level that is in reasonable agreement with experiments.

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

Evaluating conducting network based transparent electrodes from geometrical considerations

Conducting nanowire networks have been developed as viable alternative to existing indium tin oxide based transparent electrode (TE). The nature of electrical conduction and process optimization for electrodes have gained much from the theoretical models based on percolation transport using Monte Carlo approach and applying Kirchhoff's law on individual junctions and loops. While most of the literature work pertaining to theoretical analysis is focussed on networks obtained from conducting rods (mostly considering only junction resistance), hardly any attention has been paid to those made using template based methods, wherein the structure of network is neither similar to network obtained from conducting rods nor similar to well periodic geometry. Here, we have attempted an analytical treatment based on geometrical arguments and applied image analysis on practical networks to gain deeper insight into conducting networked structure particularly in relation to sheet resistance and transmittance. Many literature examples reporting networks with straight or curvilinear wires with distributions in wire width and length have been analysed by treating the networks as two dimensional graphs and evaluating the sheet resistance based on wire density and wire width. The sheet resistance values from our analysis compare well with the experimental values. Our analysis on various examples has revealed that low sheet resistance is achieved with high wire density and compactness with straight rather than curvilinear wires and with narrower wire width distribution. Similarly, higher transmittance for given sheet resistance is possible with narrower wire width but of higher thickness, minimal curvilinearity, and maximum connectivity. For the purpose of evaluating active fraction of the network, the algorithm was made to distinguish and quantify current carrying backbone regions as against regions containing only dangling or isolated wires. The treatment can be helpful in predicting the properties of a network simply from image analysis and will be helpful in improvisation and comparison of various TEs and better understanding of electrical percolation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

The evolution of acidic and ionic aggregates in ionomers during microsecond simulations

We performed microsecond-long, atomistic molecular dynamics (MD) simulations on a series of precise poly(ethylene-co-acrylic acid) ionomers neutralized with lithium, with three different spacer lengths between acid groups on the ionomers and at two temperatures. Ionic aggregates form in these systems with a variety of shapes ranging from isolated aggregates to percolated aggregates. At the lower temperature of 423 K, the ionic aggregate morphologies to do not reach a steady-state distribution over the course of the simulations. At the higher temperature of 600 K, the aggregates are sufficiently mobile that they rearrange and reach steady state after hundreds of nanoseconds. For systems that are 100% neutralized with lithium, the ions form percolated aggregates that span the simulation box in three directions, for all three spacer lengths (9, 15, 21). In the partially neutralized systems, the morphology includes lithium ion aggregates that may also include some unneutralized acid groups, along with a coexisting population of acid group aggregates that form through hydrogen bonding. In the lithium ion aggregates, unneutralized acid groups tend to be found on the ends or sides of the aggregates.

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

Evolution of structure and dynamics of thermo-reversible nanoparticle gels—A combined XPCS and rheology study

A combined X-ray photon correlation spectroscopy and rheology report is carried out to capture the evolution of structure, fast particle-scale dynamics, and moduli (elastic and loss) at early times of gel formation near the fluid-gel boundary of a suspension of nanoparticles. The system is made up of moderately concentrated suspensions of octadecyl silica in decalin (φ = 0.2) undergoing thermoreversible gelation. Near the gel boundary, the rate of gel formation is very sensitive to changes in attraction strength. Yet, we find that at different attraction strengths, the system goes through identical intermediate states of microscopic and macroscopic behavior, even though the absolute time needed to form a gel varies by orders of magnitude. We identify a single dimensionless time parameter, t w /t g , where t w is the wait time following the quench and t g is the rheologically determined gel time, that captures the similarity in gel formation at a range of attraction strengths. Following a temperature quench below the gel boundary, the system is initially fluidlike and forms diffusive clusters (similar to 8.5 times the particle diameter). After a lag-time, t L , clusters aggregate to form a network like structure which is characterized by the onset of mechanical rigidity and a rapid growth in microscopic relaxation times. At t g , the Baxter parameter obtained from adhesive hard sphere fits of the structure factor attains a constant value corresponding to the theoretical percolation boundary, thus demonstrating that gelation is percolation driven.

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

Morphology and proton diffusion in a coarse-grained model of sulfonated poly(phenylenes)

A coarse-grained model previously used to simulate Nafion using dissipative particle dynamics (DPD) is modified to describe sulfonated Diels-Alder poly(phenylene) (SDAPP) polymers. The model includes a proton-hopping mechanism similar to the Grotthuss mechanism. The intramolecular parameters for SDAPP are derived from atomistic molecular dynamics (MD) simulation using the iterative Boltzmann inversion. The polymer radii of gyration, domain morphologies, and cluster distributions obtained from our DPD model are in good agreement with previous atomistic MD simulations. As found in the atomistic simulations, the DPD simulations predict that the SDAPP nanophase separates into hydrophobic polymer domains and hydrophilic domains that percolate through the system at sufficiently high sulfonation and hydration levels. Increasing sulfonation and/or hydration leads to larger proton and water diffusion constants, in agreement with experimental measurements in SDAPP. In the DPD simulations, the proton hopping (Grotthuss) mechanism becomes important as sulfonation and hydration increase, in qualitative agreement with experiment. The turning on of the hopping mechanism also roughly correlates with the point at which the DPD simulations exhibit clear percolated, hydrophilic domains, demonstrating the important effects of morphology on proton transport.

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

Constrained model calibration of grain structure dependent spall dynamics in shock-loaded tantalum

Here, we perform a gas gun experiment by shock loading tantalum samples of varying grain structures to assess the suitability of a numerical model for simulating spall behavior. The observed differences in spall strength, as well spallation and re-compression history, are not captured in uncalibrated hydrodynamic simulations. An optimization is performed on the Johnson spall model to determine the best parameters that fit the observed trends. Linear stability analysis is employed to motivate bounds on those parameters. Herein, optimized simulations agree well with the experimental results, reproducing pullback depth and recompression timescales across the different samples tested. Further, the observed pullback time of the single crystal sample was found to imply, via the stability analysis, a percolation threshold in good agreement with the theoretical value for a body centered cubic lattice. Therefore, the combined linear stability and percolation analysis shows promise and may be applied to other materials with diverse microstructures. Collectively, the findings demonstrate that the model is suitable for reproducing spall-induced free surface behavior across various microstructures, but also points to caution in using model coefficients for uncalibrated microstructures.

58 GEOSCIENCES↗