Electroosmotic Flow Research Articles

Dynamic behavior of electro-osmosis in variable charge soils: Insights on termination and direction reversal

Electro-osmosis offers an effective method for dewatering and remediating low permeability soil. Long-term observations on nonlinear behavior of electro-osmosis and the influencing factors are not commonly reported. Connection between cessation and direction reversal of electro-osmotic flow (EOF), and the evolution of electro-chemical parameters inside of the soil mass thus remains unclear. The dynamic response of EOF in variable charge soil could be significant, whereas the investigations on which are currently lacking. A series of electro-osmotic experiments were performed with two natural variable charge soils. The results indicated that initial electro-osmotic rate was positively proportional to electric current and initial electrical conductivity of the pore fluid, which could be explained by the ion migration model. The dynamic evolution of electro-osmotic rate and electro-chemical parameters corresponding to the solute and pH conditionings at the electrode compartments demonstrated that: 1) coupling effects of non-uniform distribution of voltage gradient and pH determined the magnitude and direction of EOF rate; 2) compared to the final pHIEP value, the bigger, close and smaller values of the novel index “voltage gradient weighed mean of spatial pH″ represented the forward, terminated and reversed EOF respectively; 3) the classical Helmholtz-Smoluchowski model are proved to be more applicable interpreting the coupled nonlinearity of electro-osmosis during the later steady phase. This work would facilitate future research for a comprehensive electro-osmotic model, and provide guidance to condition the initial and boundary conditions in application of electro-osmotic dewatering and electrokinetic remediation.

Effects of multiple relaxation times in the annular flow of pulsatile electro-osmotic flow of a complex biological fluid: blood with low and high cholesterol

This study investigates the electro-osmotic flow of a biological fluid (blood with varying cholesterol levels) in annular flow to simulate a first approximation to arterial occlusion. The fluid´s rheology is characterized by a multi-modal convected Maxwell model equation. The charge density follows the Boltzmann distribution, governing the electrical field. Mathematically, this scenario can be modeled by the Poisson–Boltzmann partial differential equation. Assuming a small zeta potential (less than 25 mV) using the Debye–Huckel approximation and considering a pulsatile electrical field, analytical solutions are derived using the Fourier transform formalism. These solutions, expressed in terms of the modified Bessel function, provide transfer functions for axial velocity and volumetric flow as functions of material parameters represented by characteristic dimensionless numbers. This study further analyzes thermal, electric, inertial, viscoelastic, and various interactions within the plasma, hematocrit, hematocrit–cholesterol, and cholesterol–cholesterol as well as weight concentration through numerical simulations. Finally, the flow and rheology predictions are validated using experimental data on human blood with varying cholesterol levels. The obtained transfer functions reveal that the electric–thermal–viscoelastic effects and the multiple geometric relationships contribute to the dynamic response of the interactions between the input electrical field and output volumetric flow and shear stress functions, leading to and evolution of resonance curves. It is noteworthy that electro-osmotic flow in blood with pathologies associated with low and high cholesterol has been scarcely reported in the literature on rheology. Thus, this work represents a significant contribution to the field.

Electrokinetically controlled mixed convective heat flow in a slit microchannel

AbstractThis study performs a time‐dependent analysis of mixed convection of an incompressible fluid and heat sink/source factor in an upstanding slit superhydrophobic (SHO) microchannel in the involvement of temperature jump and electroosmotic flow conditions. The internal wall of one of the sides in the microchannel is intentionally modified to demonstrate SHO slip and temperature jump conditions. A transverse magnetic effect is introduced in the path of the flow. The steady‐state solutions of the modeled problem have been analytically derived for temperature, velocity, pressure gradient, sheer stress, and heat transfer rate. The derived results are expounded thoroughly with the use of several plots. It is deduced that the elevating mixed convection (Gre), heat source/sink (Qs), Debye–Hückel (K), and nonlinear parameters (N) are observed to increase the fluid flow as time rises, and these effects are all higher when the velocity slip and temperature jump impacts are present. Further, the application of heat‐generating parameters is viewed to encourage the fluid temperature in the microchannel. Finally, the comparison between the current investigation with the previously published findings demonstrates a very good consistency for the limiting cases.

Electroosmotic strengthening of soft clay under different electrification modes

AbstractElectrification mode is one of the important factors affecting the energy consumption and drainage effect of the electroosmosis method. In this work, a series of laboratory tests were performed to explore the impact of continuous changes in the electric field on the migration of ions and electrons, as well as the migration of electroosmotic flow. The results show that different electrification modes had a significant impact on the current and water discharge. An electrode reversed from anode to cathode, from one circuit to another, is beneficial to drainage, while an electrode reversed from cathode to anode, from one circuit to another, is unfavorable to drainage. The drainage effect of the specific cyclic and progressive electroosmosis (CPE) is not as good as the conventional electroosmosis method. However, combining the two electrification modes can further enhance the shear strength of the soil near a designated electrode. The generation of transverse cracks in soil is mainly due to the drainage difference caused by the electroosmotic flow, which is caused by the adjacent electric field at that location. When the electroosmotic flow is strong, longitudinal cracks will occur in the soil along the strongest electroosmotic flow path. In practical in situ engineering applications, the electrification mode of CPE will be more energy‐efficient than conventional electrification modes with the same processing time.

Analytical Solution for Transient Electroosmotic and Pressure-Driven Flows in Microtubes

This study focuses on deriving and presenting an infinite series as the analytical solution for transient electroosmotic and pressure-driven flows in microtubes. Such a mathematical presentation of fluid dynamics under simultaneous electric field and pressure gradients leverages governing equations derived from the generalized continuity and momentum equations simplified for laminar and axisymmetric flow. Velocity profile developments, apparent slip-induced flow rates, and shear stress distributions were analyzed by varying values of the ratio of microtube radius to Debye length and the electroosmotic slip velocity. Additionally, the “retarded time” in terms of hydraulic diameter, kinematic viscosity, and slip-induced flow rate was derived. A simpler polynomial series approximation for steady electroosmotic flow is also proposed for engineering convenience. The analytical solutions obtained in this study not only enhance the fundamental understanding of the electroosmotic flow characteristics within microtubes, emphasizing the interplay between electroosmotic and pressure-driven mechanisms, but also serve as a benchmark for validating computational fluid dynamics models for electroosmotic flow simulations in more complex flow domains. Moreover, the analytical approach aids in the parametric analysis, providing deeper insights into the impact of physical parameters on electroosmotic and pressure-driven flow behavior, which is critical for optimizing device performance in practical applications. These findings also offer insightful implications for diagnostic and therapeutic strategies in healthcare, particularly enhancing the capabilities of lab-on-a-chip technologies and paving the way for future research in the development and optimization of microfluidic systems.

Controlling Electroosmosis in Nanopores Without Altering the Nanopore Sensing Region.

Nanopores are powerful tools for single-molecule sensing of biomolecules and nanoparticles. The signal coming from the molecule to be analyzed strongly depends on its interaction with the narrower sectionof the nanopore (constriction) that may be tailored to increase sensing accuracy. Modifications of nanopore constriction have also been commonly used to induce electroosmosis, that favors the capture of molecules in the nanopore under a voltage bias and independently of their charge. However, engineering nanopores for increasing both electroosmosis and sensing accuracy is challenging. Here it is shown that large electroosmotic flows can be achieved without altering the nanopore constriction. Using continuum electrohydrodynamic simulations, it is found that an external charged ring generates strong electroosmosis in cylindrical nanopores. Similarly, for conical nanopores it is shown that moving charges away from the cone tip still results in an electroosmotic flow (EOF), whose intensity reduces increasing the diameter of the nanopore sectionwhere charges are placed. This paradigm is applied to engineered biological nanopores showing, via atomistic simulations and experiments, that mutations outside the constriction induce a relatively intense electroosmosis. This strategy provides much more flexibility in nanopore design since electroosmosis can be controlled independently from the constriction, which can be optimized to improve sensingaccuracy.

Electroosmotic permeability in kaolinite and CaCO3 poultice mixtures

Electrokinetic treatment of masonry for desalination or electroosmotic dewatering depends on a poultice, in which the electrodes are placed, which fulfills several purposes. Poultice composed of kaolinite and CaCO3 have been shown to have good workability and high pH buffering capacity. In this work, the electroosmotic (EO) permeability is studied in different kaolinite - CaCO3 mixtures. In addition, the effect on EO of using NaCl as a mixing solution is investigated. A special cell is used to test the EO in the specimens. A phenomenological approach, based on the potential gradient and the flux of solution, was used to calculate the EO flow rate and EO permeability coefficient. Results showed that by increasing the concentration of CaCO3 in the poultice mixture, the EO flow rate decreased and the poultice with 80 % CaCO3 and more did not have any EO flow. Furthermore, the ionic strength in the mixing solution decreases the EO flow rate.

Efficient separation of organic anions in beverages using aminosilane-functionalized capillary electrophoresis with contactless conductivity detection

BackgroundCapillary electrophoresis (CE) has the advantage of rapid anion analysis, when employing a reverse electroosmotic flow (EOF). The conventional CE method utilizes dynamic coatings with surfactants like cetyltrimethylammonium bromide (CTAB) in the run buffer to reverse the EOF. However, this method suffers from very slow equilibration leading to drifting effective migration times of the analyte anions, which adversely affects the identification and quantification of peaks. Permanent coating of the capillary surface may obviate this problem but has been relatively little explored. Thus, permanent capillary surface modification by the covalent binding of 3-aminopropyltriethoxysilane (APTES) was studied as an alternative. ResultsThis study investigates the effect of APTES concentration for surface functionalization on EOF mobility, separation efficiency, and reproducibility of anion separation. The performance data was complemented by X-ray photoelectron spectroscopy (XPS) and contact angle (CA) measurements. The XPS measurements showed that the coverage with APTES was dependent on its concentration in the coating solution. The XPS measurements correlated well with the EOF values determined for the capillaries tested. A standard mixture of 21 anions could be baseline separated within 10 min in the capillaries with lower EOF, but not in the capillary with the highest EOF as the residence time of the analytes was too short in this case. Compared to conventional dynamic coating with CTAB, APTES-functionalized capillaries provide faster equilibration and long-term EOF stability. The application of APTES-functionalized capillaries in analyzing different beverages demonstrates the precision, reliability, and specificity in determining organic anions, providing valuable insights of their compositions. SignificanceAPTES coating on capillaries provides a facile approach to achieve a permanent reversal of the stable EOF to determine anions. The control of the coverage via the concentration of the reagent solution allows the tailoring of the EOF to different needs, a faster EOF for less complex samples where resolution is not challenging, while a lower EOF for higher complex samples where the focus is on separation efficiency. This enhancement in efficiency and sensitivity has been applied to analyzing organic acids in several beverages.

Electro-osmotic Flow Generation via a Sticky Ion Action.

Selective transport of ions through nanometer-sized pores is fundamental to cell biology and central to many technological processes such as water desalination and electrical energy storage. Conventional methods for generating ion selectivity include placement of fixed electrical charges at the inner surface of a nanopore through either point mutations in a protein pore or chemical treatment of a solid-state nanopore surface, with each nanopore type requiring a custom approach. Here, we describe a general method for transforming a nanoscale pore into a highly selective, anion-conducting channel capable of generating a giant electro-osmotic effect. Our molecular dynamics simulations and reverse potential measurements show that exposure of a biological nanopore to high concentrations of guanidinium chloride renders the nanopore surface positively charged due to transient binding of guanidinium cations to the protein surface. A comparison of four biological nanopores reveals the relationship between ion selectivity, nanopore shape, composition of the nanopore surface, and electro-osmotic flow. Guanidinium ions are also found to produce anion selectivity and a giant electro-osmotic flow in solid-state nanopores via the same mechanism. Our sticky-ion approach to generate electro-osmotic flow can have numerous applications in controlling molecular transport at the nanoscale and for detection, identification, and sequencing of individual proteins.

Highly alkaline electrokinetic extraction: Characteristics of chromium mobilization, conversion and transport in high alkalinity soil

In site chromium (Cr) contaminated soil characterized by high alkalinity and carbonate content, protons are not effectively targeted for Cr(III) mobilization but rather accelerate the reduction of easily transportable Cr(VI) within the acidification electrokinetic (EK) system. As an alternative, the highly alkaline extraction conditions (HAECs) maintained by anolyte regulation are explored owing to the ability to desorb strong binding Cr(VI) and form anionic Cr(III)-hydroxides (Cr(OH)4−, Cr(OH)52−). The results demonstrate that HAECs were more efficient in mobilizing ions in severe alkalinity and electrical conductivity soil compared to organic acid acidifying extraction conditions (OAECs). Simultaneously, a limited amount of soluble Cr(III) was produced; however, its transportation was hindered and more noticeable in the case of Cr(VI), displaying a distinct retention phase within the intermediate soil chamber. The antagonistic interplay between electromigration and electroosmotic flow was considered the main responsible factor. The conversion intensity of Cr(VI) to Cr(III) was inhibited at HAECs. The promising mobilization and low conversion intensity contributed to total Cr removal. At HAECs, enhanced electromigration and electroosmotic flow combined with a favorable oxidation environment may facilitate in situ delivery of oxidants, offering practical implications for the EK detoxification of high alkalinity site soil contaminated with Cr. The practicability of HAECs is likely to be enhanced when the cost-benefit balance of providing a simultaneous energy supply during site treatment is resolved.

Revealing molecular insights into surface charge and local viscosity in electroosmotic flows

The limitations of the continuum theory in predicting osmotic response at the nanoscale stem from its lack of molecular-level insight into local fluid properties and the interfacial structure of fluid and electrolyte solutions. To overcome this challenge, our study integrates molecular dynamics (MD) simulation with the continuum framework to explore how surface charge and various hydrodynamic properties impact electroosmotic flow (EOF). The failure of continuum theories to account for molecular interactions and geometric boundaries leads to significant disparities between MD simulations and continuum predictions, influenced by local fluid properties and the electric field. Emphasizing the importance of incorporating appropriate local hydrodynamic properties and atomic interface boundary conditions, our findings bridge the gap between MD simulations and continuum EOF predictions. Our computational results and theoretical model, considering surface charge, atomic interface boundaries, and dynamic structure-based hydrodynamic properties, provide crucial insights and guidance for EOF investigations.

Effect of boundary slip on electroosmotic flow in a curved rectangular microchannel

Abstract The aim of this study is to numerically investigate the impact of boundary slip on electroosmotic flow (EOF) in curved rectangular microchannels. Navier slip boundary conditions were employed at the curved microchannel walls. The electric potential distribution was governed by the Poisson–Boltzmann equation, whereas the velocity distribution was determined by the Navier–Stokes equation. The finite-difference method was employed to solve these two equations. The detailed discussion focuses on the impact of the curvature ratio, electrokinetic width, aspect ratio and slip length on the velocity. The results indicate that the present problem is strongly dependent on these parameters. The results demonstrate that by varying the dimensionless slip length from 0.001 to 0.01 while maintaining a curvature ratio of 0.5 there is a twofold increase in the maximum velocity. Moreover, this increase becomes more pronounced at higher curvature ratios. In addition, the velocity difference between the inner and outer radial regions increases with increasing slip length. Therefore, the incorporation of the slip boundary condition results in an augmented velocity and a more non-uniform velocity distribution. The findings presented here offer valuable insights into the design and optimization of EOF performance in curved hydrophobic microchannels featuring rectangular cross-sections.

Slip effects on electroosmotic flow in a microchannel with squeezing wall motion

Slip effects on electroosmotic flow in a microchannel with squeezing wall motion

Progress of highly reproducible capillary electrophoresis

Following rapid developments in capillary electrophoresis (CE), this technology has become an established analytical technique owing to its microscale characteristics, high speed, high efficiency, and versatility. However, the challenges of poor peak stability and/or reproducibility have consistently hindered its wider applications. CE has long been used as a measurement tool for plotting signal intensities versus the migration time; however, the migration time is not an independent variable in CE, but is affected by many direct and indirect parameters, including capillary (length, diameter, and inner surface properties), electric field (or voltage, current, and/or power), temperature, and running buffer (electrolytes, additives, solvents, and their concentration, buffering pH, etc.). These intricacies render the acquisition of reproducible electropherograms difficult. Various studies ranging from those on the early stages of CE development to those on the exploration of three important strategies have been conducted to address this issue. In the first strategy, the CE conditions, especially those parameters that can maintain a stable electro-osmotic flow, are strictly controlled and stabilized to significantly improve peak repeatability. In the second strategy, either the peak position is corrected using internal standards or the peak time is converted into other variables, such as electrophoretic mobility, to offset or eliminate some unstable factors, thereby improving the repeatability and even reproducibility of the peaks; this strategy is useful when plotting signals versus the migration time ratio, correlated migration time, effective mobility, or temperature-correlated mobility. In the third strategy, a new methodology called highly reproducible CE (HRCE) is established using theoretical studies to explore better principles for real-time CE with the aim of the complete removal of the challenge from the root. This strategy includes the development of novel methods that plot electropherograms based on weighted mobility, migrated charge, charge density, or partial differential molar charge density. Similar to ordinary CE approaches, this strategy can also draw electropherograms based on the ratios of these properties. As theoretically predicted, these novel methods can offset or resist changes in critical CE conditions (mainly electric field strength, capillary length and diameter, and/or some buffer parameters such as concentration). Our experimental results demonstrate that given certain prerequisites, a new set of methods can produce highly reproducible electropherograms. This review focuses on the theoretical basis and advancements of HRCE, and elucidates the link between electrophoretic migration/peak expression theories and their impact on reproducibility. Studies on the transformation of time-scale electropherograms in the CE literature are summarized and analyzed in general. However, this review does not directly discuss research on and progress in improving CE repeatability or reproducibility through instrument upgrades, parameter optimization, or practical method refinements.

Unraveling the transformative impact of ternary hybrid nanoparticles on overlapped stenosis with electroosmotic vascular flow kinetics and heat transfer

This study examines how ternary hybrid nanoparticles affect electroosmotic vascular flow kinetics and heat transfer. Through a meticulous exploration of their intricate interplay, this research unveils unprecedented insights into their transformative impact. By comprehensively analyzing the dynamics of vascular flow and thermal behavior under the influence of ternary hybrid nanoparticles, novel advancements are revealed. The assessment is innovative because it incorporates electroosmotic force on blood flow that contains three different nanoparticles (Ti O2, Al2O3 and Si O2). To evaluate the numerical solution, an unraveled approach using the finite element method is employed, ensuring both stability and convergence of the solution. The computed numerical results are presented in graphs and tables, showcasing the relationship between key factors. The comparative analysis uncovers the unparalleled performance and remarkable efficacy of these nanoparticles in enhancing the electroosmotic vascular flow and optimizing heat transfer. The electric field due to the electroosmosis flow interacts with flow pattern and influence the potential flow and vortex formation. This research presents a paradigm shift in the understanding of biomedical engineering and fluid dynamics, offering promising prospects for revolutionizing healthcare technologies and achieving unprecedented levels of thermal management efficiency across diverse applications.

Enhanced migration and degradation of nitrobenzene in heterogeneous porous media using pulsed direct current electrical resistance heating with hydraulic circulation

Electrical resistance heating (ERH) is a promising in-situ technology for heterogeneous organic contaminated site remediation, yet may have low efficiency when treating semi-volatile organic contaminant (SVOC) of relatively high boiling point. Herein, we chose nitrobenzene as a representative SVOC, and proposed an ERH system powered by pulsed direct current (PDC) with simple hydraulic circulation for improved remediation efficiency in heterogeneous media. The proposed PDC-ERH with hydraulic circulation showed overall improvement in heating performance and energy efficiency, as well as migration and removal of nitrobenzene. This new system improved the uniformity of PDC heating and achieved a temperature increase of ~15°C compared to that using conventional alternating current (AC) of same voltage. Nitrobenzene migration out of the low permeability zone (LPZ) was intensified by the dual effects of heat-induced diffusion enhancement and electric field-induced electroosmotic flow, while subsequent removal was enhanced by electrochemical degradation and volatilization. After 96h, the proposed system has a higher nitrobenzene removal from LPZ (> 97.1%) compared to that using AC (84.0%–95.9%). These results suggest PDC heating coupled with hydraulic circulation was a promising approach for heterogeneous site remediation.

Impacts of activation energy and electroosmosis on peristaltic motion of micropolar Newtonian nanofluid inside a microchannel

This study investigates the impact of electroosmosis on the peristaltic flow of unsteady micropolar nanofluid with heat transfer. The findings could enhance the design of peristaltic pumps, potentially improving drug delivery systems, simulations of blood flow in medical devices, and cancer treatments. The fluid under investigation adheres to a micropolar model and flows through a microchannel that exhibits peristalsis along its walls. Moreover, the system is subjected to various external effects, including a uniform magnetic field, the electroosmotic phenomenon, heat absorption, and a chemical reaction with activation energy. Consequently, the problem is mathematically modulated by a system of nonlinear partial differential equations governing the velocity, temperature, and nanoparticle concentration. By employing wave transformation, these governing equations are reduced to ordinary differential equations (ODEs). The reduced equations were solved both analytically, using the homotopy perturbation method, and numerically, using the Runge–Kutta–Merson method. A comparison was made between the solutions, which were found to be closely aligned. Furthermore, a series of figures were employed to provide visual representation and discussion of the implications of the physical properties. The calculations reveal that the electroosmotic flow (EOF) enhances the axial flow of the micropolar fluid along the direction of the applied electric field. It is also observed that the increase in the activation energy (which indicates a low reaction rate) increases the concentration profile whereas the increase in the reaction rate parameter reduces the concentration profile. Additionally, the spin velocity of the particles is diminished by either an increase in the magnetic parameter or the coupling parameter.

Paper Spray Mass Spectrometry with On‐Paper Electrokinetic Manipulations: Part‐Per‐Trillion Detection of Per/Polyfluoroalkyl Substances in Water and Opioids in Urine

AbstractWe developed a simple, paper‐based device that enables sensitive detection by mass spectrometry (MS) without solid phase extraction or other sample preparation. Using glass fiber filter papers within a 3D printed holder, the device employs electrokinetic manipulations to stack, separate, and desalt charged molecules on paper prior to spray into the MS. Due to counter‐balanced electroosmotic flow and electrophoresis, charged analytes stack on the paper and desalting occurs in minutes. One end of the paper strip was cut into a sharp point and positioned near the inlet of a MS. The stacked analyte bands move toward the paper tip with the EOF where they are ionized by paper spray. The device was applied to analysis of PFAS in tap water with sub part‐per‐trillion detection limits in less than ten minutes with no sample pretreatment. Analysis of opioids in urine also occurs in minutes. The crucial parameters to enable stacking, separation, and MS ionization of both positively and negatively charged analytes were determined and optimized. Experimental and computational modeling studies confirm the electrokinetic stacking and analyte transport mechanisms. On‐paper separations were carried out by stacking analyte bands at different locations depending on their electrophoretic mobility, achieving baseline separation in some cases.

Paper Spray Mass Spectrometry with On-Paper Electrokinetic Manipulations: Part-Per-Trillion Detection of Per/Polyfluoroalkyl Substances in Water and Opioids in Urine.

We developed a simple, paper-based device that enables sensitive detection by mass spectrometry (MS) without solid phase extraction or other sample preparation. Using glass fiber filter papers within a 3D printed holder, the device employs electrokinetic manipulations to stack, separate, and desalt charged molecules on paper prior to spray into the MS. Due to counter-balanced electroosmotic flow and electrophoresis, charged analytes stack on the paper and desalting occurs in minutes. One end of the paper strip was cut into a sharp point and positioned near the inlet of a MS. The stacked analyte bands move toward the paper tip with the EOF where they are ionized by paper spray. The device was applied to analysis of PFAS in tap water with sub part-per-trillion detection limits in less than ten minutes with no sample pretreatment. Analysis of opioids in urine also occurs in minutes. The crucial parameters to enable stacking, separation, and MS ionization of both positively and negatively charged analytes were determined and optimized. Experimental and computational modeling studies confirm the electrokinetic stacking and analyte transport mechanisms. On-paper separations were carried out by stacking analyte bands at different locations depending on their electrophoretic mobility, achieving baseline separation in some cases.

Bilaterally Aligned Electroosmotic Flow Generated by Porous Microneedle Device for Dual-Mode Delivery.

Here, a novel porous microneedle (PMN) device with bilaterally aligned electroosmotic flow (EOF) enabling controllable dual-mode delivery of molecules is developed. The PMNs placed at anode and cathode compartments are modified with anionic poly-2-acrylamido-2-methyl-1-propanesulfonic acid and cationic poly-(3-acrylamidopropyl) trimethylammonium, respectively. The direction of EOF generated by PMN at the cathode compartment is, therefore, reversed from cathode to anode, countering the unwanted cathodal suctioning of interstitial fluid caused by reverse iontophoresis. With the bilateral alignment of EOF, the versatility of the proposed device is evaluated by delivering molecules with different charges and sizes using Franz cell. In addition, a 3D printed probe device is developed to ease practical handling and minimize electrical stimulation by integrating two PMNs in closed proximity. Finally, the performance of the integrated probe device is demonstrated by dual delivery of a variety of molecules (methylene blue, rhodamine B, and fluorescein isothiocyanate-dextran) using pig skin and vaccination using mice with delivered ovalbumin.