Rigid Particles Research Articles

Pore structure characteristics and impact factors of laminated shale oil reservoir in Chang 73 sub-member of Ordos Basin, China

This research paper focuses on the laminated shale oil reservoir in the third sub-member of the seventh member of Yanchang Formation (Chang 73 sub-member) in the Ordos Basin. The study aims to comprehensively analyze the lithofacies type, micro qualitative and quantitative pore structure parameters, and the main controlling factors of the pore structure in laminated shale. The analysis involves various techniques, including a comprehensive analysis of TOC, rock-eval, X-ray diffraction, polarized light and fluorescence microscope observation, field emission scanning electron microscope observation, and low-pressure N2 adsorption analysis. Based on the sedimentation characteristics, geochemistry, and mineral composition differences, the Chang 73 laminated shale can be classified into three lithofacies types: tuffaceous-organic matter binary laminated shale, clayey-organic matter binary laminated shale, and felsic-clayey binary laminated shale. The pore network consists primarily of organic hydrocarbon generation pressurization fractures, clay mineral felsic intergranular composite pores, and felsic plasmid intergranular pore fracture systems. Mesopores are the most developed pore type. The pore volume and specific surface area increase in the order of “tuffaceous-organic matter”, “clayey-organic matter”, and “felsic-clay”, while the heterogeneity of the pore network gradually weakens, and the roughness of the pore surface enhances. The overall development of organic matter pores is limited, with organic matter-pyrite-clay mineral composite pores being the main components of micropores. The primary intergranular pore system, associated with rigid quartz particles, dominates the mesopores and macropores, and acts as the main framework of the entire pore network. The development of feldspar dissolution pores is limited and contributes minimally to the pore network.

Numerical investigation of the dynamics of a rigid spherical particle in a vortical cross-slot flow at moderate inertia

The study of flow and particle dynamics in microfluidic cross-slot channels is of high relevance for lab-on-a-chip applications. In this work, we investigate the dynamics of a rigid spherical particle in a cross-slot junction for a channel height-to-width ratio of 0.6 and at a Reynolds number of 120 for which a steady vortex exists in the junction area. Using an in-house immersed-boundary-lattice-Boltzmann code, we analyse the effect of the entry position of the particle in the junction and the particle size on the dynamics and trajectory shape of the particle. We find that the dynamics of the particle depend strongly on its lateral entry position in the junction and weakly on its vertical entry position; particles that enter close to the centre show trajectory oscillations. Larger particles have longer residence times in the junction and tend to oscillate less due to their confinement. Our work contributes to the understanding of particle dynamics in intersecting flows and enables the design of optimised geometries for cytometry and particle manipulation.

Numerical simulation study on erosion collapse failure mechanism of cemented conglomerate accumulation

PurposeCemented conglomerate accumulation is a weak and heterogeneous medium that occurs in western China. It consists mainly of argillaceous cement that loses strength rapidly upon contact with water, leading to collapse instability failure. Its deformation failure mechanism is complex and poorly understood. In this paper, the erosion failure mechanism of cemented conglomerate accumulation is investigated.Design/methodology/approachThe collapse failure process after erosion of the slope foot for typical cemented conglomerate accumulation is studied based on field investigation using the particle discrete element method. And how the medium composition, slope angle and cementation degree influence the failure mode and process of the cemented conglomerate accumulation is examined.FindingsThe foot erosion of slope induces a tensile failure that typically manifests as “erosion at the foot of slope – tensile cracking at the back edge of slope top – integral collapse.” The collapse failure is more likely to occur when the cemented conglomerate accumulation has a higher rock content, a steeper slope angle or a weaker cementation degree.Originality/valueA model based on rigid blocks and disk particles to simulate the cemented conglomerate accumulation is developed. It shows that the hydraulic erosion at the foot of the slope resulted in a different failure mechanism than that of general slopes. The results can inform the stability management, disaster prevention and mitigation of similar slopes.

Numerical estimation via remeshing and analytical modeling of nonlinear elastic composites comprising a large volume fraction of randomly distributed spherical particles or voids

We investigate numerically via finite element (FE) simulations and analytically the nonlinear behavior of a soft matrix comprising a large volume fraction (up to 55%) of monodisperse spherical inclusions, which can be either rigid particles or voids. To address the issue of severe mesh distortion at large strains, we employ a successive remeshing and solution mapping technique that allows us to achieve a large macroscopic deformation, with the maximum attained stretch being at least two to three times when compared to that without remeshing. A general algorithm allowing to read complex arbitrary particle geometries is proposed and implemented allowing to remesh a finitely strained domain with randomly distributed inclusions of arbitrary shape. The present simulation results for a neo-Hookean elastic matrix filled with rigid particles show that the nonlinear stress–stretch uniaxial tension response can be well approximated by a neo-Hookean material with an effective shear modulus that depends on the volume fraction of the inclusions. Additionally, a study of the number of particles demonstrates that sixty four monodisperse spheres is a good compromise to reach accurate results for large deformations and maintain a reasonable CPU-time. The FE data serve as a powerful assessment tool for analytical homogenization models. Newly developed and existing models for small and large strains are proposed and assessed, respectively, for both rigid particles and voids. These results fill the gap for large volume fractions of inclusions and serve as a reference for the homogenization of the microstructures of interest.

Electrophoresis of a dielectric droplet with constant surface charge density.

Electrophoresis of a dielectric fluid droplet with constant surface charge density is investigated theoretically in this study. A pseudo-spectral method based on Chebyshev polynomials is adopted to solve the governing electrokinetic equations. It is found, among other things, that the larger the electrolyte strength in the ambient solution is, the slower the droplet moves in general. This is due to the strong screening effect of the large amount of indifferent counterions in the neighborhood of the droplet, with no reinforcement of potential-determining ions adsorbing to the droplet surface. The droplet comes to a complete halt eventually. Critical points are discovered for highly charged droplets, at which the droplet surface becomes immobile and the interior fluid stops recirculating. The droplet moves like a rigid particle with constant mobility regardless of its viscosity, a situation referred to as the "solidification phenomenon." The deadlock between the spinning motions on the charged droplet surface induced by the electric driving force and the hydrodynamic driving force respectively is responsible for this peculiar phenomenon. This is also observed for a dielectric droplet with constant surface electric potential. We demonstrate here that it occurs in the constant surface charge density situation as well.

Collision Statistics of Droplets in Turbulence Considering Lubrication Interactions, Mobility of Interfaces, and Non-continuum Molecular Effects

Collision statistics of same-size aerodynamically interacting droplets in homogeneous isotropic turbulence is examined via hybrid direct numerical simulations combined with Lagrangian particle tracking under a point-particle assumption and a one-way momentum coupling. The simulations are performed both in the absence and presence of gravity for various droplet Stokes numbers (0.1–2) over a wide range of liquid water contents (LWCs = 1–30 g/m3\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\ ext {g/m}^3$$\\end{document}). Also, the effects of using different representations of the lubrication forces, i.e. a continuum and a non-continuum one, have been investigated. The results have highlighted the importance of considering the aerodynamic interaction, especially in the presence of gravity for the systems that have high LWCs. Likewise, taking the lubrication force into account shows a notable change in statistics, but a non-continuum representation yields results that are close to the continuum one. In addition, the impact of modeling droplets as fluid drops instead of rigid particles has been examined. Obtaining quite close statistics under both cases demonstrates that a rigid particle assumption for water droplets in air is sufficiently accurate owing to the high water-to-air viscosity ratio. Moreover, the collision efficiency is shown to be in the range 60–100% in the absence of gravity, which approaches 100% as the LWC is enlarged. In the presence of gravity, however, the efficiency grows, with both the Stokes number and the LWC, to values as high as 300%. Finally, for settling droplet systems, the enhancement factor due to turbulence is quantified, exhibiting a growth with the LWC and a reduction with the Stokes number.

Study of red blood cells and particles in stenosed microvessels using coupled discrete and continuous forcing immersed boundary methods

A computational study of the blood flow in a stenosed microvessel is presented using coupled discrete ghost-cell and continuous-forcing immersed boundary methods. This study focuses on studying platelet behaviors near the stenosis with deformable red blood cells (RBCs). The influence of varying hematocrit, area blockage, stenosis shape, and driving force on flow characteristics, RBCs, and particle behaviors is considered. Distinct flow characteristics are observed in stenosed microvessels in the presence of RBCs. The motion of RBCs is the major cause of time-dependent oscillations in flow rates, while the contribution of particles to the fluctuations is negligible. However, this effect decreases when the stenosis is elongated in the axial direction. Interestingly, as the hematocrit level increases, downstream particles move closer to the vessel wall due to the enhanced shear-induced lift force resulting from the interaction among RBCs and particles. Furthermore, it is observed that geometrical changes in the stenosis have a more significant impact on the axial profile of particle concentration compared to changes in hematocrit or driving force. An asymmetric stenosis leads to asymmetric profiles in the flow velocity and the distribution of cells and particles due to the geometric focusing effect of the stenosis. There is no significant change in flow rates until a blockage of 0%–50%, but a sudden increase in the root mean square of flow rates occurs at an 80% blockage. This study contributes to our understanding of the rheological behaviors of RBCs and rigid particles in a stenosed microvessel under various hemodynamic conditions.

Internal Friction Angle of Cohesionless Binary Mixture Sand–Granular Rubber Using Experimental Study and Machine Learning

This study aimed to examine the shear strength characteristics of sand–granular rubber mixtures in direct shear tests. Two different sizes of rubber and one of sand were used in the experiment, with the sand being mixed with various percentages of rubber (0%, 10%, 20%, 30%, and 50%). The mixtures were prepared at three different densities (loose, slightly dense, and dense), and shear stress was tested at four normal stresses (30, 55, 105, and 200 kPa). The results of 80 direct shear tests were used to calculate the peak and residual internal friction angles of the mixtures, and it was found that the normal stress had a significant effect on the internal friction angle, with an increase in normal stress leading to a decrease in the internal friction angle. These results indicated that the Mohr–Coulomb theory, which applies to rigid particles only, is not applicable in sand–rubber mixtures, where stiff particles (sand) and soft particles (rubber) are mixed. The shear strength of the mixtures was also influenced by multiple factors, including particle morphology (size ratio, shape, and gradation), mixture density, and normal stress. For the first time in the literature, genetic programming, classification and regression random forests, and multiple linear regression were used to predict the peak and residual internal friction angles. The genetic programming resulted in the creation of two new equations based on mixture unit weight, normal stress, and rubber content. Both artificial intelligence models were found to be capable of accurately predicting the peak and residual internal friction angles of sand–rubber mixtures.

Derivation of the viscoelastic stress in Stokes flows induced by nonspherical Brownian rigid particles through homogenization

Derivation of the viscoelastic stress in Stokes flows induced by nonspherical Brownian rigid particles through homogenization

Incorporation of Large-Scale Prepared Poly(ethylene oxalate) into Biodegradable Poly(butylene adipate-co-terephthalate) Blown Films with Enhanced Mechanical and Barrier Performance

With the breakthrough of large-scale production of poly(ethylene oxalate) (PEOx) engineering technology, the promising sustainable polyoxalate material with high mechanical strength and marine degradability processes the feasibility of commercialization. By blending PEOx produced in 150 L pilot polymerization systems with commercial poly(butylene adipate-co-terephthalate) (PBAT) polyester, followed by a film-blowing process, the compatibility, mechanical properties, barrier properties, particularly the marine degradability, were systematically studied. Results revealed that adding the appropriate amount of PEOx as rigid particles could significantly enhance the comprehensive performance of the PBAT/PEOx blend films. More appealingly, the preferential hydrolysis of the PEOx component caused the surface of the blend film to break more easily in the marine environments, creating conditions for further degradation of the blend matrix. Additionally, the subsequent bag-making and food preservation process also proved that the film material is a potential candidate for broader commercial applications. Accordingly, for commercial applications of plastics, the PBAT/PEOx blend film is a valuable marine-degradable material in terms of environmental and economic sustainability.

Pore types and microstructure of the black shale in the Lower Silurian, Southeast of Chongqing, China

The Lower Silurian black shale in the Southeast of Chongqing is of great potential in the accumulation of shale gas in South China. In order to understand the pore system of these black shales, mercury porosimetry, surface area study, mineralogy and image analyses by scanning electronic microscopy were performed. Four major pore types (interparticle pores, intraparticle pores, organic-matter pores and microfractures) were classified under the observation of the samples. Interparticle pores, which associate with mineral particles and mostly in micro-scale, can be subdivided into pores among rigid particles, ductile particles, and rigid-ductile particles. Intraparticle pores, which can be subdivided into cleavage-plane pores within clay minerals, intercrystalline pores within pyrite framboids and pores formed by partial or complete dissolution, are in nanometer scale and more common exist in unweathered rocks. Organic-matter pores are one type of intraparticle pores which are found within organic matter and appear to be prone to form effective networks and main permeable pathway. The rocks are abundant with microfractures where absorbed gas can be transformed to free gas. Microfractures which are partly cemented by residual organic residue are of great importance. Comparing the two sets of samples from Lujiao section and Well Yuye-1, we note that the outcrop samples from Lujiao section have much more interparticle pores than core samples, while the latter own more intraparticle pores. The results of quantitative and visual qualitative analyses about the gas shale reservoir pore system are similar. With the increasing of quartz, the number of interparticle pores and macropores increase. And intraparticle pores and micro-mesopores increase with clay minerals.

Core-shell ZrO2@GO hybrid for effective interfacial adhesion improvement of carbon fiber/epoxy composites

A core-shell structured ZrO2@GO hybrid was fabricated and incorporated onto carbon fiber via dip-coating for the first time for effectively enhancing the interfacial adhesion of epoxy composites. ZrO2@GO functionalized fiber composite exhibits the optimal interfacial shear strength, with 47.5%, 21.8% and 16.8% augmentations over that of untreated fiber, only epoxy sizing coated fiber and silane coupling agent-encapsulated ZrO2 decorated fiber composites, suggesting that GO could serve as a new coupling agent for enhancing the affinity of inorganic nanoparticles with resin component instead of the traditional silane coupling agent. The reinforcing mechanisms are ascribed to the improved resin wettability, abundant functional groups offered by GO for increasing multiple interactions with epoxy resin and the formation of a mobile network absorbing energy. It is believed that the core-shell ZrO2@GO structure, with inner rigid particles and outer flexible coupling agent, can fully utilize the synergistic enhancement effect of different components and provide new insights into designing nano-structured interphase for the improvement in interfacial adhesion of fiber-polymer composites.

Effect of Spent Coffee Grounds on the Crystallinity and Viscoelastic Behavior of Polylactic Acid Composites.

This work investigated the addition of spent coffee grounds (SCG) as a valuable resource to produce biocomposites based on polylactic acid (PLA). PLA has a positive biodegradation effect but generates poor proprieties, depending on its molecular structure. The PLA and SCG (0, 10, 20 and 30 wt.%) were mixed via twin-screw extrusion and molded by compression to determine the effect of composition on several properties, including mechanical (impact strength), physical (density and porosity), thermal (crystallinity and transition temperature) and rheological (melt and solid state). The PLA crystallinity was found to increase after processing and filler addition (34-70% in the 1st heating) due to a heterogeneous nucleation effect, leading to composites with lower glass transition temperature (1-3 °C) and higher stiffness (~15%). Moreover, the composites had lower density (1.29, 1.24 and 1.16 g/cm3) and toughness (30.2, 26.8 and 19.2 J/m) as the filler content increased, which is associated with the presence of rigid particles and residual extractives from SCG. In the melt state, polymeric chain mobility was enhanced, and composites with a higher filler content became less viscous. Overall, the composite with 20 wt.% SCG provided the most balanced properties being similar to or better than neat PLA but at a lower cost. This composite could be applied not only to replace conventional PLA products, such as packaging and 3D printing, but also to other applications requiring lower density and higher stiffness.

Discrete element modeling of granular hopper flow of irregular-shaped deformable particles

Many natural and engineered granular materials have relatively deformable particles. Besides particle size and shape, particle deformability is another salient factor that significantly impacts the material’s flow behavior. In this work, the flow of irregular-shaped deformable particles in a wedge-shaped hopper is investigated using discrete element simulations. A bonded-sphere model is developed to simultaneously capture irregular particle shapes and particle-wise deformations (e.g., compression, deflection, and distortion). Quantitative analysis of the effects of irregular shapes and particle deformations shows that the increase in particle stiffness tends to increase initial packing porosity and decrease the flow rate in the hopper. Rigid particles tend to have clogging issues, whereas deformable particles have less chance to, indicating particle deformation reduces the critical bridging width in the hopper flow. Detailed analysis of stress fields is also conducted to provide insights into the mechanism of particle flow and clogging. Stresses and discharge rates calculated from numerical simulations are compared and show good agreement with Walker’s theory and the extended Beverloo formula. Simulations with various particle shape combinations are also performed and show that the initial packing porosity decreases with an increasing percentage of fibers while the discharge rate has a complex dependency on particle shapes.

Discrete element modelling of the effect of aspect ratio on compaction and shear behaviour of aggregates

Coarse aggregates are essential for constructing railways and highways, as they affect structural performance through their compacted state and shear strength. Particle geometry influences these properties, but little research has quantitatively investigated how shape indexes impact mechanical indicators. This study addresses this gap by exploring the effect of particle aspect ratio on coarse aggregate compaction and shear behavior using the discrete element method (DEM). The study's key contribution is the development of novel equations that describe the correlation between packing porosity, inter-particle contact state, and aspect ratio control thresholds in a practical format. To do so, first, a point cloud database is compiled from laser scans of 700 natural aggregate particles. Each particle is classified using three aspect-ratio parameters: elongation, flatness, and flatness-elongation index (FEI). They are then digitized as rigid block particle models and compared against actual point clouds. Next, nine aggregates of representative shapes are assessed for their macro and mesoscopic mechanical responses using DEM-based vibro-compaction and direct shear tests. The compaction simulations indicate an exponential relationship between porosity and inter-particle coefficient of friction and show that aggregates with FEI ≈ 0.5 are a good choice. The direct shear simulations consider the same initial porosity, and show that the relationship between peak shear strength and vertical stress follows a power function. Finally, the analysis of contact forces and fabric tensors implies that FEI ≈ 0.5 may also facilitate particle fragmentation control.

Discrete-element method simulations of the seismic response of tunnels in deep granular deposit

The seismic response of tunnels constructed in a deposit of granular soil of depth greater than 20 m is examined in this paper. Three-dimensional discrete-element method (DEM) simulations were used to model the soil assembly and the tunnel lining. Soil particles were simulated as rigid spherical particles that are allowed to overlap at the contact points. The tunnel lining was modeled as a flexible system using a collection of particles that were glued together using cementitious bonds to mimic the physical properties of an actual tunnel. Free-field boundaries were utilized at the model’s lateral sides to prevent the propagating wave’s reflections back to the domain and enforce the free-field motion. The DEM computational framework was validated using published centrifuge test results. DEM simulation results reveal good agreement with the experimental centrifuge test results when subjected to the same input motion. Further simulations were carried out to investigate the seismic response of the system under different input motion characteristics. The dynamic earth pressure was monitored at certain points on the tunnel lining. It was found that the dynamic earth pressure and tunnel lining forces increased during seismic excitation and the final residual values were, in most cases, considerably larger than the initial static values. The results also showed that the maximum amplification factor, ground settlement and tunnel lining internal forces occurred close to the resonant frequency of the system.

Assessing Differential Particle Deformability under Microfluidic Flow Conditions.

Assessing the mechanical behavior of nano- and micron-scale particles with complex shapes is fundamental in drug delivery. Although different techniques are available to quantify the bulk stiffness in static conditions, there is still uncertainty in assessing particle deformability in dynamic conditions. Here, a microfluidic chip is designed, engineered, and validated as a platform to assess the mechanical behavior of fluid-borne particles. Specifically, potassium hydroxide (KOH) wet etching was used to realize a channel incorporating a series of micropillars (filtering modules) with different geometries and openings, acting as microfilters in the direction of the flow. These filtering modules were designed with progressively decreasing openings, ranging in size from about 5 down to 1 μm. Discoidal polymeric nanoconstructs (DPNs), with a diameter of 5.5 μm and a height of 400 nm, were realized with different poly(lactic-co-glycolic acid) (PLGA) and poly(ethylene glycol) (PEG) ratios (PLGA/PEG), namely, 5:1 and 1:0, resulting in soft and rigid particles, respectively. Given the peculiar geometry of DPNs, the channel height was kept to 5 μm to limit particle tumbling or flipping along the flow. After thorough physicochemical and morphological characterization, DPNs were tested within the microfluidic chip to investigate their behavior under flow. As expected, most rigid DPNs were trapped in the first series of pillars, whereas soft DPNs were observed to cross multiple filtering modules and reach the micropillars with the smallest opening (1 μm). This experimental evidence was also supported by computational tools, where DPNs were modeled as a network of springs and beads immersed in a Newtonian fluid using the smoothed particle hydrodynamics (SPH) method. This preliminary study presents a combined experimental-computational framework to quantify, compare, and analyze the characteristics of particles having complex geometrical and mechanical attributes under flow conditions.

Modeling and Experiments of Droplet Evaporation with Micro or Nano Particles in Coffee Ring or Coffee Splat.

Experimental and numerical experiments were carried out to study the coffee rings or coffee splats formed by droplet evaporation with micro or nano polystyrene sphere particles (Dp = 10 μm or 100 nm). Particle image velocimetry (PIV) and a high-resolution camera were used in this experiment, along with a temperature-controlled heater and a data-acquisition computer. The results showed that a nano particle could form a homogeneous coffee splat, instead of the common coffee ring formed when using micro particles. In order to account for this phenomenon, this paper developed a complex multiphase model, one which included the smooth particle hydrodynamics (SPH) fluid model coupled with the van der Waals equation of state for droplet evaporation, the rigid particle model of finite-size micro particles, and the point-particle model of the nanometer particles. The numerical simulation was operated on a GPU-based algorithm and tested by four validation cases. A GPU could calculate 533 times the speed of a single-core CPU for about 300,000 particles. The results showed that, for rigid solid particles, the forms emerged spontaneously on the wall, and their structure was mainly affected by the boundary wettability, and less affected by the fluid flow and thermal condition. When the wall temperature was low, it was easier for the particles to be deposited on the contact line. At high wall temperature, the coffee ring effect would be weakened, and the particles were more likely to be deposited in the droplet center. The hydrophilic surface produced a larger coffee ring compared to the hydrophobic surface. The experimental and numerical results proved that particle size could play a significant role during the particle deposition, which may be a possible route for producing uniform-distributed and nano-structure coatings.

Development of a toughened and antibacterial Poly(lactide acid) (PLA) with preserved strength by elemental sulfur-based bio-renewable dynamically crosslinked elastomers

The toughening of PLA is usually realized by blending with elastomers or plasticizers, but often conflicted with lowered mechanical strength and complicated preparation steps. In the current study, inverse vulcanization was used to synthesize a series of dynamically cross-linked polysulfide (SOA) elastomers based on refining by-product elemental sulfur, soybean oil (SO) and allyl glyceryl ether (AGE), which provides dynamically reversible S-S bonds, flexible triglyceride segments and reactive epoxy functions.Melt blended with PLA, the SOAs exhibit great toughening effects and produce droplet-matrix typed blends with good interfacial compatibilization verified by rheology. SOA composition was decisive in dictating toughening efficiency. A high AGE content induces self-polymerization of epoxy groups, resulting in dominance of permanent crosslinking in SOA, giving rise to rigid SOA particles with bad processability and dispersity in PLA, hence poor toughener cavitations under external forces with unsatisfying toughening effect. A high sulfur and low AGE content results in dominance of dynamic crosslinking, providing soft SOA particles in favor of processing and dispersion in PLA matrix with small domains and thin ligament thickness, thus favoring internal cavitation of SOA particles to achieve good toughness. However, the presence of unreacted sulfur with a too high sulfur content can mitigate the toughening effect. In particular, a highly toughened blend with optimized SOA composition (PLA-S50O25A25) was obtained with excellent toughness (elongation at break 235.7% and tensile toughness 76.5 MJ/m3) balanced by only slightly impaired strength (tensile strength 50–60 MPa and modulus ∼ 900 Mpa). Moreover, the PLA-SOA exhibits excellent broad-spectrum anti-bacterial effects. This study reported the unprecedented use of fully sustainable dynamically-crosslinked elastomers as tougheners of PLA, paving new ways for PLA modification with balanced stiffness-toughness and antibacterial activity.

Flow field disturbance due to point viscosity variations in a heterogeneous fluid

A fluid can be heterogeneous in either its density or viscosity - its two basic material properties. In this Letter, we ask how an ambient flow field is disturbed due to a point viscosity variation in the fluid. If we are able to answer this question, we can create an arbitrary spatial distribution of viscosity discretely by placing viscosity sources and sinks of appropriate strength, allowing us to model complex fluids. As an example, we show how the hydrodynamics of a rigid particle is altered due to the presence of point viscosity variations.