Hi Eric and John,<br>Thanks for the speedy replies. What I have plotted here is not normalized by the surface area, it is just the total heat flux, so they should be equal in steady state. The ratio between the two number is not anything obvious like the surface area ratios or anything like that (at least so far as I can tell). I think it is fine that Aspect outputs the total flux rather than flux per unit area, by the way.<br>
Best,<br>Ian<br><br><div class="gmail_extra"><br><br><div class="gmail_quote">On Thu, Nov 15, 2012 at 8:43 AM, <span dir="ltr"><<a href="mailto:aspect-devel-request@geodynamics.org" target="_blank">aspect-devel-request@geodynamics.org</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Send Aspect-devel mailing list submissions to<br>
<a href="mailto:aspect-devel@geodynamics.org">aspect-devel@geodynamics.org</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
<a href="http://geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel" target="_blank">http://geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel</a><br>
or, via email, send a message with subject or body 'help' to<br>
<a href="mailto:aspect-devel-request@geodynamics.org">aspect-devel-request@geodynamics.org</a><br>
<br>
You can reach the person managing the list at<br>
<a href="mailto:aspect-devel-owner@geodynamics.org">aspect-devel-owner@geodynamics.org</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of Aspect-devel digest..."<br>
<br>
<br>
Today's Topics:<br>
<br>
1. REMINDER - Webinar TODAY - Timo Heister on Aspect<br>
(Lorraine Hwang)<br>
2. Re: Question regarding heat fluxes (Eric Heien)<br>
3. Re: Question regarding heat fluxes (John Naliboff)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Thu, 15 Nov 2012 08:38:00 -0800<br>
From: "Lorraine Hwang" <<a href="mailto:ljhwang@ucdavis.edu">ljhwang@ucdavis.edu</a>><br>
Subject: [aspect-devel] REMINDER - Webinar TODAY - Timo Heister on<br>
Aspect<br>
To: <<a href="mailto:aspect-devel@geodynamics.org">aspect-devel@geodynamics.org</a>>, <<a href="mailto:cig-mc@geodynamics.org">cig-mc@geodynamics.org</a>><br>
Message-ID: <015001cdc34f$938dd130$baa97390$@<a href="http://ucdavis.edu" target="_blank">ucdavis.edu</a>><br>
Content-Type: text/plain; charset="us-ascii"<br>
<br>
Dear Colleagues,<br>
<br>
Don't forget this month's CIG Webinar *** TODAY *** Thursday November 15,<br>
2012 @ 2pm PT featuring Timo Heister from TAMU. Timo's talk, "Modern<br>
Numerical Methods for Modeling Convection in the Earth's Mantle" builds on<br>
the previous seminar applying the deal.II library to ASPECT to model<br>
convection in the earth's mantle. In this webinar, Timo will step through<br>
the many ingredients necessary for a modern simulation code: time<br>
discretization, nonlinear methods, linear solvers, stabilization,<br>
adaptivity, and parallelization. Finally, the effectiveness is shown in<br>
several numerical examples.<br>
<br>
You may connect to the meeting through the link provided below or through<br>
the home page at <a href="http://geodynamics.org" target="_blank">geodynamics.org</a>. More complete instructions can also be<br>
found at:<br>
<a href="http://www.geodynamics.org/cig/community/Webinar/webinar" target="_blank">http://www.geodynamics.org/cig/community/Webinar/webinar</a><br>
<br>
Remember to test your connection before the beginning the webinar.<br>
<br>
All webinars will be recorded for later viewing at the above url and on our<br>
YouTube Channel, CIG Geodynamics<br>
(<a href="http://www.youtube.com/user/cigeodynamics" target="_blank">http://www.youtube.com/user/cigeodynamics</a>). Navigate to either to view<br>
last month's webinar by Wolfgang Bangerth, "Using Existing Libraries to<br>
Improve and Solve Computational Problems"<br>
<br>
Please forward to friends and colleagues.<br>
<br>
Regards,<br>
-Lorraine<br>
<br>
****************************************************************************<br>
Please join me in an Adobe Connect Meeting.<br>
<br>
Meeting Name: Webinar: Timo Heister<br>
Summary: Thursday November 15, 2012<br>
Modern Numerical Methods for Modeling Convection in the Earth's Mantle<br>
Timo Heister, Ph.D.<br>
Texas A&M University<br>
Invited By: Lorraine Hwang (<a href="mailto:ljhwang@ucdavis.edu">ljhwang@ucdavis.edu</a>)<br>
When: 11/15/2012 2:00 PM - 3:00 PM<br>
Time Zone: (GMT-08:00) Pacific Time (US and Canada); Tijuana<br>
<br>
To join the meeting:<br>
<a href="http://uc-d.adobeconnect.com/r28i3av93ti/" target="_blank">http://uc-d.adobeconnect.com/r28i3av93ti/</a><br>
<br>
----------------<br>
If you have never attended an Adobe Connect meeting before:<br>
<br>
Test your connection:<br>
<a href="http://uc-d.adobeconnect.com/common/help/en/support/meeting_test.htm" target="_blank">http://uc-d.adobeconnect.com/common/help/en/support/meeting_test.htm</a><br>
<br>
Get a quick overview: <a href="http://www.adobe.com/go/connectpro_overview" target="_blank">http://www.adobe.com/go/connectpro_overview</a><br>
<br>
Adobe, the Adobe logo, Acrobat and Adobe Connect are either registered<br>
trademarks or trademarks of Adobe Systems Incorporated in the United States<br>
and/or other countries.<br>
<br>
<br>
<br>
<br>
<br>
************************************<br>
<br>
Lorraine Hwang, Ph.D.<br>
<br>
Associate Director<br>
<br>
Computational Infrastructure for Geodynamics<br>
<br>
2215A EPS<br>
<br>
One Shields Avenue<br>
<br>
Davis, CA 95616<br>
<br>
w: <a href="tel:530.752.3656" value="+15307523656">530.752.3656</a><br>
<br>
e: <mailto:<a href="mailto:ljhwang@ucdavis.edu">ljhwang@ucdavis.edu</a>> <a href="mailto:ljhwang@ucdavis.edu">ljhwang@ucdavis.edu</a><br>
<br>
<<a href="http://www.geodynamics.org/cig" target="_blank">http://www.geodynamics.org/cig</a>> <a href="http://www.geodynamics.org/cig" target="_blank">http://www.geodynamics.org/cig</a><br>
<br>
Conserve P Print Sparingly<br>
<br>
<br>
<br>
<br>
<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <a href="http://geodynamics.org/pipermail/aspect-devel/attachments/20121115/721558cd/attachment-0001.htm" target="_blank">http://geodynamics.org/pipermail/aspect-devel/attachments/20121115/721558cd/attachment-0001.htm</a><br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Thu, 15 Nov 2012 09:38:30 -0700<br>
From: Eric Heien <<a href="mailto:emheien@ucdavis.edu">emheien@ucdavis.edu</a>><br>
Subject: Re: [aspect-devel] Question regarding heat fluxes<br>
To: <a href="mailto:aspect-devel@geodynamics.org">aspect-devel@geodynamics.org</a><br>
Message-ID: <<a href="mailto:6465B414-1C10-4D44-8479-757AF0307A3C@ucdavis.edu">6465B414-1C10-4D44-8479-757AF0307A3C@ucdavis.edu</a>><br>
Content-Type: text/plain; charset=iso-8859-1<br>
<br>
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>
> Hi all,<br>
><br>
> I have been doing some tests with aspect involving heat fluxes out of the<br>
> top and bottom boundaries, and am getting some behavior I don't<br>
> understand.<br>
><br>
> When the convection reaches steady state, the heat fluxes out of the top<br>
> and bottom should be equal and opposite, otherwise, of course, it wouldn't<br>
> be steady. However, if I look at a simple convection model in with 2d<br>
> shells, I find that the flux out the top is considerably lower than the<br>
> flux out of the bottom at steady state. I am attaching a prm file and a<br>
> plot of heat fluxes as demonstration. The qualitative change at ~10 Gyr is<br>
> due to a transition from organized to chaotic convection. The basic heat<br>
> flux story is unchanged, though.<br>
><br>
> Looking at the code for heat flux statistics, nothing seems obviously<br>
> wrong. I should note that doing the same test for a box geometry where the<br>
> top and bottom have the same area produces the expected results.<br>
><br>
> Am I missing something about how this should work? Any insights as to what<br>
> is going on would be appreciated.<br>
> Best,<br>
> Ian<br>
><br>
> PS All calculations are done with r1353<br>
> <test_flux.prm><<a href="http://flux.ps" target="_blank">flux.ps</a>>_______________________________________________<br>
> Aspect-devel mailing list<br>
> <a href="mailto:Aspect-devel@geodynamics.org">Aspect-devel@geodynamics.org</a><br>
> <a href="http://geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel" target="_blank">http://geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel</a><br>
<br>
<br>
<br>
------------------------------<br>
<br>
Message: 3<br>
Date: Thu, 15 Nov 2012 11:43:39 -0500<br>
From: John Naliboff <<a href="mailto:jbnaliboff@ucdavis.edu">jbnaliboff@ucdavis.edu</a>><br>
Subject: Re: [aspect-devel] Question regarding heat fluxes<br>
To: <a href="mailto:aspect-devel@geodynamics.org">aspect-devel@geodynamics.org</a><br>
Message-ID: <<a href="mailto:7F5868C7-29DD-455D-8DEE-2C5FEA3C5497@ucdavis.edu">7F5868C7-29DD-455D-8DEE-2C5FEA3C5497@ucdavis.edu</a>><br>
Content-Type: text/plain; charset="us-ascii"<br>
<br>
Hi Ian,<br>
<br>
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.<br>
<br>
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.<br>
<br>
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.<br>
<br>
Cheers,<br>
John<br>
<br>
<br>
<br>
<br>
On Nov 15, 2012, at 11:38 AM, Eric Heien wrote:<br>
<br>
> 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>
>> Hi all,<br>
>><br>
>> I have been doing some tests with aspect involving heat fluxes out of the<br>
>> top and bottom boundaries, and am getting some behavior I don't<br>
>> understand.<br>
>><br>
>> When the convection reaches steady state, the heat fluxes out of the top<br>
>> and bottom should be equal and opposite, otherwise, of course, it wouldn't<br>
>> be steady. However, if I look at a simple convection model in with 2d<br>
>> shells, I find that the flux out the top is considerably lower than the<br>
>> flux out of the bottom at steady state. I am attaching a prm file and a<br>
>> plot of heat fluxes as demonstration. The qualitative change at ~10 Gyr is<br>
>> due to a transition from organized to chaotic convection. The basic heat<br>
>> flux story is unchanged, though.<br>
>><br>
>> Looking at the code for heat flux statistics, nothing seems obviously<br>
>> wrong. I should note that doing the same test for a box geometry where the<br>
>> top and bottom have the same area produces the expected results.<br>
>><br>
>> Am I missing something about how this should work? Any insights as to what<br>
>> is going on would be appreciated.<br>
>> Best,<br>
>> Ian<br>
>><br>
>> PS All calculations are done with r1353<br>
>> <test_flux.prm><<a href="http://flux.ps" target="_blank">flux.ps</a>>_______________________________________________<br>
>> Aspect-devel mailing list<br>
>> <a href="mailto:Aspect-devel@geodynamics.org">Aspect-devel@geodynamics.org</a><br>
>> <a href="http://geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel" target="_blank">http://geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel</a><br>
><br>
> _______________________________________________<br>
> Aspect-devel mailing list<br>
> <a href="mailto:Aspect-devel@geodynamics.org">Aspect-devel@geodynamics.org</a><br>
> <a href="http://geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel" target="_blank">http://geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel</a><br>
<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <a href="http://geodynamics.org/pipermail/aspect-devel/attachments/20121115/e6dc0410/attachment.htm" target="_blank">http://geodynamics.org/pipermail/aspect-devel/attachments/20121115/e6dc0410/attachment.htm</a><br>
<br>
------------------------------<br>
<br>
_______________________________________________<br>
Aspect-devel mailing list<br>
<a href="mailto:Aspect-devel@geodynamics.org">Aspect-devel@geodynamics.org</a><br>
<a href="http://geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel" target="_blank">http://geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel</a><br>
<br>
<br>
End of Aspect-devel Digest, Vol 12, Issue 4<br>
*******************************************<br>
</blockquote></div><br></div>