Cavity and flow measurements of reproducible bubble entrainment following drop impacts

High-speed motion pictures of air–water interface dynamics of drop impacts that reproducibly make bubbles are presented. The pictures show previously unobserved details of the phenomenon. Measurements are compared with available computational methods. Experimental and numerical results agree with each other on the overall shape of the interface and the occurrence of bubble detachment. Measurements, however, show that the cavity depth stagnates before bubble entrapment. This behavior is not predicted by simulation. Also discussed are the presence of a jet that strikes the new bubble after formation and the possible effect of droplet surface oscillations on bubble entrainment. Received: 25 April 2000 / Accepted: 26 April 2001     

On the dynamics of a rolling ball actuated by internal point masses

The motion of a rolling ball actuated by internal point masses that move inside the ball’s frame of reference is considered. The equations of motion are derived by applying Euler–Poincaré’s symmetry reduction method in concert with Lagrange–d’Alembert’s principle, which accounts for the presence of the nonholonomic rolling constraint. As a particular example, we consider the case when the masses move along internal rails, or trajectories, of arbitrary shape and fixed within the ball’s frame of reference. Our system of equations can treat most possible methods of actuating the rolling ball with internal moving masses encountered in the literature, such as circular motion of the masses mimicking swinging pendula or straight line motion of the masses mimicking magnets sliding inside linear tubes embedded within a solenoid. Moreover, our method can model arbitrary rail shapes and an arbitrary number of rails such as several ellipses and/or figure eights, which may be important for future designs of rolling ball robots. For further analytical study, we also reduce the system to a single differential equation when the motion is planar, that is, considering the motion of the rolling disk actuated by internal point masses, in which case we show that the results obtained from the variational derivation coincide with those obtained from Newton’s second law. Finally, the equations of motion are solved numerically, illustrating a wealth of complex behaviors exhibited by the system’s dynamics. Our results are relevant to the dynamics of nonholonomic systems containing internal degrees of freedom and to further studies of control of such systems actuated by internal masses.     

A continuous loading solution for spherical elastic-plastic waves of small amplitude

A solution is constructed for the wave motion induced by the application of a continuously increasing uniform pressure on the surface of a spherical cavity in an infinite elastic-plastic medium. The small amplitude linearization in conjunction with parabolic work-hardening leads to constant wave speeds in both elastically and plastically deforming regions, so that linear wave function representations can be used. The induced wave pattern is an elastic (continuous) front propagating into the medium followed by an expanding region of continuous elastic loading, in turn followed by a region of continuous plastic loading which persists to the cavity boundary. The motion of the continuous elastic-plastic interface is part of the solution to be determined, in contrast to discontinuous interfaces which move with a pre-determined wave speed. Validity of this wave pattern requires appropriate pressure increase on the boundary. Several examples are presented to illustrate different forms of interface motion.     

Capillary effects on a particle rolling on a plane surface in the presence of a thin liquid film

An experimental study has been performed of the effects of a liquid film on a particle rolling on a planar surface using a combination of laser-induced fluorescence and particle-image velocimetry. Contact angle hysteresis leads to asymmetry of the liquid meniscus, resulting in a difference in contact angle between the front and rear sections of the meniscus relative to the rolling particle. This asymmetry results in a capillary torque that resists the rolling motion of the particle. The particle rolling motion also induces a viscous transport of fluid from the front to the rear of the particle, which acts to shift the location of the contact point. The laser-induced fluorescence method is used to characterize the meniscus asymmetry and the resulting change in contact angle on the two sides of the particle. Particle-image velocimetry in various horizontal and vertical cross-sectional planes is used to examine the flow trajectories and velocity magnitude within the meniscus in the presence of rolling. All experiments are conducted at small capillary number, so that the meniscus is approximately circular in shape.     

Buoyancy-driven motion of a two-dimensional bubble or drop through a viscous liquid in the presence of a vertical electric field

The effect of an electric field on the buoyancy-driven motion of a two-dimensional gas bubble rising through a quiescent liquid is studied computationally. The dynamics of the bubble is simulated numerically by tracking the gas–liquid interface when an electrostatic field is generated in the vertical gap of the rectangular enclosure. The two phases of the system are assumed to be perfect dielectrics with constant but different permittivities, and in the absence of impressed charges, there is no free charge in the fluid bulk regions or at the interface. Electric stresses are supported at the bubble interface but absent in the bulk and one of the objectives of our computations is to quantify the effect of these Maxwell stresses on the overall bubble dynamics. The numerical algorithm to solve the free-boundary problem relies on the level-set technique coupled with a finite-volume discretization of the Navier–Stokes equations. The sharp interface is numerically approximated by a finite-thickness transition zone over which the material properties vary smoothly, and surface tension and electric field effects are accounted for by employing a continuous surface force approach. A multi-grid solver is applied to the Poisson equation describing the pressure field and the Laplace equation governing the electric field potential. Computational results are presented that address the combined effects of viscosity, surface tension, and electric fields on the dynamics of the bubble motion as a function of the Reynolds number, gravitational Bond number, electric Bond number, density ratio, and viscosity ratio. It is established through extensive computations that the presence of the electric field can have an important effect on the dynamics. We present results that show a substantial increase in the bubble’s rise velocity in the electrified system as compared with the corresponding non-electrified one. In addition, for the electrified system, the bubble shape deformations and oscillations are smaller, and there is a reduction in the propensity of the bubble to break up through increasingly larger oscillations.     

Phase-field modeling of interfacial dynamics in emulsion flows: Nonequilibrium surface tension

A weakly nonlocal phase-field model is used to define the surface tension in liquid binary mixtures in terms of the composition gradient in the interfacial region so that, at equilibrium, it depends linearly on the characteristic length that defines the interfacial width. Contrary to previous works suggesting that the surface tension in a phase-field model is fixed, we define the surface tension for a curved interface and far-from-equilibrium conditions as the integral of the free energy excess (i.e., above the thermodynamic component of the free energy) across the interface profile in a direction parallel to the composition gradient. Consequently, the nonequilibrium surface tension can be widely different from its equilibrium value under dynamic conditions, while it reduces to its thermodynamic value for a flat interface at local equilibrium. In nonequilibrium conditions, the surface tension changes with time: during mixing, it decreases as the inverse square root of time, while in the linear regime of spinodal decomposition, it increases exponentially to its equilibrium value, as nonlinear effects saturate the exponential growth. In addition, since temperature gradients modify the steepness of the concentration profile in the interfacial region, they induce gradients in the nonequilibrium surface tension, leading to the Marangoni thermocapillary migration of an isolated drop. Similarly, Marangoni stresses are induced in a composition gradient, leading to diffusiophoresis. We also review results on the nonequilibrium surface tension for a wall-bound pendant drop near detachment, which help to explain a discrepancy between our numerically determined static contact angle dependence of the critical Bond number and its sharp-interface counterpart from a static stability analysis of equilibrium shapes after numerical integration of the Young-Laplace equation. Finally, we present new results from phase-field simulations of the motion of an isolated droplet down an incline in gravity, showing that dynamic contact angle hysteresis can be explained in terms of the nonequilibrium surface tension.     

Experimental investigation of circular free-surface jet impingement quenching: Transient hydrodynamics and heat transfer

An experimental study of heat transfer during quenching of a cylindrical stainless steel test specimen has been performed. A subcooled water jet is directed onto the upward facing flat face of the cylinder. The test specimen is heated to an initial temperature slightly above 900 °C and then quenched. The resulting boiling curves and heat transfer distributions are presented for impingement velocities of 2.85 and 6.4 m/s (Re = 7900 and 18,900). High-speed imaging shows that three distinct regions on the quenched surface can be identified: an expanding circular wetted region surrounding the impinging point, annular transition zone just outside the wetting front, and a unwetted region outside this zone. The free-surface of the liquid in the wetted region is smooth in the nucleate and transition boiling regimes. The annular transition zone or the wetting front region outside the wetted region is characterized by a highly disturbed liquid-gas interface, which can be attributed to intense vapor generation. At the outer edge of the transition zone, the liquid is deflected away from the surface. The velocity of the wetting front significantly increases with the jet impact velocity, which indicates that the wetting front position is governed by the ability of the flowing liquid to transport the bubbles radially outwards from the wetted region.     

Numerical simulation of the horizontal Bridgman growth. Part III: Calculation of the interface

We study the transient motion of the solidification front during the growth of semiconductor crystals in the horizontal Bridgman geometry. The calculation is based on a two-dimensional flow. We use finite elements which deform with the motion of the interface. The energy equation is coupled with the isothermal constraint of the interface in an implicit transient algorithm. Several examples show the oscillatory motion of the interface caused by the periodic flow of the melt, and they reveal the importance of the growth rate on the shape of the interface.     

Piecewise analysis of oblique vibro-impacting systems

A great number of studies have shown the complex nonlinear dynamics of mechanical systems with repeated normal impacts. An oblique frictional impact introduces even more complicated dynamics such as stick-slip motions to those systems. Hence, the dynamics of oblique vibro-impacting systems with possible sliding motion is an open problem. Based on a hybrid analysis of vibro-impact dynamics, kinematics and complementary conditions, a piecewise analysis method is developed in the paper to describe the sliding motion during an oblique impact. Thereby, a parametrically excited planar pendulum between two parallel rigid walls is studied as an illustrative example. The example, together with the corresponding numerical results, shows that the sliding impacts occur in such a system with a set of properly selected parameters. The project supported by the National Natural Science Foundation of China (59905010) and by the Scientific Research Foundation for the Returned Overseas Chinese Scholars, Ministry of Education, China     

Steady motion of a rigid disk of finite thickness on a horizontal plane

The article discusses the steady motion of a rigid disk of finite thickness rolling on its edge on a horizontal plane under the influence of gravity. The governing equations are presented and two cases allowing for a steady-state solution are considered: rolling on consistently rough ground and rolling on perfectly smooth ground. The conditions of steady motion are derived for both kinds of ground and it is shown that the possible steady motion of a disk is either on a straight line or in a circle. Oscillations about steady state are discussed and conditions for stable motion established. The bifurcations of steady motions on a smooth surface are also considered.     

Hybrid dynamics of two coupled oscillators that can impact a fixed stop

We consider two linearly coupled masses, where one mass can have inelastic impacts with a fixed, rigid stop. This leads to the study of a two degree of freedom, piecewise linear, frictionless, unforced, constrained mechanical system. The system is governed by three types of dynamics: coupled harmonic oscillation, simple harmonic motion and discrete rebounds. Energy is dissipated discontinuously in discrete amounts, through impacts with the stop. We prove the existence of a non-zero measure set of orbits that lead to infinite impacts with the stop in a finite time. We show how to modify the mathematical model so that forward existence and uniqueness of solutions for all time is guaranteed. Existence of hybrid periodic orbits is shown. A geometrical interpretation of the dynamics based on action coordinates is used to visualize numerical simulation results for the asymptotic dynamics.     

A complementarity-based rolling friction model for rigid contacts

In this work (also, preprint ANL/MCS-P3020-0812, Argonne National Laboratory) we introduce a complementarity-based rolling friction model to characterize dissipative phenomena at the interface between moving parts. Since the formulation is based on differential inclusions, the model fits well in the context of nonsmooth dynamics, and it does not require short integration timesteps. The method encompasses a rolling resistance limit for static cases, similar to what happens for sliding friction; this is a simple yet efficient approach to problems involving transitions from rolling to resting, and vice-versa. We propose a convex relaxation of the formulation in order to achieve algorithmic robustness and stability; moreover, we show the side effects of the convexification. A natural application of the model is the dynamics of granular materials, because of the high computational efficiency and the need for only a small set of parameters. In particular, when used as a micromechanical model for rolling resistance between granular particles, the model can provide an alternative way to capture the effect of irregular shapes. Other applications can be related to real-time simulations of rolling parts in bearing and guideways, as shown in examples.     

Level‐set based numerical simulation of a migrating and dissolving liquid drop in a cylindrical cavity

In the present paper the dissolution of a binary liquid drop having a miscibility gap and migrating due to thermo‐solutal capillary convection in a cylindrical cavity is studied numerically. The interest in studying this problem is twofold. From a side, in the absence of gravity, capillary migration is one of the main physical mechanisms to set into motion dispersed liquid phases and from the other side, phase equilibria of multi‐component liquid systems, ubiquitous in applications, often exhibit a miscibility gap. The drop capillary migration is due to an imposed temperature gradient between the cavity top and bottom walls. The drop dissolution is due to the fact that initial composition and volume values, and thermal boundary conditions are only compatible with a final single phase equilibrium state. In order to study the drop migration along the cavity and the coupling with dissolution, a previously developed planar two‐dimensional code is extended to treat axis‐symmetric geometries. The code is based on a finite volume formulation. A level‐set technique is used for describing the dynamics of the interface separating the different phases and for mollifying the interface discontinuities between them. The level‐set related tools of redistancing and off‐interface extension are used to enhance code resolution in the critical interface region. Migration speeds and volume variations are determined for different drop radii. Copyright © 2004 John Wiley & Sons, Ltd.     

Multiscale coupling of molecular dynamics and peridynamics

We propose a multiscale computational model to couple molecular dynamics and peridynamics. The multiscale coupling model is based on a previously developed multiscale micromorphic molecular dynamics (MMMD) theory, which has three dynamics equations at three different scales, namely, microscale, mesoscale, and macroscale. In the proposed multiscale coupling approach, we divide the simulation domain into atomistic region and macroscale region. Molecular dynamics is used to simulate atom motions in atomistic region, and peridynamics is used to simulate macroscale material point motions in macroscale region, and both methods are nonlocal particle methods. A transition zone is introduced as a messenger to pass the information between the two regions or scales. We employ the “supercell” developed in the MMMD theory as the transition element, which is named as the adaptive multiscale element due to its ability of passing information from different scales, because the adaptive multiscale element can realize both top-down and bottom-up communications. We introduce the Cauchy–Born rule based stress evaluation into state-based peridynamics formulation to formulate atomistic-enriched constitutive relations. To mitigate the issue of wave reflection on the interface, a filter is constructed by switching on and off the MMMD dynamic equations at different scales. Benchmark tests of one-dimensional (1-D) and two-dimensional (2-D) wave propagations from atomistic region to macro region are presented. The mechanical wave can transit through the interface smoothly without spurious wave deflections, and the filtering process is proven to be efficient.     

The Interface Between Fresh and Salt Groundwater in Horizontal Aquifers: The Dupuit–Forchheimer Approximation Revisited

We analyze the motion of a sharp interface between fresh and salt groundwater in horizontal, confined aquifers of infinite extend. The analysis is based on earlier results of De Josselin de Jong (Proc Euromech 143:75–82, 1981). Parameterizing the height of the interface along the horizontal base of the aquifer and assuming the validity of the Dupuit–Forchheimer approximation in both the fresh and saltwater, he derived an approximate interface motion equation. This equation is a nonlinear doubly degenerate diffusion equation in terms of the height of the interface. In that paper, he also developed a stream function-based formulation for the dynamics of a two-fluid interface. By replacing the two fluids by one hypothetical fluid, with a distribution of vortices along the interface, the exact discharge field throughout the flow domain can be determined. Starting point for our analysis is the stream function formulation. We derive an exact integro-differential equation for the movement of the interface. We show that the pointwise differential terms are identical to the approximate Dupuit–Forchheimer interface motion equation as derived by De Josselin de Jong. We analyze (mathematical) properties of the additional integral term in the exact interface motion formulation to validate the approximate Dupuit–Forchheimer interface motion equation. We also consider the case of flat interfaces, and we study the behavior of the toe of the interface. In particular, we give a criterion for finite or infinite speed of propagation.     

Spatially resolved observations of strain fields at necking and fracture of anisotropic hardened steel sheet material

In this work plastic strain localization, also referred to as necking, of press-hardened ultra-high strength steel is observed using digital speckle correlation. The region of the neck is studied during tensile tests of specimens specially designed to facilitate strain localization at an inner point of the material, thus avoiding edge effects on localization and fracture. By using measurements with a length scale small enough to properly resolve the neck, its growth and shape can be studied. Furthermore, the anisotropy of the material is investigated by examining specimens cut out at different angles to the rolling direction. It is seen that the local fracture strain of specimens cut out along the rolling direction is approximately twice as high as it is for specimens cut out perpendicular to the rolling direction.     

Motion of a solitary wave over a submerged ridge

The two- and three-dimensional motion of a stationary wave in a layer of liquid with a variable depth was studied earlier by Grimshaw [1, 2], who showed that if the unperturbed state of the liquid were a state of rest, then a certain integrated quantity characterizing the energy of the wave (referred to one element on the leading edge of the wave) was conserved during the motion. In the present investigation (which is based on this property) we shall establish the shape of the wave front for the case of steady motion over an infinite submerged cylindrical ridge; we shall present a model describing a wave of limiting amplitude and shall consider the transient perturbation of the steady-state solution. We should note that one particular case (the motion of a straight front, orthogonal to the axis of the ridge, in a direction parallel to this axis) was considered by another method in [3].     

Simulations of fibre orientation in dilute suspensions with front moving in the filling process of a rectangular channel using level-set method

The simulation of fibre orientation in dilute suspension with front moving is carried out using the projection and level-set methods. The motion of fibres is described using the Jeffery equation, and the contribution of fibres to the flow is accounted for by the configuration-field method. The dilute suspension of short fibres in Newtonian fluids is considered. The governing Navier–Stokes equation for the fluid flow is solved using the projection method with finite difference scheme, while the fibre-related equations are directly solved with the Runge–Kutta method. In the present study for fibres in dilute suspension flow for injection molding, the effects of various flow and material parameters on the fibre orientation, the velocity distributions and the shapes of the leading flow front are found and discussed. Our findings indicate that the presence of fibre motion has little influence on the front shape in the ranges of fibre parameters studied at the fixed Reynolds number. Influence of changing fibre parameters only causes variation of front shape in the region near the wall, and the front shape in the central core area does not vary much with the fibre parameters. On the other hand, the fibre motion has strong influence on the distributions of the streamwise and transverse velocities in the fountain flow. Fibre motion produces strong normal stress near the wall which leads to the reduction of transversal velocity as compared to the Newtonian flow without fibres, which in turn, leads to the increased streamwise velocity near the wall. Thus, the fibre addition to the flow weakens the strength of the fountain flow. The Reynolds number has also displayed significant influence on the distribution of the streamwise velocity behind the flow front for a given fibre concentration. It is also found that the fibre orientation is not always along the direction of the velocity vector in the process of mold filling. In the region of the fountain flow, the fibre near the centreline is more oriented across the streamwise direction compared to that in the region far behind the flow front. This leads to the fact that the fibre near the centreline in the region of fountain flow is more extended along the transverse direction. As the fibre orientation in the suspension flow and the shape of the flow front have great bearing on the quality of the product made from injection molding, this study has much implications for engineering applications. These results can also be useful in other fields dealing with fibre suspensions.     

Points singuliers d'une ligne de contact mobile

It is proposed to represent the dynamics of a moving contact line by an Onsager like mobility relation between the contact angle and the speed of the moving line, including an Arrhenius factor small enough in many physical situations to be the limiting factor for the motion. The liquid-vapor interface is then in quasiequilibrium, which allows one to analyse a dynamical wetting transition. This approach predicts well the formation of angular points on the rear edge of droplets sliding on a tilted plane.