SUMMARYSmoothed particle hydrodynamics (SPH) is becoming increasingly common in the numerical simulation of complex fluid flows and an understanding of the errors is necessary. Recent advances have established techniques for ensuring completeness conditions (low-order polynomials are interpolated exactly) are enforced when estimating property gradients, but the consequences on errors have not been investigated. Here, we present an expression for the error in an SPH estimate, accounting for completeness, an expression that applies to SPH generally. We revisit the derivation of the SPH equations for fluids, paying particular attention to the conservation principles. We find that a common method for enforcing completeness violates a property required of the kernel gradients, namely that gradients with respect the two position variables be equal and opposite. In such models this means conservation principles are not enforced and we present results that show this. As an aside we show the summation interpolant for density is a solution of, and may be used in the place of, the discretized, symmetrized continuity equation. Finally, we examine two examples of discretization errors, namely numerical boundary layers and the existence of crystallized states.

Recent basin‐scale changes to the compactness and thickness of Arctic sea ice foreshadow that encroaching swells and locally generated waves will exert more influence there in the future. Indeed, it is conceivable that waves may have already hastened the adjustments observed by breaking up ice floes. Yet waves advancing in sea ice attenuate due to being scattered from ice thickness variations and damped by ice inelasticity, turbulence and friction. While past research focuses on scattering by unnaturally perfect features in the ice, the model reported herein assimilates realistic basin‐scale swathes of heterogeneous ice and parameterizes damping. By way of example, we show how an ocean wave train evolves during its passage in an 1670‐km‐long Arctic sea ice profile obtained from submarine.

Flexural oscillations of floating sea ice sheets induced by ocean waves travelling at the boundary between the ice and the water below can propagate great distances. But, by virtue of scattering, changes of ice thickness and other properties encountered during the journey affect their passage, notwithstanding attenuation arising from several other naturally occurring agencies. We describe here a two-dimensional model that can simulate wave scattering by long (approx. 50 km) stretches of inelastic sea ice, the goal being to replicate heterogeneity accurately while also assimilating supplementary processes that lead to energy loss in sea ice at scales that are amenable to experimental validation. In work concerned with scattering from solitary or juxtaposed stylized features in the sea ice canopy, reflection and transmission coefficients are commonly used to quantify scattering, but on this occasion, we use the attenuation coefficient as we consider that it provides a more helpful description when dealing with long sequences of adjoining scatterers. Results show that scattering and viscosity both induce exponential decay and we observe three distinct regimes: (i) low period, where scattering dominates, (ii) high period, where viscosity dominates, and (iii) a transition regime. Each regime's period range depends on the sea ice properties including viscosity, which must be included for the correct identification of decay rate.

No abstract

Arctic sea-ice contains imperfections such as cracks, leads and pressure ridges that scatter flexuralgravity waves. Models for predicting scattering have been described in the literature, concentrating mainly on singular isolated features with simplified shapes or on arrays of such features. In reality ridges are seldom simple and leads are rarely entirely free of ice. Here we describe a model in which the scattering by a sheet of arbitrary thickness can be simulated. Linear wave theory and Green's functions are used to derive the governing equations for a numerical model of a two-dimensional (in the vertical) system. We examine wave scattering by random ice sheets, identifying trends in behavior as the wave period and the length, median thickness and variance of the sheet are changed. It has been suggested that wave scattering could be used to identify sea-ice thickness, a task which is difficult or expensive by other methods, and here we examine a technique by which this could potentially be achieved. However, a large data base is necessary for this to work and this may limit the practicality of the approach.