Hi David,

Ah got it - thanks a lot! I was attempting to
use write_discontinuous_equation_systems() and your help totally solved the
problem. Appreciate it!

Best,
Shawn

On Mon, Jan 14, 2019 at 6:46 PM David Knezevic <[email protected]>
wrote:

> On Mon, Jan 14, 2019 at 9:40 PM Yuxiang Wang <[email protected]> wrote:
>
>> Dear all,
>>
>> I am using the L2_LAGRANGE elements for stress recovery, because the
>> stress
>> are not continuous over the domain. I noticed that my
>> stress_system.solution->size() indicate that I have more than one value
>> for
>> each node and each stress component. This makes sense - I should have as
>> many values as the number of elements connecting to that node.
>>
>> However, when I used write_equation_system() functionality and visualize
>> in
>> ParaView, only one value is shown per node. I suspect that it is an
>> average
>> of all element DOFs on this node but am not sure - could you please help
>> share a hint?
>
>
> I'm not sure what paraview plots in this case, but you can plot a
> discontinuous solution by using write_discontinuous_exodusII, for example.
> If I recall correctly, that will create a new mesh with disconnected
> elements so that paraview will correctly plot it as discontinuous. The
> downside of this is the mesh no longer matches the mesh that you solved on
> (since the elements are disconnected) but that may be fine in your case.
>
> Best,
> David
>


-- 
Yuxiang "Shawn" Wang, PhD
[email protected]
+1 (434) 284-0836

_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users

Reply via email to