two-phase-generic.h
📄 View in API Reference (Doxygen) · View on basilisk.fr
Requires: fractions.h
double rho1 = 1., mu1 = 0., rho2 = 1., mu2 = 0.;Auxilliary fields are necessary to define the (variable) specific volume $\alpha=1/\rho$ as well as the cell-centered density.
face vector alphav[];
scalar rhov[];
event defaults (i = 0)
{
alpha = alphav;
rho = rhov;If the viscosity is non-zero, we need to allocate the face-centered viscosity field.
if (mu1 || mu2) {
mu = new face vector;
reset ((scalar *){mu}, 0);
}We add the interface to the default display.
display ("draw_vof (c = 'f');");
}The density and viscosity are defined using arithmetic averages by default. The user can overload these definitions to use other types of averages (i.e. harmonic).
#ifndef rho
# define rho(f) (clamp(f,0.,1.)*(rho1 - rho2) + rho2)
#endif
#ifndef mu
# define mu(f) (clamp(f,0.,1.)*(mu1 - mu2) + mu2)
#endifWe have the option of using some "smearing" of the density/viscosity jump.
#if FILTERED
scalar sf[];
#else
# define sf f
#endif
event tracer_advection (i++)
{When using smearing of the density jump, we initialise *sf* with the vertex-average of *f*.
#ifndef sf
#if dimension <= 2
foreach()
sf[] = (4.*f[] +
2.*(f[0,1] + f[0,-1] + f[1,0] + f[-1,0]) +
f[-1,-1] + f[1,-1] + f[1,1] + f[-1,1])/16.;
#else // dimension == 3
foreach()
sf[] = (8.*f[] +
4.*(f[-1] + f[1] + f[0,1] + f[0,-1] + f[0,0,1] + f[0,0,-1]) +
2.*(f[-1,1] + f[-1,0,1] + f[-1,0,-1] + f[-1,-1] +
f[0,1,1] + f[0,1,-1] + f[0,-1,1] + f[0,-1,-1] +
f[1,1] + f[1,0,1] + f[1,-1] + f[1,0,-1]) +
f[1,-1,1] + f[-1,1,1] + f[-1,1,-1] + f[1,1,1] +
f[1,1,-1] + f[-1,-1,-1] + f[1,-1,-1] + f[-1,-1,1])/64.;
#endif
#endif // !sf
#if TREE
set_prolongation (sf, refine_bilinear);
#endif
}
#include "fractions.h" [api]
event properties (i++)
{
foreach_face() {
double ff = (sf[] + sf[-1])/2.;
alphav.x[] = fm.x[]/rho(ff);
if (mu1 || mu2) {
face vector muv = mu;
muv.x[] = fm.x[]*mu(ff);
}
}
foreach()
rhov[] = cm[]*rho(sf[]);
#if TREE
set_prolongation (sf, fraction_refine);
#endif
}