For my PhD I have been actively involved in the development, extension, optimization and parallelization of software for simulating turbulent
phenomena. On this page I provide a summary of the research codes I have been involved with during my time at UCSD.
- Explicit RK3 for free shear flows
- Mixed RK3-ADI with the immersed boundary method for wall bounded flows with complex geometry
- Parallel geometric multigrid pressure solver
The research modules described below were created by members of the CFD lab under the supervision of professor Sutanu Sarkar. They are under
continual development and the descriptions below may not reflect the most recent versions available.
Explicit RK3 for free shear flows
- Designed for study of free shear flows (wakes, jets, shear layers)
- 3D Incompressible Navier-Stokes equations with the Boussinesq approximation
- Finite volume formulation
- Direct numerical simulation
- Staggered grid formulation: velocity at cell faces, pressure and density at cell centers
- Low storage third order Runge-Kutta for temporal advancement
- 2nd order centered differences for spatial terms
- Multigrid pressure solver
- Has been used on up to 2 billion grid points and 2048 processors
- Method of manufactured solutions available for code verification
- 3D domain decomposition, parallelization with MPICH-II
- Written in Fortran 90
Mixed RK3-ADI with the immersed boundary method for wall bounded flows with complex geometry
- Designed for study of flows with complex geometry
- 3D Incompressible Navier-Stokes equations with the Boussinesq approximation
- Finite volume formulation
- Direct numerical simulation and large eddy simulation
- LES models: Standard Smagorinsky Model, Dynamic Smagorinsky model with pathline averaging of Meneveau et. al., JFM 1996
- Sharp interface immersed boundary method implementation of Roman et. al., Computers & Fluids 2009
- Collocated grid formulation with pressure-correction algorithm
- Low storage third order Runge-Kutta for temporal advancement of convective terms
- Second order alternating direction implicit method for time advancement of viscous terms
- Efficient pipelined Thomas algorithm for ADI implementation
- 2nd order centered differences for spatial derivatives
- 5th order biased upwinding based on Rai & Moin JCP 1991 to handle thin shear layers if necessary
- Multigrid pressure solver
- Method of manufactured solutions available for code verification
- 3D domain decomposition, parallelization with MPICH-II
- Written in Fortran 90
Parallel geometric multigrid pressure solver
- V cycle and full multigrid (FMG) cycling
- For strongly stretched grids: Flexible semi-coarsened cycling based on the MG-S variant of Piquet & Vasseur, Numerical Algorithms, 2000
- Suitable for uniform and stretched grids
- Designed for all boundary conditions: Dirichlet, Neumann, and periodic cases implemented
- Cell centered coarsening
- Full volume weighting for restriction
- Full volume weighting (trilinear interpolation) for prolongation
- Smoothers: Lexicographic and Red-Black Gauss-Seidel. Block point and block line versions available. Successive over relaxation (SOR)
available for all smoothers.
- Collection of coarse sub-problem to master processor for coarse grid solution in parallel
- Method of manufactured solutions available for code verification
- 3D domain decomposition, parallelization with MPICH-II
- Written in Fortran 90
© 2012 Matt de Stadler