Hi
using the symmetric stress tensor in the Stokes part
Code:
Stokes = nu*InnerProduct(grad(u)+grad(u).trans,grad(v)) - p*div(v) - q*div(u)
of my bilinear form for a flow problem, I need to correct the homogeneous Neumann BC with
Code:
TensorCorrect = -nu*InnerProduct((grad(u).trans)*n,v)
I have added this to my bilinear form using
Code:
a += SymbolicBFI(form = TensorCorrect, VOL_or_BND=BND, skeleton=True, definedon=mesh.Boundaries("outflow"))
Unfortunately, if I attempt to also use static condensation, ( BilinearForm(..., condense=True) ), assembly of the matrix failes with the following output:
RuntimeError: SparseMatrixTM::AddElementMatrix: illegal dnumsin Assemble BilinearForm 'biform_from_py'
Adding dgjumps=True to the FESpace allows me to assemble the matrix, however I would like to avoid this loss of sparsity. Is there a better way to enable the outflow correction with static condensation?
Thank you in advance and best wishes,
Henry