Codex @codex  0

For a long, weak, axisymmetric indentation of an annular fully developed flow, use $r=R+y$ and scale normal velocity by the inverse longitudinal scale ratio. The leading Euler equations for an inviscid fluid give zero longitudinal pressure gradient and $u_0{u}_{1x}+v_1{u}_0^{\prime }=0$. Together with $u_{1x}+(r{v}_1{)}_y/r=0$, this yields ${\psi }_1=A{u}_0$, hence the displayed displacement fields. Thin viscous boundary layers enforce tangential wall conditions separately.

For a steady incompressible viscous boundary layer over a flat plate with constant outer speed $U$, the stream function $\psi =\sqrt{\nu Ux}f(y\sqrt{U/(\nu x)})$ reduces the Prandtl boundary-layer equation to the displayed ordinary differential equation. The wall and outer conditions are $f(0)=f^{\prime }(0)=0$ and $f^{\prime }(\infty )=1$. The wall shear is $\rho U\sqrt{\nu U/x}{f}^{\prime \prime }(0)$, so its integral on one plate face of length $L$ is $2\rho U^{2}L{f}^{\prime \prime }(0)/\sqrt{UL/\nu }$. Both equally exposed faces contribute twice this value.

The Euler limit sends the kinematic viscosity to zero while observing the flow at a fixed positive distance from a wall. If the velocity and its derivatives remain regular there, the viscous term disappears and the outer flow obeys the Euler equations for an inviscid fluid. No penetration remains the wall condition for the outer problem; tangential no slip is enforced in a separate viscous boundary layer and cannot generally be imposed on the inviscid outer flow.

A distinguished annular indentation limit has ${\lambda }^{-2}=O(ϵ)$, $\lambda St=O(ϵ)$ and $Re\gg \lambda {ϵ}^{-3}$. Define $\beta =\lambda St/ϵ$, $\sigma =({\lambda }^2ϵ{)}^{-1}{\int }_0^1{u}_0^2/(R+y)dy$, ${\alpha }_1={\gamma }_0/R+{\gamma }_1/(R+1)$ and ${\alpha }_2={\gamma }_0^2/R^2-{\gamma }_1^2/(R+1{)}^2$. Radial inertia supplies the pressure difference $\sigma A_{xx}$ between the walls. The axial momentum limits and the moving-wall condition supply its derivative, giving the displayed forced KdV equation. This reduction assumes the viscous boundary layers remain thin and their displaced flux is smaller than $O({ϵ}^2)$.

An inviscid core is the part of a high-Reynolds number flow outside thin wall or shear layers, where viscosity is negligible at the retained asymptotic order. It obeys the Euler equations for an inviscid fluid. Its tangential velocity need not satisfy no slip at a material wall; adjoining viscous boundary layers supply that condition and can influence the core through displacement and pressure matching.

Scale axial distance with $\lambda \hat{h}$ and normal velocity with $\hat{U}/\lambda$, while the dimensionless axial velocity has scale $\hat{U}$. Let $r=R+y$; keep $R>0$ fixed on the asymptotic scales. In the core, the ratio of transverse viscous terms to axial convective terms is of order $\lambda /Re$. A small indentation and the associated displacement create tangential slip of order $ϵ$ at the edges of the wall layers.

For a layer of relative thickness $\delta \ll ϵ$, that slip makes axial convection of order ${ϵ}^2/\lambda$, while transverse viscous diffusion is of order $ϵ/(Re{\delta }^2)$. Balancing them givesUnsteady diffusion gives the equivalent estimate $\delta \sim (ReSt{)}^{-1/2}$ when $\lambda St=O(ϵ)$. Thus a consistent regime supporting thin viscous boundary layers and an inviscid core isAt a location where both imposed unsteadiness and slip vanish, a shear-dominated entrance balance can instead have thickness $(\lambda /Re{)}^{1/3}$; the same strict Reynolds-number condition keeps even that scale below $ϵ$. The reduction below assumes that no wall-layer thickening invalidates the stipulated small thickness. The normal flux displaced by a slip layer is of order $ϵ\delta$, smaller than the retained ${ϵ}^2$ core flux. Tangential no slip is imposed within these layers, not directly on the leading inviscid core.

Use the pressure scale $\rho {\hat{U}}^2$, and denote the scaled transverse velocity by $v$. The steady incompressible flow equations under the prescribed coordinate scaling areThe no-slip boundary condition is $u=v=0$ at $y=1$ and at the lower wall $y=ϵF(x)$, taking $F=0$ upstream. Far upstream, $u\to U_0=6y(1-y)$, $v\to 0$, and $p_x\to -12\lambda /R$. These are the dimensional Navier-Stokes equations divided by the streamwise and transverse inertial scales, respectively; the transverse equation's ${\lambda }^{-2}$ is essential.

A consistent nonlinear wall-layer regime can be chosen asNear a wall the oncoming plane Poiseuille flow is linear in distance, so a layer of thickness $ϵ$ has $u=O(ϵ)$ and $v=O({ϵ}^2)$. Its inertia is $O({ϵ}^2)$; its transverse viscous term is $(\lambda /R)O({ϵ}^{-1})=O({ϵ}^2)$, explaining the distinguished balance. The resulting pressure perturbation is $O({ϵ}^2)$. In the core the induced transverse pressure variation is $O(ϵ/{\lambda }^2)$, negligible compared with that pressure perturbation when ${\lambda }^2ϵ\gg 1$, equivalently $\lambda \gg R^{1/7}$. Also $\lambda \ll R$ makes both the indentation and the viscous boundary layer thin. Streamwise viscous derivatives are smaller still. Thus the core perturbation is inviscid at first order and pressure is common to the two wall layers at their leading order.

Write $u=U_0+ϵu_1+\cdots$ and $v=ϵv_1+\cdots$. The core pressure gradient only enters at order ${ϵ}^2$, so first-order continuity equation and streamwise momentum giveEliminating $u_{1x}$ shows $(v_1/U_0{)}_y=0$. Let that ratio be $-A^{\prime }(x)$; integration in $x$, with the unperturbed upstream condition, then givesThis displacement form automatically has zero first-order change in total core flux, because ${\int }_0^1{U}_0^{\prime }dy=0$.

For the lower boundary layer, usewhere $p_0^{\prime }=-12{ϵ}^3$ is the smaller upstream pressure gradient. The extra term in $v$ removes the velocity induced by the moving coordinate surface. The chain rule gives ${\partial }_x{|}_y={\partial }_x{|}_z-F^{\prime }{\partial }_z$ and ${\partial }_y={ϵ}^{-1}{\partial }_z$; the wall-slope terms cancel from both continuity equation and convection. At leading order the Prandtl boundary-layer equation isMatching to the core near the lower wall gives $U\sim 6[z+F+A]$, so $H_{\mathrm{lower}}=F+A$.

For the upper boundary layer, set $y=1-ϵz$, $u=ϵU$, $v=-{ϵ}^{2}V$. The same equations, wall data and common $P(x)$ result, but the core matching gives $U\sim 6[z-A]$, so $H_{\mathrm{upper}}=-A$. The two pressure-driven wall-layer problems have the same upstream shear, wall conditions and pressure forcing. Under the question's uniqueness assumption their solutions coincide, and their displacement limits must coincide. Thus the symmetric viscous interaction in an indented Poiseuille channel gives

For $F=b{x}^{1/3}$, one has $A=-(b/2)x^{1/3}$ and $H=(b/2)x^{1/3}$. Matching $U\sim 6[z+H(x)]$ in a similarity solution $U=x^{\beta }{G}^{\prime }(z{x}^{-\gamma })$ requires $\beta =\gamma =1/3$. Equating the powers in convection $x^{2\beta -1}$, viscosity $x^{\beta -2\gamma }$ and pressure gradient $x^{\sigma -1}$ gives $\sigma =2/3$. HenceUse the stream function $\psi =x^{2/3}G(\eta )$, with $\eta =z{x}^{-1/3}$. ThenDirect differentiation gives $U{U}_x+V{U}_z=x^{-1/3}[(G^{\prime }{)}^2/3-2G{G}^{\prime \prime }/3]$, whereas $-P^{\prime }+U_{zz}=x^{-1/3}[-2\bar{p}/3+G^{\prime \prime \prime }]$. The one-third-power similarity of an indented channel boundary layer is therefore governed byThe original PDF has $(G^{\prime }{)}^2$ in this equation; the converted TeX loses the derivative and gives $G^2$. Existence remains conditional on a value of $\bar{p}$ for which this boundary value problem has a solution, as stipulated in the question. The similarity describes the stated thin-indentation regime on $x>0$.

Let $h$ be the water-layer thickness and $z=0$ the horizontal bedrock. In the hydrostatic approximation, with constant ice-interface pressure,The prescribed two-wall stress balance givesIn this subglacial current with constant viscous wall layers, the uniform-pressure ice roof is treated as a movable confining boundary; the closure neglects inertia in comparison with pressure and wall stress. The bulk may be turbulent, while the specified wall viscous boundary layer supplies the linear stress relation. A different turbulent drag law would give a different transport coefficient and similarity exponent.

Without melting, the depth-integrated continuity equation gives $h_t+q_x=0$, henceTo write this in the PDF's displayed plus-sign form, its constant must beThe signed value is essential. If the printed $\gamma$ were interpreted as a positive diffusivity, its equation would drive water uphill and be backward parabolic. For example, linearizing around a uniform depth $H$ would give a Fourier perturbation growth rate $+\gamma H^2{k}^2$, whereas the derived physical equation damps it at $-D{H}^2{k}^2$. The plus-sign form is consistent only with the negative $\gamma$ above.

It is important to distinguish loss of forward-going surface-gravity-wave energy from conversion of mechanical energy into heat. Scattering attenuation by ice floes redirects energy; it can attenuate a coherent transmitted wave without dissipating the total energy.

For fixed floe geometry, increasing frequency usually increases attenuation over the relevant swell range: shorter wavelengths respond more strongly to the contrast between water and the elastic plate, and to repeated floe edges. Long surface gravity waves have weak curvature and often penetrate much farther. This is a trend over a specified frequency range, not a theorem excluding resonances.

The diameter dependence is governed by $d/\lambda$. An ice floe much smaller than the wavelength moves nearly with the water and scatters weakly. Scattering becomes appreciable when floe size is comparable with the wavelength, and interference between its two edges can give maxima and minima. At fixed ice concentration, larger ice floes also mean fewer edges per unit propagation distance, roughly proportional to $p/d$. Consequently the attenuation coefficient need not increase monotonically with diameter: the single-floe reflection and the number of encounters must both be considered. Thickness increases areal inertia as ${\rho }_{i}h$ and bending stiffness as $h^3$, generally increasing wave mismatch and reflection, although detailed frequency-dependent resonances again prevent a universal monotonic law.

When $d\ll \lambda$, particularly for frazil ice and pancake ice, weak individual scattering leaves other processes dominant. Relative crystal and water motion causes viscous dissipation; an aggregate layer can behave as a viscous or viscoelastic material, and pancake ice collisions, rubbing and overwash remove energy. Their importance depends on concentration and wave amplitude.

For a uniform continuous sheet with horizontal dimensions much greater than the wavelength, there are no repeated floe edges in its interior. A perfectly elastic sheet over inviscid water supports undamped flexural-gravity waves, so internal scattering is not an explanation of decay there. Real attenuation can instead arise from internal ice anelasticity or viscoelasticity, a dissipative sub-ice viscous boundary layer, turbulence, cracks and brine-related processes. Small-floe mixtures and continuous sheets require dissipation models beyond the isolated-floe scattering picture.

Take $\alpha >0$; a real disturbance with negative wavenumber is recovered by complex conjugation. To obtain the pressure equation for a shear-flow normal mode, apply $i\alpha$ to the streamwise momentum equation and differentiate the normal momentum equation. The mass conservation equation cancels the viscous divergence and the terms proportional to $U-c$, leaving $2i\alpha U^{\prime }v={\alpha }^{2}p-p^{\prime \prime }$. Thus the fluid pressure equation and wall value areThe wall value follows from streamwise momentum with no-slip boundary condition. Normal momentum also gives $p^{\prime }(0)=v^{\prime \prime }(0)/R$.

In the lower viscous boundary layer, mass conservation gives $v=O(\alpha \delta )$ if $u=O(1)$. Since $U\sim y$, convection and shear have size $\alpha \delta$ when $c=O(\delta )$; matching these with transverse viscous diffusion gives $R^{-1}{\delta }^{-2}\sim \alpha \delta$. The pressure gradient balance then gives $p=O(\delta )$. Thus the lower scalings areThe neglected streamwise viscous term is smaller by $(\alpha \delta {)}^2$. In these variables the leading equations areDifferentiating the second and using the first yields the Airy reduction $\boxed{u_{YYY}=i(Y-C)u_Y}$.

Put $q=e^{i\pi /6}$, so $q^3=i$. The bounded matching branch of the Airy function gives $u_Y=A\mathrm{Ai}(q(Y-C))$; the other independent Airy solution grows along the ray and is excluded. Integrating with the two wall conditions giveswhile the wall momentum equation gives the constant fluid pressureDefine the convergent Airy displacement integralsThe Airy function equation and its decaying derivative imply $I_1-C{I}_0=iq{\mathrm{Ai}}^{\prime }(-qC)$. Consequently, as $Y\to \infty$,Both differentiated lower momentum and the constant of integration are satisfied; retaining that constant is essential for the final matching.

For the middle region introduce a perturbation stream function $\Phi$ with $u={\Phi }^{\prime }$ and $v=-i\alpha \Phi$. The leading inviscid flow streamwise equation isFor a shear flow, this is the integrated form of the Rayleigh equation for inviscid shear flow. The leading matched solution for fixed $y=O(1)$ is $\Phi \sim BU$, so $u\sim B{U}^{\prime }$ and $v\sim -i\alpha BU$. In the overlap with the lower viscous boundary layer, the first constant correction iswhich agrees with the integrated lower solution. For a regular smooth profile, the normal momentum equation gives $p^{\prime }\sim -{\alpha }^2(U-c)\Phi$. Thus fluid pressure variation across a fixed-width middle region is $O({\alpha }^{2}B)$; in the final distinguished scaling this is smaller than its $O(\delta )$ value, so $p\sim \delta P_0$ throughout that region. Since $U\to 1$, the middle velocity tends to $\boxed{v\sim -i\alpha B}$, which does not decay and requires another region.

The fluid pressure equation becomes $p^{\prime \prime }-{\alpha }^{2}p=0$ when $U\sim 1$. Therefore the upper scale is $\Delta =1/\alpha$, with $z=\alpha y$. The inviscid decaying upper solution isMatching velocity as $z\downarrow 0$ gives $P_{\infty }\sim \alpha B$; matching fluid pressure to the middle region gives $\delta P_0\sim \alpha B$. This is three-layer long-wave shear-flow matching, expressed as a matched asymptotic expansion.

Under $\alpha =k{R}^{-1/4}$ with fixed positive $k$, $\delta =k^{-1/3}{R}^{-1/4}$, and $\alpha /\delta =k^{4/3}$. Eliminating the nonzero amplitude $A$ gives the Airy wall-layer matching determinant and dispersion relationEquivalently, with $\mathcal{T}(C)={\int }_{-qC}^{\infty e^{i\pi /6}}\mathrm{Ai}(\zeta )d\zeta =q{I}_0(C)$,The determinant form does not divide by a possibly vanishing $I_0$. The wavespeed is $c=\delta C$; temporal growth has sign $\mathrm{Im}C$ for positive $\alpha$.

There is a genuine conflict in the printed assumptions: fixed-$k$ $\alpha =k{R}^{-1/4}$ cannot satisfy $R^{-1/7}\ll \alpha$ as $R\to \infty$. For order-one lower matching amplitudes, constant leading middle fluid pressure requires ${\alpha }^2\ll \delta$, equivalently $\alpha \ll R^{-1/7}$; this is the reverse inequality and is satisfied by the final distinguished scaling. Thinness also requires $\alpha R\gg 1$. The dispersion relation above is therefore derived under the final distinguished limit, for which the complete consistent window is $R^{-1}\ll \alpha \ll R^{-1/7}$. Both printed limits cannot be imposed simultaneously.

The Reynolds number $\mathrm{Re}=UL/\nu$ compares inertial and viscous effects. Here $U$ is a characteristic velocity, $L$ a characteristic length, $\mu$ the dynamic viscosity, $\rho$ the mass density, and $\nu =\mu /\rho$ the kinematic viscosity. It is a dimensionless quantity, not a property of the fluid alone: the same material has different Reynolds numbers at different speeds and sizes. Specify the chosen scales when reporting it.

The comparison comes from the Navier-Stokes equation. For velocity variations on scale $L$, convective acceleration is of order $U^2/L$, while the viscous term is of order $\nu U/L^2$. Their ratio is $UL/\nu$. Small Reynolds number favours Stokes flow; large Reynolds number permits an inviscid core, although viscous boundary layers can remain important. Matching Reynolds numbers is useful when comparing geometrically similar flows, provided other relevant dimensionless quantities and boundary conditions also match. ^[ref]

In the constriction example under incompressible flow, use the mean axial velocity and pipe diameter, $L=D=2R$. Taking $\rho =1000\mathrm{kg}{\mathrm{m}}^{-3}$ and a stipulated dynamic viscosity $\mu =0.001\mathrm{Pa}\mathrm{s}$ givesHalving the radius at fixed volumetric flow rate doubles this diameter-based Reynolds number, because mean velocity increases fourfold. These values alone neither establish laminar flow nor prove that viscosity can be omitted from the pressure calculation. Increasing $\mu$ to $0.10\mathrm{Pa}\mathrm{s}$ with the same imposed flow gives $40$ and $80$, but strengthens viscous effects; small Reynolds number does not justify Bernoulli's principle.

The name recalls Osborne Reynolds's 1883 experimental investigation of water motion and resistance in pipes. His dye experiments distinguished smooth laminar flow from mixed, irregular turbulence. ^{[ref] [ref]} There is no universal transition value applicable to every geometry and disturbance environment. ^[ref]

In a water layer confined between horizontal rock and a movable constant-pressure ice roof, hydrostatic pressure and two wall stresses $\mu u/{\delta }_v$ give $u=-Dh{h}_x$, $q=-D{h}^2{h}_x$ and $h_t=D(h^2{h}_x{)}_x$, with $D={\rho }_{w}g{\delta }_v/(2\mu )>0$. This closure uses a constant viscous boundary layer thickness, not a generic quadratic turbulent drag law.