SEARCH · Search NASA

Results for “LOADS AND STRESSES, STRUCTURAL - TRANSIENT DYNAMIC”

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 55 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Phase Transitions in (Mg,Fe) 2 SiO 4 Olivine under Shock Compression

The response of forsterite, Mg 2 SiO 4 , under dynamic compression is of fundamental importance for understanding its phase transformations and high-pressure behavior. Here, we have carried out an in situ X-ray diffraction study of laser-shocked polycrystalline and single-crystal forsterite from 19 to 122 GPa using the Matter in Extreme Conditions end-station of the Linac Coherent Light Source. Under laser-based shock loading, forsterite does not transform to the high-pressure equilibrium assemblage of MgSiO 3 bridgmanite and MgO periclase, as has been suggested previously. Instead, we observe forsterite and forsterite III, a metastable polymorph of Mg 2 SiO 4 , coexisting in a mixed-phase region from 33 to 75 GPa for both polycrystalline and single-crystal samples. Densities inferred from X-ray diffraction are consistent with earlier gas-gun shock data. At higher stress, the response is sample-dependent. Polycrystalline samples undergo amorphization above 79 GPa. For [010]- and [001]-oriented crystals, a mixture of crystalline and amorphous material is observed to 108 GPa, whereas the [100]-oriented forsterite adopts an unknown phase at 122 GPa. The first two sharp diffraction peaks of amorphous Mg 2 SiO 4 show a similar trend with compression as those observed for MgSiO 3 in both recent static- and laser-driven shock experiments. This study provides new insight into the transformation of forsterite under nanosecond-duration shock loading. This work emphasizes the importance of formation of metastable phases along the Hugoniot and adds to evidence that the 300-K single-crystal diamond anvil cell experiments have relevance for understanding structures formed under shock compression. In particular, the metastable phase forsterite III has now been shown to form under dynamic compression from 10s to 100s of nanoseconds as well as under 300-K static compression. Upon compression to higher pressures, Mg 2 SiO 4 transforms to an amorphous phase. These results have broad relevance for understanding the behavior of silicates under dynamic compression.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Dynamic crushing of metal lattice metamaterials: Shock mode diagrams and transition to topology-independent compaction regime

Additively manufactured lattice metamaterials offer design versatility in strength and energy absorption and provide an additional degree of freedom through the selection of the lattice topology. Under quasistatic loading, the unit cell structure can strongly affect the stiffness, yield, and post-yield behavior, but whether and to what degree the effect of lattice topology persists into dynamic loading scenarios, up to the compaction shock regime, has not been established. LLNL’ s ALE3D hydrocode was used to perform a computational investigation of dynamic loading in multiple lattice types, including the gyroid, octet, Schwarz D, and rhombic dodecahedron, under impact velocities from 0.25 to 2.25 km/s. Shock Hugoniots for each lattice topology are generated and compared, suggesting that above a critical velocity, distinctions between architectures may not persevere and compacted lattices behave similarly. Here, to investigate the transition between topology-dependent quasistatic compression and the topology-independent regime above the critical velocity, a one-dimensional elastic-linear hardening plasticity-densified solid (E-LHP-DS) shock model for lattice materials was developed that relies upon confined compression to link the quasistatic and shock mechanics. Unlike similar works, the model does not assume rigid behavior prior to yield or locking behavior at densification, allowing a richer exploration of lattice mechanics. With only six parameters, the analytical model simultaneously fit quasistatic confined compression simulations for relative densities 0.1 $≤ \bar{ρ} ≤$ 0.9 and predicted dynamic compaction behavior to traverse several distinct shock modes, each defined by a critical impact speed (equivalently, critical stresses). Comparing the numerical results to the one-dimensional E-LHP-DS shock model predictions suggests that the topology-independence under strong shocks is linked to the onset of densification, which can be predicted based on quasistatic confined compression results.

Cellular material↗

DOE OSTI2026

Mechanics of pore array collapse and interaction in shock-compressed polymethyl methacrylate (PMMA)

Recent studies on dynamic pore collapse have revealed significant development of shear localization, which can lead to material failure in porous structures and hot spot generation in energetic materials. These findings have dramatically improved the understanding of failure mechanisms during pore collapse but also prompt further investigation of realistic porous materials. In particular, porous media consist of many pores and porous networks. Even in low-porosity materials, pores can form in close proximity during the manufacturing process, leading to the critical question of pore–pore interaction during collapse under dynamic loading conditions. This study investigates, via plate impact experiments coupled with high-speed internal digital image correlation and shadowgraphy techniques, the collapse of two pores in shock-compressed PMMA at stresses between 0.4 and 1 GPa. The results of these experiments provide new insights into shear localization in pore collapse, in addition to distinct interactions between pores. Shadowgraphy measurements reveal novel, direct visualization of shear band development and crack evolution from pore surfaces. Spacing between adiabatic shear bands is measured over a range of impact stresses and is predicted accurately by the Grady–Kipp model. Pore interactions are found to effect a transition in the impact stress threshold at which different failure mechanisms initiate and are also found to possibly influence preferential sites for shear cracking. Throughout the study, numerical and theoretical models are leveraged to understand shear localization behavior. The role of baroclinicity and wave interactions between the pores is used to elucidate interaction mechanisms between pores.

Lawlor, Barry P. [California Institute of Technolo↗

DOE OSTI2025

Multi-objective automatic discovery of optimized metamaterials for varying velocity impact protection

Mechanical metamaterials have demonstrated exceptional impact performance while remaining lightweight. Impact resistance has traditionally been investigated using quasi-static simulations, often with the assumption that performance will translate to high-velocity impact scenarios. However, critical crash protection parameters—such as peak stress and absorbed energy—are highly sensitive to impact velocity, leading to inconsistent performance under dynamic loading. To address this, we introduce a strain-rate-aware, active deep learning framework that enables multi-objective optimization of impact protection metrics across a wide range of impact velocities. Our framework captures the strain-rate sensitivity of architected lattices by learning to control spatial gradation in cellular metamaterials, resulting in over 200 % enhancement in impact protection relative to state-of-the-art designs such as Voronoi and re-entrant lattices. We demonstrate its practical utility by designing next-generation lattice structures for automotive bumper systems that satisfy multiple, velocity-specific safety criteria—capabilities beyond those of conventional designs. More than just a predictive tool, this framework marks the first step towards enabling adaptable impact-resistant structures across dynamic regimes.

Deep learning↗

DOE OSTI2021

Strongly Anisotropic Thermomechanical Response to Shock Wave Loading in Oriented Samples of the Triclinic Molecular Crystal 1,3,5-Triamino-2,4,6-trinitrobenzene

All-atom molecular dynamics (MD) simulations were used to study shock wave loading in oriented single crystals of the highly anisotropic triclinic molecular crystal 1,3,5-triamino-2,4,6-trinitrobenzene (TATB). The crystal structure consists of planar hydrogen-bonded sheets of individually planar TATB molecules that stack into graphitic-like layers. Shocks were studied for seven systematically prepared crystal orientations with limiting cases that correspond to shock propagation exactly perpendicular and exactly parallel to the graphitic-like layers. The simulations were performed for initially defect-free crystals using a reverse-ballistic configuration that generates explicit, supported shocks. Final longitudinal stress components are between ≈8.5 and ≈10.5 GPa for the 1.0 km s –1 impact speed studied. Orientation-dependent properties are reported including shock speeds, stresses, temperatures, compression ratios, and local material strain rates. Spatiotemporal maps of the temperature, stress tensor, material flow, and molecular orientations reveal complicated processes that arise for specific shock directions. Furthermore, the results indicate that TATB shock response is highly sensitive to crystal orientation, with significant qualitative differences for the time evolution of the stress tensor and temperature, elastic/inelastic compression response, defect formation and growth, critical von Mises stress, and strain rates during shock rise that span nearly an order of magnitude. A variety of inelastic deformation mechanisms are identified, ranging from crumpling of graphitic-like layers to dislocation-mediated plasticity to intense shear strain localization. To our knowledge, these are the first systematic MD simulations and analysis of explicit shock wave propagation along nontrivial crystal directions in a triclinic molecular crystal.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Multiphysics Analyses of the Bottom Components of the 3D Printed Transformational Challenge Reactor

This research presents multi-physics analyses on the bottom components of the Transformational Challenge Reactor (TCR) facility. These components include the bottom axial reflector and the steel exit cone. The bottom axial reflector is made of pure silicon carbide elements hosting helium cooling channels These elements are 3D printed and therefore can host any arbitrary shape of the helium cooling channels. The design of the bottom reflector considers the neutronics and thermo-fluid dynamics performances as well as the manufacturing process optimization. More precisely, the best design of the bottom reflector reduces neutron leakage by avoiding straight cylindrical helium channels that facilitate neutron leakage, minimizes the helium flow pressure drop, and reduces the number of 3D printed silicon carbide pieces. The exit cone steel structure collects the hot helium from the bottom fuel assemblies and channels the cold helium to the top of the fuel assemblies. The steel simultaneous contact with hot and cold helium flows sets a large thermal gradient. Different designs of the exit cone are proposed to reduce the steel equivalent stress from the helium thermal load. The multi-physics analyses have been performed using Ansys Fluent, Ansys Mechanical, STAR-CCM+, and Serpent computer programs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

Anisotropic strength behavior of single-crystal TATB

High-rate strength behavior plays an important role in the shock initiation of high explosives, with plastic deformation serving to localize heat into hot spots and as a mechanochemical means to enhance reactivity. Recent simulations predict that detonation-like shocks produce highly reactive nanoscale shear bands in the layered crystalline explosive TATB (1,3,5-triamino-2,4,6-trinitrobenzene), but the thresholds leading to this response are poorly understood. We utilize molecular dynamics to simulate the high-rate compressive stress–strain response of TATB, with a focus on understanding flow behavior. Here, the dependence of strength on pressure and loading axis (crystal orientation) is explored. The deformation mechanisms fall broadly into two categories, with compression along crystal layers activating a buckling/twinning mode and compression normal to the layers producing nanoscale shear bands. Despite the complexity of the underlying mechanisms, the crystal exhibits relatively straightforward stress–strain curves. Most of the crystal orientations studied show rapid strain softening following the onset of yielding, which settles to a steady flow state. Trajectories are analyzed using five metrics for local states and structural order, but most of these metrics yield similar distributions for these deformation mechanisms. On the other hand, a recently proposed measure of intramolecular strain energy is found to most cleanly distinguish between these mechanisms, while also providing a plausible connection with mechanochemically accelerated decomposition kinetics. Localization of intramolecular strain energy is found to depend strongly on crystal orientation and pressure.

36 MATERIALS SCIENCE↗

DOE OSTI2021

On the computational solution of vector-density based continuum dislocation dynamics models: A comparison of two plastic distortion and stress update algorithms

Continuum dislocation dynamics models of mesoscale plasticity consist of dislocation transport-reaction equations coupled with crystal mechanics equations. The coupling between these two sets of equations is such that dislocation transport gives rise to the evolution of plastic distortion (strain), while the evolution of the latter fixes the stress from which the dislocation velocity field is found via a mobility law. Earlier solutions of these equations employed a staggered solution scheme for the two sets of equations in which the plastic distortion was updated via time integration of its rate, as found from Orowan's law. In this work, we show that such a direct time integration scheme can suffer from accumulation of numerical errors. We introduce an alternative scheme based on field dislocation mechanics that ensures consistency between the plastic distortion and the dislocation content in the crystal. The new scheme is based on calculating the compatible and incompatible parts of the plastic distortion separately, and the incompatible part is calculated from the current dislocation density field. Stress field and dislocation transport calculations were implemented within a finite element based discretization of the governing equations, with the crystal mechanics part solved by a conventional Galerkin method and the dislocation transport equations by the least squares method. A simple test was first performed to show the accuracy of the two schemes for updating the plastic distortion, which shows that the solution method based on field dislocation mechanics is more accurate. This method then was used to simulate an austenitic steel crystal under uniaxial loading and multiple slip conditions. By considering dislocation interactions caused by junctions, a hardening rate similar to discrete dislocation dynamics simulation results was obtained. Finally, the simulations show that dislocations exhibit some self-organized structures as the strain is increased.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Water Weakening of Calcium Oxide

Chemomechanical processes such as water weakening can control the permeability and deformation of rocks and manmade materials. Here, atomistic modeling and nanomechanical experiments were used to identify molecular origins of chemomechanical effects in calcium oxide (CaO) and its effect on observed elastic, plastic, and brittle deformation. Classical molecular dynamics simulations with the bond order-based reactive force-field ReaxFF were used to assess brittle fracture. In the presence of water, CaO fractured earlier and more often during quasi-static loading, with a calculated reduction in fracture toughness of ~80% associated with changes in the stress distribution around the crack tip. Experimentally, elastic and plastic deformation of CaO surfaces exposed to water was assessed experimentally using in situ liquid nanoindentation. Nanoindentation showed that following exposure to water, the contact hardness decreased by 1–2 orders of magnitude and decreased the modulus by 2–3 orders of magnitude due to surface hydroxylation. Furthermore, the strong chemomechanical effects on the mechanical processes in CaO suggests that minerals with similar structures may exhibit comparable effects, influencing the stability of cements and geomaterials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Biocompatible Co–P Metallic Glasses with Superior Degradation Tolerance in Physiological Environments

Metallic glasses represent a class of metallic alloys with a fully amorphous structure and attractive properties, making them promising in bioimplant applications. Here, the degradation tolerance of biocompatible cobalt–phosphorus (Co–P) metallic glasses was studied in a simulated physiological environment. The metallic glasses were synthesized in the form of coatings through a facile electrodeposition approach. This method utilizes their outstanding surface characteristics and bypasses the size limitations usually associated with their bulk counterparts. Further, the Co–P alloys showed exceptional tribological response with ~14% lower coefficient of friction and 2 orders of magnitude lesser wear rate compared to SS316 stainless steel. In addition, the Co–P alloys showed a 3 times higher hardness and 4 times higher hardness/modulus ratio compared to SS316, indicating better elastic recovery under dynamic shear stresses that are common in load-bearing bioimplants. The Co–P metallic glasses exhibited excellent hemocompatibility and cytocompatibility in terms of lower platelet adhesion, spreading, and aggregation, a hemolysis ratio lower than 1%, and enhanced surface wettability, suggesting a superlative performance in bioimplant applications.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Understanding the phase transformation mechanisms that affect the dynamic response of Fe-based microstructures at the atomic scales

Large-scale molecular dynamics (MD) simulations were carried out to investigate the shock-induced evolution of microstructure in Fe-based systems comprising single-crystal and layered Cu/Fe alloys with a distribution of interfaces. The shock compression of pure single-crystal Fe oriented along [110] above a threshold pressure results in a BCC (α)→HCP (ε) phase transformation behavior that generates a distribution of ε phase variants in the phase transformed region of the microstructure behind the shock front. The propagation of the release wave through a phase transformed ε phase causes a reverse ε→α phase transformation and renders a distribution of twins for the [110] oriented Fe that serve as void nucleation sites during spall failure. The simulations reveal that the α→ε→α transformation-induced twinning for shock loading along the [110] direction is due to a dominant ε phase variant formed during compression that rotates on the arrival of the release wave followed by a reverse phase transformation to twins in the α phase. The modifications in the evolution of the ε phase variants and twins in Fe behavior are also studied for Cu–Fe layered microstructures due to the shock wave interactions with the Cu/Fe interfaces using a newly constructed Cu–Fe alloy potential. Here, the MD simulations suggest that interfaces affect the observed variants during shock compression and, hence, distributions of twins during shock release that affects the void nucleation stresses in the Fe phase of Cu/Fe microstructures.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Pseudoelastic deformation in Mo-based refractory multi-principal element alloys

Phase diagrams supported by density functional theory methods can be crucial for designing high-entropy alloys that are subset of multi-principal-element alloys. We present phase and property analysis of quinary (MoW) x Zr y (TaTi) 1-x-y refractory high-entropy alloys from combined Calculation of Phase Diagram (CALPHAD) and density-functional theory results, supplemented by molecular dynamics simulations. Both CALPHAD and density-functional theory analysis of phase stability indicates a Mo-W-rich region of this quinary has a stable single-phase body-centered-cubic structure. We report first quinary composition from Mo-W-Ta-Ti-Zr family of alloy with pseudo-elastic behavior, i.e., hysteresis in stress-strain. Our analysis shows that only Mo-W-rich compositions of Mo-W-Ta-Ti-Zr, i.e., Mo+W, show reproducible hysteresis in stress-strain responsible for pseudo-elastic behavior. The (MoW) 85 Zr 7.5 (TaTi) 7.5 was down-selected based on temperature-dependent phase diagram analysis and molecular dynamics simulations predicted elastic behavior that reveals twinning-assisted pseudoelastic behavior. While mostly unexplored in body-centered-cubic crystals, twinning is a fundamental deformation mechanism that competes against dislocation slip in crystalline solids. Furthermore, this alloy shows identical cyclic deformation characteristics during uniaxial < 100 > loading, i.e., the pseudoelasticity is isotropic in loading direction. Additionally, a temperature increase from 77 to 1,500 K enhances the elastic strain recovery in load-unload cycles, offering possibly control to tune the pseudoelastic behavior.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Polymorphism of gold under laser-based ramp compression to 690 GPa

The high-pressure behavior and polymorphism of gold has attracted strong recent interest from theory as well as from both static and dynamic experiments. In this study, gold was dynamically ramp-compressed to 690 GPa via laser ablation using the Omega Laser Facility (U. Rochester) in order to investigate the stability regions of this noble metal's polymorphs. Here, stresses were determined from measurements of surface or interface velocities by laser interferometry. In situ X-ray diffraction was carried out to constrain the structure and density. The ambient face-centered-cubic (fcc) phase was observed at 162 GPa, followed by a mixed phase region consisting of the fcc and body-centeredcubic (bcc). The presence of only the bcc phase was observed between 377 and 690 GPa. Our results are consistent with recent shock compression experiments and demonstrate that gold transforms to the bcc phase under conditions of both ramp and shock loading.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Artificial neural network potentials for mechanics and fracture dynamics of two-dimensional crystals **

Understanding the mechanics and failure of materials at the nanoscale is critical for their engineering and applications. The accurate atomistic modeling of brittle failure with crack propagation in covalent crystals requires a quantum mechanics-based description of individual bond-breaking events. Artificial neural network potentials (NNPs) have emerged to overcome the traditional, physics-based modeling tradeoff between accuracy and accessible time and length scales. Previous studies have shown successful applications of NNPs for describing the structure and dynamics of molecular systems and amorphous or liquid phases of materials. However, their application to deformation and failure processes in materials is still uncommon. In this study, we discuss the apparent limitations of NNPs for the description of deformation and fracture under loadings and propose a way to generate and select training data for their employment in simulations of deformation and fracture simulations of crystals. We applied the proposed approach to 2D crystalline graphene, utilizing the density-functional tight-binding method for more efficient and extensive data generation in place of density functional theory. Then, we explored how the data selection affects the accuracy of the developed artificial NNPs. It revealed that NNP’s reliability should not only be measured based on the total energy and atomic force comparisons for reference structures but also utilize comparisons for physical properties, e.g. stress–strain curves and geometric deformation. In sharp contrast to popular reactive bond order potentials, our optimized NNP predicts straight crack propagation in graphene along both armchair and zigzag (ZZ) lattice directions, as well as higher fracture toughness of ZZ edge direction. Our study provides significant insight into crack propagation mechanisms on atomic scales and highlights strategies for NNP developments of broader materials.

2D materials↗

DOE OSTI2021

Analyzing Mechanical integrity of Microwave resonant Cavity flowmeter

The microwave resonant cavity flow sensor has been proposed to address the difficulties of measuring high-temperature molten salt flow in reactor vessels. The flow sensor consists of a metallic cylindrical cavity coupled to a microwave waveguide. The cylindrical cavity has one of the flat surfaces thin enough to undergo microscopic deflection due to dynamic pressure of the exerted by the flowing fluid. Deflection of the membrane causes a shift in the microwave resonant frequency. Fluid velocity can be estimated from the calibration curve. The objective of this report is to investigate mechanical integrity of the flowmeter’s membrane. The membrane is a thin stainless-steel circular plate with a clamped edge, subjected to uniform load. By calculating the stress on the plate due to deflection and, comparing the stress to the material ultimate tensile strength and yield strength, it can be estimated if the plate will fail. If the maximum stress is greater than the yield strength, the plate will undergo plastic deformation, and the sensor will be disabled. A stress greater than the ultimate tensile strength will cause the plate to fracture. The stress on the plate was calculated with an analytic closed form solution model, and with COMSOL Structural Mechanics Module which does not involve any approximations. Stresses in both the analytic model and COMSOL numerical model were calculated for stainless steel 316 thin circular plates under room temperature conditions, where the fluid is water. The diameter of the circular plate is approximately one inch, and thickness is 10mil. Fluid flow velocity is in the range from 0.5m/s to 1m/s. The results from the two models are similar. Both the analytic model and COMSOL model showed that maximum stresses on the plate, which are at the radial boundary of the plate, are three orders of magnitude smaller than the yield strength and ultimate tensile strength. This indicates that the sensor is at a low risk of mechanical failure.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2023

Investigations of correlation and coherence in turbulence from a large-eddy simulation

Abstract. Microscale flow descriptions are often given in terms of mean quantities, turbulent kinetic energy, and/or stresses. Those metrics, while valuable, give limited information about turbulent eddies and coherent turbulent structures. This work investigates the structure of an atmospheric boundary layer using coherence and correlation in space and time with a range of separation distances. We calculate spatial correlations over entire planes of velocity fluctuations, from which we can evaluate the correlation along different directions at different spacings. Similarly, coherence of the three velocity components over separations in the three directions is also investigated. We apply these analyses to a mesoscale–microscale coupled scenario with time-varying conditions and examine nuances in spatial correlations that are often overlooked. Through these analyses and results, this work highlights important differences observed in terms of coherence when comparing large-eddy simulation data to simpler models and suggests ways to improve these simpler models. We note that such differences are important for disciplines like wind energy structural dynamic analysis, in which blade loading and fatigue depend strongly on the structure of the turbulence. We emphasize the additional wealth of data that can be provided by typical atmospheric boundary layer large-eddy simulation when correlation and coherence analysis is included, and we also state the limitations of large-eddy simulation data, which inherently truncate the smaller scales of turbulence.

17 WIND ENERGY↗

DOE OSTI2024

Near wake evolution of a tidal stream turbine due to asymmetric sheared turbulent inflow with different integral length scales

Tidal stream turbines deployed at highly energetic open water sites are subjected to sheared inflow in the rotor plane. The inflow shear is expected to cause asymmetric loading on the rotor blades and affect the downstream wake. In the current study, two different turbulent inflow conditions, static-high shear and dynamic shear, were generated via an active-grid turbulence generator. A 1:20 scaled three-bladed horizontal axis tidal turbine model was tested in those conditions. The results were compared to a quasi-laminar case with no imposed turbulence or shear. The results show that the high shear reduces the average performance, with a drop of up to 16% in the optimal power coefficient. Besides, the shear profiles increase torque fluctuations and induce significant differences in wake hydrodynamics between the high-speed (upper) and low-speed (lower) regions. The large integral length scales further enhance the load fluctuations perceived by the rotor but have a negligible effect on the mean wake field quantities and the wake recovery. The lower half region featured a faster breakdown of tip vortex structure and a rapid drop of swirl number, a phenomenon conjectured to be a consequence of the strong turbulence intensities and Reynolds stresses in the lower half region. Furthermore, the sheared turbulent inflow also results in a very intensive energy redistribution process towards large-scale, low-frequency motions, which is important to the downstream turbines.

16 TIDAL AND WAVE POWER↗

DOE OSTI2025

Mesoscale Modeling of Hydrogels Under Frictional Shear Stress

Hydrogels are three-dimensional networks of hydrophilic polymers often used as a simplified model of hydrated biological materials, from cartilaginous joints to the ocular tear film. However, the lubrication mechanisms of hydrogels remain poorly understood, partly due to their complex polymeric structure, which creates blurred interfaces during sliding that are challenging to study experimentally. In this study, we employ dissipative particle dynamics (DPD) to investigate the frictional behavior of a polymeric hydrogel network sliding against a solid wall in an explicit viscous solvent. This computational approach enables us to model hydrodynamic interactions and mesoscale polymer dynamics, capturing key aspects of hydrogel friction. Our simulations reveal that hydrogel friction is governed by the interplay between polymer relaxation and viscous shear, characterized by the Weissenberg number (Wi). At low Wi, friction coefficient remain nearly constant, dominated by polymer relaxation. However, at higher Wi, friction is dominated by viscous drag within a near-wall solvent layer, leading to a linear increase in friction coefficient with Wi. Furthermore, our results demonstrate an inverse relationship between the friction coefficient and the applied normal load, consistent with experimental observations. This work provides new insights into the fundamental tribological properties of hydrogels, shedding light on the micromechanics of hydrogel friction. Improving our understanding of hydrogel structure and dynamics under friction advances our knowledge of the mechanisms regulating biological lubrication in health and disease.

36 MATERIALS SCIENCE↗