Demo problem: The (axisymmetric) static free surface bounding a layer of viscous fluid.

We consider an open circular container of unit radius containing a still viscous fluid of prescribed volume that meets the wall of the container at a contact angle . The problem is extremely similar to that described in the two-dimensional static cap tutorial. The exact solution corresponds to a free surface of constant curvature that is the arc of a circle rotated about the axis to give a section of a sphere. The mean curvature of the interface in this problem is , which differs from the two-dimensional problem in which it was .

The only differences between the axisymmetric and the two-dimensional driver codes are that:

- two-dimensional elements are replaced by the equivalent axisymmetric elements;

- specified volume is different;

- swirl (theta) component of velocity is pinned;

- analytic pressure drop is changed.

In other words we make the following changes:

Two-dimensional problem | Axisymmetric problem | |

Bulk Fluid Element | QCrouzeixRaviartElement<2> | AxisymmetricQCrouzeixRaviartElement |

Pseudo-Solid Free Surface Face Element | ElasticLineFluidInterfaceElement | ElasticAxisymmetricFluidInterfaceElement |

Spine Free Surface Face Element | SpineLineFluidInterfaceElement | SpineAxisymmetricFluidInterfaceElement |

Pseudo-Solid Volume Constraint Face Element | ElasticLineVolumeConstraintBoundingElement | ElasticAxisymmetricVolumeConstraintBoundingElement |

Spine Volume Constraint Face Element | SpineLineVolumeConstraintBoundingElement | SpineAxisymmetricVolumeConstraintBoundingElement |

Specific Volume | ||

Analytic pressure drop |

- The formula for the specified volume in the axisymmetric case is the true volume divided by because all the axisymmetric equations are divided by the common factor of . The
`AxisymmetricVolumeConstraintBoundingElement`

class must be used so that the volume is correctly calculated. - The swirl velocity is pinned on the boundaries by including the additional code Bulk_mesh_pt->boundary_node_pt(b,n)->pin(2);

- Confirm that the computed pressure differences agree with the analytic expression. Verify that the interface shape is unaffected by the capillary number, but that the pressure difference across the interface changes in inverse proportion to it. Check that the pressure difference is unaffected by the choice of reference pressure.

- Investigate what happens when the two-dimensional volume constraint elements are used. Explain your result.

- Source files for this tutorial are located in the directory:

demo_drivers/axisym_navier_stokes/axi_static_cap/

- The driver code is:

demo_drivers/axisym_navier_stokes/axi_static_cap/axi_static_cap.cc

A pdf version of this document is available.