SOLAVOF Sample Test Problem
I am rewriting the original SOLAVOF Fortran program in C++, but encouter a problem. I try to simulate the same sample test problem in the report (undular bore calculation) in page 44. All parameters are the same (NMAT=1).
At t=0 (page 46), why U = 0.2? Since F is 0, why U is set for an _empty_ cell at t=0? The SETUP routine initialize U [i][j] only when F[I][J] > 0, so it's unlikely U is set by SETUP. I've checked BC subroutine too, but it does not set U either.
The input FLHT (horizontal free surface height) is 1.0f, and del_y of mesh is 0.2, the fluid should occupy the bottom five rows: I=2~6.
At t=0 (page 46), before anything started, why F is not 1 but 9.99999E-01?
B. D. Nichols, C. W. Hirt, and R. S. Hotchkiss, SOLA-VOF: A Solution Algorithm for Transient Fluid Flow with Multiple Free Boundaries, LANL report LA-8355, 1980.
|All times are GMT -4. The time now is 08:22.|