arxiv: v1 [math.ap] 25 Jan 2018

Transcription

1 rspa.royalsocietypublishing.org Research Small-amplitude static periodic patterns at a fluid-ferrofluid interface M. D. Groves and J. Horn arxiv: v [math.ap] 5 Jan 08 Article submitted to journal Subject Areas: differential equations fluid mechanics Keywords: ferrofluids doubly periodic patterns bifurcation theory Author for correspondence: M. D. Groves groves@math.uni-sb.de Fachrichtung Mathematik Universität des Saarlandes Postfach Saarbrücken Germany Department of Mathematical Sciences Loughborough University Loughborough LE 3TU UK We establish the existence of static doubly periodic patterns in particular rolls rectangles and hexagons) on the free surface of a ferrofluid near onset of the Rosensweig instability assuming a general nonlinear) magnetisation law. A novel formulation of the ferrohydrostatic equations in terms of Dirichlet- Neumann operators for nonlinear elliptic boundaryvalue problems is presented. We demonstrate the analyticity of these operators in suitable function spaces and solve the ferrohydrostatic problem using an analytic version of Crandall-Rabinowitz local bifurcation theory. Criteria are derived for the bifurcations to be sub- super- or transcritical with respect to a dimensionless physical parameter.. Introduction Consider two static immiscible perfect fluids in the regions Ω := {x y z) : ηx z) < y < d} Ω : = {x y z) : d < y < ηx z)} separated by the free surface {y = ηx z)} where gravity acts in the negative y direction. The upper fluid is nonmagnetisable while the lower is a ferrofluid with a general nonlinear magnetisation law M = MH) = m H ) H H expressing the relationship between the magnetisation M of the ferrofluid and the strength of the magnetic field H. Subjecting the fluids to a vertically directed magnetic field of sufficient strength leads to the emergence of interfacial patterns see Figure ). c The Authors. Published by the Royal Society under the terms of the Creative Commons Attribution License by/4.0/ which permits unrestricted use provided the original author and source are credited.

2 {y = ηx z)} d Ω Figure. A rectangular periodic pattern at the interface {y = ηx z)}. In this article we present an existence theory for small amplitude doubly periodic patterns with ηx + l) = ηx) for every l L where x = x z) and L is the lattice given by Ω L = {ml + nl : m n Z} with l = l. We are especially interested in three patterns which are observed in experiments Figure ) namely rolls rectangles and hexagons see Figure 3). i) For rolls we seek functions that are independent of the z-direction and choose l = π ω 0) so that the periodic base cell is given by { x : x < ω π }. ii) For rectangles we choose l = π ω 0) l = 0 π ω ) so that the periodic base cell is given by { x z) : x z < ω π }. iii) For hexagons we choose l = π ω ) l 3 = π ω 0 3 ) so that we obtain an additional periodic direction l 3 = l + l = π ω 3 ) and the periodic base cell is given by { x z) : x < π ω x 3z < 4π ω x + 3z < 4π ω }. Notice that each of these patterns exhibits a rotational symmetry: the shape of the free surface is invariant under a rotation of the x z)-plane through respectively i) π ii) π 3 and iii) π. Figure. Experimental observation of static patterns at the fluid-ferrofluid interface: A) rolls; B) rectangles; C) hexagons Fachrichtung Physik Universität des Saarlandes). In Section a) we derive the mathematical formulation of this problem from physical principles. The governing equations equations.5).)) are formulated in terms of perturbations φ and φ of magnetic potentials corresponding to a uniform vertically directed magnetic field of strength hµh) in the upper fluid and h in the lower fluid where µ is obtained from the magnetisation law by the formula µs) = + ms)/s; the potentials are horizontally doubly periodic satisfying φ x + l y) = φ x y) φx + l y) = φx y) for every l L with a slight abuse of notation).

3 This problem was first studied by Cowley and Rosensweig []. Using a linear stability analysis they found that as the strength h of the magnetic field exceeds a critical value h c the flat surface destabilises and a hexagonal pattern of peaks appears. This phenomenon is known as the Rosensweig instability. A mathematically rigorous treatment of the problem was given by Twombly and Thomas [] who used coordinate transformations to flatten the free surface by transforming the a priori unknown domains Ω and Ω into fixed strips. Applying Lyapunov-Schmidt reduction reduces these transformed equations for rotationally symmetric patterns see below) to a locally equivalent one-dimensional equation which is solved using the implicit-function theorem; the result is the existence for values of h near h c of rolls and rectangles in addition to the hexagonal pattern. Twombly and Thomas s work is however flawed by some miscalculations and mathematical inconsistencies and is also restricted to linear magnetisation laws. In this article we present a more systematic approach which is motivated by the corresponding study of doubly periodic travelling water waves by Craig and Nicholls [3]; we also consider general nonlinear magnetisation laws. We work with dimensionless variables in terms of which the problem depends upon two dimensionless parameters β whose value is fixed) and γ see equation.)) and flatten the equations using Dirichlet-Neumann formalism. The Dirichlet-Neumann operator G for the upper fluid domain given by {ηx z) < y < } in dimensionless variables) is defined as follows. Fix Φ = Φ x z) solve the linear boundary-value problem φ xx + φ yy + φ zz = 0 η < y < φ = Φ y = η φ y = 0 y = 3 and define G η Φ ) = + ηx + ηz) φ n y=η = φ y η xφ x η zφ z) y=η..) The Dirichlet-Neumann operator G for the lower fluid domain { < y < ηx z)} is similarly defined as Gη Φ) = + η x + η z) µ gradφ + y) )φ n y=η = µ gradφ + y) )φ y η xφ x η zφ z) y=η.) where φ is the solution of the in general nonlinear) boundary-value problem divµ gradφ + y) ) gradφ + y)) = 0 < y < η.3) φ = Φ y = η.4) µ gradφ + y) )φ y + ) = µ) y =..5) The nonlinearity of.3).5) is inherited from that of the magnetisation law M = MH) for a linear magnetisation law the value of µ is constant and.3).5) are replaced by respectively Laplace s equation and a linear Neumann boundary condition). In Sections b) and c) we show that G and G are analytic functions of respectively η Φ ) and η Φ) in suitable function spaces and use these operators to recast the governing equations in terms of the variables Φ = φ y=η and Φ = φ y=η. The mathematical problem is thus to solve a system of equations of the form Gγ η Φ Φ)) = 0 where G : R X 0 Y 0 is given explicitly by the left-hand sides of equations.).4) and the function spaces X 0 Y 0 are specified in equation.5). Observe that this problem exhibits π rotational symmetry: it is invariant under rotations through respectively π and π 3 for rolls rectangles and hexagons and one may therefore replace X 0 and Y 0 by their subspaces of functions that are invariant under these rotations denoted by X sym and Y sym).

4 In Section 3 we discuss the existence of small-amplitude solutions to.6) within the framework of analytic Crandall-Rabinowitz local bifurcation theory see Buffoni and Toland [4 Chapter 8]) using γ as a bifurcation parameter. According to that theory values γ 0 of γ at which non-trivial solutions bifurcate from zero clearly Gγ 0) = 0 for all values of γ) necessarily have the property that the kernel of the linear operator L 0 := d G[γ 0 0] : X 0 Y 0 is non-trivial. We show that ker L 0 is non-trivial if and only if γ 0 = r k ) := µ µ ) µ k coth k + S k coth S ) k ) k for some k L \ {0} where µ = µ) µ = µ) and S = µ /µ + µ )) /. Choosing < µ µ ) /µ + ) and ω so that γ 0 = µ µ ) µ ω coth ω + S ω coth S ) ω ) ω is the unique maximum of the mapping k r k ) we find that ker L 0 = { v sink x) v cosk x) : k L with k = ω } where v R 3 is given by equation 3.5) see the discussion to Figure 4); this value γ 0 ) of β γ) corresponds to the Rosensweig instability. The dimension of ker L 0 is therefore determined by the number of vectors in L with length ω; for rolls rectangles and hexagons we find that dim ker L 0 is respectively 4 and 6 see Figure 5 and Sattinger [5 Section ] for a general discussion of this point). Because the kernel of L 0 is multidimensional one can not use Crandall-Rabinowitz local bifurcation theory directly. To overcome this problem we replace X 0 and Y 0 by X sym and Y sym thus restricting to solutions that are invariant under rotations through respectively π π and π 3 for rolls rectangles and hexagons. These restrictions ensure that dim ker L 0 = with ker L 0 = v 0 where v 0 = v e x z) and cos ωx rolls) e x z) = cos ωx + cos ωz rectangles) cos ωx + cos ω ) ) x + 3z + cos ω x 3z hexagons). 4 Verifying the remaining conditions in the analytic Crandall-Rabinowitz local bifurcation theorem yields the following result. Theorem.. The point γ 0 0) is a local bifurcation point for.6) that is there exist ε > 0 open neighbourhoods W sym of γ 0 0) in R X sym and V sym of 0 in X sym and analytic functions w : ε ε) V sym γ : ε ε) R with γ0) = γ 0 w0) = v 0 such that Gγs) sws)) = 0 for every s ε ε). Furthermore W sym N = {γs) sws)) : 0 < s < ε} where N = {γ v) R V sym \ {0}) : Gγ v) = 0}. In Section 4 we examine the bifurcating branches identified in Theorem.. Theorem.. Branches of small-amplitude doubly periodic solutions to the ferrohydrostatic problem bifurcate from the trivial solution at γ = γ 0. The bifurcation is i) transcritical in the case of hexagons ii) super- or subcritical in the case of rolls and rectangles depending upon the sign of a coefficient γ which is determined by µ and ω/. Explicit formulae for the coefficient γ are given in some special cases in Section 4 such formulae are unwieldy and it appears in general more appropriate to calculate them numerically

5 for a specific choice of µ). We note in particular that for constant µ corresponding to a linear magnetisation law) and very deep fluids rolls bifurcate subcritically for µ < µ c and supercritically for µ > µ c while rectangles bifurcate subcritically for µ < µ c and supercritically for µ > µ c where µ c = µ c = The same values were obtained by Silber and Knobloch [6] in a discussion of normal forms for this bifurcation problem and confirmed by Lloyd Gollwitzer Rehberg and Richter [7] as part of a wider numerical and experimental investigation. Finally we note that supercritical bifurcation of rolls is associated with supercritical) bifurcation of spatially localised patterns whose existence has been established by dynamicalsystems arguments by Groves Lloyd & Stylianou [8].. Mathematical formulation a) The physical problem We consider two static immiscible perfect fluids in the regions Ω := {x y z) : ηx z) < y < d} Ω := {x y z) : d < y < ηx z)} separated by the free interface {y = ηx z)}. The upper non-magnetisable fluid has unit relative permeability and density ρ while the lower is a ferrofluid with density ρ. The relations between the magnetic fields H H and the induction fields B B are given by the identities 5 µ 0 H = B µ 0 H + MH)) = B where µ 0 is the vacuum permeability and M is the magnetic intensity of the ferrofluid. Here and in the remainder of this paper equations for primed and non-primed quantities are supposed to hold in respectively Ω and Ω.) We suppose that MH) = m H ) H H where m is a nonnegative function so that in particular M and H are collinear. According to Maxwell s equations the magnetic and induction fields are respectively irrotational and solenoidal and introducing magnetic potential functions φ φ with grad φ = H grad φ = H one finds that these potentials satisfy the equations in which divgrad φ ) = 0 divµ grad φ ) grad φ) = 0.) µs) = + ms) ; s we assume that µ : 0 ) R is analytic and satisfies µ) + µ) > 0 so that the linearised version of the equation for φ is elliptic). Observe that φ is harmonic while φ satisfies a nonlinear elliptic partial differential equation; this nonlinearity is inherited from that of the magnetisation law M = MH) for a linear magnetisation law the value of µ is constant and the equation for φ reduces to Laplace s equation). At the interface we have the magnetic conditions where H t = H t H t = H t B n = B n t = ηx 0)T 0 ηz )T t = + η x + η z n = ηx ηz)t + η x + η z

6 are the tangent and normal vectors to the interface; it follows that φ φ y=η = 0 φ n µ grad φ )φ n y=η = 0..) The ferrohydrostatic Euler equations are given by grad p 0 + ρ gy) = 0 µ 0 MH) x z y) T) H gradp + ρgy) = 0 Rosensweig [9 Section 5.]) in which g is the acceleration due to gravity p 0 is the pressure in the upper fluid and p is a composite of the magnetostrictive and fluid-magnetic pressures. The calculation ) MH) x y z) T) H H = m H ) grad H = grad mt) dt 0 shows that these equations are equivalent to p 0 + ρ gy) = C 0 H µ 0 mt) dt p + ρgy) = C 0.3) 0 where C 0 C 0 are constants. The ferrohydrostatic boundary condition is given by p + µ 0 MH) n) = p 0 + κσ for y = ηx z) Rosensweig [9 Section 5.]) in which σ > 0 is the coefficient of surface tension and κ = ηxx + η z) + η zz + η x) η xη zη xz + η x + η z) 3/ is the mean curvature of the interface. Using.3) we find that 6 µ 0 MH) n) + µ 0 where C = C 0 C 0 or equivalently for y = ηx z) with H mt) dt + C + ρ ρ)gη σκ = 0 0 C + µ 0 + ηx + ηzφ yφ n µ grad φ )φ yφ n) ) ρ ρ )gη σκ µ 0 grad φ M grad φ ) = 0.4) s Ms) = tµt) dt. 0 The requirement that a uniform magnetic field and flat interface solves the physical problem that is η 0 φ 0 φ 0 ) = 0 µh)hy hy) is a solution to.).) and.4) leads us to choose C = µ 0 Mh) µ 0 µh)µh) )h /; we write η φ φ) = η φ φ) + η 0 φ 0 φ 0 ) so that η φ φ) = 0 0 0) is the trivial solution) and drop the tildes for notational simplicity. The next step is to introduce dimensionless variables and functions We find that ˆx ŷ ẑ) = µ 0h σ x y z) ˆη := µ 0h ˆµs) := µhs) σ η ˆφ := µ 0 h σ φ ˆMs) := Mhs). ˆφ := µ 0 h σ φ φ xx + φ yy + φ zz = 0 ηx z) < y <.5) divµ gradφ + y) ) gradφ + y)) = 0 < y < ηx z).6)

7 with boundary conditions and where φ y = 0 y =.7) µ gradφ + y) )φ y + ) µ) = 0 y =.8) φ φ + µ) ) η = 0 y = ηx z).9) φ + µ)y) n µ gradφ + y) )φ + y) n = 0 y = ηx z).0) ) M) µ) µ) γη σκ gradφ + µ)y) + M gradφ + y) ) + + ηx + ηzφ y + µ))φ + µ)y) n + ηx + ηzµ gradφ + y) )φ y + )φ + y) n = 0 y = ηx z).) α = ρ ρ )gd µ 0 h β = σ µ 0 h γ = αβ.) d and the hats have been dropped for notational simplicity. We seek periodic solutions to.5).) satisfying 7 ηx + l) = ηx) φ x + l y) = φ x y) φx + l y) = φx y) for every l L with a slight abuse of notation) where x = x z) and L is the lattice given by L = {ml + nl : m n Z} with l = l. Choose k k with k i l j = πδ ij for i j = and define the dual lattice L to L by L = {mk + nk : m n Z} so that our periodic functions can be written as ηx) = η k e ik x φ x y) = φ y)eik x k φx y) = φ k y)e ik x k L k L k L where η k = η k φ k = φ k φ k = φ k. We are especially interested in three periodic patterns namely rolls rectangles and hexagons see Figure 3). i) For rolls we seek functions that are independent of the z-direction and we choose l = π ω 0) so that the dual lattice L is generated by k = ω 0) and the periodic base cell is given by { x : x < π }. ω Furthermore the z-independent versions of equations.5).) are invariant under the reflection x x which corresponds to a rotation through π in the x z)-plane). ii) For rectangles we choose l = π ω 0) and l = 0 π ω ) so that the dual lattice L is generated by k = ω 0) and k = 0 ω) and the periodic base cell is given by { x z) : x z < π }. ω Furthermore equations.5).) are invariant under rotations through π/ in the x z)-plane.

8 iii) For hexagons we choose l = π ω 3 ) and l = π ω 0 3 ) so that we obtain an additional periodic direction l 3 = l + l = π ω 3 ). The dual lattice L is generated by k = ω 0) and k = ω 3 ) and the periodic base cell is given by { x z) : x < π x ω 4π x 3z < ω and } 4π + 3z <. ω Furthermore equations.5).) are invariant under rotations through π/3 in the x z)-plane. The mathematical problem is thus to solve.5).) for periodic functions η φ and φ in the domains Γ Ω per and Ω per where Ω per := {x y z) : x z) Γ } Ω Ω per := {x y z) : x z) Γ } Ω and Γ is the parallelogram defined by l and l or by l in the case of rolls). l π ω l π ω l π ω 4π 3ω l l 3 l 4π 3ω 4π 3ω 8 Figure 3. The lattice L and periodic base cell for rolls left) rectangles centre) and hexagons right). b) Dirichlet-Neumann formalism The Dirichlet-Neumann operator G for the upper fluid domain {ηx z) < y < β 0 } is defined as follows. Fix Φ = Φ x z) solve the linear boundary-value problem φ xx + φ yy + φ zz = 0 η < y < β 0.3) φ = Φ y = η.4) φ y = 0 y = β 0..5) and define G η Φ ) = + ηx + ηz) φ n y=η = φ y η xφ x η zφ z) y=η..6) The Dirichlet-Neumann operator G for the lower fluid domain { < y < ηx z)} is similarly defined as Gη Φ) = + η x + η z) µ gradφ + y) )φ n y=η = µ gradφ + y) )φ y η xφ x η zφ z) y=η.7) where φ is the solution of the in general nonlinear) boundary-value problem divµ gradφ + y) ) gradφ + y)) = 0 < y < η.8) φ = Φ y = η.9) µ gradφ + y) )φ y + ) = µ) y =..0)

9 It is also convenient to introduce auxiliary operators H and H given by H η Φ ) = φ y y=η Hη Φ) = φ y y=η.) where φ and φ are the solutions to the boundary-value problems.3).5) and.8).0). Using this Dirichlet-Neumann formalism we can recast the governing equations.5).) in terms of the variables η Φ = φ y=η and Φ = φ y=η as and ) η γη η Φ Φ + µ) )η = 0.) G η Φ ) + Gη Φ) + µ µ)) = 0.3) µ)g η Φ ) Gη Φ) µ µ)) + η ) H η Φ ) Φ y=η + M M)) µ Hη Φ) Hη Φ)Gη Φ) = 0.4) s in which = x z) T Ms) = tµt) dt 0 ) ) / µ = µ Φ + η Φ)Hη Φ) + + η )Hη Φ) + ) ) / M = M Φ + η Φ)Hη Φ) + + η )Hη Φ) +. 9 We study equations.).4) in the standard Sobolev spaces HperΓ r ) = ζ = ζ k e ik x : ζ k = ζ k ζ r < ζ r := CΓ ) k L k L where L is the dual lattice to L and their subspaces H r perγ ) = { ζ H r perγ ) : ζ 0 = 0 } + k ) r ζk consisting of functions with zero mean the value of the normalisation constant CΓ ) is π/ω for rolls 4π/ω) for rectangles and 8/ 3π/ω) for hexagons). In Section c) below we establish the following theorem. A function is analytic at the origin if is defined and analytic in a neighbourhood of the origin; in particular it has a convergent Maclaurin series.) Theorem.. Suppose that s > 5/. Formulae.6).7) and.) define mappings G G : HperΓ s ) Hper s / Γ ) Hper s 3/ Γ ) and H H : HperΓ s ) Hper s / Γ ) Hper s 3/ Γ ) which are analytic at the origin. Define X 0 := Hper s+/ Γ ) H perγ s ) HperΓ s ) Y 0 := HperΓ s s ) H per Γ ) Hper s 3/ Γ ).5) for s > 5/. Using Theorem. and the fact that H r Γ ) is a Banach algebra for r > we find that the left-hand sides of equations.).4) define a function G : R X 0 Y 0 which is analytic at the origin. A straightforward calculation shows that G η Φ ) + Gη Φ) + µ ) µ)) = 0 Γ and explains the choice of functions with zero mean in the second component of Y 0 ; using functions with zero mean in the second component of X 0 on the other hand ensures that the kernel of the linear operator d G[γ 0] does not contain any constant terms for any γ R). The mathematical problem is thus to solve

10 Gγ η Φ Φ)) = 0.6) for γ η Φ Φ)) R V 0 where V 0 is a neighbourhood of the origin in X 0 and Gγ 0) = 0 for all γ R. Observe that this problem exhibits rotational symmetry: it is invariant under rotations through respectively π π and π 3 for rolls rectangles and hexagons and one may therefore replace X 0 and Y 0 by their subspaces of functions that are invariant under these rotations denoted by X sym and Y sym). c) Analyticity of the Dirichlet-Neumann operators We study the boundary-value problems.8).0) and.3).5) by transforming them into equivalent problems in fixed domains cf. Nicholls and Reitich [0] and Twombly and Thomas []). The change of variable ỹ = y η + η ux ỹ z) = φx y z).7) { } transforms the variable domain Ω per into the fixed domain Σ = x y z) : x z) Γ y β 0 0) and the boundary-value problem.8).0) into div µ gradu + y) F η u) F 3 η u) F η u)) T ) ) = 0 < y < 0.8) u Φ = 0 y = 0.9) µ gradu + y).0 0) T F 3 η u) ) µ) = 0 y =.30) 0 where and F η u) = ηu x + + y)η xu y F η u) = ηu z + + y)η zu y F 3 η u) = ηu y + η + + y)η xu x + η zu z) + y) + η η x + η z)u y µ = µ grad u F η u) 0 F η u)) T) ) + 0 0) T + η we have again dropped the tildes for notational simplicity). Theorem.. Suppose that s > 5/. There exist open neighbourhoods V and U of the origin in respectively HperΓ s ) Hper s / Γ ) and HperΣ) s such that the boundary-value problem.8).30) has a unique solution u = uη Φ) in U for each η Φ) V. Furthermore uη Φ) depends analytically upon η and Φ. Proof. Write the left-hand sides of equations.8).30) as Hu η Φ) = 0 and observe that H : HperΣ) s HperΓ s ) Hper s / Γ ) is analytic at the origin with H0 0 0) = 0. Furthermore the calculation µ u xx + S uyy + uzz) d H[0 0 0]u) = u y=0 µ S uy y= where µ = µ) µ = µ) and S = µ /µ + µ )) / and standard existence and regularity theory for elliptic linear boundary-value problems show that d H[0 0 0] : HperΣ) s Hper s Hper s / Hper s 3/ Γ ) is an isomorphism. The stated result now follows from the analytic implicit-function theorem.

11 into The transformation.7) converts the formulae Gη Φ) = µ gradφ + y) )φ y η xφ x η zφ z) y=η Hη Φ) = φ y y=η Gη Φ) = µ u y F 3 η u)) u y y=0 Hη Φ) = + η.3) y=0 where u = uη Φ) is the unique solution to.8).30). Theorem.3. Hper s 3/ Γ ). Suppose that s > 5/. The formulae.3) define analytic functions G H : V To compute the Taylor-series representations of u and G we begin with the function ν : Hper s Σ)) 3 Hper s Σ) defined by νt ) = µ T + 0 0) T ). Observing that ν is analytic at the origin we write its Taylor series as νt ) = ν j {T } j) ).3) where ν j L j sh s per Σ)) 3 H s per Σ)) is given by ν j T... T j ) = j! dj ν[0]t... T j ) and may be computed explicitly from µ note in particular that ν 0 = µ ). The functions u n L n s HperΓ s ) Hper s / Γ ) HperΣ)) s with u 0 = 0) in the corresponding series uη Φ) = u n {η Φ)} n)).33) n=0 may be computed recursively by substituting the Ansatz.3).33) into equations.8).30). Consistently abbreviating m n {η Φ)} n)) to m n for notational simplicity one finds after a lengthy but straightforward calculation that divl grad u ) = 0 divl grad u n F n ) = 0 < y < 0 u Φ = 0 u n = 0 y = 0 u y = 0 L grad u n F n ) 0 0) T = 0 y = for n where 0 0 L = 0 S µ F n = µ F n F3 n F n ) T µ n η) j uy n j 0 0) T µ j= and n ν T j )grad u n j F n j F n j 3 F n j ) T ) R n 0 0) T n R j grad u n j F n j F n j 3 F n j ) T )

12 F n = ηu n x F3 n n = η η) j u n j y T n = + + y)η xu n y F n = ηu n z + + y)η zu n y + + y)η xu n x + y) ηx + ηz) n η) j u n j y n + η zu n z ) η) j grad u n j F n j 0 F n j ) T) R n = The Taylor-series representations of G and H are thus given by n=0 j n h +...+h j=n Gη Φ) = G n {η Φ)} n) ) Hη Φ) = H n {η Φ)} n) ) n=0 ν j T h... T hj ). where n n G n = µ I n + ν T j )I n j + R j I n j n H n = η) j u n j y y=0 y=0 and I n = n η) j u n j n y + η) j ηx + ηz)u n j y η xu n x + η zu n z ). For later use we record the formulae G = µ u y y=0 G = µ η) j u j y η xu x + η zu z) + µ u y) y=0 3 G 3 = µ η) j u 3 j y + ηx + ηz)u y η xu x + η zu z) + µ η) j u j y η xu x + η zu z) u y + µ u y) 3 + µ u x ) + u z) ) u y y=0 where µ = µ) and H = u y y=0 H = η) j u j y 3 H 3 = y=0 η) j u 3 j y y=0 for the first few terms in these series. The boundary-value problem.3).5) is handled in a similar fashion. The change of variable ỹ = y η η u x ỹ z) = φ x y z) { } transforms the variable domain Ω per into the fixed domain Σ = x y z) : x z) Γ y 0 β 0 ) and.3).5) into div grad u F η u ) F 3η u ) F η u )) T) = 0 0 < y <.34) u Φ = 0 y = 0.35) u y = 0 y =.36)

13 where F η u ) = ηu x + y)η xu y F η u ) = ηu z + y)η zu y 3 F 3η u ) = ηu y η + y)η xu x + η zu z) y) η η x + η z)u y we have again dropped the tildes for notational simplicity). The formulae G η Φ ) = φ y η xφ x η zφ z) y=η H η Φ ) = φ y y=η are converted into G η Φ ) = u y + F 3η u ) y=0 H η Φ u y ) = η y=0.37) where u = u η Φ ) is the unique solution to.34).36). Theorem.4. Suppose that s > 5/. i) There exist open neighbourhoods V and U of the origin in respectively HperΓ s ) Hper s / Γ ) and HperΣ s ) such that the boundary-value problem.34).36) has a unique solution u = u η Φ ) in U for each η Φ ) V. Furthermore u η Φ ) depends analytically upon η and Φ. ii) The formulae.37) define analytic functions G H : HperΓ s ) Hper s / Γ ) Hper s 3/ Γ ). The functions u n L n s HperΓ s ) Hper s / Γ ) HperΣ s )) with u 0 = 0) in the Taylor series u η Φ ) = u n {η Φ )} n)) n= are computed recursively by substituting this Ansatz into equations.34).36); one finds that divgrad u ) = 0 divgrad u n F n F n 3 F n ) T ) = 0 0 < y < for n where u Φ = 0 u n = 0 y = 0 u y = 0 u n y = 0 y = F n = ηu n x + y)η xu n y F n = ηu n z + y)η zu n y F 3 n n = η η) j u n j y + y)η xu n x y) ηx + ηz) n η) j u y n j. The Taylor-series representations of G and H are given by with G η Φ ) = G n {η Φ )} n) ) H η Φ ) = n= G n n = η) j u n j y H n n = η) j u n j y y=0 n 3 η) j η x + η z)u n j y n= + η zu n z ) H n {η Φ )} n) ) + η xu n x + η zu n z ) y=0

14 and in particular we find that and G = u y y=0 G 3 = G = η) j u j y + η xu x + η zu z ) y=0 η) j u y 3 j ηx + ηz)u y + η xu x + η zu z ). y=0 H = u y y=0 H = η) j u j y H 3 = η) j u 3 j y. y=0 y=0 3. Existence theory Next we introduce the Crandall-Rabinowitz theorem cf. Buffoni and Toland [4 Theorem 8.3.]) an application of which yields a local bifurcation point of the equation Gγ η Φ Φ)) = 0 3.) where the components of G : R V 0 Y 0 are given by the left-hand sides of.).4). Theorem 3. Crandall-Rabinowitz theorem). Let X and Y be Banach spaces V be an open neighbourhood of the origin in X and F : R V Y be an analytic function with Fλ v) = 0 for all λ R. Suppose also that i) L := d F[λ 0 0] : X Y is a Fredholm operator of index zero ii) ker L = v 0 for some v 0 X iii) the transversality condition P d d F[λ 0 0] v 0 )) 0 holds where P : Y Y is a projection with Im L = ker P. The point λ 0 0) is a local bifurcation point that is there exist ε > 0 an open neighbourhood W of λ 0 0) in R X and analytic functions w : ε ε) V λ : ε ε) R with λ0) = λ 0 w0) = v 0 such that Fλs) sws)) = 0 for every s ε ε). Furthermore 4 W N = {λs) sws)) : 0 < s < ε} where N = {λ v) R V \ {0}) : Fλ v) = 0}. The first step is to determine the maximal positive value γ 0 of the parameter γ for which the kernel of the linear operator L 0 := d G[γ 0 0] : X 0 Y 0 which is given by the explicit formula η Φ Φ + µ )η L 0 Φ = G η Φ ) + µg η Φ) + µ H η Φ) 3.) Φ η xx + η zz γ 0 η µ G η Φ ) + G η Φ) + µ H η Φ)) with G η Φ ) = k tanh k Φ eik x G k η Φ) = µ S k tanh S k Φ β k e ik x 0 k L k L and G η Φ) = µ H η Φ) is non-trivial. Writing v X 0 as vx) = k L v k e ik x with v k = η k Φ k Φ k) T and v k = v k we find that L 0 v = 3.3) k L L 0 k )v k e ik x 3.4)

15 where for k 0 and µ L 0 k ) = 0 k tanh k k γ 0 µ k tanh k µ L 0 0) = γ µ S k tanh S k µ S k tanh S k where we have identified the subspace {η 0 Φ 0 Φ 0) T : Φ 0 = 0} of R3 with R ). From this observation it follows that ker L 0 is non-trivial if that is det L 0 k ) = µ µ ) S k tanh k tanh S k γ 0 = r k ) := k + γ 0 ) µ k tanh k + S k tanh S k ) = 0 µ µ ) µ k coth k + S k coth S ) k ) k 5 for some k L \ {0}. The function k r k ) which satisfies r0) = 0 and r k ) as k takes only negative values for > µ µ ) /µ + ) while for < µ µ ) /µ + ) it has a unique maximum ω with rω) > 0 see Figure 4); we choose γ 0 = rω) and note the relationships = µ µ ) ω h ω) ω ḣ ω) h ω) where ω = ω/ and h ω) = µ coth ω + S coth S ω. ) µ µ γ 0 = ) ) ω ωh ω) r k ) r k ) γ 0 k ω k Figure 4. The graph of the function k r k ) for > µ µ ) /µ +) left) and < µ µ ) /µ +) right). Noting that ker L 0 ω) = v where v = µ we find that see Figure 5). ) µ S tanh S ω coth ω + µ S tanh S ω coth ω ker L 0 = { v sink x) v cosk x) : k L with k = ω } 3.5)

16 i) For rolls the dual lattice L is generated by k = ω 0) so that k k = ω and hence dim ker L 0 =. ii) For rectangles the dual lattice L is generated by k = ω 0) and k = 0 ω) so that k k k k = ω and hence dim ker L 0 = 4. iii) For hexagons the dual lattice L is generated by k = ω 0) and k = ω 3 ) so that k k k k k 3 k 3 = ω where k 3 = k k and hence dim ker L 0 = 6. k k 3 k = ω k = ω k = ω k k k k k k k k k 3 Figure 5. The vectors generated by L with length ω in the case of rolls left) rectangles centre) and hexagons right). Define the projection P ω : X 0 X 0 Y 0 Y 0 by P ωvx) = where v is given by formula 3.3) so that k L k =ω v k e ik x k 6 X 0 = X ω X r Y 0 = Y ω Y r with X ω = P ω[x 0 ] X r = I P ω)[x 0 ] Y ω = P ω[y 0 ] Y r = I P ω)[y 0 ]. Using 3.4) one finds that Im L 0 Xω Y ω and Im L 0 Xr Y r and we prove that L 0 : X 0 Y 0 is a Fredholm operator of index zero in two steps. Lemma 3.. The mapping L 0 is an isomorphism X r Y r. Proof. The mapping L Xr 0 is formally invertible on Y r with L 0 χ Ψ Ψ) T = k L k =ω L 0 k ) χ k Ψ Ψ k k) T e ik x 3.6) where L 0 k ) = µ S k tanh S k det L 0 k ) µ ) k tanh k γ 0 + k ) µ ) µ ) k tanh k det L 0 k ) 0 µ γ 0 + k µ µ ) µ γ 0 + k det L 0 k )

17 for k L with k > ω and χ 0 L 0 0) 0 = Ψ 0 γ 0 Ψ 0 0 χ 0 µ )γ 0 Ψ 0 Denoting the right-hand side of equation 3.6) by η Φ Φ) T and using the estimates for k > ω one finds that η s+/ χ 0 + Ψ 0 ) + S k tanh S k k tanh k χ 0 + Ψ 0 + k L k =ω0 k L k =ω0 k s+ χ s + Ψ s + Ψ s 3/ ; k tanh k ) det L 0 k )) k tanh k k det L 0 k ) k s 3 k χ k + Ψ k + Ψ k ) k tanh k ) χ β k + Ψ + Ψ k k 0 7 similar calculations yield Φ s χ s + Ψ s + Ψ s 3/ Φ s χ s + Ψ s + Ψ s 3/. We conclude that L 0 : Y r X r exists and is continuous. Corollary 3.. The operator L 0 : X 0 Y 0 is a Fredholm operator of index zero. Proof. A straightforward calculation shows that v / Im L 0 ω) so that kerl 0 ω)) = ker L 0 ω)) and hence Y ω = Im L0 ω) ker L 0 ω) ) sink x) Im L 0 ω) ker L 0 ω) ) cos x)) k L k =ω = Im L 0 Xω ker L 0. Using this decomposition and Lemma 3. we find that ) Y 0 = Y ω Y r = Im L Xω 0 ker L 0 Im L Ur 0 = Im L 0 ker L 0. It follows that Im L 0 is closed and codim Im L 0 = dim ker L 0 so that L 0 : X 0 Y 0 is a Fredholm operator of index zero. Because the kernel of L 0 is multidimensional we can not use Theorem 3. directly. To overcome this problem we recall that G and hence L 0 ) is invariant under certain rotations see below) and seek solutions to 3.) in X 0 that have this rotational symmetry denoting the relevant subspaces of H r perγ ) H r perγ ) X 0 and Y 0 by H r symγ ) H r symγ ) X sym and Y sym so that X sym = Hsym s+/ Γ ) H symγ s ) HsymΓ s ) Y sym = HsymΓ s s ) H sym Γ ) Hsym s 3/ Γ ) for s > 5/. Note that X sym and Y sym are invariant under P ω so that according to the above analysis L 0 : X sym Y sym is a Fredholm operator of index zero and Y sym = Im L 0 ker L 0.

18 i) For rolls we consider functions that are independent of the z-coordinate and lie in the subspace H r symγ ) = { ζ H r perγ ) : ζ x) = ζ x) } of functions which are invariant under rotations through π. ii) For rectangles we work with the subspace H r symγ ) = { ζ H r perγ ) : ζ x z) = ζ z x) } of functions which are invariant under rotations through π/. iii) For hexagons we work with the subspace { HsymΓ r ) = ζ HperΓ r ) : ζ x z) = ζ x ) 3z ) )} 3x + z of functions which are invariant under rotations through π/3. These restrictions ensure that dim ker L 0 = with ker L 0 = v 0 where v 0 = ve x z) with cos ωx rolls) e x z) = cos ωx + cos ωz rectangles) cos ωx + cos ω ) ) x + 3z + cos ω x 3z hexagons). The projection P : Y sym Y sym onto ker L 0 along Im L 0 is given by 8 P χ Ψ Ψ) T = C χ ω0) Ψ ω0) Ψ ω0) )T v )v e x z) 3.7) where γ 0 + ω µ v = γ 0 + ω ) µ ) S ω coth S ω + coth ω C = µ ωs tanh S ω µ S tanh S ω coth ω + S coth S ω ) tanh ω µ coth ω + S coth S ω ) + v R 3 solves the equation L 0 ω) T v = 0 and C = v v ) ). Lemma 3.. The transversality condition P d d G[γ 0 0] v 0 )) 0 is satisfied. Proof. It follows from the calculation d d G[γ 0 0] η Φ Φ) T ) = 0 0 η) T and the formula 3.7) for P that P d d G[γ 0 0] v 0 )) = C µ ) S tanh S ωµ coth ω + S coth S ω)v e x z). The facts established above confirm that the hypotheses of Theorem 3. are satisfied an application of which yields Theorem.. 4. The bifurcating solution branches In this section we examine the bifurcating solution branches identified in Theorem. by applying the following supplement to the Crandall-Rabinowitz theorem.

19 Theorem 4.. Suppose that the hypotheses of Theorem 3. hold. In the notation of that theorem let Q : X X be a projection with Im Q = ker L and the Taylor series of the functions w : ε ε) V λ : ε ε) R be given by λs) = λ 0 + sλ + s λ +... ws) = v 0 + sw +... where λ λ... R and w w... ker Q. i) The coefficient λ satisfies the equation ) P! d F[λ 0 0]v 0 v 0 ) + λ P d d F[λ 0 0] v 0 )) = 0 and the bifurcation is transcritical if λ is non-zero see Figure 6). s λ 0 λs) Figure 6. Transcritical Crandall-Rabinowitz bifurcation for λ < 0 left) and λ > 0 right). s λ 0 λs) 9 ii) Suppose that λ is zero. The coefficient λ satisfies the equation P d F[λ 0 0]v 0 w ) + ) 3! d3 F[λ 0 0]v 0 v 0 v 0 ) + λ P d d F[λ 0 0] v 0 )) = 0 where w ker Q solves the equation d F[λ 0 0]w ) =! d F[λ 0 0]v 0 v 0 ). The bifurcation is supercritical for λ > 0 and subcritical for λ < 0 see Figure 7). s s λ 0 λs) λ 0 λs) Figure 7. Crandall-Rabinowitz bifurcation for λ = 0 λ > 0 supercritical left) and λ = 0 λ < 0 subcriticalright). To apply this theorem we write the Taylor series of the analytic functions w : ε ε) V sym γ : ε ε) R given in Theorem. as γs) = γ 0 + sγ + s γ +... ws) = v 0 + sw +... with γ γ... R and w w... keri P ) and introduce the operators L = d d G[γ 0 0] ) Q 0 =! d G[γ 0 0] C 0 = 3! d3 G[γ 0 0].

20 One finds that L v) = 0 0 η) T and Q 0 v v) = 0 G + G + µ H + µ H + µ Φ ) Φ + µ µ )H + µ Φ G H µ H µ Φ ) µ G G µ H G C 0 v v v) 0 G 3 + G 3 + µ H 3 + µ H H + µ µ ) Φ H µ η Φ H µ H 3 = µ G 3 G 3 µ H 3 µ H H µ µ ) Φ H + µ η Φ H... 6 µ H 3 +G G η Φ ) H G + µ η Φ) + µ µ )H H + 3 µ µ )H 3 G H + ηxη zz + ηzη xx η xη zη xz 3 η η where v = η Φ Φ) T and... µ =... µ). Theorem 4. shows that γ = [Q 0v 0 v 0 )] v [L v 0 ] v 0 where [ζ] = ζ ω0) = ζe Γ with componentwise extension) and γ = [Q 0v 0 w ) + C 0 v 0 v 0 v 0 )] v [L v 0 ] v for γ = 0 where w keri P ) solves the equation L 0 w = Q 0 v 0 v 0 ). A straightforward calculation shows Q 0 v 0 v 0 ) can be written as a sum in which each summand is a constant vector multiplied by either e or e. For hexagons we find that γ generally does not vanish while for rolls and for rectangles [ ] [ e = + cos ωx [ ] e = [ ] = 0 [ e ] [ = ] cos ωx = 0 + cos ωx + z) + cos ωx z) + cos ωx + cos ωz) [ e ] [ = ] cos ωx + cos ωz) = 0 ] = 0 so that in both cases γ = 0. Attempting to compute explicit general expressions for γ leads to unwieldy formulae it appears more appropriate to calculate them numerically for a specific choice of µ that is a specific magnetisation law). Here we confine ourselves to stating the values of the coefficients for two particular special cases.

21 i) Constant relative permeability µ corresponding to a linear magnetisation law): We find that γ = µ µ µ + ω µc 8µ + )ω 4ω ) + γ 0 4µ ) ω 3 t t ) t for rolls and γ = µ µ µ + ω µc + ω 3 µ ) 4 + t ) t t ) 4 + t ) )) µω µ ) 3 4γ 0 µ + t ) µ + ) ω 3 ω t µ 8γ 0 + ω ) 3 ) + t t 8µ + )ω ω + γ 0 t µ ) ω 3 ) t t ) ω3 µ ) 4 t 4 4t )) t t ) µc 8µ + )ω 4ω ) + γ 0 4µ ) ω 3 t t ) t + ω 3 µ ) 4 + t ) t t ) 4 + t ) )) ωµ + ) µ ) t 5ω 4γ 0 + ω ) 9+ t t +4 ) t t ) ω µµ ) 3 γ 0 µ + ) t ) for rectangles where t = tanh ω t = tanh ω t = tanh ω ω = ω/ and ω + γ0 µ ) t ) γ 0 ω )ω t C = µ γ 0 + ω ) t ) γ 0 4ω )ω t C = The sign of γ clearly depends upon µ and ω see Figure 8). ω 0 ω µ µ Figure 8. The sign of the coefficient γ as a function of µ and ω for a linear magnetisation law for rolls left) and rectangles right). The shaded and white areas show the regions in which the bifurcation is respectively super- and subcritical.

22 ii) Small values of corresponding to deep fluids): Abbreviating µ) µ) µ) µ)... to respectively µ µ µ µ... one finds that γ = t µ 3... µ + µ... µ µ µ µ µ µ µ µ µ + µ µ + 4µ3 µ + 49µ µ µ 8µ µ + 9µ µ µ 0µ µ + 5µ µ µ + 3 µ µ + 8 µµ3 + 6µ µ + 8µ µ3 6µ µ 8 ) µ3 + 6µ 4 µ + 3µ3 µ µ + 6µ µ µ + 30µ µ µ + 5µ µ + 0µ µ µ 6µ µ 46µ µ µ ) 5µ µ 0 µ ) 7 µ µ µ 5t 6 µ + µ ) + t 96µ µ5 µ + 9µ4 µ + 850µ3 µ3 + 7µ µ4 5µ5 70µ4 µ 80µ3 µ + o) 740µ µ3 50µ µ4 + 96µ4 + 0 µµ3 + 73µ µ 4µ µ3 7 ) µ4 + 88µ µ7 µ + 880µ6 µ + 880µ5 µ µ4 µ4 + 88µ3 µ5 56µ7 98µ6 µ 3µ5 µ 864µ 4 µ3 4µ3 µ4 80µ6 46µ5 µ 64µ4 µ 40µ3 µ3 0µ µ4 +56µ5 +67µ4 µ ) 8µ µ3 µ +04µ µ3 0µ µ4 +7 µµ3 +94µ µ +84µ µ3 +0 µ ) 6 µ µ4 04t 6 µ + µ ) as 0 for rolls and γ = t )µ µ µ 4 60)µ... µ µ )µ µ ) µ µ µ ) µ µ µ ) µ µ µ ) µ )µ µ )µ µ ) µ µ ) µ µ ) µ µ +4 60)µ... µ 4 89)µ µ 68 40)µ µ ) µ µ )µ µ 4 60)µ 3... ) µ )µ 4 µ ) µ µ 3 µ )µ µ µ )µ µ )µ µ µ 44 00)µ µ µ 445 6)µ µ )µ µ )µ µ µ +44 ) 00) µ µ ) 7 µ µ 50 7 )t 6 µ + µ ) + t ) µ µ )µ µ )µ µ )µ 5 µ )µ 4 µ )µ 3 µ )µ µ )µ 4 µ )µ 3 µ )µ µ )µ µ )µ )µ ) µ )µ 4 ) )µ 3 µ )µ 7 µ )µ 6 µ )µ 5 µ )µ 4 µ )µ 6 µ )µ 5 µ )µ 4 µ )µ 3 µ ) µ µ )µ µ )µ µ )µ )µ )µ 5 µ )µ 4 µ )µ 3 µ )µ µ )µ 4 µ )µ 3 µ )µ µ )µ µ )µ )µ ) µ ) 608)µ 4 µ ) 6 µ )t 6 µ + µ + o) ) as 0 for rectangles where t = µ µ + µ ) +. Note that ω = µ µ ) µ + S ) + o) γ 0 = µ µ ) 4 4µ + S ) + o)

23 as 0.) We note in particular that for constant µ corresponding to a linear magnetisation law) rolls bifurcate subcritically for µ < µ c and supercritically for µ > µ c while rectangles bifurcate subcritically for µ < µ c and supercritically for µ > µ c where µ c = µ c = Figure 9 shows the sign of γ for the Langevin magnetisation law µs) = + M cothγs) ) s γs in the limit 0 where M and χ 0 are respectively the magnetic saturation and initial susceptibility of the ferrofluid and γ = 3χ 0 /M. M M ) γ 4 6 γ Figure 9. The sign of the coefficient γ as a function of M and γ for the Langevin magnetisation law 4.) for rolls left) and rectangles right) in a ferrofluid of great depth. The shaded and white areas show the regions in which the bifurcation is respectively super- and subcritical. Data Accessibility. This paper has no additional data. Authors Contributions. M.D.G. and J.H. collaborated equally on each aspect of this work. Competing Interests. We declare we have no competing interests. Funding. No external funding was involved in this work. References. Cowley MD Rosensweig RE. 967 The interfacial stability of a ferromagnetic fluid. J. Fluid Mech Twombly EE Thomas JW. 983 Bifurcating instability of the free surface of a ferrofluid. SIAM J. Math. Anal Craig W Nicholls DP. 000 Traveling two and three dimensional capillary gravity water waves. SIAM J. Math. Anal Buffoni B Toland JF. 003 Analytic Theory of Global Bifurcation. Princeton N. J.: Princeton University Press. 5. Sattinger DH. 978 Group representation theory bifurcation theory and pattern formation. J. Func. Anal Silber M Knobloch E. 988 Pattern selection in ferrofluids. Physica D Lloyd DJB Gollwitzer C Rehberg I Richter R. 05 Homoclinc snaking near the surface instability of a polarisable fluid. J. Fluid Mech Groves MD Lloyd DJB Stylianou A. 07 Pattern formation on the free surface of a ferrofluid: spatial dynamics and homoclinic bifurcation. Physica D Rosensweig RE. 997 Ferrohydrodynamics. New York: Dover. 0. Nicholls DP Reitich F. 00 A new approach to analyticity of Dirichlet-Neumann operators. Proc. Roy. Soc. Edin. A