This is a read only copy of the old FEniCS QA forum. Please visit the new QA forum to ask questions

How to Apply Pessure on Normal Direction of Boundary

0 votes

Usually, to apply the pressure P boundary condition, the pressure P is decomposed into x, y, z directions as Px, Py, Pz. Now I have a curvy boundary to apply pressure. In this case, the normal direction of the boundary need to be calculated by use: n=FacetNormal(mesh).

But I can't use P*n to define the traction boundary conditions for defining Px, Py, Pz.

Anyone knows how to define a pressure boundary condition which is normal to the boundary.

asked May 12, 2014 by wmy32552316wmy FEniCS Novice (230 points)

Pressure is a scalar and doesn't have direction. Do you want to apply a traction that is normal to the surface?

Which equation are you solving?

1 Answer

0 votes
T = Expression(('p*x[0]/R', 'p*x[1]/R'), t=10000., R = R)

thus, the surface integral in variational formular:

L = inner(T, v)*ds

in which $p$ is the pressure imposed on the boundary of a circle of radius $R$. Tell me if you need sample code (linear elasticity).

answered Oct 7, 2015 by newuser FEniCS Novice (650 points)
edited Oct 9, 2015 by newuser

But what if the boundary is not circular?

...