dolfin team mailing list archive

[Xavier Raynaud <question151961@xxxxxxxxxxxxxxxxxxxxx>]

New question #151961 on DOLFIN:
https://answers.launchpad.net/dolfin/+question/151961

I want to solve a mixed problem where we have one of the unknown functions which is defined in the interior of the domain and the other defined on the boundary. The two are coupled. My plan was to use mixed elements and I started to implement the following dummy problem:

from dolfin import *


mesh = UnitSquare(32,32)
bcmesh=BoundaryMesh(mesh)

# Choose finite element vector space
U = FunctionSpace(mesh, "Lagrange", 1)
P = FunctionSpace(bcmesh, "Lagrange", 1)

W=U*P

(u, p) = TrialFunctions(W)
(v, q) = TestFunctions(W)

# Construct bilinear and linear form
(u,p) = TrialFunction(W)
(v,q) = TestFunction(W)

alpha=k2=gamma=1
dt=0.1

a = u*v*dx+alpha*inner(grad(u), grad(v))*dx-k2*gamma*p*vds+(1+dt*k2)*p*q*ds


I get an error message when executing W=U*P

I wonder if Fenics is able to solve such problem. My question is actually very close to question sent 

Thank you very much for your help!

Xavier Raynaud.



-- 
You received this question notification because you are a member of
DOLFIN Team, which is an answer contact for DOLFIN.