[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: poroelastic 1d benchmark

From: Konstantinos Poulios
Subject: Re: poroelastic 1d benchmark
Date: Wed, 10 Nov 2021 09:14:54 +0100

Dear Anne-Cecile

This is just a way to implement a homogeneous BC on p. If you know that p has to be 0 at the top boundary (T_RG) you can use a reduced fem to remove the corresponding degrees of freedom. It is the simplest way to implement homogeneous Dirichlet BCs.


On Wed, Nov 10, 2021 at 1:05 AM Lesage,Anne Cecile J <> wrote:

Dear all


Why shall I use a reduced fem for the pressure to solve the biot equations?


mfp_ = gf.MeshFem(mesh, 1)


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







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.

reply via email to

[Prev in Thread] Current Thread [Next in Thread]