[aspect-devel] Internal heating in aspect

[Scott King]

Yes.  If you take stresses or any gradient quantities at the integration points and project to the surface it is much better.  

Scott

Sent from my iPhone

> On Aug 27, 2018, at 4:49 PM, Wolfgang Bangerth <bangerth at colostate.edu> wrote:
> 
>> On 08/25/2018 02:17 PM, Ian Rose wrote:
>> I found that using CBF in calculating the dynamic topography significantly improved the convergence of the solution (quadratic rather than linear, if memory serves). The mathematics for computing the heat flux is essentially identical, so a lot of the dynamic topography code could be copied whole-cloth into the heat flux postprocessor, in case anyone wants to give it a shot.
> 
> Jacky and I found similar things with the dynamic topography postprocessor which is another example of computing something along the boundary that contains the gradient of the solution. There, instead of taking the stress at the surface, we take the stress at the center of the top-most cell (if I recall correctly) or some average over the cell, and found that this *substantially* improved the accuracy of the result. Something like this would be easy to implement in the temperature flux postprocessor as well.
> 
> Best
> W.
> 
> -- 
> ------------------------------------------------------------------------
> Wolfgang Bangerth          email:                 bangerth at colostate.edu
>                           www: http://www.math.colostate.edu/~bangerth/
> _______________________________________________
> Aspect-devel mailing list
> Aspect-devel at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel