On 9/28/24 01:37, Praveen C wrote:

I want to map a fixed Cartesian mesh to some other domain, e.g., as in the work of Thomas Wick. So I dont want to modify the Cartesian mesh, as it is supposed to describe the reference domain.

I am constructing a MappingFEField to achieve this mapping by solving some elliptic equation.

To describe the displacement of the boundary of Cartesian mesh, I was hoping to attach a Manifold. But Manifold fixes vertex locations,

The manifold you would attach this way applies to the triangulation, which you don't want to change, not to the computed transformation. You need to prescribe your boundary displacement as the boundary conditions for the elliptic equation you describe. (In fact, you can compute this via the GridTools::laplace_transform() function.)


So I think I should just fill boundary values of the euler_vector of MappingFEField myself and proceed.

That's essentially right, except that the euler_vector is going to be the output of the elliptic solver that takes the boundary values as inputs.

Best
 W.

--
------------------------------------------------------------------------
Wolfgang Bangerth          email:                 bange...@colostate.edu
                           www: http://www.math.colostate.edu/~bangerth/


--
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.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/1f5e7d83-1571-4f14-86a7-a3595db29302%40colostate.edu.

Reply via email to