Your Wind Farm Simulation Might Be Wrong: Compressibility Matters More Than You Think for Atmospheric Gravity Waves

Boussinesq approximation underestimates wave displacements at altitude and mischaracterizes wave angles at low Froude numbers — exactly the regimes that matter for wind farm siting.

Imagine you're simulating ocean waves in a bathtub. The water in your tub is incompressible — push one end down, the other end rises instantly and proportionally. Now imagine the real ocean, where water density changes with depth, temperature, and salinity, and where deep pressure actually compresses the fluid. Your bathtub model gets the wave shapes roughly right, but it misses how much bigger the waves get in deep water and how they bend at the interface between layers. That's the core of this paper: the atmospheric equivalent of realizing your bathtub isn't the ocean. The committed claim: solving the fully compressible Navier-Stokes equations for stratified flow over topography produces materially different gravity wave characteristics compared to the incompressible Boussinesq approximation — specifically larger parcel displacements at altitude, greater variation in apparent wavelengths, and different wave angles, especially in low Froude number (strongly stratified) regimes. This matters because Boussinesq is the workhorse assumption behind most atmospheric boundary layer simulations used for wind energy. The study compares three formulations head-to-head: linearized incompressible Boussinesq, full nonlinear incompressible Boussinesq, and fully compressible Navier-Stokes. The test case is stratified flow over bell-shaped hills — a classical setup in mountain wave literature that isolates gravity wave generation cleanly. All three produce qualitatively similar wave structures, so the Boussinesq approximation isn't catastrophically wrong. But quantitatively, the differences grow with altitude and with decreasing Froude number. The compressible solver predicts larger vertical parcel displacements at height, which directly affects how much kinetic energy is available or unavailable to downstream turbines. The architecture here is computational fluid dynamics at the mesoscale atmospheric level — three separate solver formulations applied to the same idealized terrain geometry. The linearized Boussinesq is analytical (closed-form mountain wave theory), the full Boussinesq is a nonlinear numerical solve, and the compressible solver adds density as a prognostic variable coupled to the energy equation. The key structural choice is the comparison methodology itself: same geometry, same inflow, three physics levels. This is a controlled numerical experiment, not a field campaign. On integrity, this is simulation-vs-simulation with no field validation. The bell-shaped hill is an idealized canonical test case, which is both a strength (clean isolation of the effect) and a limitation (real terrain is not a bell curve). The comparison across three equation sets is internally consistent, and the use of classical mountain wave theory as a baseline is appropriate. But the absence of any observational data or field measurements means we're watching models argue with each other. The authors are transparent about this — they frame the work as a sensitivity study, not a field validation. The practical stakes are real. Wind farm layout optimization, wake modeling, and power prediction all depend on getting atmospheric gravity waves right. Gravity waves triggered by terrain upstream of a wind farm, or by the farm's own blockage effect, can modulate hub-height wind speeds by several percent. If the Boussinesq approximation underestimates wave amplitude at turbine height in strongly stratified conditions, then the performance predictions feeding investment decisions are systematically biased. The low Froude number regime — where differences are largest — corresponds to stable nighttime conditions, which is precisely when gravity wave effects on wind farms are most pronounced. The gap this paper doesn't fill is the obvious one: running the same comparison with a real wind farm in the domain and checking whether the compressibility correction actually changes predicted annual energy production by an amount that matters commercially. That experiment requires coupling to a turbine model and running at sufficient resolution over realistic terrain, which is computationally expensive. The authors are almost certainly planning this — the paper reads as the foundational validation step before the applied wind-energy paper.