Re: fem order

[Konstantinos Poulios]

Dear Anne Cecile

Yes, I can confirm this. Calling set_classical_fem(2) makes tertrahedral elements into 10-node tetras, and hexahedral elements into 27-node hexas, independent of whether they have the same number of nodes in the mesh that you provide.

BR

Kostas

On Thu, Mar 24, 2022 at 4:56 PM Lesage,Anne Cecile J <AJLesage@mdanderson.org> wrote:

Dear all

 

I am writing poroelastic python script for two equations coupling pressure and displacement

If I write the following, does it mean that I am using a second order tetrahedra approximation for my displacement?

 

##########################################################################################################

#

# Definition of finite elements methods and integration method

#

#########################################################################################################

 

press_fem_order = 1      # pressure finite element order

disp_fem_order  = 2       # displacements finite element order

 

 

mfub = gf.MeshFem(meshb, 3)

mfub.set_classical_fem(disp_fem_order)

 

 

 

mimu3b=gf.MeshIm(meshb, 5)   # integration displacement on tetrahedra brain

 

 

mfpb= gf.MeshFem(meshb, 1)

mfpb.set_classical_fem(press_fem_order)

The information contained in this e-mail message may be privileged, confidential, and/or protected from disclosure. This e-mail message may contain protected health information (PHI); dissemination of PHI should comply with applicable federal and state laws. If you are not the intended recipient, or an authorized representative of the intended recipient, any further review, disclosure, use, dissemination, distribution, or copying of this message or any attachment (or the information contained therein) is strictly prohibited. If you think that you have received this e-mail message in error, please notify the sender by return e-mail and delete all references to it and its contents from your systems.