Homepage - Theoreticalibus - Software - Experimentation - Junk - About me



∇ Hyberbolic equations, hypersonics and kinetic theory

Hyper2D - [Fortran90/03, Jan 2022] - Hyper2D is a 1D or 2D finite-volume solver for hyperbolic PDEs and hypersonic flows. I'm now in the process of uploading it to github as Open Source (a wiki here!).

This is my main software contribution from the PhD and postdoc. What is cool about it? Simplicity.
Seriously, it is simple and easy to learn. This makes it ideal for anyone wishing to learn CFD from scratch. Also, it is very flexible, and you can easily implement your own sets of equations and solve them, as long as they are hyperbolic. Extending the method to parabolic equations will be quite trivial.

There are different versions of Hyper2D, with increasing complexity.

  • A basic version to show the flow of the program;
  • Some better implentations, with van Leer's MUSCL reconstruction (2nd order in space), or 5th order WENO for instance;
  • A CUDA Fortran version for running serious cases, with 1st order or 3rd order low-storage TVD Runge-Kutta explicit time integrator.
What does it have to do with my PhD?
Well, I had to solve some strange hyperbolic system of equations (moment equations) for non-equilibrium hypersonic flows & plasmas. Hyper2D makes it easy to implement your own equations and is not limited to traditional hydrodynamics!


LARSEN - [C++, May 2019] - LAgrangian Reactor for StrEams in Nonequilibrium, is a Lagrangian solver able to recompute the chemical composition of a multi-temperature reacting mixture, from a low-fidelity baseline simulation. Result of my MSc thesis and Research Master project work at VKI. Only the single-streamline mode is available in this repository. For the multi-streamline mode, drop me an email.



FVKin - [C++, November 2018] - This is a minimal C++ Finite Volume solver for the 1D1V collisionless Vlasov/Boltzmann equation with forcing and source terms. External force field (assumed to be electrostatic in nature) and chemical source are imported from a file. Small and undocumented, the code should be self-explanatory.

The solver returns the distribution function "f", depicted here at the right for a testcase. The scheme is second order in space and first order explicit in time.


∇ Parallel Computing

Linux Clustering Step-by-Step plus something extra - [PDF document, October 2015] - This is a step by step tutorial on making a Beowulf cluster out of computers using Linux and Open MPI. Part I: basic cluster tutorial - Part II: diskless cluster tutorial. Actually, that's not just a tutorial since I've tried to fill it with various appendixes: from the bootstrap process to a brief history of UNIX shells. Hope you enjoy!



∇ FEM

CheeseFEM - [Octave, July 2014] - Here's a little FEM code I wrote under Octave to analyze the thermal transitory into a piece of Emmental cheese! The goal is to solve parabolic (no convection/reaction) equations into a heterogeneous domain. The code is composed by a mesh generator, a preprocessor and the solver / postprocessor. CheeseFEM presentation (pdf).

In the picture, a simulation of a material with very-little-conductive particles inside


NSSSS - Navier Stokes Stability and Structural Sensitivity - [MATLAB and FreeFem++, Feb 2015] - This is a collection of Matlab and FreeFem++ scripts to perform stability and structural sensitivity analysis of an incompressible fluid flow. Click here for a pdf document with some theory and many colored plots.

In the picture, the structural sensitivity parameter for a flow past a triangular bluff-body


∇ Aerodynamics

VorDyn2.0 - [MATLAB, May 2014] - VorDyn2.0 is a vortex lattice method for computation of forces, moments and stability derivatives of aerodynamic surfaces. You can create the geometry by editing an input file, export it / run the simulation etc etc. VorDyn2 is basically a collection of MATLAB scripts. Click on the title to download a .zip or click here for the manual: VorDyn2 Manual. You can click on the image for bigger resolution.

VorDyn1.0 - [MATLAB, Feb 2014] - VorDyn1.0 is simply the previous version of VorDyn2, supporting only aircraft-like geometries. The advantage of this old version is that is gives an extimation of the stability derivative Cl_alphadot and Cm_alphadot, not supported anymorein VorDyn2. VorDyn1.0 Manual.



∇ Electromagnetism

Qsci - [MATLAB, April 2015] - Qsci is a MATLAB script that plots the electrostatic field generated by a set of conductors. You can assemble some geometries, impose a Voltage or total Charge on conductors and Qsci will plot surface charge density, the electrostatic potential and compute Capacitances. With a single click on the name you can download the marvellous script in compressed .tar.gz format, or click here to download a simple manual with colored pictures. Qsci Manual. You can click on the image for bigger resolution.