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