Brandon Talamini
Brandon Talamini
Use MacPorts.
This PR was spurred by #754. The main change is to add more rigorous tests of the solid mechanics module. These are patch tests, which ensure that the numerical approximation...
The functional-based solid mechanics operator uses a total Lagrangian setting, so it uses the first Piola-Kirchhoff stress in the residual calculation. However, in the kernel [here](https://github.com/LLNL/serac/blob/5aeeea93f029d3b3afcbe30e9ce39253c1ea77e8/src/serac/numerics/functional/integral_utilities.hpp#L523) the stress is contracted...
The request here is to be able to grab the single component of a vector-valued field and treat it as a scalar, say to compute norms. I suppose this would...
Inside of the control flow in the functional paradigm, we currently have the following steps: 1. Interpolate, which goes from L-vector to Q-vector (DOFS to primal fields and their gradients...
There are some core tensor operations that are needed, as well as some that would be nice to have. In my opinion, these are the needs. They're all for 2D...
The reaction force calculation in the solid mechanics module assumes that the accelerations are zero: https://github.com/LLNL/serac/blob/ef413804eda27556cec2b456ac45b44c830b23c1/src/serac/physics/solid_mechanics.hpp#L669-L673 This is fine for quasi-statics, but for dynamics, we must pass the updated accelerations...
We need to include quadrature data in restart files. We can't use the `mfem::QuadratureFunction` interface because this only works for struct-of-arrays, while serac uses array-of-structs. Instead, we should use a...
We're currently checking that all elements can reproduce affine functions (i.e., patch tests). We should go further and verify that higher-order elements can reproduce polynomial solutions up to their order.
@kswartz92 wants to solve a series of problems on the same mesh with different essential boundary conditions (for concreteness, this is for homogenization). It would be helpful if one could...