State space matrices including dirichlet BC

Asked by Praveen C

I want to compute linearized state space matrices of a non-linear system, like incompressible Navier-Stokes discretized with say Taylor-Hood elements. I know how to generate the matrices to obtain a linearized system like

M*du/dt = A*u + B*p
0 = B^T *u

but I want to completely eliminate the dirichlet degrees of freedom for the velocity. If I write velocity as

u = uf + ud


uf = free dofs
ud = dirichlet dofs

then I want a system like this

Mf*d(uf)/dt = Af*uf + Ad*ud + B*p

where I simply neglect the term d(ud)/dt. So here Mf, Af are smaller matrices from which dirichlet dofs have been completely removed.

Is there a simple way to construct these matrices ?


Question information

English Edit question
DOLFIN Edit question
No assignee Edit question
Last query:
Last reply:
Revision history for this message
Johannes Ring (johannr) said :

FEniCS no longer uses Launchpad for Questions & Answers. Please consult the documentation on the FEniCS web page for where and how to (re)post your question:

Can you help with this problem?

Provide an answer of your own, or ask Praveen C for more information if necessary.

To post a message you must log in.