Re: [Libmesh-users] Question about the nodal solution value outputted via write_equation_systems() from L2_LAGRANGE elements

[Yuxiang W. <yw...@vi...>]

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 <dav...@ak...>
wrote:

> On Mon, Jan 14, 2019 at 9:40 PM Yuxiang Wang <yw...@vi...> 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
yw...@vi...
+1 (434) 284-0836