On 6/6/21 7:05 AM, TGCMr. Anonymous wrote:
From figure b, we can see the values on the refined grid are very lower than that before refine. I am thinking, the distribution on refined grid should look same as before, whatever the grid is refined or not. I was confused that the distribution after the refine (figure b) is quite different with that before (figure a). In addition, when I assembled the source term from local cell to global vector, the globe vector of r.h.s has some problem.It looks like the values on some nodes are accounted more than one time.

The terms on the right hand side are of the form

  F_i  =  \int_\Omega  phi_i(x)  f(x)  dx

but because phi_i only lives on the cells adjacent to node i, what you really have is

  F_i  =  \int_{cells around node i}  phi_i(x)  f(x)  dx

If i is a node on a refined cell, then the area of the cells around node i is quite small, so you have

  F_i  =  \int_{small area}  phi_i(x)  f(x)  dx

whereas if i is on a coarse cell, then the area of the cells around node i is quite large, so you have

  F_i  =  \int_{large area}  phi_i(x)  f(x)  dx

This explains why the values of F_i scale with the size/area/volume of cells around node i.

Best
 W>

--
------------------------------------------------------------------------
Wolfgang Bangerth          email:                 [email protected]
                           www: http://www.math.colostate.edu/~bangerth/

--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see 
https://groups.google.com/d/forum/dealii?hl=en
--- You received this message because you are subscribed to the Google Groups "deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to [email protected].
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/8e6b9c3f-5a69-6c30-8d5a-b78e74111913%40colostate.edu.

Reply via email to