skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Estimating viscoelastic, soft material properties using a modified Rayleigh cavitation bubble collapse time
Accurate determination of high strain rate (>103 1/s) constitutive properties of soft materials remains a formidable challenge. Albeit recent advancements among experimental techniques, in particular inertial microcavitation rheometry (IMR), the intrinsic requirement to visualize the bubble cavitation dynamics has limited its application to nominally transparent materials. Here, in an effort to address this challenge and to expand the experimental capability of IMR to optically opaque materials, we investigated whether one could use the acoustic signature of the time interval between the bubble's maximum radius and first collapse time point, characterized as the bubble collapse time, to infer the viscoelastic material properties without being able to image the bubble directly in the tissue. By introducing a modified Rayleigh collapse time for soft materials, which is strongly dependent on the stiffness of the material at hand, we show that, in principle, one can obtain an order of magnitude or better estimate of the viscoelastic material properties of the soft material under investigation. Using a newly developed energy-based theoretical framework, we show that for materials stiffer than 10 kPa the bubble collapse time during a single bubble cavitation event can provide quantitative and meaningful information about the constitutive properties of the material at hand. For very soft materials (i.e., shear modulus less than 10 kPa), our theory shows that unless the collapse time measurement has very high precision and low uncertainties, the material property estimates based on the bubble collapse time only will not be accurate and require visual resolution of the full cavitation kinematics.  more » « less
Award ID(s):
2232427 2232428
PAR ID:
10509387
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
AIP Publishing
Date Published:
Journal Name:
Physics of Fluids
Volume:
36
Issue:
1
ISSN:
1070-6631
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract BackgroundHydrogels are one of the most ubiquitous polymeric materials. Among them gelatin, agarose and polyacrylamide-based formulations have been effectively utilized in a variety of biomedical and defense-related applications including ultrasound-based therapies and soft tissue injury investigations stemming from ballistic and blast exposures. Interestingly, while in most cases accurate prediction of the mechanical response of these surrogate gels requires knowledge of the underlying finite deformation, high-strain rate material properties, it is these properties that have remained scarce in the literature. ObjectiveBuilding on our prior works using Inertial Microcavitation Rheometry (IMR), here we present a comprehensive list of the high-strain rate (> 10$$^3$$ 3 1/s) mechanical properties of these three popular classes of hydrogel materials characterized via laser-based IMR, further showing that the choice in finite-deformation, rate-dependent constitutive model can be informed directly by the type of crosslinking mechanism and resultant network structure of the hydrogel, thus providing a chemophysical basis of the the choice of phenomenological constitutive model. MethodsWe analyze existing experimental gelatin IMR datasets and compare the results with prior data on polyacrylamide. ResultsWe show that a Neo-Hookean Kelvin-Voigt (NHKV) model can suitably simulate the high-rate material response of dynamic, physically crosslinked hydrogels like gelatin, while the introduction of a strain-stiffening parameter through the use of the quadratic Kelvin-Voigt (qKV) model was necessary to appropriately model chemically crosslinked hydrogels such as polyacrylamide due to the nature of the static,covalent bonds that comprise their structure. ConclusionsIn this brief we show that knowledge of the type of underlying polymer structure, including its bond mobility, can directly inform the appropriate finite deformation, time-dependent viscoelastic material model for commonly employed tissue surrogate hydrogels undergoing high strain rate loading within the ballistic and blast regimes. 
    more » « less
  2. Unlike traditional structural materials, soft solids can often sustain very large deformation before failure, and many exhibit nonlinear viscoelastic behavior. Modeling nonlinear viscoelasticity is a challenging problem for a number of reasons. In particular, a large number of material parameters are needed to capture material response and validation of models can be hindered by limited amounts of experimental data available. We have developed a Gaussian Process (GP) approach to determine the material parameters of a constitutive model describing the mechanical behavior of a soft, viscoelastic PVA hydrogel. A large number of stress histories generated by the constitutive model constitute the training sets. The low-rank representations of stress histories by Singular Value Decomposition (SVD) are taken to be random variables which can be modeled via Gaussian Processes with respect to the material parameters of the constitutive model. We obtain optimal material parameters by minimizing an objective function over the input set. We find that there are many good sets of parameters. Further the process reveals relationships between the model parameters. Results so far show that GP has great potential in fitting constitutive models. 
    more » « less
  3. Erosion of solid surfaces due to cavitation has been studied for decades. However, it has been a long debate that which mechanism, namely shockwaves, microjets towards the surface, or both, during the cavitation bubble collapse is the primary factor responsible for that erosion. In this project we investigate the small-scale mechanisms of material erosion induced by the collapse of a single cavitation bubble close to a wall. More specifically, our experimental setup includes modification of the initial nucleus size, the maximum bubble radius, the stand-off distance to the wall, the material softness, and the initial flow temperature. We record the evolution of the bubble using high speed cameras as well as the local impacts on the materials. With the help of specifically designed cold-wires, we also measure the temperature in the liquid and in the bubble. Two different methods are used to generate the bubble: (i) an acoustic shockwave of variable intensity, (ii) a YAG laser, which may introduce a high temperature at the start. We also combine the two methods in which the laser initially creates a nucleus, then the shockwave triggers the expansion of the bubble. The objectives of the project are included in this paper, while some first results will be presented at the CAV2018 conference. 
    more » « less
  4. null (Ed.)
    The fluid dynamics of a bubble collapsing near an elastic or viscoelastic material is coupled with the mechanical response of the material. We apply a multiphase fluid–solid coupled computational model to simulate the collapse of an air bubble in water induced by an ultrasound shock wave, near different types of materials including metals (e.g. aluminium), polymers (e.g. polyurea), minerals (e.g. gypsum), glass and foams. We characterize the two-way fluid–material interaction by examining the fluid pressure and velocity fields, the time history of bubble shape and volume and the maximum tensile and shear stresses produced in the material. We show that the ratio of the longitudinal acoustic impedance of the material compared to that of the ambient fluid, $$Z/Z_0$$ , plays a significant role. When $$Z/Z_0<1$$ , the material reflects the compressive front of the incident shock into a tensile wave. The reflected tensile wave impinges on the bubble and decelerates its collapse. As a result, the collapse produces a liquid jet, but not necessarily a shock wave. When $$Z/Z_0>1$$ , the reflected wave is compressive and accelerates the bubble's collapse, leading to the emission of a shock wave whose amplitude increases linearly with $$\log (Z/Z_0)$$ , and can be much higher than the amplitude of the incident shock. The reflection of this emitted shock wave impinges on the bubble during its rebound. It reduces the speed of the bubble's rebound and the velocity of the liquid jet. Furthermore, we show that, for a set of materials with $$Z/Z_0\in [0.04, 10.8]$$ , the effect of acoustic impedance on the bubble's collapse time and minimum volume can be captured using phenomenological models constructed based on the solution of Rayleigh–Plesset equation. 
    more » « less
  5. Abstract Traumatic brain injury (TBI), particularly from explosive blasts, is a major cause of casualties in modern military conflicts. Computational models are an important tool in understanding the underlying biomechanics of TBI but are highly dependent on the mechanical properties of soft tissue to produce accurate results. Reported material properties of brain tissue can vary by several orders of magnitude between studies, and no published set of material parameters exists for porcine brain tissue at strain rates relevant to blast. In this work, brain tissue from the brainstem, cerebellum, and cerebrum of freshly euthanized adolescent male Göttingen minipigs was tested in simple shear and unconfined compression at strain rates ranging from quasi-static (QS) to 300 s−1. Brain tissue showed significant strain rate stiffening in both shear and compression. Minimal differences were seen between different regions of the brain. Both hyperelastic and hyper-viscoelastic constitutive models were fit to experimental stress, considering data from either a single loading mode (unidirectional) or two loading modes together (bidirectional). The unidirectional hyper-viscoelastic models with an Ogden hyperelastic representation and a one-term Prony series best captured the response of brain tissue in all regions and rates. The bidirectional models were generally able to capture the response of the tissue in high-rate shear and all compression modes, but not the QS shear. Our constitutive models describe the first set of material parameters for porcine brain tissue relevant to loading modes and rates seen in blast injury. 
    more » « less