Fibrin network adaptation to cell-generated forces

Fibrin promotes wound healing by serving as provisional extracellular matrix for fibroblasts that realign and degrade fibrin fibers, and sense and respond to surrounding substrate in a mechanical-feedback loop. We aimed to study mechanical adaptation of fibrin networks due to cell-generated forces at the micron-scale. Fibroblasts were elongated-shaped in networks with ≤?2 mg/ml fibrinogen, or cobblestone-shaped with 3 mg/ml fibrinogen at 24 h. At frequencies f?<?10² Hz, G′ of fibroblast-seeded fibrin networks with ≥?1 mg/ml fibrinogen increased compared to that of fibrin networks. At frequencies f?>?10³ Hz, G″ of fibrin networks decreased with increasing concentration following the power-law in frequency with exponents ranging from 0.75?±?0.03 to 0.43?±?0.03 at 3 h, and of fibroblast-seeded fibrin networks with exponents ranging from 0.56?±?0.08 to 0.28?±?0.06. In conclusion, fibroblasts actively contributed to a change in viscoelastic properties of fibrin networks at the micron-scale, suggesting that the cells and fibrin network mechanically interact. This provides better understanding of, e.g., cellular migration in wound healing.

Open image in new window

Graphical abstract

     

A bootstrap mechanism for non-colloidal suspension viscosity

The role of friction in non-colloidal suspensions is examined with a model which splits the viscosity into a frictionless component (τ*) plus a frictional component which depends on the ratio of the particle pressure (P) to the shear stress (τ). The model needs the input by computation of τ* and P and a suitable choice of particle friction coefficient (μ). It can be extended to elongational flows and cases where sphere roughness is important; volume fractions up to 0.5 are considered. It is shown that friction acts in a feedback or “bootstrap” manner to increase the suspension viscosity. The analysis is also useful for deducing the friction coefficient in suspensions from experimental data. It was applied to several sets of experimental data and reasonable correlations of the viscosities were demonstrated. An example of the correlation for spheres in a silicone oil is shown for volume fractions 0.1–0.5.

Open image in new window

Graphical abstract

     

Re-examination of terminal relaxation behavior of high-molecular-weight ring polystyrene melts

For high-molecular-weight (M) ring polymers with low contamination of linear chains, recent viscoelastic tests revealed broad terminal relaxation associated with no clear entanglement plateau. This relaxation behavior is qualitatively similar to that deduced from molecular models (double-folded lattice-animal model and the fractal loopy globule model) for entangled ring polymers, but quantitatively important differences are also noted: For example, the full terminal relaxation of those polymers is slower than the model prediction. This study re-examined the viscoelastic data of entangled high-M ring polystyrene (PS) samples (coded as R-240; M = 244×10³) specifically for two points: the purity of the ring samples after the viscoelastic tests and the molecular origin of the stress. For the first point, the R-240 samples contaminated with linear chains at low but different levels were prepared by tuning either the purification efficiency or the retention time of the sample at high temperature (T) before/during the viscoelastic test. The fraction w _L of the linear contaminant, determined after the viscoelastic measurement, was ranging from 0.7 to 4.9%, and the extrapolation of the modulus data to w _L = 0 gave the data for the ideally pure ring melt. This pure ring melt exhibited broad terminal relaxation that started faster but completed slower compared to the model prediction, indicating that the ring relaxation is not well described by the current model(s) even in the absence of linear contaminant. For the second point, dynamic birefringence measurements were conducted for the R-240 samples with w _L = 4.6 and 1.0%. These samples obeyed the stress-optical rule, and their stress-optical coefficient was indistinguishable from that for linear PS samples, revealing that the stress of the ring PS chains reflects the orientational anisotropy of the chains (as is the case also for linear chains). The relaxation behavior of pure ring PS melt is discussed on the basis of these findings, with the focus being placed on the ring-ring threading not considered in the models.     

Inviscid Models Generalizing the Two-dimensional Euler and the Surface Quasi-geostrophic Equations

Any classical solution of the two-dimensional incompressible Euler equation is global in time. However, it remains an outstanding open problem whether classical solutions of the surface quasi-geostrophic (SQG) equation preserve their regularity for all time. This paper studies solutions of a family of active scalar equations in which each component u _j of the velocity field u is determined by the scalar θ through \({u_j =\mathcal{R}\Lambda^{-1}P(\Lambda) \theta}\) , where \({\mathcal{R}}\) is a Riesz transform and Λ = (?Δ)^1/2. The two-dimensional Euler vorticity equation corresponds to the special case P(Λ) = I while the SQG equation corresponds to the case P(Λ) = Λ. We develop tools to bound \({\|\nabla u||_{L^\infty}}\) for a general class of operators P and establish the global regularity for the Loglog-Euler equation for which P(Λ) = (log(I + log(I ? Δ))) ^γ with 0 ≦ γ ≦ 1. In addition, a regularity criterion for the model corresponding to P(Λ) = Λ ^β with 0 ≦ β ≦ 1 is also obtained.     

Rheological and micro structural studies of methylcellulose solutions in <Emphasis Type="Italic">N</Emphasis>,<Emphasis Type="Italic">N</Emphasis>-dimethylformamide for preparation of extruded beads: controlled urea release

The rheological properties of methylcellulose in N,N-dimethylformamide (MC-DMF) gel are investigated to prepare extruded beads. The temperature scan under dynamic compression for various concentrations of MC in DMF is performed to investigate the rapture of MC gel at a constant frequency of 1 Hz. The morphological studies are performed using a scanning electron microscope (SEM) to analyze the size and shape of dried bead. However, during swelling studies, the MC beads have the capability to swell and retain a large amount of water >?9150% by weight and 9192.6% by volume. The mechanism of swelling is thermodynamically verified, where the enthalpy of hydration of initial layer of MC bead is negative. The newly defined electrostatic penta-pole model explains the anomalous behavior of urea release, where urea is assumed to be electrostatically bounded with the MC molecules.

Open image in new window

Grapichal abstract ?

     

Rheological behavior of composites made from linear medium-density polyethylene and hemp fibers treated by surface-initiated catalytic polymerization

This work presents different rheological methods to determine the effect of fiber surface treatment on their interaction with a polymer matrix. In particular, surface-initiated catalytic polymerization was investigated on hemp fibers to improve their adhesion with linear medium-density polyethylene (LMDPE). The selected rheological tests (creep-recovery (solid state), small and large amplitude oscillation shear, and transient rheology (melt state)) were used to compare the treated and untreated fiber composites with the neat matrix. The results showed a significant improvement of the treated hemp composite (LPHC) creep modulus with respect to its untreated counterpart (LNHC) leading to a reduction of the creep strain, especially as temperature increases. The transient viscosity was modeled using a modified Kohlrausch-Williams-Watt (KWW) equation showing an increase in the transient viscosity (\( {\eta}_0^{+} \)) and relaxation time (τ) with fiber addition and surface treatment. These results were confirmed by large amplitude oscillatory shear (LAOS) through the reduction of the relative third harmonic (I_3/1), intrinsic nonlinearity parameter (Q₀), and nonlinear viscoelastic ratio (NRL). The results clearly show that catalytic polymerization is a good surface modification technique to increase the compatibility between natural fibers and polymer matrices as to improve all their final properties.

Open image in new window

Graphical Abstract ?

     

Comparative study of the electrorheological effect in suspensions of needle-like and isotropic cerium dioxide nanoparticles

A novel approach to the analysis of the electrorheological effect is proposed, based on the expansion of dimensionless relative shear stress as function of electric field strength in the power series \( {\tau}_{\mathrm{rel}}=\frac{\tau_E}{\tau }=1+\frac{\alpha }{\tau }E+\frac{\beta }{\tau }{E}^n \). The application of this approach to investigation of the electrorheological effect in suspensions of isotropic and needle-like CeO₂ nanoparticles in polydimethylsiloxane has revealed that the polynomial coefficients can be judged as a measure of the efficiency of transformation of electrical energy into mechanical energy. The values of α and β coefficients depend on the shape and concentration of filler particles, as well as on the shear rate. The value and the sign of these coefficients determine both the magnitude of the electrorheological effect and the type of dependence of the shear stress (linear or power law) on the strength of the electric field. It has been shown that the values of α and β coefficients for the electrorheological fluids with needle-like particles are greater than for fluids with isotropic particles (at the same concentration of suspensions), which is associated with the different polarization of particles in the applied electric field.

Open image in new window

Graphical abstract A novel approach to the analysis of the electrorheological effect is proposed.

     

Mechanical Properties of Hard W-C Coating on Steel Substrate Deduced from Nanoindentation and Finite Element Modeling

Mechanical properties of a hard and stiff W-C coating on steel substrate have been investigated using nanoindentation combined with finite element modeling (FEM) and extended FEM (XFEM). The significant pile-up observed around the indents in steel substrate caused an overestimation of hardness and indentation modulus. A simple geometrical model, considering the additional contact surfaces due to pile-up, has been proposed to reduce this overestimation. The presence of W-C coating suppressed the pile-up in the steel substrate and a transition to sink-in behavior occurred. The FEM simulations adequately reproduced the surface topography of the indents in the substrate and coating/substrate systems as well. The maximum principal stresses of the indented W-C/steel coated system were tensile; they were always located in the coating and evolved in 3 stages. Cohesive cracking occurred during loading in the sink-in zone (stage III) when the ultimate tensile strength (σ _max ) of the coating was reached. The obtained hardness (H _c ), indentation modulus (E _c ), yield stress (Y) and strength (σ _max ) of the W-C coating were H _c ? =?20 GPa, E _c ? =?250 GPa, Y?=?9.0 GPa and σ _max ? =?9.35 GPa, respectively. XFEM resulted in fracture energy of the W-C coating of G?=?38.1 J?·?m^-2 and fracture toughness of K _IC ? =?3.5 MPa?·?m^0.5.     

Linearized Stability for Semiflows Generated by a Class of Neutral Equations,with Applications to State-Dependent Delays

We prove a principle of linearized stability for semiflows generated by neutral functional differential equations of the form x′(t) = g(? x _t , x _t ). The state space is a closed subset in a manifold of C ²-functions. Applications include equations with state-dependent delay, as for example x′(t) = a x′(t + d(x(t))) + f (x(t + r(x(t)))) with \({a\in\mathbb{R}, d:\mathbb{R}\to(-h,0), f:\mathbb{R}\to\mathbb{R}, r:\mathbb{R}\to[-h,0]}\).     

An iterative updating method for undamped structural systems

Finite element model updating is a procedure to minimize the differences between analytical and experimental results and can be mathematically reduced to solving the following problem. Problem P: Let M _a ∈SR ^n×n and K _a ∈SR ^n×n be the analytical mass and stiffness matrices and Λ=diag{λ ₁,…,λ _p }∈R ^p×p and X=[x ₁,…,x _p ]∈R ^n×p be the measured eigenvalue and eigenvector matrices, respectively. Find \((\hat{M}, \hat{K}) \in \mathcal{S}_{MK}\) such that \(\| \hat{M}-M_{a} \|^{2}+\| \hat{K}-K_{a}\|^{2}= \min_{(M,K) \in {\mathcal{S}}_{MK}} (\| M-M_{a} \|^{2}+\|K-K_{a}\|^{2})\), where \(\mathcal{S}_{MK}=\{(M,K)| X^{T}MX=I_{p}, MX \varLambda=K X \}\) and ∥?∥ is the Frobenius norm. This paper presents an iterative method to solve Problem P. By the method, the optimal approximation solution \((\hat{M}, \hat{K})\) of Problem P can be obtained within finite iteration steps in the absence of roundoff errors by choosing a special kind of initial matrix pair. A numerical example shows that the introduced iterative algorithm is quite efficient.     

Regularity and Singularities of Optimal Convex Shapes in the Plane

We focus here on the analysis of the regularity or singularity of solutions Ω ₀ to shape optimization problems among convex planar sets, namely:

$J(\Omega_{0})={\rm min} \{J(\Omega), \Omega \quad {\rm convex},\Omega \in \mathcal{S}_{\rm ad}\},$

where \({\mathcal{S}_{\rm ad}}\) is a set of 2-dimensional admissible shapes and \({J:\mathcal{S}_{\rm ad}\rightarrow\mathbb{R}}\) is a shape functional. Our main goal is to obtain qualitative properties of these optimal shapes by using first and second order optimality conditions, including the infinite dimensional Lagrange multiplier due to the convexity constraint. We prove two types of results:

1. i)
   under a suitable convexity property of the functional J, we prove that Ω ₀ is a W ^2,p -set, \({p\in[1, \infty]}\). This result applies, for instance, with p = ∞ when the shape functional can be written as J(Ω) = R(Ω) + P(Ω), where R(Ω) = F(|Ω|, E _f (Ω), λ₁(Ω)) involves the area |Ω|, the Dirichlet energy E _f (Ω) or the first eigenvalue of the Laplace–Dirichlet operator λ₁(Ω), and P(Ω) is the perimeter of Ω;
    

1. ii)
   under a suitable concavity assumption on the functional J, we prove that Ω ₀ is a polygon. This result applies, for instance, when the functional is now written as J(Ω) = R(Ω) ? P(Ω), with the same notations as above.
    

     

Measuring and assessing first and second normal stress differences of polymeric fluids with a modular cone-partitioned plate geometry

We propose a simple, robust method to measure both the first and second normal stress differences of polymers, hence obtaining the full set of viscometric material functions in nonlinear shear flow. The method is based on the use of a modular cone-partitioned plate (CPP) setup with two different diameters of the inner plate, mounted on a rotational strain-controlled rheometer. The use of CPP allows extending the measured range of shear rates without edge fracture problems. The main advantage of such a protocol is that it overcomes limitations of previous approaches based on CPP (moderate temperatures not exceeding 120 °C, multiple measurements of samples with different volume) and yields data over a wide temperature range by performing a two-step measurement on two different samples with the same volume. The method was tested with two entangled polystyrene solutions at elevated temperatures, and the results were favorably compared with both the limited literature data on the second normal stress difference and the predictions obtained with a recent tube-based model of entangled polymers accounting for shear flow-induced molecular tumbling. Limitations and possible improvements of the proposed simple experimental protocol are also discussed.

Open image in new window

Graphical abstract The effects of edge fracture in start-up shear experiments can be circumvented with the use of a cone-partitioned plate (CPP) geometry. Such a device consists of an inner measuring plate surrounded by an outer nonmeasuring corona. The radius of the sample exceeds that of the measuring plate so that the measured volume is not affected by edge instability. However, the measured first normal stress difference is an apparent one (N_app,1), owing to the contribution of the nonmeasured part of the sample. The figure depicts a schematic design of a modular CPP geometry. Such a fixture is built in a way that the inner tool and the outer partition can be easily replaced, in order to have different measuring diameters (i.e., 6 and 10 mm). From the corresponding signals of the N_app,1, the effective first and second normal stress differences can be calculated.

     

Evaluation of the Spectral Element Dynamic Model for Large-Eddy Simulation on Unstructured,Deformed Meshes

In this paper, the Spectral-Element Dynamic Model (SEDM), suited for Large-Eddy Simulation (LES) using Discontinuous Finite Element Methods (DFEM), is assessed using unstructured meshes. Five test cases of increasing complexity are considered, namely, the Taylor-Green vortex at Re =?5000, the turbulent channel flow at Re_τ =?587, the circular cylinder in cross-flow at Re_D =?3900, the square cylinder in cross-flow at Re_D =?22400 and the channel with periodic constrictions at Re_h =?10595. Various discretization parameters such as the grid spacing, polynomial degree and numerical flux are assessed and very accurate results are reported in all cases. This consistency in the results demonstrates the versatility of the SEDM approach and its ability to gage the actual resolution and quality of the mesh and, accordingly, to introduce an amount of sub-grid dissipation which is adapted to the spatial discretization considered.     

Suspensions of titania nanoparticle networks in nematic liquid crystals: rheology and microstructure

We study the influence of confinement on the rheology and structure of nematic liquid crystals (NLCs). NLCs get confined in networks of titania (TiO₂, primary particle size = 21 nm) nanoparticles in suspensions of TiO₂ and NLC, N-(4-methoxybenzylidene)-4-butylaniline (MBBA). Suspensions with TiO₂ nanoparticle volume fraction (?) of 0.006–0.017, form viscoelastic solids with low elastic modulus (G′) of 10¹ Pa–10² Pa and short relaxation times. Increase in TiO₂ nanoparticle ? leads to a rise in G′ with TiO₂ nanoparticles forming a percolating network at a critical volume fraction (? _c) = 0.023, and G′ of ~10³ Pa. TiO₂/MBBA NLC suspensions at and above ? _c = 0.023 show G′ ~ ω ^x?1 scaling, where ω is the angular frequency and the minimum in loss modulus (G′′) with ω. The effective noise temperature, x decreases and approaches 1 with the increase in the TiO₂ nanoparticle ? from 0.023–0.035, is indicative of an increase in the glassy dynamics. Through the polarized light microscopy and differential scanning calorimetry experiments, we propose that the progressive addition of TiO₂ nanoparticles introduces a quenched random disorder (QRD) in the NLC medium which disturbs the nematic order. This results in metastable TiO₂/MBBA NLC suspensions in which NLC domains get confined in the network of flocs of TiO₂ nanoparticles. We also show that the salient rheological signatures of soft glassy rheology develop only in the presence of NLC MBBA and are absent in the isotropic phase of MBBA.     

Multipulse Phases in k-Mixtures of Bose–Einstein Condensates

For the system

$-\Delta U_i+ U_i=U_i^3-\beta U_i\sum_{j\neq i}U_j^2,\quad i=1,\dots,k,$

(with k ≧ 3), we prove the existence for β large of positive radial solutions on \({\mathbb R^N}\) . We show that as β →  + ∞, the profile of each component U _i separates, in many pulses, from the others. Moreover, we can prescribe the location of such pulses in terms of the oscillations of the changing-sign solutions of the scalar equation  ? ΔW  +  W  =  W³. Within an Hartree–Fock approximation, this provides a theoretical indication of phase separation into many nodal domains for the k-mixtures of Bose–Einstein condensates.

     

Effects of Fuel Lewis Number on Localised Forced Ignition of Globally Stoichiometric Stratified Mixtures: a Numerical Investigation

The influences of fuel Lewis number Le_F on localised forced ignition of globally stoichiometric stratified mixtures have been analysed using three-dimensional compressible Direct Numerical Simulations (DNS) for cases with Le_F ranging from 0.8 to 1.2. The globally stoichiometric stratified mixtures with different values of root-mean-square (rms) equivalence ratio fluctuation (i.e. ?^′= 0.2, 0.4 and 0.6) and the Taylor micro-scale l_? of equivalence ratio ? variation (i.e. l_?/l_f= 2.1, 5.5 and 8.3 with l_f being the Zel’dovich flame thickness of the stoichiometric laminar premixed flame) have been considered for different initial rms values of turbulent velocity u^′. A pseudo-spectral method is used to initialise the equivalence ratio variation following a presumed bi-modal distribution for prescribed values of ?^′ and l_?/l_f for global mean equivalence ratio 〈?〉=1.0. The localised ignition is accounted for by a source term in the energy transport equation that deposits energy for a stipulated time interval. It has been observed that the maximum values of temperature and the fuel reaction rate magnitude increase with decreasing Le_F during the period of external energy deposition. The initial values of Le_F, u^′/S_b(?=1), ?^′ and l_?/l_f have been found to have significant effects on the extent of burning of the stratified mixtures following localised ignition. For a given value of u^′/S_b(?=1), the extent of burning decreases with increasing Le_F. An increase in u^′ leads to a monotonic reduction in the burned gas mass for all values of Le_F in all stratified mixture cases but an opposite trend is observed for the Le_F=0.8 homogeneous mixture. It has been found that an increase in ?^′ has adverse effects on the burned gas mass, whereas the effects of l_?/l_f on the extent of burning are non-monotonic and dependent on ?^′ and Le_F. Detailed physical explanations have been provided for the observed Le_F, u^′/S_b(?=1), ?^′ and l_?/l_f dependences.     

Turbulent Drag Reduction by a Near Wall Surface Tension Active Interface

In this work we study the turbulence modulation in a viscosity-stratified two-phase flow using Direct Numerical Simulation (DNS) of turbulence and the Phase Field Method (PFM) to simulate the interfacial phenomena. Specifically we consider the case of two immiscible fluid layers driven in a closed rectangular channel by an imposed mean pressure gradient. The present problem, which may mimic the behaviour of an oil flowing under a thin layer of different oil, thickness ratio h₂/h₁ =?9, is described by three main flow parameters: the shear Reynolds number Re _τ (which quantifies the importance of inertia compared to viscous effects), the Weber number We (which quantifies surface tension effects) and the viscosity ratio λ = ν₁/ν₂ between the two fluids. For this first study, the density ratio of the two fluid layers is the same (ρ₂ = ρ₁), we keep Re _τ and We constant, but we consider three different values for the viscosity ratio: λ =?1, λ =?0.875 and λ =?0.75. Compared to a single phase flow at the same shear Reynolds number (Re _τ =?100), in the two phase flow case we observe a decrease of the wall-shear stress and a strong turbulence modulation in particular in the proximity of the interface. Interestingly, we observe that the modulation of turbulence by the liquid-liquid interface extends up to the top wall (i.e. the closest to the interface) and produces local shear stress inversions and flow recirculation regions. The observed results depend primarily on the interface deformability and on the viscosity ratio between the two fluids (λ).     

Growths of mechanical elasticity and electrical conductance of graphene nanoplatelet/poly(lactic acid) composites under strong electric field: correlation with time evolution of higher order structure of graphene nanoplatelets

For a composite of poly(lactic acid) containing graphene nanoplatelets (GNPs) at a low concentration (0.34 vol%), this study examined growths of mechanical and electrical properties under an alternating current (AC) electric field, focusing on field-induced GNP structures governing those properties. The composite was subjected to the AC field (60 Hz) of various intensities E for various times t _E at 190 °C. A fraction of randomly oriented GNPs was aligned by the field and then connected into columns, as suggested from optical microscopy. This structural evolution led to qualitatively similar growths of low-frequency storage modulus and static electrical conductivity. The key quantity for understanding this growth was a time t _E* for occurrence of short circuit that detected formation of GNP columns conductively bridging the electrodes. The growths of both modulus and conductivity for various E were summarized as functions of a reduced variable, t _E/t _E*, confirming the growths commonly reflected the evolution of the GNP columns. However, the modulus grew fast and leveled off by t _E/t _E* ~ 1, whereas the conductivity kept growing gradually even at t _E/t _E* > 1. This difference was discussed in relation to the matrix chains and leftover GNPs out the column.     

Effects of Fuel Lewis Number on Localised Forced Ignition of Turbulent Mixing Layers

The influences of fuel Lewis number Le _F on localised forced ignition of inhomogeneous mixtures are analysed using three-dimensional compressible Direct Numerical Simulations (DNS) of turbulent mixing layers for Le _F  = 0.8, 1.0 and 1.2 and a range of different root-mean-square turbulent velocity fluctuation u′ values. For all Le _F cases a tribrachial flame has been observed in case of successful ignition. However, the lean premixed branch tends to merge with the diffusion flame on the stoichiometric mixture fraction isosurface at later stages of the flame evolution. It has been observed that the maximum values of temperature and reaction rate increase with decreasing Le _F during the period of external energy addition. Moreover, Le _F is found to have a significant effect on the behaviours of mean temperature and fuel reaction rate magnitude conditional on mixture fraction values. It is also found that reaction rate and mixture fraction gradient magnitude \(\vert \nabla \xi \vert \) are negatively correlated at the most reactive region for all values of Le _F explored. The probability of finding high values of \(\vert \nabla \xi \vert \) increases with increasing Le _F . For a given value of u′, the extent of burning decreases with increasing Le _F . A moderate increase in u′ gives rise to an increase in the extent of burning for Le _F  = 0.8 and 1.0, which starts to decrease with further increases in u′. For Le _F  = 1.2, the extent of burning decreases monotonically with increasing u′. The extent of edge flame propagation on the stoichiometric mixture fraction ξ = ξ _st isosurface is characterised by the probability of finding burned gas on this isosurface, which decreases with increasing u′ and Le _F . It has been found that it is easier to obtain self-sustained combustion following localised forced ignition in case of inhomogeneous mixtures than that in the case of homogeneous mixtures with the same energy input, energy deposition duration when the ignition centre is placed at the stoichiometric mixture. The difficultly to sustain combustion unaided by external energy addition in homogeneous mixture is particularly prevalent in the case of Le _F  = 1.2.     

Brinkman Flow of a Viscous Fluid Through a Spherical Porous Medium Embedded in Another Porous Medium

An analytical investigation for a two-dimensional steady, viscous, and incompressible flow past a permeable sphere embedded in another porous medium is presented using the Brinkman model, assuming a uniform shear flow far away from the sphere. Semi-analytical solutions of the problem are derived and relevant quantities such as velocities and shearing stresses on the surface of the sphere are obtained. The streamlines inside and outside the sphere and the radial velocity are shown in several graphs for different values of the porous parameters \({\sigma _1 =(\mu /\tilde {\mu }) (a/\sqrt{K_1 })}\) and \({\sigma _2 =(\mu /\tilde {\mu }) (a/\sqrt{K_2 })}\) , where a is the radius of the sphere, μ is the dynamic viscosity of the fluid, \({\tilde {\mu }}\) is an effective or Brinkman viscosity, while K ₁ and K ₂ are the permeabilities of the two porous media. It is shown that the dimensionless shearing stress on the sphere is periodic in nature and its absolute value increases with an increase of both porous parameters σ ₁ and σ ₂.