Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix”

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,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Rapid isotopic analysis of uranium microparticles via SP-ICP-TOF-MS

Inductively coupled plasma – time-of-flight – mass spectrometry (ICP-TOF-MS) was employed for the isotopic analysis of uranium particles of varying 235 U enrichment levels. Here, a single particle (SP)-based introduction scheme was employed such that individual particles, in a suspension, were analyzed. The uranium oxide microparticles were comprised of depleted uranium (DU, 235 U/ 238 U of 0.0017316(14)), natural uranium (NU, 235 U/ 238 U of 0.0072614(39)), and low enriched uranium (LEU, 235 U/ 238 U of 0.051025(15)). The percent relative difference of the SP-ICP-TOF-MS measured isotopic ratios compared to the expected values for the DU, NU, and LEU particle populations were 8.75, 0.12, and 1.23 %, respectively. After characterization, the DU and NU particles were doped within a complex sample matrix (Arizona Test Dust) containing Fe, Ti, Al, and Si particles, among others. Then, the suspension was analyzed via SP-ICP-TOF-MS and the detected particles were classified as DU or NU based on their measured 235 U/ 238 U ratio. In the same analysis, the matrix particles (i.e., Al, Fe, and Ti) were detected, demonstrating the simultaneous nuclide detection provided by the measurement platform. The presented SP-ICP-TOF-MS methodology for uranium particle characterization proved to be a high throughput method for detecting and isotopically discerning uranium particles with varying enrichment levels, in a complex matrix.

Stanberry, Jordan S. [Oak Ridge National Laborator↗
DOE OSTI2021

Decellularized ECM derived from normal bone involved in the viability and chemo-sensitivity in multiple myeloma cells

Multiple myeloma (MM) is an incurable plasma cell malignancy. The progression of MM is closely related to the bone microenvironment. Bone matrix proteins are remodeled and manipulated to govern cancer growth during the process of MM. However the role of normal bone extracellular matrix in MM is still unclear. In this study the decellularized extracellular matrix derived from normal SD rats’ skulls (N-dECM) was prepared by decellularization technology. The CCK 8 assay and the dead-live cell kit assay were used to determine the viability of MM cells and the sensitivity to bortezomib. The Realtime PCR and Western blot assay were used to assay the mRNA and protein related to MM. Under the treatment of N-dECM, we found that the viability of MM cells was inhibited and the sensitivity of MM cells to bortezomib was increased. Additionally, the expression levels of APRIL and TACI, which participated in the progression of MM, were significantly decreased in MM cells. It suggested that N-dECM might inhibit the development of MM via APRIL-TACI axis, and our study may provide a novel and potential biomaterial for MM therapy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Netrin-like domain of sFRP4, a Wnt antagonist inhibits stemness, metastatic and invasive properties by specifically blocking MMP-2 in cancer stem cells from human glioma cell line U87MG

Rapid proliferation, high stemness potential, high invasiveness and apoptotic evasion are the distinctive hallmarks of glioma malignancy. The dysregulation of the Wnt/β-catenin pathway is the key factor regulating glioma malignancy. Wnt antagonist, secreted frizzled-related protein 4 (sFRP4), which has a prominent pro-apoptotic role in glioma stem cells, has two functional domains, the netrin-like domain (NLD), and cysteine-rich domain (CRD) both of which contribute to apoptotic properties of the whole protein. However, there are no reports elucidating the specific effects of individual domains of sFRP4 in inhibiting the invasive properties of glioma. This study explores the efficacy of the domains of sFRP4 in inhibiting the key hallmarks of glioblastoma such as invasion, metastasis, and stemness. We overexpressed sFRP4 and its domains in the glioblastoma cell line, U87MG cells and observed that both CRD and NLD domains played prominent roles in attenuating cancer stem cell properties. Significantly, we could demonstrate for the first time that both NLD and CRD domains negatively impacted the key driver of metastasis and migration, the matrix metalloproteinase-2 (MMP-2). Mechanistically, compared to CRD, NLD domain suppressed MMP-2 mediated invasion more effectively in glioma cells as observed in matrigel invasion assay and a function-blocking antibody assay. Fluorescent matrix degradation assay further revealed that NLD reduces matrix degradation. NLD also significantly disrupted fibronectin assembly and decreased cell adhesion in another glioma cell line LN229. In conclusion, the NLD peptide of sFRP4 could be a potent short peptide therapeutic candidate for targeting MMP-2-mediated invasion in the highly malignant glioblastoma multiforme.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Alpha Shape Analysis (ASA) Framework for Post- Clustering Property Determination in Atom Probe Tomographic Data

While application of clustering algorithms to atom probe tomography data have enabled quantification of solute clusters in terms of number density, size, and subcomposition there exist other properties (e.g., volume, surface area, and composition) that are better determined by defining an interface between the cluster and the surrounding matrix. The limitation in composition results from an ion selection step where the expected matrix ion types are omitted from the cluster search algorithm to enhance the contrast between the matrix and cluster and to reduce the complexity of the search. Previously, composition determination within solute clusters has utilized a secondary envelopment and erosion step on top of conventional methods such as maximum separation. In this work, we present a novel stochastic method that combines the particle identification fidelity of a conventional clustering algorithm with the analytical flexibility of mesh-based approaches through the generation of alpha shapes for each identified cluster. Furthermore, the corresponding mesh accounts for concave components of the clusters and determines the volume and surface area of the clusters; additionally, the mesh boundary is utilized to update the total composition according to the internal ions.

42 ENGINEERING↗
DOE OSTI2024

Assessing Gas-Phase Ion Reactivity of 50 Elements with NO and the Direct Application for 239 Pu in Complex Matrices Using ICP-MS/MS

Understanding the reactivity of metal cations with various reaction gases in ICP-MS/MS is important to determine the best gas to use for a given analyte/interference pair. In this study, nitric oxide (NO) was investigated as the reaction gas following previous experimental designs. The reactions with 50 elements were investigated to examine periodic trends in reactivity, validate theoretical modeling of reaction enthalpies as a method to screen reactant gases, and provide a baseline for potential in-line gas separation methods. ICP-MS/MS studies involving actinides are typically limited to Th, U and Pu, with analyses of Np and Am rarely reported in the literature. To date, only two previous methods have investigated the use of NO in ICP-MS/MS analyses. To showcase the utility of NO, a method was developed to measure 239 Pu in the presence of environmental matrix constituent and other actinides, like what could be expected from post-detonation debris, with no chemical separation prior to analysis. 239 Pu + was reacted to form 239 Pu 16 O + , eliminating interferences derived from the sample matrix by measuring 239 Pu + intensity at m/z = 255 ( 239 Pu 16 O + ). To validate NO for 238 U 1 H + interference removal in environmental matrices, standard reference materials were diluted to 1 mg per gram of solution and spiked to 0.05 pg/g of 239 Pu and 1 µg/g 238 U (Pu/U = 5 x 10 -8 ). Measured 239 Pu concentrations were within 6% of the spiked value. In conclusion, these results demonstrate that reliable 239 Pu measurements can be made at levels relevant to nuclear forensics without the need for extensive chemical matrix separation prior to analysis.

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

An Evaluation of Actinide Reactivity with CO 2 , O 2 , and O 2 /He Gases using Inductively Coupled Plasma Tandem Mass Spectrometry: Application to Simultaneous Measurement of 241 Am/ 241 Pu Ratios in Unseparated Complex Matrices

Accurate actinide measurements are critical within the field of nuclear science. Traditional methods for actinide quantification require time-consuming sample processing prior to analysis. There is a need for rapid analytical techniques that still maintain a high degree of accuracy. In this work, actinide reactivity was assessed for multiple oxygen-containing reaction gases using quadrupole inductively coupled plasma tandem mass spectrometry (Q-ICP-MS/MS) to evaluate actinide analysis in complex sample matrices without analyte-matrix separation. A novel method was developed to measure 241 Am/ 241 Pu in complex sample matrices using O 2 /He reaction gas with no matrix removal or analyte pre-concentration. This inline method reduces matrix-derived polyatomic interferences that complicate traditional ICP-MS analyses by mass-shifting to 241 Am 16 O + and 241 Pu 16 O 2 + , allowing Am and Pu to be mass separated for simultaneous analysis. While mass shifting is efficient, a small portion of Am + (<1.3%) and Pu + (<1.4%) react to from AmO 2 + and PuO + , respectively. Therefore, a mass balance approach was used, in combination with reactivity determined from 242 Pu and 243 Am standard solutions, to correct for residual 241 PuO + and 241 AmO 2 + . The method was validated by measuring 241 Am/ 241 Pu in Pu isotope standards CRM-136 and CRM-137 (separated in March/April 1970 and February 2022, respectively) in both neat solutions and complex matrices containing diluted soil (NIST SRM 2711a, >1000 µg·g -1 ). Method detection limits of 15.9 and 9.6 fg·g -1 were determined for 241 Am and 241 Pu, respectively, and 241 Am/ 241 Pu ratios were measured with accuracies within <3.5%. In conclusion, this work presents the first direct analysis of 241 Am/ 241 Pu in unseparated complex matrices, advancing capabilities for rapid actinide measurements.

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

Surfactants are Ineffective for Reducing Imbibition of Water-Based Fracturing Fluids in Deep Gas Reservoirs

Minimizing loss of injected hydraulic fracturing fluids into shale along fracture-matrix boundaries is desired because imbibed water restricts gas production and wastes valuable water resources. This problem has motivated the addition of surfactants into water-based hydraulic fracturing fluids in order to reduce the capillary driving force for imbibition. Here, we show that reduction in interfacial tension and wettability alteration has negligible ability to reduce imbibition in deep gas reservoirs. The effectiveness of altering capillary forces acting at the wetting front also depends on the injection pressure acting at the fracture-matrix boundary. The pressure at the interface between the fracture and the shale matrix is constrained between the reservoir pore pressure and formation pressure (rock fracture pressure, also known as breakdown pressure of the rock) and increases with depth to magnitudes that greatly exceed that of capillary pressures. Furthermore, the analyses presented here show that even maximum alteration of interfacial properties that result in strongly hydrophobic interactions between the fracturing fluid and reservoir rock is incapable of significantly reducing imbibition in deep reservoirs. Instead of using surfactants, this analysis points to decreases in wellbore shut-in pressures and shut-in times as practical options for reducing imbibition losses of water-based fluids.

03 NATURAL GAS↗
DOE OSTI2022

Impact of Concurrent Solubilization and Fines Migration on Fracture Aperture Growth in Shales during Acidized Brine Injection

Acidic hydraulic fracturing fluid chemically and physically alters shale rock fabric during injection and shut-in, creating a “reaction-altered zone” along the fracture faces. To better characterize the variable thickness and composition of this reaction-altered zone under advective flow, we take a coupled experimental and modeling approach. A fluidic cell, with six fiducial markers, is first fabricated to keep the rock sample in place during the core floods and to allow image alignment of acquired images. Then, we conduct a series of reactive core floods in a clay-rich siliceous Wolfcamp shale sample with 10 wt % carbonate, using a synthetic fracturing fluid under no confining stress and at room temperature. High-resolution computed tomography (CT) scans are periodically conducted to observe the spatial alteration of the fracture network. We then perform scanning electron microscopy (SEM-EDS) on the two orthogonal surfaces (fracture surface and freshly cut profile face) to generate high-resolution elemental maps that show the change in mineralogy, both with distance along a given flow path along a fracture surface and with depth from the fracture surface into the shale matrix. These results are contrasted against a two-dimensional (2D) advection-diffusion reaction model developed previously for batch reactions between shale and synthetic fracturing fluids. The model simulates the geochemical interaction occurring at the fracture/matrix interface and penetrating into the shale matrix during the reactive core flood. Both model and experimental results show that the acidic brine is neutralized during the core flood, corresponding to an increase in fracture aperture as a function of fluid volume injected with the greatest change near the inlet. SEM-EDS scans reveal significant dissolution of carbonates on the fracture surface without pyrite oxidation. The reactive transport model indicates that carbonate depletion into the shale interior should be observable, yet SEM-EDS shows no discernible loss of carbonate in the orthogonal profile face. The combination of these observations suggests an additional fracture evolution mechanism in the reactive system, i.e., fines migration. We show that fines migration enhances the access of fracturing fluid to the matrix resulting in a more pronounced fracture widening. Finally, we conclude that coupled mineral dissolution and fines migration govern fracture aperture growth during acidized brine injection. In this work, we effectively show the underlying risk of relying solely on models that do not include an important (transport) process that can alter the system significantly and propose a combined chemomechanical mechanism for fracture evolution appropriate for this shale mineralogy.

58 GEOSCIENCES↗
DOE OSTI2022

Fe Oxidation and Species Distribution at the Rock–Fluid Interface of Marcellus Shale Reacted with Hydraulic Fracturing Fluid

Hydraulic fracturing of shale reservoirs resulted in significant opportunity for increased oil and gas production in the United States. Rock-fluid interactions can cause mineral dissolution and precipitation reactions that lead to permeability changes in the shale matrix, which ultimately may affect transport pathways and hydrocarbon production. Understanding the distribution of secondary precipitates, such as barite and Fe(III) (hydro)oxides, and cation leaching at the rock-fluid interface is an important step to further investigate how these geochemical processes can change permeability and transport pathways. In this study, thin sections of the fracture-matrix interface were made from reacted Marcellus shale cores. The thin sections were characterized using synchrotron X-ray fluorescence imaging and synchrotron X-ray absorption spectroscopy. Fe species with different oxidation states were identified in the maps, together with barite and Ca distribution. The results show that ferrihydrite, as newly formed Fe(III)-bearing precipitates, aligned well with the border of the Ca (e.g., calcite) leaching region in the reaction front. Some Fe-containing clay also dissolved, but the dissolution region for the clay was not as deep as the calcite. Further, the reaction front is about three times deeper in the direction parallel to the shale bedding than that perpendicular to the bedding. The Ca leaching region can be an index for reaction front detection for Marcellus shale. Reactive transport modeling was conducted and the predicted Ca leaching boarder align well with ferrihydrite precipitation, consistent with the experimental observation. The carbonate mineral dissolution can be crucial to promote fluid access into the shale matrix. Together with our previous study on the shale reactive surface, this follow-up study showed similar Ca leaching region and Fe(III) precipitates distribution in the reaction front regardless of barite precipitation on the surface, indicating that the barite coatings on the surface may not pose a significant impact on reactive transport at the shale-fluid interface.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2025

Strain-Based Assessment of Shale Caprock during Cyclic Underground Hydrogen Storage

Successful large-scale underground hydrogen storage (UHS) in depleted gas reservoirs depends on the integrity of the overlying caprock to prevent hydrogen loss during cyclic injection and depletion. Prior studies on crushed Marcellus shale, a potential caprock, indicate that cyclic hydrogen injection and depletion induces microstructural changes, increasing porosity and permeability. However, the extent of these changes in intact shale remains unclear. This study presents a strain-based experimental approach to quantify volumetric strain evolution in intact Marcellus shale matrix under unconstrained stress conditions. A quadrant-shaped shale sample without visible fractures underwent eight hydrostatic pore-pressure cycles (injection to 1500 psi and depletion to 500 psi in 250 psi steps). Linear strain gauges measured strain in three orthogonal directions. Results indicate progressive plastic strain accumulation, leading to an ∼12% increase in matrix porosity after eight cycles, with an estimated 19% increase after 30 cycles. This porosity increase follows a logarithmic trend, suggesting a diminishing effect in later cycles. Additionally, permeability and diffusive mass flux are projected to rise by ∼70% over 30 cycles, enhancing hydrogen migration risk. The shale matrix also exhibited mechanical stiffening over successive cycles, limiting large-scale deformation but not preventing porosity enhancement. A new parameter, α, was introduced to characterize shale sensitivity to cyclic loading, aiding UHS caprock assessments. These findings underscore the necessity of incorporating cyclic loading effects in UHS site selection and operational strategies to ensure long-term storage integrity.

08 HYDROGEN↗
DOE OSTI2022

Quantum Perturbation Theory Using Tensor Cores and a Deep Neural Network

In this work, time-independent quantum response calculations are performed using Tensor cores. This is achieved by mapping density matrix perturbation theory onto the computational structure of a deep neural network. The main computational cost of each deep layer is dominated by tensor contractions, i.e., dense matrix–matrix multiplications, in mixed-precision arithmetics, which achieves close to peak performance. Quantum response calculations are demonstrated and analyzed using self-consistent charge density-functional tight-binding theory as well as coupled-perturbed Hartree–Fock theory. For linear response calculations, a novel parameter-free convergence criterion is presented that is well-suited for numerically noisy low-precision floating point operations and we demonstrate a peak performance of almost 200 Tflops using the Tensor cores of two Nvidia A100 GPUs.

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

A High-Efficiency Delayed Update Algorithm for Evaluating Slater Determinants in Quantum Monte Carlo

For quantum Monte Carlo simulations of molecular systems or supercells with thousands of electrons, matrix operations related to Slater determinants lead the computational cost. McDaniel et al. [J. Chem. Phys. 2017, 147, 174107] proposed a delayed update algorithm to increase computational efficiency by using matrix–matrix multiplication when updating the inverse matrices of Slater determinants. However, preparing intermediate matrices for applying the Sherman–Morrison–Woodbury formula remained a bottleneck. Here, in this work, we introduce an improved algorithm for CPUs and GPUs that (1) reduces this bottleneck by iteratively updating the intermediate matrices and (2) is efficient at any acceptance ratio, with no cost for rejected moves on CPUs and minimal cost on GPUs. We show the full scheme of integrating the delayed update algorithm into a single-electron move. The high efficiency of our algorithm is demonstrated on CPUs and GPUs for a 512 atom/6144 valence electron calculation, with 12× and 2× overall speed-up compared to traditional rank-1 update schemes in diffusion quantum Monte Carlo, respectively.

Luo, Ye [Argonne National Laboratory (ANL), Argonn↗
DOE OSTI2021

Effects of Strain Rate and Temperature on the Mechanical Properties of Simulated Silica Ionogels

Ionogels are hybrid materials formed by impregnating the pore space of a solid matrix with a conducting ionic liquid. By combining the properties of both component materials, ionogels can act as self-supporting electrolytes in Li batteries. In this study, molecular dynamics simulations are used to investigate the dependence of mechanical properties of silica ionogels on solid fraction, temperature, and pore width. Comparisons are made with corresponding aerogels. We find that the solid matrix fraction increases the moduli and strength of the ionogel. This varies nonlinearly with temperature and strain rate, according to the contribution of the viscous ionic liquid to resisting deformation. Additionally, owing to the temperature and strain sensitivity of the ionic liquid viscosity, the mechanical properties approach a linear mixing law at high temperature and low strain rates. The median pore width of the solid matrix plays a complex role, with its influence varying qualitatively with deformation mode. Narrower pores increase the relevant elastic modulus under shear and uniaxial compression but reduce the modulus obtained under uniaxial tension. Conversely, shear and tensile strength are increased by narrowing the pore width. All of these pore size effects become more pronounced as the silica fraction increases. Pore size effects, similar to the effects of temperature and strain rate, are linked to the ease of fluid redistribution within the pore space during deformation-induced changes in the geometry of the pores.

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

Structure–Thermodynamic Relationship of a Polysaccharide Gel (Alginate) as a Function of Water Content and Counterion Type (Na vs Ca)

Biofilms are the predominant mode of microbial life on Earth, and so a deep understanding of microbial communities-and their impacts on environmental processes-requires a firm understanding of biofilm properties. Because of the importance of biofilms to their microbial inhabitants, microbes have evolved different ways of engineering and reconfiguring the matrix of extracellular polymeric substances (EPS) that constitute the main non-living component of biofilms. This ability makes it difficult to distinguish between the biotic and abiotic origins of biofilm properties. An important route toward establishing this distinction has been the study of simplified models of the EPS matrix. This study builds on such efforts by using atomistic simulations to predict the nanoscale (≤10 nm scale) structure of a model EPS matrix and the sensitivity of this structure to interpolymer interactions and water content. To accomplish this, we use replica exchange molecular dynamics (REMD) simulations to generate all-atom configurations of ten 3.4 kDa alginate polymers at a range of water contents and Ca–Na ratios. Simulated systems are solvated with explicitly modeled water molecules, which allows us to capture the discrete structure of the hydrating water and to examine the thermodynamic stability of water in the gels as they are progressively dehydrated. Our primary findings are that (i) the structure of the hydrogels is highly sensitive to the identity of the charge-compensating cations, (ii) the thermodynamics of water within the gels (specific enthalpy and free energy) are, surprisingly, only weakly sensitive to cation identity, and (iii) predictions of the differential enthalpy and free energy of hydration include a short-ranged enthalpic term that promotes hydration and a longer-ranged (presumably entropic) term that promotes dehydration, where short and long ranges refer to distances shorter or longer than ~0.6 nm between alginate strands.

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

Exploiting a Shortcoming of Coupled-Cluster Theory: The Extent of Non-Hermiticity as a Diagnostic Indicator of Computational Accuracy

The fundamental non-Hermitian nature of the forms of the coupled-cluster (CC) theory widely used in quantum chemistry has usually been viewed as a negative, but the present paper shows how this can be used to an advantage. Specifically, the non-symmetric nature of the reduced one-particle density matrix (in the molecular orbital basis) is advocated as a diagnostic indicator of computational quality. In the limit of the full coupled-cluster theory [which is equivalent to full configuration interaction (FCI)], the electronic wave function and correlation energy are exact within a given one-particle basis set, and the symmetric character of the exact density matrix is recovered. The extent of the density matrix asymmetry is shown to provide a measure of “how difficult the problem is” (like the well-known T 1 diagnostic), but its variation with the level of theory also gives information about “how well this particular method works”, irrespective of the difficulty of the problem at hand. The proposed diagnostic is described and applied to a select group of small molecules, and an example of its overall utility for the practicing quantum chemist is illustrated through its application to the beryllium dimer (Be 2 ). Future application of this idea to excited states, open-shell systems, and symmetry-breaking problems and an extension of the method to the two-particle density are then proposed.

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

How Topological Polymer Loops on the Nanoparticle Surface Control the Mechanical Properties of Nanocomposites

Carbon black (CB) and silica (SiO 2 ) filled elastomers are known to be the most successful polymer nanocomposites (PNCs) in industry, where “bound rubber (BR)” (i.e., polymer chains that are physically or chemically adsorbed on the nanofiller surface) plays a critical role in their reinforcement. Here, we report a molecular-scale mechanism underlying the “BR-induced reinforcement” by integrating neutron scattering experiments and molecular dynamics simulations. Simplified non-cross-linked SiO 2 -filled polybutadiene (PB) and CB-filled PB reveal the critical role of topological polymer loops in the BR for the enhanced mechanical performance. The average loop size on the SiO 2 surface modified with a silane coupling agent is much smaller than that on the CB surface and the loops on the SiO 2 surface are densely formed, preventing interdigitation with the matrix chains. On the other hand, the larger, uncrowded loops formed on the CB surface facilitate the interdigitation with the matrix polymer chains even near the filler surface. In this way, a strong connectivity is established between a matrix and a nanofiller, resulting in an adhesive filler–polymer interface. Furthermore, our findings shed light on rich and complex physics and materials design problems in PNCs, where the topological polymer structure on the nanofiller surface directly controls the macroscopic mechanical properties.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Elastic Properties of Plasma-Exposed Tungsten Predicted by Molecular-Dynamics Simulations

We report results of systematic molecular-dynamics computations of the elastic properties of single-crystalline tungsten containing structural defects, voids and overpressurized He nanobubbles, related to plasma exposure of tungsten serving as a plasma- facing component (PFC) in nuclear fusion devices. Our computations reveal that the empty voids are centers of dilatation resulting in development of tensile stress in the tungsten matrix, whereas He-filled voids (nanobubbles) introduce compressive stress in the plasma-exposed tungsten. We find that the dependence of the elastic moduli of plasma-exposed tungsten, namely, the bulk, Young, and shear modulus, on its void fraction follows a universal exponential scaling relation. We also find that the elastic moduli of plasma-exposed tungsten soften substantially as a function of He content in the tungsten matrix, following an exponential scaling relation; this He-induced exponential softening is in addition to the softening caused in the matrix with increasing temperature. A systematic characterization of the dependence of the elastic moduli on the He bubble size reveals that He bubble growth affects significantly both the bulk modulus and the Poisson ratio of plasma-exposed tungsten, while its effect on the Young and shear moduli of the plasma-exposed material is weak. Our findings contribute directly to the development of a structure-properties database that is required for the predictive modeling of the dynamical response of PFCs in nuclear fusion devices.

36 MATERIALS SCIENCE↗
DOE OSTI2020

3D Printing of Liquid Crystal Elastomer Foams for Enhanced Energy Dissipation Under Mechanical Insult

Polymer foams are an essential class of lightweight materials used to protect assets against mechanical insults, such as shock and vibration. Two features are important to enhance their energy absorption characteristics: the foam structure and the matrix phase mechanical behavior. This study investigates novel approaches to control both of these features to enhance the energy absorption capability of flexible lattice foams. First, we consider 3D printing via digital light processing (DLP) as a method to control the foam mesostructure across a suite of periodic unit cells. Second, we introduce an additional energy dissipation mechanism in the solid matrix phase material by 3D printing the lattice foams with polydomain liquid crystal elastomer (LCE), which undergo a mechanically induced phase transition under large strains. This phase transition is associated with LC mesogen rotation and alignment and provides a second mechanism for mechanical energy dissipation in addition to the viscoelastic relaxation of the polymer network. Here, we contrast the 3D printed LCE lattices with conventional, thermomechanically near-equivalent elastomer lattice foams to quantify the energy-absorbing enhancement the LCE matrix phase provides. Under cyclic quasi-static uniaxial compression conditions, the LCE lattices show dramatically enhanced energy dissipation in uniaxial compression compared to the non-LCE equivalent foams printed with a commercially available photocurable elastomer resin. The lattice geometry also plays a prominent role in determining the energy dissipation ratio between the LCE and non-LCE foams. We show that when increasing the lattice connectivity, the foam deformation transitions from bending-dominated to stretching-dominated deformations, which generates higher axial strains in the struts and higher energy dissipation in the lattice foam, as stretching allows greater mesogen rotation than bending. The LCE foams demonstrate superior energy absorption during the repeated dynamic loading during drop testing compared with the non-LCE equivalent foams, demonstrating the potential of LCEs to enhance physical protection systems against mechanical impact.

3D printing↗