[aspect-devel] Help needed with setting up an EBA benchmark

[Katrina Arredondo]

Magali and I agree with Ian. I was confusing myself because I didn't read
the pdf very carefully. Good luck with the benchmark.


On Fri, Sep 5, 2014 at 2:40 PM, Wolfgang Bangerth <bangerth at tamu.edu> wrote:

>
> Ian,
>
>  Excuse the overly verbose answer, I'm talking this out as much for me as
>> anyone.  So, thermodynamically, the adiabatic heating term is sensitive
>> to the
>> total temperature rather than differences in temperature (as in, for
>> instance,
>> the buoyancy).  I see that at the top of page three they have defined T_0
>> to
>> be the dimensional surface temperature over the temperature scale, \Delta
>> T.
>>   In the anelastic case there would be an adiabatic profile starting from
>> T_0,
>> but in this case, because it is incompressible, the reference temperature
>> is
>> T_0 all the way down.
>>
>> The way it shakes out in these equations, T_0 becomes a free parameter
>> that
>> can increase or decrease the effect of adiabatic heating on the model.
>> For
>> the benchmark case, they use T_0 = 273/3000, like you said, but I could
>> put in
>> anything, really.  All this is a long way of saying that in your parameter
>> file you have set the boundary conditions for temperature to go from zero
>> to
>> one at the bottom.  To my eyes, it looks like you can reproduce the
>> conditions
>> of the benchmark by simply changing the temperature boundary conditions to
>> 0.091 - 1.091.
>>
>
> Excellent suggestion. Now I just need to figure out why I get exactly the
> same results as if I didn't have friction and adiabatic heating switched
> off...
>
> Thanks for talking me through this!
>
> Best
>  Wolfgang
>
> --
> ------------------------------------------------------------------------
> Wolfgang Bangerth               email:            bangerth at math.tamu.edu
>                                 www: http://www.math.tamu.edu/~bangerth/
>
> _______________________________________________
> Aspect-devel mailing list
> Aspect-devel at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel
>



-- 

Katrina Arredondo
Ph.D. Candidate
Geodynamics/Subduction Zones
University of California, Davis
Department of Geology
One Shield Avenue, Davis CA 95616
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20140905/e90dd707/attachment.html>