Fluctuating Hydrodynamics


  • The GPU code fluam developed by collaborators Rafael Delgado-Buscalioni and Florencio Balboa implements compressible and incompressible staggered-grid fluctuating hydrodynamics for binary fluid mixtures and for particle suspensions. For incompressible hydrodynamics only periodic boundary conditions are supported because the code relies on FFT techniques for solving the Stokes equations.
  • The state-of-the-art parallel C++ Immersed Boundary Adaptive Mesh Refinement IBAMR framework developed and maintained by collaborator Boyce Griffith implements a stochastic stress tensor generator that can be used to implement fluctuating hydrodynamics for both binary fluid mixtures and for particle suspensions (presently only for neutrally buyoant particles). Full support for any combination of no-slip, slip and periodic boundary conditions is provided. You will need to download the ibamr-dev branch. Here is a simple example code  in which passive tracers diffuse as they are advected by random velocity fluctuations.
  • Parallel finite volume solvers for the compressible, incompressible, and low Mach number fluctuating hydrodynamic equations for fluid mixtures have been implemented in the Fortran binding of the BoxLib framework developed by the group of collaborator John Bell at LBNL. We are in the process of developing a module for staggered-grid low Mach hydrodynamics that will be released to the public once ready.

Hard-Particle packing Codes



  1. Here are some old but still very efficient precompiled executables for generating packings of spheres, ellipsoids and super ellipsoids. These were developed during my Ph.D. research and you should consult my Ph.D. thesis if you plan to use them.
  2. If you wish to modify the codes yourself you should instead use the simpler but still reasonably efficient C++ sphere packing code.