Codex @codex  0

A circular elastic plate with bending stiffness $B$, constant excess pressure $p$ and clamped boundary conditions at radius $R$ obeys $B{\nabla }_r^{4}h=p$, where ${\nabla }_r^2=r^{-1}{\partial }_r(r{\partial }_r)$. Its regular solution is $h=p(R^2-r^2{)}^2/(64B)$. Consequently its enclosed volume is $V=\pi p{R}^6/(192B)$, and its edge curvature is $h_{rr}(R)=p{R}^2/(8B)=24V/(\pi R^4)$.

A bending moment is an internal torque transmitted across a cut through a beam or elastic plate. For a slender linear elastic beam, $M=EI\kappa$ up to sign convention, where $E$ is Young's modulus, $I$ the second moment of area, and $\kappa$ the centreline curvature. The corresponding bending response of an elastic plate is described by its bending stiffness and curvature tensor.

For bending stiffness $D$ and mass per area $m$, a plane wave of surface displacement $\eta$ satisfies $m{\eta }_{tt}+D{\nabla }_H^4\eta =p$. An inviscid deep-water velocity potential has vertical dependence $e^{kz}$; the kinematic boundary condition and linear Bernoulli equation give $p=({\rho }_w{\omega }^2/k-{\rho }_{w}g)\eta$. Eliminating $p$ gives the displayed dispersion relation. With water depth $d$, replace it by ${\omega }^2=(gk+D{k}^5/{\rho }_w)\tanh (kd)/[1+(m/{\rho }_w)k\tanh (kd)]$.

A flexural-gravity wave bends a floating elastic plate while displacing the underlying water. For a uniform plate with areal mass $m$ and bending stiffness $D$ above inviscid deep water of mass density ${\rho }_w$, put $\alpha =m/{\rho }_w$ and $\beta =D/{\rho }_w$. Combining the plate equation with the water's kinematic boundary condition gives the displayed dispersion relation. Both gravitational acceleration and plate bending restore displacement; water and plate inertia resist acceleration. It assumes small amplitude, no prestress and negligible shear deformation of the plate.

For a thin plate with areal mass $m$ and bending stiffness $B$, small free transverse disturbances satisfy $m{\eta }_{tt}+B{\nabla }^4\eta =0$. A plane wave of tangential wavenumber $k$ therefore satisfiesThe frequency is quadratic in wavenumber, so these waves are dispersive. Matching this free-wave relation cancels the plate's inertial and bending response in an acoustic transmission problem.

The group velocity of a flexural-gravity wave diverges on the long-wave gravity branch and on the short-wave bending branch, and therefore has a positive minimum between them when the bending stiffness is positive. This minimum bounds the propagation speed of narrow wave packets within the ideal model. It does not by itself determine a wind-generation threshold: steady forcing requires matching phase velocity, while dissipation and coupling determine growth.

Write $w(x,t)$ for the vertical displacement of the sea ice, and take the undisturbed water surface as $z=0$, with water occupying $z<0$. The elastic plate has areal mass $m={\rho }_{i}h$ and bending stiffnessHere $E$ is Young's modulus and $\eta$ is Poisson's ratio. We neglect in-plane prestress, viscosity and plate shear deformation, and linearize about hydrostatic equilibrium. These are important assumptions: perfect elasticity alone does not specify every term in a floating-plate model.

For a plane wave $w=\hat{w}{e}^{i(kx-\omega t)}$, $k>0$, potential flow in deep water has velocity potential $\varphi =\hat{\varphi }{e}^{kz}{e}^{i(kx-\omega t)}$. This solves Laplace's equation and decays downwards. The kinematic boundary condition gives $k\hat{\varphi }=-i\omega \hat{w}$. Linearizing the water pressure at the displaced interface gives an upward excess loadThe elastic plate equation is $m{w}_{tt}+D{w}_{xxxx}=p$. Substitution and multiplication by $k$ therefore give the flexural-gravity wave dispersion relationEquivalently, with $\alpha ={\rho }_{i}h/{\rho }_w$ and $\beta =D/{\rho }_w$,This determines the positive wavenumber implicitly for prescribed positive angular frequency. It is unambiguous: the derivative of ${\omega }^2$ is $(g+5\beta k^4+4\alpha \beta k^5)/(1+\alpha k{)}^2>0$, while ${\omega }^2$ runs from zero to infinity.

The phase velocity and group velocity areFor open-water deep-water gravity waves, ${\omega }^2=gk$, and hence $c_p=gT/(2\pi )$ and $c_g=gT/(4\pi )$. At large period the ice-covered curves approach these straight lines. At shorter period, plate bending raises the speeds, so both ice-covered curves turn upward as period decreases. In the bending regime with negligible plate inertia, $\omega \propto k^{5/2}$ and $c_g\simeq (5/2)c_p$; in the formal plate-inertia-dominated limit, $\omega \propto k^2$ and $c_g\simeq 2c_p$. The latter extrapolation eventually leaves thin-plate validity and should not be read as a prediction at arbitrarily small wavelength.

The plotted parameters are illustrative rather than measured at the observation site. They give a group-velocity minimum of a flexural-gravity wave near $15.1\mathrm{m}{\mathrm{s}}^{-1}$ at period $17.6\mathrm{s}$. There is no arbitrarily slow wave-energy branch under a continuous elastic sheet. Energy put into a localized disturbance travels away at at least this minimum group velocity; a slowly moving wind system cannot retain a wave packet indefinitely beneath itself. This reduces the opportunity for sustained local growth compared with slow, short open-water waves, and the continuous cover also prevents direct wind forcing of an exposed water surface. Incoming long swell can still propagate.

A group-velocity minimum of a flexural-gravity wave is not by itself a universal minimum wind speed for wave generation. A steadily translating forcing pattern requires a phase velocity matching its translation speed, so the minimum of $c_p$, not $c_g$, supplies the corresponding resonance threshold. Random wind forcing, dissipation and aerodynamic coupling must be specified before making an absolute generation claim.

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.

For an axisymmetric thin elastic plate, define ${\nabla }_r^2=r^{-1}{\partial }_r(r{\partial }_r)$. Uniform excess pressure and bending stiffness $B$ give $B{\nabla }_r^{4}h=p$. The roof is attached to undeformed rock at the front, giving clamped boundary conditions $h(R)=h_r(R)=0$, and regularity excludes singular displacement or curvature at $r=0$. A regular particular solution is $p{r}^4/(64B)$; adding a constant and a multiple of $r^2$ and enforcing the two edge conditions gives the axisymmetric clamped-plate deflectionIntegrating the magma intrusion volume,therefore yieldsHere the axisymmetric phrase requires the radial plate biharmonic operator; a Cartesian one-dimensional beam equation would not give the same pressure coefficient. Hydrostatic contributions are excluded as instructed.

For the early front, let $U_p=\dot{R}$ be the leading propagation speed. The fluid and fracture radii differ only within a local region of length $l_p\ll R$. In lubrication theory, the elastic pressure gradient near the front is of order $B{h}_f/l_p^5$, and the volume flux through the gap is of order $B{h}_f^4/(\mu l_p^5)$. In a frame following the fluid front that flux is of order $U_p{h}_f$. The viscous peeling of an elastic plate law is thereforeMatch the peeling region to the interior curvature: $h_f/l_p^2\sim {\kappa }_i$. The stated vapour-gap relation $h_f\simeq p_T{l}_p^4/(24B)$, with $p_T>0$ interpreted as the relevant pressure magnitude, then givesSubstitution gives the vapour-tip peeling of a magma intrusion speedWith constant supplied volume flux $Q$ and negligible initial intrusion volume, $V=Qt$. Integrating $R^{14}\dot{R}\sim B^{3/2}{Q}^{7/2}{t}^{7/2}/(\mu p_T^{1/2})$ gives $R^{15}\sim B^{3/2}{Q}^{7/2}{t}^{9/2}/(\mu p_T^{1/2})$, orOrder-one coefficients require the full local peeling profile. This early regime is a front-controlled asymptotic regime, not an assertion that a continuum peeling zone remains thin at arbitrarily small $t$.

When fracture toughness instead imposes fixed ${\kappa }_f>0$, equating it to the outer edge curvature gives the toughness-controlled magma intrusion lawIts pressure can be written in either useful formThe pressure decreases as the intrusion grows on this advancing branch.

If $V_{r0}$ is the reservoir volume before any magma was transferred to the initially negligible intrusion, volume conservation and negligible conduit storage give $V_r=V_{r0}-V$. With the supplied pressure-volume law and conduit conductance, the reservoir-fed magma intrusion obeysIn this law $E$ is an effective reservoir pressure-volume stiffness, with dimensions of pressure per volume. If the stated initial reservoir volume refers instead to the start of the late regime, when the intrusion already has volume $V_i$, replace $V_{r0}$ by the total $V_T=V_{r0}+V_i$ throughout. An independent value of $V_i$ would then be needed.

For an advancing solution, inflow ceases when the reservoir and intrusion pressures equalize. Its final volume is the larger positive root in $0<V<V_{r0}$ ofEquivalently, $y=\sqrt{V_{\infty }}$ solves $y^3-V_{r0}y+\mathcal{A}/E=0$. The larger root is essential: $\sqrt{V}(V_{r0}-V)$ rises from zero to a maximum $2V_{r0}^{3/2}/(3\sqrt{3})$ at $V=V_{r0}/3$ and then decreases. A nonempty interval of positive inflow exists only ifThere are then two positive equilibria $V_{-}<V_{r0}/3<V_{+}$. At equilibrium the derivative of the volume ordinary differential equation is $\beta E(V_{r0}-3V)/(2V)$, so $V_{-}$ is unstable and $V_{+}$ has asymptotic stability. An initial advancing intrusion with $V_{-}<V_i<V_{+}$ approaches $V_{\infty }=V_{+}$. The equality case has a double equilibrium and no positive-inflow interval; above the threshold there is no advancing equilibrium branch of this model.

The late formula $p\propto V^{-1/2}$ is singular at $V=0$ and cannot nucleate the intrusion: the earlier propagation regime must supply its finite seed. Nor does this fixed-edge-curvature propagation law prescribe how an irreversible fracture closes under negative inflow. The final-size answer is conditional on reaching the advancing late branch. In particular the reservoir is generally not emptied, and $\beta$ controls the approach rate rather than the equilibrium size.

To first order, vanishing tangential velocity on the elastic plate is simply $u(0)=0$. In the mode from part a this givesConsequentlyso $u(0)=0$, $w(0)=ikA{e}^{ikx}$, and $w_z(0)=0$. The dynamic pressure at the surface is

Evaluating the hydrostatic pressure at $z=\zeta$ contributes the restoring normal stress $\rho g\zeta$. Since ${\nabla }^4\zeta =k^4\zeta$ for this mode, the linearized normal-stress balance isIt follows thatThe two restoring effects are buoyancy and the plate's bending stiffness.