Computational Fluid Dynamics - Simulation of Dye Mixing in a Fluid using Matlab
This project was a part of the Graduate Coursework - MAE 561 - Computational Fluid Dynamics. It included
- Simulating the 2D flow described by the Navier Stokes Equations using the Fractional Step Method on a Staggered Mesh.
- The code utilizes Adams-Bashforth Crank-Nicholson method was used for the predictor step and the V Cycle Multi Grid method for the pressure Poisson Equation.
- The dye mixing part was coded using the Second order Lax Wendroff method. The dye was described using the Viscosity and Schmidt numbers.
- Performed a Grid Convergence Analysis by mesh refinement to estimate velocity of the fluid within 0.25% accuracy.
- Visualized the flow using velocity vectors, vorticity, kinetic energy, segregation intensity plots with time and made a small animation of the whole mixing using Matlab.
The fluid is modeled as a rectangular plate with some source terms kept at the walls. The dye were kept as two separate circular blobs in the fluid at t= 0. The mixing is shown in the simulations for a period of 10 seconds. The contour indicates the segregation intensity with 1 being the highest.
A simulation of the Dye mixing for Re = 66 and Sc = 0.1
A slower simulation for Re = 66 and Sc = 0.01