Search

Hydraulic fracturing of unconventional reservoirs has seen a boom in the last century, as a means to fulfill the growing energy demand in the world. The fracturing fluid used in the process plays a substantial role in determining the results. Hence, several research and development efforts have been geared towards developing more sustainable, efficient, and improved fracturing fluids. Herein, we present a dynamic binary complex (DBC) solution, with potential to be useful in the hydraulic fracturing domain. It has a supramolecular structure formed by the self-assembly of low molecular weight viscosifiers (LMWVs) oleic acid and diethylenetriamine into an elongated entangled network under alkaline conditions. With less than 2 wt% constituents dispersed in aqueous solution, a viscous gel that exhibits high viscosities even under shear was formed. Key features include responsiveness to pH and salinity, and a zero-shear viscosity that could be tuned by a factor of ~280 by changing the pH. Furthermore, its viscous properties were more pronounced in the presence of salt. Sand settling tests revealed its potential to hold up sand particles for extended periods of time. In conclusion, this DBC solution system has potential to be utilized as a smart salt-responsive, pH-switchable hydraulic fracturing fluid.

The rheological properties parameters of paddy soil affect the interaction between the tillage tools and soil, thus influencing the operation quality and power consumption. In order to study the effects of tillage methods and moisture content on the rheological properties parameters of paddy soil in the middle and lower reaches of the Yangtze River, uniaxial compression creep tests of paddy soils with four moisture contents under no tillage (moisture contents: 26.71%, 24.52%, 23.26%, 21.28%) and plough tillage (moisture contents: 26.77%, 25.55%, 23.40%, 20.56%) were carried out using a TMS-PRO texture analyzer. The creep properties curves obtained from the tests, and the rheological constitutive equation of paddy soil under compression was established by Burgers viscoelastic model. Respectively, the quantitative change rules of creep properties of paddy soil with different moisture contents under different tillage methods and the correlation between these parameters were explored. The results showed that the moisture content under the three-year plough tillage and no tillage methods had significant influence on the rheological properties parameters of paddy soil (P < 0.05). The instantaneous elastic modulus, delay elastic modulus, and viscosity coefficient of the two paddy soils (no tillage and plough tillage soils) decreased with the increase of moisture content. However, the variation rules of relaxation time and delay viscosity coefficient with moisture content differed between these two paddy soils. Specifically, the strain rate of the two paddy soils decreased as moisture content decreased, where the total strain combines elastic strain, viscous strain, and viscoelastic strain. The initial strain rate and steady strain rate of the plough tillage paddy soils were lower than that of the no tillage paddy soils. The established creep model equation could be used to obtain viscoelastic rheological parameters of paddy soil in a wide range. The fitting equations between rheological parameters and moisture content were introduced into Burgers model, and the coupling equations between creep deformation and moisture content and time were derived, which could be used to predict the creep properties and deformation behavior of paddy soil in a certain range of no tillage and ploughed field. Overall, this study has a certain theoretical significance for the development and improvement of paddy soil rheology theory, and can also provide theoretical basis and technical support for the research of agricultural machinery design optimization, field water, soil conservation, soil tillage and compaction related simulation analysis in the middle and lower reaches of the Yangtze River.

Recent studies showed that wave intensity analysis (WIA) provides clinically valuable information about local and global cardiovascular function. Wave intensity (WI) is computed as the product of the pressure change and the velocity change during short time intervals. The major limitation of WIA in clinical practice is the need for invasive pressure measurement. Since vessel wall displacement can be measured non-invasively, the usage of WI will be expanded if the vessel wall dilation is used instead of pressure in derivation of WI waveform. Our goal in this study is to investigate the agreement between wall displacement-based WI and the pressure-based WI for different vessel wall models including linear elastic, nonlinear and viscoelastic cases. The arbitrary Eulerian Lagrangian finite element method is employed to solve the coupled fluid-structure interaction (FSI). Our computational models also include two types of vascular disease-related cases with geometrical irregularities, aneurysm and stenosis. Our results show that for vessels with linear elastic wall, the displacement-based WI is almost identical to the pressure-based WI. The existence of vessel irregularities does not impact the accuracy of displacement-based WI. However, in a viscoelastic wall where there is a phase difference between pressure and vessel wall dilation, displacement-based WI deviated from pressure-based WI. The error associated with this phase difference increased nonlinearly with increasing viscosity. This results in a maximum error of 6.8% and 7.13% for a regular viscoelastic vessel wall and an irregular viscoelastic vessel wall, respectively. A separate analysis has also been performed on the agreement of backward and forward running waves extracted from a decomposition of the displacement-based and pressure-based WI. Our findings suggest that displacement-based WI is a reliable method of WIA for large central arteries that do not show viscoelastic behaviors. This can be clinically significant since the required information can be measured non-invasively.

Cross-linked actin networks are the primary component of the cell cytoskeleton and have been the subject of numerous experimental and modeling studies. While these studies have demonstrated that the networks are viscoelastic materials, evolving from elastic solids on short timescales to viscous fluids on long ones, questions remain about the duration of each asymptotic regime, the role of the surrounding fluid, and the behavior of the networks on intermediate timescales. Here we perform detailed simulations of passively cross-linked non-Brownian actin networks to quantify the principal timescales involved in the elastoviscous behavior, study the role of nonlocal hydrodynamic interactions, and parameterize continuum models from discrete stochastic simulations. To do this, we extend our recent computational framework for semiflexible filament suspensions, which is based on nonlocal slender body theory, to actin networks with dynamic cross linkers and finite filament lifetime. We introduce a model where the cross linkers are elastic springs with sticky ends stochastically binding to and unbinding from the elastic filaments, which randomly turn over at a characteristic rate. We show that, depending on the parameters, the network evolves to a steady state morphology that is either an isotropic actin mesh or a mesh with embedded actin bundles. For different degrees of bundling, we numerically apply small-amplitude oscillatory shear deformation to extract three timescales from networks of hundreds of filaments and cross linkers. We analyze the dependence of these timescales, which range from the order of hundredths of a second to the actin turnover time of several seconds, on the dynamic nature of the links, solvent viscosity, and filament bending stiffness. We show that the network is mostly elastic on the short time scale, with the elasticity coming mainly from the cross links, and viscous on the long time scale, with the effective viscosity originating primarily from stretching and breaking of the cross links. We show that the influence of nonlocal hydrodynamic interactions depends on the network morphology: for homogeneous meshworks, nonlocal hydrodynamics gives only a small correction to the viscous behavior, but for bundled networks it both hinders the formation of bundles and significantly lowers the resistance to shear once bundles are formed. We use our results to construct three-timescale generalized Maxwell models of the networks., Author summary The cytoskeleton is composed of semiflexible, inextensible actin filaments, dynamic cross linkers, and molecular motors, and makes the primary contribution to the structural properties of the cell. Despite its being so fundamental to cell biology, the biological complexity of the cytoskeleton hinders our ability to understand its mechanical properties through in vitro experiments. In this paper, we perform microscopic simulations of actin fibers and transient cross linkers to quantify the principle timescales involved in the network, study how these timescales influence the morphology and rheology of the system, and examine the role of hydrodynamic interactions in cytoskeletal networks. We find three principle timescales which we associate with fiber flexibility, cross linker detachment, and network remodeling, respectively. We show that the morphology of the network is more important on longer timescales, where the viscosity of links inside of fiber bundles is enhanced. We also show that hydrodynamic interactions reduce the stress inside of bundles because of entrainment flows. Finally, we propose a continuum model which can be used to coarse-grain our agent-based simulations and enable modeling of larger systems.

Hyaluronic acid injections have been a mainstay of arthritis treatment for decades. However, much controversy remains about their clinical efficacy and their potential mechanism of action. This approach to arthritis therapy is often called viscosupplementation, a term which is rooted in the elevated viscosity of the injected solutions. This terminology also suggests a mechanical pathway of action and further implies that their efficacy is dependent on viscosity. Notably, previous studies of the relationship between viscous properties of hyaluronic acid solutions and their clinical efficacy have not been definitive. Recently we developed an experimental and analytical framework for studying cartilage lubrication that captures the Stribeck-like behavior of cartilage in an elastoviscous transition curve. Here we apply this framework to study the lubricating behavior of six hyaluronan products currently used for injectable arthritis therapy in the US. Despite the fact that the source and chemical modifications endow these products with a range of lubricating properties, we show that the lubricating effect of all of these materials can be described by this Stribeck-like elastoviscous transition. Fitting this data to the elastoviscous transition model enables the calculation of effective lubricating viscosities for each material, which differ substantially from the viscosities measured using standard rheometry. Further we show that while data from standard rheometry are poor predictors of clinical performance of these materials, measurements of friction coefficient and effective lubricating viscosity correlate well (R2 = 0.77; p < 0.005) with assessments of improved clinical function reported previously. This approach offers both a novel method that can be used to evaluate potential clinical efficacy of hyaluronic acid formulations and provide new insight on their mode of action.

Wheat (Triticum spp.) gluten consists mainly of intrinsincally disordered storage proteins (glutenins and gliadins) that can form megadalton-sized networks. These networks are responsible for the unique viscoelastic properties of wheat dough and affect the quality of bread. These properties have not yet been studied by molecular level simulations. Here, we use a newly developed α-C-based coarse-grained model to study ∼ 4000-residue systems. The corresponding time-dependent properties are studied through shear and axial deformations. We measure the response force to the deformation, the number of entanglements and cavities, the mobility of residues, the number of the inter-chain bonds, etc. Glutenins are shown to influence the mechanics of gluten much more than gliadins. Our simulations are consistent with the existing ideas about gluten elasticity and emphasize the role of entanglements and hydrogen bonding. We also demonstrate that the storage proteins in maize and rice lead to weaker elasticity which points to the unique properties of wheat gluten., Author summary During the breadmaking process, expanding gas bubbles cause the dough to increase volume. Gluten proteins act as an elastic scaffold in that process, allowing the wheat dough to grow more than other kinds of dough. Thus, explaining the unique viscoelastic properties of gluten at the molecular level may be of great interest to the baking industry. Assessing the structural properties of gluten is difficult because its proteins are disordered. We provide the first molecular dynamics model of gluten elasticity, that is able to distinguish gluten and proteins from different plants, like maize and rice. Our model shows the structural changes the gluten proteins undergo during their deformation, which mimics the mixing of dough during kneading. It also allows for a determination of the force required to extend gluten proteins, as during baking. The data confirms existing theories about gluten, but it also provides molecular-level information about the extraordinary elasticity of gluten.

Dairy products exhibit several physical properties that are crucial to define whether we like the food or not: firmness, creaminess, thickness, or lightness. Viscosity changes the flow properties of food and influences the appearance and the consistency of a product; this control variable is important in most production stages-manufacture, processing, and storage. Viscosity of heterogeneous products at a given temperature depends on its composition and physical state of its substances. Although rheology provides a method to access the product viscosity, it lacks non-contact full-field monitoring. We apply a simple correlation analysis of laser speckle images to evaluate viscosity properties of dairy products. Our approach ensures robust measurements with high degree of detectability.

Knowledge of the biomechanical properties of the human cornea is crucial for understanding the development of corneal diseases and impact of surgical treatments (e.g., corneal laser surgery, corneal cross-linking). Using a Surface Force Apparatus we investigated the transient viscous response of the anterior cornea from donor human eyes compressed between macroscopic crossed cylinders. Corneal biomechanics was analyzed using linear viscoelastic theory and interpreted in the framework of a biphasic model of soft hydrated porous tissues, including a significant contribution from the pressurization and viscous flow of fluid within the corneal tissue. Time-resolved measurements of tissue deformation and careful determination of the relaxation time provided an elastic modulus in the range between 0.17 and 1.43 MPa, and fluid permeability of the order of 10(-13) m(4)/(N.s). The permeability decreased as the deformation was increased above a strain level of about 10%, indicating that the interstitial space between fibrils of the corneal stromal matrix was reduced under the effect of strong compression. This effect may play a major role in determining the observed rate-dependent non-linear stress-strain response of the anterior cornea, which underlies the shape and optical properties of the tissue.

Shear wave elastography is emerging as a clinically valuable diagnostic tool to differentiate between benign and malignant breast masses. Elastography techniques assume that soft tissue can be modelled as a purely elastic medium. However, this assumption is often violated as soft tissue exhibits viscoelastic properties. In order to explore the role of viscoelastic parameters in suspicious breast masses, a study was conducted on a group of patients using shear wave dispersion ultrasound vibrometry in the frequency range of 50-400 Hz. A total of 43 female patients with suspicious breast masses were recruited before their scheduled biopsy. Of those, 15 patients did not meet the data selection criteria. Voigt model based shear elasticity showed a significantly (p = 7.88x10(-6)) higher median value for the 13 malignant masses (16.76±13.10 kPa) compared to 15 benign masses (1.40±1.12 kPa). Voigt model based shear viscosity was significantly different (p = 4.13x10(-5)) between malignant (8.22±3.36 Pa-s) and benign masses (2.83±1.47 Pa-s). Moreover, the estimated time constant from the Voigt model, which is dependent on both shear elasticity and viscosity, differed significantly (p = 6.13x10(-5)) between malignant (0.68±0.33 ms) and benign masses (3.05±1.95 ms). Results suggest that besides elasticity, viscosity based parameters like shear viscosity and time constant can also be used to differentiate between malignant and benign breast masses.

Actomyosin-based cortical flow is a fundamental engine for cellular morphogenesis. Cortical flows are generated by cross-linked networks of actin filaments and myosin motors, in which active stress produced by motor activity is opposed by passive resistance to network deformation. Continuous flow requires local remodeling through crosslink unbinding and and/or filament disassembly. But how local remodeling tunes stress production and dissipation, and how this in turn shapes long range flow, remains poorly understood. Here, we study a computational model for a cross-linked network with active motors based on minimal requirements for production and dissipation of contractile stress: Asymmetric filament compliance, spatial heterogeneity of motor activity, reversible cross-links and filament turnover. We characterize how the production and dissipation of network stress depend, individually, on cross-link dynamics and filament turnover, and how these dependencies combine to determine overall rates of cortical flow. Our analysis predicts that filament turnover is required to maintain active stress against external resistance and steady state flow in response to external stress. Steady state stress increases with filament lifetime up to a characteristic time τm, then decreases with lifetime above τm. Effective viscosity increases with filament lifetime up to a characteristic time τc, and then becomes independent of filament lifetime and sharply dependent on crosslink dynamics. These individual dependencies of active stress and effective viscosity define multiple regimes of steady state flow. In particular our model predicts that when filament lifetimes are shorter than both τc and τm, the dependencies of effective viscosity and steady state stress on filament turnover cancel one another, such that flow speed is insensitive to filament turnover, and shows a simple dependence on motor activity and crosslink dynamics. These results provide a framework for understanding how animal cells tune cortical flow through local control of network remodeling., Author summary In this paper, we develop and analyze a minimal model for a 2D network of cross-linked actin filaments and myosin motors, representing the cortical cytoskeleton of eukaryotic cells. We implement coarse-grained representations of force production by myosin motors and stress dissipation through an effective cross-link friction and filament turnover. We use this model to characterize how the sustained production of active stress, and the steady dissipation of elastic stress, depend individually on motor activity, effective cross-link friction and filament turnover. Then we combine these results to gain insights into how microscopic network parameters control steady state flow produced by asymmetric distributions of motor activity. Our results provide a framework for understanding how local modulation of microscopic interactions within contractile networks control macroscopic quantities like active stress and effective viscosity to control cortical deformation and flow at cellular scales.

Cell mechanics has proven to be important in many biological processes. Although there is a number of experimental techniques which allow us to study mechanical properties of cell, there is still a lack of understanding of the role each sub-cellular component plays during cell deformations. We present a new mesoscopic particle-based eukaryotic cell model which explicitly describes cell membrane, nucleus and cytoskeleton. We employ Dissipative Particle Dynamics (DPD) method that provides us with the unified framework for modeling of a cell and its interactions in the flow. Data from micropipette aspiration experiments were used to define model parameters. The model was validated using data from microfluidic experiments. The validated model was then applied to study the impact of the sub-cellular components on the cell viscoelastic response in micropipette aspiration and microfluidic experiments., Author summary Predictive simulations of cell flow in microfluidic devices and capillary networks may help to quantify the impact of different cell components on its behavior. Cells have complex mechanical properties and can undergo significant deformations, requiring detailed models to give an insight into the cell rheology. We developed computational model for simulations of cells with nucleus and cytoskeleton in flows in complex domains such as capillary networks and microfluidic devices. We validated the model using experimental data and used it to quantify the effects of cell components on its behavior. We envision that the proposed model will allow to study in silico numerous problems related to the cell biomechanics in flows.

An orb web's prey capture thread features tiny glue droplets, each formed of an adhesive glycoprotein core surrounded by an aqueous layer. Small molecules in the aqueous layer confer droplet hygroscopicity and maintain glycoprotein viscoelasticity, causing droplet volume and glycoprotein performance to track changes in environmental humidity. Droplet extension combines with that of a thread's supporting flagelliform fibers to sum the adhesive forces of multiple droplets, creating an effective adhesive system. We combined measurements of the force on an extending droplet, as gauged by the deflection of its support line, with measurements of glycoprotein volume and droplet extension to determine the Young's modulus (E) and toughness of three species' glycoproteins. We did this at five relative humidities between 20 - 90% to assess the effect of humidity on these properties. When droplets of a thread span extend, their extensions are constrained and their glycoprotein filaments remain covered by aqueous material. This was also the case during the first extension phase of the individual droplets that we examined. However, as extension progressed, the aqueous layer was progresses disrupted, exposing the glycoprotein. During the first extension phase E ranged from 0.00003 GPa, a value similar to that of fibronectin, a glycoprotein that anchors cells in the extracellular matrix, to 0.00292 GPa, a value similar to that of resilin in insect ligaments. Second phase E increased 4.7 - 19.4-fold. When compared at the same humidity the E of each species' glycoprotein was less than 5% of the value reported for its flagelliform fibers. This difference may facilitate the coordinated extension of these two capture thread components that is responsible for summing the thread's adhesive forces. National Science Foundation grant [IOS-1257719] National Science Foundation grant IOS-1257719 (BO) supported this study. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Thin liquid films are central to everyday life. They are ubiquitous in modern technology (pharmaceuticals, coatings), consumer products (foams, emulsions) and also serve vital biological functions (tear film of the eye, pulmonary surfactants in the lung). A common feature in all these examples is the presence of surface-active molecules at the air-liquid interface. Though they form only molecular-thin layers, these surfactants produce complex surface stresses on the free surface, which have important consequences for the dynamics and stability of the underlying thin liquid film. Here we conduct simple thinning experiments to explore the fundamental mechanisms that allow the surfactant molecules to slow the gravity-driven drainage of the underlying film. We present a simple model that works for both soluble and insoluble surfactant systems in the limit of negligible adsorption-desorption dynamics. We show that surfactants with finite surface rheology influence bulk flow through viscoelastic interfacial stresses, while surfactants with inviscid surfaces achieve stability through opposing surface-tension induced Marangoni flows.