A Robust Truncated-Domain Approach for Cone--Jet Simulations in Electrospinning and Electrospraying

This paper presents a robust, parameter-free truncated-domain framework for electrohydrodynamic simulations of cone-jet flows that leverages inexpensive full-domain electrostatic data to impose accurate boundary conditions, thereby significantly reducing computational costs while maintaining high predictive accuracy without requiring empirical tuning.

The Gist

Imagine you are trying to paint a tiny, perfect masterpiece on a single grain of sand, but you are standing in the middle of a massive football stadium.

That is essentially the challenge scientists face when studying electrospinning and electrospraying. These are high-tech processes used to make ultra-thin fibers (for things like medical filters or drug delivery) or tiny particles. They work by shooting a stream of liquid out of a tiny needle using a strong electric charge.

The Problem: The "Stadium" vs. The "Grain of Sand"

In a real experiment, the needle is tiny (about the width of a human hair), but the target collector plate is far away (like the size of a room).

• The Needle: Needs to be simulated with extreme detail, like looking at it under a microscope.
• The Distance: Needs to be simulated as a huge empty space.

If you try to simulate the whole thing on a computer, you need a grid so massive that it would take supercomputers years to crunch the numbers. It's like trying to map every single blade of grass in a stadium just to study the texture of one specific blade.

The Old Solution: The "Guess-and-Check" Map

To save time, scientists used to cut the simulation short. They would simulate just the area near the needle (the "grain of sand") and pretend the rest of the stadium didn't exist. To do this, they used a mathematical formula (a "map") to guess what the electric field looked like at the edge of their cut-off zone.

The Flaw: This old map was like a tourist's sketch of a city. It was close, but not accurate.

• It often got the "traffic" (electric field) wrong right near the needle.
• To make the simulation work, scientists had to tweak the numbers manually until the result looked right. This is like adjusting a GPS until it tells you the route you want to take, rather than the route that actually exists.
• The Catch: If you didn't already know the answer from an experiment, you couldn't use this method to predict new results. It was a tool for checking, not for discovering.

The New Solution: The "High-Res Snapshot"

The authors of this paper, Ghanashyam K. C. and his team, came up with a clever, two-step trick to fix this.

Step 1: The Cheap Snapshot
First, they run a very simple, fast computer simulation of just the electricity (ignoring the liquid flow for a moment) across the whole "stadium." Because electricity is easier to calculate than moving liquid, this takes very little time and computer power.

• Analogy: Imagine taking a high-resolution photo of the wind patterns in the stadium before you even start the race.

Step 2: The Perfect Boundary
They take that photo (the exact electric field data) and paste it onto the edges of their small, cut-off simulation.

• Analogy: Instead of guessing what the wind is doing at the edge of your small garden, you look at the photo of the whole stadium and say, "Okay, at this exact spot, the wind is blowing at 5 mph from the north."

Now, when they simulate the liquid jet in their small, fast domain, it behaves exactly as if it were in the huge stadium, because the "wind" (electric field) pushing on it is perfectly accurate.

Why This is a Big Deal

1. No More Guessing: You don't need to know the answer beforehand. The computer figures out the physics naturally.
2. Super Fast: It runs on a much smaller computer domain, saving massive amounts of time and money.
3. Predictive Power: Because it doesn't rely on "tweaking" numbers to match old experiments, scientists can now use it to design new machines and processes they've never tried before.

The Bottom Line

Think of the old method as trying to draw a map of a forest by looking at a blurry sketch and guessing where the trees are. The new method is like taking a drone photo of the whole forest, zooming in on the specific tree you care about, and using that perfect detail to study it.

This new approach allows scientists to design better medical filters, more efficient drug sprays, and advanced materials much faster and more accurately than ever before.

Technical Summary

1. Problem Statement

Electrospinning and electrospraying are critical techniques for producing sub-micron fibers and particles. The core physical phenomenon governing these processes is the cone-jet mode, where a pendant liquid drop deforms into a Taylor cone under high electric fields, ejecting a thin jet.

• Computational Challenge: Direct Numerical Simulations (DNS) of these processes are computationally prohibitive due to a massive disparity in length scales. The needle-to-collector distance ($H'$) is typically 15–20 cm, while the cone-jet region is only 0.6–1 cm, and the needle inner diameter is ~250–600 $\mu$m. Resolving the fine interfacial features within a full domain requires grid resolutions of $O(1-5)$ $\mu$m, making full-domain 3D simulations extremely expensive.
• Limitations of Existing Truncated-Domain Methods: To reduce cost, researchers often use Truncated-Domain Simulations (TDS), simulating only the local cone-jet region. However, existing TDS approaches rely on analytical expressions (specifically the Jones-Thong model) to approximate far-field electric boundary conditions.
  • These analytical models systematically underestimate the electric field near the needle tip.
  • They require empirical tuning of parameters (e.g., coefficient $K_V$) based on prior experimental data or full-domain results to match the interface shape.
  • This reliance on prior knowledge destroys the predictive capability of the simulations, making them unsuitable for exploring uncharted parameter regimes.

2. Methodology

The authors propose a new, general framework for TDS that eliminates empirical tuning by leveraging the low computational cost of purely electrostatic simulations.

A. Governing Equations & Setup

• Fluid Model: The study uses the Taylor-Melcher leaky-dielectric model coupled with the Navier-Stokes equations.
• Numerical Solver: The simulations are performed using an extended OpenFOAM solver (interFoam) with the MULES scheme for Volume-of-Fluid (VOF) interface tracking.
• Physics: The model accounts for surface tension, viscous stresses, inertia, and electrostatic body forces (Maxwell stress). It employs a specific interpolation scheme for conductivity and permittivity at the interface to minimize spurious charge leakage.
• Test Case: The methodology is demonstrated using 1-octanol (Newtonian fluid) in an axisymmetric configuration, validated against experimental data from Herrada et al.

B. The Proposed Truncated-Domain Strategy
Instead of using analytical approximations, the authors propose a two-step workflow:

1. Full-Domain Electrostatic Simulation (FDS-ES): A computationally inexpensive simulation is run on the full geometry (needle to collector) without fluid dynamics. This captures the exact electric potential ($\Phi$) and electric field ($E$) distributions surrounding the needle, accounting for the specific geometry and boundary conditions.
2. Data Extraction and Fitting: The electric field and potential data are extracted from the boundaries of a proposed truncated domain (top, lateral, and bottom). These data are smoothed and fitted using three-term Gaussian profiles to ensure continuity and differentiability.
3. Truncated EHD Simulation: The fitted profiles are imposed as Dirichlet (potential) and Neumann (field gradient) boundary conditions on the truncated domain. A full electrohydrodynamic simulation is then run on this small domain.

C. Convergence Studies

• Lateral Boundary: Pure electrostatic simulations determined that the lateral boundary must be placed at $L \geq 1.5 H'$ to eliminate artificial confinement effects on the electric field.
• Truncation Size: The proposed method achieves convergence for a truncated domain size as small as $6R_i \times 6R_i$ (where $R_i$ is the needle radius), significantly smaller than the sizes required by analytical methods.

3. Key Contributions

1. Elimination of Tunable Parameters: The proposed method removes the need for empirical coefficients (like $K_V$) or prior knowledge of the cone-jet shape. It is purely physics-based.
2. Exact Far-Field Boundary Conditions: By deriving boundary conditions from exact electrostatic solutions rather than analytical approximations, the method accurately reproduces the electric field strength near the needle tip, which is critical for jet stability.
3. Predictive Capability: The framework allows for the simulation of cone-jet configurations in unexplored parameter regimes without requiring experimental calibration.
4. Generality: While demonstrated in 2D axisymmetry, the method is geometry-agnostic and can be extended to 3D and non-Newtonian fluids without modifying the electrostatic reconstruction procedure.

4. Results and Validation

The authors compared three approaches:

• FDS-EHD: Full-domain electrohydrodynamic simulation (Reference).
• TDS-EHD-JT: Truncated-domain using Jones-Thong analytical expressions (Existing method).
• TDS-EHD-P: Truncated-domain using the Proposed method.

Key Findings:

• Interface Morphology: The proposed method (TDS-EHD-P) accurately reproduces the cone-jet shape, matching both full-domain simulations and experimental data. In contrast, the Jones-Thong approach (TDS-EHD-JT) showed high sensitivity to domain size and required specific parameter tuning ($K_V=0.56$ for $H'/R_i=40$) to achieve a stable cone, even when the physical $H'/R_i$ was different (10).
• Physical Quantities: The proposed method significantly outperformed the analytical approach in predicting secondary quantities:
  • Currents: Captured the spatial variation of conduction and convection currents more accurately.
  • Charge Density: The charge distribution profile closely matched the full-domain results, whereas the analytical method showed deviations.
  • Maxwell Stresses & Velocity: The proposed method yielded velocity fields and Maxwell stress profiles with markedly reduced discrepancies compared to the analytical truncation.
• Computational Efficiency: The proposed method achieves convergence at a domain size of $6R_i \times 6R_i$, whereas analytical methods often required larger domains ($10R_i \times 15R_i$) and still failed to capture accurate physics without tuning. This leads to a substantial reduction in computational cost.

5. Significance

This work provides a robust, predictive, and cost-effective framework for studying electrohydrodynamic flows. By decoupling the expensive fluid dynamics from the electrostatics and using the latter to inform the boundary conditions of the former, the authors have solved a long-standing issue in the field: the inability to simulate cone-jet dynamics in truncated domains without empirical calibration.

This approach enables researchers to:

• Perform parametric studies and optimization of electrospinning/electrospraying systems efficiently.
• Study complex instabilities and transition regimes without the prohibitive cost of full-domain simulations.
• Extend simulations to 3D geometries and complex fluid rheologies with confidence in the boundary condition accuracy.

In summary, the paper replaces "approximate analytical boundary conditions" with "exact numerical boundary conditions derived from cheap electrostatic solves," offering a superior alternative for both fundamental research and industrial application design.