[Top][All Lists]

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

adding buoyancy forces to Biot poroelastic equations

From: Lesage,Anne Cecile J
Subject: adding buoyancy forces to Biot poroelastic equations
Date: Tue, 16 Nov 2021 22:12:27 +0000

Dear all


To implement the building of my fem matrix for the mechanical equilibrium, I presently write

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')


How can i add an additional buoyancy terms to the equation (see right-hand side attached equation picture)?

gravity is vector g = - 9180 N along the x axis for my mesh

rhot (density tissue is constant) but rhof = rho water for x<0 and rhof = rho air for x>0


Thank you

Anne-Cecile Lesage




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.

Attachment: mechanicaleqbuoyancy.png
Description: mechanicaleqbuoyancy.png

reply via email to

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