How do instabilities form and evolve in the turbulent waters of the ocean? This project explores the question by numerically solving the 3D incompressible Navier-Stokes equations (with mass equaiton evolved using either density or buoyancy at the user's discretion) under the Boussinesq approximation for stratified fluids using a user customisable choice of either finite spatial differences or Fast Fourier Transform (FFT) spectral methods advanced in time through an RK4 time integrator. Incompressibility is maintained using an FFT based spectral divergence cleaner. Rayleigh-Taylor instabilities and Kevin-Hemmholtz instabilities are simulated. The solver has been benchmarked against the Taylor-Green Vortex with a maximum error of 2% in energy decay rate. (Mar 2025, Year 2)
Github link and PDF coming soon.