SUPG implementation
Hello everyone,
I have already posted this question in the Dolfin forum, but I guess it is rather pertinent to this forum.
I am trying to add the SUPG terms to a steady NS solver form (PSPG already working, from CBC.PDESys package).
The input objects are
u_= Previous step velocity
u = Unknown actual velocity
p = Unknown actual pressure
v = Velocity Test Function
esu = SUPG stabilization parameter (a constant)
So I define my SUPG stabilization like this
def supg_stabilizat
"""Add SUPG stabilization to Navier Stokes solver. Bilinear shapefunctions"""
esu = esu*u_
F = - esu*inner(grad(q), grad(p))*dx \
- esu*inner(grad(v), dot(grad(u), u_))*dx
for bc in self.bcs:
if bc.type() in ('ConstantPress
F = F + esu*inner(v, dot(grad(p), n))*ds(bc.bid)
return F
I am having trouble with the inner products. I know tensor dimensions are not coherent, but I do not know how to fix it.
Any ideas will be appreciated!
Thanks,
Marcos
Question information
- Language:
- English Edit question
- Status:
- Solved
- For:
- CBC.PDESys Edit question
- Assignee:
- No assignee Edit question
- Solved by:
- Marcos Samudio
- Solved:
- Last query:
- Last reply: