<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Hi Wolfgang,<div><br></div><div>Happy Thanksgiving!</div><div><br></div><div>That was Ian who sent that input file over. The spherical shell models I tested actually all crashed a little bit after I last emailed you, with all of them reporting very high (> 1e20) residuals at there final time step. The crash occurred over a range (~ 10^8 - 10^4).</div><div><br></div><div>I'm a bit surprised at the high resolutions, considering my mesh setup: 4 initial global refinements, 3 initial adaptive refinements, refinement fraction of 0.3 and coarsening fraction of 0.05. I'll try modifying the mesh to see if that helps with convergence.</div><div><br></div><div>In terms of resolving basal verse surface heat flux, would it be difficult to include a parameter that modifies the spatial resolution as a function of depth, in addition to temperature, cp, etc?</div><div><br></div><div><div apple-content-edited="true">
<span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; font-size: medium; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; font-size: medium; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; font-size: medium; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; font-size: medium; "><div>Cheers,</div><div>John</div><div><br></div></span><br class="Apple-interchange-newline"></span><br class="Apple-interchange-newline"></div></span></span>
</div>
<br><div><div>On Nov 19, 2012, at 3:27 PM, Wolfgang Bangerth wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div><br>All,<br>I know it's terrible -- I don't recall who exactly sent me this input file (attached) and whether I modified it, but in any case I ran it over the weekend to see what happens to the boundary fluxes in the shell geometry. Attached is what I get from simply running this on 16 processors for a *very* long time.<br><br>The two heat fluxes never reach a steady state -- rather, the jump up and down every time a plume detaches from the inner or outer boundaries. Looking at the solution picture (also attached), I'm not surprised that the heat fluxes don't coincide: it looks like the boundary layers are less than one cell width, and so we can not expect to resolve them accurately. Rather, the heat flux given by these numerical approximations is determined by the smallest length scale that can be resolved by the mesh -- and because this length scale is different for the inner and outer boundaries, the heat fluxes is different. (Of course, had we chosen a mesh where the radial resolution is the same at the two boundaries, we may have gotten results where the two heat fluxes match -- but are both completely wrong.)<br><br>Best<br> W.<br><br>-- <br>------------------------------------------------------------------------<br>Wolfgang Bangerth email: <a href="mailto:bangerth@math.tamu.edu">bangerth@math.tamu.edu</a><br> www: <a href="http://www.math.tamu.edu/~bangerth/">http://www.math.tamu.edu/~bangerth/</a><br><br><span><test_flux.prm></span><span><heatflux.png></span><span><visit0001.png></span></div></blockquote></div><br></div></body></html>