NUMERICAL SOLUTION OF WAVY-STRATIFIED FLUID-FLUID FLOW WITH THE ONE-DIMENSIONAL TWO-FLUID MODEL: STABILITY, BOUNDEDNESS, CONVERGENCE AND CHAOS

W FULLMER; A CLAUSSE; A VAIDHEESWARAN; M LÓPEZ DE BERTODANO

Chicago

Congreso; ASME 2014 4th Joint US-European Fluids Engineering Division Summer; 2014

ASME

In this paper the one-dimensional two-fluid model is used to dynamically simulate slightly inclined fluid-fluid flow in a rectangular channel. By that, it is specifically meant that the solutions exhibit a wavy pattern arising from the inherent instability of the model. The conditions and experimental data of Thorpe are used for comparison. The linear instability of the model is regularized, i.e., made well-posed, with surface tension and axial turbulent stress with a simple turbulent viscosity model. Nonlinear analysis in an infinite domain demonstrates for the first time one-dimensional two-fluid model chaotic behavior in addition to limit cycle behavior and asymptotic stability. The chaotic behavior is a consequence of the linear instability (the long wavelength energy source) the nonlinearity (the energy transfer mechanism) and the viscous dissipation (the short wavelength energy sink). Since the model is chaotic, solutions exhibit sensitive dependence on initial conditions which results in non-convergence of particular solutions with grid refinement. However, even chaotic problems have invariants and the ensemble averaged water void fraction amplitude spectrum is used to demonstrate convergence and make comparisons to the experimental data.