{'fraction_non_alphanumeric': 0.03710493046776232, 'fraction_numerical': 0.015518331226295827, 'mean_word_length': 4.216122650840751, 'pattern_counts': {'":': 0, '<': 7, '<?xml version=': 0, '>': 6, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 2, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
On the reflection of Alfv\'en waves and its implication for Earth's core modeling

Alfv\'en waves propagate in electrically conducting fluids in the presence of a magnetic field. Their reflection properties depend on the ratio between the kinematic viscosity and the magnetic diffusivity of the fluid, also known as the magnetic Prandtl number Pm. In the special case Pm=1, there is no reflection on an insulating, no-slip boundary, and the wave energy is entirely dissipated in the boundary layer. We investigate the consequences of this remarkable behaviour for the numerical modeling of torsional Alfv\'en waves (also known as torsional oscillations), which represent a special class of Alfv\'en waves, in rapidly rotating spherical shells. They consist of geostrophic motions and are thought to exist in the fluid cores of planets with internal magnetic field. In the geophysical limit Pm<<1, these waves are reflected at the core equator, where they are entirely absorbed for Pm=1. Our numerical calculations show that the reflection coefficient at the equator of these waves remains below 0.2 for Pm>0.3, which is the range of values for which geodynamo numerical models operate. As a result, geodynamo models with no-slip boundary conditions cannot exhibit torsional oscillation normal modes.

Alfvén waves propagate in electrically conducting fluids in the presence of a magnetic field. Their reflection properties depend on the ratio between the kinematic viscosity and the magnetic diffusivity of the fluid, also known as the magnetic Prandtl number P m. In the special case P m = 1, there is no reflection on an insulating, no-slip boundary, and the incoming wave energy is entirely dissipated in the boundary layer.
We investigate the consequences of this remarkable behaviour for the numerical modeling of torsional Alfvén waves (also known as torsional oscillations), which represent a special class of Alfvén waves, in rapidly rotating spherical shells. They consist of geostrophic motions and are thought to exist in the fluid cores of planets with internal magnetic field. In the geophysical limit P m 1, these waves are reflected at the core equator, but they are entirely absorbed for P m = 1. Our numerical calculations show that the reflection coefficient at the equator of these waves remains below 0.2 for P m ≥ 0.3, which is the range of values for which geodynamo numerical models operate. As a result, geodynamo models with no-slip boundary conditions cannot exhibit torsional oscillation normal modes.

I. INTRODUCTION
Hannes Alfvén first showed the theoretical existence, in an inviscid fluid of infinite electrical conductivity, of hydromagnetic waves that couple fluid motion and magnetic field [2]. The propagation of torsional Alfvén waves in the Earth's fluid core was thereafter predicted by Braginsky [4]. Such waves arise in rapidly rotating spheres or spherical shells in the presence of a magnetic field. In torsional Alfvén waves, the motions are geostrophic and consist in the rotation ω g (s) of nested cylinders centered on the rotation axis. They thus depend only on the distance s to the rotation axis. The period of the fundamental modes of torsional Alvén waves in the Earth's fluid core was first estimated to be about 60 years. This timescale was inferred from the analysis of the decadal length of day changes since the first half of the 19th century [17] and of the geomagnetic secular variation after 1900 [3]. With hindsight, these time series were not long enough to show convincingly variations with 60 years periodicity. Torsional waves with much shorter periods have now been extracted from time series of core surface flows for the time interval 1955-1985 [11]. If this discovery is confirmed, the period of the fundamental modes is of the order of 6 years and, as such, is much shorter than initially calculated.
Several authors have searched for torsional Alfvén waves in geodynamo simulations. Using stress-free boundary conditions, Dumberry & Bloxham [10] and Busse & Simitev [7] illustrated some parts of the torsional wave mechanism. Dumberry & Bloxham [10] found that the whole length of the geostrophic cylinders accelerates azimuthally as if they were rigid. The inertial forces, in their simulation, are however so influential that they dominate the Lorentz forces. Torsional Alfvén waves have finally been detected in a set of numerical simulations of the geodynamo with no-slip boundary conditions, for 0.5 ≤ P m ≤ 10, by Wicht & Christensen [27] (the magnetic Prandtl number P m is the ratio of kinematic viscosity over magnetic diffusivity). In both the geophysical (P m ∼ 10 −5 ) and the numerical studies, there seems to be no reflection of the torsional Alfvén waves upon their arrival at the equator. However, experimental studies in liquid metals have shown resonance effects on Alfvén normal modes [14] as well as reflection of wave packets [1].
In this paper we elaborate on the remark that reflection of Alfvén waves is controlled not only by the boundary condition, but also by the magnetic Prandtl number of the fluid in which they propagate [see 13, p. 23,24]. In the next section, we discuss the governing equations for one dimensional Alfvén waves and the associated boundary conditions for a solid and electrically insulating wall. We remark that for P m = 1 all the energy of the incident Alfvén wave is dissipated in a boundary layer, resulting in no reflected wave. In the following section, we change geometry to further emphasize our point and briefly present a direct numerical simulation of propagation and reflection of Alfvén wave in a non rotating spherical shell. That introduces the section devoted to the geophysical application, where we investigate torsional Alfvén waves in the Earth's core, modeled as a rapidly rotating spherical shell, calculating the energy loss on reflection at the Equator as a function of P m. Finally, we discuss the implications concerning the ability of geodynamo simulations to produce torsional eigenmodes and waves which are expected in the Earth's core.

II. REFLECTION OF ONE-DIMENSIONAL ALFVÉN WAVES
We introduce the problem through the example of Alfvén waves, transverse to a uniform magnetic field in an homogeneous and electrically conducting fluid, hitting a solid wall perpendicular to the imposed magnetic field [19]. The imposed uniform magnetic field B 0 is along the x-axis, while the induced magnetic field b(x, t) and the velocity field u(x, t) are transverse to this field, along y. Assuming invariance along y and z axes, the problem reduce to a 1-dimensional problem, u and b depending only on x. Projecting the Navier-Stokes equation and the induction equation on the y direction (on which the pressure gradient and the non-linear terms do not contribute), one obtains the following equations: where µ 0 is the magnetic permeability, ρ is the fluid density, ν the kinematic viscosity, and σ the electrical conductivity.

A. Elsasser variables
Introducing the two Elsasser variables h ± = u ± b/ √ µ 0 ρ, the equation of momentum (1) and the equation of magnetic induction (2) can be combined into where V A = B 0 / √ µ 0 ρ is the Alfvén wave speed, and η = (µ 0 σ) −1 is the magnetic diffusivity. It is already apparent that when ν = η, the right hand side of the previous equation vanishes, in which case h + and h − are fully decoupled. One can also show that h − travels in the direction of the imposed magnetic field, while h + travels in the opposite direction.
Introducing a length scale L and the time-scale L/V A , the previous equations take the following non-dimensional form: where the Lundquist number S and the magnetic Prandtl number P m are defined as: The propagation of Alfvén wave requires that the dissipation is small enough for the wave to propagate. This is ensured by S 1. The fact that (P m − 1)/(P m + 1) = −(P m −1 − 1)/(P m −1 + 1) establishes a fundamental symmetry of these equations: when changing P m into P m −1 , only the sign of the coupling term (right hand side of equations 4) changes.

B. Physical boundary conditions and reflection of Alfvén waves
These equations must be completed by boundary conditions. We assume that the wall is electrically insulating, and that the fluid velocity vanishes at the solid boundary (no-slip boundary condition), which translate to b = 0 and u = 0, leading to h ± = 0.
For P m = 1 the equations for h + and h − are fully decoupled, regardless of the value of S: In addition, for an insulating solid wall, the boundary condition h ± = 0 does not couple h + and h − either. As a result, reflection is not allowed at an insulating boundary when P m = 1, because reflection requires change of traveling direction, and thus transformation of h + into h − and vice versa. The energy carried by the wave has to be dissipated in the boundary layer. For P m = 1 the equations are coupled: for very small diffusivities (that is large Lundquist number S), the coupling will be effective only in a thin boundary layer. In addition the coupling will be more efficient as P m is further from 1. This gives a mechanism for reflection of Alfvén waves on an insulating boundary when P m = 1. Before giving a numerical illustration, it is instructive to consider the boundary conditions in the two limits P m = 0 and P m = ∞, with S 1 (dissipationless interior). In the limit P m = 0, there is no viscous term and the boundary condition, at the wall x = x 0 , reduces to There is perfect reflection. The incident (+) and reflected (-) waves have equal velocities and opposite magnetic fields. This also corresponds to a stress-free boundary condition for the velocity field in combination with an insulating wall (infinitely small vorticity sheet at the wall), leading to perfect reflection regardless of the value of P m used in equation 4. In this case the boundary condition for the velocity field is ∂ x u = 0, which translates into ∂ x (h + + h − ) = 0 and h + − h − = 0, effectively coupling h + and h − . In the limit P m = ∞, the boundary condition , at the wall x = x 0 , reduces instead to 1e-06 1e-05 0,0001 0,001 0,01 0,1 1 10 100 1e+03 Pm, magnetic Prandtl number The incident and reflected waves have opposite velocities and equal magnetic fields. This also corresponds to a no-slip boundary condition for the velocity field in combination with a perfectly conducting wall (infinitely small current sheet at the wall), leading to perfect reflection regardless of the value of P m used in equation 4. In this case the boundary condition for the magnetic field is ∂ x b = 0, which couples h + and h − . Another combination of boundary conditions inhibits reflection for P m = 1: for a stress-free (∂ x u = 0) and perfectly conducting wall (∂ x b = 0), which translates into ∂ x h + = 0 and ∂ x h − = 0, the fields h + and h − are decoupled, as for a no-slip insulating wall. Note finally that a wall with finite conductivity will allow some weak reflection, as illustrated by figure 4h.

C. Numerical simulations
We have performed a numerical simulation in a channel 0 ≤ x ≤ x 0 with a one-dimensional finite difference scheme. The Lundquist number is chosen large enough so that dissipation can be neglected in the interior. The boundary conditions were set to be electrically insulating and no-slip. The grid is refined next to the boundaries, in order to have at least 4 points in each boundary layer, which are Hartmann layers of thickness δ = √ νη/V A (see appendix A).
From the simulation of the traveling wave, we compute the transmission coefficient as the ratio of the velocity amplitude of the reflected and incident waves for different values of P m and S. The results are reported on figure 1.
As expected, there is full dissipation for P m = 1 and energy conservation for P m 1 or P m 1. Furthermore, the reflection coefficient R is independent of S, and exhibits the expected symmetry R(P m −1 ) = −R(P m). The

III. REFLECTION OF A LOCALIZED ALFVÉN WAVE PACKET ON A SPHERICAL BOUNDARY
The peculiar case where no reflection occurs is not specific to the planar, one-dimensional ideal experiment. Here, we run an axisymmetric simulation in a spherical shell permeated by a non-uniform magnetic field, without global rotation. The imposed magnetic field is the same as in Jault [16], and is represented by the dashed field lines of figure 2. Contrary to the simplest case of the previous section, it is a non-uniform magnetic field, which is not perpendicular to the boundaries. The observed behaviour of Alfvén wave packets hitting the curved boundaries should therefore apply to many systems.
The numerical pseudo-spectral code is the one used in Gillet et al. [12], but restrained to axisymmetry. It uses the SHTns library [23] for spherical harmonic expansion (Legendre polynomials) in the latitudinal direction, and second order finite differences in radius with many points concentrated near the boundaries. It time-steps both induction and momentum equation in the sphercial shell using a semi-implicit Crank-Nicholson scheme for the diffusive terms, while the coupling and (negligible) non-linear terms are handled by an Adams-Bashforth scheme (second order in time). The number of radial grid points is set to 500 and the maximum degree of Legendre polynomials to 120.
The Alfvén wave packets are generated mechanically by spinning the conducting inner core for a very short duration (compared to the Alfvén propagation time). Since the imposed magnetic field strength is not uniform, the wave front deforms as it propagates along the field lines. When the wave packet hits the outer insulating spherical shell, it does reflect and propagates back towards the inner shell for P m = 0.1 and P m = 10 but there is no reflection for P m = 1. This is illustrated by the snapshots of figure 2.

IV. REFLECTION OF TORSIONAL ALFVÉN WAVES
Finding evidence of propagation of Torsional Alfvén Waves (TAW) in the Earth's fluid core may open a window on the core interior. Properties of TAW in the Earth's core have thus been thoroughly investigated after the initial study of Braginsky [4]. They have been recently reviewed by Jault [15] and Roberts & Aurnou [20].

A. Model of Torsional Alfvén waves
In order to model TAW, magnetic diffusion and viscous dissipation are neglected in the interior of the fluid. The Earth's fluid core is modeled as a spherical shell of inner radius r i , outer radius r o and rotation rate Ω. Rapid rotation introduces an asymmetry between the velocity and magnetic fields and makes the velocity geostrophic, provided that λ ≡ V A /Ωr o 1 [16]. Note that the Lehnert number λ is about 10 −4 in the Earth's core. Geostrophic velocity in a spherical shell consists of the rotation ω g (s) of nested cylinders centered on the rotation axis. It thus depends only on the distance s from the rotation axis. A one-dimensional wave equation for the geostrophic velocity sω g (s) is obtained after elimination of the magnetic field b: with L = s 3 H(s) and H(s) the half-height of the geostrophic cylinders, andṼ 2 A involves only the z-average of the squared s-component of the imposed magnetic field. Braginsky [4] derived (8) rigorously in the geophysical case for which the viscous Ekman layer is thin compared to the magnetic diffusion layer located at the top and bottom rims of the geostrophic cylinders. This condition amounts to P mλ 1. Then, the velocity remains geostrophic in the magnetic diffusion layer. We have written the equation (8) in its simplest form, when the imposed magnetic field is axisymmetric, the mantle is insulating and Ekman friction at the rims of the geostrophic cylinders is neglected. The equation (8) needs to be completed by two boundary conditions, which can be derived when either P m 1 or P m 1. Interestingly, the equation (8) may be valid in the limit P m 1 but also in the limit P m 1 (provided P mλ 1). In the specific case P m 1, the appropriate boundary condition on the geostrophic velocity at the equator (on the inner edge of the Hartmann boundary layer) can be inferred from the boundary condition on the magnetic field. For an insulating outer sphere, it yields ∂ s ω g = 0 which corresponds to a stress-free boundary, as in the one dimensional wave case with P m → 0. In the case P m 1, the appropriate boundary condition is ω g = 0 as the angular velocity of the outermost geostrophic cylinder is immediately synchronized with the rotation of the solid outer sphere in the course of a spin-up experiment. This is equivalent to a no-slip boundary, as for the one dimensional wave case with P m → ∞.

B. Normal modes
Assuming that ω g varies with time as e ict , the equation (8) can be transformed into a normal mode equation: Transmission and reflection of TAW on the geostrophic cylinder tangent to the inner core set a special problem that we do not address here. As an intermediate step, we simply illustrate our discussion with results for the full sphere case, imposing ∂ s ω g | s=ε = 0, with ε 1 (we have checked the convergence of the numerical results as ε → 0). It is of interest to write the solution of this equation in the case c = 0 andṼ A uniform: A non-zero solution (uniform rotation ω g (s) = α 2 ) exists for the boundary condition ∂ s ω g | s=ro = 0 but not for the condition ω g | s=ro = 0 that applies when P m 1. We are interested in this latter case, despite its lack of geophysical realism, as contrasting the two boundary conditions sheds light on the nature of the constraint ∂ s ω g | s=ro = 0 that has always been used in TAW studies.
In the general case (c = 0, non-uniformṼ A ), it remains easy to calculate numerically a solution of (9) for 0 < s < r o . We have successfully checked our numerical results against the eigenvalues listed in the table C1 of Roberts & Aurnou [20], that have been obtained analytically for ∂ s ω g | s=ro = 0 andṼ A = 1. Then, the first eigenvalues are (0, 5.28, 8.63, 11.87, 15.07, ..), while in the caseṼ A = 1 and ω g | s=ro = 0 they are (2.94, 6.35, 9.58, 12.78, 15.95, ..). In the latter case, we recover our previous observation that 0 is not an eigenvalue.
In contrast with an often-made statement [5,15,20], the study of the case P m 1 shows that it is not required to have ∂ s ω g | s=ro = 0 to obtain solutions with bounded values of ω g for s ≤ r o . On the other hand, the singularity of ∂ s L at s = r o implies a singularity of ∂ s ω g (which is O((1 − s) −1/2 ) as s → 1) . That points to significant viscous dissipation once the viscous term is reintroduced.
When P m is neither very small nor very large, it is not possible to separate the interior region (where (9) applies) and the Hartmann boundary layer.
We can conclude the discussion of normal modes by noting that the solutions for the two cases P m 1 and P m 1 differ in a significant way at the equator. In both cases, solutions with bounded values of ω g in the interval

C. Numerical experiments
In order to determine the reflection coefficient of TAW at the equator of the outer shell, we use a set-up that resembles the Earth's core. The code is the same as the one described in section III, but this time with imposed global rotation. The total number of radial points is typically 1200 and the maximum degree of Legendre polynomials is set to 360.
For reflection to occur, there must be a non-zero imposed magnetic field B s at the equator. Hence we set the simplest potential quadrupolar field (generated from outside the sphere): B s = B 0 s, B z = 2B 0 z and B φ = 0. This ensures a local traveling speed V A (s) = B s (s)/ √ µ 0 ρ that is large near the reflection point (s = 1). The Lehnert number is small and always set to λ = V A /(Ωr o ) = 5 × 10 −4 , so that λP m is also small. The initial velocity field is along the azimutal direction φ and depends only on the cylindrical radius s: u φ (s) = sω g (s) = u 0 s exp(−(s − s 0 ) 2 / 2 ) with s 0 = 0.675. We used two different width = 0.02 and = 0.063. This initial velocity field splits into a torsional Alfvén wave packet propagating inwards that we do not consider here, and another traveling outwards that we carefully follow and we focus on the reflection of this wave packet at the equator of the outer shell (s = 1). The Lundquist number S based on the size of the spherical shell ranges from 6 × 10 2 to 8 × 10 4 and the Ekman number E = ν/Ωr 2 o and magnetic Ekman number Em = η/Ωr 2 o are both always very low and range from 5 × 10 −10 to 5 × 10 −7 over a wide range of magnetic Prandtl number: from P m = 10 −3 to P m = 10 2 .
We measure the extremum of the velocity field in the wave packet before and after the reflection, a i and a r respectively, at a fixed radius (s = 0.925 for = 0.02 and s = 0.75 for = 0.063), from which we compute the corresponding reflection coefficient R = a r /a i , reported in figure 3 for an insulating outer shell. We found no significant dependence with the Lundquist number S or the width of the initial pulse (R varies by less than 0.03).
As expected from the discussion of Alfvén waves equations, the combination P m = 1, no-slip boundary condition and insulating wall corresponds to a special case whereby no reflection at all occurs at the equator (see also fig. 4g).
However, there are differences with the planar case. First, the reflection coefficient is not symmetric with respect to P m = 1, as expected from our discussion of torsional eigenmodes in spherical geometry in the previous section. For large P m there is high dissipation and very little reflection compared to low P m. Second, the reflection coefficient is not as large.
Space-time diagrams of the reflection of the wave at the equator are presented in figure 4 for a few representative cases. The highest reflection coefficient occurs for the stress-free insulating case at P m = 1: from R = 0.86 at S = 1000 to R = 0.88 at S = 1.5 × very near the boundary, as the magnetic field must vanish, doing so by producing the reflected wave, just as in the planar case. This is not a boundary layer, but simply the superposition of the incident and reflected wave (see also appendix A). The Hartmann boundary layer is too small to be seen on these plots, but we checked that its size and relative amplitude for velocity and magnetic fields do match the analytic theory developed in appendix A.
For P m = 0.1, the reflected wave carries only 16% of the energy, the remaining being dissipated in the boundary layer. The magnetic field changes sign at the reflection, while the velocity keeps the same sign ( fig. 4cd). For P m = 10, the reflected energy drops to 3% and the small reflected velocity field has opposite sign, while the magnetic field (barely visible on figure 4) keeps the same sign ( fig. 4ef). During its propagation, the incoming wave is also much more damped than for P m = 0.1, even in the case where S or E have comparable values. This is due to strong dissipation at the top and bottom boundaries, which increases as the wave propagates toward the equator (visible in figure 4e) for P m > 1. This may not be unrelated to the previously discussed singularity for normal modes in the case P m > 1. A consequence of this large dissipation, is the difficulty to clearly identify the reflected wave, and to properly define a reflection coefficient. The values reported in figure 3 are thus not very precise for P m > 1.
It may also be worth noting that changing the magnetic boundary from insulating to a thin conducting shell allows weak reflection for P m = 1 and no-slip velocity ( fig. 4h), in agreement with the analysis of the governing equations (section II B).

D. Energy dissipation and normal modes
We want to emphasize that when no reflection occurs, the energy of the wave is dissipated very quickly. However, for liquid metals (P m 1), only a small amount of the wave energy is absorbed in the event of a reflection, but many successive reflections can lead to significant dissipation. Using the theoretical reflection coefficient, we can estimate the time-scale of dissipation of an Alfvén wave due to its reflections at the boundaries. In the case of an Alfvén wave turbulence (many wave packets) in a spherical shell of radius L with homogeneous mean energy e, permeated by a magnetic field of rms intensity B 0 , any wave packet will reach the outer insulating boundary once (on average) in the time interval L/V A . When it reflects on the boundary, it loses the fraction 1 − R 2 (P m) of its energy, where R(P m) is the reflection coefficient (in amplitude). We can then estimate the dissipation rate of energy e due to this process: Hence, the time-scale of dissipation at the boundaries which is inversely proportional to the strength of the magnetic field, and depends on the diffusivities only through P m. We can compare this to the dissipation of Alfvén waves of length scale in the bulk of the fluid: τ v = 2 2 /(η + ν). It appears that the length scale where surface and bulk dissipation are comparable is such that Replacing R by its theoretical expression A11 and assuming P m 1, we find for liquid metals Hence for the Earth's core with S ∼ 10 4 , and P m ∼ 10 −5 , the dissipation of Alfvén waves is dominated by the partial absorption at the boundaries for length scales larger than L/10. For numerical simulations of the geodynamo with S ∼ 10 3 ans R ∼ 0.2, we have L/ ∼ 30. These time-scales are also relevant for torsional normal modes. In one dimension, normal modes are a superposition of waves propagating in opposite directions. Hence, if the dissipation of waves is dominated by their reflection, so will it be for the normal modes. From the previous estimation of L/ in the Earth's core, we expect the dissipation of large wave-length torsional Alfvén waves (the ones that can be observed) to be dominated by the effect of reflection. Furthermore, in order to detect a normal mode, its dissipation time must be much larger than its period T = 2πL(cV A ) −1 . The pulsation c of the first torsional normal modes are given in section IV B in Alfvén frequency units, and their dissipation time can be estimated by τ s for the large-scale normal modes. We define a quality factor for torsional normal modes by Presence of normal modes requires Q 1. In the Earth's core, we find Q E 13 c and for numerical simulations of the geodynamo we find Q sim 0.16 c. Considering the largest modes (with c 5), the torsional oscillations should therefore be present in the Earth's core, but are completely absent even from the best current geodynamo simulations.

V. DISCUSSION: IMPLICATION FOR NUMERICAL GEODYNAMO MODELS AND THE EARTH-CORE
We showed that numerical simulations conducted for P m ∼ 1 cannot adequately reproduce the boundary conditions for torsional Alfvén waves in the Earth's core (where P m 1). The small reflection coefficient observed for TAW (figure 3) means that it is hard to observe TAW reflection at the equator in numerical simulations of the geodynamo which currently operate with 0.1 < P m < 10 [e.g. 22,26], where the waves are moreover mixed with thermal convection.
As for possible torsional eigenmodes, it is almost impossible to observe them with such low reflection coefficients. Unfortunately, that severely limits the ability of geodynamo simulations to exhibit torsional oscillation normal modes, because normal modes require a large reflection coefficient to be observable: their period (of order L/V A ) must be much larger than the energy dissipation time τ s (see expression 12). A few studies have tried to pin down torsional eigenmodes [10,21,27] but even though they report waves propagating with the appropriate speed, they report neither reflection of these waves, nor eigenmodes.
Another issue for geodynamo models with very low diffusivities, is that the part of the energy carried by Alfvén waves (regular or torsional) is dissipated very quickly (on an Alfvén wave time-scale), so that an Alfvén wave turbulence would be damped much faster, and the turbulent state may be far from what we would expect in the Earth's core.
Changing the boundary condition to stress-free simulates the case P m = 0 with a high reflection coefficient (R = 0.88), but still lower than the planar case. Even though this may still be problematic to observe eigenmodes, numerical models that use stress-free boundaries [e.g. 8,10,18,25] are intrinsically much more suited for the study of torsional normal modes. Quasi-geostrophic dynamo models that can compute dynamo models at very low magnetic Prandtl numbers [P m < 10 −2 in 24], could also provide an interesting tool to study torsional oscillations.
In the case of the Earth's core, a recent study [11] found no clear evidence for reflection at the equator, although this has yet to be confirmed. One might want to invoke turbulent viscosity (see the contrasted views of Deleplace & Cardin [9] and Buffett & Christensen [6] in a different context) to explain this fact, leading to an effective P m close to 1 and inhibiting reflection of torsional Alfvén waves. This would make numerical models more relevant, but is rather speculative. A solid conductive layer at the top of the core can also have a damping effect on the propagation and reflection of torsional waves, and we plan to investigate these matters in a forthcoming study.
where is the reciprocal Lundquist number based on the frequency: In the regime where Alfvén waves do propagate, we have 1 and also ω √ νη/V 2 A 1 so we can approximate the square root by its first order Taylor expansion, which leads to two solutions k 2 1 and k 2 2 : The solutions k = ±k 1 = ±ω/V A (1 − i ), correspond to the propagation in both directions of an Alfvén wave at the speed V A and with attenuation on a length scale V A /( ω). The solutions k = ±k 2 ±i/δ correspond to a Hartmann boundary layer of thickness δ ≡ √ νη/V A .
Finally, from equation A1 and A2 we know that U and B are related for each k by: and for the solutions k = ±k 1 and k = ±k 2 , it reduces to This means that for the travelling wave solution, U and B have always the same amplitude and the same phase when propagating in the direction opposite to the imposed magnetic field, or opposite phase when propagating in the same direction. For the boundary layers, in the limit P m 1 they involve the velocity field alone, whereas for P m 1 they involve only the magnetic field.

Reflection coefficient at an insulating wall
In order to derive the reflection coefficient, we consider an insulating wall at x = 0 with an incoming Afvén wave from the x > 0 region (k = +k 1 ), giving rise to a reflected wave (k = −k 1 ). The boundary conditions are matched by a boundary layer (k = +k 2 ) localized near x = 0 (the solution k = −k 2 is growing exponentially for x > 0 and has to be rejected for this problem). The solution to this problem reads where we have taken into account the fact that α −k1 = −α k1 (see eq. A8). The boundary conditions u = 0 and b = 0 at x = 0 lead to: from which we find the amplitude β of the velocity boundary layer contribution, and the reflection coefficient R of the amplitude of the velocity component: We are left to evaluate α k2 /α k1 using equations A8, which gives α k2 /α k1 = ν/η at leading order in , and thus which is independent of ω and V A . In the case P m = 1, we then have R = 0 and β = −1 which means that no reflection occurs and that the amplitude of the incoming wave is canceled by the boundary layer alone. It may be worth emphasizing that, although the boundary layer has the same thickness δ in the velocity and the magnetic field components, in the limit P m → 0, we have β → −2 and α k2 β → 0, so that the boundary layer is