Dear all
Why shall I use a reduced fem for the pressure to solve the biot equations?
mfp_ = gf.MeshFem(mesh, 1)
mfp_.set_classical_fem(press_fem_order)
keptdofs = np.arange(mfp_.nbdof())
keptdofs = np.setdiff1d(keptdofs, mfp_.basic_dof_on_region(T_RG))
mfp = gf.MeshFem("partial", mfp_, keptdofs)
.........................
#Equation 1
md.add_linear_term(mim9,
'G*Grad(u):Grad(Test_u)+G/(1-2*nu)*Div(u)*Div(Test_u)+alpha*Grad(p).Test_u')
#Equation 2
md.add_linear_term(mim4,
'alpha/dt*(Div(u)-Div(u_prev))*Test_p+k*Grad(p).Grad(Test_p)')
Thank you
Regards
Anne-Cecile
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.
