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.
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.
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
.
Put , so . The power series coefficients of the hyperbolic sine preprocessing are
Thus for every , and the series converge on the full interval. Applying coefficient-dominated entrywise positivity gives .
Every diagonal entry belongs to one of the diagonal blocks and equals . Because and , this is . Therefore
The matrix lies in the elliptope and admits a Gram matrix representation by unit vectors. This preprocessing is the Krivine rounding scheme; the equality of absolute coefficients is what preserves positive semidefiniteness even though the cross-block sine coefficients alternate in sign.
The three fixed- estimates below cease to be uniform as the minimum of the phase approaches the endpoint. For the cubic endpoint-to-saddle transition, write
The Taylor series of the hyperbolic sine gives, for bounded and fixed ,
The cubic term controls the tail, so localization of the Laplace integral gives the uniform leading formula
Here is the cubic Laplace transition integral, equal to in terms of the Scorer Hi function.
To recover the endpoint regime, let . Scale in ; the cubic term becomes negligible and . Therefore
which agrees with part (i) in the overlap . At , the Gamma integral gives , recovering part (iii).
For , set . The exponent has its maximum at , with second derivative . Thus Laplace's method gives
and the transition formula becomes
For with , part (ii) has
Putting reproduces both the exponential and its prefactor. For relative agreement of these leading exponentials, one may use the overlap , which makes . Thus the same transition integral connects all three regimes.
Figure 1.
The cubic endpoint-to-saddle transition
. Direct numerical integration of the original phase approaches the same cubic transition function as the large parameter increases. Negative transition parameter places the minimum inside the interval; positive parameter leaves an ordinary endpoint minimum.
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