SFEMaNS
version 4.1 (work in progress)
Reference documentation for SFEMaNS

It can be used by the function init_level_set
to initialize the level set functions.
This functions defines a level set for one given interface, one given Fourier mode, one given component (cosine or sine) on a number of nodes, denoted nb_node, of the finite element mesh. We denote by level_set_mesh the finite element mesh used to approximate the level set.
The inputs of this function are the following:
interface_nb
is the label of the interface considered. It is an integer between one and the number of fluids minus one. TYPE
is the component of the level set is computed (cosine or sine). It is an integer between one and two. rr
is a real valued tabular that contains two columns with dimensions (2,nb_node). The tabular rr(1,:) contains the radial cylindrical coordinate of each nodes considered. Respectively, rr(2,:) contains the vertical coordinates of these nodes. We note that the integer nb_node is generally equal to the total number of node level_set_mesh%np. m
is the Fourier mode \(m\) considered. It is an integer. t
is the time at which this term is computed. It is a real number. The output of this function is a real valued tabular vv of dimension SIZE(rr,2).
Here is an exemple where we consider a problem with one level set (so two fluids). The initial interface is \(z=1\). We define a level set \(\varphi\) that depends of the vertical coordinates as follows:
\begin{align} \varphi(z)= \frac{1}{2}(1+\tanh(\frac{z1}{0.05})). \end{align}
One can note that \(\varphi\) is close to zero for z smaller than 0.9, close to one for z larger than 1.1 and equal to 0.5 for z equal to 1.
The corresponding code lines are written as follows.
We refer to the sections Examples with manufactured solutions and Examples on physical problems for more examples.