Formation of Nanopillar Arrays in Ultrathin Viscous Films: The Critical Role of Thermocapillary Stresses

Experiments by several groups during the past decade have shown that a molten polymer nanofilm subject to a large transverse thermal gradient undergoes spontaneous formation of periodic nanopillar arrays. The prevailing explanation is that coherent reflections of acoustic phonons within the film cause a periodic modulation of the radiation pressure which enhances pillar growth. By exploring a deformational instability of particular relevance to nanofilms, we demonstrate that thermocapillary forces play a crucial role in the formation process. Analytic and numerical predictions show good agreement with the pillar spacings obtained in experiment. Simulations of the interface equation further determine the rate of pillar growth of importance to technological applications.Comment: 5 pages, 4 figure

Model of the meniscus of an ionic liquid ion source.

A simple model of the transfer of charge and ion evaporation in the meniscus of an ionic-liquid ion source working in the purely ionic regime is proposed on the basis of order-of-magnitude estimates which show that, in this regime, _i_ the flow in the meniscus is dominated by the viscosity of the liquid and is affected very little by the mass flux accompanying ion evaporation, and _ii_ the effect of the space charge around the evaporating surface is negligible and the evaporation current is controlled by the finite electrical conductivity of the liquid. The model predicts that a stationary meniscus of a very polar liquid undergoing ion evaporation is nearly hydrostatic and can exist only below a certain value of the applied electric field, at which the meniscus attains its maximum elongation but stays smooth. The electric current vs applied electric field characteristic displays a frozen regime of negligible ion evaporation at low fields and a conduction-controlled regime at higher fields, with a sharp transition between the two regimes owing to the high sensitivity of the ion evaporation rate to the electric field. A simplified treatment of the flow in the capillary or liquid layer through which liquid is delivered to the meniscus shows that the size of the meniscus decreases and the maximum attainable current increases when the feeding pressure is decreased, and that appropriate combinations of feeding pressure and pressure drop may lead to high maximum currents

Transport coefficients for electrolytes in arbitrarily shaped nano and micro-fluidic channels

We consider laminar flow of incompressible electrolytes in long, straight channels driven by pressure and electro-osmosis. We use a Hilbert space eigenfunction expansion to address the general problem of an arbitrary cross section and obtain general results in linear-response theory for the hydraulic and electrical transport coefficients which satisfy Onsager relations. In the limit of non-overlapping Debye layers the transport coefficients are simply expressed in terms of parameters of the electrolyte as well as the geometrical correction factor for the Hagen-Poiseuille part of the problem. In particular, we consider the limits of thin non-overlapping as well as strongly overlapping Debye layers, respectively, and calculate the corrections to the hydraulic resistance due to electro-hydrodynamic interactions.Comment: 13 pages including 4 figures and 1 table. Typos corrected. Accepted for NJ

Entrance effects in concentration-gradient-driven flow through an ultrathin porous membrane

Published Online: 29 July 2019Transport of liquid mixtures through porous membranes is central to processes such as desalination, chemical separations, and energy harvesting, with ultrathin membranes made from novel 2D nanomaterials showing exceptional promise. Here, we derive, for the first time, general equations for the solution and solute fluxes through a circular pore in an ultrathin planar membrane induced by a solute concentration gradient. We show that the equations accurately capture the fluid fluxes measured in finite-element numerical simulations for weak solute-membrane interactions. We also derive scaling laws for these fluxes as a function of the pore size and the strength and range of solute-membrane interactions. These scaling relationships differ markedly from those for concentration-gradient-driven flow through a long cylindrical pore or for flow induced by a pressure gradient or an electric field through a pore in an ultrathin membrane. These results have broad implications for transport of liquid mixtures through membranes with thickness on the order of the characteristic pore size.Daniel J. Rankin, Lydéric Bocquet, and David M. Huan

Electrohydrodynamics within electrical double layer in a pressure-driven flow in presence of finite temperature gradients

A wide spectrum of electrokinetic studies is modelled as isothermal ones to expedite analysis even when such conditions may be extremely difficult to realize in practice. As a clear and novel departure from this trend, we address the case of flow-induced electrohydrodynamics, commonly referred to as streaming potential, in a situation where finite temperature gradients do indeed exist. By way of analysing a model problem of flow through a narrow parallel plate channel, we show that the temperature gradients have a significant effect on the streaming potential, and, consequently, on the flow itself. We incorporate thermoelectric effects in our model by a full-fledged coupling among the electric potential, the ionic species distribution, the fluid velocity and the local fluid temperature fields without resorting to ad hoc simplifications. We expect this expository study to contribute towards more sophisticated future inquiries into practical micro-/nano-fluidic applications coupling thermal field focusing with electrokinetic effects.Comment: 13 pages, 5 figure

Aerothermodynamic Analysis of a Reentry Brazilian Satellite

This work deals with a computational investigation on the small ballistic reentry Brazilian vehicle SARA (acronyms for SAt\'elite de Reentrada Atmosf\'erica). Hypersonic flows over the vehicle SARA at zero-degree angle of attack in a chemical equilibrium and thermal non-equilibrium are modeled by the Direct Simulation Monte Carlo (DSMC) method, which has become the main technique for studying complex multidimensional rarefied flows, and that properly accounts for the non-equilibrium aspects of the flows. The emphasis of this paper is to examine the behavior of the primary properties during the high altitude portion of SARA reentry. In this way, velocity, density, pressure and temperature field are investigated for altitudes of 100, 95, 90, 85 and 80 km. In addition, comparisons based on geometry are made between axisymmetric and planar two-dimensional configurations. Some significant differences between these configurations were noted on the flowfield structure in the reentry trajectory. The analysis showed that the flow disturbances have different influence on velocity, density, pressure and temperature along the stagnation streamline ahead of the capsule nose. It was found that the stagnation region is a thermally stressed zone. It was also found that the stagnation region is a zone of strong compression, high wall pressure. Wall pressure distributions are compared with those of available experimental data and good agreement is found along the spherical nose for the altitude range investigated.Comment: The paper will be published in Vol. 42 of the Brazilian Journal of Physic

Simulating (electro)hydrodynamic effects in colloidal dispersions: smoothed profile method

Previously, we have proposed a direct simulation scheme for colloidal dispersions in a Newtonian solvent [Phys.Rev.E 71,036707 (2005)]. An improved formulation called the ``Smoothed Profile (SP) method'' is presented here in which simultaneous time-marching is used for the host fluid and colloids. The SP method is a direct numerical simulation of particulate flows and provides a coupling scheme between the continuum fluid dynamics and rigid-body dynamics through utilization of a smoothed profile for the colloidal particles. Moreover, the improved formulation includes an extension to incorporate multi-component fluids, allowing systems such as charged colloids in electrolyte solutions to be studied. The dynamics of the colloidal dispersions are solved with the same computational cost as required for solving non-particulate flows. Numerical results which assess the hydrodynamic interactions of colloidal dispersions are presented to validate the SP method. The SP method is not restricted to particular constitutive models of the host fluids and can hence be applied to colloidal dispersions in complex fluids

THE SPIRAL WAVE INSTABILITY INDUCED BY A GIANT PLANET. I. PARTICLE STIRRING IN THE INNER REGIONS OF PROTOPLANETARY DISKS

We have recently shown that spiral density waves propagating in accretion disks can undergo a parametric instability by resonantly coupling with and transferring energy into pairs of inertial waves (or inertial-gravity waves when buoyancy is important). In this paper, we perform inviscid three-dimensional global hydrodynamic simulations to examine the growth and consequence of this instability operating on the spiral waves driven by a Jupiter-mass planet in a protoplanetary disk. We find that the spiral waves are destabilized via the spiral wave instability (SWI), generating hydrodynamic turbulence and sustained radially-alternating vertical flows that appear to be associated with long wavelength inertial modes. In the interval $0.3~R_p \leq R \leq 0.7~R_p$, where $R_p$ denotes the semi-major axis of the planetary orbit (assumed to be 5~au), the estimated vertical diffusion rate associated with the turbulence is characterized by $\alpha_{\rm diff} \sim (0.2-1.2) \times 10^{-2}$. For the disk model considered here, the diffusion rate is such that particles with sizes up to several centimeters are vertically mixed within the first pressure scale height. This suggests that the instability of spiral waves launched by a giant planet can significantly disperse solid particles and trace chemical species from the midplane. In planet formation models where the continuous local production of chondrules/pebbles occurs over Myr time scales to provide a feedstock for pebble accretion onto these bodies, this stirring of solid particles may add a time constraint: planetary embryos and large asteroids have to form before a gas giant forms in the outer disk, otherwise the SWI will significantly decrease the chondrule/pebble accretion efficiency.Comment: Accepted for publication in the The Astrophysical Journal, 19 pages, 12 figures, 1 tabl

Diffuse charge and Faradaic reactions in porous electrodes

Porous electrodes instead of flat electrodes are widely used in electrochemical systems to boost storage capacities for ions and electrons, to improve the transport of mass and charge, and to enhance reaction rates. Existing porous electrode theories make a number of simplifying assumptions: (i) The charge-transfer rate is assumed to depend only on the local electrostatic potential difference between the electrode matrix and the pore solution, without considering the structure of the double layer (DL) formed in between; (ii) the charge-transfer rate is generally equated with the salt-transfer rate not only at the nanoscale of the matrix-pore interface, but also at the macroscopic scale of transport through the electrode pores. In this paper, we extend porous electrode theory by including the generalized Frumkin-Butler-Volmer model of Faradaic reaction kinetics, which postulates charge transfer across the molecular Stern layer located in between the electron-conducting matrix phase and the plane of closest approach for the ions in the diffuse part of the DL. This is an elegant and purely local description of the charge-transfer rate, which self-consistently determines the surface charge and does not require consideration of reference electrodes or comparison with a global equilibrium. For the description of the DLs, we consider the two natural limits: (i) the classical Gouy-Chapman-Stern model for thin DLs compared to the macroscopic pore dimensions, e.g., for high-porosity metallic foams (macropores >50 nm) and (ii) a modified Donnan model for strongly overlapping DLs, e.g., for porous activated carbon particles (micropores <2 nm). Our theory is valid for electrolytes where both ions are mobile, and it accounts for voltage and concentration differences not only on the macroscopic scale of the full electrode, but also on the local scale of the DL. The model is simple enough to allow us to derive analytical approximations for the steady-state and early transients. We also present numerical solutions to validate the analysis and to illustrate the evolution of ion densities, pore potential, surface charge, and reaction rates in response to an applied voltage