This issue documents work to implement a simple first order finite volume scheme, against analytic solutions. This is an academic exercise to further my understanding of the implementation of finite volume schemes.

The issue will follow work in E. Toro, validating the scheme against simple analytic problems - either Burger's equation or the Sod problems, for which the analytic solution has been generated in #1.

This issue covers the implementation of particle motion solver for analytic uniform field solutions. These simulations will be used to validate very simple PIC schemes in the future.

Currently the code requires each simulation to be generated independently of others, leading to code repetition etc. This issue aims to improve the code with a Controller class to improve modularity. The flux calculator may also be separated from the boundary conditions although this may be left for a later stage.

This problem type is being written to ensure the reflecting boundary conditions are working properly, and also to further test the Riemann solver implementation.

This issue covers changes to the code in order to implement random choice on the shock tube 1D problems.

This issue documents to the finite volume fluid solver to add the capability to model multiple fluids in a single simulation.

This involves recording the density of each individual fluid component, as well as its material properties, in order to combine ideal gas equations of state for a particular mixture in a single cell.

The expected solution generated by this issue will do nothing to track the interface.

This issue documents the generation of a random fluid simulation to test the finite volume codes. In theory the simulations should evolve into a uniform flow field if left for sufficient time, which is what I aim to test.

This issue documents the implementation of BEM theory for the generation of power coefficient curves for an arbitrary airfoil.

This issue documents the implementation of the HLLC approximate Riemann solver for ideal gases.

This issue aims to implement some form of WENO scheme for an HLLC finite volume method. The aim is to further my understanding of the state of the art in higher order finite volume methods.

The architecture could do with some improvements to help modularity of different aspects of the simulations. This issue aims to resolve these issues and then also run Noh in 2D to test 2D simulations further.

This issue documents changes made to the Riemann solver code to handle a second order method for better accuracy.

This issue documents the implementation of a simple Lax Friedrichs scheme for Euler's equations. The aim is simply to go back to the basic methods and get a better understanding of artificial viscosity.

A underlying grid to store the degrees of freedom on the DG grid is written.

This issue documents the cleaning up of the directory structure to make sure the work in this repository is cleaner.

This issue documents the implementation of Boris method for analysing the motion of particles in arbitrary magnetic and electric fields.

This issue is an architectural change to the repository, in which I'm going to make the riemann solver directory more general, so as to accommodate the new issues I'm working on in non-riemann solver methods within the same framework.

This issue documents the re-implementation of ESM1, a 1D PIC code for resolving particle motion in electromagnetic fields, taken from C. K. Birdsall - Plasma Physics via Computer Simulation, as a learning exercise.

This issue documents the implementation of a FORCE method for solving Euler's equation. This centred scheme is considered to be very simple, but is implemented so that I can go back to the basic building blocks of some more complex methods, which might rival Riemann solvers.

This issue is a simple stepping stone to implementing similar artificial viscosity schemes like these for the full Euler's equations.

This issue documents the use of boost python wrapping to make use of the BEM solver from a python interface. This will allow easier use of the code to explore what the power coefficient curves I am generating look like.

Currently the VTK writer I am using is broken. The implementation was already quite basic, and could do with some overhaul. For that reason, I should spend some time writing the data stucture properly, including the mesh, and making sure the size of the actualy mesh is representative of the sim etc. etc.

This issue documents the implementation.

This is an experiment to see if such an approach can even be done. The aim would be to recover a HLLC riemann solver using machine learning methods on data generated from an exact riemann solver. I already have the exact method to generate the necessary data, but need to figure out a way to correctly fit the state.

The Double Mach is one of the most complex 2D 1 Fluid Simulations. This issue documents the generation of a double mach simulation, commenting on the results obtained. Currently, I am only going to run this simulation with Godunov. It's already established that this won't resolve the complex flow structures, which requires a MUSCL implementation. Given that I'm working in python and that would probably be a bit too expensive computationally, this issue will only try to use Godunov.

I have implemented previously a C++ version of the Riemann solver code I already have in python in this repository. This issue documents the addition of the serial C code to this repository. In this issue, the old, rusty code needs to be made to compile and output results from VTK as it used to.

This issue covers the generation of a idealised wake rotation model for a wind turbine. This is being written to get a better understanding of the idealised aerodynamics of wind turbines before moving onto BEM theory.

Write tests for gaussian quadrature rules to show order of accuracy for different basis functions in DG grid.

This issue documents the work I undertook to write my own Riemann Solver based on the work in Toro. This was an academic exercise to help me further my understanding of the initial value problem and how it can be solved in hyperbolic differential equations and applied to Euler's equations to solve CFD problems.

Much of the work follows that in Chapter 4 of the book Riemann solvers and Numerical Methods for Fluid Dynamics by E. Toro

This issue documents the addition of a simple profiling tool to the Riemann solver code, to see if I can optimise the python version a little, by using cython or numba.

This numerical method aims to solve for the density behind a shock, given the equation of state, the state ahead of the shock, and the pressure behind the shock.

This issue documents the implementation of a 1D electrostatic field solver. The implementation is taken from C. K. Birdsall - Plasma Physics via Computer Simulation, and is used as a learning exercise to better understand simple PIC codes.