[Getdp] Rectangular shell transformation

[François Henrotte]

Le 13 janv. 2015 à 10:09, Tobias Nähring <i_inbox at tn-home.de> a écrit :

A small additional piece of advice.
Personally I would completely disconnect the solid structure
(your parallelepiped and cylinder) from the infinite shell.
This means leaving a small distance between them and the inner boundary shell. 
This is not strictly necessary, but I find it logical
to keep geometrical objects and modelling artifacts (the shell) independent of each other.
This will also make it easier to you to reuse the shell description in another problem. 

Fr.


> Dear Getdp-users/developers, 
> I am afraid I do not completely understand the meaning of VolRectShell. 
> I assume it transforms a shell (constructed as set difference of two co-centric axi-parallel cubes, the smaller one contained in the interior of the larger one) into the (infinitely large) complement of a cube. 
> 
> The only hint I found so far for the rectangular shell transformation is from the paper of Abert et. al.: "magnum.fe: A micromagnetic finite-element simulation code based on FEniCS" (http://arxiv.org/abs/1301.3061). 
> 
> I have attached a geometry file CuboidCyl.geo and a problem description CuboidCyl.pro. Both of them include the attached parameter file CuboidCyl.par. Note, that I am exploiting symmetry in all three space directions. 
> 
> In GuboidCyl.pro I have tried to apply VolRectShell without luck. 
> When I run getdp I get the error 
> Error   : GetDP - Bad parameters for transformation Jacobian: 0.0163227 not in [0.023,0.047571] 
> 
> Perhaps my misunderstanding of the radius parameters of VolRectShell is the origin of this error. I already tried several variants in this respect. One of those variants is implemented in CuboidCyl.pro. 
> 
> Could you please tell me what I am doing wrong? 
> 
> The problem description works (of course with a worse approximation) if I replace VolRectShell by Vol in domInf. 
> 
> With best regards, 
> Tobias Zawada 
> 
> <CuboidCyl.geo><CuboidCyl.par><CuboidCyl.pro>_______________________________________________
> getdp mailing list
> getdp at geuz.org
> http://www.geuz.org/mailman/listinfo/getdp

---
François Henrotte Dr Ir
UCL - Bâtiment Euler - Av. G. Lemaître 4-6 
B-1348 Louvain-la-Neuve
Tél: +32(0)10 47 23 64  Fax: id. 21 80