[aspect-devel] Aspect-devel Digest, Vol 12, Issue 3

[John Naliboff]

Hi Wolfgang,

Thanks for the reply.

I've attached two images that show the temperature (whole model) and mesh (upper section) corresponding to the end of an updated heat flux vs time plot (also attached).

The model is indeed  in a similar initial phase as the model posted in your link, so good to know the heat flux discrepancy levels out eventually.

 I'll report back if the model does not transition to the chaotic stage, but I don't think there is any reason to think it will not.

Thanks again,
John
 



On Nov 17, 2012, at 8:17 AM, Wolfgang Bangerth wrote:

> 
> John,
> 
>> I've attached a plot showing the basal and surface heat flow for a
>> spherical shell model I'm currently running.  The Ra number is ~ 7.66e5.
>>  The model has not completely reached steady state yet, but if the
>> trend continues there will be a distinct gap between the surface and
>> basal heat flow.  The .prm file I used to run the model is also attached.
> 
> Can you post a picture of the solution at the time corresponding to the end of the plot you show? If this setup is anything similar to the ones I've made in this geometry then you shouldn't expect to run into any steady-state limit at all. Rather, you get one set of blobs rising at the very beginning and then a long phase where hot material streams up to the surface. The heat flux will then look like in your graph. But it's not a steady state -- it simply moves heat from the bottom into the domain but not out of the top, like in your graph. This takes a long time until the system finally converts to some more chaotic movement. See for example here:
>  http://www.math.tamu.edu/~bangerth/step-32.2d.convection.gif
> So if you're still in the initial phase where nothing happens for a long time, then I wouldn't expect the heat flux at the top and bottom to equal.
> 
> 
>> In terms of resolution, I did 4 initial global refinement steps, 3
>> initial adaptive refinement steps, and the coarsening and refinement
>> fractions are, respectively, 0.05 and 0.3.
>> 
>> I'll try increasing the refinement fraction and initial refinements to
>> see if this reduces the gap between the global and surface heat fluxes.
>> I am bit surprised that the mesh parameters above did not adequately
>> resolve the basal boundary layer, if this is also the issue in my models.
> 
> Can you post a picture of the mesh as well?
> 
> Best
> W.
> 
> -- 
> ------------------------------------------------------------------------
> Wolfgang Bangerth               email:            bangerth at math.tamu.edu
>                                www: http://www.math.tamu.edu/~bangerth/
> 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://geodynamics.org/pipermail/aspect-devel/attachments/20121117/b1301df9/attachment-0004.htm 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: heat_flow_shell_15Nov17.pdf
Type: application/pdf
Size: 432876 bytes
Desc: not available
Url : http://geodynamics.org/pipermail/aspect-devel/attachments/20121117/b1301df9/attachment-0002.pdf 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://geodynamics.org/pipermail/aspect-devel/attachments/20121117/b1301df9/attachment-0005.htm 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: b3_2ds_phi360_temp.pdf
Type: application/pdf
Size: 54127 bytes
Desc: not available
Url : http://geodynamics.org/pipermail/aspect-devel/attachments/20121117/b1301df9/attachment-0003.pdf 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://geodynamics.org/pipermail/aspect-devel/attachments/20121117/b1301df9/attachment-0006.htm 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: b3_2ds_phi360_mesh.png
Type: image/png
Size: 64865 bytes
Desc: not available
Url : http://geodynamics.org/pipermail/aspect-devel/attachments/20121117/b1301df9/attachment-0001.png 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://geodynamics.org/pipermail/aspect-devel/attachments/20121117/b1301df9/attachment-0007.htm