Field-only surface integral equations: scattering from a perfect electric conductor.

A field-only boundary integral formulation of electromagnetics is derived without the use of surface currents that appear in the Stratton-Chu formulation. For scattering by a perfect electrical conductor (PEC), the components of the electric field are obtained directly from surface integral equation solutions of three scalar Helmholtz equations for the field components. The divergence-free condition is enforced via a boundary condition on the normal component of the field and its normal derivative. Field values and their normal derivatives at the surface of the PEC are obtained directly from surface integral equations that do not contain divergent kernels. Consequently, high-order elements with fewer degrees of freedom can be used to represent surface features to a higher precision than the traditional planar elements. This theoretical framework is illustrated with numerical examples that provide further physical insight into the role of the surface curvature in scattering problems.

Field-only surface integral equations: scattering from a dielectric body.

An efficient field-only nonsingular surface integral method to solve Maxwell's equations for the components of the electric field on the surface of a dielectric scatterer is introduced. In this method, both the vector wave equation and the divergence-free constraint are satisfied inside and outside the scatterer. The divergence-free condition is replaced by an equivalent boundary condition that relates the normal derivatives of the electric field across the surface of the scatterer. Also, the continuity and jump conditions on the electric and magnetic fields are expressed in terms of the electric field across the surface of the scatterer. Together with these boundary conditions, the scalar Helmholtz equation for the components of the electric field inside and outside the scatterer is solved by a fully desingularized surface integral method. Compared with the most popular surface integral methods based on the Stratton-Chu formulation or the Poggio-Miller-Chew-Harrington-Wu-Tsai (PMCHWT) formulation, our method is conceptually simpler and numerically straightforward because there is no need to introduce intermediate quantities such as surface currents, and the use of complicated vector basis functions can be avoided altogether. Also, our method is not affected by numerical issues such as the zero-frequency catastrophe and does not contain integrals with (strong) singularities. To illustrate the robustness and versatility of our method, we show examples in the Rayleigh, Mie, and geometrical optics scattering regimes. Given the symmetry between the electric field and the magnetic field, our theoretical framework can also be used to solve for the magnetic field.

On the derivation of the Smoluchowski result of electrophoretic mobility.

One of the most enduring, broadly applicable and widely used theoretical results of electrokinetic theory is the Smoluchowski expression for the electrophoretic mobility. It is a limiting form that holds for any solid particle of arbitrary shape in an electrolyte of any composition provided the thickness of the electrical double layer is "infinitely" thin compared to the particle size and the particle has uniform surface potential. The familiar derivation of this result that is a simplified version of the original Smoluchowski analysis in 1903, considers the motion of the electrolyte adjacent to a planar surface. The theory is deceptively simple but as a result much of the interesting physics and characteristic hydrodynamic behavior around the particle have been obscured, leading to a significantly incorrect picture of the fluid velocity profile near the particle surface. This paper provides a derivation of this key theoretical result by starting from Smoluchowski's original 1903 analysis but brings out overlooked details of the hydrodynamic features near and far from the particle that have not been canvassed in detail. The objective is to draw together all the key physical features of the electrophoretic problem in the thin double layer regime to provide an accessible and complete exposition of this important result in colloid science.

Mobile-surface bubbles and droplets coalesce faster but bounce stronger.

Enhancing the hydrodynamic interfacial mobility of bubbles and droplets in multiphase systems is expected to reduce the characteristic coalescence times and thereby affect the stability of gas or liquid emulsions that are of wide industrial and biological importance. However, by comparing the controlled collision of bubbles or water droplets with mobile or immobile liquid interfaces, in a pure fluorocarbon liquid, we demonstrate that collisions involving mobile surfaces result in a significantly stronger series of rebounds before the rapid coalescence event. The stronger rebound is explained by the lower viscous dissipation during collisions involving mobile surfaces. We present direct numerical simulations to confirm that the observed rebound is enhanced with increased surface mobility. These observations require a reassessment of the role of surface mobility for controlling the dynamic stability of gas or liquid emulsion systems relevant to a wide range of processes, from microfluidics and pharmaceuticals to food and crude oil processing.

The unusual fluid dynamics of particle electrophoresis.

The classical problem of the electrophoretic motion of a spherical particle has been treated theoretically by Overbeek in his 1941 PhD thesis and almost 40 years later by O'Brien & White. Although both approaches used identical assumptions, the details are quite different. Overbeek solved for the pressure, velocity fields as well as the electrostatic potential, whereas O'Brien & White obtained the electrophoretic mobility without the need to consider the pressure and velocity explicitly. In this paper, we establish the equivalence of these two approaches that allow us to show that the tangential component of the fluid velocity has a maximum near the surface of the particle and outside the double layer, the velocity decays as 1/r3, where r is the distance from the sphere, instead of 1/r in normal Stokes flow. Associated with this behavior is that of an irrotational outer flow field. This is consistent with the fact that a sphere moving with a constant electrophoretic velocity experiences zero net force. A study of the forces on the particle also provides a physical explanation of the independence of the electrophoretic mobility on the electrostatic boundary conditions or dielectric permittivity of the particle. These results are important in situations where inter-particle interaction is considered, for instance, in electrokinetic deposition.

Coalescence Dynamics of Mobile and Immobile Fluid Interfaces.

Coalescence dynamics between deformable bubbles and droplets can be dramatically affected by the mobility of the interfaces with fully tangentially mobile bubble-liquid or droplet-liquid interfaces expected to accelerate the coalescence by orders of magnitude. However, there is a lack of systematic experimental investigations that quantify this effect. By using high speed camera imaging we examine the free rise and coalescence of small air-bubbles (100 to 1300 µm in diameter) with a liquid interface. A perfluorocarbon liquid, PP11, is used as a model liquid to investigate coalescence dynamics between fully mobile and immobile deformable interfaces. The mobility of the bubble surface was determined by measuring the terminal rise velocity of small bubbles rising at Reynolds numbers, Re, less than 0.1 and the mobility of free PP11 surface by measuring the deceleration kinetics of the small bubble toward the interface. Induction or film drainage times of a bubble at the mobile PP11-air surface were found to be more than 2 orders of magnitude shorter compared to the case of bubble and an immobile PP11-water interface. A theoretical model is used to illustrate the effect of hydrodynamics and interfacial mobility on the induction time or film drainage time. The results of this study are expected to stimulate the development of a comprehensive theoretical model for coalescence dynamics between two fully or partially mobile fluid interfaces.

Field-only integral equation method for time domain scattering of electromagnetic pulses.

The scattering of electromagnetic pulses is described using a non-singular boundary integral method to solve directly for the field components in the frequency domain, and Fourier transform is then used to obtain the complete space-time behavior. This approach is stable for wavelengths both small and large relative to characteristic length scales. Amplitudes and phases of field values can be obtained accurately on or near material boundaries. Local field enhancement effects due to multiple scattering of interest to applications in microphotonics are demonstrated.

Self-determined shapes and velocities of giant near-zero drag gas cavities.

Minimizing the retarding force on a solid moving in liquid is the canonical problem in the quest for energy saving by friction and drag reduction. For an ideal object that cannot sustain any shear stress on its surface, theory predicts that drag force will fall to zero as its speed becomes large. However, experimental verification of this prediction has been challenging. We report the construction of a class of self-determined streamlined structures with this free-slip surface, made up of a teardrop-shaped giant gas cavity that completely encloses a metal sphere. This stable gas cavity is formed around the sphere as it plunges at a sufficiently high speed into the liquid in a deep tank, provided that the sphere is either heated initially to above the Leidenfrost temperature of the liquid or rendered superhydrophobic in water at room temperature. These sphere-in-cavity structures have residual drag coefficients that are typically less than [Formula: see text] those of solid objects of the same dimensions, which indicates that they experienced very small drag forces. The self-determined shapes of the gas cavities are shown to be consistent with the Bernoulli equation of potential flow applied on the cavity surface. The cavity fall velocity is not arbitrary but is uniquely predicted by the sphere density and cavity volume, so larger cavities have higher characteristic velocities.

Space-time domain solutions of the wave equation by a non-singular boundary integral method and Fourier transform.

The general space-time evolution of the scattering of an incident acoustic plane wave pulse by an arbitrary configuration of targets is treated by employing a recently developed non-singular boundary integral method to solve the Helmholtz equation in the frequency domain from which the space-time solution of the wave equation is obtained using the fast Fourier transform. The non-singular boundary integral solution can enforce the radiation boundary condition at infinity exactly and can account for multiple scattering effects at all spacings between scatterers without adverse effects on the numerical precision. More generally, the absence of singular kernels in the non-singular integral equation confers high numerical stability and precision for smaller numbers of degrees of freedom. The use of fast Fourier transform to obtain the time dependence is not constrained to discrete time steps and is particularly efficient for studying the response to different incident pulses by the same configuration of scatterers. The precision that can be attained using a smaller number of Fourier components is also quantified.

Leidenfrost Vapor Layers Reduce Drag without the Crisis in High Viscosity Liquids.

The drag coefficient C_{D} of a solid smooth sphere moving in a fluid is known to be only a function of the Reynolds number Re and diminishes rapidly at the drag crisis around Re∼3×10^{5}. A Leidenfrost vapor layer on a hot sphere surface can trigger the onset of the drag crisis at a lower Re. By using a range of high viscosity perfluorocarbon liquids, we show that the drag reduction effect can occur over a wide range of Re, from as low as ∼600 to 10^{5}. The Navier slip model with a viscosity dependent slip length can fit the observed drag reduction and wake shape.

A robust and accurate formulation of molecular and colloidal electrostatics.

This paper presents a re-formulation of the boundary integral method for the Debye-Hückel model of molecular and colloidal electrostatics that removes the mathematical singularities that have to date been accepted as an intrinsic part of the conventional boundary integral equation method. The essence of the present boundary regularized integral equation formulation consists of subtracting a known solution from the conventional boundary integral method in such a way as to cancel out the singularities associated with the Green's function. This approach better reflects the non-singular physical behavior of the systems on boundaries with the benefits of the following: (i) the surface integrals can be evaluated accurately using quadrature without any need to devise special numerical integration procedures, (ii) being able to use quadratic or spline function surface elements to represent the surface more accurately and the variation of the functions within each element is represented to a consistent level of precision by appropriate interpolation functions, (iii) being able to calculate electric fields, even at boundaries, accurately and directly from the potential without having to solve hypersingular integral equations and this imparts high precision in calculating the Maxwell stress tensor and consequently, intermolecular or colloidal forces, (iv) a reliable way to handle geometric configurations in which different parts of the boundary can be very close together without being affected by numerical instabilities, therefore potentials, fields, and forces between surfaces can be found accurately at surface separations down to near contact, and (v) having the simplicity of a formulation that does not require complex algorithms to handle singularities will result in significant savings in coding effort and in the reduction of opportunities for coding errors. These advantages are illustrated using examples drawn from molecular and colloidal electrostatics.

The hydrodynamics of bubble rise and impact with solid surfaces.

A bubble smaller than 1mm in radius rises along a straight path in water and attains a constant speed due to the balance between buoyancy and drag force. Depending on the purity of the system, within the two extreme limits of tangentially immobile or mobile boundary conditions at the air-water interface considerably different terminal speeds are possible. When such a bubble impacts on a horizontal solid surface and bounces, interesting physics can be observed. We study this physical phenomenon in terms of forces, which can be of colloidal, inertial, elastic, surface tension and viscous origins. Recent advances in high-speed photography allow for the observation of phenomena on the millisecond scale. Simultaneous use of such cameras to visualize both rise/deformation and the dynamics of the thin film drainage through interferometry are now possible. These experiments confirm that the drainage process obeys lubrication theory for the spectrum of micrometre to millimetre-sized bubbles that are covered in this review. We aim to bridge the colloidal perspective at low Reynolds numbers where surface forces are important to high Reynolds number fluid dynamics where the effect of the surrounding flow becomes important. A model that combines a force balance with lubrication theory allows for the quantitative comparison with experimental data under different conditions without any fitting parameter.

The impact and bounce of air bubbles at a flat fluid interface.

The rise and impact of bubbles at an initially flat but deformable liquid-air interface in ultraclean liquid systems are modelled by taking into account the buoyancy force, hydrodynamic drag, inertial added mass effect and drainage of the thin film between the bubble and the interface. The bubble-surface interaction is analyzed using lubrication theory that allows for both bubble and surface deformation under a balance of normal stresses and surface tension as well as the long-range nature of deformation along the interface. The quantitative result for collision and bounce is sensitive to the impact velocity of the rising bubble. This velocity is controlled by the combined effects of interfacial tension via the Young-Laplace equation and hydrodynamic stress on the surface, which determine the deformation of the bubble. The drag force that arises from the hydrodynamic stress in turn depends on the hydrodynamic boundary conditions on the bubble surface and its shape. These interrelated factors are accounted for in a consistent manner. The model can predict the rise velocity and shape of millimeter-size bubbles in ultra-clean water, in two silicone oils of different densities and viscosities and in ethanol without any adjustable parameters. The collision and bounce of such bubbles with a flat water/air, silicone oil/air and ethanol/air interface can then be predicted with excellent agreement when compared to experimental observations.

Anomalous Sinking of Spheres due to Local Fluidization of Apparently Fixed Powder Beds.

The sinking of an intruder sphere into a powder bed in the apparently fixed bed regime exhibits complex behavior in the sinking rate and the final depth when the sphere density is close to the powder bed density. Evidence is adduced that the intruder sphere locally fluidizes the apparently fixed powder bed, allowing the formation of voids and percolation bubbles that facilitates spheres to sink slower but deeper than expected. By adjusting the air injection rate and the sphere-to-powder bed density ratio, this phenomenon provides the basis of a sensitive large particle separation mechanism.

Electric Double-Layer Interaction between Dissimilar Charge-Conserved Conducting Plates.

Small metallic particles used in forming nanostructured to impart novel optical, catalytic, or tribo-rheological can be modeled as conducting particles with equipotential surfaces that carry a net surface charge. The value of the surface potential will vary with the separation between interacting particles, and in the absence of charge-transfer or electrochemical reactions across the particle surface, the total charge of each particle must also remain constant. These two physical conditions require the electrostatic boundary condition for metallic nanoparticles to satisfy an equipotential whole-of-particle charge conservation constraint that has not been studied previously. This constraint gives rise to a global charge conserved constant potential boundary condition that results in multibody effects in the electric double-layer interaction that are either absent or are very small in the familiar constant potential or constant charge or surface electrochemical equilibrium condition.

Drag Moderation by the Melting of an Ice Surface in Contact with Water.

We report measurements of the effects of a melting ice surface on the hydrodynamic drag of ice-shell-metal-core spheres free falling in water at a Reynolds of number Re~2×10^{4}-3×10^{5} and demonstrate that the melting surface induces the early onset of the drag crisis, thus reducing the hydrodynamic drag by more than 50%. Direct visualization of the flow pattern demonstrates the key role of surface melting. Our observations support the hypothesis that the drag reduction is due to the disturbance of the viscous boundary layer by the mass transfer from the melting ice surface.

Force Balance Model for Bubble Rise, Impact, and Bounce from Solid Surfaces.

A force balance model for the rise and impact of air bubbles in a liquid against rigid horizontal surfaces that takes into account effects of buoyancy and hydrodynamic drag forces, bubble deformation, inertia of the fluid via an added mass force, and a film force between the bubble and the rigid surface is proposed. Numerical solution of the governing equations for the position and velocity of the center of mass of the bubbles is compared against experimental data taken with ultraclean water. The boundary condition at the air-water interface is taken to be stress free, which is consistent for bubbles in clean water systems. Features that are compared include bubble terminal velocity, bubbles accelerating from rest to terminal speed, and bubbles impacting and bouncing off different solid surfaces for bubbles that have already or are yet to attain terminal speed. Excellent agreement between theory and experiments indicates that the forces included in the model constitute the main physical ingredients to describe the bouncing phenomenon.

Boundary regularized integral equation formulation of the Helmholtz equation in acoustics.

A boundary integral formulation for the solution of the Helmholtz equation is developed in which all traditional singular behaviour in the boundary integrals is removed analytically. The numerical precision of this approach is illustrated with calculation of the pressure field owing to radiating bodies in acoustic wave problems. This method facilitates the use of higher order surface elements to represent boundaries, resulting in a significant reduction in the problem size with improved precision. Problems with extreme geometric aspect ratios can also be handled without diminished precision. When combined with the CHIEF method, uniqueness of the solution of the exterior acoustic problem is assured without the need to solve hypersingular integrals.

Measurement of surface and interfacial tension using pendant drop tensiometry.

Pendant drop tensiometry offers a simple and elegant solution to determining surface and interfacial tension - a central parameter in many colloidal systems including emulsions, foams and wetting phenomena. The technique involves the acquisition of a silhouette of an axisymmetric fluid droplet, and iterative fitting of the Young-Laplace equation that balances gravitational deformation of the drop with the restorative interfacial tension. Since the advent of high-quality digital cameras and desktop computers, this process has been automated with high speed and precision. However, despite its beguiling simplicity, there are complications and limitations that accompany pendant drop tensiometry connected with both Bond number (the balance between interfacial tension and gravitational forces) and drop volume. Here, we discuss the process involved with going from a captured experimental image to a fitted interfacial tension value, highlighting pertinent features and limitations along the way. We introduce a new parameter, the Worthington number, Wo, to characterise the measurement precision. A fully functional, open-source acquisition and fitting software is provided to enable the reader to test and develop the technique further.

Measurement and modeling on hydrodynamic forces and deformation of an air bubble approaching a solid sphere in liquids.

The interaction between bubbles and solid surfaces is central to a broad range of industrial and biological processes. Various experimental techniques have been developed to measure the interactions of bubbles approaching solids in a liquid. A main challenge is to accurately and reliably control the relative motion over a wide range of hydrodynamic conditions and at the same time to determine the interaction forces, bubble-solid separation and bubble deformation. Existing experimental methods are able to focus only on one of the aspects of this problem, mostly for bubbles and particles with characteristic dimensions either below 100 µm or above 1 cm. As a result, either the interfacial deformations are measured directly with the forces being inferred from a model, or the forces are measured directly with the deformations to be deduced from the theory. The recently developed integrated thin film drainage apparatus (ITFDA) filled the gap of intermediate bubble/particle size ranges that are commonly encountered in mineral and oil recovery applications. Equipped with side-view digital cameras along with a bimorph cantilever as force sensor and speaker diaphragm as the driver for bubble to approach a solid sphere, the ITFDA has the capacity to measure simultaneously and independently the forces and interfacial deformations as a bubble approaches a solid sphere in a liquid. Coupled with the thin liquid film drainage modeling, the ITFDA measurement allows the critical role of surface tension, fluid viscosity and bubble approach speed in determining bubble deformation (profile) and hydrodynamic forces to be elucidated. Here we compare the available methods of studying bubble-solid interactions and demonstrate unique features and advantages of the ITFDA for measuring both forces and bubble deformations in systems of Reynolds numbers as high as 10. The consistency and accuracy of such measurement are tested against the well established Stokes-Reynolds-Young-Laplace model. The potential to use the design principles of the ITFDA for fundamental and developmental research is demonstrated.