[aspect-devel] Question about spherical shell subset

[FELIPE ORELLANA ROVIROSA]

 Hi Wolfgang,

  Thanks a lot for replying.

   I guess there may two different relevant cases for implementation:

  1) Construct the spherical sector (solid angle) given two (plane) angles
(a1,a2), which would be measured as apertures about perpendicular great
circles. The first aperture would be given about one pole, doing
(-a1/2,+a1/2); and the second aperture would be given about a perpendicular
pole, doing (-a2/2,+a2/2).

  2) Construct the spherical sector (solid angle) given intervals in
spherical coordinates (tet1,t1t2);(phi1,phi2)

..thinking about this..

cheers,
Felipe


On Tue, Mar 20, 2018 at 8:53 AM, Wolfgang Bangerth <bangerth at colostate.edu>
wrote:

>
> Felipe,
>
>       I am trying to implement a model in spherical coordinates, but I
>> would like to model only a sector of a spherical shell, say, a solid angle.
>> I can't find references for this.
>>
>>       I have seen this stuff in 2D, using a variable named Opening angle.
>>
>>       Can a sector be implemented in 3D?
>>
>>        it would require 2 plane opening angles..
>>
>
> It is of course possible, but it's not implemented right now either in
> ASPECT or in deal.II. You'd have to think about how you want to subdivide
> the sector into hexahedra and then write the code for it. You could look at
> how this is done for the current hyper-shell in 3d here:
>
> https://github.com/dealii/dealii/blob/master/source/grid/
> grid_generator.cc#L272
>
> or in line 3389 for the half hyper shell. You'd have to implement
> something similar for the sector. The ASPECT part of it would actually be
> quite simple then -- basically just this:
>
> https://github.com/geodynamics/aspect/blob/master/source/
> geometry_model/spherical_shell.cc#L56
>
> Best
>  W.
>
> --
> ------------------------------------------------------------------------
> Wolfgang Bangerth          email:                 bangerth at colostate.edu
>                            www: http://www.math.colostate.edu/~bangerth/
> _______________________________________________
> Aspect-devel mailing list
> Aspect-devel at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20180320/e1409a0e/attachment.html>