This repository contains the codes used for simulating the cases discussed in the manuscript: Lifting a sessile oil drop with an impacting one. We investigate the dynamics of an oil drop impacting an identical sessile drop sitting on a superamphiphobic surface. On this page, I am presenting the code that we used to simulate the process shown in the above video. The results presented here are currently under review in Science Advances. For a detailed documentation on the code included in the manuscript, please visit my Basilisk sandbox.
In the manuscript, offset parameter
Weber number is based on the impact velocity,
Navier Stokes equation for this case:
$$
\partial_tU_i+\nabla\cdot(U_iU_j) =
\frac{1}{\hat{\rho}}\left(-\nabla p + Oh\nabla\cdot(2\hat{\mu}D_{ij}) + \kappa\delta_sn_i\right) + Bog_i
$$
The
Ohnesorge number
Bond number
Note: The subscript
Velocity scale as the intertial-capillary velocity, $$ U_\gamma = \sqrt{\frac{\gamma}{\rho_l R}} $$
For Bond numbers
We refine based on curvatures of the two drops along with the generally used VoF and velocity fields. This ensures that the refinement level along the interface is MAXlevel.
Use the following procedure:
Step 1: Importing the stl file and generating the first dump file
#!/bin/bash
mkdir intermediate
qcc -fopenmp -O2 -Wall dropOnDropImpact.c -o dropOnDropImpact -lm
export OMP_NUM_THREADS=8
./dropOnDropImpact
Step 2: Follow the method described (here). Do not forget to use the dump file generated in the previous step.