Navneet,
> I am solving cardiac tissue problem, i.e. a time-dependent > electro-mechanical problem. In this problem electrical potential generated > from the pacemaker is taken care by a reaction-diffusion equation, the > mechanical response is an elastic contraction. A picture is attached below > showing how a plane wave of the electrical signal passes the tissue from left > to right and how contraction follows. > > > Till now I have done only the electrical part of the problem i.e. a > reaction-diffusion problem using adaptively refined meshes. > > > For the full electro-mechanics problem, I am now coupling the electrical > potential to the mechanical response of the tissue. My doubts are following, > > > Q1. Is it possible/numerically correct to adaptively refine adeforming mesh? Yes, you can do that. You have to make sure that you move hanging nodes in such a way that they continue to be at the midpoint of the longer edge. An easy way to make sure that you do is if you deform a mesh by a displacement that is computed using a continuous finite element. (I'll note that the easier way to do this is to actually use the MappingQEulerian class -- it allows you to not actually move vertices, but just uses a finite element field as the displacement vector.) > For reaction-diffusion(electrical) part, a highly refined mesh is needed(due > to constraints on time steps) but for solving the mechanical deformation > problem coarser mesh will do, because doing mechanics calculations on the > highly refined mesh will be expensive, for efficiency I wanted to solve the > mechanics problem on a coarser version of the same mesh used for solving only > the electrical part at each time step. > > > Q2. How do I get a mesh and quadrature points of a mesh that is one or two > levels less refined than the one used for the solving purely electrical part? > how to get the quadrature points of such mesh? How do I interpolate > electrical > potential solution on this coarser mesh? Using different meshes is possible, but awkward. Have you though of just using a higher-order element for the electrical part, and a lower-order element for the displacement? This way, both variables live on the same mesh, making data exchange simpler. Best W. -- ------------------------------------------------------------------------ Wolfgang Bangerth email: [email protected] 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 [email protected]. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/ef98f8ba-a420-6654-09a4-a100429d67e0%40colostate.edu.
