fem_hughes VS dolfinx

I was lucky enough to have taken two FEA classes from Hughes in the 80s. I worked with his DLEARN [1] code many times back then. He was a great teacher.

I never had the pleasure of taking classes from Juan Carlos Simo, but he was known to have outstanding classes. His was a very brilliant light & life cut too short by cancer at the young age of 42.

[1]. DLEARN is a linear static and dynamic finite element code written in Fortran. https://github.com/fit087/fem_hughes

One thing that is a pain when learning this stuff is that actually performing the optimization requires a good understanding of the numerical discretization of PDEs. Finite elements are a natural choice because it is very easy to characterize the adjoint with this formulation. There are some good free tools that you can use to actually learn and do some computations yourself. The first is hIPPYlib (paper, code), which is built on top of FEniCS (link), for which there are many good tutorials. Beware trying to install this on Windows though. You will need to work in Docker or in Ubuntu via Windows Linux Subsystem.

If you're ok with coding, fenics is a solid place to start. Also if you're comfortable with coding, openfoam is FVM, rather than FEM, but it can handle solidmechanics.

> FEniCs made FEM so easy

https://fenicsproject.org/

Indeed, was blown away when I saw it for the first time over a decade ago, compared to the convoluted C++ FEM libraries I had seen before that.

Have you looked at FEniCS? Pretty much everything else I'm aware of is probably overkill (e.g., MOOSE in C++, HYPRE's Python bindings, etc.)

FEniCSx is quite good.

You mean Python? For many research tasks it's fine. High level libraries let you define your computation in a minimal amount of code. FEniCS is a great example of this - underneath it compiles the abstracted high level stuff to calls to low-level libraries that do the heavy lifting. For many applications you can just write vectorized code with Numpy that performs well, or use Numba to JIT what you can't vectorize. For some tasks, however, you need interfaces that don't exist in the high level libraries, and that was the case for me.

I just came across FEniCSX. I’m not sure if it’s what you want but here’s the description:

Fenics: Very popular finite element framework “UseTab: Never” https://github.com/FEniCS/dolfinx/blob/main/.clang-format