Sedimentation of a Charged Soft Sphere within a Charged Spherical Cavity

The sedimentation of a soft particle composed of an uncharged hard sphere core and a charged porous surface layer inside a concentric charged spherical cavity full of a symmetric electrolyte solution is analyzed in a quasi-steady state. By using a regular perturbation method with small fixed charge densities of the soft sphere and cavity wall, a set of linearized electrokinetic equations relevant to the fluid velocity field, electrical potential profile, and ionic electrochemical potential energy distributions are solved. A closed-form formula for the sedimentation velocity of the soft sphere is obtained as a function of the ratios of core-to-particle radii, particle-to-cavity radii, particle radius-to-Debye screening length, and particle radius-to-porous layer permeation length. The existence of the surface charge on the cavity wall increases the settling velocity of the charged soft sphere, principally because of the electroosmotic enhancement of fluid recirculation within the cavity induced by the sedimentation potential gradient. When the porous layer space charge and cavity wall surface charge have the same sign, the particle velocity is generally enhanced by the presence of the cavity. When these fixed charges have opposite signs, the particle velocity will be enhanced/reduced by the presence of the cavity if the wall surface charge density is sufficiently large/small relative to the porous layer space charge density in magnitude. The effect of the wall surface charge on the sedimentation of the soft sphere increases with decreases in the ratios of core-to-particle radii, particle-to-cavity radii, and particle radius-to-porous layer permeation length but is not a monotonic function of the ratio of particle radius-to-Debye length.


Introduction
The sedimentation of charged particles in ionic fluids under gravity is a common phenomenon in various fields of colloidal science, as well as biomedical, mechanical, chemical, civil, and environmental engineering [1,2].This phenomenon is more complex than the migration of uncharged particles because the ambient fluid flow distorts the electrical double layer surrounding each charged particle and induces a sedimentation potential gradient [3,4].This induced electric field alters the velocity distribution of the ionic fluid via electrostatic interaction and diminishes the settling velocity of the charged particle by an electrophoretic effect.
Using a regular perturbation method with the particle zeta potential as a small perturbation parameter, Booth [5] first obtained analytical formulas for the sedimentation velocity and potential in suspensions of hard (impermeable to the ionic fluid) spheres with arbitrary double layer thickness as power expansions in their low zeta potential.Numerical results relaxing the assumption of low zeta potential in this analysis were calculated by Stigter [6].Ohshima et al. [7] obtained analytical formulas and numerical results of the sedimentation velocity and potential in suspensions of hard spheres for a wide range of zeta potential and double-layer thickness.On the other hand, the regular perturbation analyses have been extended to obtain the sedimentation velocity and potential in concentrated suspensions of interacting hard spheres [8][9][10][11][12], porous (permeable to the ionic fluid) spheres [13,14], and soft (hard core covered with porous surface layer) spheres [15,16].The equations that govern the ionic concentration, electric potential, and fluid flow fields around each charged particle migrating in an electrolyte solution are linearized, assuming that the system is only slightly distorted from the equilibrium.These linearized equations are solved with the fixed charge densities of the particles as the small perturbation parameters.A unit cell model that allows for the overlap of the electric double layers of adjacent particles is used to account for the particle-interaction effects.The sedimentation velocity and potential are not necessarily monotonic functions of the volume fraction of the particles.The particle concentration effects are significant, even in dilute suspensions.
In various applications of sedimentation, colloid particles are rarely unbounded but often settle close to solid boundaries [17,18].In the past, the boundary effect on sedimentation of uncharged particles were studied extensively [19][20][21][22][23], but not much information of this effect on charged particles was conveyed.Pujar and Zydney [24] calculated the settling velocity of a charged spherical particle in a concentric uncharged spherical cavity using perturbation expansions of small zeta potentials and Peclet number.This study has been extended to numerical calculations in the case of arbitrary zeta potentials [25].More recently, the settling of a charged hard [26] or porous [27] spherical particle within a concentric charged spherical cavity with arbitrary electric double layers has been analytically investigated for the case of small fixed charge densities of the particle and cavity wall.
In this paper, the sedimentation of a charged soft spherical particle inside a concentric charged spherical cavity with arbitrary double-layer thickness is analytically studied.The fluid velocity, electric potential, and ionic electrochemical potential energy distributions satisfying the linearized electrokinetic equations are determined as a power series of the small fixed charge densities of the soft sphere and cavity wall.A closed-form formula for the settling velocity of the soft sphere is obtained as a function of relevant parameters.

Electrokinetic Equations
As illustrated in Figure 1, we consider the quasi-steady sedimentation of a soft sphere of radius a, consisting of an uncharged hard sphere core of radius r 0 and a charged porous surface layer of thickness a − r 0 , inside a concentric charged spherical cavity of radius b full of the fluid solution of a symmetric electrolyte.The porous surface layer is permeable to the electrolyte solution and has fixed charges distributed at a uniform space density.The gravitational acceleration field ge z (e z is the unit vector in the z direction) is acting on the system, and the sedimentation velocity Ue z of the soft sphere will be determined.The origin of the spherical coordinates (r, θ, ϕ) is attached to the particle/cavity center (at z = 0), and the system is independent of ϕ (symmetric about the z axis).

Differential Equations
The system is assumed to deviate slightly from the equilibrium.Therefore, the pressure profile Figure 1.Geometric sketch for the sedimentation of a soft sphere inside a concentric spherical cavity.
The small perturbations δp, δψ, δn ± , and fluid velocity field u(r, θ) satisfy the continuity equation of the incompressible fluid (∇ • u = 0) following linearized electrokinetic equations [13]: resulting from the modified Stokes/Brinkman equation, Poisson's equation, and continuity equations of the ionic species, respectively.In Equations ( 2)-( 4), δµ ± (r, θ) are the perturbed quantities of the ionic electrochemical potential energies defined in terms of δψ and δn ± , η and ε are the viscosity and dielectric permittivity of the fluid, respectively, n ∞ and Z are the bulk concentration and valence, respectively, of the symmetric electrolyte, D ± is the diffusivity of the ionic species, 1/λ is the square root of the permeability of the fluid in the porous surface layer of the soft sphere, h(r) equals unity if r 0 ≤ r ≤ a and zero, otherwise k is Boltzmann's constant, T is the absolute temperature, and e is the charge of a proton.For porous particles made of plastic foam slab [28] and steel wool [29], experimental data for 1/λ can be 400 microns, while in the surface layers of grafted polymer microcapsules [30], rat lymphocytes [31], and human erythrocytes [32] in electrolyte solutions, 1/λ was found as low as 3 nm.

Boundary Conditions
The boundary conditions of the perturbations u, δψ, and δµ ± at the interface between the hard sphere core and the porous surface layer, as well as at the particle surface, are [15,33,34] r = r 0 : u = 0, e r • ∇δψ = 0, e r • ∇δµ ± = 0, (6) r = a : u, e r • τ, δψ, ∇δψ, δµ ± , and ∇δµ ± are continuous, (7) where τ is the viscous stress and e r is the unit vector in the r direction.Equation (7) represents the continuity requirements of the fluid velocity and stress, electrical potential and field, as well as ionic concentrations and fluxes at the interface.Various boundary conditions describing fluid flow at interfaces between porous media and surrounding fluids are discussed in detail in the literature [35,36] related to Darcy's law and Brinkman's equation, and Equation ( 7) is physically true and mathematically consistent with Equation (2).The boundary conditions of the small perturbations at the cavity wall are [26,27] Equations ( 6) and ( 8) take a reference frame moving with the particle and show that the Gauss condition holds at the surfaces of the hard sphere core and cavity wall, while no ions can penetrate into these surfaces.
Substituting Equations ( 9)-( 12) into the Gauss condition at r = b, we obtain the following relation between the zeta potential ζ and the surface charge density σ of the cavity wall confining the soft sphere: Namely, ψ (eq) expressed by Equation ( 9) after substituting Equation ( 13) is still valid for the case of constant zeta potential at the cavity wall.

Small Perturbations
To solve the small perturbations u, δp, δψ, and δµ ± in terms of the particle velocity U for small dimensionless charge densities Q and σ, these variables are expressed as perturbed expansions in powers of Q and σ up to the second orders, such as where the coefficients U ij with i and j equal to 0, 1, or 2. They are determined to be independent of Q and σ but are functions of the ratios of core-to-particle radii r 0 /a, particleto-cavity radii a/b, particle radius-to-Debye length κa, and particle radius-to-porous-layerpermeation length λa.In the expansions of δψ and δµ ± , there is no zeroth-order term of Q and σ because no ionic concentration gradient or electric field is applied.Substituting the expansions of u, δp, δψ, δµ ± , and U in the form of Equation ( 14) and Equation ( 9) for ψ (eq) into Equations ( 2)-( 8), we obtain the following solution for the components of u in spherical coordinates δp (to the orders Q 2 , Qσ, and σ 2 ), δψ, and δµ ± (to the orders Q and σ): Here, the functions ψ eq01 (r) and ψ eq10 (r) have been provided by Equations ( 10) and (11), and the functions F ijr (r), F pij (r), F ψij (r), and F ij± (r) are provided by Equations (A1)-(A4), (A9), and (A10) in Appendix A. Because F ψ01 (r), F ψ10 (r), F 01± (r), and F 10± (r) are influenced by the zeroth-order fluid flow field, the leading orders of the relaxation effect on the electrical double layers adjacent to the particle and cavity wall are contained in the solutions for δψ and δµ ± to the first orders Q and σ (which are sufficient for calculations of the settling velocity to the second orders Q 2 , Qσ, and σ 2 ).

Forces on the Particle
The net force exerted on the soft spherical particle undergoing sedimentation includes the gravity, electrical, and hydrodynamic forces.The gravity, which has nothing to do with the fixed and mobile electric charges, is where ρ p and ε p are the mass density and porosity, respectively, of the surface layer of the particle, ρ and ρ c are the mass densities of the fluid and hard core, respectively, and ge z is the gravitational acceleration.
The electrical force acting on the charged soft particle is [26,37] Substitution of Equations ( 9) and ( 16) into Equation ( 19) results in Only the second orders Q 2 , Qσ, and σ 2 contribute to F e , since both δψ (sedimentationinduced electric potential) and ψ (eq) are of the first orders Q and σ.
The hydrodynamic drag force exerted on the soft particle is provided by Substitution of Equation ( 15) into the previous equation leads to where the coefficients C 006 , C 026 , C 116 , and C 206 are provided in Equations (A1)-(A4).

Sedimentation Velocity
At the quasi-steady state, the net force acting on the charged soft sphere disappears.Using this constraint after adding Equations ( 18), (20), and (22), we obtain the sedimentation velocity of the particle inside the charged cavity in the expansion form of Equation ( 14) with the coefficients as where U 00 is the settling velocity of an uncharged soft spherical particle inside a concentric uncharged spherical cavity [40], which is positive but smaller than that in an unbounded fluid.As expected, the correction to the sedimentation velocity of the charged soft sphere begins with the second order of the fixed charge densities.
The sedimentation velocity of a charged soft spherical particle within a charged cavity provided by Equations ( 14) and ( 23)-( 25) can also be expressed as where U 00 (κa) 3 , (27b) Note that (κa)σ and (κa) 2 Q do not depend on κ or n ∞ .For a soft sphere undergoing sedimentation inside a cavity filled with a symmetric electrolyte solution, the dimensionless second-order coefficients H 1 , H 2 , and H 3 are functions of the ratios of core-to-particle radii r 0 /a, particle-to-cavity radii a/b, particle radius-to-Debye length κa, and particle radius-toporous-layer-permeation length λa.
The terms involving the coefficients H 1 and H 3 can be viewed as the corrections to the sedimentation of a charged soft spherical particle inside a concentric uncharged spherical cavity (σ = 0) and an uncharged soft particle (Q = 0) inside a concentric charged cavity, respectively.The surface charges of the cavity wall alter the settling velocity of the soft sphere by means of an electroosmotic recirculating flow developed from interactions of the electrical potential gradient caused by the sedimentation with the electrical double layer adjoining the cavity wall and a wall-induced electrical potential on the soft sphere.In the limiting case of r 0 /a = 0, Equation (26) reduces to the sedimentation velocity formula for a charged porous spherical particle within a concentric charged cavity [27].

Results and Discussion
The sedimentation velocity of a charged soft spherical particle within a concentric charged spherical cavity full of the solution of a symmetric electrolyte is expressed by Equations ( 26) and ( 27) as a power expansion of the fixed charge densities of its porous surface layer Q and the cavity wall σ up to the second orders Q 2 , Qσ, and σ 2 .

The Coefficients H 1 , H 2 , and H 3
For the settling of a soft particle within a cavity filled with an aqueous solution of potassium chloride (KCl, with εk 2 T 2 /ηD ± Z 2 e 2 = 0.26 at room temperature), the coefficients H 1 , H 2 , and H 3 , calculated from Equations ( 25) and ( 27), are plotted versus the ratios of particle-to-cavity radii a/b, particle radius-to-Debye length κa, particle radius-toporous-layer-permeation length λa, and core-to-particle radii r 0 /a in Figures 2-4 for a wide range.The coefficient H 1 (and U 02 ) is negative, so the presence of stationary space charges in the porous surface layer of the soft sphere reduces its sedimentation velocity.This retarding effect is caused by the electrophoresis in the opposite direction generated by the sedimentation-induced electrical potential field, as provided by Equation (16).For specified values of λa, κa, and r 0 /a, as shown in Figure 2b,d, −H 1 first increases with an increase in the particle-to-cavity radius ratio a/b, reaches a maximum, and then decreases with a further increase in a/b.It is understood that both −U 02 and U 00 , provided by Equations ( 23) and (25a), are monotonic-decreasing functions of a/b, due to the electrophoretic and viscous retardation effects, respectively, of the cavity wall on the moving soft sphere.Compared to U 00 , −U 02 decreases slow for a low value of a/b and fast for a high value of a/b, causing −H 1 to reach a maximum according to Equation (27a).As κa increases or λa decreases, the position of this maximum moves to larger a/b.U ) is negative, so the presence of stationary space charges in the porous surface layer of the soft sphere reduces its sedimentation velocity.This retarding effect is caused by the electrophoresis in the opposite direction generated by the sedimentation-induced electrical potential field, as provided by Equation (16).For specified values of a λ , a κ , and 0 / r a , as shown in Figure 2b,d,  Both coefficients H 2 and H 3 are positive, so the contribution of H 3 from σ 2 increases the settling velocity of an uncharged soft sphere in the cavity.The counterion concentration in the electric double layer adjoining the cavity wall near the front side of the settling particle increases because of the squeezing of the fluid, while the counterion concentration near the back side of the soft sphere decreases due to fluid compensation, generating a sedimentation-induced electrical potential gradient and electroosmosis along the cavity wall to enhance the recirculating flow and particle sedimentation.The contribution of H 2 of interactions between the stationary space charge of the porous surface layer of the soft sphere and the surface charge of the cavity wall increases the sedimentation velocity if these charges have the same sign (Qσ > 0), while counterions near the particle surface and cavity wall are reduced by mutual competition, thereby dropping the sedimentation-induced electric field and electrophoretic retardation and decreases the particle velocity if these charges are opposite in sign (Qσ < 0), and the enrichment of the counterions in the double layers enhances the induced potential gradient and electrophoretic retardation.For the provided values of λa, κa, and r 0 /a, as shown in Figures 3b,d and 4b,d, the coefficients H 2 and H 3 are decreasing functions of the particle-to-cavity radius ratio a/b (because the increase of a/b decreases the surface area of the cavity for a particular particle and, therefore, reduces the effect of cavity surface charge on the settling particle).
For the fixed values of λa, a/b, and r 0 /a, as shown in Figures 2a, 3a, and 4a, the coefficients −H 1 , H 2 , and H 3 have maxima at some finite values of the ratio of particle radius to Debye length κa and gradually fade as κa becomes smaller or larger.In general, as a/b increases, the location of these maxima shifts to larger κa but is insensitive to changes in λa.In the limiting case of very thick double layers ( κa → 0 ), the counterions in them are negligible and the particle sedimentation is not affected by the interaction between the space charge of the porous surface layer of the soft sphere and the surface charge of the cavity wall.In the case of very thin double layers ( κa → ∞ ), the charge density vanishes everywhere in fluid and the interaction between fixed charges disappears.
For specified values of a/b, r 0 /a, and κa, as shown in Figures 2a,c, 3a,c and 4a,c, the coefficients −H 1 , H 2 , and H 3 are decreasing functions of the ratio of particle radius-toporous-layer-permeation length λa as expected.Although H 3 is only a relatively weak function of λa (the charged cavity wall effect on the settling of the soft particle principally through the electroosmotic recirculating flow has not much to do with the relative permeability of its porous surface layer), H 1 and H 2 can strongly depend on λa.
For provided values of κa, λa, and a/b, as shown in Figures 2c,d and 3c,d, the coefficients −H 1 and H 2 decrease monotonically and substantially with an increase in the core-to-particle radius ratio r 0 /a (i.e., a decrease in the relative volume, and thus total space charge, of the porous surface layer of the soft sphere), as expected.On the other hand, as shown in Figure 4c,d, the coefficient H 3 , due to the surface charge of the cavity wall (it has little to do with the space charge of the porous surface layer), can be a relatively weak nonmonotonic function of r 0 /a.It is understood that, as r 0 /a increases, both U 00 and U 20 , provided by Equations ( 23) and (25c), decrease due to the effects of viscous retardation and mainly electroosmotic recirculation, respectively, of the cavity wall on the moving soft sphere.Compared to U 00 , U 20 decreases slower for low values of r 0 /a and faster for high values of r 0 /a, causing H 3 to reach a maximum according to Equation (27c).

The Normalized Sedimentation Velocity
Our results in Figures 2-4 show that the coefficient H 2 in Equation ( 26) is an order of magnitude larger than the other coefficients −H 1 and H 3 .Consequently, for a charged soft spherical particle migrating inside a spherical cavity with a surface charge of the same sign (Qσ > 0), the net effect of these coefficients is generally an increase of the sedimentation velocity of the particle (with exceptions occurring as σ is much smaller than Q in magnitude).For the case of the cavity with a surface charge of the opposite sign (Qσ < 0), if |σ| is sufficiently larger/smaller than Q , the net effect is to enhance/reduce the settling velocity.This trend is manifested in Figure 5, which is a plot of the normalized particle velocity U/U 00 of a soft spherical particle with a space charge density in its porous surface layer Q = 1 inside a spherical cavity full of an aqueous solution of potassium chloride calculated using Equation (26) versus the surface charge density σ of the cavity wall.The enhancement and reduction of the particle velocity caused by the cavity wall can be substantial if |σ| is relatively large, especially in cases of small core-to-particle radius ratio r 0 /a, ratio of particle radius-to-porous-layer-permeation length λa, and particle-tocavity radius ratio a/b, with moderate values of the ratio of the particle radius to Debye length κa.density σ of the cavity wall.The enhancement and reduction of the particle velocity caused by the cavity wall can be substantial if σ is relatively large, especially in cases of small core-to-particle radius ratio 0 / r a , ratio of particle radius-to-porous-layer-  ,d, 4b,d and 5a show that the influence of the charged cavity wall on the sedimentation of the soft sphere increases as b a / decreases, with other parameters remaining constant.This result reflects the fact that for a fixed size particle, the reduction of b a / increases the surface area of the cavity and, therefore, enhances the effect of its surface charge on the settling particle.Figure 5b,c,d illustrates that the effect of the cavity wall surface charge on the sedimentation of the soft sphere increases with decreases in 0 / r a and a λ but is not a monotonic function of a κ .

Conclusions
This paper presents an analysis of the sedimentation of a soft spherical particle inside a concentric spherical cavity filled with a symmetric electrolyte solution.By using a regular perturbation method with small fixed charge densities of the porous surface layer of the soft sphere Q and cavity wall σ , a set of linearized electrokinetic equations Figures 3b,d, 4b,d and 5a show that the influence of the charged cavity wall on the sedimentation of the soft sphere increases as a/b decreases, with other parameters remaining constant.This result reflects the fact that for a fixed size particle, the reduction of a/b increases the surface area of the cavity and, therefore, enhances the effect of its surface charge on the settling particle.Figure 5b,c,d illustrates that the effect of the cavity wall surface charge on the sedimentation of the soft sphere increases with decreases in r 0 /a and λa but is not a monotonic function of κa.

Conclusions
This paper presents an analysis of the sedimentation of a soft spherical particle inside a concentric spherical cavity filled with a symmetric electrolyte solution.By using a regular perturbation method with small fixed charge densities of the porous surface layer of the soft sphere Q and cavity wall σ, a set of linearized electrokinetic equations related to the fluid velocity field, electrical potential profile, and ionic electrochemical potential energy distributions are solved with the relaxation effect of the electric double layers.A closedform formula for the settling velocity of the soft sphere is provided by Equation (26) as a function of the ratios of particle-to-cavity radii a/b, particle radius-to-Debye length κa, particle radius-to-porous-layer-permeation length λa, and core-to-particle radii r 0 /a.Corrections due to the effect of the fixed charges Q and σ on the sedimentation velocity start from the second orders Q 2 , Qσ, and σ 2 .The existence of the surface charge on the cavity wall increases the settling velocity of the charged soft sphere, mainly because of the electroosmotic enhancement of recirculating flow within the cavity induced by the sedimentation potential gradient.When Q and σ have the same sign, the particle velocity is generally enhanced by the presence of the cavity.When these fixed charges have opposite signs and σ is sufficiently large/small relative to Q in magnitude, the particle velocity will be enhanced/reduced by the presence of the charged cavity.The effect of σ on the sedimentation of the soft sphere increases with decreasing a/b, λa, and r 0 /a but has a maximum at some finite values of κa and disappears as κa approaches zero and infinity.Our results might be used to optimize sedimentation processes, design materials with desired settling properties, or study efficient flotation of low-rank coal with poor hydrophobicity [41,42].

Figure 1 .
Figure 1.Geometric sketch for the sedimentation of a soft sphere inside a concentric spherical cavity.

Figure 2 .Figure 3 .coefficient 2 H
Figure 2. The coefficient H 1 (accounting for the effect of the particle charge) in Equation (26) for the sedimentation of a soft particle in a cavity full of aqueous KCl solution: (a) versus the ratio of particle radius to Debye length κa with r 0 /a = 0 and a/b = 0.5; (b) versus the particle-to-cavity radius ratio a/b with r 0 /a = 0 and λa = 1; (c) versus the ratio of particle radius to porous layer permeation length λa with κa = 1 and a/b = 0.5; (d) versus the core-to-particle radius ratio r 0 /a with κa = 1 and λa = 1.

Figure 3 .
Figure 3.The coefficient H 2 (accounting for the coupling effect of the fixed charges) in Equation (26) for the sedimentation of a soft particle in a cavity full of aqueous KCl solution: (a) versus the ratio of particle radius-to-Debye length κa with r 0 /a = 0 and a/b = 0.5; (b) versus the particle-to-cavity radius ratio a/b with r 0 /a = 0 and λa = 1; (c) versus the ratio of particle radius-to-porous-layerpermeation length λa with κa = 1 and a/b = 0.5; (d) versus the core-to-particle radius ratio r 0 /a with κa = 1 and λa = 1.

Figure 4 .coefficient 3 H 1 H
Figure 4.The coefficient 3 H (accounting for the effect of the cavity charge) in Equation (26) for the sedimentation of a soft particle in a cavity full of aqueous KCl solution: (a) versus the ratio of particle radius-to-Debye length a κ with 0 / 0 r a = and 5 .0 / = b a ; (b) versus the particle-to-cavity

1 H
− first increases with an increase in the particle-to-cavity radius ratio b a / , reaches a maximum, and then decreases with a further increase in b a / .It is understood that both 02 U − and 00 U , provided by Equations (23) and (25a), are monotonic-decreasing functions of b a / , due to the electrophoretic and viscous retardation effects, respectively, of the cavity wall on the moving soft sphere.Compared to 00 U ,

Figure 4 .
Figure 4.The coefficient H 3 (accounting for the effect of the cavity charge) in Equation (26) for the sedimentation of a soft particle in a cavity full of aqueous KCl solution: (a) versus the ratio of particle radius-to-Debye length κa with r 0 /a = 0 and a/b = 0.5; (b) versus the particle-to-cavity radius ratio a/b with r 0 /a = 0 and λa = 1; (c) versus the ratio of particle radius-to-porous layer permeation length λa with κa = 1 and a/b = 0.5; (d) versus the core-to-particle radius ratio r 0 /a with κa = 1 and λa = 1.

Figure 5 .
Figure 5.The normalized sedimentation velocity U/U 00 of a soft particle with the dimensionless space charge density Q = 1 in a cavity full of aqueous KCl solution versus the dimensionless surface charge density σ: (a) for several values of the particle-to-cavity radius ratio a/b with r 0 /a = 0.5, λa = 1, and κa = 1; (b) for several values of the ratio of particle radius to Debye length κa with r 0 /a = 0.5, a/b = 0.5, and λa = 1; (c) for several values of the ratio of particle radius to porous layer permeation length λa with r 0 /a = 0.5, a/b = 0.5, and κa = 1; (d) for several values of the core-to-particle radius ratio r 0 /a with a/b = 0.5, λa = 1, and κa = 1.