jtoumey/2d_fvm

Hi there 👋

2d_fvm's Issues

Document the Code

Add a LaTeX document that contains instructions for use and explains the theory of the code.

Develop General Method for Handling Boundaries

Currently, the code solves boundaries separately from interior cells (within each function). This seems inefficient, and maybe some OOP techniques such as polymorphism or operator overloading will help.

Implement Iterative Solvers

The current solution method is a direct method -- the Thomas algorithm for tri-diagonal matrices. This works by solving along N-S lines and treating the W and E coefficients explicitly.

Add the option to choose iterative solvers (Jacobi, Gauss-Seidel, etc.) as an exercise.

Additionally, I can modularize all of the components of the Thomas algorithm procedure (e.g., bounds calculation, update_explicit, and thomas).

Make FVM Coefficients Directional

Currently, the code only accepts flow from W - > E and from S -> N. Add functionality to handle any direction.

Increase Precision of Output Data

Currently, the write_results function includes only two significant figures beyond the decimal point. This is a huge problem because for fine grids (e.g., 1000 x 1000), just two significant figures cannot represent the difference in grid spacing.

This is likely the reason that plotting the diagonal values of phi results in discontinuities.

Improve read_input Function

• Test if the input file exists
• Make the input file a general file (not only grid data)
• Allow the inputs to be in any order