Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion behavior”

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 685 records · Page 38

Nuclear energy density functionals grounded in ab initio calculations

Here, we discuss the construction of a nuclear energy density functional (EDF) from ab initio computations and advocate the need for a methodical approach that is free from ad hoc assumptions. The equations of state (EoSs) of symmetric nuclear and pure neutron matter are computed using the chiral NNLO sat and the phenomenological AV4' + UIX c Hamiltonians as inputs to self-consistent Green's function (SCGF) and auxiliary field diffusion Monte Carlo (AFDMC) methods. We propose a convenient parametrization of the EoS as a function of the Fermi momentum and fit it on the SCGF and AFDMC calculations. We apply the ab initio based EDF to carry out an analysis of the binding energies and charge radii of different nuclei in the local density approximation. The NNLO sat -based EDF produces encouraging results, whereas the AV4' + UIX c -based one is farther from experiment. Possible explanations of these different behaviors are suggested, and the importance of gradient and spin-orbit terms is analyzed. Our paper paves the way for a practical and systematic way to merge ab initio nuclear theory and density functional theory, while shedding light on some critical aspects of this procedure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The Method of Finite Averages: A rigorous upscaling methodology for heterogeneous porous media

Rigorous upscaling techniques offer accurate and computationally-efficient strategies for modeling the average behaviors of multi-physical, multiscale phenomena in geological porous media. However, such techniques often rely on a variety of methodological assumptions that prohibit their rigorous application to practical systems (e.g., systems involving heterogeneous porous media, system-scale boundary conditions, and fine-scale dynamics that are not diffusion-dominant). In this work, we aim to formulate an upscaling methodology with few methodological assumptions to provide high levels of model generality and foster the utilization of rigorously-derived upscaled models in practice. In particular, we introduce the Method of Finite Averages (MoFA), a novel upscaling methodology for rigorously modeling heterogeneous porous media and system-scale boundary conditions. We then detail MoFA’s implementation for the advective–diffusive transport of a single species and compare the methodology with classic numerical techniques, as well as other rigorous upscaling techniques, to highlight MoFA’s unique combination of rigor and generality. We then validate the derived model while demonstrating its benefits in three numerical experiments. The results suggest that (1.) the applicability and a priori error guarantees of MoFA models do not directly depend on system geometry, (2.) a model’s applicability and error guarantees can be can arbitrarily expanded and reduced, respectively, with further computational expense, and (3.) downscaling with MoFA provides an efficient strategy for generating accurate pore-scale solutions from upscaled results. Ultimately, the results evidence that upscaled models can be rigorously derived for heterogeneous porous media systems and resolved in a fraction of the time it takes to perform the equivalent pore-scale simulations.

58 GEOSCIENCES↗

Explosive nucleation and growth of Pb islands on Ge(111) below room temperature via collective diffusion

Nucleation is a fundamental process in nature controlling phase transitions and pattern formation. Classical nucleation is based on the gradual aggregation of diffusing atoms after a critical cluster size is reached. Contrary to these expectations, we have observed in real time the formation of perfect Pb(111) islands with more than 10 5 atoms emerging out of the compressed wetting layer within a few seconds during growth of Pb on Ge(111) below room temperature. Pb deposited on Ge(111) in a temperature range between −30 °C and 10 °C exhibited explosive nucleation of height-selected islands upon reaching a critical coverage of 1.33 ± 0.07 ML with respect to Ge(111). Island nucleation was fueled by the 2% compression of the Pb wetting layer at the critical coverage. The island areas grew linearly with time, exhibiting collective diffusion, where thousands of atoms follow correlated, non-random walk diffusion. The total growth rate was higher with increasing temperature. Density functional theory simulations of the chemical potential and binding sites of the Pb/Ge(111) system give insight into the role of the compression of the wetting layer and the strain within the growing Pb layers, which explain this nonclassical behavior of the system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

The Role of Surface Hydrophobicity on the Structure and Dynamics of CO2 and CH4 Confined in Silica Nanopores

Advancing a portfolio of technologies that range from the storage of excess renewable natural gas for distributed use to the capture and storage of CO 2 in geological formation are essential for meeting our energy needs while responding to challenges associated with climate change. Delineating the surface interactions and the organization of these gases in nanoporous environments is one of the less explored approaches to ground advances in novel materials for gas storage or predict the fate of stored gases in subsurface environments. To this end, the molecular scale interactions underlying the organization and transport behavior of CO 2 and CH 4 molecules in silica nanopores need to be investigated. To probe the influence of hydrophobic surfaces, a series of classical molecular dynamics (MD) simulations are performed to investigate the structure and dynamics of CO 2 and CH 4 confined in OH-terminated and CH 3 -terminated silica pores with diameters of 2, 4, 6, 8, and 10 nm at 298 K and 10 MPa. Higher adsorption extents of CO 2 compared to CH 4 are noted on OH-terminated and CH 3 -terminated pores. The adsorbed extents increase with the pore diameter. Further, the interfacial CO 2 and CH 4 molecules reside closer to the surface of OH-terminated pores compared to CH 3 -terminated pores. The lower adsorption extents of CH 4 on OH-terminated and CH 3 -terminated pores result in higher diffusion coefficients compared to CO 2 molecules. The diffusivities of both gases in OH-terminated and CH 3 -terminated pores increase systematically with the pore diameter. The higher adsorption extents of CO 2 on OH-terminated and CH 3 -terminated pores are driven by higher van der Waals and electrostatic interactions with the pore surfaces, while CH 4 adsorption is mainly due to van der Waals interactions with the pore walls. These findings provide the interfacial chemical basis underlying the organization and transport behavior of pressurized CO 2 and CH 4 gases in confinement.

Mohammed, Sohaib↗

Regulating Li Nucleation and Growth Heterogeneities via Near-Surface Lithium-Ion Irrigation for Stable Anode-Less Lithium Metal Batteries

The inhomogeneous nucleation and growth of Li dendrite combined with the spontaneous side reactions with the electrolytes dramatically challenge the stability and safety of Li metal anode (LMA). Despite tremendous endeavors, current success relies on the use of significant excess of Li to compensate the loss of active Li during cycling. Herein, a near-surface Li + irrigation strategy is developed to regulate the inhomogeneous Li deposition behavior and suppress the consequent side reactions under limited Li excess condition. The conformal polypyrrole (PPy) coating layer on Cu surface via oxidative chemical vapor deposition technique can induce the migration of Li + to the interregional space between PPy and Cu, creating a near-surface Li + -rich region to smooth diffusion of ion flux and uniform the deposition. Moreover, as evidenced by multiscale characterizations including synchrotron high-energy X-ray diffraction scanning, a robust N-rich solid-electrolyte interface (SEI) is formed on the PPy skeleton to effectively suppress the undesired SEI formation/dissolution process. Strikingly, stable Li metal cycling performance under a high areal capacity of 10 mAh cm -2 at 2.0 mA cm -2 with merely 0.5 × Li excess is achieved. The findings not only resolve the long-standing poor LMA stability/safety issues, but also deepen the mechanism understanding of Li deposition process.

25 ENERGY STORAGE↗

Evolution of PTCDA-derived seeds prior to graphene nanoribbon growth on Ge(001)

The synthesis of graphene nanoribbons (GNRs) can be realized via CH4 chemical vapor deposition (CVD) on substrates such as Ge(001) that promote highly anisotropic growth kinetics. Small polycyclic aromatic hydrocarbons (PAHs) are first sublimed onto Ge near room temperature to form graphene-like seeds that subsequently initiate GNR growth with CH4 exposure at 1173 K. The behaviors of PAHs in the time in between their sublimation onto Ge and GNR growth are unclear. Here, we study an archetypical PAH, perylene-3,4,9,10-tetracarboxyl acid dianhydride (PTCDA), on Ge(001) using both scanning tunneling microscopy (STM) and density functional theory (DFT) to characterize PAH conformation, surface diffusivity, and clustering, versus temperature. PTCDA becomes mobile above 673 K (consistent with a DFT calculated diffusion barrier of 2.01 eV). The mobile PTCDA molecules meet, cluster, and fuse at Ge step edges. These clusters have an apparent height of 0.4 nm, similar to small graphene islands, and grow laterally with increasing temperature – reaching 1.2 – 2.1 nm in width at 1173 K, consistent with the seed size extrapolated from CVD experiments. These results provide a plausible picture for how PTCDA forms seeds for anisotropic GNR CVD and show that PAHs with reduced surface diffusivity and inter-molecular reactivity are needed to enable more monodisperse PAH-seeded GNR synthesis.

36 MATERIALS SCIENCE↗

A family of independent Variable Eddington Factor methods with efficient preconditioned iterative solvers

We present a family of discretizations for the Variable Eddington Factor (VEF) equations that have high-order accuracy on curved meshes and efficient preconditioned iterative solvers. The VEF discretizations are combined with the Discontinuous Galerkin transport discretization from to form effective high-order, linear transport methods. The VEF discretizations are derived by extending the unified analysis of Discontinuous Galerkin methods for elliptic problems presented by Arnold et al. to the VEF equations. This framework is used to define analogs of the interior penalty, second method of Bassi and Rebay, minimal dissipation local Discontinuous Galerkin, and continuous finite element methods. The analysis of subspace correction preconditioners, which use a continuous operator to iteratively precondition the discontinuous discretization, is extended to the case of the non-symmetric VEF system. Numerical results demonstrate that the VEF discretizations have arbitrary-order accuracy on curved meshes, preserve the thick diffusion limit, and are effective on a proxy problem from thermal radiative transfer in both outer transport iterations and inner preconditioned linear solver iterations. We demonstrate that the VEF solution converges to the S N transport solution as the mesh is refined on both problems with smooth and non-smooth behavior in angle. Parallel performance studies show that the interior penalty VEF discretization's linear solve weak scales out to 1024 processors and strong scales well on a single node. Particular attention is paid to the parallel performance of the VEF algorithm when used in combination with a parallel block Jacobi transport sweep.

97 MATHEMATICS AND COMPUTING↗

A high-dimensional neural network potential for molecular dynamics simulations of condensed phase nickel and phase transitions

A high-dimensional neural network interatomic potential was developed and used in molecular dynamics simulations of condensed phase Ni and Ni systems with liquid–solid phase coexistence. The reference data set was generated by sampling the potential energy surface over a broad temperature-pressure domain using ab initio MD simulations to train a unified potential. Excellent agreement was achieved between bulk face-centered cubic nickel thermal expansion simulations and relevant experimental data. The same potential also yields accurate structures and diffusivities in the liquid state. The phase transition between liquid and solid phases was simulated using the two-phase interface method. The predicted melting point temperature is within a few kelvins of the literature value. Here, the general methodology could be applied to describe crystals with much more complex phase behaviors.

74 ATOMIC AND MOLECULAR PHYSICS↗

Extended linear-in- T resistivity due to electron-phason scattering in moiré superlattices

Due to its incommensurate nature, moiré superlattices host not only acoustic phonons but also another type of soft collective modes called phasons. Here, we investigate the impact of electron-phason scattering on the transport properties of moiré systems. We show that the resistivity can scale linearly with temperature down to temperatures much lower than the Bloch-Grüneisen scale defined by electron kinematics on the Fermi surface. This result stems from the friction between layers, which transfers phason spectral weight to a broad diffusive low-energy peak in the mechanical response of the system. As a result, phason scattering becomes a very efficient channel for entropy production at low temperatures. We also consider the contributions of phasons to thermodynamic properties at low temperatures and find a “metallic-like” linear-in-T behavior for the specific heat, despite the fact that this behavior is due to mechanical and not electronic degrees of freedom. Finally, we discuss the implications of this finding to reports of linear-in-T resistivity in the phase diagram of twisted bilayer graphene.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Radionuclide Waste Disposal: Development of Multi-scale Experimental and Modeling Capabilities (Final Report)

The DOE EPSCoR Implementation project “Radioactive waste management: Development of multi-scale experimental and modeling capabilities” helped to develop a team of scientists and engineers from Clemson University, South Carolina State University, and the University of South Carolina to address the disposition of nuclear wastes and study the transport of radioisotopes from a waste repository in the near and far field. The project involved 20 faculty from the three institutions as well as 13 postdoctoral fellows, 32 graduate students, and 32 undergraduate students and was active from 2014-2019. Additionally, we forged new collaborations with nine researchers from DOE laboratories SRNL, LLNL, and ANL as well as the University of Manchester and the China Academy of Engineering Physics. The overarching goal of the project was to understand the conditions under which important classes of co-reactants, ranging from counter ions in crystal lattices to dissolved oxygen in pores, control the chemistry and transport characteristics of radionuclides in engineered waste forms and natural soils. Our approach was to characterize the time and length scales over which non-equilibrium states are maintained by rate-limiting, or rate-enhancing, reactions between radionuclides and co-reactants due to interactions between physical mass-transfer processes (i.e., advection, diffusion) and (biogeo) chemical reactions. We have focused our project on three specific classes of reactions relevant to radionuclide transport at DOE legacy sites: ion exchange/substitution, ligand complexation, and redox-mediated reactions. Understanding radionuclide migration requires detailed knowledge of how changes to a system – whether engineered or natural – drive the behavior of co-reactants, which in turn provide the geochemical context controlling radionuclide transport. Student engagement and training were a primary focus of the project in order to create a pipeline of researchers who could work in the area of nuclear waste disposition to support the state and the nation. Over the duration of the project we worked with 32 undergraduate, graduated 18 M.S. students and 14 Ph.D. students, and advised 13 postdoctoral fellows. The Ph.D students and postdocs have primarily taken positions at DOE laboratories, academia, and industry. Through our collaborative team, numerous follow on projects have been started with over $5M in sponsored research. Additionally, thus far the team has published 53 peer reviewed papers and given over 75 technical presentations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Intermediate time sub-diffusion and stress relaxation in ring polymer melts

The slow dynamics of non-concatenated ring melts remains a frontier problem in polymer science with implications for many soft material environments including cellular biophysics. Here, in this work, we report large-scale simulations of model ring melts that analyze the monomer and center-of-mass (CM) mean square displacements (MSD) and stress relaxation function on intermediate time and length scales. The degree of dynamical slowing down is characterized by the maximally sub-diffusive fractional time scaling exponents. The data span an exceptionally wide range of ring degrees of polymerization and stiffnesses and are not successfully organized based on the classic measure linear chain entanglement, N/N e . Rather, we find that the crossover degree of polymerization, N D , based on ring macromolecular caging that successfully allows master curves to be constructed for the long-time CM self-diffusion constant also collapses these temporal dynamic scaling exponents. Different properties display different exponents and exhibit one or two regimes of linear variation with the logarithm of N D / N . A distinct crossover of the CM-MSD and stress relaxation exponents emerges at sufficiently large N or stiffness that is not found for the monomer MSD, indicating a novel form of dynamic decoupling. This crossover aligns with the predicted critical degree of polymerization for transitioning from a weak to strong caging regime, indicative of activated transport. The latter may reflect the emergence of an intermolecular collective contribution to stress in analogy with dense soft colloidal matter. Suggestions are made for future theoretical work to address the rich patterns of behavior discovered.

Anomalous diffusion↗

On the Flow of a Cement Suspension: The Effects of Nano-Silica and Fly Ash Particles

Additives such as nano-silica and fly ash are widely used in cement and concrete materials to improve the rheology of fresh cement and concrete and the performance of hardened materials and increase the sustainability of the cement and concrete industry by reducing the usage of Portland cement. Therefore, it is important to study the effect of these additives on the rheological behavior of fresh cement. In this paper, we study the pulsating Poiseuille flow of fresh cement in a horizontal pipe by considering two different additives and when they are combined (nano-silica, fly ash, combined nano-silica, and fly ash). To model the fresh cement suspension, we used a modified form of the power-law model to demonstrate the dependency of the cement viscosity on the shear rate and volume fraction of cement and the additive particles. The convection–diffusion equation was used to solve for the volume fraction. After solving the equations in the dimensionless forms, we conducted a parametric study to analyze the effects of nano-silica, fly ash, and combined nano-silica and fly ash additives on the velocity and volume fraction profiles of the cement suspension. According to the parametric study presented here, larger nano-silica content results in lower centerline velocity of the cement suspension and larger non-uniformity of the volume fraction. Compared to nano-silica, fly ash exhibits an opposite effect on the velocity. Larger fly ash content results in higher centerline velocity, while the effect of the fly ash on the volume fraction is not obvious. For cement suspension containing combined nano-silica and fly ash additives, nano-silica plays a dominant role in the flow behavior of the suspension. The findings of the study can help the design and operation of the pulsating flow of fresh cement mortars and concrete in the 3D printing industry.

36 MATERIALS SCIENCE↗

Numerical simulation of flow and mixing in fracture intersections

Fluid transport through fractured geological formations is strongly influenced by the redistribution of solutes at fracture intersections. In this study, we perform detailed numerical simulations of flow and scalar transport within the intersection of two smooth, planar fractures. The analysis focuses on the mixing ratio, the proportion of solute flux exiting along the outlet branch aligned with the primary inlet flow direction, relative to the total solute flux at the outlets. We systematically investigate how the mixing ratio varies with four key parameters: Peclet number, Reynolds number, flow rate ratio between outlet branches, and fracture intersection angle. Results show that the mixing ratio decreases with increasing Peclet number and outlet flow rate ratio, consistent with reduced diffusive spreading and enhanced streamline routing. While low Reynolds numbers have minimal impact, inertial effects at higher Reynolds numbers significantly increase the mixing ratio. Additionally, acute and obtuse intersection angles alter flow partitioning and modify the solute distribution at the outlets. These findings provide a quantitative basis for incorporating physically realistic mixing behavior—intermediate between complete mixing and streamline-following assumptions—into network-scale transport models. The results have direct relevance to subsurface energy systems, including geothermal energy production, carbon sequestration, and contaminant remediation.

58 GEOSCIENCES↗

Complex Dependence of Calcite Crack Kinetics on Salinity: The Role of DLVO and Hydration Forces

Abstract Subcritical crack growth (SCG) plays an important role in many geological processes such as delayed earth rupture and rock weathering. The complex dependency of SCG on the in‐crack fluid chemistry, however, is still poorly understood. In this study, we utilize the newly developed surface force‐based fracture theory (SFFT) to elucidate the relative contributions of surface forces and solute transport to the crack growth kinetics of calcite in NaCl solutions. Expanding on Barenblatt's cohesive crack model, SFFT introduces an effective stress intensity at the crack tip that encompasses all the relevant intermolecular forces across the crack in addition to the external far‐field stresses. The nonlinear system of equations portraying the crack opening profile, the solute distribution in a propagating crack, and the crack growth velocity are numerically solved via an implicit scheme. After carefully calibrating the model for calcite‐water systems, the SFFT is used to predict the SCG response of calcite at different NaCl concentrations, based on various hypotheses. These predictions are then compared to existing SCG data from the literature. We demonstrate that the experimentally observed variation of SCG rate with NaCl concentration cannot be explained solely by DLVO forces (electrostatic and Van der Waals interactions). This can be remediated by introducing an exponentially decaying hydration force with a nonlinear, nonmonotonic dependence on NaCl concentration. Furthermore, we demonstrate that accounting for both diffusive and advective transport of ions is important in explaining the absence of a stage‐II SCG response for calcite in electrolyte solutions. Plain Language Summary Subcritical crack growth (SCG) refers to the slow propagation of cracks in materials under a stress below the threshold for catastrophic failure. SCG is a key process in many geological events, for example, delayed earth ruptures and rock weathering. New initiatives such as underground CO 2 and H 2 storage in carbonate reservoirs further call for better understanding of SCG in carbonate minerals subjected to varying fluid chemistry. This study examines the SCG of calcite, a key mineral found in carbonate rocks, intergranular cement in sandstones, and filling material in mineral veins and faults, determining their deformation and strength. A mathematical model is developed to describe how the crack opens and propagates, how solutes (like salts) distribute within the crack, and how the crack surfaces interact with each other. We used the model to predict calcite SCG in water at different salt concentrations and compared it with experimental data. Our results revealed that the hydration force is the dominating factor in determining the complex, non‐linear dependency of SCG on salinity. We also found that both the movement of ions by diffusion and by bulk water flow are crucial for explaining the SCG rates, especially when the cracks grow quickly. Key Points Surface Force‐Based Fracture Theory predicts the complex subcritical crack growth patterns of calcite crystals immersed in NaCl solutions Results highlight the dominant role of hydration forces in altering the fracture behavior of calcite compared to VdW and electric double‐layer forces Advective solute transport explains the absence of stages‐II and ‐III subcritical crack growth responses in solid‐liquid systems

DLVO↗

Neutronics Calculation Advances at Los Alamos: Manhattan Project to Monte Carlo

The history and advances of neutronics calculations at Los Alamos during the Manhattan Project through the present are reviewed. Substantial improvements to neutron diffusion methods and the invention of both the Monte Carlo neutron transport methods in 1947 and deterministic discrete ordinates Sn in 1953 were all made at Los Alamos just after the Manhattan Project. We briefly summarize early simpler and more approximate neutronics methods and then describe the need to better predict neutronics behavior through consideration of theoretical equations, models and algorithms, experimental measurements, and available computing capabilities and their limitations. This paper briefly covers key advances in deterministic methods during the Manhattan Project. These capabilities, coupled with increasing postwar defense needs and the invention of electronic computing with the Electronic Numeric Integrator and Computer, known as ENIAC, and the Mathematical Analyzer Numerical Integrator and Automatic Computer Model, known as MANIAC, led to the creation of Monte Carlo and deterministic discrete ordinates neutronics transport methods. We note the important role that the scientific comradery between the Los Alamos scientists played in the process. This paper briefly covers the early methods, algorithms, computers, and electronic and women pioneers that enabled Monte Carlo to spread to all areas of science. We focus heavily on these early developments and the subsequent creation of the MCNP® code, advances in its associated nuclear data, and its applications to problems of national defense at Los Alamos.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Role of Intradomain Heterogeneity on Ion and Polymer Dynamics in Block Polymer Electrolytes: Investigating Interfacial Mobility and Ion-Specific Dynamics and Transport

Block polymers show promise as solid-state battery electrolytes due to the optimization of conductive and mechanical properties enabled via tuning of block chemistry and length. We investigate a polystyrene-block-poly(oligo-oxyethylene methacrylate) (PS-b-POEM) electrolyte doped with various lithium salts to investigate the role of molecular structure on ion transport properties and on local ion dynamics and associations. Anion charge becomes more delocalized with increasing size, reducing the coupling between salt ions while increasing coupling between ion and polymer chain motions and creating a more mobile overall environment. We observe support for this ion-polymer coupling via 1 H, 7 Li and 19 F NMR spectroscopy, from which we obtain ion-specific mobility transition temperatures that differ from the polymer glass transition temperature. We also note faster transport and weaker local energetic interactions with anion size using temperature-dependent NMR diffusometry. 1 H NMR spectroscopy further elucidates polymer chain dynamics and enables quantification of the temperature-dependent fraction of the conducting block that is immobile near the PS-POEM domain interface. As a result, NMR thus represents a species-specific and timescale-specific platform to quantify phase and interface behavior, and to correlate ion-specific transport with polymer chain dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Your Clean Graphene is Still Not Clean

Researchers working with thin samples, such as monolayer graphene, are consistently struggling against contamination. Indeed, the problem of hydrocarbon contamination is known from the earliest days of electron microscopy and efforts to reduce this problem are ubiquitous to almost all high-vacuum experiments. Accurate knowledge of the behavior of such contamination is essential for electron beam (e-beam) based atomic fabrication, where it is aspired to select and control matter on an atom-by-atom basis. Here, the vexing question of hydrocarbon contamination on graphene is taken up. Image intensity is used to directly reveal the presence of diffusing hydrocarbons on ostensibly clean graphene. These diffusing hydrocarbons are previously inferred but not directly observed. Surprising dynamic variations of the concentration of these hydrocarbons impels questions about their origin. Here, some possible explanations are presented and some tentative conclusions are drawn. This work updates the conceptual model of “clean graphene” and offers refinements to the description of e-beam induced hydrocarbon deposition.

atomic fabrication↗

Macromolecular gelatin properties affect fibrin microarchitecture and tumor spheroid behavior in fibrin-gelatin gels

The biophysical properties of extracellular matrices (ECM) are known to regulate cell behavior, however decoupling cell behavior changes due to the relative contributions of material microstructure versus biomechanics or nutrient permeability remains challenging, especially within complex, multi-material matrices. We developed four gelatin-fibrin interpenetrating network (IPN) formulations which are identical in composition but possess variable gelatin molecular weight distributions, and display differences in microstructure, biomechanics, and diffusivity. In this work we interrogate the response of multicellular tumor spheroids to these IPN formulations and found that a high stiffness, gelatin-network dominated IPNs impeded remodeling and invasion of multicellular tumor spheroids; whereas relatively lower stiffness, fibrin-network dominated IPNs permitted protease-dependent remodeling and spheroid invasion. Cell proliferation correlated to nutrient diffusivity across tested IPN formulations. These findings demonstrate the complexity of ECM IPNs, relative to single polymer matrices, and highlight that cell response does not derive from a single aspect of the ECM, but rather from the interplay of multiple biomechanical properties. The methodology developed in this study represents a framework for future studies which aim to characterize cellular phenotypic responses to biophysical cues present within complex, multi-material matrices.

36 MATERIALS SCIENCE↗