True, with . The hyperbolic cosine identity and the initial values give by mathematical induction, using the same Chebyshev polynomial recurrence.
Use a power-series solution of a differential equation, . Comparing coefficients gives
In particular , while and are free. The two prescribed derivative pairs select and . The next coefficients are , , , . Hence
These are the first three nonzero terms when the respective parameter is nonzero; if or , the corresponding solution vanishes identically.
For closed forms, set . Cancellation of the first-derivative terms reduces the equation, for , to . The analytic solutions extend across zero, and their hyperbolic cosine and hyperbolic sine expansions identify
This is the hyperbolic reduction of a regular-singular differential equation: although the original leading coefficient vanishes at zero, both selected power-series solutions of a differential equation are entire.
A positive-eigenvalue normal mode satisfies . Solving the constant-coefficient equation gives the sine, cosine, hyperbolic sine and hyperbolic cosine family. To avoid confusing modal coefficients with the filament bending modulus, write those coefficients as .
The free conditions at zero require and . With , the remaining conditions are
Its determinant is . Therefore the free-free bending spectrum is
A convenient corresponding eigenfunction is
where gives unit norm. The stable characteristic equation for free-free bending modes is . Intersections or numerical bracketing give
For the entire positive sequence, let . Expanding near its zero and using gives
There are also two zero-energy filament modes, not captured by substituting into the finite-coefficient trigonometric expression. Solving with free endpoints gives the affine kernel. An orthonormal basis of that kernel is
The positive normal modes are orthogonal to both. Including these two rigid-motion modes is essential for completeness; the regular self-adjoint operator on a finite interval has a complete discrete eigenbasis.
Figure 1.
Characteristic roots and bending shapes of a filament with two free ends
. Details.
Using the exponential definition of the hyperbolic cosine and the complex exponential function,
Thus
For real , the hyperbolic sine vanishes only at . Also , so the only zero of sine available to is . The product therefore vanishes exactly when or :
There are three possible parameter values. At the roots are and is real; at the roots themselves are real.
For the simple symmetric random walk, independence and the Rademacher distribution give the moment-generating function
The function is even; for its derivative is , since the derivative of is at most one and . Here and are the hyperbolic cosine and hyperbolic tangent. Hence the exponential-moment estimate is
The supplied exponential maximal bound for a symmetric random walk now gives . For , choose ; for , use the elementary bound by one. Thus
Fix , choose , and set . For large these form an increasing sequence. Apply the preceding exponential maximal bound for a symmetric random walk at with threshold . The resulting probability is at most
We have and . Choose strictly between and ; the displayed probabilities are . By the Borel-Cantelli first lemma, almost surely, eventually . For every , monotonicity of for large therefore gives
Intersect the probability-one events for to obtain the upper law of the iterated logarithm:
Integrating the resistive induction equation once introduces a constant :
Because , integrating this relation again over the channel gives , where is the volumetric flow rate per unit span. Put and , taking without loss of generality. Reversing the imposed field reverses but leaves and unchanged. Eliminating from the magnetohydrodynamic momentum equation gives
Its two zero wall values remove the odd homogeneous solution. Write . Integrating the first relation with the two magnetic wall values yields
Therefore . Substitution and integration produce the velocity and induced field:
Here is the Hartmann number. The hyperbolic cosine makes even, while the hyperbolic sine makes odd. Direct differentiation verifies both coupled equations and all four wall values.
Integrating the velocity gives the flux
The apparent singularity at is removable. The small-field limit is plane Poiseuille flow:
The cubic coefficient of is its first-order response in ; itself vanishes at zero imposed field. For a fixed positive , the flux decreases monotonically as increases. One exact way to see the sketch's monotonicity is the partial-fraction expansion of the hyperbolic cotangent:
whose derivative is strictly negative for . The flux has a horizontal tangent at and approaches zero as
The flux is even if signed is used.
For , put . Away from the walls,
The velocity has a nearly uniform core, with thin Hartmann layers of thickness enforcing the no-slip boundary condition. Near the upper wall, with fixed, and ; the lower wall follows by even/odd symmetry. Thus is negative for , positive for , and returns rapidly to zero at both walls. The required sketches show a decreasing flux, a flat velocity core, and an odd induced field with magnetic wall layers.
Figure 1.
Hartmann-flow flux versus magnetic field strength, and velocity and induced-field profiles at Hartmann number 20
. Details.
Write . The Taylor series of the cosine is
After expanding each power by the multinomial theorem, the coefficient of has absolute value when is even and is zero when is odd. The hyperbolic cosine
therefore has exactly the absolute values of the coefficients of . Thus is an entire majorant series for .
Write
The first column has unit length for the indefinite quadratic form:
Because , there is a unique such that
The two columns are orthogonal for the same form and the second has squared length :
The vectors with these properties are
Using the hyperbolic cosine and hyperbolic sine, the two possible families are therefore
Write , so the two circular boundaries are and . The zero angular mode must interpolate between and , giving .
For the th cosine mode, the radial factor has equal value at both boundaries. The unique harmonic function with those data is
Hence the solution of the annular Dirichlet problem is
At and , the hyperbolic cosine quotient equals one, so the boundary conditions are satisfied term by term.
Parametrize the finite cylindrical surface patch by
The surface area element is . The seam has zero surface area, so it causes no overcounting. Therefore
Unrolling the unit circular cylinder preserves surface area, which explains the equality with the area under the hyperbolic cosine in the figure. The patch consists of the cylindrical surface alone; no flat base is included.
For the given wrapped hyperbolic cosine curve,
so . The curve makes one full turn around the circular cylinder, descending from to and ascending to its initial point. Its two parameter endpoints coincide; their tangent vectors and do not agree, so the resulting closed curve has a corner there.
Figure 1.
The hyperbolic cosine curve wrapped around a unit circular cylinder, and the same surface unrolled
. The shaded finite patch extends from to . On the unrolled cylinder, its horizontal coordinate is ; the two vertical edges are identified. Details.
For separation of variables, set , with . The side Neumann boundary conditions give . As in the standard Sturm-Liouville problem, integration by parts shows . The zero eigenvalue gives constant , and together with makes constant too.
For positive eigenvalues, the side conditions give , , . The top condition chooses the hyperbolic cosine factor . Normalize it to one at , obtaining
The orthogonality of integer Fourier modes and the permitted interchange of summation and integration give
For the remaining Fourier coefficients, integrate by parts twice. The sine boundary term vanishes, and with ,
The last integral vanishes because and is a positive integer. Thus : only even Fourier modes survive, as expected from reflection symmetry about . Writing gives
The factor involving hyperbolic cosine lies between zero and one on the closed rectangle. Comparison with proves uniform convergence, so the solution is continuous on the entire boundary and takes the prescribed bottom values. On every strip , the Fourier coefficients have additional exponential decay, allowing differentiation and verifying Laplace's equation and all three Neumann boundary conditions.
To justify the endpoint calculation in b, differentiate at , in the interior-edge sense:
Subtracting from the hyperbolic tangent leaves an absolutely convergent series. For the remaining part, the power series gives, by taking real parts at and then ,
The same formula follows as the limit of the differentiated interior solution, since its leading multiplier is . This establishes the boundary normal derivative used in b. It also explains why the failure of the corner compatibility condition for mixed boundary data affects first derivatives at the bottom corners without spoiling continuity of itself.
Set in Laplace's equation. Separation of variables gives
The Sturm-Liouville problem in has the constant mode with , and the eigenfunctions with , . Indeed, integration by parts makes nonnegative; for the first Neumann boundary condition removes the sine term and the second requires to be an integer. For , both derivative conditions make constant.
Expand the top data using a Fourier cosine series:
The constant coefficient is
For , the constant part has zero integral against , and two integrations by parts give
Since the bottom value is zero, a convenient factor depending on for the nonconstant modes is , which vanishes at . Its value at is . The zero mode is linear in , so matching both constant boundary values gives
To obtain exactly the paper's form, use the hyperbolic-function identity
Thus
Each product of hyperbolic sine or hyperbolic cosine with has . The zero mode is also a harmonic function. On , the factor lies between zero and one. Since converges, the Weierstrass M-test gives uniform convergence on the closed rectangle. The series therefore gives zero at the bottom and the continuous Fourier cosine series of at the top, including its endpoints. On compact subsets with , the factors decay exponentially with , justifying termwise partial derivatives. In particular on both open vertical edges because every differentiated cosine factor vanishes there. The solution has finite Dirichlet energy: Fourier orthogonality gives the th nonconstant mode's energy as , whose sum is finite. Here is the hyperbolic cotangent. Thus a proves uniqueness in that class.
There is a corner compatibility condition for mixed boundary data in the printed wording. The top data have , but the paper writes even at . The series solves the interior equation and all the open-edge boundary conditions; the literal derivative condition at is incompatible with the top data. A function differentiable in along the top edge would necessarily have . Hence it cannot simultaneously satisfy the printed zero derivative at that corner, and there is no solution with all those literal corner requirements. The standard interpretation of this mixed boundary condition imposes the Neumann boundary condition on the open sides and the Dirichlet boundary condition continuously on the horizontal edges. Differentiating the series term by term on the top edge would be invalid at this corner; the uniform convergence of the undifferentiated series does not license that operation.
Substitute the exponential definitions of hyperbolic cosine and hyperbolic sine. The cross terms cancel:
This proves the identity for real arguments, and the same algebra proves it for complex arguments.
Differentiating the hyperbolic cosine with the chain rule gives . For real , , so the only zero of the derivative is , namely . Thus
Since , this stationary point is the global minimum.
Using the exponential definition of hyperbolic cosine, multiply by the nonzero quantity . This yields , hence .
If , its complex modulus requires , so ; its phase then requires . Thus all complex solutions are
There are no real solutions, since real .