An official website of the United States government
A comparative study is presented to solve the inverse problem in elasticity for the shear modulus (stiffness) distribution utilizing two constitutive equations: (1) linear elasticity assuming small strain theory, and (2) finite elasticity with a hyperelastic neo-Hookean material model. Assuming that a material undergoes large deformations and material nonlinearity is assumed negligible, the inverse solution using (2) is anticipated to yield better results than (1). Given the fact that solving a linear elastic model is significantly faster than a nonlinear model and more robust numerically, we posed the following question: How accurately could we map the shear modulus distribution with a linear elastic model using small strain theory for a specimen undergoing large deformations? To this end, experimental displacement data of a silicone composite sample containing two stiff inclusions of different sizes under uniaxial displacement controlled extension were acquired using a digital image correlation system. The silicone based composite was modeled both as a linear elastic solid under infinitesimal strains and as a neo-Hookean hyperelastic solid that takes into account geometrically nonlinear finite deformations. We observed that the mapped shear modulus contrast, determined by solving an inverse problem, between inclusion and background was higher for the linear elastic model as compared to that of the hyperelastic one. A similar trend was observed for simulated experiments, where synthetically computed displacement data were produced and the inverse problem solved using both, the linear elastic model and the neo-Hookean material model. In addition, it was observed that the inverse problem solution was inclusion size-sensitive. Consequently, an 1-D model was introduced to broaden our understanding of this issue. This 1-D analysis revealed that by using a linear elastic approach, the overestimation of the shear modulus contrast between inclusion and background increases with the increase of external loads and target shear modulus contrast. Finally, this investigation provides valuable information on the validity of the assumption for utilizing linear elasticity in solving inverse problems for the spatial distribution of shear modulus associated with soft solids undergoing large deformations. Thus, this work could be of importance to characterize mechanical property variations of polymer based materials such as rubbers or in elasticity imaging of tissues for pathology.
more »
« less
The effective shear modulus of a random isotropic suspension of monodisperse liquid n -spheres: from the dilute limit to the percolation threshold
A numerical and analytical study is made of the macroscopic or homogenized mechanical response of a random isotropic suspension of liquid n -spherical inclusions ( n = 2, 3), each having identical initial radius A , in an elastomer subjected to small quasistatic deformations. Attention is restricted to the basic case when the elastomer is an isotropic incompressible linear elastic solid, the liquid making up the inclusions is an incompressible linear elastic fluid, and the interfaces separating the solid elastomer from the liquid inclusions feature a constant initial surface tension γ . For such a class of suspensions, it has been recently established that the homogenized mechanical response is that of a standard linear elastic solid and hence, for the specific type of isotropic incompressible suspension of interest here, one that can be characterized solely by an effective shear modulus n in terms of the shear modulus μ of the elastomer, the initial elasto-capillary number eCa = γ /2 μA , the volume fraction c of inclusions, and the space dimension n . This paper presents numerical solutions—generated by means of a recently introduced finite-element scheme—for n over a wide range of elasto-capillary numbers eCa and volume fractions of inclusions c . Complementary to these, a formula is also introduced for n that is in quantitative agreement with all the numerical solutions, as well as with the asymptotic results for n in the limit of dilute volume fraction of inclusions and at percolation . The proposed formula has the added theoretical merit of being an iterated-homogenization solution.
more »
« less
- Award ID(s):
- 1922371
- PAR ID:
- 10461010
- Date Published:
- Journal Name:
- Soft Matter
- Volume:
- 19
- Issue:
- 2
- ISSN:
- 1744-683X
- Page Range / eLocation ID:
- 208 to 224
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
A microfluidic technique recently proposed in the literature to measure the interfacial tension between a liquid droplet and an immiscible suspending liquid [Hudson et al. , Appl. Phys. Lett. , 2005, 87 , 081905], [Cabral and Hudson, Lab Chip , 2006, 6 , 427] is suitably adapted to the characterization of the elastic modulus of soft particles in a continuous-flow process. A microfluidic device consisting of a cylindrical pipe with a reduction in cross-section is designed, and the deformation and velocity of incompressible elastic particles suspended in a Newtonian liquid are tracked as they move along the centerline through the constriction. Kinematic and shape information is exploited to calculate the particle's elastic modulus by means of the theory of elastic particle deformation in extensional flow. The approach is validated for different orders of magnitude of the elastic capillary number through experiments and numerical simulations.more » « less
-
The displacement of a suspension of particles by an immiscible fluid in a capillary tube or in porous media is a canonical configuration that finds application in a large number of natural and industrial applications, including water purification, dispersion of colloids and microplastics, coating and functionalization of tubings. The influence of particles dispersed in the fluid on the interfacial dynamics and on the properties of the liquid film left behind remain poorly understood. Here, we study the deposition of a coating film on the walls of a capillary tube induced by the translation of a suspension plug pushed by air. We identify the different deposition regimes as a function of the translation speed of the plug, the particle size, and the volume fraction of the suspension. The thickness of the coating film is characterized, and we show that similarly to dip coating, three coating regimes are observed, liquid only, heterogeneous, and thick films. We also show that, at first order, the thickness of films thicker than the particle diameter can be predicted using the effective viscosity of the suspension. Nevertheless, we also report that for large particles and concentrated suspensions, a shear-induced migration mechanism leads to local variations in volume fraction and modifies the deposited film thickness and composition.more » « less
-
We present a numerical study of non-colloidal spherical and rigid particles suspended in Newtonian, shear thinning and shear thickening fluids employing an immersed boundary method. We consider a linear Couette configuration to explore a wide range of solid volume fractions ( $$0.1\leqslant \unicode[STIX]{x1D6F7}\leqslant 0.4$$ ) and particle Reynolds numbers ( $$0.1\leqslant Re_{p}\leqslant 10$$ ). We report the distribution of solid and fluid phase velocity and solid volume fraction and show that close to the boundaries inertial effects result in a significant slip velocity between the solid and fluid phase. The local solid volume fraction profiles indicate particle layering close to the walls, which increases with the nominal $$\unicode[STIX]{x1D6F7}$$ . This feature is associated with the confinement effects. We calculate the probability density function of local strain rates and compare the latter’s mean value with the values estimated from the homogenisation theory of Chateau et al. ( J. Rheol. , vol. 52, 2008, pp. 489–506), indicating a reasonable agreement in the Stokesian regime. Both the mean value and standard deviation of the local strain rates increase primarily with the solid volume fraction and secondarily with the $$Re_{p}$$ . The wide spectrum of the local shear rate and its dependency on $$\unicode[STIX]{x1D6F7}$$ and $$Re_{p}$$ point to the deficiencies of the mean value of the local shear rates in estimating the rheology of these non-colloidal complex suspensions. Finally, we show that in the presence of inertia, the effective viscosity of these non-colloidal suspensions deviates from that of Stokesian suspensions. We discuss how inertia affects the microstructure and provide a scaling argument to give a closure for the suspension shear stress for both Newtonian and power-law suspending fluids. The stress closure is valid for moderate particle Reynolds numbers, $$O(Re_{p})\sim 10$$ .more » « less
-
Hydrodynamic interactions generate a diffusive motion in particulates in a shear flow, which plays seminal roles in overall particulate rheology and its microstructure. Here we investigate the shear induced diffusion in a red-blood cell (RBC) suspension using a numerical simulation resolving individual motion and deformation of RBCs. The non-spherical resting shape of RBCs gives rise to qualitatively different regimes of cell dynamics in a shear flow such as tank-treading, breathing, tumbling and swinging, depending on the cell flexibility determined by the elastic capillary number. We show that the transition from tumbling to tank-treading causes a reduction in the gradient diffusivity. The diffusivity is computed using a continuum approach from the evolution of a randomly packed cell-layer width with time as well as by the dynamic structure factor of the suspension. Both approaches, although operationally different, match and show that for intermediate capillary numbers RBCs cease tumbling accompanied by a drop in the coefficient of gradient diffusivity. A further increase of capillary number increases the diffusivity due to increased deformation. The effects of bending modulus and viscosity ratio variations are also briefly investigated. The computed shear induced diffusivity was compared with values in the literature. Apart from its effects in margination of cells in blood flow and use in medical diagnostics, the phenomenon broadly offers important insights into suspensions of deformable particles with non-spherical equilibrium shapes, which also could play a critical role in using particle flexibility for applications such as label free separation or material processing.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1039/d2sm01219g
