 momentum and kinetic-energy budgets. Further, the drag force for the perpendicular motion shows quantitative agreement with the infinite series solution of [@bar1968] with an error of less than $0.043\%$. The drag force for the parallel motion approaches the wide-gap solution of [@mag2003] with increasing $\varepsilon$. To make sure of numerical stability and accuracy, we set the clearance parameter in a range of $10^{-3}\leq \varepsilon\leq 9$. We performed the convergence tests with varying the size $R_{\rm max}$ of the computational domain, the number $N_\xi$ of nodes in the gap between the bubble and the wall, and the number $N_\eta$ of nodes describing the bubble surface. We confirmed that the relative errors in the migration velocity $V_{B3}$ and the drag forces for the perpendicular and parallel bubble motions to those obtained on the base meshes $N_\xi\times N_\eta = 200 \times 200$ for various $\varepsilon$ decrease with increasing the size $R_{\rm max}$ and the resolutions $N_\xi$ and $N_\eta$. From the convergence behaviors, we deduce to obtain the migration velocities $V_{B3}$ with errors of much less than $1\%$ using the present base meshes, which are accurate enough for the subsequent discussion.

Migration velocity
------------------

As shown in the inset of figure \[fig:migvel\_comp\_exp\], the ratio of the simulation to the analytical migration velocity becomes close to unity as $\kappa (=(1+\varepsilon)^{-1})$ approaches zero, indicating the simulation result is consistent with the analytical solution ([@mag2003]) for small $\kappa (< \sim 0.3)$. Such a consistency between the different approaches may deny the erroneous prediction of the migration velocity in [@mag2003], of which the possibility was pointed out by [@wan2006]. By contrast, the simulated migration velocity for the bubble closer to the wall with the clearance shorter than the bubble radius (i.e., for $\kappa > \sim 0.5$ presumably beyond the theoretical precondition $\kappa\ll 1$) is considerably higher than the analytical solution. Although the simulation result reveals lower velocity than the experimental one, the tendency of the higher velocity than the theoretical one for the narrow gap is qualitatively similar to that in the experiment. Thus, the present simulation also indicates the presence of the additional narrow-gap repulsive force.

For an undeformed spherical bubble at small but non-zero Reynolds numbers $0
momentum and kinetic - energy budgets. Further, the puff power for the perpendicular motion testify quantitative agreement with the countless series solution of [ @bar1968 ] with an error of less than $ 0.043\%$. The drag force out for the parallel motion approaches the across-the-board - gap solution of [ @mag2003 ] with increasing $ \varepsilon$. To take sure of numeral constancy and accuracy, we set the clearance argument in a range of $ 10^{-3}\leq \varepsilon\leq 9$. We perform the overlap tests with varying the size $ R_{\rm max}$ of the computational domain, the phone number $ N_\xi$ of nodes in the gap between the house of cards and the wall, and the phone number $ N_\eta$ of nodes identify the bubble surface. We confirmed that the proportional errors in the migration velocity $ V_{B3}$ and the drag forces for the perpendicular and parallel bubble apparent motion to those obtained on the base meshes $ N_\xi\times N_\eta = 200 \times 200 $ for various $ \varepsilon$ decrease with increasing the size $ R_{\rm max}$ and the resolutions $ N_\xi$ and $ N_\eta$. From the convergence behaviors, we deduce to obtain the migration velocities $ V_{B3}$ with mistake of much less than $ 1\%$ using the present infrastructure mesh topology, which are accurate enough for the subsequent discussion. 

 Migration velocity 
 ------------------ 

 As shown in the inset of digit \[fig: migvel\_comp\_exp\ ], the proportion of the simulation to the analytical migration velocity become close to unity as $ \kappa (= (1+\varepsilon)^{-1})$ approaches zero, indicating the simulation result is consistent with the analytical solution ([ @mag2003 ]) for small $ \kappa (< \sim 0.3)$. Such a consistency between the different approaches may traverse the erroneous prediction of the migration velocity in [ @mag2003 ], of which the hypothesis was pointed out by [ @wan2006 ]. By contrast, the simulated migration velocity for the bubble closer to the rampart with the clearance shorter than the bubble radius (i.e., for $ \kappa > \sim 0.5 $ presumably beyond the theoretical precondition $ \kappa\ll 1 $) is considerably higher than the analytical solution. Although the simulation solution reveals lower speed than the experimental one, the tendency of the higher velocity than the theoretical one for the narrow opening is qualitatively similar to that in the experiment. therefore, the present simulation also indicates the bearing of the additional narrow - gap repulsive force. 

 For an undeformed spherical bubble at small but non - zero Reynolds numbers $ 0
 molentum and kinetic-energy budgets. Furthet, rhe dreg forcs for thd perpendicular motion shows qyantiuctive agreement with the infijite serues wilution of [@bar1968] witm an svror mh less than $0.043\%$. Thg drag force for the paralneu lotion approaches the wide-gap solutyon of [@kah2003] with increasyng $\nawepsjlon$. To make sure of numerical stagility end accuracy, we set the clearance parametfr ij a range of $10^{-3}\leq \vwrepsilon\lew 9$. Wq performed tfe convergtnee tests wijh varying the size $R_{\rm max}$ of tfe cokputationao eomwhn, the numbxr $N_\xi$ of nodes in the gap betweem the bubble akd thx waol, and the number $N_\ete$ of nodes describind the bubtlz surface. We confirmee rhat jhe ralatkce drrkrx jn the mijration velkcity $V_{B3}$ ane the drag forces fpr nne perpendichlar agd parallel bubble motions to those obtaited on the base meshes $N_\xi\rimes N_\eta = 200 \times 200$ flr variouf $\varepsilon$ decrease with increasing the size $R_{\rk max}$ xnd tme rdwoputions $N_\xi$ and $N_\eta$. From the convergence behwbiprx, we deduce to obtain the mlgtwtion velocitkes $V_{B3}$ wifh errors of much pess thwn $1\%$ uwing the kresemt base meshes, which are acxurate enougk fir the subsequent biscussion.

Miyratiom velpcity
------------------

As shown in the inret kf figure \[flg:migvel\_ckop\_exp\], the ratio uf nhe vimulation to the analyticwl migratmon vzlocity cecokes cljse to unihy as $\kappa (=(1+\varepsilon)^{-1})$ aoproaehes sero, indicwting the simulation result is rpnsistent wijh dhe analytieal sokution ([@mag2003]) fow small $\kappa (< \sim 0.3)$. Sbch a zonsistencj between the differqnt approachev may deny thx erroneots peedixtion ow the migration velocity in [@mag2003], of which the possibikiti sas pointed out by [@wan2006]. By contrasy, tfe filukaead migration velucigu for the bubble zlosrr to the wall with dhe dlearance shorter yhwn the bobble radyus (i.e., for $\kalpa > \sim 0.5$ presumabpy bebond tie theprejical precondition $\kappa\ll 1$) is donsiderahly higher than the analytical xolution. Although the simulation result reveals lower velicity than the expetimtntal one, the tendqncy of tve higher velocity tyan the theoreticsl one for the narrow gap iv quapitatively similar to that in the experiment. Thus, the present simulation qlso iidycates the lresrnce mf che addytioiau narrow-gap repuksive force.

For an undeformed spierical butbje at small but non-zero Reynplas numbers $0
momentum and kinetic-energy budgets. Further, the drag the motion shows agreement with the with error of less $0.043\%$. The drag for the parallel motion approaches the solution of [@mag2003] with increasing $\varepsilon$. To make sure of numerical stability and we set the clearance parameter in a range of $10^{-3}\leq \varepsilon\leq 9$. We the tests varying size $R_{\rm max}$ of the computational domain, the number $N_\xi$ of nodes in the gap between bubble and the wall, and the number $N_\eta$ nodes describing the bubble We confirmed that the relative in migration velocity and drag for the perpendicular parallel bubble motions to those obtained on the base meshes $N_\xi\times N_\eta = 200 \times 200$ for $\varepsilon$ decrease the size max}$ the $N_\xi$ and $N_\eta$. convergence behaviors, we deduce to obtain $V_{B3}$ with errors of much less than $1\%$ the present meshes, which are accurate enough for subsequent discussion. Migration velocity ------------------ As shown in inset of figure \[fig:migvel\_comp\_exp\], the ratio of the simulation to the analytical migration velocity becomes unity as $\kappa (=(1+\varepsilon)^{-1})$ zero, indicating the result consistent the solution ([@mag2003]) small $\kappa (< \sim 0.3)$. Such a consistency between the different may deny the erroneous prediction of the migration velocity in which possibility was pointed by [@wan2006]. By contrast, simulated velocity for the bubble the with than bubble (i.e., for $\kappa > 0.5$ presumably beyond the theoretical $\kappa\ll 1$) is considerably Although the simulation result reveals lower velocity than experimental one, the tendency of the higher than the theoretical one for the narrow gap is qualitatively similar to in the the present simulation also indicates the presence of additional narrow-gap repulsive force. an undeformed spherical bubble at small but non-zero Reynolds $0
 momentum and kinetic-energy bUdgets. FurtHer, thE drAg fOrCe foR the Perpendicular mOTion Shows quantitative agreeMent wItH The iNFiNite sEries soLUtION of [@BaR1968] wIth An ERrOr of lEss Than $0.043\%$. The Drag force fOr tHe Parallel motiON aPproaches tHe wIde-gap solutiOn oF [@mag2003] wiTh IncREasinG $\vaRepsiLon$. To mAKe sure Of numericAl STabiliTY and accURAcY, we sEt the clearance parAMeTEr in a range of $10^{-3}\leQ \varepSiLOn\LEQ 9$. We PerFormed the cOnVergeNCe tests WItH VARyiNG the size $R_{\rm maX}$ of the compuTAtiOnal doMaIn, tHE numbeR $N_\xi$ oF nODes In the gap betWeen The bubble And the WAll, and tHE number $n_\eta$ of NodEs dEscrIBiNg The BuBBle SUrFacE. we cOnfirmed ThAt The reLatiVE ERRors In tHe miGratiOn velocity $V_{B3}$ aNd tHe drAG foRces fOr the PerpEnDiculAr and pArallEl Bubble motions to ThosE obtained On tHe BasE mEshes $n_\Xi\timeS N_\eTa = 200 \tImes 200$ for Various $\VArePsILON$ dEcrease with increasInG THe Size $R_{\rm mAx}$ and tHE rEsOLutions $N_\Xi$ And $n_\eta$. fROm the ConvERgEnce behaViors, wE DeDuCe to obtAiN the miGrAtiOn vElociTIes $V_{b3}$ with eRrors of mUch leSS than $1\%$ using the pREsent base meshES, wHICh ARe acCurAte enough foR the SUbseQuenT DiScuSSion.

MIgratIoN VeLOcity
------------------

As shown in the inSeT of figUre \[fiG:migvel\_comp\_exP\], the ratio oF THE simulatIon tO ThE Analytical migrAtion Velocity beCOmes closE to unIty as $\kapPa (=(1+\varepsiLON)^{-1})$ approacHes ZerO, inDicATInG the simulatioN REsulT iS consisTenT with thE anAlyTicAl sOlUtion ([@mag2003]) fOr small $\kApPa (< \SiM 0.3)$. SUch A consIStency beTwEen ThE diFfereNT approAches May dEnY tHE erRoneous PReDICtioN oF tHe miGraTiOn velOcitY In [@mAg2003], of whiCh the possIbiLIty wAs PoInted ouT by [@wan2006]. By contrAsT, the simulaTeD miGratioN VElocity fOr the bubble closer to the wALl with tHe cLearaNce sHorter thaN thE bubblE raDIus (i.e., fOr $\kappA > \sim 0.5$ pReSumABLy beyOND tHe tHeOretical prECOndItion $\KaPpa\lL 1$) is consIderably higher than THe aNalytical soluTioN. AltHOUgH thE SiMUlaTiON reSULt reveals lower vElocity thaN tHE eXperimentaL One, ThE tendenCy of the HigheR VelocitY than the tHeoreticaL oNe foR THe nArrow gap is QualitatIvely simiLAr to tHAt In the ExpErimenT. THus, The prEsent sIMulAtion Also inDiCates tHe preSeNce of the Additional narrow-gap repuLsive fOrce.

FOr aN undeformEd sPHerIcal bubblE at sMall but non-ZerO ReYnoldS nuMBers $0
 momentum and kinetic-ener gy budgets . Fur the r,th e dr ag f orce for the p e rpen dicular motion shows q uanti ta t ivea gr eemen t witht he i nfi ni te  se ri e ssolut ion  of [@b ar1968] wi than  error of le s sthan $0.04 3\% $. The dragfor ce for  t hep arall elmotio n appr o achesthe wide- ga p  solut i on of [ @ m ag 2003 ] with increasing $\ v arepsilon$. To  makesu r eo f  nu mer ical stabi li ty an d  accura c y, w e  se t  the clearanc e parameter ina rang eof$ 10^{-3 }\leq  \ v are psilon\leq9$.We perfor med th e  conver g ence te sts wi thvar ying th esiz e$ R_{ \ rm  ma x }$of the c om pu tatio nald o m a in,the  num ber $ N_\xi$ of nod esin t h e g ap be tween  the  b ubble  and t he wa ll , and the numbe r $N _\eta$ of  no de s d es cribi n g thebub ble  surfac e. We c o nfi rm e d th at the relative er ro r s  i n the mi gratio n  v el o city $V_ {B 3}$  and t he dr ag f o rc es for t he per p en di cular a nd  paral le l b ubb le mo t ions  to th ose obta inedo n the base mes h es $N_\xi\tim e sN _ \e t a =200  \times 200 $ fo r  var ious $\ var e psilo n$ de cr e as e  with increasing th esize $ R_{\r m max}$ and t he resolut i o n s $N_\xi $ an d  $ N _\eta$. From t he co nvergenceb ehaviors , wededuce t o obtaint h e migrat ion  ve loc iti e s  $ V_{B3}$ withe r rors  o f muchles s than$1\ %$usi ngth e present  base me sh es ,wh ich  area ccurateen oug hfor  thes ubsequ ent d iscu ss io n .
Migrati o nv e loci ty 
- ---- --- -- ----- ---Asshown i n the ins eto f fi gu re  \[fig: migvel\_comp\ _e xp\], thera tio  of th e simulati on to the analytical mi g rationvel ocity  bec omes clos e t o unit y a s  $\kap pa (=( 1+\va re psi l o n)^{- 1 } )$  ap pr oaches zer o ,  in dicat in g th e simul ation result is co n sis tent with the  an alyt i c al  so l ut i on([ @ mag 2 0 03]) for small$\kappa (<  \ s im  0.3)$. Su c h a  c onsiste ncy bet weent he diff erent app roaches m ay  den y the  erroneous  predict ion of th e  migr a ti on ve loc ity in  [ @ma g2003 ], ofw hic h the  possi bi lity w as po in ted outby [@wan2006]. By contr ast, t he si mul ated migr ati o n v elocity f or t he bubbleclo ser  to t hew all w itht he  cl e aranc e sh o rter than th e b u b bl e radius (i . e . , f or $\ kap p a > \s im 0 .5$ presumably be y ond the theore tica l pre con d itio n$\kappa\ll 1$)  is  c o n siderabl yhigher than  the ana ly t icalsoluti on. Al thought h es imulat ionres ult revea lslo w er velo ci ty than t he e xp erimen tal on e , th e tendency of thehighe r veloc i tythanth e theor e tica l one forthe narrowgap is  qua litat ively s im ilar t o t ha t in the e x periment.  Thus , the p re sent  si mulati on a l s o ind icat es  th e presenc e of th ea ddi tion al na rr ow-g ap repuls i ve force .
F or an u nd efo r m ed sph e ri c a l bubble a t s mallb u t non-zero Reyn o ld s  numb ers $0 
 momentum_and kinetic-energy_budgets. Further, the drag_force for_the_perpendicular motion_shows_quantitative agreement with_the infinite series_solution of [@bar1968] with_an error of_less_than $0.043\%$. The drag force for the parallel motion approaches the wide-gap solution of_[@mag2003]_with increasing_$\varepsilon$._To_make sure of numerical stability_and accuracy, we set the_clearance parameter_in a range of $10^{-3}\leq \varepsilon\leq 9$. We_performed_the convergence tests_with varying the size $R_{\rm max}$ of the computational_domain, the number $N_\xi$ of nodes_in the gap_between_the_bubble and the wall,_and the number $N_\eta$ of nodes_describing the bubble surface. We confirmed_that the relative errors in the migration_velocity $V_{B3}$ and the drag forces_for the perpendicular and parallel_bubble motions_to those obtained on the_base meshes $N_\xi\times_N_\eta =_200 \times 200$_for various $\varepsilon$ decrease with increasing_the size $R_{\rm_max}$ and the resolutions $N_\xi$ and_$N_\eta$._From the convergence_behaviors,_we_deduce to_obtain the migration_velocities_$V_{B3}$ with_errors_of much less than $1\%$ using_the_present base meshes, which are accurate enough_for the subsequent discussion.

Migration_velocity
------------------

As_shown in the inset_of figure \[fig:migvel\_comp\_exp\], the ratio_of the simulation to the analytical_migration velocity_becomes close_to unity as $\kappa (=(1+\varepsilon)^{-1})$ approaches zero, indicating the simulation result_is consistent with the analytical solution_([@mag2003]) for small $\kappa_(< \sim_0.3)$._Such a consistency_between_the different_approaches may deny the erroneous prediction of_the migration_velocity in [@mag2003], of which the_possibility was pointed out_by_[@wan2006]. By contrast, the simulated migration_velocity for the bubble closer to_the wall with the clearance_shorter_than_the bubble radius (i.e., for_$\kappa > \sim 0.5$ presumably beyond_the theoretical precondition_$\kappa\ll 1$) is considerably higher than the_analytical_solution. Although the simulation result reveals_lower_velocity than the experimental one, the_tendency_of_the higher velocity than the_theoretical one for the narrow gap_is qualitatively similar to that in the experiment. Thus,_the present simulation_also indicates the presence of_the_additional_narrow-gap repulsive force.

For an undeformed spherical bubble at small but_non-zero Reynolds_numbers $0