arxiv:physics/ v1 [physics.ed-ph] 25 Apr 2000

Transcription

1 Matching Conditions on Capillar Ripples: Polariation Estudio de las condiciones de empalme para las oscilaciones de intercara entre dos fluidos. Polariación. arxiv:phsics/ v1 [phsics.ed-ph] 25 Apr 2000 Arek H. Rodrígue, J. Marín-Antuña and H. Rodrígue-Coppola Dpto. de Física Teórica, Fac. de Física, Universidad de la Habana, C. de la Habana 10400, Cuba Februar 28, 2018) The matching conditions at the interface between two non-mixed fluids at rest are obtained directl using the equation of movement of the whole media. This is a non-usual point of view in hdrodnamics courses and our aim is to fix ideas about the intrinsic information contained in the matching conditions, on fluids in this case. Afterward, it is analed the polariation of the normal modes at the interface and it is shown that this information can be achieved through a phsical analsis and reinforced later b the matching conditions. A detailed analsis of the matching conditions is given to understand the role that plas the continuit of the stress tensor through the interface on the phsics of the surface particle movement. The main importance of the viscosit of each medium is deduced. En el presente trabajo se utilian las ecuaciones de movimiento de todo el sistema compuesto por dos fluidos no miscibles en reposo relativo para obtener las condiciones de empalme en la intercara. Este procedimiento es inusual el objetivo de hacerlo en esa forma es fijar ideas sobre la información que está contenida en dichas condiciones de empalme, en este caso en fluidos. Posteriormente es analiada la polariación de los modos normales en la interface se demuestra que estas características pueden ser obtenidas directamente de análisis físico de las condiciones de empalme g, e I. INTRODUCTION The boundar conditions constitute the ke feature in an theor of surface waves. It is through them that one introduces in the analsis the phsical consequences of specific surface effects, such as local changes of mass or stresses, thus going beond the simple case in which one merel matches two semi-infinite bulk media at an interface which, in particular, can be a free surface). In the last ears, increased attention has been paid to the properties of capillar waves b phsicists and chemists [1 8]. Ripples represent one of the few cases in which the relation between the dnamical properties of a surface and liquid flows can be predicted completel. The stud of capillar ripples has clarified which properties of a liquid surface determine the surface s resistance against deformation. The boundar conditions for the stress at the interface are derived from the principle that the forces acting upon such an interfacial element result not onl from viscous stresses in the liquid but also from stresses existing in the deformed interface. The cause of the difference is that an interface, unlike a three dimensional liquid, can not enjo the propert of incompressibilit [1]. Work done on an element of liquid is partiall degrade into heat b viscous friction and partiall transformed into kinetic energ which is transmitted to adjoining elements. Work done on an interfacial element leads, at least partiall, to an increment of the surface potential energ. It is this potential energ of the deformed interface which enables the whole sstem, including the interface, to carr out an oscillator motion. In this article, first, it will be shown a non-traditional wa to obtain the matching conditions at the interface between two non-mixed fluid at rest, considering the equation of motion of the whole media directl. Usuall, universit courses do not use this approach and state the boundar conditions from outside the constitutive equations governing the studied problem. It is important from a pedagogical point of view to evidence that in fact, the matching conditions at the interfaces are contained, almost in all cases, in the equation of movement for the whole sstem taken as the composition of all media. The arek@ff.oc.uh.cu 1

2 other cases arise when the fluid interfaces have intrinsic properties not included on the equations of the media in that cases mentioned before. This point of view was alread used in [9 13] to elaborate a formalism, based on Surface Green Functions, which establish an isomorphism between solids and fluidsrelatedtotheinterfacenormalmodes. Theaimofthis paper is to emphasie the fact that this wa to establish the matching conditions allows us to introduce the phsical characteristics related with the boundar conditions in a firmer floating. II. EQUATIONS FOR THE BULK IN A VISCOUS INCOMPRESSIBLE FLUID. To achieve a sstem of equations which describes the oscillator motion of a viscous incompressible fluid, the starting point is the linearied Navier-Stokes equation for a viscous fluid [14], which reads: ρ i = τ ij x j ρf i,ext where ρ is the densit of equilibrium, V i the components of the fluid particle, F i,ext are the external forces and τ ij is the stress tensor which, for a viscous incompressible fluid, has the following form [14] i τ ij = pδ ij η ) j 2) x j x i Here p are the small variations of pressure, η is the viscosit parameter and δ ij is the unit matrix elements. The subscripts i and j take the values and identicall. In the sum over repeated subscripts is understood. It had been assumed that the sstem is smmetric with respect to the x direction. Putting 2) in it is obtained ρ ρ [ η p2η ) [ η )] = 0 3) )] p2η ) = 0 4) where we have neglected the external forces. Also the continuit equation is needed, which expresses that the volume of an element of the incompressible fluid does not change during the motion. It has the form [14] = 0 5) The first step is to obtain the equation of motion and its solution of the each medium taken as infinite. So, the parameters as densit and viscosit are considered constants in the whole medium and it leads to transform equations 3) and 4) to: ρ ρ p η 2 V = 0 6) p η 2 V = 0 7) where 2 = 2 / 2 2 / 2 ). The solution of the sstem 5), 6) and 7) written as a vector field velocit, can be putted as the sum of an irrotational field related with the longitudinal mode) and a divergence free field related with the transverse modes) [1], i.e.: which satisf: V = V 1 V 2 8) V 1 = 0 9) V 2 = 0 10) An irrotational field is characteried b a scalar function, the potential function ϕ,, t), such that V 1 = ϕ 11) and the divergence free field can be described b a vector function [14], the stream or vorticit function ψ,, t), such that V 2 = ψ, ψ ) 12) The velocit components can thus be written in terms of the potential and the stream functions, as: V = ϕ ψ V = ϕ ψ 13) 14) Substitution of eqns. 13) and 14) in the continuit condition 5) gives 2 ϕ = 0 15) and substitution into eqns. 6) and 7) leads to { ρ ϕ } p { ρ ψ { ρ ϕ } p } η 2 ψ } { ρ ψ η 2 ψ = 0 16) = 0 17) Equations 16) and 17) are simultaneousl satisfied if one considers: 2

3 ρ ϕ p = C 1 18) ρ ψ η 2 ψ = C 2 19) The constants C 1 and C 2 are obtained from the condition at eroflow, and it givesraiseto C 1 = p o and C 2 = 0 where p o is the reference atmospheric pressure. The solution of equations 15), 18) and 19) are looking for as the following ϕ,,t) = Φ)e iκ ωt) 20) ψ,,t) = Ψ)e iκ ωt) 21) where the parameters κ and ω are the wavevector and frequenc of the wave respectivel. Substituting eqns. 20) and 21) in 15) and 19) gives the -dependence of the functions ϕ and ψ, which satisf: d 2 Φ) d 2 κ 2 Φ) = 0 22) d 2 Ψ) d 2 q 2 tψ) = 0 23) with q 2 t = κ2 iρω/η Equations 22) and 23) lead to the solutions: Φ) = C 1 e κ C 2 e κ 24) Ψ) = C 3 e q t C4 e q t 25) and in combination with eqns. 20) and 21) the give a solution of the form: ϕ,,t) = C 1 e κ C 2 e κ) e iκ ωt) 26) ψ,,t) = C 3 e q t C4 e q t ) e iκ ωt) 27) where C 1,C 2,C 3 and C 4 are constants to be determined b boundar and matching conditions. Then, the expressions for the varing velocit components and the pressure are finall obtained b substitution of eqns. 26) and 27) into 13), 14) and 18), respectivel. III. INTERFACE PROBLEM: MATCHING CONDITIONS. Now we will match the two media. We consider a surface which at rest coincides with the plane = 0 and it separates medium M 1 at < 0 from medium M 2 at > 0. Each one is viscous and incompressible. First of all, an solution has to fulfill the continuit of the velocit field across the surface according to V = V 2) at = 0 28) V = V 2) at = 0 29) Superscript 1 denotes medium M 1 and superscript 2 denotesmedium M 2. Therestofthematchingconditions at the interface are derived from the sstem of equations which govern the whole sstem. This is no usuall done in the normal program courses at the Universities and we consider that it is important to state that the matching conditions are, almost in all of the cases, contained in the equation of movement for the whole sstem taken as the composition of each one. These are eqns. 3) and 4) where the parameters η and ρ are constant, but taking different values on each medium. Integrating these equations through the surface about = 0 from ǫ to ǫ and later taking ǫ 0, it is obtained from eq. 3) [η 1 and from eq. 4) )] 2) = [η 2 = 0 = )] 2) 30) = 0 [ ] p 1 2η 1 = [ = 0 ] 2) = p 2 2η 2 p γ 31) = 0 wherep γ is the jump due to the surfacetension according with the Laplace Law [14]. The other terms in eqns. 3) and 4) vanish when ǫ 0 because the are continuous or have a finite jump at = 0. It is easil seen according to eq. 2) that eqns. 30) and 31) are the conditions of the continuit of the stress tensor components through the interface, as expected. On the other hand, the boundar conditions of the problem are regularit at ±, which leads to eliminate 4 of the 8 constants appearing in 26) and 27) for the two media. Using the resulting functions ϕ and ψ for each medium in conditions28)-31) the following sstem is met: iκ iκ q t1 q t2 C i i C 2iκ 2 η 1 2iκ 2 2 = 0 η 2 η 1 Q t1 η 2 Q t2 C 3 d 41 η 2 Q t2 d 43 2iκη 2 q t2 C 4 32) to determine the remaining constants, with the definitionsd 41 = η 1 Q t1 iυ, d 43 = 2iκη 1 q t1 Υ, Υ = γκ 3 /ω, Q t1 = 2κ 2 iρ 1 ω/η 1 and Q t2 = 2κ 2 iρ 2 ω/η 2. In deducing the expression in sstem 32) which comes from 31) it was considered that p γ = γ ) 3

4 on the surface, but as we are dealing with the velocit field, then it is obtained that: p γ which finall leads to: = γ 2 V 2 34) p γ = γκ2 iω κc 1 iκc 3 )e iκ ωt) 35) according to the substitution / iω and / κ 2 The vanishing of the determinant of 32) leads to the equation for the dispersion relation of the existing surface modes. It can be written as: ω 2[ ρ 1 ρ 2 )ρ 1 q t2 ρ 2 q t1 ) κρ 1 ρ 2 ) 2] γκ 3 [ρ 1 κ q t2 )ρ 2 κ q t1 )] 4κ 3 η 2 η 1 ) 2 κ q t1 )κ q t2 ) 4iκ 2 ωη 2 η 1 )ρ 1 κ ρ 2 κ ρ 1 q t2 ρ 2 q t1 ) = 0 36) This is the dispersion relation for capillar waves for the interface of two viscous incompressible fluids. In order to obtain the constants C 1,C 2,C 3 and C 4, an initial stimulus is needed according to an initial value problem [15] but this method, simple at the beginning, becomes rather complicated later and it is not good for a quite general stud. As was said on the introduction, the aim of this paper is not to get inside the dispersion relation of the normal modes at the interface of two viscous fluids at rest. For a better stud of this subject we recommend paper [16]. We get here to show a wa of solution also different from the usuall taken as onl the velocit vector as a function of an stream function. From now on, we will put our attention on the matching conditions and we will show that more that a mathematical information of the matching can be found on it, but also the phsics of the polariation can be deduced and how the interface moves in its oscillation. From eqns. 30) and 31) it can be obtained more information about the polariation of the modes on the interface. This will be done in the next section. IV. MATCHING CONDITIONS AND POLARIZATION. Therearetwopossiblemodesonthefluid: oneinwhich the fluid particle movesin the direction of the wavepropagation called longitudinal with notation LV,0) and another transverse to the direction of the wave propagation and normal to the interface with notation TN0,V). The longitudinal mode is related to the fluid compressibilit because this motion of the fluid particles is onl possible when its volume changes [14]. This analsis also holdswhenthemode isonthe interface, butthis doesnot mean that there are no longitudinal modes of oscillation on the interface when the fluids involved are incompressible. It had been shown in [1] that the interface, when oscillating, must be considered as a compressible one, because precisel its change in area is responsable for the increasing of its potential energ and therefore, for its oscillation. It is important to state that this is a fundamental argument to understand the movement of an interface in hdrodnamics. Nevertheless, now it can be shown that in spite of the compressibilit of the interface, there does not exist pure longitudinal mode if we are dealing with incompressible media. Let us demonstrate this. If a point o on the interface is considered moving with velocit, sa V S o in the -axis direction, then, according to the continuit of velocit, the point o, ǫ) in M 1 and the point o,ǫ) in M 2 must have the same velocit V S o if ǫ is small enough. As the interface is compressible, at the point 1 near enough o the velocit can be, for instance, V S 1 different ingeneralfromvs o but asthemedia are incompressible, at the point 1, ǫ) and 1,ǫ) the velocit must be Vo. S See Fig. 1. This is not in conformit with the continuit of the velocit through the surface and hence the pure longitudinal mode is not possible and onl the TN mode seems to be valid when the media are incompressible. After these considerations during the above demonstration, the student can keep the idea that the interface oscillations can onl occur in the -axis. This is the accurate moment to show to the student that things not alwas are as the apparentl seem to be, because that assumption does not take into account the different properties of each medium, whose response depends on its fundamental parameters, as densit and viscosit, which are different for each medium. Hence, it is evident that it must be analed, precisel, the interface matching conditions. It is useful to compare and to support the previous qualitative analsis with a quantitative and more profound one regarding the interface matching conditions. Recalling carefull eqns. 30) and 31) and supposing that such a wave propagates in direction with movement onl in direction TN mode), then V = V 2) = 0 and eq. 30) becomes η 1 = 0 = η 2) 2 = 0 37) ItisknownthatV iscontinuousalongtheinterfacefor all points. Then, the derivative with respect to is also the same in both hands of 37) and this expression onl holds if η 1 = η 2, i.e., if the media have the same viscosit. It does not mean for the interface to disappear because the densit of each medium can be different. Then, if the viscosities of the media have not the same value, the velocit component V along the interface must be different from ero to compensate the inequalit 37) and to 4

5 fulfill the continuit of the stress tensor in the direction given b eq. 30) ielding to a component of movement along direction. This mode will be called Sagittal mode or SV,V). The above analsis was done for the general case. Now we are able to take the particularcasein which one ofthe media is vacuum, for instance, M 2, with η 2 = 0. Then, condition 30) becomes [η 1 )] = 0 38) = 0 and it can be seen that also V must be non ero on the surface to hold eq. 38) with the corresponding Sagittal polariation movement. With this analsis on the conditions of stress component continuit in direction along the interface, it can be seen that if the two media are viscous at least one of them), the fluid particle of the interface movesin a Sagittal mode which combines movement in both directions: along the wave propagation in direction, and normal to the interface in direction. Then, it is qualitativel clear that the viscosit of each media plas a fundamental role in the coupling of modes evenforincompressiblefluids. Inspiteofthat, it couldbe amistaketosathatthe modesdecoupleifthe viscosities are equal. It should not be forgotten that eq. 31) is also important in the characteriation of the interface particle behaviour and it includes the pressure on each side of the surface. According to eq. 18), the pressure is associated with the longitudinal mode and the inertial effect of the fluid particle according with the densit of the media. Then, it contains the information of each components of the velocit and also of the densit and according to eq. 26) the pressure has a jump through the interface. This result, in combination with the analsis of eq. 30) make difficult to understand the role plaed b the densities of each media on the surface polariation movement, and it can not be reached from this onl analsis. This point is still a matter of investigation. V. CONCLUSIONS The present workis an attempt to givean example, using the hdrodnamics, of how the stud of the interface matching conditions allows us to make a plentiful and rich in details discussion. Moreover, of how the interface matching conditions content a sufficient information to conclude that the interface oscillation must be with a Sagittal mode and not with neither a pure longitudinal, nor a pure transversal one. This movement has been shown to be close related to the phsical properties of the media such as viscosities and that fact allows us to establish rigorousl that those are the parameter which characterie the interface movement and the response of each medium to an stimulus coming from the other one. It was seen that viscosit is the main parameter in the coupling of the two modes to achieve a Sagittal one, nevertheless within the framework of this formalism it is difficult to determine the role of viscosit and of the densit ratios in the coupling of modes. This aspect needs further investigation. [1] E. H. Lucassen-Renders and J. Lucassen, Advan. Colloid Interface Sci. 2, ). [2] J. Lucassen, Trans. Farada Soc. 64, ). [3] M. vandentempelandr.p. vanderiet, J. Chem. Phs. 42, ) [4] R. S. Hansen and J. A. Mann, J. Appl. Phs. 35, ) [5] C-S. Yih, J. Fluid Mech. 27, ) [6] A. P. Hooper and W. G. C. Bod, J. Fluid Mech. 128, ) [7] X. Li and R. S. Tankin, J. Fluid Mech. 226, ) [8] J. Cousin and C. Dumouchel, Atomiation and Spras 6, 1996) [9] García-Moliner, F., Ann. Phsique 2, ). [10] Velasco, V. R. and García-Moliner, F., Surface Sci. 67, ). [11] Velasco, V. R. and García-Moliner, F., Phsica Scripta 20, ). [12] Platero, G., Velasco, V. R. and García-Moliner, F., Phsica Scripta 23, ). [13] García-Moliner, F. and Velasco, V. R., Theor of Single and Multiple Interfaces, World Scientific, Singapore, 1992). [14] L. D. Landau and E. M. Lifshit, Fluid Mechanics Butterworth-Heinemann, 1987). [15] P. G. Drain and W. H. Reid, Hdrodnamics Stabilit Cambridge Universit Press, 1981). [16] A. Hernánde Rodrígue, J. Marín Antuña, H. Rodrígue Coppola and C. Dopao, Phs. Scripta 59, ). FIG. 1. Relation of velocities on the interface between two viscous fluids. 5

6 v o 2) ε v 1 2) v o S v 1 S o 1 v o -ε v 1