I wanted to add that ChartManifold along with MappingQ seems to be a very general way of parametrizing geometries and work for more general shapes than circles and ellipses (see attached .png). However, the corner points of hyper_ball proved to be very important at the time of performing the transformation. The figure I present has not been optimized (initial vertexes), but serves as motivation for illustrating the point.
I believe this thread is answered now, however I cannot mark it as answered ... maybe an administrator can? Thanks again, El viernes, 7 de abril de 2017, 15:24:51 (UTC+2), Juan Carlos Araujo Cabarcas escribió: > > In the same spirit I would like to share a code that uses ChartManifold, > for generating a mapping of an ellipse applied to HyperBall. > > El miércoles, 5 de abril de 2017, 14:47:17 (UTC+2), Juan Carlos Araujo > Cabarcas escribió: >> >> Dear all, >> I recently found the thread: Something wrong with ChartManifold when >> chartdim=1: >> https://groups.google.com/forum/#!topic/dealii/Sfo9xKoeRpw >> a more straight answer to this thread. >> >> I took the liberty to copy David's code and modify it, so that the upper >> face of a square follows: y(x)=0.25*sin(PI*x) + 1. >> The function is not the same as I first proposed, but from this example >> it is pretty clear how to proceed. >> The result is plotted in the attached .png, along with a variation from >> step-10 in the deal.ii tutorial. >> >> I hope this may be of any use to somebody! >> >> Juan Carlos Araújo, >> Umeå Universitet >> El martes, 9 de febrero de 2016, 17:10:51 (UTC+1), Juan Carlos Araujo >> Cabarcas escribió: >>> >>> Dear all, >>> >>> I am working with the wave equation with variable refractive indexes >>> within the domain. >>> I receive meshes that come from a shape optimization routine, and my >>> task is to >>> run on arbitrarily curved elements resulting from the optimization. >>> Step-10 and others show how to use mappings and Manifolds to curve >>> elements following basic shapes like circles, cylinders etc. >>> >>> My guess is that it should be possible in deal.II to define my own >>> Manifolds based on how I want to bend my edges (optimization routine). >>> >>> I have prepared a very basic example on the matlab code attached that >>> has the info to plot the attached figure. There, the edges follow the >>> equation |x|^p + |y|^p=1, with p=6. >>> >>> It would be very illustrative to apply a mapping like what is done in >>> step-10 on a circle, but for the presented case. I would appreciate any >>> advise on how to achieve this, hopefully wth this example =) >>> >>> Thanks in advance! >>> >>> Juan Carlos Araújo-Cabarcas >>> Umeå Universitet. >>> >>> >>> -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.
