What is Harbor Wave Diffraction?
🙋
What exactly is wave diffraction, and why does it matter for harbors?
🎓
Basically, it's how waves bend and spread out after passing through a gap, like an opening in a breakwater. In practice, this bending determines if the water inside a harbor stays calm or gets choppy. For instance, a ship trying to dock needs calm water, which depends on how much wave energy sneaks in. Try moving the "Gap Width (B)" slider above to see how a wider opening lets more wave energy through.
🙋
Wait, really? So the wave height inside isn't just the same as outside? What's this K_D value the simulator shows?
🎓
Exactly right! The wave height changes. K_D is the diffraction coefficient. It's the ratio of the wave height at a point inside the harbor to the incident wave height outside. A common case is K_D = 0.5, meaning the wave is half as high inside. When you change the "Incident Angle (θ)" parameter, you'll see K_D change in real-time because waves hitting the gap at an angle diffract differently.
🙋
That makes sense. But the math looks intense with "Fresnel integrals". What's the physical idea behind the equation?
🎓
In plain language, the theory treats the wave crests like light waves. The gap acts like a new source, and waves from different parts of this source interfere with each other, creating a pattern of high and low spots. The Fresnel integral sums up all these tiny wave contributions. The simulator does this complex math instantly. For instance, adjust the "Water Depth (h)" and see how it affects the wavelength and thus the diffraction pattern—shallower water makes waves "feel" the bottom and slow down, changing everything.
Physical Model & Key Equations
The core of this simulator is Sommerfeld's theory for diffraction by a semi-infinite barrier (a very long breakwater with one opening). The wave height at any point behind the barrier is found by summing the contributions from the wave front at the gap edge.
$$K_D(\xi) = \frac{1}{2}\left|F(\xi_+) + F(\xi_-)\right|$$
Here, $K_D$ is the diffraction coefficient (relative wave height). $F(\xi)$ is the complex Fresnel integral. The parameters $\xi_+$ and $\xi_-$ are dimensionless coordinates that depend on the observer's position $(x, y)$, the wavelength $L$, and the gap geometry. They essentially tell us "how far" the point is from the geometric shadow boundaries.
The wave power, or energy flux, is what engineers really care about for harbor tranquility. It's directly related to the square of the wave height.
$$P = \frac{\rho g^2 H^2 T}{32 \pi}$$
Where $P$ is wave power per unit crest width (W/m), $\rho$ is seawater density (~1025 kg/m³), $g$ is gravity (9.81 m/s²), $H$ is wave height, and $T$ is wave period. Inside the harbor, the diffracted wave power is $P_{harbor}= K_D^2 \times P_{incident}$. This is why a small $K_D$ leads to a dramatic reduction in wave energy.
Real-World Applications
Harbor Design & Breakwater Layout: Engineers use these exact calculations to optimize the length, orientation, and gap width of breakwaters. A common case is placing a main breakwater with a small, angled gap to allow ship entry while minimizing wave penetration into the mooring basin.
Marina Tranquility Analysis: For yacht marinas, even small waves can damage boats. Diffraction analysis predicts "hot spots" of wave energy behind breakwaters, helping to position floating docks in the calmest zones. Changing the wave period (T) in the simulator mimics different storm conditions.
Validation for Numerical Wave Models: Advanced CAE tools like MIKE 21 or OpenFOAM solve the full fluid equations. The Sommerfeld solution provides a precise, analytical benchmark to validate these complex numerical simulations in their far-field results.
Coastal Structure Retrofit: When an existing harbor needs upgrading, engineers model diffraction to see if extending a breakwater or narrowing an entrance will effectively reduce wave heights at critical infrastructure like cargo terminals or ferry berths.
Common Misunderstandings and Points to Note
When you start using this tool, there are a few points that are easy to misunderstand. First, the diffraction coefficient K_D is a ratio of wave heights, not of energy itself. If K_D=0.5, the wave height is halved, but since wave energy (an indicator of wave force or overtopping) is proportional to the square of the wave height, the energy is actually reduced to 0.25 times. In design, don't be complacent thinking "the wave height is small, so it's OK"; make sure to properly estimate the forces acting on the structure.
Next, water depth h and wave period T are not independent parameters. While you can change them separately in the tool, real waves are governed by the dispersion relation $L = (gT^2/2\pi) \tanh(2\pi h/L)$, which links water depth and wavelength L. For example, when a wave with a 10-second period enters a shallow area with a 5m depth, its wavelength becomes about 50m (compared to about 150m in deep water). This "shoaling transformation" also changes how diffraction spreads, so when using real data, pay attention to the combination of water depth and period.
Finally, understand the fundamental limitation that this calculation assumes a "steady, single wave". The actual sea consists of irregular waves with a directional spectrum. Therefore, even in an area calculated by the tool to have K_D=0.3 (a calm area), wave components from other directions might intrude. In practice, the standard workflow is to use these results as a basis, then verify with more advanced numerical simulations that consider irregular waves and directional distribution.
Worked Example
Harbor with 8 m incident swell, 12 s period, 15 m depth, 40 m breakwater gap. Wavelength L = 9.81×144/6.28 ≈ 225 m. At harbor center 120 m from gap: K_D ≈ 0.35 (Sommerfeld solution). Harbor wave height = 0.35 × 8 = 2.8 m. Sheltering = (1 − 0.35) × 100 = 65%. Energy reduction from 7.85 kW/m incident to 0.96 kW/m sheltered (proportional to K_D²). H/L steepness = 8/225 = 0.036 (mild swell regime).
