Source author record

Steffen Hardt

Steffen Hardt appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

25works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

25 published item(s)

preprint2022arXiv

Controlling the electrostatic Faraday instability using superposed electric fields

When the interface between a dielectric and a conducting liquid is excited by an oscillatory electric field, electrostatic Faraday waves can be induced. Here, we study the response of the interface to an AC electric field, which is superposed by either a second AC field of different frequency, or by a DC field. An algorithm based on light refraction at the fluidic interface is used to obtain the spatio-temporal response of the Faraday waves, and the critical voltage corresponding to onset of instability, the interfacial oscillation frequency and the dominant wavelengths are determined. The influence of the mixing ratio, which denotes the relative amplitudes of the different components of the driving signal, is analyzed, and the experimental results are compared with theoretical predictions. For AC/AC driving, gradual variations of the mixing ratio can induce a jump of the pattern wavelength, which is a result of the transition from harmonic to a subharmonic oscillation. For AC/DC driving, the interface oscillates either harmonically or subharmonically, and the response wavelength can be tuned continuously by adjusting the admixture of the DC component. For both driving modes, the experiments show good agreement with theory.

preprint2022arXiv

Deformation modes of an oil-water interface under a local electric field: From Taylor cones to surface dimples

Fluidic interfaces disintegrate under sufficiently strong electric fields, leading to electrohydrodynamic (EHD) tip streaming. Taylor cones, which emit charged droplets from the tip of a conical cusp, are among the most prominent and well-studied examples of EHD instabilities. In liquid-liquid systems, more complex interface deformation modes than simple Taylor cones can be observed, with the interface being pushed away from the electrode, and additional cone structures emerging from the rim of the dimple. In this article, we investigate the mechanisms behind these deformation modes experimentally and numerically, and demonstrate that the presence of droplets triggers the dimple at the interface. In order to characterize the underlying processes, we replace the pin electrode by a hollow metallic needle with a prescribed electrolyte volume flow. The submerged electrospray introduces droplets of an aqueous KCl solution with varying ion concentrations into silicone oils with varying viscosities. By measuring the corresponding electric current and by optical investigation of the interface deformation, we study the system response to variations of the ionic concentration, viscosity, applied voltage as well as flow rate. In addition to the experiments, we use a finite element solver and compute the charge transport due to the droplets in the oil phase. Further, we compute the electric potential distribution, flow field and interface deformation. After calibration of our model with particle tracking velocimetry data of the flow inside the oil phase, we reproduce the experimentally observed dimple at the liquid-liquid interface. In summary, this work highlights the importance of charged droplets for the complex dynamic modes observed when a liquid-liquid interface is exposed to a local electric field.

preprint2022arXiv

Resonantly-driven nanopores can serve as nanopumps

Inducing transport in electrolyte-filled nanopores with dc fields has led to influential applications ranging from nanosensors to DNA sequencing. Here we use the Poisson-Nernst-Planck and Navier-Stokes equations to show that unbiased ac fields can induce comparable directional flows in gated conical nanopores. This flow exclusively occurs at intermediate driving frequencies and hinges on the resonance of two competing timescales, representing space charge development at the ends and in the interior of the pore. We summarize the physics of resonant nanopumping in an analytical model that reproduces the results of numerical simulations. Our findings provide a generic route towards real-time controllable flow patterns, which might find applications in controlling the translocation of particles such as small molecules or nanocolloids.

preprint2020arXiv

Flow and streaming potential of an electrolyte in a channel with an axial temperature gradient

The effect of an axial temperature gradient on the flow profile and the induced streaming potential of a pressure-driven symmetric electrolyte in a slit channel is investigated. Based on the non-isothermal Nernst-Planck equations as well as the Poisson equation in the lubrication approximation, expressions for the ion distribution in the electric double layer (EDL) are derived. It is found that thermophoretic ion motion and a temperature-dependent electrophoretic ion mobility increase the local EDL thickness with temperature, whereas a temperature-dependent permittivity shrinks the EDL. Within the Debye-Hückel approximation, the Navier-Stokes equation with the corresponding electric body force terms is solved. Analytical expressions for the flow profile and the induced (streaming) field under non-isothermal conditions are derived. It is shown that for such a situation the induced electric field is the linear superposition of at least seven individual contributions. For very wide channels, only the thermoelectric field typically present in bulk electrolytes when subjected to a temperature gradient (Soret equilibrium) as well as the conventional pressure-induced streaming field are of importance. Under extreme confinement, selective thermo-electro-migration driven by the interplay between the temperature-dependent electrophoretic ion mobility and the interaction of the ions with the surface wall charge causes a thermoelectric field of non-advective origin. For wider channels and besides the well-known thermoosmosis due to the temperature dependence of the dielectric permittivity, it is demonstrated that a temperature gradient renders the ion cloud in the EDL out of mechanical equilibrium. This leads to a thermoosmotic flow, and the ion advection affiliated with it may induce a thermoelectric field of similar order of magnitude as the one caused by more conventional thermal effects.

preprint2016arXiv

Conjugated liquid layers driven by the short-wavelength Bénard-Marangoni instability: experiment and numerical simulation

The coupled dynamics of two conjugated liquid layers of disparate thicknesses, which coat a solid substrate and are subjected to a transverse temperature gradient, is investigated. The upper liquid layer evolves under the short-wavelength (SW) Bénard-Marangoni (BM) instability, whereas the lower, much thinner film undergoes a shear-driven long-wavelength deformation. Although the lubricating film should reduce the viscous stresses acting on the up to one hundred times thicker upper layer by only 10%, it is found that the critical Marangoni number of marginal stability may be as low as if a stress-free boundary condition were applied at the bottom of the upper layer, i.e. much lower than the classical value of 79.6 known for a single film. Furthermore, it is experimentally verified that the deformation of the liquid-liquid interface, albeit small, has a non-negligible effect on the temperature distribution along the liquid-gas interface of the upper layer. This stabilizes the hexagonal pattern symmetry towards external disturbances and indicates a two-way coupling of the different layers. The experiments also demonstrate how convection patterns formed in a liquid film can be used to pattern a second conjugated film. The experimental findings are verified by a numerical model of the coupled layers.

preprint2016arXiv

Interfacial instability of liquid films coating the walls of a parallel-plate channel and sheared by a gas flow

The stability of liquid films coating the walls of a parallel-plate channel and sheared by a pressure-driven gas flow along the channel centre is studied. The films are susceptible to a long-wavelength instability, whose dynamic behaviour is found - for sufficiently low Reynolds numbers and thick gas layers - to be described by two coupled non-linear partial differential equations. To the best of our knowledge, such coupled fully non-linear equations for the film thicknesses have not been derived previously. A linear stability analysis conducted under the condition that the material properties and the initial undisturbed liquid film thicknesses are equal can be utilized to determine whether the interfaces are predominantly destabilized by the variations of the shear stress or by the pressure gradient acting upon them. The analysis of the weakly non-linear equations performed for this case shows that instabilities corresponding to a vanishing Reynolds number are absent from the system. Moreover, for this configuration, the patterns emerging along the two interfaces are found to be identical in the long-time limit, implying that the films are fully synchronized. A different setup, where the liquid films have identical material properties but their undisturbed thicknesses differ, is studied numerically. The results show that even for this configuration the interfacial waves remain phase-synchronized and closely correlated for a broad range of parameters. These findings are particularly relevant for multiphase flow in narrow ducts, for example in the respiratory system or in microfluidic channels.

preprint2015arXiv

Brownian dynamics of rigid particles in an incompressible fluctuating fluid by a meshfree method

A meshfree Lagrangian method for the fluctuating hydrodynamic equations (FHEs) with fluid-structure interactions is presented. Brownian motion of the particle is investigated by direct numerical simulation of the fluctuating hydrodynamic equations. In this framework a bidirectional coupling has been introduced between the fluctuating fluid and the solid object. The force governing the motion of the solid object is solely due to the surrounding fluid particles. Since a meshfree formulation is used, the method can be extended to many real applications involving complex fluid flows. A three-dimensional implementation is presented. In particular, we observe the short and long-time behaviour of the velocity autocorrelation function (VACF) of Brownian particles and compare it with the analytical expression. Moreover, the Stokes-Einstein relation is reproduced to ensure the correct long-time behaviour of Brownian dynamics.

preprint2015arXiv

Drag and diffusion coefficients of a spherical particle attached to a fluid interface

Explicit analytical expressions for the drag and diffusion coefficients of a spherical particle attached to the interface between two immiscible fluids are constructed for the case of a small viscosity ratio between the fluid phases. The model is designed to explicitly account for the dependence on the contact angle between the two fluids and the solid surface. The Lorentz reciprocal theorem is applied in the context of a geometric perturbation approach, which is based on the deviation of the contact angle from a 90°-value. By testing the model against experimental and numerical data from the literature, good agreement is found within the entire range of contact angles below 90°. As an advantage of the method reported, the drag and diffusion coefficients can be calculated up to second order in the perturbation parameter, while it is sufficient to know the velocity and pressure fields only up to first order. Extensions to other particle shapes with known velocity and pressure fields are straightforward.

preprint2015arXiv

Driven particles at fluid interfaces acting as capillary dipoles

The dynamics of spherical particles driven along an interface between two immiscible fluids is investigated asymptotically. Under the assumptions of a pinned three-phase contact line and very different viscosities of the two fluids, a particle assumes a tilted orientation. As it moves, it causes a deformation of the fluid interface which is also computed. The case of two interacting driven particles is studied via the Linear Superposition Approximation. It is shown that the capillary interaction force resulting from the particle motion is dipolar in terms of the azimuthal angle and decays with the fifth power of the inter-particle separation, similar to a capillary quadrupole originating from undulations of the three-phase contact line. The dipolar interaction is demonstrated to exceed the quadrupolar interaction at moderate particle velocities.

preprint2014arXiv

Coupled self-organization: Thermal interaction between two liquid films undergoing long-wavelength instabilities

The effects of thermal coupling between two thin liquid layers, separated by a gas layer, are discussed. The liquid layers undergo long-wavelength instabilities driven by gravitational and thermocapillary stresses. To study the dynamics, both a linear stability analysis and a full numerical solution of the thin-film equations are performed. The results demonstrate that the stability properties of the combined system differ substantially from the case where both layers evolve independently from each other. Most prominently, oscillatory instabilities, not present in single-liquid layer configurations, may occur.

preprint2014arXiv

Droplet Formation via Solvent Shifting in a Microfluidic Device

Solvent shifting is a process in which a non-solvent is added to a solvent/solute mixture and extracts the solvent. The solvent and the non-solvent are miscible. Because of solution supersaturation a portion of the solute transforms to droplets. In this paper, based on this process, we present an investigation on droplet formation and their radial motion in a microfluidic device in which a jet is injected in a co-flowing liquid stream. Thanks to the laminar flow, the microfluidic setup enables studying diffusion mass transfer in radial direction and obtaining well-defined concentration distributions. Such profiles together with Ternary Phase Diagram (TPD) give detailed information about the conditions for droplet formation condition as well as their radial migration in the channel. The ternary system is composed of ethanol (solvent), de-ionized water (non-solvent) and divinyle benzene (solute). We employ analytical/numerical solutions of the diffusion equation to obtain concentration profiles of the components. We show that in the system under study droplets are formed in a region of the phase diagram between the binodal and the spinodal, i.e. via a thermally activated process. The droplets are driven to the channel centerline by the solutal Marangoni effect but are not able to significantly penetrate into the single-phase region, where they get rapidly dissolved. Therefore, the radial motion of the binodal surface carries the droplets to the centerline where they get collected.

preprint2014arXiv

Energy harvesting through gas dynamics in the free molecular flow regime between structured surfaces at different temperatures

For a gas confined between surfaces held at different temperatures the velocity distribution shows a significant deviation from the Maxwell distribution when the mean free path of the molecules is comparable to or larger than the channel dimensions. If one of the surfaces is suitably structured, this non-equilibrium distribution can be exploited for momentum transfer in tangential direction between the two surfaces. This opens up the possibility to extract work from the system which operates as a heat engine. Since both surfaces are held at constant temperatures, the mode of momentum transfer is different from thermal creep flow that has gained more attention so far. This situation is studied in the limit of free-molecular flow for the case that an unstructured surface is allowed to move tangentially with respect to a structured surface. Parameter studies are conducted, and configurations with maximum thermodynamic efficiency are identified. Overall, it is shown that significant efficiencies can be obtained by tangential momentum transfer between structured surfaces.

preprint2014arXiv

Line tension and reduction of apparent contact angle associated with electric double layers

The line tension of an electrolyte wetting a non-polar substrate is computed analytically and numerically. The results show that, depending on the value of the apparent contact angle, positive or negative line tension values may be obtained. Furthermore, a significant difference between Young's contact angle and the apparent contact angle measured several Debye lengths remote from the three-phase contact line occurs. When applying the results to water wetting highly charged surfaces, line tension values of the same order of magnitude as found in recent experiments can be achieved. Therefore, the theory presented may contribute to the understanding of line tension measurements and points to the importance of the electrostatic line tension. Being strongly dependent on the interfacial charge density, electrostatic line tension is found to be tunable via the pH value of the involved electrolyte. As a practical consequence, the stability of nanoparticles adsorbed at fluid-fluid interfaces is predicted to be dependent on the pH value. The theory is suited for future incorporation of effects due to surfactants where even larger line tension values can be expected.

preprint2014arXiv

Numerical simulation of moving rigid body in rarefied gases

In this paper we present a numerical scheme to simulate a moving rigid body with arbitrary shape suspended in a rarefied gas. The rarefied gas is simulated by solving the Boltzmann equation using a DSMC particle method. The motion of the rigid body is governed by the Newton-Euler equations, where the force and the torque on the rigid body is computed from the momentum transfer of the gas molecules colliding with the body. On the other hand, the motion of the rigid body influences the gas flow in its surroundings. We validate the numerical results by testing the Einstein relation for Brownian motion of the suspended particle. The translational as well as the rotational degrees of freedom are taken into account. It is shown that the numerically computed translational and rotational diffusion coefficients converge to the theoretical values.

preprint2014arXiv

Numerical simulation of wetting phenomena by a meshfree particle method

Simulations of wetting phenomena by a meshfree particle method are presented. The incompressible Navier-Stokes equations are used to model the two-phase flow. The continuous surface force model is used to incorporate the surface tension force. Chorin's projection method is applied to discretize the Navier-Stokes equations. The different fluid phases are identified by assigning different colors and different material properties (density, viscosity) to the particles that remain unchanged throughout a simulation. Two-phase flow is captured by a one-fluid model via using weighted averages of the density and viscosity in a region around the fluid-fluid interface. The differential operators at each particle are computed from the surrounding cloud of particles with the help of the least-squares method. The numerical results are compared with specific analytical solutions, but also with previously considered test cases involving wetting of a container and sessile drops. A good overall agreement is found.

preprint2013arXiv

Influence of the enclosed fluid on the flow over a microstructured surface in the Cassie state

Analytical expressions for the flow field as well as for the effective slip length of a shear flow over a surface with periodic rectangular grooves are derived. The primary fluid is in the Cassie state with the grooves being filled with a secondary immiscible fluid. The coupling of both fluids is reflected in a locally varying slip distribution along the fluid-fluid interface, which models the effect of the secondary fluid on the outer flow. The obtained closed-form analytical expressions for the flow field and effective slip length of the primary fluid explicitly contain the influence of the viscosities of the two fluids as well as the magnitude of the local slip, which is a function of the surface geometry. They agree well with results from numerical computations of the full geometry. The analytical expressions allow investigating the influence of the viscous stresses inside the secondary fluid for arbitrary geometries of the rectangular grooves. For classic superhydrophobic surfaces, the deviations in the effective slip length compared to the case of inviscid gas flow are are pointed out. Another important finding with respect to an accurate modeling of flow over microstructured surfaces is that the local slip length of a grooved surface is anisotropic.

preprint2013arXiv

Propulsion Mechanisms for Leidenfrost Solids on Ratchets

We propose a model for the propulsion of Leidenfrost solids on ratchets based on viscous drag due to the flow of evaporating vapor. The model assumes pressure-driven flow described by the Navier-Stokes equations and is mainly studied in lubrication approximation. A scaling expression is derived for the dependence of the propulsive force on geometric parameters of the ratchet surface and properties of the sublimating solid. We show that the model results as well as the scaling law compare favorably with experiments and are able to reproduce the experimentally observed scaling with the size of the solid.

preprint2013arXiv

Sample dispersion in isotachophoresis with Poiseuille counterflow

A particular mode of isotachophoresis (ITP) employs a pressure-driven flow opposite to the sample electromigration direction in order to anchor a sample zone at a specific position along a channel or capillary. We investigate this situation using a two-dimensional finite-volume model based on the Nernst-Planck equation. The imposed Poiseuille flow profile leads to a significant dispersion of the sample zone. This effect is detrimental for the resolution in analytical applications of ITP. We investigate the impact of convective dispersion, characterized by the area-averaged width of a sample zone, for various values of the sample Péclet-number, as well as the relative mobilities of the sample and the adjacent electrolytes. A one-dimensional model for the area-averaged concentrations based on a Taylor-Aris-type effective axial diffusivity is shown to yield good agreement with the finite-volume calculations. This justifies the use of such simple models and opens the door for the rapid simulation of ITP protocols with Poiseuille counterflow.

preprint2012arXiv

Electric double layer structure close to the three-phase contact line in an electrolyte wetting a solid substrate

The electric double layer structure in an electrolyte close to a solid substrate near the three-phase contact line is approximated by considering the linearized Poisson-Boltzmann equation in a wedge geometry. The mathematical approach complements the semi-analytical solutions reported in the literature by providing easily available characteristic information on the double layer structure. In particular, the model contains a length scale that quantifies the distance from the fluid-fluid interface over which this boundary influences the electric double layer. The analysis is based on an approximation for the equipotential lines. Excellent agreement between the model predictions and numerical results is achieved for a significant range of contact angles. The length scale quantifying the influence of the fluid-fluid interface is proportional to the Debye length and depends on the wall contact angle. It is shown that for contact angles approaching 90° there is a finite range of boundary influence.

preprint2012arXiv

Longitudinal and transversal flow over a cavity containing a second immiscible fluid

An analytical solution for the flow field of a shear flow over a rectangular cavity containing a second immiscible fluid is derived. While flow of a single-phase fluid over a cavity is a standard case investigated in fluid dynamics, flow over a cavity which is filled with a second immiscible fluid, has received little attention. The flow filed inside the cavity is considered to define a boundary condition for the outer flow which takes the form of a Navier slip condition with locally varying slip length. The slip-length function is determined from the related problem of lid-driven cavity flow. Based on the Stokes equations and complex analysis it is then possible to derive a closed analytical expression for the flow field over the cavity for both the transversal and the longitudinal case. The result is a comparatively simple function, which displays the dependence of the flow field on the cavity geometry and the medium filling the cavity. The analytically computed flow field agrees well with results obtained from a numerical solution of the Navier-Stokes equations. The studies presented in this article are of considerable practical relevance, for example for the flow over superhydrophobic surfaces.

preprint2012arXiv

Simulation of a moving liquid droplet inside a rarefied gas region

We study the dynamics of a liquid droplet inside a gas over a large range of the Knudsen numbers. The moving liquid droplet is modeled by the incompressible Navier-Stokes equations, the surrounding rarefied gas by the Boltzmann equation. The interface boundary conditions between the gas and liquid phases are derived. The incompressible Navier-Stokes equations are solved by a meshfree Lagrangian particle method called Finite Pointset Method (FPM), and the Boltzmann equation by a DSMC type of particle method. To validiate the coupled solutions of the Boltzmann and the incompressible Navier-Stokes equations we have further solved the compressible and the incompressible Navier-Stokes equations in the gas and liquid phases, respectively. In the latter case both the compressible and the incompressible Navier-Stokes equations are also solved by the FPM. In the continuum regime the coupled solutions obtained from the Boltzmann and the incompressible Navier-Stokes equations match with the solutions obtained from the compressible and the incompressible Navier-Stokes equations. In this paper, we presented solutions in one-dimensional physical space.

preprint2012arXiv

Thermally-driven flows between a Leidenfrost solid and a ratchet surface

The significance of thermally-driven flows for the propulsion of Leidenfrost solids on a ratchet surface is studied based on a numerical solution of the Boltzmann equation. In contrast to a previous analysis, it is found that no significant thermal creep flow is established. Instead, the flow pattern is dominated by thermal edge and thermal-stress slip flow, the latter being directed opposite to thermal creep flow. However, in total thermally-induced flows only make a minor contribution to the propulsion of Leidenfrost solids on ratchet surfaces which is dominated by the pressure-driven flow due to the sublimating solid.

preprint2011arXiv

Momentum and Mass Fluxes in a Gas Confined between Periodically Structured Surfaces at Different Temperatures

It is well known that in a gas-filled duct or channel along which a temperature gradient is applied, a thermal creep flow is created. Here we show that a mass and momentum flux can also be induced in a gas confined between two parallel structured surfaces at different temperatures, i.e. \textit{orthogonal} to the temperature gradient. We use both analytical and numerical methods to compute the resulting fluxes. The momentum flux assumes its maximum value in the free-molecular flow regime, the (normalized) mass flux in the transition flow regime. The discovered phenomena could find applications in novel methods for energy-conversion and thermal pumping of gases.

preprint2011arXiv

Non-equilibrium electric double layers at the interface between two electrolytes

The charge distribution at the interface between two electrolytes is studied for the case of non-vanishing ion fluxes. The analysis is an extension of the established Verwey-Niessen theory to non-equilibrium situations. Applying matched asymptotic expansions for the region of the electric double layer and the bulk electrolytes, analytical expressions for the ion concentrations and the electrostatic potential are derived. It is found that the electric double layer is surprisingly stable. Even in cases where the applied electric field is so strong that it completely cancels the field of the charge clouds at the interface, two opposing space charge regions are formed. This phenomenon is qualitatively explained using a simple model incorporating ion partitioning and diffusive ion transport. Another notable phenomenon is that the interface acquires a net charge in the presence of an applied electric field, being of relevance for the hydrodynamic stability of interfacial flows.

preprint2010arXiv

Thermocapillary Flow on Superhydrophobic Surfaces

A liquid in Cassie-Baxter state above a structured superhydrophobic surface is ideally suited for surface driven transport due to its large free surface fraction in close contact to a solid. We investigate thermal Marangoni flow over a superhydrophobic array of fins oriented parallel or perpendicular to an applied temperature gradient. In the Stokes limit we derive an analytical expression for the bulk flow velocity above the surface and compare it with numerical solutions of the Navier-Stokes equation. Even for moderate temperature gradients comparatively large flow velocities are induced, suggesting to utilize this principle for microfluidic pumping.