“…Simulations were carried out using in-house code BLUE, ,, which considers a three-dimensional single-field formulation of the Navier–Stokes equations in a Cartesian domain normal∇ · boldu = 0 ρ ( ∂ u ∂ t + u · ∇ u ) = prefix− normal∇ P + ρ boldg + normal∇ · [ μ false( normal∇ boldu + normal∇ u T false) ] + boldF where t , P , u , g , and F denote time, pressure, velocity, gravity, and a local surface force. For clean systems, this term follows a hybrid formulation of the form boldF = ∫ A σ κ boldn δ ( x − boldx f ) .25em normald A where σ is a constant interfacial tension coefficient, κ is the interface curvature, and the 3D Dirac delta function δ( x – x f ) is set to 0 everywhere and unity at the interface, which is located at x = x f . A numerical Heaviside function, scriptH ( x , t ) , is used to define density and viscosity fields throughout the domain, which are respectively given by .25ex2ex ρ ( x , t ) = ρ a + ( ρ normalo<...…”