Closed samuelpmishLLNL closed 1 month ago
I am seeing a numerical inconsistency in parallel:
I'll investigate tomorrow
edit: this has been fixed
I believe the parallel issue is fixed now, and this is ready for review.
Is the L2 field equivalent to the volume averaging trick I've been using for a while in LiDO? That is, integrating the stress in each element and dividing by element volume
Is the L2 field equivalent to the volume averaging trick I've been using for a while in LiDO? That is, integrating the stress in each element and dividing by element volume
Yes, the volume averaging trick is equivalent to specifying L2<0>
for the approximation space. However, now you can try higher order fits too, if you like (it might also work for fitting vector-valued quantities as well, I can't remember). In theory, the same feature also supports fitting functions to Hcurl
spaces as well, although I don't have a test for that.
Is the L2 field equivalent to the volume averaging trick I've been using for a while in LiDO? That is, integrating the stress in each element and dividing by element volume
Yes, the volume averaging trick is equivalent to specifying
L2<0>
for the approximation space. However, now you can try higher order fits too, if you like (it might also work for fitting vector-valued quantities as well, I can't remember). In theory, the same feature also supports fitting functions toHcurl
spaces as well, although I don't have a test for that.
Awesome, once this propagates to LiDO I'll update our serac visualization operator to use it so we can plot smoother fields. Thanks @samuelpmish!
Finally revisiting this PR.
Here's an example of how to fit a field to the J2 stress resulting from a displacement field (in 2D), the user writes:
where
input_space
is the function space for the displacement, andoutput_space
is the function space of the fitted approximation, andu
is the displacement field to interpolate for the supplied q-function.Here are examples of the J2 stress on a fracture specimen in plane strain, being approximated by both a piecewise constant field (
L2<0>
) and a continuous, piecewise linear field (H1<1>
):And here are some closeups of the
L2<0>
andH1<1>
approximations of the J2 stress near the crack tip: