Jeroen Molemaker

Low viscosity flow simulations for animation

We present a combination of techniques to simulate turbulent fluid flows in 3D. Flow in a complex domain is modeled using a regular rectilinear grid with a finite-difference solution to the incompressible Navier-Stokes equations. We propose the use of the QUICK advection algorithm over a globally high resolution grid. To calculate pressure over the grid, we introduce the Iterated Orthogonal Projection (IOP) framework. In IOP a series of orthogonal projections ensures that multiple conditions such as non-divergence and boundary conditions arising through complex domains shapes or moving objects will be satisfied simultaneously to specified accuracy. This framework allows us to use a simple and highly efficient multigrid method to enforce non-divergence in combination with complex domain boundary conditions. IOP is amenable to GPU implementation, resulting in over an order of magnitude improvement over a CPU-based solver. We analyze the impact of these algorithms on the turbulent energy cascade in simulated fluid flows and the resulting visual quality.

Small‐scale open ocean currents have large effects on wind wave heights

Tidal currents and large-scale oceanic currents are known to modify ocean wave properties, causing extreme sea states that are a hazard to navigation. Recent advances in the understanding and modeling capability of open ocean currents have revealed the ubiquitous presence of eddies, fronts, and filaments at scales 10–100 km. Based on realistic numerical models, we show that these structures can be the main source of variability in significant wave heights at scales less than 200 km, including important variations down to 10 km. Model results are consistent with wave height variations along satellite altimeter tracks, resolved at scales larger than 50 km. The spectrum of significant wave heights is found to be of the order of 70〈Hs〉2/(g2〈Tm0,−1〉2) times the current spectrum, where 〈Hs〉 is the spatially averaged significant wave height, 〈Tm0,−1〉 is the energy-averaged period, and g is the gravity acceleration. This variability induced by currents has been largely overlooked in spite of its relevance for extreme wave heights and remote sensing.

Submesoscale Dynamics in the Northern Gulf of Mexico. Part III: Lagrangian Implications

AbstractFour numerical simulations are used to characterize the impact of submesoscale circulations on surface Lagrangian statistics in the northern Gulf of Mexico over 2 months, February and August, representative of winter and summer. The role of resolution and riverine forcing is explored focusing on surface waters in regions where the water column is deeper than 50 m. Whenever submesoscale circulations are present, the probability density functions (PDFs) of dynamical quantities such as vorticity and horizontal velocity divergence for Eulerian and Lagrangian fields differ, with particles preferentially mapping areas of elevated negative divergence and positive vorticity. The stronger the submesoscale circulations are, the more skewed the Lagrangian distributions become, with greater differences between Eulerian and Lagrangian PDFs. In winter, Lagrangian distributions are modestly impacted by the presence of the riverine outflow, while increasing the model resolution from submesoscale permitting to sub...

Ocean convergence and the dispersion of flotsam

Significance Ocean currents move material released on the ocean surface away from the release point and, over time, spread it over an increasingly large area. However, observations also show high concentrations of the material even after significant spreading. This work examines a mechanism for creating such concentrations: downwelling of water at the boundaries of different water masses concentrates floating material at this boundary. Hundreds of satellite-tracked drifters were released near the site of the 2010 Deepwater Horizon oil spill. Surprisingly, most of these gathered into a single cluster less than 100 m in size, dramatically demonstrating the strength of this mechanism. Floating oil, plastics, and marine organisms are continually redistributed by ocean surface currents. Prediction of their resulting distribution on the surface is a fundamental, long-standing, and practically important problem. The dominant paradigm is dispersion within the dynamical context of a nondivergent flow: objects initially close together will on average spread apart but the area of surface patches of material does not change. Although this paradigm is likely valid at mesoscales, larger than 100 km in horizontal scale, recent theoretical studies of submesoscales (less than ∼10 km) predict strong surface convergences and downwelling associated with horizontal density fronts and cyclonic vortices. Here we show that such structures can dramatically concentrate floating material. More than half of an array of ∼200 surface drifters covering ∼20 × 20 km2 converged into a 60 × 60 m region within a week, a factor of more than 105 decrease in area, before slowly dispersing. As predicted, the convergence occurred at density fronts and with cyclonic vorticity. A zipperlike structure may play an important role. Cyclonic vorticity and vertical velocity reached 0.001 s−1 and 0.01 ms−1, respectively, which is much larger than usually inferred. This suggests a paradigm in which nearby objects form submesoscale clusters, and these clusters then spread apart. Together, these effects set both the overall extent and the finescale texture of a patch of floating material. Material concentrated at submesoscale convergences can create unique communities of organisms, amplify impacts of toxic material, and create opportunities to more efficiently recover such material.

Observations of Near-Surface Current Shear Help Describe Oceanic Oil and Plastic Transport

Plastics and spilled oil pose a critical threat to marine life and human health. As a result of wind forcing and wave motions, theoretical and laboratory studies predict very strong velocity variation with depth over the upper few centimeters of the water column, an observational blind spot in the real ocean. Here we present the first-ever ocean measurements of the current vector profile defined to within 1 cm of the free surface. In our illustrative example, the current magnitude averaged over the upper 1 cm of the ocean is shown to be nearly four times the average over the upper 10 m, even for mild forcing. Our findings indicate that this shear will rapidly separate pieces of marine debris which vary in size or buoyancy, making consideration of these dynamics essential to an improved understanding of the pathways along which marine plastics and oil are transported.

Topographic and Mixed Layer Submesoscale Currents in the Near-Surface Southwestern Tropical Pacific

AbstractThe distribution and strength of submesoscale (SM) surface layer fronts and filaments generated through mixed layer baroclinic energy conversion and submesoscale coherent vortices (SCVs) generated by topographic drag are analyzed in numerical simulations of the near-surface southwestern Pacific, north of 16°S. In the Coral Sea a strong seasonal cycle in the surface heat flux leads to a winter SM “soup” consisting of baroclinic mixed layer eddies (MLEs), fronts, and filaments similar to those seen in other regions farther away from the equator. However, a strong wind stress seasonal cycle, largely in sync with the surface heat flux cycle, is also a source of SM processes. SM restratification fluxes show distinctive signatures corresponding to both surface cooling and wind stress. The winter peak in SM activity in the Coral Sea is not in phase with the summer dominance of the mesoscale eddy kinetic energy in the region, implying that local surface layer forcing effects are more important for SM gene...