Modelling approach of a near-far-field model for bubble formation and transport

The first modelling approach is given with a convection-diffusion equation in cylindrical coordinate as:

where we assume $u_{near}$ is the solution of the bubble-formation in the near-field and we assume to have Dirichlet-boundary conditions.

Here, we have the benefit and drawbacks of the modelling approach:

• •

  Benefits:

  • –

    The model is simple and fast to compute.

  • –

    The model also allows to discuss a dynamical shape.

• •

  Drawbacks:

  • –

    The shape of the bubble is not preserved, while we assume a static shape.

  • –

    The influence of the speed of motion in the outer normal direction is not possible.

We apply an improved model, that allows to follow the shapes of the bubble, see [12].

The convection-diffusion equation is reformulated in the notation of a level-set equation, which is given as

where ${\bf v}$ is the convection vector and $F_{0}$ is the speed of motion in the outer normal direction. Further $\Omega$ is the computational domain and $T$ is the end time. The initialisation $u({\bf x},0)$ is the results of the near-field computations.

Such equations are wel-known as level-set equations and can be solved like convection-diffusion equations, see [12].

In the following, we apply the explicit different discretization methods in space, while we apply the level-set equation with the explicit time-discretisation and apply upwind methods for the advection and outer normal direction term only in the $x$-direction, the same is also done with the $y$-direction.

We have the following terms:

where we assume $v_{x},v_{y},F_{0}\geq 0$.

We discretize the Level-set equation with the explicit time-discretisation and apply upwind methods for the advection and outer normal direction term in the $x$- and $y$-direction and have the following terms:

where we assume $v_{x},v_{y},F_{0}\geq 0$.

In the following, we discuss the near-field approach.

The near-field model is discussed with respect to the formation of a drop or bubble, see [15] and [13].

The basic modelling idea is based on the so called Young-Laplace equation, see [3] and deals with the following simplied shape of the bubble, see Figure 3.

We deal with the following parameterisation, see [15]:

where $\beta$ is the Bond number, $\sigma$ is the surface tension, $\rho$ is the liquid density, $g$ is the gravity and $R_{t}$ is the curvature of the drop.

The near-field equations are given as:

where $s$ is the arc length along the curve and $\theta$ the angle of elevation for its slope and $\alpha$ is the mono-layer surface tension.

We have the conditions:

where $L$ is the arc length of the bubble which is a-priori unknown and is numerically computed via the boundary value problem. Further, we apply $r=a\rightarrow 0$.

The near-field bubble is computed with the ODE’s given in (18)-(20).

We estimate the characteristic parameters of the ellipse in the Figure 4.

Based on the estimation of the elliptic-parameters, we obtain the curvature of the ellipse:

The ellipse is the curvature of the far-field, which is computed with the level-set method.

Based on the electric-field on smaller time-scale, we have to update the near-field and far-field coupling in a numerical cycle,

The near-field bubble is computed with the ODE’s and with the underlying influence of the E-field. Further the far-field is computed by the transport equation with level-set methods. We couple the coupling via an interpolation between the mesh-free space of the near-field and the grid-space of the far-field with the give diameters of the bubbles and the trajectories of the transport field

The oscillation of the spherical bubble, see [9], [14] and [1], is given as:

where $r_{0}$ is the mean radius, $\omega_{0}$ is the resonance frequency and $r_{\epsilon_{0}}$ amplitude of maximal oscillations of the bubble radius.

The characteristic frequency of the breathing mode with respect to the electric stress (pressure) is given as:

where $\sigma$ is the gas-liquid surface tension, $\rho$ is the liquid density, $k$ is the polytropic exponent ($k=1.4$ for air). $p$ is the pressure, where $p=|p_{0}-p_{E}|$ with $p_{0}$ is the hydrostatic pressure in the fluid far from the bubble and $p_{E}$ is the electric stress (pressure) given of the electric field with

the pressure maximum is given at $\theta=\pi/2$, $\epsilon$ is the peremitivity and ${\bf E}_{0}$ is the uniform electric field.

The numerical cycle is given in Figure 5.

In the following, we discuss the detailed coupling of the near-field and far-field, see Figure 6.

The near-field equations are given as:

or in vectorial form:

where $\tilde{A}(y(t))=\left(\begin{array}[]{c}\cos(\theta)\\ \sin(\theta)\\ -\frac{\sin(\theta)}{r}\end{array}\right)$ and $\tilde{q}(t)=\left(\begin{array}[]{c}0\\ 0\\ \frac{\Delta p}{\alpha}\end{array}\right)$. Further, $y(t)=(r(t),z(t),\theta(t))^{t}$ and $t$ is the arc length along the curve and $r$ is the radius, $z$ is the vertical distance and $\theta$ the angle of elevation for its slope and $\alpha$ is the mono-layer surface tension.

We have to linearise the equations, which are given as:

where $A(t)=\frac{\partial\tilde{A}}{\partial y}|_{y(t_{i})}y(t)$, $q(t)=\tilde{A}(y(t_{i}))-\frac{\partial\tilde{A}}{\partial y}|_{y(t_{i})}y(t_{i})+\tilde{q}(t)$. Further, we have $y=(r,z,\theta)^{t}$ and

We have the conditions:

where $a$ is the radius of the bubble, $L$ is the arc length of the bubble which is a-priori unknown, where we start with $L=L_{0}=\frac{2\pi a^{2}}{4}$ and we go on with $L=L_{1}\leq L_{2}\leq\ldots\leq L_{stop}$.

We have the following linear equation system:

4 where $S_{i}=R_{i}=-h_{i}^{-1}I-\frac{1}{2}A(t_{i+1/2})$, $q_{i}=q(t_{i+1/2}$ with $h_{i}=t_{i+1}-t_{i}$, $t_{i+1/2}=t_{i}+\frac{1}{2}h_{i}$.

Further, we have the matrices and vectors:

• •

  Test example 1:
  We apply the following test-example with the a quarter of a circle, means $L=2\pi a^{2}/4$, where $a=1,2,3,3.5$ and we assume $N=L/\Delta t$, where we assume $\Delta t=2\pi\;0.001$ and $\frac{\Delta p}{\alpha}=0$.

• •

  Test example 2:
  We have $N=L/\Delta t$, where we assume $\Delta t=0.001$ and $L=2$ and $\frac{\Delta p}{\alpha}=0,0.8,1.4$.

The numerical results of bubble-formation is given in Figure 7.

In the following, we couple the near-field and far-field computations.

We have the following setting, see Figure 8.

We apply the following parameters:

• •

  Input-parameters of the near-field bubbles:

  • –

    $r_{1}(0)=0.01,\;z_{1}(0)=0,\;r_{1}(L)=0,\;L=1,\;\frac{\Delta p_{1}}{\alpha_{1}}=0.8$,

  • –

    $r_{2}(0)=0.01,\;z_{2}(0)=0,\;r_{2}(L)=0,\;L=2,\;\frac{\Delta p_{2}}{\alpha_{2}}=0.8$.

• •

  We compute the bubbles based on the near-field code and obtain the ellipse-diameters $a_{bubble},\;b_{bubble}$.

• •

  We initialise the two ellipses for the far-field computations given as:

  • –

    $(x-20)^{2}+((y-100)*a_{bubble,1}/b_{bubble,1})^{2}-a_{bubble,1}^{2}$,

  • –

    $(x-70)^{2}+((y-90)*a_{bubble,2}/b_{bubble,2})^{2}-a_{bubble,2}^{2}$.

The numerical results of near-field bubble-formation is given in Figure 9.

The numerical results of the near-far-field coupled bubble-transport code, which is given in Figure 10.

In the following, we extend the near-field and far-field computations with 10 bubbles. We also apply the decomposition of near-field and far-field computations as given in Figure 8.

We apply the following parameters:

• •

  Computation of the near-field bubbles (a representing bubble is computed)

  • –

    Input-parameters of the near-field bubbles computation:

    • *

      Bubble 1:
      $r_{1}(0)=0.01,\;z_{1}(0)=0,\;r_{1}(L)=0,\;L=2,\;\frac{\Delta p_{1}}{\alpha_{1}}=0.2$,

    • *

      Bubble 2:
      $r_{2}(0)=0.01,\;z_{2}(0)=0,\;r_{2}(L)=0,\;L=2,\;\frac{\Delta p_{2}}{\alpha_{2}}=0.4$.

    • *

      Bubble 3:
      $r_{3}(0)=0.01,\;z_{3}(0)=0,\;r_{3}(L)=0,\;L=2,\;\frac{\Delta p_{3}}{\alpha_{3}}=0.6$.

    • *

      Bubble 4:
      $r_{4}(0)=0.01,\;z_{4}(0)=0,\;r_{4}(L)=0,\;L=2,\;\frac{\Delta p_{4}}{\alpha_{4}}=0.8$.

    • *

      Bubble 5:
      $r_{5}(0)=0.01,\;z_{5}(0)=0,\;r_{5}(L)=0,\;L=2,\;\frac{\Delta p_{5}}{\alpha_{5}}=1.0$.

    • *

      Bubble 6:
      $r_{6}(0)=0.01,\;z_{6}(0)=0,\;r_{6}(L)=0,\;L=2,\;\frac{\Delta p_{6}}{\alpha_{6}}=1.2$.

    • *

      Bubble 7:
      $r_{7}(0)=0.01,\;z_{7}(0)=0,\;r_{7}(L)=0,\;L=2,\;\frac{\Delta p_{7}}{\alpha_{7}}=1.4$.

    • *

      Bubble 8:
      $r_{8}(0)=0.01,\;z_{8}(0)=0,\;r_{8}(L)=0,\;L=2,\;\frac{\Delta p_{8}}{\alpha_{8}}=1.6$.

    • *

      Bubble 9:
      $r_{9}(0)=0.01,\;z_{9}(0)=0,\;r_{9}(L)=0,\;L=2,\;\frac{\Delta p_{9}}{\alpha_{9}}=1.8$.

    • *

      Bubble 10:
      $r_{10}(0)=0.01,\;z_{10}(0)=0,\;r_{10}(L)=0,\;L=2,\;\frac{\Delta p_{10}}{\alpha_{10}}=2.0$.

  • –

    Output-parameters of the near-field bubble computation:

    • *

      Bubble 1:
      $a_{bubble_{1}}=0.5141,\;b_{bubble_{1}}=0.9926$.

    • *

      Bubble 2:
      $a_{bubble_{2}}=0.5219,\;b_{bubble_{2}}=0.9718$.

    • *

      Bubble 3:
      $a_{bubble_{3}}=0.5295,\;b_{bubble_{3}}=0.9506$.

    • *

      Bubble 4:
      $a_{bubble_{4}}=0.5369,\;b_{bubble_{4}}=0.9289$.

    • *

      Bubble 5:
      $a_{bubble_{5}}=0.5443,\;b_{bubble_{5}}=0.9067$.

    • *

      Bubble 6:
      $a_{bubble_{6}}=0.5514,\;b_{bubble_{6}}=0.8841$.

    • *

      Bubble 7:
      $a_{bubble_{7}}=0.5584,\;b_{bubble_{7}}=0.8612$.

    • *

      Bubble 8:
      $a_{bubble_{8}}=0.5651,\;b_{bubble_{8}}=0.8382$.

    • *

      Bubble 9:
      $a_{bubble_{9}}=0.5717,\;b_{bubble_{9}}=0.8154$.

    • *

      Bubble 10:
      $a_{bubble_{10}}=0.5780,\;b_{bubble_{10}}=0.7928$.

  • –

    Ellipse: $(x-x_{bubble_{i}})^{2}+((y-y_{bubble_{i}})*a_{bubble}/b_{bubble})^{2}-a_{bubble}^{2}$,
    where $(x_{bubble_{i}},y_{bubble_{i}})$ is the origin of the $i$-th bubble.

• •

  Computation of the far-field bubbles (level-set initialisation):

  • –

    Parameterisation of the level-set initial-function, e.g., two bubbles:

    $\displaystyle\phi_{0}(x,y)=\left\{\begin{array}[]{cc}(x-x_{bubble_{1}})^{2}+((y-y_{bubble_{1}})\frac{a_{bubble_{1}}}{b_{bubble_{1}}})^{2}-a_{bubble_{1}}^{2},&\\ a_{x}\leq x\leq 50,\;a_{y}\leq x\leq b_{y},&\\ (x-x_{bubble_{2}})^{2}+((y-y_{bubble_{2}})\frac{a_{bubble_{2}}}{b_{bubble_{2}}})^{2}-a_{bubble_{2}}^{2},&\\ 50\leq x\leq b_{x},\;a_{y}\leq x\leq b_{y},&\\ \end{array}\right.$

    where $(x_{bubble_{1}},y_{bubble_{1}})=(20,50)$, $(x_{bubble_{2}},y_{bubble_{2}})=(80,50)$ with the coordinates of the grid $(a_{x},a_{y})=(0,0)$ and $b_{x},b_{y})=(100,200)$.

The numerical results of the formation of the bubbles are given in Figures 11.

The numerical results of the near-far-field coupled bubble-transport code, which is given in Figures 12 and 13.

In the following, we simulate a bubble filled with air in a electrical field, see [14].

We have the following setting of the influenced bubble, see Figure 14.

We apply the following near-field equation based on an extension of the pressure-term with an E-field.

The near-field equations are given as:

where $s$ is the arc length along the curve and $\theta$ the angle of elevation for its slope and $\alpha$ is the mono-layer surface tension.

Further, we have

where $p_{tube}=p_{E}=\frac{9}{8}|{\bf E}_{0}|^{2}\;\sin(\theta)$ is the pressure based on the electrical field of the tube and we obtain:

We apply the following parameters:

• •

  Computation of the near-field bubbles (a representing bubble is computed)

  • –

    Input-parameters of the near-field bubbles computation:

    • *

      Bubble 1:
      $r_{1}(0)=0.01,\;z_{1}(0)=0,\;r_{1}(L)=0,\;L=2,\;\frac{\Delta p_{1}}{\alpha_{1}}=0.2$,

    • *

      Bubble 2:
      $r_{2}(0)=0.01,\;z_{2}(0)=0,\;r_{2}(L)=0,\;L=2,\;\frac{\Delta p_{2}}{\alpha_{2}}=0.4$.

    • *

      Bubble 3:
      $r_{3}(0)=0.01,\;z_{3}(0)=0,\;r_{3}(L)=0,\;L=2,\;\frac{\Delta p_{3}}{\alpha_{3}}=0.6$.

    • *

      Bubble 4:
      $r_{4}(0)=0.01,\;z_{4}(0)=0,\;r_{4}(L)=0,\;L=2,\;\frac{\Delta p_{4}}{\alpha_{4}}=0.8$.

    • *

      Bubble 5:
      $r_{5}(0)=0.01,\;z_{5}(0)=0,\;r_{5}(L)=0,\;L=2,\;\frac{\Delta p_{5}}{\alpha_{5}}=1.0$.

    • *

      Bubble 6:
      $r_{6}(0)=0.01,\;z_{6}(0)=0,\;r_{6}(L)=0,\;L=2,\;\frac{\Delta p_{6}}{\alpha_{6}}=1.2$.

    • *

      Bubble 7:
      $r_{7}(0)=0.01,\;z_{7}(0)=0,\;r_{7}(L)=0,\;L=2,\;\frac{\Delta p_{7}}{\alpha_{7}}=1.4$.

    • *

      Bubble 8:
      $r_{8}(0)=0.01,\;z_{8}(0)=0,\;r_{8}(L)=0,\;L=2,\;\frac{\Delta p_{8}}{\alpha_{8}}=1.6$.

    • *

      Bubble 9:
      $r_{9}(0)=0.01,\;z_{9}(0)=0,\;r_{9}(L)=0,\;L=2,\;\frac{\Delta p_{9}}{\alpha_{9}}=1.8$.

    • *

      Bubble 10:
      $r_{10}(0)=0.01,\;z_{10}(0)=0,\;r_{10}(L)=0,\;L=2,\;\frac{\Delta p_{10}}{\alpha_{10}}=2.0$.

  • –

    Electrical field parameters are given as:

    $\alpha=0.1$, $\rho=0.1,\;g=9.81,\;|{\bf E}_{0}|^{2}=0.1$.

  • –

    Output-parameters of the near-field bubble computation (formation):

    • *

      Bubble 1:
      $a_{bubble_{1}}=0.5141,\;b_{bubble_{1}}=0.9926$.

    • *

      Bubble 2:
      $a_{bubble_{2}}=0.5219,\;b_{bubble_{2}}=0.9718$.

    • *

      Bubble 3:
      $a_{bubble_{3}}=0.5295,\;b_{bubble_{3}}=0.9506$.

    • *

      Bubble 4:
      $a_{bubble_{4}}=0.5369,\;b_{bubble_{4}}=0.9289$.

    • *

      Bubble 5:
      $a_{bubble_{5}}=0.5443,\;b_{bubble_{5}}=0.9067$.

    • *

      Bubble 6:
      $a_{bubble_{6}}=0.5514,\;b_{bubble_{6}}=0.8841$.

    • *

      Bubble 7:
      $a_{bubble_{7}}=0.5584,\;b_{bubble_{7}}=0.8612$.

    • *

      Bubble 8:
      $a_{bubble_{8}}=0.5651,\;b_{bubble_{8}}=0.8382$.

    • *

      Bubble 9:
      $a_{bubble_{9}}=0.5717,\;b_{bubble_{9}}=0.8154$.

    • *

      Bubble 10:
      $a_{bubble_{10}}=0.5780,\;b_{bubble_{10}}=0.7928$.

  • –

    Output-parameters of the near-field bubble computation (in the E-field):

    • *

      Bubble 1:
      $a_{bubble_{1}}=0.5558,\;b_{bubble_{1}}=0.8698$.

    • *

      Bubble 2:
      $a_{bubble_{2}}=0.5626,\;b_{bubble_{2}}=0.8468$.

    • *

      Bubble 3:
      $a_{bubble_{3}}=0.5693,\;b_{bubble_{3}}=0.8239$.

    • *

      Bubble 4:
      $a_{bubble_{4}}=0.5757,\;b_{bubble_{4}}=0.8013$.

    • *

      Bubble 5:
      $a_{bubble_{5}}=0.5819,\;b_{bubble_{5}}=0.7789$.

    • *

      Bubble 6:
      $a_{bubble_{6}}=0.5878,\;b_{bubble_{6}}=0.7571$.

    • *

      Bubble 7:
      $a_{bubble_{7}}=0.5935,\;b_{bubble_{7}}=0.7359$.

    • *

      Bubble 8:
      $a_{bubble_{8}}=0.5990,\;b_{bubble_{8}}=0.7153$.

    • *

      Bubble 9:
      $a_{bubble_{9}}=0.6042,\;b_{bubble_{9}}=0.6953$.

    • *

      Bubble 10:
      $a_{bubble_{10}}=0.6091,\;b_{bubble_{10}}=0.6760$.

  • –

    Ellipse: $(x-x_{bubble_{i}})^{2}+((y-y_{bubble_{i}})*a_{bubble}/b_{bubble})^{2}-a_{bubble}^{2}$,
    where $(x_{bubble_{i}},y_{bubble_{i}})$ is the origin of the $i$-th bubble.

• •

  Computation of the far-field bubbles (level-set initialisation):

  • –

    Parameterisation of the level-set initial-function, e.g., two bubbles:

    $\displaystyle\phi_{0}(x,y)=\left\{\begin{array}[]{cc}(x-x_{bubble_{1}})^{2}+((y-y_{bubble_{1}})\frac{a_{bubble_{1}}}{b_{bubble_{1}}})^{2}-a_{bubble_{1}}^{2},&\\ a_{x}\leq x\leq 50,\;a_{y}\leq x\leq b_{y},&\\ (x-x_{bubble_{2}})^{2}+((y-y_{bubble_{2}})\frac{a_{bubble_{2}}}{b_{bubble_{2}}})^{2}-a_{bubble_{2}}^{2},&\\ 50\leq x\leq b_{x},\;a_{y}\leq x\leq b_{y},&\\ \end{array}\right.$

    where $(x_{bubble_{1}},y_{bubble_{1}})=(20,50)$, $(x_{bubble_{2}},y_{bubble_{2}})=(80,50)$ with the coordinates of the grid $(a_{x},a_{y})=(0,0)$ and $b_{x},b_{y})=(100,200)$.