Laplace Law for Surface Tension

Problem Description

Laplace’s law dictates the pressure difference between inside and the outside of a bubble due to surface tension forces. The expression for a 2D problem is:
Δ p = α / r
Where,
α MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaabaaaaaaaaape GaeqySdegaaa@37A8@
Surface tension coefficient
r MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaabaaaaaaaaape GaamOCaaaa@3700@
Final radius of the bubble

The expression for a 3D problem is:

Δ p = 2 α / r

The initial particle arrangement for the 2D simulation as well as the problem geometry are shown in Figure 1, where the bubble and bulk fluid are shown in blue and red, respectively. For 3D simulation, the side lengths of the bubble fluid and bulk fluid regions are the same as the 2D case.
Figure 1. Initial Configuration of the 2D Laplace Law Problem. The blue particles are the bubble fluid phase, while the red is the bulk fluid.


Numerical Setup

Here, a particle spacing of dx = 2*10-3 m is used for both 2D and 3D cases. Both simulations run for five physical seconds. All boundaries are periodic.

Results

The pressure inside the droplet is estimated by averaging values of all the particles belonging to the droplet and comparing them to the average pressure outside of the droplet. The estimated pressure difference Δ p for 2D and 3D cases and their comparison with the analytic predictions are given in Table 1.
Table 1. Comparison of the Analytic Results Against the nanoFluidX
Δ p N F X Δ p t h e o r y Relative error
2-Dimensional 0.452 0.443 2%
3-Dimensional 0.774 0.806 4%
Figure 2 shows the final particle arrangement at t = 5 s for both cases. As expected, the bubble has reached a circular shape in 2D and a spherical shape in 3D.
Figure 2. Final Phase Configuration of the 2D Case at t = 5 s. The colors show bubble and bulk fluid phases.


Figure 3. Final Phase Configuration of the 3D Case at t = 5 s. The colors show particle pressure.