investigations of dynamical properties of fluid  fluid dispersions ( e.g. , emulsions @xcite and polymer blends @xcite ) require detailed understanding of the behavior of viscous drops in creeping flows .
such understanding is also crucial in development of new drop - based microfluidic systems @xcite .
therefore , drop dynamics at small reynolds numbers have been extensively studied experimentally @xcite , computationally @xcite and theoretically @xcite .
these investigations revealed complex nonlinear drop dynamics resulting from the coupling between the drop shape and fluid flow .
examples of nonlinear phenomena that occur under creeping - flow conditions include formation of self - similar neck regions during a drop breakup process @xcite , universal slow evolution of drops near the critical flow strength above which there are no stationary drop shapes @xcite , and existence of two branches of stable stationary shapes of highly viscous drops in two - dimensional stokes flows with nonzero vorticity @xcite .
as revealed by the analysis presented by @xcite , there exists a flow - parameter range where a high - viscosity drop can either adopt a nearly spherical shape stabilized primarily by the rotational flow component or an elongated shape stabilized primarily by capillary forces .
abrupt changes of the drop shape from one state to the other can be used in manipulation of emulsion micro - structure and for controlling the behavior of highly viscous drops in microfluidic devices . due to discontinuous changes of emulsion microstructure ,
the bistable drop behavior may also significantly affect emulsion rheology .
moreover , the mechanism of the bistability is of fundamental interest because of close analogies to the dynamics of vesicles @xcite and macromolecules @xcite in external flows .    due to the fundamental significance and because of potential applications ( such as those mentioned above ) it is important to explore the dynamics of highly viscous drops in linear flows with nonzero vorticity .
however , the investigations presented by @xcite were limited to stationary drop shapes and stationary external flows . in the present study we focus on drop behavior in time - dependent flows .
we elucidate the physical mechanism that give rise to the bistable drop behavior and examine the consequence of these mechanism for drop response to time variation of the fluid vorticity .
the system dynamics is investigated via direct boundary - integral simulations @xcite and by using a small - deformation approach @xcite .
in particular we show that the small - deformation equations with only several essential terms retained reproduce complex dynamical features of drop evolution that are associated with drop bistability .    to emphasize important aspects of the drop dynamics we consider three flow variation protocols . in the first protocol ,
the vorticity is slowly increased and then decreased .
we find that such quasistatic vorticity ramping gives access to both the elongated and compact , nearly spherical stationary drop shapes .
the drop exhibits a hysteretic behavior , with transitions between the compact shape ( rotationally stabilized ) and elongated shape ( stabilized by capillary forces ) occurring at different values of the vorticity when it is slowly ramped up or down .    in the second protocol
, the vorticity undergoes finite - frequency periodic oscillations .
as expected , at low frequencies the drop behavior is quasistatic , with a hysteresis loop analogous to the one observed for linear ramping . at high frequencies the oscillations average out , and the drop undergoes small oscillations around the stationary shape corresponding to the average flow .
however , at intermediate frequencies we find a much more complex behavior . in particular we show that there exists a frequency and amplitude domain where the drop response to the periodic forcing is chaotic .
since in the creeping flow regime fluid motion is governed by the linear stokes equations , the nonlinear chaotic behavior of the drop stems entirely from the coupling of the drop shape to the fluid velocity .
an analysis of drop motion in the chaotic domain indicates that the transition to chaos is associated with the existence of two stationary states observed in steady flow .    in our third flow - variation protocol
, the vorticity of the imposed flow undergoes random changes .
we observe , that the resulting statistical distribution of the drop length is bimodal in a certain regime of flow parameters , with two peaks around the drop length corresponding to the short and long stationary solutions .
hence , we find that also in this problem the existence of two stationary states underpins drop behavior in a time - dependent flow .    in  [ viscous drop in linear flows ] the system considered in our paper is defined . the quasistatic hysteretic drop behavior is analyzed in
[ hysteretic behavior ] , our results for chaotic drop dynamics are presented in
[ sec : sinusoidal_chaos ] , and drop motion in linear flows with randomly varying vorticity is discussed in  [ drop statistics in a linear flow with stochastic vorticity ] .
our findings are summarized in
[ sec : conclusion ] .
we consider a viscous drop suspended in an incompressible newtonian fluid of a constant viscosity @xmath0 .
the viscosity of the drop fluid is @xmath1 , and the interfacial tension between the two phases is @xmath2 .
the drop is surfactant free , and no marangoni stresses are present .
there are also no buoyancy forces .
we focus here on nonlinear effects that stem entirely from the coupling of the fluid flow to the drop shape ( but not from the fluid inertia ) .
therefore , the creeping - flow conditions are assumed .    in the creeping - flow regime
the fluid motion in the regions inside ( @xmath3 ) and outside ( @xmath4 ) the drop is governed by the stokes equations @xmath5 the fluid velocity @xmath6 is continuous at the drop interface @xmath7 . due to
the absence of the marangoni stresses the tangential viscous traction is also continuous .
the jump in the normal viscous traction across @xmath7 is equal to the capillary pressure @xmath8 = 2\kappa\sigma,\ ] ] where @xmath9 is the viscous stress tensor , @xmath10 is the outward normal unit vector , and @xmath11 is the local curvature .
( 300,8)(0,0 ) ( 75,0)@xmath12 ( 176,0)@xmath13 ( 290,0)@xmath14    the drop is subject to two - dimensional linear incident flow @xmath15 where @xmath16 is the strain rate , @xmath17 is the dimensionless vorticity parameter , @xmath18 is the position , and @xmath19 and @xmath20 are the symmetric and antisymmetric parts of the normalized velocity - gradient tensor . in an appropriately adopted coordinate system
we have @xmath21 without a loss of generality .
accordingly , @xmath22 corresponds to a purely straining flow with the extensional axis @xmath23 and the compressional axis @xmath24 , and @xmath25 corresponds to shear flow in the @xmath26 direction with the velocity gradient in the @xmath27 direction .
the tensor @xmath20 in equation describes the rigid - body rotation in the anti - clockwise direction with the angular velocity @xmath28 the decomposition of the flow field into the straining and rotational components is sketched in figure [ flow decomposition ] .
the dynamics of our system is characterized by three dimensionless parameters : there is the viscosity ratio @xmath29 that characterizes dissipative forces in the drop- and continuous - phase fluids .
the capillary number @xmath30 ( where @xmath31 is the radius of an undeformed drop ) characterizes the ratio between the deforming viscous forces produced by the imposed flow and the capillary forces that resist drop deformation and drive the drop towards the equilibrium spherical shape .
finally the vorticity parameter @xmath17 characterizes the magnitude of the rotational component of the external flow relative to the extensional component .    in this paper
we focus on the parameter regime where the drop deformation may be significant , which requires that @xmath32 .
( however , the flow is not strong enough to cause drop breakup . )
we also assume that the drop - phase fluid is much more viscous than the continuous - phase fluid , @xmath33 drop deformation process is hindered at large drop - phase viscosities , while drop rotation is only weakly affected by the viscous stresses inside the drops .
therefore the relative effect of the drop rotation is amplified in the regime , and the rotational component of the external flow produces nontrivial qualitative effects .     and capillary number @xmath34 .
normalized drop length @xmath35 and drop angle @xmath36 versus vorticity parameter @xmath17 .
arrows indicate the direction of increasing time , and inset shows dependence of @xmath17 on time .
( results from boundary - integral simulations . ) , title="fig:",width=240 ]   and capillary number @xmath34
. normalized drop length @xmath35 and drop angle @xmath36 versus vorticity parameter @xmath17 .
arrows indicate the direction of increasing time , and inset shows dependence of @xmath17 on time .
( results from boundary - integral simulations . ) , title="fig:",width=240 ]
to illustrate the effect of the rotational component of the flow on the dynamics of a highly viscous drop we consider a system where the parameter @xmath17 is first slowly increased and then slowly decreased .
we adopt here a linear ramping protocol where the vorticity is slowly ramped up from @xmath22 to @xmath37 and then ramped down back to zero .
before the ramping occurs the flow is maintained at @xmath22 ( for 5% of the total ramping time ) to allow the drop to relax to the stationary shape in purely straining flow .
this vorticity variation protocol is represented in the inset of figure [ fig:1 ] .
the evolution of the drop shape in this time - dependent flow is depicted in figure [ fig:1 ] .
figure [ fig:1 ] represents the drop length @xmath35 , and figure [ fig:1 ] shows the drop angle @xmath38 ( measured anticlockwise from the axis @xmath26 ) ; both quantities are plotted versus the vorticity parameter @xmath17 . in our example , the drop viscosity is @xmath39 .
the capillary number @xmath40 is below the critical value for drop breakup ( @xmath41 in 2d straining flow ) but it is sufficiently large to allow for a significant flow - induced drop deformation .
the total ramping time is @xmath42 , so that the drop response to the flow variation is nearly quasistatic .
the calculations were performed using the adaptive boundary - integral procedure developed by @xcite .
the results shown in figure [ fig:1 ] indicate that the drop response to the vorticity variation is hysteretic . at @xmath22 drop
is elongated , and it is aligned with the extensional axis of the straining component of the flow ( @xmath43 ) .
with increasing @xmath17 , the drop orientation slowly changes towards the symmetry axis @xmath26 ( @xmath44 ) , and the drop length slowly decreases . at a critical value of the vorticity parameter , @xmath45 ,
a discontinuous change occurs : the drop length and the angle suddenly decrease .
afterwards the drop is almost spherical and nearly aligned with the axis @xmath26 .
when the direction of the vorticity change is reversed , the drop initially retraces its trajectory .
however , the drop does not jump back to the elongated shape until the vorticity reaches the lower critical value @xmath46 .    the bistable drop behavior and the associated hysteretic shape evolution stem from the existence of two mechanisms that can stabilize a viscous drop in linear flows with rotation .
namely , the drop can be stabilized by the capillary stresses ( which drive the drop towards the equilibrium spherical shape ) or by the vorticity flow component ( which rotates the drop out of the extensional axis of the straining component of the flow ) .    in a purely straining flow
the drop assumes the interfacial - tension - stabilized elongated shape @xcite .
the form of the drop results from the balance between drop deformation by the flow and drop relaxation due to the capillary forces .
the deformation and relaxation occur on the respective time scales @xmath47 @xmath48 both of which are proportional the viscosity ratio @xmath49 .
the drop deformation @xmath50 is determined by the time scale ratio @xmath51 and therefore it is independent of the viscosity ratio in the limit @xmath52 . in purely staining flows , the drop is oriented along the extensional axis @xmath23 .        for small values of @xmath17 ,
the vorticity flow component produces an @xmath53 perturbation of the drop orientation .
the corresponding decrease of the drop length is @xmath54 .
however , a further drop rotation is arrested because the straining component of the flow produces hydrodynamic stresses that pull the elongated drop back towards the straining axis ( as illustrated in figure [ mechanism ] ) .
since the @xmath55 internal circulation inside an elongated high - viscosity drop is weak , the drop in its stationary state behaves analogously to a rigid body whose equilibrium orientation results from the balance of the torques produced by the straining and rotational components of the external flow .
the transition to the compact drop shape occurs when the vorticity flow component becomes too strong to be balanced . under such conditions ,
a rigid body would undergo a transition to a periodic motion with continuous rotation in the clockwise direction .
similarly , a drop also starts to continuously rotate when @xmath17 achieves the upper critical value @xmath56 .
however , during the rotation the drop length decreases because the drop becomes misaligned with the extensional axis of the flow . as a result
, the drop relaxes to a nearly spherical shape .    in this new , compact shape the fluid inside
the drop circulates with the angular velocity @xmath57 that is nearly equal to the angular velocity of the external flow . within each period of rotation
the drop undergoes a small deformation produced by the straining component of the external flow ( as schematically illustrated in figure [ mechanism ] ) .
however , the deformation does not accumulate because it is constantly convected away by the rotational component of the flow .    since the rotation occurs on the time scale @xmath58 and the drop deforms on the much longer timescale , we find that the drop deformation @xmath59 is small for @xmath49 , consistent with the results shown in figure [ fig:1 ] .
relation indicates that the deformation of the rotationally stabilized drop increases with the decreasing parameter @xmath17 .
when @xmath17 falls below the lower critical value @xmath60 , the hydrodynamic torque associated with the straining component of the flow acting on a slightly elongated drop becomes strong enough to reorient the drop along the straining axis and arrest further drop rotation .
deformation thus starts to accumulate , the drop is stretched , and a transition to the interfacial - tension stabilized elongated state takes place .    as shown in figure [ fig:1 ] , a drop in the compact , rotationally stabilized stationary state is nearly aligned with the symmetry axis of the applied flow @xmath26 .
this behavior stems from the flow - reflection symmetry of stokes equations and the fact that the drop is stabilized by rotation rather than the capillary forces . in the absence of the capillary forces ( or in the limit of infinitely strong flow )
the symmetry of stokes equations implies that the stationary drop shape is invariant with respect to flow reflection .
hence , the shape is also invariant with respect to the corresponding transformation @xmath61 of the spatial coordinates ( and this symmetry corresponds to drop alignment in the @xmath26 direction ) .
a perturbation due to the presence of the capillary stresses produces only a small asymmetry because the effect of capillary forces is insignificant for a nearly spherical drop .
with vorticity parameter @xmath17 for @xmath39 , 100 , 50 , and different values of capillary number ( as labeled ) .
( results from boundary - integral simulations . ) , title="fig:",width=307 ] +     with vorticity parameter @xmath17 for @xmath39 , 100 , 50 , and different values of capillary number ( as labeled ) .
( results from boundary - integral simulations . ) , title="fig:",width=240 ]   with vorticity parameter @xmath17 for @xmath39 , 100 , 50 , and different values of capillary number ( as labeled ) .
( results from boundary - integral simulations . ) , title="fig:",width=240 ]      the quasistatic response of the drop length to the variation of the vorticity is illustrated in figure [ drop response for different ca and lambda ] for different values of the capillary number @xmath62 and viscosity ratio @xmath29 .
the results indicate that the size of the hysteresis loop is the largest for large values of @xmath29 and @xmath62 .
when the capillary number is decreased , the upper critical vorticity parameter @xmath56 ( corresponding to the transition from the elongated to the compact drop ) decreases but the lower critical parameter @xmath63 remains nearly unaffected .
in contrast , the viscosity ratio @xmath29 affects primarily the position of the lower critical parameter @xmath60 ( corresponding to the transition from the compact to the elongated drop ) .
this behavior is consistent with the scaling relations and and the mechanism of drop bistability explained in  [ capillary and rotational stabilizing mechanisms ] .
according to these relations the deformation of an elongated drop scales with @xmath62 and is independent of @xmath29 , while the deformation of a compact drop is independent of @xmath62 and scales with @xmath64 .
we recall that the critical states corresponding to the transitions between the elongated and compact drop shapes correspond to the points where the maximal torque @xmath65 exerted by the extensional component of the flow on the drop just before the transition marginally balances the torque @xmath66 exerted by the vorticity flow component .
since @xmath67 whereas @xmath66 is approximately independent of @xmath68 , we conclude that @xmath60 varies with @xmath64 and @xmath56 with the @xmath62 .
this conclusion is consistent with our simulation results .
the plots shown in figure [ drop response for different ca and lambda ] indicate that for a given viscosity ratio , there exists a critical value of the capillary number @xmath69 below which the drop response to the changes of the flow vorticity does not exhibit a hysteretic loop .
the bifurcation point occurs at the critical value of the vorticity parameter @xmath70 that corresponds to the position of the infinitesimal hysteresis loop for @xmath62 slightly above the critical value @xmath69 .
it can be shown from the results derived by @xcite that @xmath71 for @xmath49 .
the plots shown in figure [ drop response for different ca and lambda ] are consistent with the above expressions .
we note that for @xmath72 ( see figure [ drop response for different ca and lambda ] ) the hysteretic drop behavior is not observed .
the reason is that in this viscosity range @xmath69 exceeds the critical value of the capillary parameter for drop breakup .
hence , the elongated stationary drop shape does not exist for @xmath73 .
crucial features of the evolution of a highly viscous drop in two - dimensional flows with nonzero vorticity are captured by small - deformation equations . in our approach @xcite , the position @xmath74 of the drop interface
is expanded in spherical harmonics    @xmath75,\ ] ]    where @xmath76 and @xmath77 denote the spherical - harmonic order , @xmath78 and @xmath79 are the expansion coefficients , and the parameter @xmath80 is given by the drop - volume constraint . since for @xmath81 all spherical harmonics are real , we set @xmath82 in expansion .
moreover , flow - induced drop deformation preserves the symmetry of the incident flow .
therefore only even values of @xmath83 and @xmath84 need to be included in the analysis .
evolution equation for the expansion coefficients @xmath78 and @xmath79 are obtained by inserting the series into the boundary - value problem  , performing a boundary - perturbation analysis , and reexpanding resulting products of spherical harmonics using appropriate clebsch
gordan coupling coefficients .
the detailed analysis and explicit expressions for the evolution equations at different truncation levels are presented elsewhere @xcite .    for simplicity ,
our small - deformation calculations are performed with the expansion truncated at the lowest spherical - harmonic order @xmath85 , which leaves us with three independent drop - shape components : @xmath86 , and @xmath87 .
noting that @xmath88 we find that the shape parameters @xmath89 and @xmath90 correspond to the drop deformation along the symmetry axis @xmath26 and the straining axis @xmath23 , respectively .
the parameter @xmath87 described an axisymmetric deformation along the axis @xmath91 .
the evolution equations for the shape parameters @xmath86 , and @xmath87 , truncated at the second - order in the drop deformation , can be represented in the following form    [ small - deformation equations ] @xmath92,\\ \label{equation f'22 } \dot{f'_{22 } } & = & -2\omega f''_{22}+\lambda^{-1}\left[d_{21}f'_{22}f''_{22 }      -{{\rm ca}}^{-1}(d_1 + 2d_2f'_{20})f'_{22}\right],\\ \label{equation f''22 } \dot{f''_{22 } } & = &   2\omega f'_{22}+\lambda^{-1 }      \left[(d_{31}+d_{32}f'_{20 }        + d_{33}f^{\prime\hspace{0.5pt}2}_{20 }        + d_{34}f^{\prime\hspace{0.5pt}2}_{22 }        + d_{35}f^{\prime\hspace{0.5pt}2}_{22})-                    { { \rm ca}}^{-1}(d_1 + 2d_2f'_{20})f''_{22}\right],\nonumber\\\end{aligned}\ ] ]    where the dot denotes the time derivative ( normalized by @xmath93 ) . the terms involving the coefficients @xmath94 correspond to drop deformation by the external flow , and the terms involving @xmath95 describe the capillary relaxation .
all these terms are @xmath55 in the large - viscosity - ratio regime .
explicit expressions for the coefficients @xmath94 and @xmath95 are given in @xcite ; here we only note that these coefficients are functions of the viscosity ratio @xmath29 and have finite limits for @xmath52 .
the two remaining terms on the right - hand side of equations and ( the terms proportional to @xmath96 ) are viscosity independent .
these terms represent the rigid - body rotation of the drop , with the angular velocity @xmath97 where @xmath98 .
consistent with our qualitative physical picture described in  [ capillary and rotational stabilizing mechanisms ] ( and illustrated in figure [ mechanism ] ) , the rotational velocity involves two terms .
the first term corresponds to the rotation of the drop by the vorticity component of the flow . the second term , which is proportional to the shape parameter @xmath89 that described deformation in the @xmath26 direction , corresponds to the rotation of a deformed drop by the straining component of the external flow towards the straining axis @xmath23 .
it has been shown by @xcite that for @xmath49 the drop behavior near the bifurcation point can be described by simplified asymptotic equations    [ simplified evolution equation ] @xmath99    where @xmath100 and @xmath101 are the high - viscosity limits of @xmath102 and @xmath103 .
the asymptotic result is obtained form on assumption that near the bifurcation point there is a balance between drop deformation and rotation ( which corresponds to @xmath104 ) and the balance between capillary relaxation and rotation ( which yields @xmath105 and @xmath106 ) . moreover , the two contributions to the angular velocity are of the same order but do not cancel ( i.e. , @xmath107 ) .
equations are rescaled accordingly , and only the leading order terms are retained .
equations have all necessary ingredients that are needed to describe the hysteretic drop behavior .
we have terms representing drop rotation by the straining and vorticity components of the external flow ( i.e. , the terms proportional to @xmath96 ) , drop relaxation due to the presence of the capillary forces ( the terms proportional to @xmath108 ) , and stretching of the drop along the straining axis of the straining component of the external flow ( the last term in equation ) . neglecting any of these terms
would qualitatively alter the solution structure , which no longer would manifest the key features of the drop evolution in the parameter range considered herein .      in the regime @xmath109 and
@xmath110 the stationary solutions of equations can be obtained by a singular - perturbation analysis .
to the leading order in the small parameters , we find that the elongated drop is described by the relations    [ equations for elongated drop ] @xmath111 @xmath112    the first of the above relations correspond to the fact that an elongated drop does not rotate ( there is only a weak fluid circulation inside it , as predicted by our qualitative analysis ) .
the second relation describes the balance between drop deformation by the external flow and relaxation due to the capillary forces . recalling the definition of the angular velocity @xmath96 we find that relations yield    [ solution for elongated drop ] @xmath113
@xmath114    by inserting the above relations back into one can verify that they constitute a consistent leading - order asymptotic stationary solution .    according to equation drop elongation along the straining axis
@xmath23 ( i.e. , @xmath43 ) scales with the capillary number , which is consistent with the scaling result .
the drop angle @xmath115 only slightly deviates from @xmath43 because @xmath116 ( assuming that @xmath117 and @xmath110 ) .
the leading - order stationary solution corresponding to the compact drop is obtained from the following relations    [ equations for compact drop ] @xmath118 @xmath119    which are obtained by dropping from the evolution equations the @xmath55 capillary - relaxation terms . taking into account the definition of @xmath96 we thus obtain    [ solution for compact drop ] @xmath120 and @xmath121    ( the solution with
the plus sign in front of the square root is unstable ) .
since the shape parameter @xmath90 vanishes according to equation , the drop is oriented in the @xmath26 direction .
for @xmath122 we find @xmath123 which is consistent with our scaling estimate .
for @xmath124 , where @xmath125 the solution does not exist ; the drop thus undergoes a transition to the long solution . the critical vorticity parameter decreases with increasing @xmath29 , in agreement with our numerical results presented in figure [ drop response for different ca and lambda ] .
the solution  of the simplified small - deformation equations is perturbative .
we note , however , the exact stationary solution can also be found @xcite .
our analytical solution quantitatively agrees with the results of numerical simulations , provided that the drop deformation is not too large .     with vorticity parameter @xmath17 for @xmath39 and different values of capillary number ( as labeled ) .
solution of small - deformation equations ( [ small - deformation equations ] ) ; comparison of small - deformation results ( dashed lines ) with boundary integral simulations ( solid lines ) .
, title="fig:",width=240,height=259 ]   with vorticity parameter @xmath17 for @xmath39 and different values of capillary number ( as labeled ) .
solution of small - deformation equations ( [ small - deformation equations ] ) ; comparison of small - deformation results ( dashed lines ) with boundary integral simulations ( solid lines ) . ,
title="fig:",width=240,height=259 ]      predictions of the small - deformation equations for drop behavior in a two - dimensional linear flow with slowly varying vorticity are depicted in figure [ small deformation results for lambda=200 ] for a system with the viscosity ratio @xmath39 .
figure [ small deformation results for lambda=200 ] shows the dependence of the drop length @xmath126 on the vorticity parameter @xmath17 for the same set of capillary numbers as those represented in figure [ drop response for different ca and lambda ] .
figure [ small deformation results for lambda=200 ] compares the small deformation results directly with the results of the boundary - integral simulations .
the small - deformation calculations were performed using the second - order equations because for an elongated drop they are more accurate than the simplified equations .
the results shown in figures [ drop response for different ca and lambda ] and [ small deformation results for lambda=200 ] indicate that for small and moderate capillary numbers the small - deformation theory yields accurate quantitative predictions . at high values of @xmath62 drop behavior
is also captured quantitatively , except for the upper portion of the hysteresis loop ( i.e.  when the drop is in the elongated state ) .
for all values of the capillary number the lower and upper critical vorticity parameters @xmath63 and @xmath127 are obtained within the numerical error of the boundary - integral simulations .
our additional calculations ( not shown ) indicate that a similar accuracy is obtained for @xmath72 and @xmath128 .
@xmath129 { figs / fig6a.eps } } & { \includegraphics[width=2.5in , height=4.8in , trim=20 0 0 0 ] { figs / fig6b.eps } } \end{array } $ ]     ( solid lines ) , drop deformation @xmath68 ( dashed ) , and drop angle @xmath36 ( dash - dotted ) are shown versus time @xmath130 normalized by the oscillation period @xmath131 .
viscosity ratio @xmath132 , capillary number @xmath40 , period @xmath133 , mean vorticity @xmath134 , and vorticity oscillation amplitude @xmath135 and 0.8 .
( results from small - deformation theory . ) , title="fig:",width=432 ]   ( solid lines ) , drop deformation @xmath68 ( dashed ) , and drop angle @xmath36 ( dash - dotted ) are shown versus time @xmath130 normalized by the oscillation period @xmath131 .
viscosity ratio @xmath132 , capillary number @xmath40 , period @xmath133 , mean vorticity @xmath134 , and vorticity oscillation amplitude @xmath135 and 0.8 .
( results from small - deformation theory . ) , title="fig:",width=432 ]
dynamical systems with multiple equilibrium states often exhibit novel dynamics when driven by simple forcing @xcite . thus , despite the laminar nature of the stokes flow , we expect to find interesting nonlinear dynamics of a viscous drop in a time - varying linear flow with rotation . to explore this dynamics
we will now investigate the drop response to harmonic variation of the vorticity @xmath136 where @xmath137 is the average vorticity value , @xmath138 is the oscillation amplitude , and @xmath131 is the oscillation period .
we have performed a series of small deformation calculations (  [ subsec : chaos_sdt ] ) and boundary - integral simulations (  [ subsec : chaos_bis ] ) for different values of the flow parameters @xmath137 , @xmath138 , and @xmath131 .
if the oscillation period @xmath131 is much shorter than the drop deformation and oscillation times and , we find that the drop undergoes small oscillations about a stationary shape corresponding to the mean value of @xmath17 ( which is an expected behavior ) . in the opposite regime @xmath139 , the quasistatic drop behavior described in
[ hysteretic behavior ] is recovered . in what follows we focus on the most interesting parameter domain @xmath140 and @xmath141 , in which an interaction of different timescales as well as an interplay between the short and elongated drop shapes is anticipated .       and
period @xmath133 . viscosity ratio and capillary number are the same as in figures [ fig : chaos01 ] and [ example of period doubling ] .
( results from small - deformation theory . ) , title="fig:",width=336 ]    figure [ fig : chaos01 ] illustrates the dependence of the drop evolution in linear flow with the oscillatory vorticity on the oscillation period @xmath131 .
the viscosity ratio in this example is @xmath132 , and the capillary number is @xmath40 .
the mean value of the vorticity @xmath134 is close to the lower critical value @xmath142 , and the oscillation amplitude is @xmath143 .
figure [ fig : chaos01 ] represent our results for the shortest oscillation period of the flow vorticity @xmath144 .
the drop oscillates about the compact stationary shape in this case .
both the drop lengths and the drop angle vary periodically , with the period @xmath145 equal to the period @xmath131 of the external forcing .
the drop length decreases when the drop is in the compressional quadrant @xmath146 and increases for @xmath147 . upon an increase of the period of the external forcing the amplitude of the angular drop oscillations increases . when the oscillation amplitude reaches @xmath148 the drop starts to tumble , as illustrated in figure [ fig : chaos01 ] .
@xmath129 { figs / fig9a.eps } } & { \includegraphics[width=2.5in , height=4.8 in ] { figs / fig9b.eps } } \end{array } $ ]    a further increase of the period of the external forcing results in a qualitative change of drop response .
we find that the drop still undergoes a tumbling motion ; however , the evolution is not periodic but it becomes _ chaotic _ , as shown in figure [ fig : chaos01 ] .
the chaotic motion continues up to @xmath133 , and then the drop reverts to periodic motion . for @xmath149 ( figure [ fig : chaos01 ] )
the drop oscillates about the elongated stationary shape , and in the regime @xmath150 ( figure [ fig : chaos01 ] ) the system approaches the quasistatic behavior discussed in  [ hysteretic behavior ] .
occurs right after the drop jumps from the long to the compact shape when @xmath17 increases above the upper critical value @xmath56 .
namely , before the drop settles down to the compact stationary shape it undergoes a tumbling motion with the amplitude decaying on the timescale . ] the transition to the chaotic drop motion occurs through a cascade of period doubling events , as illustrated in figures [ example of period doubling ] and [ fig : period_doubling ] .
figure [ example of period doubling ] depicts the drop evolution at three values of the amplitude of vorticity oscillations @xmath138 ( the remaining system parameters are the same as those that yield the chaotic motion depicted in figure [ fig : chaos01 ] ) .
the results shown in figure [ example of period doubling ] indicate that for a sufficiently small oscillation amplitude ( @xmath151 in our example ) the drop evolves with the period @xmath152 ( i.e. , the period equal to that of the external forcing ) . at a larger amplitude @xmath153 the drop oscillation period
is @xmath154 and for @xmath155 we find @xmath156 .
the period - doubling scenario of the transition to chaos in our system is further supported by the bifurcation diagram shown in figure [ fig : period_doubling ] , where the drop length @xmath35 at times @xmath157 ( @xmath158 ) is plotted versus the flow - oscillation amplitude @xmath138 .
an analysis of the results shown in figure [ example of period doubling ] indicates that the period doubling occurs as a result of a resonance between the drop tumbling motion and the vorticity oscillations .
namely , if the drop is relatively long and approximately aligned with the straining axis of the external flow when the vorticity parameter @xmath159 reaches a minimum , the drop rotation may be significantly slowed down or even arrested ( as in the long - drop stationary state discussed in
[ hysteretic behavior ] ) . such a temporary arrest of drop rotation corresponds to the shoulders in the plots of the angular evolution depicted in figure [ example of period doubling ] . on the other hand ,
if the drop angle exceeds @xmath160 when the vorticity goes through a minimum , the arrest of drop rotation does not occur .
as seen in figure [ example of period doubling ] this interplay of drop tumbling with oscillations of the external forcing produces the period - doubling bifurcation that leads to alternating accelerated and retarded drop - rotation cycles .    in our numerical calculations depicted in figures [ fig :
chaos01]-[fig : period_doubling ] we have used the full set of the small - deformation equations but we find that the simplified asymptotic equations yield similar results . in particular these equations correctly reproduce the cascade of the period doubling bifurcations and the chaotic - evolution domain . however , a further simplification of the evolution equations is not possible : if any of the terms in equations is removed the solutions qualitatively change and the chaotic domain disappears . we also find that chaotic drop dynamics occurs only for highly viscous drops with @xmath161 .    , @xmath162 , @xmath132 , and @xmath40 .
( results from boundary - integral simulations . ) , title="fig:",width=336 ]      figure [ fig : chaos03 ] shows examples of our simulation results for the periodic and chaotic drop evolution .
the flow parameters are similar to those used in our small - deformation calculations described in  [ subsec : chaos_sdt ] . as with the small - deformation calculations , for short periods of the external forcing @xmath131 the drop oscillates around the compact stationary shape , for moderate periods the system undergoes a transition to chaotic evolution , and for long periods
the drop motion approaches the quasistatic behavior .
consistent with the small - deformation results , the chaotic drop dynamics revealed by the boundary - integral simulations is due to a cascade of period doubling bifurcations .
a bifurcation diagram illustrating this behavior is presented in figure [ fig : period_doubling_dns ] .
the domain of chaotic motion found in the direct boundary - integral simulations somewhat differs from the corresponding domain obtained from the small - deformation theory . also , the magnitude of the chaotic fluctuations in the drop length is larger in the boundary - integral runs .
we have tested the convergence of the boundary - integral simulations , and we believe that the differences between the drop behavior obtained by the two different methods stem from the approximations involved in the small - deformation theory .
we note , however , that the evolution in the chaotic and period - doubling regimes is very sensitive to the to small changes of the system parameters ( hence , also sensitive to the approximations involved in our calculations ) .     and
@xmath40 , in linear flow with stochastic vorticity .
mean value of the vorticity parameter @xmath163 , variance @xmath164 , and correlation time @xmath165 .
inset shows vorticity probability distribution .
( results from boundary - integral simulations.),title="fig:",width=336,height=336 ]     and @xmath40 , in linear flow with stochastic vorticity . mean and variance of the vorticity parameter are the smae as in figure [ fig : cgaussdns ] ; the correlation times are @xmath166 for lines marked 15 , respectively .
( results from small - deformation theory . ) , title="fig:",width=249,height=211 ]   and @xmath40 , in linear flow with stochastic vorticity . mean and variance of the vorticity parameter are the smae as in figure [ fig : cgaussdns ] ; the correlation times are @xmath166 for lines marked 15 , respectively .
( results from small - deformation theory . ) , title="fig:",width=249,height=211 ]
in some systems ( e.g. emulsion flows through a packed bed of fibers @xcite or turbulent emulsion flows with drops that are much smaller then the kolmogorov scale @xcite ) a viscous drop undergoes deformation in a random external creeping flow . to gain some understanding of the role that drop bistability may play in such systems we consider the drop behavior in a flow with stochastic vorticity .
we assume that the time variation of the vorticity parameter @xmath17 is described by a stationary markovian gaussian process ( i.e , the ornstein  uhlenbeck process ) with the mean @xmath167 , variance @xmath168 , and correlation time @xmath169 . a standard numerical scheme @xcite for generating such a time - correlated gaussian process is used to model the time variation of @xmath17 for a given drop trajectory .
an example of drop behavior in a stochastic flow with a gaussian variation of the vorticity is presented in figure [ fig : cgaussdns ] .
the mean value of the vorticity is @xmath163 and the variance is @xmath164 .
the correlation time of the vorticity distribution @xmath170 is several times shorter than the drop relaxation time .
the figure depicts the probability distribution of the drop length , obtained using the boundary - integral simulations .
the capillary number is @xmath40 , and the drop viscosity is @xmath39 .
the results indicate that the drop - length distribution is bimodal for the above parameter values .
this behavior is expected since the vorticity undergoes random variation in the domain that includes both the lower and upper critical values @xmath60 , and @xmath56 of the vorticity parameter @xmath17 . since drop response to a flow with slow variation of vorticity
is hysteretic , a drop in the random flow tends to stay in the neighborhood of the compact and the elongated stationary states .
we note that the peak of the length probability distribution at @xmath171 is shifted towards the shorter drop lengths compared to the length of a drop in the elongated stationary shape .
this is because @xmath40 is close to the critical capillary number for drop breakup .
thus , due to the slow time scale in the drop dynamics near the critical capillary number @xcite , the drop does not have sufficient time to fully extend before the vorticity significantly changes .
figure [ fig : cgauss01 ] shows the probability density distribution for the drop length @xmath35 and drop angle @xmath36 for different values of the flow correlation time @xmath169 .
other system parameters are the same as in figure [ fig : cgaussdns ] .
the calculations were performed using the small - deformation equations .
the results indicate that at short flow correlation times the drop - length probability distribution is peaked around small values corresponding to the short - drop stationary solution , and the has a moderate - hight peak at @xmath172 .
as the flow correlation time @xmath169 increases , the length probability distribution becomes bimodal : one of its peaks corresponds to compact and the other to elongated drops .
a corresponding change occurs in the angle distribution , i.e. , its peak becomes more pronounced , and shifts towards the straining axis @xmath43 .
the shift of the typical drop length and orientation from the compact to elongated state when the flow correlation time is increased resembles the analogous shift for a system with harmonic vorticity oscillations ( see figures [ fig : chaos01 ] and figures [ fig : chaos01 ] ) .
this behavior is further illustrated in figure [ fig : cgauss02 ] which shows the average values and the variance of @xmath35 and @xmath36 versus the correlation time @xmath169 .
normalized by drop deformation time , for @xmath39 , @xmath40 , @xmath163 , and @xmath164 .
, title="fig:",width=249,height=211 ]   normalized by drop deformation time , for @xmath39 , @xmath40 , @xmath163 , and @xmath164 .
, title="fig:",width=249,height=211 ]
we have presented results of numerical and theoretical investigations of the dynamics of highly viscous drops in two - dimensional linear creeping flows with time - dependent vorticity . in our earlier publication @xcite we predicted that in stationary flows such drops exhibit bistable behavior : there is a range of system parameters where the drop may assume either an elongated shape approximately aligned with the straining axis of the external flow or a nearly spherical shape approximately aligned in the flow direction . here
we analyze the consequences of this behavior for the system dynamics .
we also elucidate the physical mechanism that leads to drop bistability .
a direct consequence of the existence of two stationary states is hysteretic drop response to a flow with slowly varying vorticity .
we have explained that the rapid transition from an elongated non - rotating drop shape to the nearly spherical compact shape occurs when the vorticity becomes strong enough to overcome the effect of the straining flow component that aligns the drop with the straining axis .
this transition is thus analogous to the behavior of an elongated rigid particle which starts to tumble when the vorticity grows above a critical value @xcite .
a viscous drop also begins to tumble at a critical vorticity magnitude @xmath56 .
however , when the drop becomes misaligned with the extensional axis it relaxes under the action of capillary forces towards a nearly spherical rotationally stabilized stationary shape .
if , in turn , the vorticity is slowly decreased , the drop returns to the elongated shape only after the vorticity magnitude reaches a lower critical value @xmath173 .
this hysteretic drop response occurs at high drop viscosities because the rotational stabilizing mechanism is more efficient in the high - viscosity regime .
a highly viscous drop deforms less within each drop revolution , so the compact shape remains stable even for small vorticity magnitudes .
the existence of two stationary states affects drop dynamics not only in the quasistatic regime but also at finite frequencies of the external forcing . at small amplitudes of harmonic vorticity oscillations
the drop simply oscillates ( with the same frequency as the external forcing ) about one of the stationary states .
however , if the vorticity - variation range includes both critical values @xmath60 and @xmath56 the dynamics of the system is much richer . we find that with an increasing magnitude of the vorticity oscillations the system undergoes a cascade of period - doubling bifurcations resulting in chaotic drop dynamics . the period doubling stems from the resonance between the periodicity of the external forcing and the tumbling motion of the drop when it jumps from a ( partially ) elongated shape towards the compact rotationally stabilized state .
chaos in our system emerges despite linearity of stokes equations
the system dynamics is nonlinear because of the coupling of the flow to the evolving fluid interface .
a detailed analysis of small - deformation equations describing drop dynamics reveals that in addition to chaos associated with the period - doubling mechanism there also exists in our system a different kind of chaotic evolution that results from manifold tangling @xcite .
our analysis of different types of chaos will be presented in a separate publication .    to our knowledge ,
chaotic drop dynamics in stokes - flow regime has never been observed before .
we note , however , that in a recent independent study @xcite has reported chaotic dynamics in co - annular stokes flow with insoluble surfactant adsorbed on the fluid interface .
chaos in their system also appears as a result of period doubling .
understanding of drop bistability and the associated dynamical phenomena is relevant for many practical problems .
for example , interpretation of rheological response of emulsions of highly viscous drops to time - varying flows requires insight into drop dynamics .
our results may also be useful in design new methods for manipulating emulsion microstructure in material processing and controlling drop behavior in microfluidic flows .
drop bistability could , e.g. , be used to construct microfluidic switches , and chaotic drop dynamics may be relevant for microfluidic mixing .
in particular , our analysis suggests that macromolecules with high degree of internal dissipation may undergo a transition between a nearly spherical and moderately elongated states ( in addition to the standard coil - stretch transition predicted by * ? ? ? * ) .
in fact , the dynamics of macromolecules can be modeled using equations analogous to , supplemented with terms representing random thermal forces .
such a simplified description correctly captures the most important features of power spectra of dna molecules evolving in linear flows with nonzero rotational component @xcite .
there are also close analogies between drop and vesicle motion .
the main difference between these two systems is that vesicles satisfy a constant - area constraint whereas the drop area can vary .
this constraint gives rise to periodic vesicle motion ( such as tank treading and tumbling ) even in stationary flows @xcite .
it would be interesting to determine if a coupling of vesicle oscillations to a harmonic variation of the external flow can lead to chaotic dynamics .
it would also be of significant interest to experimentally explore the bistable and chaotic drop dynamics ( as well as related phenomena for other deformable particles ) .
in such experiments a four - roll mill could be used to produce a linear flow with a controlled magnitude of vorticity .
the experiments could also be performed using recently developed microfluidic analogues of a four - roll mill device @xcite .
we would like to acknowledge helpful discussions with paul steen , demetrios papageorgiou , and michael loewenberg .
we thank michael loewenberg for permission to use his improved boundary - integral code , and also petia vlahovska for use of her mathematica codes for the coefficients in the small - deformation theory .
jb was supported by nsf career grant cts-0348175 , yny acknowledges a nsf / dms grant ( dms-0708977 ) and a sbr grant from njit .
the simulations were conducted on the njit computer cluster supported by nsf / mri grant dms-0420590 .
1983 nonlinear oscillations , dynamical systems , and bifurcations of vector fields . in _ nonlinear oscillations , dynamical systems , and bifurcations of vector fields _
f.  john , j.  e. marsden & l.  sirovich ) .
new york : springer - verlag .