[aspect-devel] Cylindrical coordinates

[Jonathan Perry-Houts]

Do you mean the axisymmetric case with a non-zero inner radius (i.e. 
still in (r, z) coordinates)? The perpendicular plane (r,theta) already 
exists as the 2D 'spherical shell', correct?

On 05/09/2018 11:07 PM, Maxwell Rudolph wrote:
> If you decide to pursue this it might make sense to also implement the 2D spherical annulus simultaneously.
> 
>> On May 9, 2018, at 10:16 PM, Jonathan Perry-Houts <jperryh2 at uoregon.edu> wrote:
>>
>> Thanks Wolfgang and Felipe. Sounds like there's no library magic that saves me in this case. Thanks for the advice!
>>
>> -Jonathan
>>
>>> On 05/09/2018 08:11 PM, Wolfgang Bangerth wrote:
>>>> On 05/10/2018 10:16 AM, Jonathan Perry-Houts wrote:
>>>>
>>>> Does anybody have a sense of how difficult it would be to implement 2D
>>>> cylindrical coordinates in Aspect? It doesn't look like there's support
>>>> for non-Cartesian derivatives in Deal.II, but maybe somebody has a
>>>> suggestion for hacking something together.
>>> You'd have to rewrite the assemblers to implement the bilinear form that describes the equation in cylindrical coordinates. You'd also have to hack all of the other places where we do any integrals (nullspace removal, postprocessors).
>>> It's likely going to be a substantial amount of work -- not difficult, just going through it all.
>>> Best
>>>   W.
>>
>> -- 
>> Jonathan Perry-Houts
>> Ph.D. Candidate
>> Department of Earth Sciences
>> 1272 University of Oregon
>> Eugene, OR 97403-1272
>> _______________________________________________
>> Aspect-devel mailing list
>> Aspect-devel at geodynamics.org
>> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel
> _______________________________________________
> Aspect-devel mailing list
> Aspect-devel at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel
> 

-- 
Jonathan Perry-Houts
Ph.D. Candidate
Department of Earth Sciences
1272 University of Oregon
Eugene, OR 97403-1272