<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Hi Ian,<div><br></div><div><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; ">That's definitely odd. At first glance I don't see anything in the input file that would cause that kind of difference between the surface and basal heat flow. You do have different velocity boundary conditions for the top and bottom boundaries, but that should not cause that much of a difference in the heat flow.</span></div><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><br></div><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">As Eric mentioned, the scaling of the heat flux from the output file would be different for the top and bottom boundaries. That's the only think I can think of right now as well.</div><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><br class="Apple-interchange-newline"></div><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">I'm also running some 2D shell models right now and I'll check out how the surface and basal heat fluxes compare. I don't think they have reached steady yet, but it will be good to see if they exhibit the same behavior. I'll send some figures over later this afternoon.</div><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><br></div><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Cheers,</div><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">John</div><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><br></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; "><br class="Apple-interchange-newline"></span><br class="Apple-interchange-newline"></div></span></span>
</div>
<br><div><div>On Nov 15, 2012, at 11:38 AM, Eric Heien wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div>Are you normalizing based on the size of the boundary? One thing that is not obvious and that we've had troubles with is that the printed flux is not per unit area, but over the whole boundary area. Maybe this should be changed?<br><br>-Eric<br><br>On Nov 15, 2012, at 9:27 AM, Ian Rose wrote:<br><br><blockquote type="cite">Hi all,<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">I have been doing some tests with aspect involving heat fluxes out of the<br></blockquote><blockquote type="cite">top and bottom boundaries, and am getting some behavior I don't<br></blockquote><blockquote type="cite">understand.<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">When the convection reaches steady state, the heat fluxes out of the top<br></blockquote><blockquote type="cite">and bottom should be equal and opposite, otherwise, of course, it wouldn't<br></blockquote><blockquote type="cite">be steady. However, if I look at a simple convection model in with 2d<br></blockquote><blockquote type="cite">shells, I find that the flux out the top is considerably lower than the<br></blockquote><blockquote type="cite">flux out of the bottom at steady state. I am attaching a prm file and a<br></blockquote><blockquote type="cite">plot of heat fluxes as demonstration. The qualitative change at ~10 Gyr is<br></blockquote><blockquote type="cite">due to a transition from organized to chaotic convection. The basic heat<br></blockquote><blockquote type="cite">flux story is unchanged, though.<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Looking at the code for heat flux statistics, nothing seems obviously<br></blockquote><blockquote type="cite">wrong. I should note that doing the same test for a box geometry where the<br></blockquote><blockquote type="cite">top and bottom have the same area produces the expected results.<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Am I missing something about how this should work? Any insights as to what<br></blockquote><blockquote type="cite">is going on would be appreciated.<br></blockquote><blockquote type="cite">Best,<br></blockquote><blockquote type="cite">Ian<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">PS All calculations are done with r1353<br></blockquote><blockquote type="cite"><test_flux.prm><flux.ps>_______________________________________________<br></blockquote><blockquote type="cite">Aspect-devel mailing list<br></blockquote><blockquote type="cite"><a href="mailto:Aspect-devel@geodynamics.org">Aspect-devel@geodynamics.org</a><br></blockquote><blockquote type="cite"><a href="http://geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel">http://geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel</a><br></blockquote><br>_______________________________________________<br>Aspect-devel mailing list<br><a href="mailto:Aspect-devel@geodynamics.org">Aspect-devel@geodynamics.org</a><br>http://geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel<br></div></blockquote></div><br></div></body></html>