hub.darcs.net :: zaleski -> paris-stable -> changes

The code did not compile.

EPSC is a "clipping parameter" that allows to discard all cells with C_ijk < EPSC It has a default value and a value accessible by the NAMELIST inputvof. The default value was specified in two different places, making confusions highly likely.

Originally the time-step for the fine grid flux computation was set to half of the coarse grid time-step, which was erroneous.

We can plot this velocity vs. time to check if we have a converged terminal velocity while running the Raindrop case.

A very minor change, has no bearing on code functionality.

Earlier versions of the Hybrid mode which utilizes a subgrid had some major issues regarding implementation. These have henceforth been corrected, and tested for both single proc and MPI runs.

NOTE: The 'Rudman' mode still has some asymmetry originating from faulty flux restriction operators between the sub-grid and coarse-grid. Testing and correction is under way.

Brief Description : This mode allows us to advect and reconstruct VOF in a twice finer grid i.e. a volume 8 times smaller than the original control volume. The VOF function is initialized on the fine grid, and for most procedures that require coarse grid values of VOF, appropriate restriction operators enable access to such values derived from fine grid information. This can be implemented both in momentum conserving mode or non-conserving mode. The underlying advection algorithm is based on the paper from Weymouth & Yue (JCP 2010).

How to Use : Add the line "Hybrid = T" to the input script file. Do not set Hybrid and Rudman switches to true at the same time. It can work with both momentum-conserving ("MOMCONS = T" in the 'inputvof' file) or non momentum-conserving modes.

Bugs : The fully coupled momentum-conserving method using the subgrid VOF fluxes (the switch "Rudman" in the 'input' file) STILL has an asymmetry bug (although it does not lead to code failure/blowup). Therefore, avoiding it's usage is recommended for the time being, it shall be debugged in the subsequent versions.

The capillary wave test now shows second-order convergence down to h/lambda = 1/128, the smallest grid size.

Also add the disk on the figure in Tests/Inflowdisk

  The vector image format .svg is now used for html. Unlike pdf it seems supported
  by mozilla. Paper figures are generated in pdf.