科研之友
微信
新浪微博
Facebook
分享链接
摘要
We first develop a new mathematical model for two-fluid interface motion, subjected to the Rayleigh Taylor (RT) instability in two-dimensional fluid flow, which in its simplest form, is given by htt(alpha, t) = Ag Lambda h - (rho)+(sigma)+(rho)- Lambda(3)h Ap alpha(Hh(t)h(t)), where Lambda = Hp alpha and H denotes the Hilbert transform. In this so-called h-model, A is the Atwood number, g is the acceleration, a is surface tension, and p(+/-) denotes the densities of the two fluids. We derive our h-model using asymptotic expansions in the Birkhoff Rott integral-kernel formulation for the evolution of an interface separating two incompressible and irrotational fluids. The resulting h-model equation is shown to be locally and globally well-posed in Sobolev spaces when a certain stability condition is satisfied; this stability condition requires that the product of the Atwood number and the initial velocity field be positive. The asymptotic behavior of these global solutions, as t -> Do, is also described. The h-model equation is shown to have interesting balance laws, which distinguish the stable dynamics from the unstable dynamics. Numerical simulations of the h-model show that the interface can quickly grow due to nonlinearity, and then stabilize when the lighter fluid is on top of the heavier fluid and acceleration is directed downward. In the unstable case of a heavier fluid being supported by the lighter fluid, we find good agreement for the growth of the mixing layer with experimental data in the "rocket rig" experiment of Read and Youngs. We then derive an interface model for RT instability, with a general parameterization z(alpha, t) = (z(1)(alpha,t), z(2)(alpha,t)) such that z satisfies z(tt) = Lambda[A/(broken vertical bar partial derivative alpha Z broken vertical bar 2) Eta(zt center dot (partial derivative alpha(Z))(perpendicular to)Eta((Zt) center dot (partial derivative(alpha Z))(perpendicular to))) + ([P])/(+)(-)(rho)(+rho) + Ag-Z2] (()partial derivative()perpendicular to)(alpha Zeta)/broken vertical bar(partial derivative alpha Zeta)broken vertical bar(2) +(Zeta t) (partial derivative(az)) (perpendicular to) ( ((partial derivative alpha Zeta t)perpendicular to)/broken vertical bar(partial derivative alpha Zeta)broken vertical bar(2) - ((partial derivative alpha Zeta)perpendicular to 2(partial derivative alpha Zeta.partial derivative alpha Zeta t))/(broken vertical bar partial derivative alpha Zeta broken vertical bar 4) ) . This more general RT z-model allows for interface turnover. Numerical simulations of the z-model show an even better agreement with the predicted mixing layer growth for the rocket rig experiment.
- 出版日期2017
