SEARCH · Search NASA

Results for “Triaxiality”

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.

SourceResearch connectionFrom yearThrough year

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

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Exploring pressure-dependent inelastic deformation and failure in bonded granular composites: An energetic materials perspective

In polymer-filled granular composites, damage may develop in mechanical loading prior to material failure. Damage mechanisms such as microcracking or plastic deformation in the binder phase can substantially alter the material’s mesostructure. For energetic materials, such as solid propellants and plastic bonded explosives, these mesostructural changes can have far reaching effects including degraded mechanical properties, potentially increased sensitivity to further insults, and changes in expected performance. Unfortunately, predicting damage is nontrivial due to the complex nature of these composites and the entangled interactions between inelastic mechanisms. In this work, we assess the current literature of experimental knowledge, focusing on the pressure-dependent shear response, and propose a simple simulation framework of bonded particles to study four limiting-case material formulations at both meso- and macro-scales. Further, to construct the four cases, we systematically vary the relative interfacial strength between the polymer binder and granular filler phase and also vary the polymer’s glass transition temperature relative to operating temperature which determines how much the binder can plastically deform. These simulations identify key trends in global mechanical response, such as the emergence of strain hardening or softening regimes with increasing pressure which qualitatively resemble experimental results. By quantifying the activation of different inelastic mechanisms, such as bonds breaking and plastically straining, we identify when each mechanism becomes relevant and provide insight into potential origins for changes in mechanical responses. The locations of broken bonds are also used to define larger, mesoscopic cracks to test various metrics of damage. We primarily focus on triaxial compression, but also test the opposite case of triaxial extension to highlight the impact of Lode angle on mechanical behavior.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Evolution of anisotropy with saturation and its implications for the elastoplastic responses of clay rocks

Many clay rocks have distinct bedding planes. Experimental studies have shown that their mechanical properties evolve with the degree of saturation (DOS), often with higher stiffness and strength after drying. For transversely isotropic rocks, the effects of saturation can differ between the bed-normal (BN) and bed-parallel (BP) directions, which gives rise to saturation-dependent stiffness and strength anisotropy. Accurate prediction of the mechanical behavior of clay rocks under partially saturated conditions requires numerical models that can capture the evolving elastic and plastic anisotropy with DOS. In this study, we present an anisotropy framework for coupled solid deformation-fluid flow in unsaturated elastoplastic media. Here we incorporate saturation-dependent strength anisotropy into an anisotropic modified Cam-Clay (MCC) model and consider the evolving anisotropy in both the elastic and plastic responses. The model was calibrated using experimental data from triaxial tests to demonstrate its capability in capturing strength anisotropy at various levels of saturation. Through numerical simulations, we demonstrate the role of evolving stiffness and strength anisotropy in the mechanical behavior of clay rocks. Plane strain simulations of triaxial compression tests were also conducted to demonstrate the impacts of material anisotropy and DOS on the mechanical and fluid flow responses.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Coupled Thermo-Hydrological-Mechanical-Chemical Behavior of Anisotropic Granite for Geologic Disposal of High-Level Radioactive Waste: A Core-Scale Laboratory Investigation

The coupled thermo-hydrological-mechanical-chemical (THMC) behavior of rock within an Excavation Damaged Zone (EDZ) is critical for the safety and long-term performance of a geological repository for high-level radioactive wastes. While many laboratory experiments have been conducted to investigate EDZ rocks, the flow and deformation characteristics resulting from anisotropic rock textures and microcrack distribution under triaxial loading and elevated temperatures remain poorly understood. Particularly, cracks at various scales serve as fast paths for fluid flow and solute transport and present as focal points of mechanical weakness, which complicate the coupled THMC processes in anisotropic EDZ rocks and challenge modeling predictions. Here, in this study, a series of core-scale experiments was conducted on three granite samples under repository-relevant conditions. These rock samples were obtained from the Grimsel Underground Research Laboratory (URL), featured by anisotropic minerals (represented by bedding layers) and microcrack distributions and coarse cm-scale grain sizes. During the experiments, samples were subjected to an elevated temperature at 90 °C and different triaxial loading conditions either by radial (normal to bedding layers) or axial (parallel to bedding layers) compaction. Water was injected into the samples, and the rock permeability evolutions and effluent water chemistry were monitored closely. For intact samples, thermal expansion of minerals at 90 °C resulted in a large, 75% irreversible permeability reduction and rock strengthening under radial compaction, while thermal impact was limited to a 15% permeability reduction under axial compaction. In contrast, for a sample containing open cracks, the growth of fractures during the experiment resulted in an abrupt permeability increase and fast failure at 90 °C. The effluent water chemistry indicates much more considerable mineral dissolution from large shear sliding than that in rocks dominated by mechanical compaction. These results helped better understand the coupled THMC processes in anisotropic rocks containing cracks, evaluate the behaviors of EDZ rocks, and predict the long-term evolution of EDZ for the performance of the repositories.

Coupled THMC processes↗

DOE OSTI2020

An inelastic homogenization framework for layered materials with planes of weakness

Many geologic materials have a composite structure, in which macroscopic mechanical behavior is determined by the properties, shape, and heterogeneous distribution of individual constituents. In particular, sedimentary rocks commonly exhibit a layered microstructure, with distinct bedding planes that can also form planes of weakness. In this work, we present a homogenization framework for modeling inelastic layered media. The proposed constitutive model allows for distinct micro-constitutive laws for each layer, explicit representation of layer distributions, as well as incorporation of imperfect bonding at the interface between adjacent layers. No a priori assumptions are needed regarding the specific constitutive models used for the layers and interfaces, providing significant modeling flexibility. The overall framework provides a simple and physically-motivated way of defining anisotropic material behavior as an emergent property of the layered microstructure. The model is calibrated using triaxial and true-triaxial experimental data to demonstrate its ability to describe anisotropic deformation and multiple modes of failure.

58 GEOSCIENCES↗

DOE OSTI2022

Experimental study of fracture slip due to stress perturbation in fractured geo-resources

Anthropogenic activities such as underground fluid injection/extraction, well stimulation, and dam impoundment will perturb the state-of-stress in the subsurface, which could activate critically stressed pre-existing fractures and faults. Fault and fracture activation can induce seismicity, change flow behavior, and damage asperities on the fracture surface to alter reservoir permeability. Here, in this study, triaxial shear tests were conducted on saw cut and rough fracture Barre granite specimens to better characterize how a decrease in normal stress and increase in shear stress will influence crystalline rock fracture hydromechanical properties. In a triaxial set up, confining pressure (CP) and differential stress (DS) were used to induce normal and shear stress on a 45° oriented fracture. Fracture slip and the corresponding stress relaxation was observed when the confining stress decreased. Fracture permeability increased when the fracture slipped. Moreover, new micro-cracks observed by x-ray images indicated the creation of new flow paths for fluid movement. The experimental outcomes of this study indicate that decrease in normal stress due to stress perturbations along the faults/fractures in the extensional stress regime can induce the slip along the fracture. This could have a significant influence on slope stability in quarries, stimulation and injection strategy in deep geo-resources, and potentially minimize the risk of induced seismicity. The observed stress dependency of hydraulic aperture in rough fracture during pre- and post-slip stages could provide key inputs for field-scale simulation of geo-resources.

58 GEOSCIENCES↗

DOE OSTI2021

Role of crystallographic orientation on intragranular void growth in polycrystalline FCC materials

In this study, we study the effect of crystallographic orientation and applied triaxiality on the growth of intragranular voids. Two 3D full-field micromechanics methods are used, the dilatational visco-plastic fast-Fourier transform (DVP-FFT) and the crystal plasticity Finite Elements (CP-FE), both of which incorporate a combination of crystalline plasticity and dilatational plasticity. We demonstrate with several select cases that predictions of void growth from both formulations agree qualitatively. With the more computationally efficient DVP-FFT, additional effects of polycrystalline microstructure and the influence of nearest neighborhood are investigated. Crystals bearing a single intracrystalline void are studied in three types of 3D microstructural environments: isolated single crystals, individual equal-sized grains within a regular polycrystal, and individual variable sized grains within a polycrystal with grains and voids randomly located. We show that loading type plays a significant role. In strain-rate controlled conditions, voids in the hardest [111]-crystals grow the fastest in time, whereas in stress-controlled conditions, voids in the softest [100]-crystal grow the fastest in time. The analysis reveals that on average void growth is slower for the same starting orientation in the polycrystal than in the single crystal. We find that at the highest triaxiality tested that the correlation between crystal orientation and void growth rate in the polycrystal strengthens, drawing closer to that seen in the isolated single crystals. These results and model can help guide the microstructural design of polycrystalline materials with high strength and damage-tolerance in high-rate deformation.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Coupled phase field damage and crystal plasticity analysis of intragranular fracture: The role of crystallographic orientation and voids

Damage evolution in engineering metal alloys at the grain scale exhibits significant microstructural heterogeneity and anisotropy. These heterogeneities create local hotspots for stress and strain localization, leading to void nucleation. Crystal orientation influences the active slip systems around voids, affecting lattice rotation and potentially forming discontinuities. At low triaxiality, voids may change shape due to lower stress, rotation, elongation, and coalescence. At high triaxiality, the correlation between crystal orientation and void growth rate becomes stronger, resembling the behavior observed in isolated single crystals. Therefore, understanding the effects of crystal orientation, heterogeneous strain, and defect evolution is crucial for single crystal fracture characterization. Here, in this work, a coupled phase-field damage (PFD) and crystal plasticity (CP) model is implemented within a finite element framework to analyze crystal deformation and failure. The CP method employs a dislocation density-based constitutive model, while intragranular failure is modeled using an anisotropic PFD method. The PFD model considers both the stored energy due to elastic stretching and the energy release due to defect formation and crack formation. A single crystal Al2219 with an intracrystalline spherical void is chosen to analyze fracture. The study finds that fracture propagation is strongly correlated with crystal orientations. This coupled CP-PFD model provides accurate failure prediction in crystalline materials by incorporating the effects of crystal orientations and existing voids. This study demonstrates how the local microstructure and defects influence plastic deformation and failure mechanisms in metal alloys.

Aluminum alloy↗

DOE OSTI2020

Void coalescence and ductile failure in IN718 investigated via high-energy synchrotron X-ray tomography and diffraction

Ductile failure through the growth and coalescence of voids is of particular relevance for many engineering materials. Yet, the lack of experimental measurements of the mechanical state of the material at an appropriate length scale has limited further understanding of the regime of ductile failure. In this study, local grain-scale experimental measurements are used to determine their relevance for describing and predicting ductile failure. Additive manufactured materials, due to the presence of inherent porosity, as well as the ability to tailor internal porosity, provides a promising avenue to study ductile failure. Selective laser melting is used to manufacture a specially designed specimen with two large, internal voids, in addition to the natural porosity, which is characteristic to the production process. The initial porosity and its evolution upon tensile loading are characterized via micro tomography. Several locations of void coalescence are captured in the sample, revealing activity of multiple modes of failure. Finite element simulations, with a simplified J2 plasticity model, instantiated with the initial void structure, is deployed, in which geometric localizations of stress concentration corresponded with experimentally observed sites of coalescence but was inadequate in capturing the appropriate failure mechanism. Furthermore, from the experimental results of the far-field high-energy diffraction microscopy, the heterogeneous micromechanical state is identified and tracked around voids due to the local grain interactions. These experiments determine narrow bands of low stress triaxiality, at the onset of failure, which highlight the path of coalescence through intervoid shearing. Diffraction spot spreading analysis aided characterization of intragranular plasticity and strain heterogeneity, which can be coupled with high triaxiality to form the conditions inducive for coalescence through intervoid necking. In this study, the experimentally determined grain-scale description of the micromechanical state provides a physical basis that can accurately capture both the sites and the mechanism of void coalescence at the onset of ductile failure.

36 MATERIALS SCIENCE↗

DOE OSTI2021

An improved specimen preparation method for marine shallow gas-bearing sand sediments and its validations

Instabilities of shallow gas-charged seabed are potential geological hazards in ocean engineering. In practice, the conventional field sampling techniques failed to obtain undisturbed gas-bearing sediments from the seabed for laboratory mechanical testing because of sensitive gas exsolution and escape from sediments. However, preparation of representative remoulded gas-charged specimens is a challenging issue, because it is rather difficult to quantitatively control the gas content and obtain uniform distribution of gas bubbles within the specimen. Given the above problems, this work proposes a reliable approach to reconstitute the high-saturation specimen of gas-charged sediments in the laboratory by an improved multifunction integrated triaxial apparatus (MITA). This apparatus is developed based on an advanced stress path triaxial system by introducing a temperature-controlled system and a wave-monitoring system. The temperature-controlled system is used to accurately mimic the in situ environments of sediments in the seabed. The wave-monitoring system is utilized to identify exsolution point of free gas and examine the disturbance of gas to specimens during gas exsolution. The detailed procedure of gassy specimen preparation is introduced. Then, the quality of prepared specimens using our improved apparatus is validated by the high-resolution micro-X-ray computed tomography (μCT) scanning test, from which bubble occurrence and size distribution within the gassy sand specimen can be obtained; and preliminary mechanical tests on gassy sand specimens with various initial saturation degrees are performed. The proposed specimen preparation procedure succeeds in proving the postulated occurrence state of gas bubbles in coarse-grained sediments and accurately controlling the gas content.

58 GEOSCIENCES↗

DOE OSTI2024

Impact of supercritical carbon dioxide on the frictional strength of and the transport through thin cracks in shale

Understanding the mechanical and transport behavior of thin (i.e. small aperture) cracks slipping under supercritical carbon dioxide (sc-CO 2 ) conditions is essential to evaluate the integrity of sealing formations with buoyant sc-CO 2 below and the success of waterless fracturing. The two major items of interest in this work are frictional strength and permeability change of the crack. We used a triaxial cell that permits in situ visualization to conduct and monitor slippage along the faces of narrow cracks subjected to triaxial stresses. Such cracks are analogs to small geological faults. We tested carbonate-rich, 1-inch diameter Wolfcamp shale samples that are saw cut 30° to vertical to create a thin crack. Friction coefficients ranged from about 0.6 to 0.8 consistent with expectations for brittle rocks. The sc-CO 2 generally did not alter friction coefficient over the time scale of experiments. From a transport perspective, saturating cracks with sc-CO 2 substantially decreased permeability of the crack by 26%–52%, while slip resulted in a variety of permeability responses. Overall, the combined impact of sc-CO 2 saturation and slip reduced fault permeability for all tests. Our observations support the notion that the sealing capacity of some caprocks improves when saturated with sc-CO 2 and that some slip of small fractures is not necessarily detrimental to caprock integrity.

58 GEOSCIENCES↗

DOE OSTI2023

Laboratory investigation of hydraulic fracturing in granitic rocks using active and passive seismic monitoring

SUMMARY Knowledge of the fracturing processes can be important for the optimization of pressurized fluid injection operations in the deep underground rock mass. Active and passive seismic monitoring techniques have been used in the field for tracking or mapping the propagating hydraulic fracture. Although both these monitoring techniques provide valuable information about the generated fracture network, it is difficult for either technique to comprehensibly identify the different processes associated with hydraulic fracturing. The combined active and passive monitoring has the potential for better characterization of the complex hydraulic fracturing phenomena. In this study, laboratory hydraulic fracturing experiments with combined active and passive seismic monitoring were conducted on true triaxially loaded Barre granite cubes with different fluid injection rates. The seismic inelastic fracturing was detected by 16 passive acoustic emission sensors, where 3678 and 2370 seismic source events were detected for the high and low injection rate experiments, respectively. For active monitoring, strong variations in the attributes of signals were observed which were transmitted through four source–receiver pairs, placed both perpendicular and parallel to the generated hydraulic fracture. Positive velocity changes were observed for active sensor pairs with ray paths passing through the generated hydraulic fracture indicating fluid permeation, whereas isolated dry deformation was characterized by a slight but permanent velocity decrease. Compared to velocity, the energy of the active signals was 1–2 orders of magnitude more sensitive to different hydraulic fracturing processes. However, the sensitivity and signatures of the active signal attributes were found to be dependent on the frequency range and direction of ray path with respect to the location of the generated fracture network. Using the coupled evaluation of the active and passive signals we were able to systematically identify various hydraulic fracturing processes including: (1) aseismic deformation, (2) fracture initiation and fluid permeation, (3) pressure build-up, (4) fracture propagation and (5) pressure release and leak-off. The results of this study showed that combining the respective advantages of active and passive seismic techniques and using both of them to monitor the failure processes can facilitate a more comprehensive understanding and better control of the hydraulic stimulations in subsurface operations.

Geochemistry & Geophysics↗

DOE OSTI2022

Constraining the shape of Milky Way satellites with distance gradients

We combine the Dark Energy Camera Legacy Survey DR8 photometry with Gaia photometry to study the 3D structure of Bootes I, Draco, Ursa Minor, Sextans, and Sculptor dwarf galaxies using blue horizontal branch (BHB) stars as distance indicators. We construct a new colour–absolute magnitude of BHB stars that we use to measure the distance gradients within the body of the dwarf galaxies. We detect a statistically significant non-zero gradient only in Sextans and Sculptor. Through modelling of the gradient and 2D density of the systems by triaxial Plummer models, we find that the distance gradients in both dwarf galaxies are inconsistent with prolate shape, but compatible with oblate or triaxial shapes. In order to explain the observed gradients, oblate models of Sextans and Sculptor need to have a significant intrinsic ellipticity larger than 0.47 for Sextans and 0.46 for Sculptor. Furthermore, the flattened oblate shape may imply a significant anisotropy in velocity distribution in order to be consistent with the lack of significant velocity gradients in these systems.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

First Observation of Multiphonon 𝛾-Vibrations in an Odd-Odd Nuclear System

Here, the identification of the first multiphonon 𝛾-vibrational bands in an odd-odd neutron-rich nucleus of the nuclear chart is presented. These high-spin structures of hard to access $^{104}_{41}\textrm{Nb}_{63}$, produced in fission, were studied by combining a spectrometer with isotopic resolution coupled to a 𝛾-ray tracking array and independently from high-fold 𝛾 coincidence measurements. Triaxial projected shell model calculations for the high-spin states are in good agreement with the measured observables for the yrast, one-phonon, and two-phonon 𝛾 bands. The possibility of an oblate shape of an isomeric state and coexistence of triaxial and oblate configurations is investigated from the decay of the 141 keV isomer. The present Letter illustrates the robustness of vibration excitations in the presence of an odd valence proton and neutron as well as the possibly coexisting shapes beyond the 𝑁 = 60 transitional region.

Wang, E. H. [Shandong Univ., Jinan (China); Vander↗

DOE OSTI2020

Effective field theory for deformed odd-mass nuclei

We develop an effective field theory (EFT) for deformed odd-mass nuclei. These are described as an axially symmetric core to which a nucleon is coupled. In the coordinate system fixed to the core the nucleon is subject to an axially symmetric potential. Power counting is based on the separation of scales between low-lying rotations and higher-lying states of the core. In leading order, core and nucleon are coupled by universal derivative terms. These comprise a covariant derivative and gauge potentials which account for Coriolis forces and relate to Berry-phase phenomena. At leading order, the EFT combines the particle-rotor and Nilsson models. Here, we work out the EFT up to next-to-leading order and illustrate the results in 239 Pu and 187 Os . At leading order, odd-mass nuclei with rotational bandheads that are close in energy and differ by one unit of angular momentum are triaxially deformed. For bandheads that are well separated in energy, triaxiality becomes a subleading effect. The EFT developed in this paper presents a model-independent approach to the particle-rotor system that is capable of systematic improvement.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Signature splitting of the g 7 / 2 [ 404 ] 7 / 2 + bands in Ba 131 and Ce 133

Excited states in Ba 131 and Ce 133 were studied using in-beam γ -ray spectroscopy through the Sn 122 ( C 13 , 4 n ) Ba 131 and Te 125 ( C 12 , 4 n ) Ce 133 reactions, respectively. A strongly coupled band, associated with the ν g 7 / 2 [ 404 ] 7 / 2 + configuration, was identified in Ba 131 and Ce 133 . It is the first time to observe the ν g 7 / 2 [ 404 ] 7 / 2 + bands in the N = 75 isotones. The signature partners exhibit considerable energy splitting in comparison with those in the π g 7 / 2 [ 404 ] 7 / 2 + bands in the odd- A Ta and Re isotopes. Here, extensive cranked shell model and quasiparticle-plus-triaxial-rotor model calculations reveal the origin of the signature splitting, which depends not only on the triaxiality, but also on the configuration mixing with nearby low- j orbitals.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Probing nuclear quadrupole deformation from correlation of elliptic flow and transverse momentum in heavy ion collisions

In heavy ion collisions, elliptic flow v2 and radial flow, characterized by eventwise average transverse momentum [p T ], are related to the shape and size of the overlap region, which are sensitive to the shape of colliding atomic nuclei. The Pearson correlation coefficient between v 2 and [p T ], ρ 2 , was found to be particularly sensitive to the quadrupole deformation parameter β that is traditionally measured in low energy experiments. Built on earlier insight that the prolate deformation β>0 reduces the ρ 2 in ultracentral collisions (UCC), we show that the prolate deformation β<0 enhances the value of ρ 2 . As β>0 and β<0 are the two extremes of triaxiality, the strength and sign of v 2 2 -[pT] correlation can be used to provide valuable information on the triaxiality of the nucleus. Our study provide further arguments for using the hydrodynamic flow as a precision tool to directly image the deformation of the atomic nuclei at extremely short timescale (<10 -24 s).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Skyrme–Hartree–Fock–Bogoliubov mass models on a 3D mesh: IIb. Fission properties of BSkG2

Large-scale models of nuclear structure are currently the only way to provide consistent datasets for the many properties of thousands of exotic nuclei that are required by nucleosynthesis simulations. In [W. Ryssens et al., Eur. Phys. J. A 58, 246 (2022)], we recently presented the new BSkG2 model based on an energy density functional of the Skyrme type. Relying on a flexible three-dimensional coordinate representation of the nucleus, the model takes into account both triaxial deformation and time-reversal symmetry breaking. BSkG2 achieves a state-of-the-art global description of nuclear ground state (g.s.) properties and reproduces in particular the known masses with a root-mean-square (rms) deviation of 678 keV. Moving beyond g.s. properties, the model also reproduces all empirical values for the primary and secondary barriers as well as isomer excitation energies of actinide nuclei with rms deviations below 500 keV, i.e. with unprecedented accuracy. Here we discuss in detail the extension of our framework to the calculation of the fission barriers of 45 actinide nuclei, including odd-mass and odd-odd systems. We focus in particular on the impact of symmetry breaking which is key to the accuracy of the model: we allow systematically for axial, reflection and time-reversal symmetry breaking. The effect of the latter on the fission properties of odd-mass and odd-odd nuclei is small, but we find that allowing for shapes with triaxial or octupole deformation, as well as shapes with both, is crucial to achieving this accuracy. Here, the numerical accuracy of our coordinate space approach, the variety of nuclear configurations explored and the simultaneous successful description of fission properties and known masses makes BSkG2 the tool of choice for the large-scale study of nuclear structure.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Atomistic Simulations of the Elastic Compression of Platinum Nanoparticles

Abstract The elastic behavior of nanoparticles depends strongly on particle shape, size, and crystallographic orientation. Many prior investigations have characterized the elastic modulus of nanoscale particles using experiments or simulations; however their reported values vary widely depending on the methods for measurement and calculation. To understand these discrepancies, we used classical molecular dynamics simulation to model the compression of platinum nanoparticles with two different polyhedral shapes and a range of sizes from 4 to 20 nm, loaded in two different crystal orientations. Multiple standard methods were used to calculate the elastic modulus from stress-vs-strain data for each nanoparticle. The magnitudes and particle-size dependence of the resulting moduli varied with calculation method and, even for larger nanoparticles where bulk-like behavior may be expected, the effective elastic modulus depended strongly on shape and orientation. Analysis of per-atom stress distributions indicated that the shape- and orientation-dependence arise due to stress triaxiality and inhomogeneity across the particle. When the effective elastic modulus was recalculated using a representative volume element in the center of a large nanoparticle, the elastic modulus had the expected value for each orientation and was shape independent. It is only for single-digit nanoparticles that meaningful differences emerged, where even the very center of the particle had a lower modulus due to the effect of the surface. These findings provide better understanding of the elastic properties of nanoparticles and disentangle geometric contributions (such as stress triaxiality and spatial inhomogeneity) from true changes in elastic properties of the nanoscale material.

36 MATERIALS SCIENCE↗