Basilisk CFD
Adaptive Cartesian mesh PDE framework
Loading...
Searching...
No Matches
reduced.h File Reference
#include "iforce.h"
#include "curvature.h"
Include dependency graph for reduced.h:

Go to the source code of this file.

Functions

void event_acceleration (void)
 We need the interfacial force module as well as some functions to compute the position of the interface.
 

Variables

coord G = {0.,0.,0.}
 
coord Z = {0.,0.,0.}
 

Function Documentation

◆ event_acceleration()

void event_acceleration ( void  )

We need the interfacial force module as well as some functions to compute the position of the interface.

We overload the acceleration() event to add the contribution of gravity to the interfacial potential \(\phi\).

If \(\phi\) is already allocated, we add the contribution of gravity, otherwise we allocate a new field and set it to the contribution of gravity.

Event: acceleration (i++)

Definition at line 39 of file reduced.h.

References dimension, f, G, scalar::i, phi, position, rho1, rho2, x, coord::x, and Z.

Variable Documentation

◆ G

coord G = {0.,0.,0.}

Reduced gravity

We re-express gravity in two-phase flows as an interfacial force as

\[ -\nabla p + \rho\mathbf{g} = -\nabla p' - [\rho]\mathbf{g}\cdot\mathbf{x}\mathbf{n}\delta_s \]

with \(p'= p - \rho\mathbf{g}\cdot\mathbf{x}\) the dynamic pressure and \(\rho\mathbf{g}\cdot\mathbf{x}\) the hydrostatic pressure. The corresponding potential is

\[ \phi = [\rho]\mathbf{G}\cdot(\mathbf{x} - \mathbf{Z}) \]

with \(\mathbf{G}\) the gravity vector and \(\mathbf{Z}\) an optional reference level.

Definition at line 21 of file reduced.h.

Referenced by event_acceleration().

◆ Z

coord Z = {0.,0.,0.}

Definition at line 21 of file reduced.h.

Referenced by event_acceleration().