next up previous contents
Next: Summary Up: Comparison Between Interfacial Wave Previous: Comparison of Dimensionless Parameters

Numerical Comparisons

The experimental results [10] are shown in figures 8-11 and 8-13

     figure6955
Figure 8-12: Vector plot of the dimensionless velocity tex2html_wrap_inline17479 for wave (b).
Figure 8-11: Vector plot of the dimensionless velocity tex2html_wrap_inline17479 for wave (a).

     figure6971
Figure 8-14: Vector plot of the dimensionless velocity tex2html_wrap_inline17479 for wave (b).
Figure 8-13: Vector plot of the dimensionless velocity tex2html_wrap_inline17479 for wave (a).

at two different times. The simulated results are shown in figures 8-12 and 8-14 for areas with the same dimensionless size. The separation between the simulated results plotted is 16lu. The results plotted are the dimensionless velocities, tex2html_wrap_inline17479 . The axis have been picked so that z = 0 is the mwl, as before, and x = 0 is at a trough. When comparing the simulation results and the experimental results it should be realised that the spacing between the vectors is different. The separation between the experimental results is 0.0125m, giving about 59 points per wavelength in both directions. The simulation results are plotted every 16 sites. That gives 64 vectors per wavelength in the x-direction and 74 vectors per wavelength in the z-direction. Thus, since the dimensionless areas are the same in both plots, the simulated results have slightly more vectors in the horizontal direction and considerably more in the vertical direction. This means that the length of the vectors should be compared and not the closeness of packing.

The velocity profile was found for the experimental and the simulation results. The horizontal profile is shown in figure 8-15 for the results at the wave trough and in figure 8-16 at the wave crest. The vertical velocity profile at tex2html_wrap_inline12645 , the midpoint between the crest and the trough, is shown in figure 8-17.

   figure6994
Figure: The horizontal component of the dimensionless velocity u' = u/c as a function of the dimensionless length tex2html_wrap_inline12737 for the two troughs shown in figure 8-11 tex2html_wrap_inline12845 and figure 8-13 (+). The solid line represents the simulation results.

   figure7004
Figure: The horizontal component of the dimensionless velocity u' = u/c as a function of the dimensionless length tex2html_wrap_inline12737 for the crest shown in figure 8-13 (+). The solid line represents the simulation results.

   figure7013
Figure: The vertical component of the dimensionless velocity w' = w/c as a function of the dimensionless length tex2html_wrap_inline12737 at tex2html_wrap_inline12645 shown in figure 8-11 tex2html_wrap_inline12845 and figure 8-13 (+). The solid line represents the simulation results.

The vertical velocities, displayed in figure 8-17, show good agreement between the simulation results and the experimental results. There is clearly some degree of uncertainty in the experimental results. The horizontal results, shown in figures 8-15 and 8-16, are qualitatively similar, however quantitatively there are some discrepancies. In each of the three sets of measurements the horizontal velocities in one of the fluids agree well, while in the other fluid the magnitude of the experimental results are somewhat smaller than the simulation results. As figure 8-15 demonstrates it is not always the same fluid which contains the discrepancy. This effect can also be seen in figures 8-11 and 8-13 where the magnitude of the horizontal velocity above and below the troughs and crest are different. The difference between the two sets of experimental results in figure 8-15 and the fact that theory predicts only a small difference close to the interface, suggest that the discrepancy is due to some problem with the experimental results rather than the simulations. There are a number of sources which can introduce errors into the experimental results. The most likely cause is the image shifting system which is imposing a horizontal shift velocity onto the flow [94]. This velocity is considerably larger than the flow velocity and may be slightly irregular. Despite these errors the vertical profile is seen to fit well with the experimental results. The horizontal velocity shows the same qualitative features and the comparison is as good as can be expected, given the uncertainties in the experimental results. It is interesting to note that despite the difference in the Reynolds number there is little difference is the shape of the horizontal velocity close to the interface. This could be partly due to the difference in the interface thicknesses.


next up previous contents
Next: Summary Up: Comparison Between Interfacial Wave Previous: Comparison of Dimensionless Parameters

James Buick
Tue Mar 17 17:29:36 GMT 1998