MOF test case with filaments and periodic boundary conditions
We simulate three fluids in a box with filaments of coordinates (0,0) and (1,1).
All the fluids considered have the same properties :
\(\rho\) | \(\mu\) | \(\sigma\) | C | \(\alpha\) | \(T_{0}\) |
---|---|---|---|---|---|
1 | 1 | 1 | 1 | 1 | 1 |
where \(\rho\) is the density, \(\mu\) the dynamic viscosity coefficient, \(\sigma\) the conductivity, C the specific heat capacity, \(\alpha\) the thermal expansion coefficient and \(T_{0}\) the reference temperature.
The velocity field is uniform unidirectional. It is set as: \( u=(\pm \dfrac{1}{2},\pm \dfrac{1}{2})\)
We consider periodic boundary conditions along both horizontal and vertical axis.
The numerical method considered for Phase advection is MOF with filaments.
A regular mesh with \(\dfrac{L}{\Delta x}=32\) for each direction is considered.
The time step is fixed and equals to \(0.01 \ s\). The final time is set to \(2 \ s\) and corresponds to one period.
The use_analytic variable is used to toggle the analytic reconstruction. The tolerance angle and derivative are respectively \(10^{-15}\) and \(0\).
Two variables are included: sign and direction. sign allows to change the velocity direction. direction variable concerns the positioning of the filaments.
Figure 1 shows the volume fraction of the main fluids at three different times with analytic reconstruction.
Analytic reconstruction | direction | sign | Symmetric difference area |
---|---|---|---|
enable | 1 | 1.0 | 5.40995354406858999E-016 |
enable | 1 | -1.0 | 6.77361330403222784E-016 |
enable | 2 | 1.0 | 1.49848843384206664E-013 |
enable | 2 | -1.0 | 1.49839585951506342E-013 |
disable | 1 | 1.0 | 5.40995354406858999E-016 |
disable | 1 | -1.0 | 6.77361330403222784E-016 |
disable | 2 | 1.0 | 1.49848843384206664E-013 |
disable | 2 | -1.0 | 1.49839585951506342E-013 |